Deletion on chromosome 5 associated with decreased reproductive efficiency in female cattle

A large structural variant collection in Holstein cattle and associated database for variant discovery, characterization, and application

BACKGROUND

Structural variants (SVs) such as deletions, duplications, and insertions are known to contribute to phenotypic variation but remain challenging to identify and genotype. A more complete, accessible, and assessable collection of SVs will assist efforts to study SV function in cattle and to incorporate SV genotyping into animal evaluation.

RESULTS

In this work we produced a large and deeply characterized collection of SVs in Holstein cattle using two popular SV callers (Manta and Smoove) and publicly available Illumina whole-genome sequence (WGS) read sets from 310 samples (290 male, 20 female, mean 20X coverage). Manta and Smoove identified 31 K and 68 K SVs, respectively. In total the SVs cover 5% (Manta) and 6% (Smoove) of the reference genome, in contrast to the 1% impacted by SNPs and indels. SV genotypes from each caller were confirmed to accurately recapitulate animal relationships estimated using WGS SNP genotypes from the same dataset, with Manta genotypes outperforming Smoove, and deletions outperforming duplications. To support efforts to link the SVs to phenotypic variation, overlapping and tag SNPs were identified for each SV, using genotype sets extracted from the WGS results corresponding to two bovine SNP chips (BovineSNP50 and BovineHD). 9% (Manta) and 11% (Smoove) of the SVs were found to have overlapping BovineHD panel SNPs, while 21% (Manta) and 9% (Smoove) have BovineHD panel tag SNPs. A custom interactive database ( https://svdb-dc.pslab.ca ) containing the identified sequence variants with extensive annotations, gene feature information, and BAM file content for all SVs was created to enable the evaluation and prioritization of SVs for further study. Illustrative examples involving the genes POPDC3, ORM1, G2E3, FANCI, TFB1M, FOXC2, N4BP2, GSTA3, and COPA show how this resource can be used to find well-supported genic SVs, determine SV breakpoints, design genotyping approaches, and identify processed pseudogenes masquerading as deletions.

CONCLUSIONS

The resources developed through this study can be used to explore sequence variation in Holstein cattle and to develop strategies for studying SVs of interest. The lack of overlapping and tag SNPs from commonly used SNP chips for most of the SVs suggests that other genotyping approaches will be needed (for example direct genotyping) to understand their potential contributions to phenotype. The included SV genotype assessments point to challenges in characterizing SVs, especially duplications, using short-read data and support ongoing efforts to better characterize cattle genomes through long-read sequencing. Lastly, the identification of previously known functional SVs and additional CDS-overlapping SVs supports the phenotypic relevance of this dataset.

Multi-breed genomic predictions and functional variants for fertility of tropical bulls

Worldwide, most beef breeding herds are naturally mated. As such, the ability to identify and select fertile bulls is critically important for both productivity and genetic improvement. Here, we collected ten fertility-related phenotypes for 6,063 bulls from six tropically adapted breeds. Phenotypes were comprised of four bull conformation traits and six traits directly related to the quality of the bull's semen. We also generated high-density DNA genotypes for all the animals. In total, 680,758 single nucleotide polymorphism (SNP) genotypes were analyzed. The genomic correlation of the same trait observed in different breeds was positive for scrotal circumference and sheath score on most breed comparisons, but close to zero for the percentage of normal sperm, suggesting a divergent genetic background for this trait. We confirmed the importance of a breed being present in the reference population to the generation of accurate genomic estimated breeding values (GEBV) in an across-breed validation scenario. Average GEBV accuracies varied from 0.19 to 0.44 when the breed was not included in the reference population. The range improved to 0.28 to 0.59 when the breed was in the reference population. Variants associated with the gene HDAC4, six genes from the spermatogenesis-associated (SPATA) family of proteins, and 29 transcription factors were identified as candidate genes. Collectively these results enable very early in-life selection for bull fertility traits, supporting genetic improvement strategies currently taking place within tropical beef production systems. This study also improves our understanding of the molecular basis of male fertility in mammals.

Bovine models for human ovarian diseases

During early embryonic development, late fetal growth, puberty, adult reproductive years, and advanced aging, bovine and human ovaries closely share molecular pathways and hormonal signaling mechanisms. Other similarities between these species include the size of ovaries, length of gestation, ovarian follicular and luteal dynamics, and pathophysiology of ovarian diseases. As an economically important agriculture species, cattle are a foundational species in fertility research with decades of groundwork using physiologic, genetic, and therapeutic experimental techniques. Many technologies used in modern reproductive medicine, such as ovulation induction using hormonal therapy, were first used in cows before human trials. Human ovarian diseases with naturally occurring bovine correlates include premature ovary insufficiency (POI), polycystic ovarian syndrome (PCOS), and sex-cord stromal tumors (SCSTs). This article presents an overview of bovine ovary research related to causes of infertility, ovarian diseases, diagnostics, and therapeutics, emphasizing where the bovine model can offer advantages over other lab animals for translational applications.

Distribution of Copy Number Variation in SYT11 Gene and Its Association with Growth Conformation Traits in Chinese Cattle

Simple Summary

It is known that many different breeds of cattle are widely distributed in China. However, due to a lengthy selection of draught direction, there are obvious shortcomings in Chinese cattle, such as less meat production, slow weight gain, poor meat quality, and a lack of specialized beef cattle breeds. Animal breeding heavily benefits from molecular technologies, among which molecular genetic markers were widely used to improve the economic traits of beef cattle. Because the copy number variation (CNV) involves a longer DNA sequence or even the entire functional gene, it may have a greater impact on the phenotype. Recent studies have indicated that CNVs are widespread in the Chinese cattle genome. By investigating the effects of CNVs on gene expression and cattle traits, we aim to find those genomic variations which could significantly affect cattle traits, and which could provide a basis for genetic selection and molecular breeding of local Chinese cattle.

Abstract

Currently, studies of the SYT11 gene mainly focus on neurological diseases such as schizophrenia and Parkinson’s disease. However, some studies have shown that the C2B domain of SYT11 can interact with RISC components and affect the gene regulation of miRNA, which is important for cell differentiation, proliferation, and apoptosis, and therefore has an impact on muscle growth and development in animals. The whole-genome resequencing data detected a CNV in the SYT11 gene, and this may affect cattle growth traits. In this study, CNV distribution of 672 individuals from four cattle breeds, Yunling, Pinan, Xianan, and Qinchuan, were detected by qPCR. The relationship between CNV, gene expression and growth traits was further investigated. The results showed that the proportion of multiple copy types was the largest in all cattle breeds, but there were some differences among different breeds. The normal type had higher gene expression than the abnormal copy type. The CNVs of the SYT11 gene were significantly correlated with body length, cannon circumference, chest depth, rump length, and forehead size of Yunling cattle, and was significantly correlated with the bodyweight of Xianan cattle, respectively. These data improve our understanding of the effects of CNV on cattle growth traits. Our results suggest that the CNV of SYT11 gene is a protentional molecular marker, which may be used to improve growth traits in Chinese cattle.

Identification of Loci and Pathways Associated with Heifer Conception Rate in U.S. Holsteins

Heifer conception rate (HCR) is defined as the percentage of inseminated heifers that become pregnant at each service. The genome-wide association analyses in this study focused on identifying the loci associated with Holstein heifer (n = 2013) conception rate at first service (HCR1) and the number of times bred (TBRD) to achieve a pregnancy. There were 348 unique loci associated (p < 5 × 10⁻⁸) with HCR1 and 615 unique loci associated (p < 5 × 10⁻⁸) with TBRD. The two phenotypes shared 302 loci, and 56 loci were validated in independent cattle populations. There were 52 transcription factor binding sites (TFBS) and 552 positional candidate genes identified in the HCR1- and TBRD-associated loci. The positional candidate genes and the TFBS associated with HCR1 and TBRD were used in the ingenuity pathway analysis (IPA). In the IPA, 11 pathways, 207 master regulators and 11 upstream regulators were associated (p < 1.23 × 10⁻⁵) with HCR1 and TBRD. The validated loci associated with both HCR1 and TBRD make good candidates for genomic selection and further investigations to elucidate the mechanisms associated with subfertility and infertility.

Identification of loci associated with conception rate in primiparous Holstein cows

BACKGROUND

Subfertility is a major issue facing the dairy industry as the average US Holstein cow conception rate (CCR) is approximately 35%. The genetics underlying the physiological processes responsible for CCR, the proportion of cows able to conceive and maintain a pregnancy at each breeding, are not well characterized. The objectives of this study were to identify loci, positional candidate genes, and transcription factor binding sites (TFBS) associated with CCR and determine if there was a genetic correlation between CCR and milk production in primiparous Holstein cows. Cows were bred via artificial insemination (AI) at either observed estrus or timed AI and pregnancy status was determined at day 35 post-insemination. Additive, dominant, and recessive efficient mixed model association expedited (EMMAX) models were used in two genome-wide association analyses (GWAA). One GWAA focused on CCR at first service (CCR1) comparing cows that conceived and maintained pregnancy to day 35 after the first AI (n = 494) to those that were open after the first AI (n = 538). The second GWAA investigated loci associated with the number of times bred (TBRD) required for conception in cows that either conceived after the first AI (n = 494) or repeated services (n = 472).

RESULTS

The CCR1 GWAA identified 123, 198, and 76 loci associated (P < 5 × 10- 08) in additive, dominant, and recessive models, respectively. The TBRD GWAA identified 66, 95, and 33 loci associated (P < 5 × 10- 08) in additive, dominant, and recessive models, respectively. Four of the top five loci were shared in CCR1 and TBRD for each GWAA model. Many of the associated loci harbored positional candidate genes and TFBS with putative functional relevance to fertility. Thirty-six of the loci were validated in previous GWAA studies across multiple breeds. None of the CCR1 or TBRD associated loci were associated with milk production, nor was their significance with phenotypic and genetic correlations to 305-day milk production.

CONCLUSIONS

The identification and validation of loci, positional candidate genes, and TFBS associated with CCR1 and TBRD can be utilized to improve, and further characterize the processes involved in cattle fertility.

Extracellular vesicle-coupled miRNA profiles in follicular fluid of cows with divergent post-calving metabolic status

Most high-yielding dairy cows enter a state of negative energy balance (NEB) during early lactation. This, in turn, results in changes in the level of various metabolites in the blood and follicular fluid microenvironment which contributes to disturbed fertility. Extracellular vesicles (EVs) are evolutionarily conserved communicasomes that transport cargo of miRNA, proteins and lipids. EV-coupled miRNAs have been reported in follicular fluid. However, the association between postpartum NEB and EV-coupled miRNA signatures in follicular fluid is not yet known. Energy balance analysis in lactating cows shortly after post-calving revealed that the majority of the cows exhibited transiently negative energy balance levels, whereas the remaining cows exhibited either consistently negative or consistently positive energy levels. Metabolic status was associated with EV-coupled miRNA composition in the follicular fluid. Cows experiencing NEB showed reduced expression of a large number of miRNAs while cows with positive energy balances primarily exhibited elevated expression of EV-coupled miRNAs. The miRNAs that were suppressed under NEB were found to be involved in various metabolic pathways. This is the first study to reveal the presence of an association between EV-coupled miRNA in follicular fluid and metabolic stress in dairy cows. The involvement of differentially expressed miRNAs in various pathways associated with follicular growth and oocyte maturation suggest the potential involvement of specific follicular miRNAs in oocyte developmental competence, which may partially explain reduced fertility in cows due to post-calving metabolic stress.

Validation of 46 loci associated with female fertility traits in cattle

Background

Subfertility is one challenge facing the dairy industry as the average Holstein heifer conception rate (HCR), the proportion of heifers that conceive and maintain a pregnancy per breeding, is estimated at 55–60%. Of the loci associated with HCR, few have been validated in an independent cattle population, limiting their usefulness for selection or furthering our understanding of the mechanisms involved in successful pregnancy. Therefore, the objectives here were to identify loci associated with HCR: 1) to the first artificial insemination (AI) service (HCR1), 2) to repeated AI services required for a heifer to conceive (TBRD) and 3) to validate loci previously associated with fertility. Breeding and health records from 3359 Holstein heifers were obtained after heifers were bred by AI at observed estrus, with pregnancy determined at day 35 via palpation. Heifer DNA was genotyped using the Illumina BovineHD BeadChip, and genome-wide association analyses (GWAA) were performed with additive, dominant and recessive models using the Efficient Mixed Model Association eXpedited (EMMAX) method with a relationship matrix for two phenotypes. The HCR1 GWAA compared heifers that were pregnant after the first AI service (n = 497) to heifers that were open following the first AI service (n = 405), which included those that never conceived. The TBRD GWAA compared only those heifers which did conceive, across variable numbers of AI service (n = 712). Comparison of loci previously associated with fertility, HCR1 or TBRD were considered the same locus for validation when in linkage disequilibrium (D’ > 0.7).

