Searching for homozygous haplotype deficiency in Manech Tête Rousse dairy sheep revealed a nonsense variant in the MMUT gene affecting newborn lamb viability

BACKGROUND

Recessive deleterious variants are known to segregate in livestock populations, as in humans, and some may be lethal in the homozygous state.

RESULTS

We used phased 50 k single nucleotide polymorphism (SNP) genotypes and pedigree data to scan the genome of 6845 Manech Tête Rousse dairy sheep to search for deficiency in homozygous haplotypes (DHH). Five Manech Tête Rousse deficient homozygous haplotypes (MTRDHH1 to 5) were identified, with a homozygous deficiency ranging from 84 to 100%. These haplotypes are located on Ovis aries chromosome (OAR)1 (MTRDHH2 and 3), OAR10 (MTRDHH4), OAR13 (MTRDHH5), and OAR20 (MTRDHH1), and have carrier frequencies ranging from 7.8 to 16.6%. When comparing at-risk matings between DHH carriers to safe matings between non-carriers, two DHH (MTRDHH1 and 2) were linked with decreased insemination success and/or increased stillbirth incidence. We investigated the MTRDHH1 haplotype, which substantially increased stillbirth rate, and identified a single nucleotide variant (SNV) inducing a premature stop codon (p.Gln409*) in the methylmalonyl-CoA mutase (MMUT) gene by using a whole-genome sequencing approach. We generated homozygous lambs for the MMUT mutation by at-risk mating between heterozygous carriers, and most of them died within the first 24 h after birth without any obvious clinical symptoms. Reverse transcriptase-qPCR and western blotting on post-mortem liver and kidney biological samples showed a decreased expression of MMUT mRNA in the liver and absence of a full-length MMUT protein in the mutant homozygous lambs.

CONCLUSIONS

We identified five homozygous deficient haplotypes that are likely to harbor five independent deleterious recessive variants in sheep. One of these was detected in the MMUT gene, which is associated with lamb lethality in the homozygous state. A specific management of these haplotypes/variants in the MTR dairy sheep selection program would help enhance the overall fertility and lamb survival.

Integrative Multiomics Analysis of the Heat Stress Response of Enterococcus faecium

A continuous heat-adaptation test was conducted for one Enterococcus faecium (E. faecium) strain wild-type (WT) RS047 to obtain a high-temperature-resistant strain. After domestication, the strain was screened with a significantly higher ability of heat resistance. which is named RS047-wl. Then a multi-omics analysis of transcriptomics and metabolomics was used to analyze the mechanism of the heat resistance of the mutant. A total of 98 differentially expressed genes (DEGs) and 115 differential metabolites covering multiple metabolic processes were detected in the mutant, which indicated that the tolerance of heat resistance was regulated by multiple mechanisms. The changes in AgrB, AgrC, and AgrA gene expressions were involved in quorum-sensing (QS) system pathways, which regulate biofilm formation. Second, highly soluble osmotic substances such as putrescine, spermidine, glycine betaine (GB), and trehalose-6P were accumulated for the membrane transport system. Third, organic acids metabolism and purine metabolism were down-regulated. The findings can provide target genes for subsequent genetic modification of E. faecium, and provide indications for screening heat-resistant bacteria, so as to improve the heat-resistant ability of E. faecium for production.

Disorders of purine biosynthesis metabolism

Purines are essential molecules that are components of vital biomolecules, such as nucleic acids, coenzymes, signaling molecules, as well as energy transfer molecules. The de novo biosynthesis pathway starts from phosphoribosylpyrophosphate (PRPP) and eventually leads to the synthesis of inosine monophosphate (IMP) by means of 10 sequential steps catalyzed by six different enzymes, three of which are bi-or tri-functional in nature. IMP is then converted into guanosine monophosphate (GMP) or adenosine monophosphate (AMP), which are further phosphorylated into nucleoside di- or tri-phosphates, such as GDP, GTP, ADP and ATP. This review provides an overview of inborn errors of metabolism pertaining to purine synthesis in humans, including either phosphoribosylpyrophosphate synthetase (PRS) overactivity or deficiency, as well as adenylosuccinate lyase (ADSL), 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/IMP cyclohydrolase (ATIC), phosphoribosylaminoimidazole succinocarboxamide synthetase (PAICS), and adenylosuccinate synthetase (ADSS) deficiencies. ITPase deficiency is being described as well. The clinical spectrum of these disorders is broad, including neurological impairment, such as psychomotor retardation, epilepsy, hypotonia, or microcephaly; sensory involvement, such as deafness and visual disturbances; multiple malformations, as well as muscle presentations or consequences of hyperuricemia, such as gouty arthritis or kidney stones. Clinical signs are often nonspecific and, thus, overlooked. It is to be hoped that this is likely to be gradually overcome by using sensitive biochemical investigations and next-generation sequencing technologies.

