Detection of polymorphisms in six genes and their association analysis with litter size in sheep

Litter size in sheep is a complex trait controlled by micro-effective polygenes. APAF1, CLSTN2, CTH, PLCB1, PLCB4, and CHST11 are all involved in mammalian reproduction. However, the effects of these genes on litter size in sheep are still unclear. Therefore, in this study, we used Sequenom MassARRAY® SNP assay technology to type the single nucleotide polymorphisms (SNPs) loci of six genes in five sheep breeds. The results showed that most sheep breeds contain three genotypes at each locus. Then, we conducted population genetic analysis on the SNPs of six genes and found that the polymorphic information content in all sheep breeds ranged from 0 to 0.37, and most sheep breeds were in Hardy-Weinberg equilibrium (p > 0.05). In addition, association analysis in Small Tail Han sheep indicated that the rs399534524 locus in CLSTN2 was highly associated with first parity litter size, and litter size in ewes with CT genotype was higher than that in ewes with CC genotype or TT genotype. Furthermore, the rs407142552 locus in CTH was highly associated with second parity litter size in Small Tail Han sheep, and litter size in ewes with CT genotype was higher than that in ewes with TT genotype. Finally, we predicted the CTH and CLSTN2 protein interaction network and found that HTR1E, NOM1, CCDC174 and ALPK3 interact with CLSTN2 and have been reported as candidate genes related to litter size in sheep. These results suggest that they may be useful genetic markers for increasing litter size in sheep.

Boar semen cryopreservation: State of the art, and international trade vision

Biosecurity is a major concern in the global pig production. The separation in time of semen collection, processing and insemination in the pig farm is a few days for chilled semen but it can be indefinite when using cryopreserved semen. Field fertility results of boar cryopreserved semen are close to chilled semen, which makes it a valuable resource for the establishment of semen genebanks, long-distance semen trade, and the implementation of other technologies such as the sex-sorted semen. But cryopreserved semen is far from being routine in pig farms. The most recent research efforts to facilitate its implementation include the use of additives before freezing, or in the thawing extender. Long-term preserved semen trade is a biosecurity challenge. To harmonize international trade of germplasm, the World Organization of Animal Health (WOAH) established a regulatory framework for all member countries. The present paper aims to review the latest advances of boar semen cryopreservation with special focus on the benefits of its inclusion as a routine tool in the pig industry. We also review recently reported field fertility results of cryopreserved semen, its international trade compared to chilled semen, and the regulatory framework involved. Boar cryopreserved semen is a valuable tool to control biosecurity risk, implement other technologies, and facilitate international trade. Research already demonstrated good field fertility results, but it still represents less than 0.1 % of the international trade. As boar cryopreserved semen gets closer to implementation, the correspondent authorities are reviewing the trade rules.

Identification of genetic variants affecting reproduction traits in Vrindavani cattle

Genome-wide association studies (GWAS) are one of the best ways to look into the connection between single-nucleotide polymorphisms (SNPs) and the phenotypic performance. This study aimed to identify the genetic variants that significantly affect the important reproduction traits in Vrindavani cattle using genome-wide SNP chip array data. In this study, 96 randomly chosen Vrindavani cows were genotyped using the Illumina Bovine50K BeadChip platform. A linear regression model of the genome-wide association study was fitted in the PLINK program between genome-wide SNP markers and reproduction traits, including age at first calving (AFC), inter-calving period (ICP), dry days (DD), and service period (SP) across the first three lactations. Information on different QTLs and genes, overlapping or adjacent to genomic coordinates of significant SNPs, was also mined from relevant databases in order to identify the biological pathways associated with reproductive traits in bovine. The Bonferroni correction resulted in total 39 SNP markers present on different chromosomes being identified that significantly affected the variation in AFC (6 SNPs), ICP (7 SNPs), DD (9 SNPs), and SP (17 SNPs). Novel potential candidate genes associated with reproductive traits that were identified using the GWAS methodology included UMPS, ITGB5, ADAM2, UPK1B, TEX55, bta-mir-708, TMPO, TDRD5, MAPRE2, PTER, AP3B1, DPP8, PLAT, TXN2, NDUFAF1, TGFA, DTNA, RSU1, KCNQ1, ADAM32, and CHST8. The significant SNPs and genes associated with the reproductive traits and the enriched genes may be exploited as candidate biomarkers in animal improvement programs, especially for improved reproduction performance in bovines.

Genome-wide in silico analysis leads to identification of deleterious L290V mutation in RBBP5 gene in Bos indicus

Genome-wide deleterious mutations were identified in zebu cattle (Bos indicus) using in silico approach. The ddRAD sequence data of Sahiwal cattle were annotated and aligned with the cattle reference genome (ARS-UCD1.2). A total of 279,383 SNPs were identified at Read Depth10, which were further filtered to 692 missense SNPs. These SNPs were further analyzed, for functional consequences, by using Variant Effect Predictor, PolyPhen, PROVEAN, and PANTHER tools. A total of 18 SNPs, were finally identified as deleterious, and among these, 12 SNPs were mapped on nine different genes. ERRAT, ProSA-web, Project HOPE, TM-Align, and YASSARA tools, further confirmed the protein malfunctioning of one missense (L290V) mutation of Retinoblastoma binding protein-5 (RBBP5) gene, transcribing a cell cycle regulatory protein and associated with Retinoblastoma in human. This derived bioinformatics pipeline may be useful for preliminarily identifying the deleterious DNA mutations in livestock, specifically in absence of any genetic disease records.

Role of genetic introgression in introducing mutant alleles in Bos indicus cattle and prevalence of lethal genetic disorders in Bos taurus × Bos indicus and Bos indicus cattle in India

Fertility is an important trait associated with reproductive performance and animal welfare concern. Lethal alleles affect fertility through early embryonic death, abortions, and stillbirth depending on the genetic expression of the allele. Holstein Friesian and Jersey are two major Bos taurus breeds used widely for increasing milk yield along with purebreds of Bos indicus breeds like Gir, Kankrej, Sahiwal, and Tharparkar. In the present study, prevalence of lethal mutants in crossbred Holstein Friesian (CBHF, n = 2435), crossbred Jersey (CBJY, n = 2874), Gir (n = 3288), Kankrej (n = 593), Sahiwal (n = 965), and Tharparkar (n = 18) were studied. Heterozygous carrier animals were identified for bovine leukocyte adhesion deficiency (BLAD), Citrullinemia, complex vertebral malformation (CVM), Brachyspina, Holstein Haplotype 1 (HH1), Holstein Haplotype 3 (HH3),Holstein Haplotype 4 (HH4) and Jersey Haplotype 1 (JH1). Breed purity analysis confirmed inheritance of Bos taurus genes contributing to the presence of lethal mutant alleles like BLAD, Citrullinemia, HH1, and JH1 in apparently phenotypic Bos indicus animals. Screening and elimination of heterozygous carrier bulls/cows is essential to control fertility loss associated with lethal alleles.

Case report: Investigation of genetic mutations in a case of schistosomus reflexus in a Holstein dairy cattle fetus in Korea

Schistosomus reflexus (SR) is one of the most common congenital anomalies found in cases of cattle dystocia; this disorder occurs mostly in cattle. Congenital anomalies such as SR are caused by various genetic and environmental factors, but no specific cause has been elucidated for SR. This study reports a case of SR in a Holstein dairy cattle fetus with congenital anomalies in Korea. Grossly, a distinct spine curvature was observed between the thoracic and lumbar vertebrae, accompanied by a consequential malformation from the sacrum to the occipital bone. Furthermore, the thoracic and abdominal organs were exposed. In computed tomography (CT) images, mild and severe kyphoscoliosis was observed in T1~11 and L1~6, respectively. Additionally, vertebral dysplasia was observed in S1~5 and Cd 1~5. To pinpoint the causal genes and mutations, we leveraged a custom 50K Hanwoo SNP-Chip and the Online Mendelian Inheritance in Animals (OMIA) database. As a result, we identified a nonsense mutation in apoptotic protease activating factor 1 (APAF1) within HH1 that was associated with a decrease in conception rate and an increase in abortion in Holstein dairy cattle. The genotype of the SR case was A/A, and most of the 1,142 normal Holstein dairy cattle tested as a control group had the genotype G/G. In addition, the A/A genotype did not exist in the control group. Based on the pathological, genetic, and radiological findings, the congenital abnormalities observed were diagnosed as SR.