Results

The HCR1 GWAA identified 116, 187 and 28 loci associated (P < 5 × 10^− 8) in additive, dominant and recessive models, respectively. The TBRD GWAA identified 235, 362, and 69 QTL associated (P < 5 × 10^− 8) with additive, dominant and recessive models, respectively. Loci previously associated with fertility were in linkage disequilibrium with 22 loci shared with HCR1 and TBRD, 5 HCR1 and 19 TBRD loci.

Conclusions

Loci associated with HCR1 and TBRD that have been identified and validated can be used to improve HCR through genomic selection, and to better understand possible mechanisms associated with subfertility.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5935-3) contains supplementary material, which is available to authorized users.

Copy number variation of bovine SHH gene is associated with body conformation traits in Chinese beef cattle

Sonic Hedgehog (Shh) regulates many key developmental processes during vertebrate limb development, fat formation, and skeletal tissue regeneration. Current whole genome sequencing data have identified a copy number variation mapping to bovine Sonic Hedgehog gene (SHH-CNV). The object of this study was to characterize the SHH-CNV distributions in 648 individuals from 11 Chinese cattle populations and further to investigate the associations of the copy number changes with gene expression and cattle growth traits. The SHH-CNV showed a high variance within Chinese indigenous yellow cattle. Compared to yak and dairy cattle, the beef cattle like Luxi and Xianan breed had significantly higher median copy numbers, suggesting the diversity of SHH-CNV in beef cattle selections. The negative correlation of SHH-CNV with SHH transcriptional level in adult adipose tissue (P < 0.01) indicated the dosage effects of SHH-CNV related to bovine fat formation. Association analysis of SHH-CNV and body size traits was conducted in five breeds. The results revealed that the copy number gain type of SHH-CNV exhibited significantly better chest depth in 24 months old Qinchuan cattle, and better body weight, body length, and chest girth in 18 months old Nanyang cattle, whereas the normal copy number had superior chest girth and body weight in adult Jinnan cattle (P < 0.05 or P < 0.01). In summary, this research uncovered meaningful effects of SHH-CNV on gene expression and cattle phenotypic traits, indicating its potential applications for genetic improvement of beef cattle.