Mining massive genomic data of two Swiss Braunvieh cattle populations reveals six novel candidate variants that impair reproductive success

Background

This study was carried out on the two Braunvieh populations reared in Switzerland, the dairy Brown Swiss (BS) and the dual-purpose Original Braunvieh (OB). We performed a genome-wide analysis of array data of trios (sire, dam, and offspring) from the routine genomic selection to identify candidate regions showing missing homozygosity and phenotypic associations with five fertility, ten birth, and nine growth-related traits. In addition, genome-wide single SNP regression studies based on 114,890 single nucleotide polymorphisms (SNPs) for each of the two populations were performed. Furthermore, whole-genome sequencing data of 430 cattle including 70 putative haplotype carriers were mined to identify potential candidate variants that were validated by genotyping the current population using a custom array.

Results

Using a trio-based approach, we identified 38 haplotype regions for BS and five for OB that segregated at low to moderate frequencies. For the BS population, we confirmed two known haplotypes, BH1 and BH2. Twenty-four variants that potentially explained the missing homozygosity and associated traits were detected, in addition to the previously reported TUBD1:p.His210Arg variant associated with BH2. For example, for BS we identified a stop-gain variant (p.Arg57*) in the MRPL55 gene in the haplotype region on chromosome 7. This region is associated with the ‘interval between first and last insemination’ trait in our data, and the MRPL55 gene is known to be associated with early pregnancy loss in mice. In addition, we discuss candidate missense variants in the CPT1C, MARS2, and ACSL5 genes for haplotypes mapped in BS. In OB, we highlight a haplotype region on chromosome 19, which is potentially caused by a frameshift variant (p.Lys828fs) in the LIG3 gene, which is reported to be associated with early embryonic lethality in mice. Furthermore, we propose another potential causal missense variant in the TUBGCP5 gene for a haplotype mapped in OB.

Conclusions

We describe, for the first time, several haplotype regions that segregate at low to moderate frequencies and provide evidence of causality by trait associations in the two populations of Swiss Braunvieh. We propose a list of six protein-changing variants as potentially causing missing homozygosity. These variants need to be functionally validated and incorporated in the breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00686-3.

Reverse Genetic Screen for Deleterious Recessive Variants in the Local Simmental Cattle Population of Switzerland

We herein report the result of a large-scale reverse genetic screen in the Swiss Simmental population, a local dual-purpose cattle breed. We aimed to detect possible recessively inherited variants affecting protein-coding genes, as such deleterious variants can impair fertility and rearing success significantly. We used 115,000 phased SNP data of almost 10 thousand cattle with pedigree data. This revealed evidence for 11 genomic regions of 1.17 Mb on average, with haplotypes (SH1 to SH11) showing a significant depletion in homozygosity and an allele frequency between 3.2 and 10.6%. For the proposed haplotypes, it was unfortunately not possible to evaluate associations with fertility traits as no corresponding data were available. For each haplotype region, possible candidate genes were listed based on their known function in development and disease. Subsequent mining of single-nucleotide variants and short indels in the genomes of 23 sequenced haplotype carriers allowed us to identify three perfectly linked candidate causative protein-changing variants: a SH5-related DIS3:p.Ile678fs loss-of-function variant, a SH8-related CYP2B6:p.Ile313Asn missense variant, and a SH9-related NUBPL:p.Ser143Tyr missense variant. None of these variants occurred in homozygous state in any of more than 5200 sequenced cattle of various breeds. Selection against these alleles in order to reduce reproductive failure and animal loss is recommended.

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.

Comparative transcriptome analysis reveals key insights into male sterility in Salvia miltiorrhiza Bunge

Background