Unravelling transmission ratio distortion across the bovine genome: identification of candidate regions for reproduction defects

BACKGROUND

Biological mechanisms affecting gametogenesis, embryo development and postnatal viability have the potential to alter Mendelian inheritance expectations resulting in observable transmission ratio distortion (TRD). Although the discovery of TRD cases have been around for a long time, the current widespread and growing use of DNA technologies in the livestock industry provides a valuable resource of large genomic data with parent-offspring genotyped trios, enabling the implementation of TRD approach. In this research, the objective is to investigate TRD using SNP-by-SNP and sliding windows approaches on 441,802 genotyped Holstein cattle and 132,991 (or 47,910 phased) autosomal SNPs.

RESULTS

The TRD was characterized using allelic and genotypic parameterizations. Across the whole genome a total of 604 chromosomal regions showed strong significant TRD. Most (85%) of the regions presented an allelic TRD pattern with an under-representation (reduced viability) of carrier (heterozygous) offspring or with the complete or quasi-complete absence (lethality) for homozygous individuals. On the other hand, the remaining regions with genotypic TRD patterns exhibited the classical recessive inheritance or either an excess or deficiency of heterozygote offspring. Among them, the number of most relevant novel regions with strong allelic and recessive TRD patterns were 10 and 5, respectively. In addition, functional analyses revealed candidate genes regulating key biological processes associated with embryonic development and survival, DNA repair and meiotic processes, among others, providing additional biological evidence of TRD findings.

CONCLUSIONS

Our results revealed the importance of implementing different TRD parameterizations to capture all types of distortions and to determine the corresponding inheritance pattern. Novel candidate genomic regions containing lethal alleles and genes with functional and biological consequences on fertility and pre- and post-natal viability were also identified, providing opportunities for improving breeding success in cattle.

Identification of candidate lethal haplotypes and genomic association with post-natal mortality and reproductive traits in Nellore cattle

The wide use of genomic information has enabled the identification of lethal recessive alleles that are the major genetic causes of reduced conception rates, longer calving intervals, or lower survival for live-born animals. This study was carried out to screen the Nellore cattle genome for lethal recessive haplotypes based on deviation from the expected population homozygosity, and to test SNP markers surrounding the lethal haplotypes region for association with heifer rebreeding (HR), post-natal mortality (PNM) and stayability (STAY). This approach requires genotypes only from apparently normal individuals and not from affected embryos. A total of 62,022 animals were genotyped and imputed to a high-density panel (777,962 SNP markers). Expected numbers of homozygous individuals were calculated, and the probabilities of observing 0 homozygotes was obtained. Deregressed genomic breeding values [(G)EBVs] were used in a GWAS to identify candidate genes and biological mechanisms affecting HR, STAY and PNM. In the functional analyses, genes within 100 kb down and upstream of each significant SNP marker, were researched. Thirty haplotypes had high expected frequency, while no homozygotes were observed. Most of the alleles present in these haplotypes had a negative mean effect for PNM, HR and STAY. The GWAS revealed significant SNP markers involved in different physiological mechanisms, leading to harmful effect on the three traits. The functional analysis revealed 26 genes enriched for 19 GO terms. Most of the GO terms found for biological processes, molecular functions and pathways were related to tissue development and the immune system. More phenotypes underlying these putative regions in this population could be the subject of future investigation. Tests to find putative lethal haplotype carriers could help breeders to eliminate them from the population or manage matings in order to avoid homozygous.

Mating allocations in Holstein combining genomic information and linear programming optimization at the herd level

In this study, we explored mating allocation in Holstein using genomic information for 24,333 Holstein females born in Denmark, Finland, and Sweden. We used 2 data sets of bulls: the top 50 genotyped bulls and the top 25 polled genotyped bulls on the Nordic total merit scale. We used linear programming to optimize economic scores within each herd, considering genetic level, genetic relationship, semen cost, the economic impact of genetic defects, polledness, and β-casein. We found that it was possible to reduce genetic relationships and eliminate expression of genetic defects with minimal effect on the genetic level in total merit index. Compared with maximizing only Nordic total merit index, the relative frequency of polled offspring increased from 13.5 to 22.5%, and that of offspring homozygous for β-casein (A2A2) from 66.7 to 75.0% in one generation, without any substantial negative impact on other comparison criteria. Using only semen from polled bulls, which might become necessary if dehorning is banned, considerably reduced the genetic level. We also found that animals carrying the polled allele were less likely to be homozygous for β-casein (A2A2) and more likely to be carriers of the genetic defect HH1. Hence, adding economic value to a monogenic trait in the economic score used for mating allocation sometimes negatively affected another monogenetic trait. We recommend that the comparison criteria used in this study be monitored in a modern genomic mating program.

Whole-genome resequencing reveals genetic diversity and selection characteristics of dairy goat

The dairy goat is one of the earliest dairy livestock species, which plays an important role in the economic development, especially for developing countries. With the development of agricultural civilization, dairy goats have been widely distributed across the world. However, few studies have been conducted on the specific characteristics of dairy goat. In this study, we collected the whole-genome data of 89 goat individuals by sequencing 48 goats and employing 41 publicly available goats, including five dairy goat breeds (Saanen, Nubian, Alpine, Toggenburg, and Guanzhong dairy goat; n = 24, 15, 11, 6, 6), and three goat breeds (Guishan goat, Longlin goat, Yunshang Black goat; n = 6, 15, 6). Through compared the genomes of dairy goat and non-dairy goat to analyze genetic diversity and selection characteristics of dairy goat. The results show that the eight goats could be divided into three subgroups of European, African, and Chinese indigenous goat populations, and we also found that Australian Nubian, Toggenburg, and Australian Alpine had the highest linkage disequilibrium, the lowest level of nucleotide diversity, and a higher inbreeding coefficient, indicating that they were strongly artificially selected. In addition, we identified several candidate genes related to the specificity of dairy goat, particularly genes associated with milk production traits (GHR, DGAT2, ELF5, GLYCAM1, ACSBG2, ACSS2), reproduction traits (TSHR, TSHB, PTGS2, ESR2), immunity traits (JAK1, POU2F2, LRRC66). Our results provide not only insights into the evolutionary history and breed characteristics of dairy goat, but also valuable information for the implementation and improvement of dairy goat cross breeding program.

Identification of a putative haplotype associated with recumbency in Holstein calves

Thirty-four Holstein calves from multiple farms were found recumbent during the neonatal period with no detectable neurologic, infectious, or metabolic abnormalities. Most calves did not survive beyond 6 wk of age. The objective of this study was to conduct a genome-wide association and pedigree analysis to determine if a genetic origin was plausible. There were 101,917 DNA markers for 18 affected calves and 26 unaffected family controls available for analysis. Genome-wide association, homozygosity screening, and a parental based transmission disequilibrium test were conducted in PLINK. A genomic region on the end of chromosome 16 that contained 78 markers based on a recessive inheritance model and that spanned 5.1 million bp was considered the most probable region for a genetic defect; the region was narrowed to 2.1 million bp following homozygosity screening and the transmission disequilibrium test with all affected calves homozygous in the candidate region and 1 homozygous control. A genotyped sire and 2 dams with imputed genotypes were heterozygous in the candidate region. A common sire born in 2008 was identified that was present for both paternal and maternal lineages of all affected calves; nearly all lineages traced through a prolific son born in 2010 who was genotyped and was heterozygous for the candidate region. Therefore, a possible genetic defect with incomplete penetrance on chromosome 16 that results in recumbency has been identified. Further efforts with an increase in families represented are needed to confirm a genetic basis, and identify the mutation and mode of inheritance.

