Integrating RNA-seq and GWAS reveals novel genetic mutations for buffalo reproductive traits

Genome-wide association study (GWAS) has been applied in buffalo breeding programs and been used to identify a number of candidate genes associated with buffalo reproductive traits. The genetic code of specific genes underlying buffalo reproductive traits remains unclear. Association study that measures both genetic and transcriptional variation has been applied for the investigation of complex traits. To investigate genes involved in buffalo reproductive traits, integrated RNA-seq results were investigated of buffalo granulosa cells and candidate genes which were reported to be associated with buffalo reproductive traits in a previous GWAS. A large number of variants were detected by RNA-seq, and 214 variants were located within the buffalo reproductive candidate genes identified by GWAS. A further association study in 462 Italian Mediterranean buffalo indicated that 25 SNPs distributed in 13 genes were associated with reproductive traits. Of the 13 genes, 11 were expressed in granulosa cells of all antral follicle development stages, and significant difference was found in the expression of NDUFS2 between follicles of diameter<8mm and>8 mm. These findings extend the results of GWAS by expanding the knowledge about new and potentially functional single-nucleotide polymorphisms and provide useful information about regulatory genes affecting buffalo reproductive traits.