Large-scale heterosis breeding depends upon stable, inherited male sterility lines. We accidentally discovered a male sterility line (SW-S) in the F₁progeny of a Salvia miltiorrhiza Bunge from Shandong, China (purple flowers) crossed with a S. miltiorrhiza f. alba from Sichuan, China (white flowers). We sought to provide insights into the pollen development for male sterility in S. miltiorrhiza.

Methods

The phenotypic and cytological features of the SW-S and fertile control SW-F were observed using scanning electron microscopy and paraffin sections to identify the key stage of male sterility. Transcriptome profiles were recorded for anthers at the tetrad stage of SW-S and SW-F using Illumina RNA-Seq.

Results

The paraffin sections showed that sterility mainly occurred at the tetrad stage of microspore development, during which the tapetum cells in the anther compartment completely fell off and gradually degraded in the sterile line. There was little-to-no callose deposited around the microspore cells. The tetrad microspore was shriveled and had abnormal morphology. Therefore, anthers at the tetrad stage of SW-S and fertile control SW-F were selected for comparative transcriptome analysis. In total, 266,722,270 clean reads were obtained from SW-S and SW-F, which contained 36,534 genes. There were 2,571 differentially expressed genes (DEGs) in SW-S and SW-F, of which 63.5% were downregulated. Gene Ontology (GO) enrichment analysis indicated that the differentially expressed genes were enriched in 56 functional groups (GO terms); of these, all DEGs involved in microgametogenesis and developmental maturation were downregulated in SW-S. These results were confirmed by quantitative RT-PCR. The two GO terms contained 18 DEGs, among which eight DEGs (namely: GPAT3, RHF1A, phosphatidylinositol, PFAS, MYB96, MYB78, Cals5, and LAT52) were related to gamete development. There were 10 DEGs related to development and maturation, among which three genes were directly related to pollen development (namely: ACT3, RPK2, and DRP1C). Therefore, we believe that these genes are directly or indirectly involved in the pollen abortion of SW-S. Our study provides insight into key genes related to sterility traits in S. miltiorrhiza, and the results can be further exploited in functional and mechanism studies.

Genomic Approaches Reveal Pleiotropic Effects in Crossbred Beef Cattle

Carcass and meat quality are two important attributes for the beef industry because they drive profitability and consumer demand. These traits are of even greater importance in crossbred cattle used in subtropical and tropical regions for their superior adaptability because they tend to underperform compared to their purebred counterparts. Many of these traits are challenging and expensive to measure and unavailable until late in life or after the animal is harvested, hence unrealistic to improve through traditional phenotypic selection, but perfect candidates for genomic selection. Before genomic selection can be implemented in crossbred populations, it is important to explore if pleiotropic effects exist between carcass and meat quality traits. Therefore, the objective of this study was to identify genomic regions with pleiotropic effects on carcass and meat quality traits in a multibreed Angus-Brahman population that included purebred and crossbred animals. Data included phenotypes for 10 carcass and meat quality traits from 2,384 steers, of which 1,038 were genotyped with the GGP Bovine F-250. Single-trait genome-wide association studies were first used to investigate the relevance of direct additive genetic effects on each carcass, sensory and visual meat quality traits. A second analysis for each trait included all other phenotypes as covariates to correct for direct causal effects from identified genomic regions with pure direct effects on the trait under analysis. Five genomic windows on chromosomes BTA5, BTA7, BTA18, and BTA29 explained more than 1% of additive genetic variance of two or more traits. Moreover, three suggestive pleiotropic regions were identified on BTA10 and BTA19. The 317 genes uncovered in pleiotropic regions included anchoring and cytoskeletal proteins, key players in cell growth, muscle development, lipid metabolism and fat deposition, and important factors in muscle proteolysis. A functional analysis of these genes revealed GO terms directly related to carcass quality, meat quality, and tenderness in beef cattle, including calcium-related processes, cell signaling, and modulation of cell-cell adhesion. These results contribute with novel information about the complex genetic architecture and pleiotropic effects of carcass and meat quality traits in crossbred beef cattle.

Improved dairy cattle mating plans at herd level using genomic information