Four novel candidate causal variants for deficient homozygous haplotypes in Holstein cattle

Mendelian variants can determine both insemination success and neonatal survival and thus influence fertility and rearing success of cattle. We present 24 deficient homozygous haplotype regions in the Holstein population of Switzerland and provide an overview of the previously identified haplotypes in the global Holstein breed. This study encompasses massive genotyping, whole-genome sequencing (WGS) and phenotype association analyses. We performed haplotype screenings on almost 53 thousand genotyped animals including 114 k SNP data with two different approaches. We revealed significant haplotype associations to several survival, birth and fertility traits. Within haplotype regions, we mined WGS data of hundreds of bovine genomes for candidate causal variants, which were subsequently evaluated by using a custom genotyping array in several thousand breeding animals. With this approach, we confirmed the known deleterious SMC2:p.Phe1135Ser missense variant associated with Holstein haplotype (HH) 3. For two previously reported deficient homozygous haplotypes that show negative associations to female fertility traits, we propose candidate causative loss-of-function variants: the HH13-related KIR2DS1:p.Gln159* nonsense variant and the HH21-related NOTCH3:p.Cys44del deletion. In addition, we propose the RIOX1:p.Ala133_Glu142del deletion as well as the PCDH15:p.Leu867Val missense variant to explain the unexpected low number of homozygous haplotype carriers for HH25 and HH35, respectively. In conclusion, we demonstrate that with mining massive SNP data in combination with WGS data, we can map several haplotype regions and unravel novel recessive protein-changing variants segregating at frequencies of 1 to 5%. Our findings both confirm previously identified loci and expand the spectrum of undesired alleles impairing reproduction success in Holstein cattle, the world's most important dairy breed.

Allele Frequency of APAF1 Mutation in Holstein Cattle in Brazil

APAF1 is an autosomal recessive inherited mutation, associated with Holstein haplotype 1 (HH1) and characterized by a substitution of cytosine for a thymine (c.1741C>T) in chromosome 5. The mutation causes fetal and embryonic loss, between 60 and 200 days of gestation, and reduced conception rate. The ARMS-PCR is considered a simple and low-cost method to determine single nucleotide polymorphism (SNP) with no need for genetic sequencing of the animal genome. This study aimed to verify the allelic frequency of APAF1 mutation in Brazilian Holstein cattle. A total of 248 Holstein DNA samples (210 cows and 38 bulls) were analyzed, and synthetic genes were manufactured to validate the primers developed by the authors. All animals assessed in this study were classified as wild-type for APAF1 mutation. The primers and protocol developed for the ARMS-PCR technique work with 100% specificity and efficiency since the amplicon formations are as expected according to the genotypes. In conclusion, the mutation responsible for APAF1 was not detected in the Brazilian Holstein cattle population assessed in this prevalence study, although it is not possible to affirm that APAF1 does not occur in Brazilian Holstein animals. The tetra-primer ARMS-PCR protocol for APAF1 mutation that has been validated here may be a relatively simple and economical method to determine the animals' genotype.

A Nonsense Variant in CCDC65 Gene Causes Respiratory Failure Associated with Increased Lamb Mortality in French Lacaune Dairy Sheep

We recently demonstrated that the Lacaune deficient homozygous haplotype 6 (LDHH6) potentially hosts a recessive perinatal lethal mutation in Lacaune dairy sheep mapped on OAR3. In the present study, we have analyzed the whole-genome sequences of two Lacaune ram heterozygous carriers of LDHH6. After variant calling and filtering against the variants of 86 non-carrier rams, we have identified a single nucleotide variant (SNV) in the two LDHH6 carriers whose variant allele induced a premature stop codon (p.Glu111*) in the Coiled-Coil Domain Containing 65 (CCDC65) gene. CCDC65 is involved in the assembly of the nexin-dynein regulatory complex for the formation of microtubules in ciliated cells. In order to identify the phenotype in homozygous sheep, we generated at-risk matings (n = 17) between rams and ewes heterozygous for the candidate variant in CCDC65. A total of 16 lambs were born alive with five genotyped as homozygous carriers. The homozygous lambs suffered from respiratory problems, and four of them died within the first month of life. At necropsy, we observed a broad hepatization of lung lobes possibly induced by infectious pneumonia. The management of this lethal recessive allele (frequency of 0.06) through reasoned mating in the Lacaune sheep selection schemes could reduce lamb mortality by 2%.

Analysis of Genetic Diversity in the American Standardbred Horse Utilizing Short Tandem Repeats and Single Nucleotide Polymorphisms

American Standardbreds were developed as a harness racing horse breed. The United States Trotting Association closed the studbook in 1973 and implemented a book size cap in 2009. This study aimed to investigate genetic diversity in the American Standardbred after the studbook cap was introduced using short tandem repeats (STRs) and single-nucleotide polymorphisms (SNPs). Sixteen STRs from horses foaled from 2010 to 2015 and their sires and dams (n = 50 621) were utilized to examine allelic richness (A_r), expected heterozygosity (H_E), observed heterozygosity (H_O), unbiased heterozygosity (H_U), inbreeding coefficient (F_IS), and fixation index (F_ST). These analyses found that trotting and pacing sires were less genetically diverse than dams (H_EP_Bonferroni = 0.029 and 6.3 × 10⁻⁵, respectively) and their offspring (A_rP_Bonferroni = 0.034 and 6.9 × 10^-6, respectively), and pacing offspring were significantly less diverse than their dams (H_EP_Bonferroni = 2 × 10^-3). Inbreeding coefficients for trotters (F_IS = −0.014) and pacers (F_IS = −0.012) suggest that breeding practices have maintained diversity. Moderate levels of genetic differentiation (0.066 < F_ST < 0.11) were found between pacing and trotting groups. Additionally, 10 of the most prolific trotting sires and their male offspring (n = 84) were genotyped on the 670K Axiom Equine HD Array. H_O values higher than H_E (P < 0.001), low inbreeding coefficients (mean F = −0.064), and mean F_ROH = 21% indicate relatively high levels of diversity in this cohort, further supporting the STR data. However, in contrast, H_O values were higher for trotting sires (0.41) than their offspring (0.36). This observation warrants further monitoring of diversity over time. These data provide an updated foundation of diversity indices for further, long-term analysis in the breed.

Study on the concordance between different SNP-genotyping platforms in sheep

Different SNP genotyping technologies are commonly used in multiple studies to perform QTL detection, genotype imputation, and genomic predictions. Therefore, genotyping errors cannot be ignored, as they can reduce the accuracy of different procedures applied in genomic selection, such as genomic imputation, genomic predictions, and false-positive results in genome-wide association studies. Currently, whole-genome resequencing (WGR) also offers the potential for variant calling analysis and high-throughput genotyping. WGR might overshadow array-based genotyping technologies due to the larger amount and precision of the genomic information provided; however, its comparatively higher price per individual still limits its use in larger populations. Thus, the objective of this work was to evaluate the accuracy of the two most popular SNP-chip technologies, namely, Affymetrix and Illumina, for high-throughput genotyping in sheep considering high-coverage WGR datasets as references. Analyses were performed using two reference sheep genome assemblies, the popular Oar_v3.1 reference genome and the latest available version Oar_rambouillet_v1.0. Our results demonstrate that the genotypes from both platforms are suggested to have high concordance rates with the genotypes determined from reference WGR datasets (96.59% and 99.51% for Affymetrix and Illumina technologies, respectively). The concordance results provided in the current study can pinpoint low reproducible markers across multiple platforms used for sheep genotyping data. Comparing results using two reference genome assemblies also informs how genome assembly quality can influence genotype concordance rates among different genotyping platforms. Moreover, we describe an efficient pipeline to test the reliability of markers included in sheep SNP-chip panels against WGR datasets available on public databases. This pipeline may be helpful for discarding low-reliability markers before exploiting genomic information for gene mapping analyses or genomic prediction.