Whole-genome analysis of structural variations between Xiang pigs with larger litter sizes and those with smaller litter sizes

To gain a better knowledge of structural variations (SVs) in Xiang pig, we used next-generation sequencing to analyze the Xiang pigs with larger (XL) or smaller litter sizes (XS). Our analysis yielded 28,040 putative SVs in the Xiang pig. These SVs distributed throughout all of chromosomes. Some functional regions including exons and untranslated regions were less varied than introns and intergenic regions. We detected 4637 and 4119 specific SVs, which contained 1697 and 1582 genes in XL and XS group, respectively. These genes were mainly enriched in the well-known pathways involved in development and reproduction processes. Population validation was carried out on 50 SVs candidates using PCR method in 144 Xiang pig crowds. All of 50 SVs were confirmed by PCR method and 14 SVs were associated with the litter size of Xiang pigs. These results may be helpful for the elucidation of growth and reproduction regulation in Xiang pig.

The Bos taurus-Bos indicus balance in fertility and milk related genes

Numerical approaches to high-density single nucleotide polymorphism (SNP) data are often employed independently to address individual questions. We linked independent approaches in a bioinformatics pipeline for further insight. The pipeline driven by heterozygosity and Hardy-Weinberg equilibrium (HWE) analyses was applied to characterize Bos taurus and Bos indicus ancestry. We infer a gene co-heterozygosity network that regulates bovine fertility, from data on 18,363 cattle with genotypes for 729,068 SNP. Hierarchical clustering separated populations according to Bos taurus and Bos indicus ancestry. The weights of the first principal component were subjected to Normal mixture modelling allowing the estimation of a gene’s contribution to the Bos taurus-Bos indicus axis. We used deviation from HWE, contribution to Bos indicus content and association to fertility traits to select 1,284 genes. With this set, we developed a co-heterozygosity network where the group of genes annotated as fertility-related had significantly higher Bos indicus content compared to other functional classes of genes, while the group of genes associated with milk production had significantly higher Bos taurus content. The network analysis resulted in capturing novel gene associations of relevance to bovine domestication events. We report transcription factors that are likely to regulate genes associated with cattle domestication and tropical adaptation. Our pipeline can be generalized to any scenarios where population structure requires scrutiny at the molecular level, particularly in the presence of a priori set of genes known to impact a phenotype of evolutionary interest such as fertility.

