Identification of differentially expressed genes in testis
We detected differentially expressed genes (DEGs) of testis using both DESeq2 v1.34.0 (Love et al., 2014) and edgeR v3.36.0 (Robinson et al., 2010) to compare expression levels between the BY and LJX tree sparrows. We used a threshold value of fold change of 2 and adjusted Pvalue (Benjamin-Hochberg correction in DESeq2 and FDR correction in edgeR) of 0.05 for differential expression analysis. Only genes identified by both DESeq2 and edgeR were considered to be DEGs and were used for further analysis.