Bayesian inference of the inbreeding load variance for fertility traits in Brown Swiss cattle

Our study investigated the inbreeding load for fertility traits in the Italian Brown Swiss dairy cattle breed. Fertility traits included continuous traits (i.e., interval from calving to first service, days open, and calving interval) and categorical traits (i.e., calving rate at first insemination and nonreturn date at d 56). We included only records of the first 3 parities of cows that calved between 2010 and 2018. We traced up the pedigree of the cows with records as far as possible, ending up with a total of 73,246 animals. The final data set consisted of 59,864 records from 34,921 cows. We analyzed all models using a Bayesian approach that included a covariate with total inbreeding in addition to systematic, permanent environment, additive genetic, and inbreeding load effects. We then evaluated the trends in heritabilities and ratios of the inbreeding load using a continuum of partial inbreeding coefficients from 0.001 to 0.100 as reference. Posterior estimates of heritabilities tended to decrease across the continuum, whereas ratios of the inbreeding load tended to increase, more noticeably in categorical traits (calving rate at first insemination and nonreturn date at d 56). From the results obtained, we confirmed the presence of heterogeneity in inbreeding depression. We then predicted the inbreeding load effects, which had a low reliability of prediction, explained by having only 513 ancestors generating inbreeding. However, reliability of prediction was high enough for some of the individuals, obtaining a favorable prediction of inbreeding load for a relevant percentage, which improved the phenotypic performance of their inbred descendants. These results make it feasible to implement breeding and management strategies that select ancestors with a favorable inbreeding load prediction. In addition, it opens the possibility to define a global index for the expected consequences of the inbreeding generated by each individual.

Non-additive association analysis using proxy phenotypes identifies novel cattle syndromes

Mammalian species carry ~100 loss-of-function variants per individual^1,2, where ~1-5 of these impact essential genes and cause embryonic lethality or severe disease when homozygous³. The functions of the remainder are more difficult to resolve, although the assumption is that these variants impact fitness in less manifest ways. Here we report one of the largest sequence-resolution screens of cattle to date, targeting discovery and validation of non-additive effects in 130,725 animals. We highlight six novel recessive loci with impacts generally exceeding the largest-effect variants identified from additive genome-wide association studies, presenting analogs of human diseases and hitherto-unrecognized disorders. These loci present compelling missense (PLCD4, MTRF1 and DPF2), premature stop (MUS81) and splice-disrupting (GALNT2 and FGD4) mutations, together explaining substantial proportions of inbreeding depression. These results demonstrate that the frequency distribution of deleterious alleles segregating in selected species can afford sufficient power to directly map novel disorders, presenting selection opportunities to minimize the incidence of genetic disease.

Prevalence of nine genetic defects in Chinese Holstein cattle

Worldwide use of elite sires has caused inbreeding accumulation and high frequencies of genetic defects in dairy cattle populations. In recent years, several genetic defect genes or haplotypes have been identified in Holstein cattle. A rapid and reliable microfluidic chip with Kompetitive allele-specific PCR (KASP) assay was developed in our previous study for the detection of heterozygotes at eight genetic defect loci of bovine leukocyte adhesion deficiency (BLAD), Brachyspina syndrome (BS), complex vertebral malformation (CVM), Holstein haplotype 1 (HH1), Holstein haplotype 3 (HH3), Holstein haplotype 4 (HH4), Holstein haplotype 5 (HH5) and haplotype for cholesterol deficiency (HCD). This study aimed to extend that assay to include a newly identified genetic defect of Holstein haplotype 6 (HH6) and to estimate the frequencies of carriers for each of the nine genetic defects in six Chinese Holstein herds. Of the 1633 cows, carrier frequencies of the genetic defects were 6.92%, 5.76%, 4.46%, 4.30%, 3.62%, 2.94%, 1.86% and 0.37% for HH1, HH3, CVM, HH5, HCD, BS, HH6 and BLAD, respectively. No carrier was found for HH4. Notably, 27.43% of cows carried at least one genetic defect, while 2.27% and 0.12% of cows carried double and triple genetic defect alleles, respectively. The existence of genetic defects calls for routine molecular testing and effective management of genetic defects by avoiding carrier-to-carrier mating in production herds and eliminating or at least reducing the frequency of the defective alleles through marker-assisted selection in breeding herds.

Identification of homozygous haplotypes carrying putative recessive lethal mutations that compromise fertility traits in French Lacaune dairy sheep

Background

Homozygous recessive deleterious mutations can cause embryo/fetal or neonatal lethality, or genetic defects that affect female fertility and animal welfare. In livestock populations under selection, the frequency of such lethal mutations may increase due to inbreeding, genetic drift, and/or the positive pleiotropic effects of heterozygous carriers on selected traits.

Results

By scanning the genome of 19,102 Lacaune sheep using 50 k single nucleotide polymorphism (SNP) phased genotypes and pedigree data, we identified 11 Lacaune deficient homozygous haplotypes (LDHH1 to LDHH11) showing a highly significant deficit of homozygous animals ranging from 79 to 100%. These haplotypes located on chromosomes 3, 4, 13, 17 and 18, spanned regions from 1.2 to 3.0 Mb long with a frequency of heterozygous carriers between 3.7 and 12.1%. When we compared at-risk matings (between carrier rams and daughters of carrier rams) and safe matings, seven of the 11 haplotypes were associated with a significant alteration of two fertility traits, a reduced success of artificial insemination (LDHH1, 2, 8 and 9), and/or an increased stillbirth rate (LDHH3, 6, 8, 9, and 10). The 11 haplotypes were also tested for a putative selective advantage of heterozygous carrier rams based on their daughter yield deviation for six dairy traits (milk, fat and protein yields, fat and protein contents and lactation somatic cell score). LDHH1, 3, 4, 5, 7, 9 and 11 were associated with positive effects on at least one selected dairy trait, in particular milk yield. For each haplotype, the most probable candidate genes were identified based on their roles in lethality of mouse knock-out models and in mammalian genetic disorders.

Conclusions

Based on a reverse genetic strategy, we identified at least 11 haplotypes with homozygous deficiency segregating in French Lacaune dairy sheep. This strategy represents a first tool to limit at-risk matings in the Lacaune dairy selection scheme. We assume that most of the identified LDHH are in strong linkage disequilibrium with a recessive lethal mutation that affects embryonic or juvenile survival in sheep but is yet to be identified.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00634-1.

Identification of deleterious recessive haplotypes and candidate deleterious recessive mutations in Japanese Black cattle

Intensive use of a few elite sires has increased the risk of the manifestation of deleterious recessive traits in cattle. Substantial genotyping data gathered using single-nucleotide polymorphism (SNP) arrays have identified the haplotypes with homozygous deficiency, which may compromise survival. We developed Japanese Black cattle haplotypes (JBHs) using SNP array data (4843 individuals) and identified deleterious recessive haplotypes using exome sequencing of 517 sires. We identified seven JBHs with homozygous deficiency. JBH_10 and JBH_17 were associated with the resuming of estrus after artificial insemination, indicating that these haplotypes carried deleterious mutations affecting embryonic survival. The exome data of 517 Japanese Black sires revealed that AC_000165.1:g.85341291C>G of IARS in JBH_8_2, AC_000174.1:g.74743512G>T of CDC45 in JBH_17, and a copy variation region (CNVR_27) of CLDN16 in JBH_1_1 and JBH_1_2 were the candidate mutations. A novel variant AC_000174.1:g.74743512G>T of CDC45 in JBH_17 was located in a splicing donor site at a distance of 5 bp, affecting pre-mRNA splicing. Mating between heterozygotes of JBH_17 indicated that homozygotes carrying the risk allele died around the blastocyst stage. Analysis of frequency of the CDC45 risk allele revealed that its carriers were widespread throughout the tested Japanese Black cattle population. Our approach can effectively manage the inheritance of recessive risk alleles in a breeding population.

Of rodents and ruminants: a comparison of small noncoding RNA requirements in mouse and bovine reproduction