From 2012 to 2018, 223 180 Montbéliarde females were genotyped in France and the number of newly genotyped females increased at a rate of about 33% each year. With female genotyping information, farmers have access to the genomic estimated breeding values of the females in their herd and to their carrier status for genetic defects or major genes segregating in the breed. This information, combined with genomic coancestry, can be used when planning matings in order to maximize the expected on-farm profit of future female offspring. We compared different mating allocation approaches for their capacity to maximize the expected genetic gain while limiting expected progeny inbreeding and the probability to conceive an offspring homozygous for a lethal recessive allele. Three mate allocation strategies (random mating (RAND), sequential mating (gSEQ€) and linear programing mating (gLP€)) were compared on 160 actual Montbéliarde herds using male and female genomic information. Then, we assessed the benefit of using female genomic information by comparing matings planned using only female pedigree information with the equivalent strategy using genomic information. We measured the benefit of adding genomic expected inbreeding and risk of conception of an offspring homozygous for a lethal recessive allele to Net merit in mating plans. The influence of three constraints was tested: by relaxing the constraint on availability of a particular semen type (sexed or conventional) for bulls, by adding an upper limit of 8.5% coancestry between mate pairs or by using a more stringent maximum use of a bull in a herd (5% vs 10%). The use of genomic information instead of pedigree information improved the mate allocation method in terms of progeny expected genetic merit, genetic diversity and risk to conceive an offspring homozygous for a lethal recessive allele. Optimizing mate allocation using linear programming and constraining coancestry to a maximum of 8.5% per mate pair reduced the average coancestry with a small impact on expected Net Merit. In summary, for male and female selection pathways, using genomic information is more efficient than using pedigree information to maximize genetic gain while constraining the expected inbreeding of the progeny and the risk to conceive an offspring homozygous for a lethal recessive allele. This study also underlines the key role of semen type (sexed vs conventional) and the associated constraints on the mate allocation algorithm to maximize genetic gain while maintaining genetic diversity and limiting the risk to conceive an offspring homozygous for a lethal recessive allele.

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: Lethal mutations in Vorderwald cattle through Montbéliarde incrossings

Vorderwald cattle are a dual-purpose cattle breed with high migrant contributions from Montbéliarde bulls in the recent past. Through the wide use of Montbéliarde bulls, undesirable alleles were also disseminated into the Vorderwald population. Haplotypes on bovine chromosome 19 (MH1) and 29 (MH2), supposed to harbor lethal mutations, were identified in Montbéliarde cattle. A study in French Montbéliarde cattle identified the PFAS:g.28511199C>T (rs455876205) variant as the most likely MH1 embryonic lethal mutation. The objective of the present study was to determine whether the PFAS:g.28511199C>T variant was introduced into Vorderwald cattle through Montbéliarde bulls and disseminated in this population. The present study expands on previous work on the deleterious SLC37A2 variant (ss2019324563) of the MH2 locus. Herein, we traced the ss2019324563 variant back to the Montbéliarde bull, which was the most likely source for this deleterious mutation in Vorderwald cattle. We genotyped 354 Vorderwald cattle for the PFAS variant, resulting in 41 heterozygous individuals and a T allele frequency of 0.058. An aborted fetus homozygous mutant for SLC37A2 from our previous study on the MH2 locus in Vorderwald cattle was wild type for the PFAS variant. Both lethal mutations were segregating independently of each other, and we found no indications of joint occurrence in a larger number of animals. Neither SLC37A2 nor PFAS double heterozygous mutants were lethal. The earliest animal with a heterozygous PFAS genotype was 1 of 5 migrant Montbéliarde bulls, and this bull was the most likely origin of the deleterious PFAS allele in Vorderwald cattle. All Vorderwald cattle under study born before introgression of this Montbéliarde bull were homozygous wild type. In addition, all 41 heterozygous Vorderwald cattle had genetic contributions from this Montbéliarde bull, whereas in 74 Vorderwald cattle without genes from Montbéliarde bulls, the PFAS T allele was not observed. In a sample of actual German Fleckvieh the PFAS T allele could be found at a very low frequency. Our study demonstrated the introgression of lethal variants through Montbéliarde bulls into a traditional cattle breed highly adapted to harsh local conditions. These findings underline the need to screen bulls for lethal mutations before their wide use in breeding, particularly in breeds with a focus on preservation of their genetic uniqueness.

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.

Genomic Regions Associated With Gestation Length Detected Using Whole-Genome Sequence Data Differ Between Dairy and Beef Cattle