Genomic analysis of stayability in Nellore cattle

Stayability, which can be defined as the probability of a cow calving at a certain age when given the opportunity, is an important reproductive trait in beef cattle because it is directly related to herd profitability. The objective of this study was to estimate genetic parameters and to identify possible genomic regions associated with the phenotypic expression of stayability in Nellore cows. The variance components were estimated by Bayesian inference using a threshold animal model that included the systematic effects of contemporary group and sexual precocity and the random effects of animal and residual. The SNP effects were estimated by the single-step genomic BLUP method using information of 2,838 animals (2,020 females and 930 sires) genotyped with the Illumina High-Density BeadChip Array (San Diego, CA, USA). The variance explained by windows formed by 200 consecutive SNPs was used to identify genomic regions of largest effect on the expression of stayability. The heritability was 0.11 ± 0.01 when A matrix (pedigree) was used and 0.14 ± 0.01 when H matrix (relationship matrix that combines pedigree information and SNP data) was used. A total of 147 candidate genes for stayability were identified on chromosomes 1, 2, 5, 6, 9 and 20 and on the X chromosome. New candidate regions for stayability were detected, most of them related to reproductive, immunological and central nervous system functions.

Comparative analyses across cattle genders and breeds reveal the pitfalls caused by false positive and lineage-differential copy number variations

We compared CNV region (CNVR) results derived from 1,682 Nellore cattle with equivalent results derived from our previous analysis of Bovine HapMap samples. By comparing CNV segment frequencies between different genders and groups, we identified 9 frequent, false positive CNVRs with a total length of 0.8 Mbp that were likely caused by assembly errors. Although there was a paucity of lineage specific events, we did find one 54 kb deletion on chr5 significantly enriched in Nellore cattle. A few highly frequent CNVRs present in both datasets were detected within genomic regions containing olfactory receptor, ATP-binding cassette, and major histocompatibility complex genes. We further evaluated their impacts on downstream bioinformatics and CNV association analyses. Our results revealed pitfalls caused by false positive and lineage-differential copy number variations and will increase the accuracy of future CNV studies in both taurine and indicine cattle.

Genome-Wide Study of Structural Variants in Bovine Holstein, Montbéliarde and Normande Dairy Breeds

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.