Ruminants are major producers of meat and milk, thus managing their reproductive potential is a key element in cost-effective, safe, and efficient food production. Of particular concern, defects in male germ cells and female germ cells may lead to significantly reduced live births relative to fertilization. However, the underlying molecular drivers of these defects are unclear. Small noncoding RNAs, such as piRNAs and miRNAs, are known to be important regulators of germ-cell physiology in mouse (the best-studied mammalian model organism) and emerging evidence suggests that this is also the case in a range of ruminant species, in particular bovine. Similarities exist between mouse and bovids, especially in the case of meiotic and postmeiotic male germ cells. However, fundamental differences in small RNA abundance and metabolism between these species have been observed in the female germ cell, differences that likely have profound impacts on their physiology. Further, parentally derived small noncoding RNAs are known to influence early embryos and significant species-specific differences in germ-cell born small noncoding RNAs have been observed. These findings demonstrate the mouse to be an imperfect model for understanding germ-cell small noncoding RNA biology in ruminants and highlight the need to increase research efforts in this underappreciated aspect of animal reproduction.

Mandibulofacial Dysostosis Attributed to a Recessive Mutation of CYP26C1 in Hereford Cattle

In spring 2020, six Hereford calves presented with congenital facial deformities attributed to a condition we termed mandibulofacial dysostosis (MD). Affected calves shared hallmark features of a variably shortened and/or asymmetric lower mandible and bilateral skin tags present 2–10 cm caudal to the commissure of the lips. Pedigree analysis revealed a single common ancestor shared by the sire and dam of each affected calf. Whole-genome sequencing (WGS) of 20 animals led to the discovery of a variant (Chr26 g. 14404993T>C) in Exon 3 of CYP26C1 associated with MD. This missense mutation (p.L188P), is located in an α helix of the protein, which the identified amino acid substitution is predicted to break. The implication of this mutation was further validated through genotyping 2 additional affected calves, 760 other Herefords, and by evaluation of available WGS data from over 2500 other individuals. Only the affected individuals were homozygous for the variant and all heterozygotes had at least one pedigree tie to the suspect founder. CYP26C1 plays a vital role in tissue-specific regulation of retinoic acid (RA) during embryonic development. Dysregulation of RA can result in teratogenesis by altering the endothelin-1 signaling pathway affecting the expression of Dlx genes, critical to mandibulofacial development. We postulate that this recessive missense mutation in CYP26C1 impacts the catalytic activity of the encoded enzyme, leading to excess RA resulting in the observed MD phenotype.

Breed- and trait-specific associations define the genetic architecture of calving performance traits in cattle

Reducing the incidence of both the degree of assistance required at calving, as well as the extent of perinatal mortality (PM) has both economic and societal benefits. The existence of heritable genetic variability in both traits signifies the presence of underlying genomic variability. The objective of the present study was to locate regions of the genome, and by extension putative genes and mutations, that are likely to be underpinning the genetic variability in direct calving difficulty (DCD), maternal calving difficulty (MCD), and PM. Imputed whole-genome single-nucleotide polymorphism (SNP) data on up to 8,304 Angus (AA), 17,175 Charolais (CH), 16,794 Limousin (LM), and 18,474 Holstein-Friesian (HF) sires representing 5,866,712 calving events from descendants were used. Several putative quantitative trait loci (QTL) regions associated with calving performance both within and across dairy and beef breeds were identified, although the majority were both breed- and trait-specific. QTL surrounding and encompassing the myostatin (MSTN) gene were associated (P < 5 × 10⁻⁸) with DCD and PM in both the CH and LM populations. The well-known Q204X mutation was the fifth strongest association with DCD in the CH population and accounted for 5.09% of the genetic variance in DCD. In contrast, none of the 259 segregating variants in MSTN were associated (P > × 10⁻⁶) with DCD in the LM population but a genomic region 617 kb downstream of MSTN was associated (P < 5 × 10⁻⁸). The genetic architecture for DCD differed in the HF population relative to the CH and LM, where two QTL encompassing ZNF613 on Bos taurus autosome (BTA)18 and PLAG1 on BTA14 were identified in the former. Pleiotropic SNP associated with all three calving performance traits were also identified in the three beef breeds; 5 SNP were pleiotropic in AA, 116 in LM, and 882 in CH but no SNP was associated with more than one trait within the HF population. The majority of these pleiotropic SNP were on BTA2 surrounding MSTN and were associated with both DCD and PM. Multiple previously reported, but also novel QTL, associated with calving performance were detected in this large study. These also included QTL regions harboring SNP with the same direction of allele substitution effect for both DCD and MCD thus contributing to a more effective simultaneous selection for both traits.

Symposium review: How to implement genomic selection

Genomic selection was adopted very quickly in the 10 yr after first implementation, and breeders continue to find new uses for genomic testing. Breeding values with higher reliability earlier in life are estimated by combining DNA genotypes for many thousands of loci using existing identification, pedigree, and phenotype databases for millions of animals. Quality control for both new and previous data is greatly improved by comparing genomic and pedigree relationships to correct parent-progeny conflicts and discover many additional ancestors. Many quantitative trait loci and gene tests have been added to previous assays that used only evenly spaced, highly polymorphic markers. Imputation now combines genotypes from many assays of differing marker densities. Prediction models have gradually advanced from normal or Bayesian distributions within trait and breed to single-step, multitrait, or other more complex models, such as multibreed models that may be needed for crossbred prediction. Genomic selection was initially applied to males to predict progeny performance but is now widely applied to females or even embryos to predict their own later performance. The initial focus on additive merit has expanded to include mating programs, genomic inbreeding, and recessive alleles. Many producers now use DNA testing to decide which heifers should be inseminated with elite dairy, beef, or sex-sorted semen, which should be embryo donors or recipients, or which should be sold or kept for breeding. Because some of these decisions are expensive to delay, predictions are now provided weekly instead of every few months. Predictions from international genomic databases are often more accurate and cost-effective than those from within-country databases that were previously designed for progeny testing unless local breeds, conditions, or traits differ greatly from the larger database. Selection indexes include many new traits, often with lower heritability or requiring large initial investments to obtain phenotypes, which provide further incentive to cooperate internationally. The genomic prediction methods developed for dairy cattle are now applied widely to many animal, human, and plant populations and could be applied to many more.

Novel haplotypes responsible for prenatal death in Nordic Red and Danish Jersey cattle

Haplotypes that are common in a population but not observed as homotypes in living animals may harbor lethal alleles that compromise embryo survival. In this study, we searched for homozygous-deficient haplotypes in the genomes of 19,309 Nordic Red Dairy (RDC) and 4,291 Danish Jersey (JER) cattle genotyped using the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). For statistically significant deficient haplotypes, we evaluated the effect on nonreturn rate in at-risk matings (mating between carrier bull and daughter of carrier sire) versus not-at-risk matings (mating between noncarrier bull and daughter of noncarrier sire). Next, we analyzed whole-genome sequence variants from the 1000 Bull Genomes Project to identify putative causal variants underlying these haplotypes. In RDC, we identified 3 homozygous-deficient regions (HDR) that overlapped with known recessive lethal mutations: a 662-kb deletion on chromosome 12 in RDC [Online Mendelian Inheritance in Animals (OMIA) 001901-9913), a missense mutation in TUBD1, g.11063520T>C, in Braunvieh cattle (OMIA 001939-9913), and a 525-kb deletion on chromosome 23 in RDC (OMIA 001991-9913)]. In addition, we identified 15 novel HDR and their tag haplotypes for the underlying causative variants. The tag haplotype located between 39.2 and 40.3 Mbp on chromosome 18 had a negative effect on nonreturn rate in at-risk mating, confirming embryonic lethality. In Danish Jersey, we identified 12 novel HDR and their tag haplotypes for underlying causative variants. For 3 of these 12 tag haplotypes, insemination records of at-risk mating showed a negative effect on nonreturn rate, confirming the association with early embryonic mortality. Cattle that had both genotype and whole-genome sequence data were analyzed to detect the causative variants underlying each tag haplotype. However, none of the functional variants or deletions showed concordance with carrier status of the novel tag haplotypes. Carrier status of these detected haplotypes can be used to select bulls to reduce the frequencies of lethal alleles in the population and to avoid at-risk matings.