While many association studies exist that have attempted to relate genomic markers to phenotypic performance in cattle, very few have considered gestation length as a phenotype, and of those that did, none used whole genome sequence data from multiple breeds. The objective of the present study was therefore to relate imputed whole genome sequence data to estimated breeding values for gestation length using 22,566 sires (representing 2,262,706 progeny) of multiple breeds [Angus (AA), Charolais (CH), Holstein-Friesian (HF), and Limousin (LM)]. The associations were undertaken within breed using linear mixed models that accounted for genomic relatedness among sires; a separate association analysis was undertaken with all breeds analysed together but with breed included as a fixed effect in the model. Furthermore, the genome was divided into 500 kb segments and whether or not segments harboured a single nucleotide polymorphism (SNP) with a P ≤ 1 × 10^-4 common to different combinations of breeds was determined. Putative quantitative trait loci (QTL) regions associated with gestation length were detected in all breeds; significant associations with gestation length were only detected in the HF population and in the across-breed analysis of all 22,566 sires. Twenty-five SNPs were significantly associated (P ≤ 5 × 10^-8) with gestation length in the HF population. Of the 25 significant SNPs, 18 were located within three QTLs on Bos taurus autosome number (BTA) 18, six were in two QTL on BTA19, and one was located within a QTL on BTA7. The strongest association was rs381577268, a downstream variant of ZNF613 located within a QTL spanning from 58.06 to 58.19 Mb on BTA18; it accounted for 1.37% of the genetic variance in gestation length. Overall there were 11 HF animals within the edited dataset that were homozygous for the T allele at rs381577268 and these had a 3.3 day longer (P < 0.0001) estimated breeding value (EBV) for gestation length than the heterozygous animals and a 4.7 day longer (P < 0.0001) EBV for gestation length than the homozygous CC animals. The majority of the 500 kb windows harboring a SNP with a P ≤ 1 × 10^-4 were unique to a single breed and no window was shared among all four breeds for gestation length, suggesting any QTLs identified are breed-specific associations.

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.

A Meta-Analysis Including Pre-selected Sequence Variants Associated With Seven Traits in Three French Dairy Cattle Populations

A within-breed genome-wide association study (GWAS) is useful when identifying the QTL that segregates in a breed. However, an across-breed meta-analysis can be used to increase the power of identification and precise localization of QTL that segregate in multiple breeds. Precise localization will allow including QTL information from other breeds in genomic prediction due to the persistence of the linkage phase between the causal variant and the marker. This study aimed to identify and confirm QTL detected in within-breed GWAS through a meta-analysis in three French dairy cattle breeds. A set of sequence variants selected based on their functional annotations were imputed into 50 k genotypes for 46,732 Holstein, 20,096 Montbeliarde, and 11,944 Normande cows to identify QTL for milk production, the success rate at insemination of cows (fertility) and stature. We conducted within-breed GWAS followed by across-breed meta-analysis using a weighted Z-scores model on the GWAS summary data (i.e., P-values, effect direction, and sample size). After Bonferroni correction, the GWAS result identified 21,956 significantly associated SNP (P _FWER < 0.05), while meta-analysis result identified 9,604 significant SNP (P _FWER < 0.05) associated with the phenotypes. The meta-analysis identified 36 QTL for milk yield, 48 QTL for fat yield and percentage, 29 QTL for protein yield and percentage, 13 QTL for fertility, and 16 QTL for stature. Some of these QTL were not significant in the within-breed GWAS. Some previously identified causal variants were confirmed, e.g., BTA14:1802265 (fat percentage, P = 1.5 × 10^-760; protein percentage, P = 7.61 × 10^-348) both mapping the DGAT1-K232A mutation and BTA14:25006125 (P = 8.58 × 10^-140) mapping PLAG1 gene was confirmed for stature in Montbeliarde. New QTL lead SNP shared between breeds included the intronic variant rs109205829 (NFIB gene), and the intergenic variant rs41592357 (1.38 Mb upstream of the CNTN6 gene and 0.65 Mb downstream of the CNTN4 gene). Rs110425867 (ZFAT gene) was the top variant associated with fertility, and new QTL lead SNP included rs109483390 (0.1 Mb upstream of the TNFAIP3 gene and 0.07 Mb downstream of PERP gene), and rs42412333 (0.45 Mb downstream of the RPL10L gene). An across-breed meta-analysis had greater power to detect QTL as opposed to a within breed GWAS. The QTL detected here can be incorporated in routine genomic predictions.