Genome editing approaches to augment livestock breeding programs

The prospect of genome editing offers a number of promising opportunities for livestock breeders. Firstly, these tools can be used in functional genomics to elucidate gene function, and identify causal variants underlying monogenic traits. Secondly, they can be used to precisely introduce useful genetic variation into structured livestock breeding programs. Such variation may include repair of genetic defects, the inactivation of undesired genes, and the moving of useful alleles and haplotypes between breeds in the absence of linkage drag. Editing could also be used to accelerate the rate of genetic progress by enabling the replacement of the germ cell lineage of commercial breeding animals with cells derived from genetically elite lines. In the future, editing may also provide a useful complement to evolving approaches to decrease the length of the generation interval through in vitro generation of gametes. For editing to be adopted, it will need to seamlessly integrate with livestock breeding schemes. This will likely involve introducing edits into multiple elite animals to avoid genetic bottlenecks. It will also require editing of different breeds and lines to maintain genetic diversity, and enable structured cross-breeding. This requirement is at odds with the process-based trigger and event-based regulatory approach that has been proposed for the products of genome editing by several countries. In the absence of regulatory harmony, researchers in some countries will have the ability to use genome editing in food animals, while others will not, resulting in disparate access to these tools, and ultimately the potential for global trade disruptions.

Diagnostic studies of abortion in Danish cattle 2015-2017

Background

Abortion is a major source of economic losses in cattle breeding. Abortion occurs due to a wide range of causes, but infections are the most frequently diagnosed. However, establishing an aetiological diagnosis remains challenging due to the large variety of bacteria, protozoa, viruses, and fungi that have been associated with abortion in cattle. Economic restraints limit the range of diagnostic methods available for routine diagnostics, and decomposition of the conceptus or lack of proper fetal and/or maternal samples further restrict the diagnostic success. In this study, we report recent diagnostic findings from bovine abortions in Denmark, a country that has a large dairy sector and is free from most infectious agents causing epizootic abortion in cattle. The aims of the study were: (i) to identify infectious causes of bovine abortion in Denmark, (ii) to categorise the diagnostic findings based on the level of diagnostic certainty, and (iii) to assess the diagnostic rate. Due to economic restraints, only a limited panel of routine diagnostic methods were available. Placentas and/or fetuses from mid- to late-term abortions and stillbirths (n = 162) were submitted to the Danish National Veterinary Institute between January 2015 and June 2017. The aborted materials were examined macroscopically, histologically, and by bacterial culture. Maternal blood samples were tested for bovine viral diarrhoea virus (BVDV) antibodies.

Results

The likely aetiology of the abortion was diagnosed in 52 cases, resulting in a diagnostic rate of 33%. The most common cause was protozoal infection (19%) followed by infection with Trueperella pyogenes (3%), Staphylococcus aureus (2%), and non-haemolytic Escherichia coli (2%). Lesions in fetuses with a protozoal infection were consistent with neosporosis. In many cases (38%), inflammatory changes were found in the placenta and/or fetal organs but no specific aetiology was identified. Neither infection with Brucella spp. nor maternal BVDV antibodies were detected. The majority of submitting herds (92%) were each represented by fewer than three abortion cases over the study period.

Conclusions

Protozoal infection, most likely neosporosis, was the most commonly diagnosed cause of abortion and the only one associated with potential epizootic abortion events. Despite using a reduced number of diagnostic methods in comparison to other abortion studies, the diagnostic rate of this study was within the range reported in an earlier Danish study, as well as in recent international studies. The low number of submitted cases per herd and the sparse anamnestic information provided at submission hampered conclusions on the potential epizootic character of the abortion events in question.

Fine mapping of genomic regions associated with female fertility in Nellore beef cattle based on sequence variants from segregating sires

Background

Impaired fertility in cattle limits the efficiency of livestock production systems. Unraveling the genetic architecture of fertility traits would facilitate their improvement by selection. In this study, we characterized SNP chip haplotypes at QTL blocks then used whole-genome sequencing to fine map genomic regions associated with reproduction in a population of Nellore (Bos indicus) heifers.

Methods

The dataset comprised of 1337 heifers genotyped using a GeneSeek® Genomic Profiler panel (74677 SNPs), representing the daughters from 78 sires. After performing marker quality control, 64800 SNPs were retained. Haplotypes carried by each sire at six previously identified QTL on BTAs 5, 14 and 18 for heifer pregnancy and BTAs 8, 11 and 22 for antral follicle count were constructed using findhap software. The significance of the contrasts between the effects of every two paternally-inherited haplotype alleles were used to identify sires that were heterozygous at each QTL. Whole-genome sequencing data localized to the haplotypes from six sires and 20 other ancestors were used to identify sequence variants that were concordant with the haplotype contrasts. Enrichment analyses were applied to these variants using KEGG and MeSH libraries.

Results

A total of six (BTA 5), six (BTA 14) and five (BTA 18) sires were heterozygous for heifer pregnancy QTL whereas six (BTA 8), fourteen (BTA 11), and five (BTA 22) sires were heterozygous for number of antral follicles’ QTL. Due to inadequate representation of many haplotype alleles in the sequenced animals, fine mapping analysis could only be reliably performed for the QTL on BTA 5 and 14, which had 641 and 3733 concordant candidate sequence variants, respectively. The KEGG “Circadian rhythm” and “Neurotrophin signaling pathway” were significantly associated with the genes in the QTL on BTA 5 whereas 32 MeSH terms were associated with the QTL on BTA 14. Among the concordant sequence variants, 0.2% and 0.3% were classified as missense variants for BTAs 5 and 14, respectively, highlighting the genes MTERF2, RTMB, ENSBTAG00000037306 (miRNA), ENSBTAG00000040351, PRKDC, and RGS20. The potential causal mutations found in the present study were associated with biological processes such as oocyte maturation, embryo development, placenta development and response to reproductive hormones.

Conclusions

The identification of heterozygous sires by positionally phasing SNP chip data and contrasting haplotype effects for previously detected QTL can be used for fine mapping to identify potential causal mutations and candidate genes. Genomic variants on genes MTERF2, RTBC, miRNA ENSBTAG00000037306, ENSBTAG00000040351, PRKDC, and RGS20, which are known to have influence on reproductive biological processes, were detected.

1000 Bull Genomes Project to Map Simple and Complex Genetic Traits in Cattle: Applications and Outcomes

The 1000 Bull Genomes Project is a collection of whole-genome sequences from 2,703 individuals capturing a significant proportion of the world's cattle diversity. So far, 84 million single-nucleotide polymorphisms (SNPs) and 2.5 million small insertion deletions have been identified in the collection, a very high level of genetic diversity. The project has greatly accelerated the identification of deleterious mutations for a range of genetic diseases, as well as for embryonic lethals. The rate of identification of causal mutations for complex traits has been slower, reflecting the typically small effect size of these mutations and the fact that many are likely in as-yet-unannotated regulatory regions. Both the deleterious mutations that have been identified and the mutations associated with complex trait variation have been included in low-cost SNP array designs, and these arrays are being genotyped in tens of thousands of dairy and beef cattle, enabling management of deleterious mutations in these populations as well as genomic selection.

A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves

The domestic dog serves as an excellent model to investigate the genetic basis of disease. More than 400 heritable traits analogous to human diseases have been described in dogs. To further canine medical genetics research, we established the Dog Biomedical Variant Database Consortium (DBVDC) and present a comprehensive list of functionally annotated genome variants that were identified with whole genome sequencing of 582 dogs from 126 breeds and eight wolves. The genomes used in the study have a minimum coverage of 10× and an average coverage of ~24×. In total, we identified 23 133 692 single-nucleotide variants (SNVs) and 10 048 038 short indels, including 93% undescribed variants. On average, each individual dog genome carried ∼4.1 million single-nucleotide and ~1.4 million short-indel variants with respect to the reference genome assembly. About 2% of the variants were located in coding regions of annotated genes and loci. Variant effect classification showed 247 141 SNVs and 99 562 short indels having moderate or high impact on 11 267 protein-coding genes. On average, each genome contained heterozygous loss-of-function variants in 30 potentially embryonic lethal genes and 97 genes associated with developmental disorders. More than 50 inherited disorders and traits have been unravelled using the DBVDC variant catalogue, enabling genetic testing for breeding and diagnostics. This resource of annotated variants and their corresponding genotype frequencies constitutes a highly useful tool for the identification of potential variants causative for rare inherited disorders in dogs.

Short communication: A splice site mutation in CENPU is associated with recessive embryonic lethality in Holstein cattle

A genome scan for homozygous haplotype deficiency coupled with whole-genome sequence data analysis is a very effective method to identify embryonic lethal mutations in cattle. Among other factors, the power of the approach depends on the availability of a greater amount of genotyping and sequencing data. In the present study, we analyzed the largest known panel of Illumina BovineSNP50 (Illumina Inc., San Diego, CA) genotypes, comprising 401,896 Holstein animals, and we report the mapping of a new embryonic lethal haplotype on chromosome 27, called HH7. We fine mapped the locus in a 2.0-Mb interval using an identical-by-descent approach and analyzed genome sequence data from 4 carrier and 143 noncarrier Holstein bulls to identify the causative mutation. We detected a strong candidate variant in the gene encoding centromere protein U (CENPU), a centromere component essential for proper chromosome segregation during mitosis. The mutant allele is a deletion of 4 nucleotides located at position +3 to +6 bp after the splicing donor site of exon 11. Cross-species nucleotide alignment revealed that the nucleotide at position +3 is entirely conserved among vertebrates, suggesting that it plays an important role in the regulation of CENPU splicing. For verification, we genotyped the candidate variant in 232,775 Holstein individuals and did not observe any homozygotes, whereas 16 were expected (Poisson P-value = 1.1 × 10^-7; allele frequency = 0.8%). In addition, genotyping of 250,602 animals from 19 additional breeds revealed that the mutant allele is restricted to animals of Holstein descent. Finally, we estimated the effect of the candidate variant on 2 fertility traits in at-risk mating (i.e., between carrier bulls and daughters of carrier bulls) versus non-risk mating. In agreement with a recessive lethal inheritance pattern, we observed a marked reduction in both conception rate and 56-d nonreturn rate in heifers and cows. The effect on 56-d nonreturn rate suggests that a substantial proportion of homozygous mutants die before 35 d after insemination, which is consistent with the early embryonic death previously reported in CENPU^-/- mouse embryos. In conclusion, we demonstrate that with more than 400,000 genotypes, we can map very rare recessive lethal mutations segregating at a frequency below 1% in the population. We recommend performing new analyses regularly as data are accumulating.

Technical note: Development and application of KASP assays for rapid screening of 8 genetic defects in Holstein cattle

Specific DNA mutations underlying several genetic defects associated with embryo loss or reduced calf survivability have been identified in dairy cattle, and a convenient and cost-effective platform is required for their routine screening. We developed Kompetitive allele-specific PCR (KASP) assays for discrimination of the wild-type alleles from the associated defective alleles at each of 8 common genetic defects in Holstein cattle, involving 5 SNP [HH1, HH3, HH4, bovine leukocyte adhesion deficiency (BLAD), and complex vertebral malformation (CVM)] and 3 insertion or deletion mutations [HH5, haplotype for cholesterol deficiency (HCD), and brachyspina (BS)]. A total of 390 cows from a Chinese Holstein herd were genotyped and the carriers identified at 7 of these 8 loci (except HH4), with the highest carrier frequencies found for CVM (10.5%) and HH1 (10.0%), followed by HH3 (2.6%), BS (2.1%), HCD (1.3%), HH5 (0.8%), and BLAD (0.5%). Surprisingly, 102 cows (26.2%) carried at least 1 of the 7 defective alleles. Our results demonstrate that these KASP assays are simple, rapid, and reliable for the detection of multiple genetic defects. The high carrier frequency of these genetic defects indicates an urgent need for routine molecular testing to eliminate the deleterious alleles from Chinese Holstein cattle.

Management of lethal recessive alleles in beef cattle through the use of mate selection software

BACKGROUND

Recessive loss-of-function (LOF) alleles at genes which are essential for life, can result in early embryonic mortality. Cattle producers can use the LOF carrier status of individual animals to make selection and mate allocation decisions.

METHODS

Two beef cattle breeding strategies i.e. (1) selection against LOF carriers as parents and (2) simultaneous selection and mate allocation to avoid the occurrence of homozygous offspring in three scenarios, which differed in number and frequency of LOF alleles were evaluated using the mate selection program, MateSel. Scenarios included (a) seven loci with high-frequency LOF alleles, (b) 76 loci with low-frequency LOF alleles, and (c) 50 loci with random high- and low-frequency LOF alleles. In addition, any savings resulting from the information obtained by varying the percentage (0-100%) of the herd genotyped, together with segregation analysis to cover ungenotyped animals, were calculated to determine (1) which percentage optimized net profit for a fixed cost of genotyping ($30/test), and (2) the breakeven cost for genotyping.

RESULTS

With full knowledge of the LOF alleles carried by selection candidates, the most profitable breeding strategy was always simultaneous selection and mate allocation to avoid homozygous affected offspring (aa) as compared to indiscriminate selection against carrier parents (Aa). The breakeven value of genotyping depended on the number of loci modeled, the LOF allele frequencies, and the mating/selection strategies used. Genotyping was most valuable when it was used to avoid otherwise high levels of embryonic mortalities. As the number of essential loci with LOF alleles increased, especially when some were present at relatively high minor allele frequencies, embryonic losses increased, and profit was maximized by genotyping 10 to 20% of a herd and using that information to reduce these losses.

CONCLUSIONS

Genotyping 100% of the herd was never the most profitable outcome in any scenario; however, genotyping some proportion of the herd, together with segregation analysis to cover ungenotyped animals, maximized overall profit in scenarios with large numbers of loci with LOF alleles. As more LOF alleles are identified, such a mate selection software will likely be required to optimally select and allocate matings to balance the rate of genetic gain, embryonic losses, and inbreeding.

The cholesterol-deficiency associated mutation in APOB segregates at low frequency in Chinese Holstein cattle

Cholesterol deficiency (CD) is a recessive genetic defect in Holstein cattle, which affects the rearing success of calves. We investigated whether CD has been spread into the People’s Republic of China. A total of 138 Holstein bulls and 90 cows were screened for CD using polymerase chain reaction amplification and then seven bulls (5.07%) and one cow (1.11%) were identified as CD carriers. Pedigree analyses revealed that all seven CD carrier bulls were traced back to MAUGHLIN STORM and CD may have spread into the People’s Republic of China through his daughter, Braedale Baler Twine. Our results demonstrated that the CD defect exists in the Chinese Holstein population.

Estimating the effect of the deleterious recessive haplotypes AH1 and AH2 on reproduction performance of Ayrshire cattle

The effects of 2 deleterious recessive haplotypes on reproduction performance of Ayrshire cattle, Ayrshire Haplotype 1 (AH1) and Ayrshire Haplotype 2 (AH2), were investigated in Canadian Ayrshire cattle. We calculated their phenotypic effects on stillbirth (SB) rate and 56-d nonreturn rate (NRR) by estimating the interaction of service sire carrier status with maternal grandsire carrier status using the official Canadian evaluation models for those 2 traits. The interaction term included 9 subclasses for the 3 possible statuses of each bull: haplotype carrier, noncarrier, or not genotyped. For AH1, 394 carriers and 1,433 noncarriers were available, whereas 313 carriers and 1,543 noncarriers were available for the AH2 haplotype. The number of matings considered for SB was 34,312 for heifers (first parity) and 115,935 for cows (later parities). For NRR, 49,479 matings for heifers and 160,528 for cows were used to estimate the haplotype effects. We observed a negative effect of AH1 on SB rates, which was 2.0% higher for matings of AH1-carrier sires to dams that had an AH1-carrier sire; this effect was found for both heifers and cows. However, AH1 had small, generally nonsignificant effects on NRR. The AH2 haplotype had a substantial negative effect on NRR, with 5.1% more heifers and 4.0% more cows returning to service, but the effects on SB rates were inconsistent and mostly small effects. Our results validate the harmful effects of AH1 and AH2 on reproduction traits in the Canadian Ayrshire population. This information will be of great interest for the dairy industry, allowing producers to make mating decisions that would reduce reproductive losses.

Analysis of a large dataset reveals haplotypes carrying putatively recessive lethal and semi-lethal alleles with pleiotropic effects on economically important traits in beef cattle

BACKGROUND

In livestock, deleterious recessive alleles can result in reduced economic performance of homozygous individuals in multiple ways, e.g. early embryonic death, death soon after birth, or semi-lethality with incomplete penetrance causing reduced viability. While death is an easy phenotype to score, reduced viability is not as easy to identify. However, it can sometimes be observed as reduced conception rates, longer calving intervals, or lower survival for live born animals.

METHODS

In this paper, we searched for haplotypes that carry putatively recessive lethal or semi-lethal alleles in 132,725 genotyped Irish beef cattle from five breeds: Aberdeen Angus, Charolais, Hereford, Limousin, and Simmental. We phased the genotypes in sliding windows along the genome and used five tests to identify haplotypes with absence of or reduced homozygosity. Then, we associated the identified haplotypes with 44,351 insemination records that indicated early embryonic death, and postnatal survival records. Finally, we assessed haplotype pleiotropy by estimating substitution effects on estimates of breeding value for 15 economically important traits in beef production.

RESULTS

We found support for one haplotype that carries a putatively recessive lethal (chromosome 16 in Simmental) and two haplotypes that carry semi-lethal alleles (chromosome 14 in Aberdeen Angus and chromosome 19 in Charolais), with population frequencies of 8.8, 15.2, and 14.4%, respectively. These three haplotypes showed pleiotropic effects on economically important traits for beef production. Their allele substitution effects are €2.30, €3.42, and €1.47 for the terminal index and €1.03, - €3.11, and - €0.88 for the replacement index, where the standard deviations for the terminal index are €22.52, €18.65, and €22.70 and for the replacement index they are €31.35, €29.82, and €35.79. We identified ZFAT as the candidate gene for semi-lethality in Aberdeen Angus, several candidate genes for the lethal Simmental haplotype, and no candidate genes for the semi-lethal Charolais haplotype.

CONCLUSIONS

We analysed genotype, reproduction, survival, and production data to detect haplotypes that carry putatively recessive lethal or semi-lethal alleles in Irish beef cattle and identified one lethal and two semi-lethal haplotypes, which have pleiotropic effects on economically important traits in beef production.

Determination of quantitative trait nucleotides by concordance analysis between quantitative trait loci and marker genotypes of US Holsteins

Experimental designs that exploit family information can provide substantial predictive power in quantitative trait nucleotide discovery projects. Concordance between quantitative trait locus genotype as determined by the a posteriori granddaughter design and marker genotype was determined for 30 trait-by-chromosomal segment effects segregating in the US Holstein population with probabilities of <10^-20 to accept the null hypotheses of no segregating gene affecting the trait within the chromosomal segment. Genotypes for 83 grandsires and 17,217 sons were determined by either complete sequence or imputation for 3,148,506 polymorphisms across the entire genome; 471 Holstein bulls had a complete genome sequence, including 64 of the grandsires. Complete concordance was obtained only for stature on chromosome 14 and daughter pregnancy rate on chromosome 18. For each quantitative trait locus, effects of the 30 polymorphisms with highest concordance scores for the analyzed trait were computed by stepwise regression for predicted transmitting abilities of 26,750 bulls with progeny test and imputed genotypes. Effects for stature on chromosome 11, daughter pregnancy rate on chromosome 18, and protein percentage on chromosome 20 met 3 criteria: complete or almost complete concordance, nominal significance of the polymorphism effect after correction for all other polymorphisms, and marker coefficient of determination >40% of total multiple-regression coefficient of determination for the 30 polymorphisms with highest concordance. An intronic variant marker on chromosome 5 at 93,945,738 bp explained 7% of variance for fat percentage and 74% of total multiple-marker regression variance but was concordant for only 24 of 30 families. The missense polymorphism Phe279Tyr in GHR at 31,909,478 bp on chromosome 20 was confirmed as the causative mutation for fat and protein concentration. For effect on fat percentage on chromosome 14, 12 additional missense polymorphisms were found that had almost complete concordance with the suggested causative polymorphism (missense mutation Ala232Glu in DGAT1). The only polymorphism found likely to improve predictive power for genomic evaluation of dairy cattle was on chromosome 15; that polymorphism had a frequency of 0.45 for the allele with economically positive effects on all production traits.

A survey of functional genomic variation in domesticated chickens

BACKGROUND

Deleterious genetic variation can increase in frequency as a result of mutations, genetic drift, and genetic hitchhiking. Although individual effects are often small, the cumulative effect of deleterious genetic variation can impact population fitness substantially. In this study, we examined the genome of commercial purebred chicken lines for deleterious and functional variations, combining genotype and whole-genome sequence data.

RESULTS

We analysed over 22,000 animals that were genotyped on a 60 K SNP chip from four purebred lines (two white egg and two brown egg layer lines) and two crossbred lines. We identified 79 haplotypes that showed a significant deficit in homozygous carriers. This deficit was assumed to stem from haplotypes that potentially harbour lethal recessive variations. To identify potentially deleterious mutations, a catalogue of over 10 million variants was derived from 250 whole-genome sequenced animals from three purebred white-egg layer lines. Out of 4219 putative deleterious variants, 152 mutations were identified that likely induce embryonic lethality in the homozygous state. Inferred deleterious variation showed evidence of purifying selection and deleterious alleles were generally overrepresented in regions of low recombination. Finally, we found evidence that mutations, which were inferred to be evolutionally intolerant, likely have positive effects in commercial chicken populations.

CONCLUSIONS

We present a comprehensive genomic perspective on deleterious and functional genetic variation in egg layer breeding lines, which are under intensive selection and characterized by a small effective population size. We show that deleterious variation is subject to purifying selection and that there is a positive relationship between recombination rate and purging efficiency. In addition, multiple putative functional coding variants were discovered in selective sweep regions, which are likely under positive selection. Together, this study provides a unique molecular perspective on functional and deleterious variation in commercial egg-laying chickens, which can enhance current genomic breeding practices to lower the frequency of undesirable variants in the population.

Invited review: A perspective on the future of genomic selection in dairy cattle

Genomic evaluation has been successfully implemented in the United States, Canada, Great Britain, Ireland, New Zealand, Australia, France, the Netherlands, Germany, and the Scandinavian countries. Adoption of this technology in the major dairy producing countries has led to significant changes in the worldwide dairy industry. Gradual elimination of the progeny test system has led to a reduction in the number of sires with daughter records and fewer genetic ties between years. As genotyping costs decrease, the number of cows genotyped will continue to increase, and these records will become the basic data used to compute genomic evaluations, most likely via application of "single-step" methodologies. Although genomic selection has been successful in increasing rates of genetic gain, we still know very little about the genetic architecture of quantitative variation. Apparently, a very large number of genes affect nearly all economic traits, in accordance with the infinitesimal model for quantitative traits. Less emphasis in selection goals will be placed on milk production traits, and more on health, reproduction, and efficiency traits and on environmentally friendly production with reduced waste and gas emission. Genetic variance for economic traits is maintained by the increase in frequency of rare alleles, new mutations, and changes in selection goals and management. Thus, it is unlikely that a selection plateau will be reached in the near future.

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.