A large population study to assess the magnitude of prenatal programming in dairy cattle.

The performance of an adult dairy cow may be influenced by events that occur before her birth. The present study investigated potential effects of 2 prenatal groups of factors, Assisted Reproductive Techniques (ART) and maternal characteristics (e.g., dam parity), on offspring performance during their first lactation, in populations of 2 dairy cow breeds: French Holstein and Montbéliarde. The different ART studied included the type of semen (conventional or X-sorted) used for Artificial Insemination (AI) and the technology of conception used (AI, embryo transfer, or in vitro fertilization). Three maternal characteristics were considered: (1) the dam age at first calving, (2) dam parity number, and (3) indicators of dam udder health during gestation (somatic cell score and events of clinical mastitis). First, we investigated whether heifer survival from 3d to 18 mo old was associated with any of the prenatal factors considered. We then estimated the associations of these prenatal factors with 8 traits of commercial interest: (1) stature, (2-4) milk, fat, and protein yields, (5) somatic cell score, (6) clinical mastitis, and (7-8) heifer and cow conception rate, all measured on genotyped cows. Linear models were used for this study with the prenatal factors as covariates in the model, and for the 8 traits, phenotypes were adjusted for their corresponding genomic estimated breeding value. The results indicated that the survival rate of heifers born from embryo transfer was significantly higher than that of heifers born from AI (probably due to preferential management practices), while the other prenatal factors did not explain differences in heifer survival. Among the Montbéliarde cows born from AI, those born from X-sorted semen showed a lightly but significantly lower milk yield than those born without X-sorting of the semen (-52 kg of milk in the first lactation). Among the Holstein cows, those born from embryo transfer presented significantly lower milk performance than cows born from AI. Regarding the maternal characteristics, none or very weak associations were found between the dam age at first calving and the offspring performance in both breeds. Dam parity, on the other hand, was associated with offspring performance for milk, fat, and protein yield in both breeds, however not in the same direction. In the Holstein breed, an increase in dam parity was favorable for offspring performance for milk, fat, and protein yield, whereas in the Montbéliarde breed, an increase in dam parity was associated with lower milk and protein yield and no association was found for fat yield. The udder health of the dam during gestation was not or only weakly associated with the traits studied in the offspring. Although some significant associations were identified due to the large sample size, the effects were modest, typically less than 1% of the phenotypic mean, and were not consistently observed across the 2 breeds.

Detailed analysis of mortality rates in the female progeny of 1,001 Holstein bulls allows the discovery of new dominant genetic defects.

Reducing juvenile mortality in cattle is important for both economic and animal welfare reasons. Previous studies have revealed a large variability in mortality rates between breeds and sire progeny groups, with some extreme cases due to dominant mutations causing various syndromes among the descendants of mosaic bulls. The purpose of this study was to monitor sire-family calf mortality within the French and Walloon Holstein populations, and to use this information to detect genetic defects that might have been overlooked by lack of specific symptoms. In a population of heifers born from 1,001 bulls between 2017 and 2020, the average sire-family mortality rates were of 11.8% from birth to 1 year of age and of 4.2, 2.9, 3.1, and 3.2% for the perinatal, postnatal, preweaning, and postweaning subperiods, respectively. After outlining the 5 worst bulls per category, we paid particular attention to the bulls Mo and Pa, because they were half-brothers. Using a battery of approaches, including necropsies, karyotyping, genetic mapping, and whole-genome sequencing, we described 2 new independent genetic defects in their progeny and their molecular etiology. Mo was found to carry a de novo reciprocal translocation between chromosomes BTA26 and BTA29, leading to increased embryonic and juvenile mortality because of aneuploidy. Clinical examination of 2 calves that were monosomic for a large proportion of BTA29, including an orthologous segment deleted in human Jacobsen syndrome, revealed symptoms shared between species. In contrast, Pa was found to be mosaic for a dominant de novo nonsense mutation of GATA 6 binding protein (GATA6), causing severe cardiac malformations. In conclusion, our results highlight the power of monitoring juvenile mortality to identify dominant genetic defects due to de novo mutation events.

A missense mutation in the FZD7 gene is associated with dilution of the red areas of the coat in Montbéliarde cattle.

Recently, a new genetically autosomal recessive color phenotype emerged in the red pied bovine Montbéliarde breed. It is characterized by a dilution of the red areas of the coat and was denominated 'milca'. A genome-wide homozygosity scan of 106 cases followed by haplotype analysis revealed a candidate region within BTA2 between positions 89.95 and 91.63 Mb. Analysis of whole-genome sequence data generated from milca animals identified a strong candidate variant within the coding region of the Frizzled-7 gene (FZD7). This gene encodes for a G-protein coupled receptor for Wnt signaling proteins. The variant induces a glycine to alanine substitution in the second extracellular loop, p.(Gly414Ala). Cross-species amino acid alignments revealed that this glycine is conserved among orthologs and most paralogs, suggesting that it plays an important role in FZD function. In addition, genotyping data revealed that the mutant allele is restricted to the Montbéliarde breed, at a 3.7% frequency. All homozygous cows for the mutant allele exhibited the milca phenotype whereas all heterozygotes had no coat color defects. In conclusion, this study strongly suggests that, in cattle, a mutation of FZD7 alone is sufficient to cause a coat color phenotype without any strong other adverse effect.

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.

Inferring the population structure of the Maghreb sheep breeds using a medium-density SNP chip.

North Africa has a great diversity of indigenous sheep breeds whose origin is linked to its environmental characteristics and to certain historical events that took place in the region. To date, few genome-wide studies have been conducted to investigate the population structure of North African indigenous sheep. The objective of the present study was to provide a detailed assessment of the genetic structure and admixture patterns of six Maghreb sheep populations using the Illumina 50K Ovine BeadChip and comparisons with 22 global populations of sheep and mouflon. Regardless of the method of analysis used, patterns of multiple hybridization events were observed within all North African populations, leading to a heterogeneous genetic architecture that varies according to the breed. The Barbarine population showed the lowest genetic heterogeneity and major southwest Asian ancestry, providing additional support to the Asian origin of the North African fat-tailed sheep. All other breeds presented substantial Merino introgression ranging from 15% for D'man to 31% for Black Thibar. We highlighted several signals of ancestral introgression between North African and southern European sheep. In addition, we identified two opposite gradients of ancestry, southwest Asian and central European, occurring between North Africa and central Europe. Our results provide further evidence of the weak global population structure of sheep resulting from high levels of gene flow among breeds occurring worldwide. At the regional level, signs of recent admixture among North African populations, resulting in a change of the original genomic architecture of minority breeds, were also detected.

A missense mutation (p.Tyr452Cys) in the CAD gene compromises reproductive success in French Normande cattle.

We scanned the genome of 77,815 Normande cattle with different Illumina SNP chips (Illumina Inc., San Diego, CA) to map recessive embryonic lethal mutations using homozygous haplotype deficiency. We detected 2 novel haplotypes on chromosomes 11 and 24 but did not confirm 6 previously reported haplotypes. The one on chromosome 11 showed a marked reduction in conception rates and moderate decrease in nonreturn rate in at-risk versus control mating, supporting late embryonic mortality. After fine mapping and analyzing whole-genome sequences, we prioritized a missense mutation in CAD (g.72399397T>C; p.Tyr452Cys)-a gene encoding a protein (carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase) essential for de novo pyrimidine biosynthesis-as a candidate causal variant. This transition mutation replaces a tyrosine residue, which is perfectly conserved among living organisms, with a cysteine residue in the carbamoyl-phosphate synthetase 2 domain of the protein. A single animal was confirmed to be homozygous for the mutation based on Sanger sequencing. However, large-scale genotyping of the candidate variant with the Illumina EuroG10k BeadChip revealed an absence of live homozygotes in a panel of 33,323 Normande animals and an absence of carriers in 348,593 animals from 19 other cattle breeds. These results support recessive embryonic lethality with nearly complete penetrance, as was previously reported in CAD mutants in several eukaryote species. The only homozygous cow had extremely poor udder conformation, suggesting a potential role of CAD in udder development, but no effect was detected when comparing daughter yield deviations of 250 heterozygous bulls with that of 2,912 homozygotes for the ancestral allele. Together, our results showed the importance of large-scale screening for homozygous haplotype deficiency with hundreds of thousands of animals, validating results with an independent data set, and considering unexpected live homozygotes, to avoid both false-positive and false-negative discoveries. These discoveries will be used primarily in mating decisions to avoid at-risk mating. In addition, we recommend including CAD in the breeding objectives of Normande cattle.

Use of meta-analyses and joint analyses to select variants in whole genome sequences for genomic evaluation: An application in milk production of French dairy cattle breeds.

As a result of the 1000 Bull Genome Project, it has become possible to impute millions of variants, with many of these potentially causative for traits of interest, for thousands of animals that have been genotyped with medium-density chips. This enormous source of data opens up very interesting possibilities for the inclusion of these variants in genomic evaluations. However, for computational reasons, it is not possible to include all variants in genomic evaluation procedures. One potential approach could be to select the most relevant variants based on the results of genome-wide association studies (GWAS); however, the identification of causative mutations is still difficult with this method, partly because of weak imputation accuracy for rare variants. To address this problem, this study assesses the ability of different approaches based on multi-breed GWAS (joint and meta-analyses) to identify single-nucleotide polymorphisms (SNP) for use in genomic evaluation in the 3 main French dairy cattle breeds. A total of 6,262 Holstein bulls, 2,434 Montbéliarde bulls, and 2,175 Normande bulls with daughter yield deviations for 5 milk production traits were imputed for 27 million variants. Within-breed and joint (including all 3 breeds) GWAS were performed and 3 models of meta-analysis were tested: fixed effect, random effect, and Z-score. Comparison of the results of within- and multi-breed GWAS showed that most of the quantitative trait loci identified using within-breed approaches were also found with multi-breed methods. However, the most significant variants identified in each region differed depending on the method used. To determine which approach highlighted the most predictive SNP for each trait, we used a marker-assisted best unbiased linear prediction model to evaluate lists of SNP generated by the different GWAS methods; each list contained between 25 and 2,000 candidate variants per trait, which were identified using a single within- or multi-breed GWAS approach. Among all the multi-breed methods tested in this study, variant selection based on meta-analysis (fixed effect) resulted in the most-accurate genomic evaluation (+1 to +3 points compared with other multi-breed approaches). However, the accuracies of genomic evaluation were always better when variants were selected using the results of within-breed GWAS. As has generally been found in studies of quantitative trait loci, these results suggest that part of the genetic variance of milk production traits is breed specific in Holstein, Montbéliarde, and Normande cattle.

KIRREL is differentially expressed in adipose tissue from 'fertil+' and 'fertil-' cows: in vitro role in ovary?

We have previously shown that dairy cows carrying the 'fertil-' haplotype for one quantitative trait locus affecting female fertility located on the bovine chromosome three (QTL-F-Fert-BTA3) have a significantly lower conception rate and body weight after calving than cows carrying the 'fertil+' haplotype. Here, we compared by Tiling Array the expression of genes included in the QTL-F-Fert-BTA3 in 'fertil+' and 'fertil-' adipose tissue one week after calving when plasma non-esterified fatty acid concentrations were greater in 'fertil-' animals. We observed that thirty-one genes were overexpressed whereas twelve were under-expressed in 'fertil+' as compared to 'fertil-' cows (P < 0.05). By quantitative PCR and immunoblot we confirmed that adipose tissue KIRREL mRNA and protein were significantly greater expressed in 'fertil+' than in 'fertil-'. KIRREL mRNA is abundant in bovine kidney, adipose tissue, pituitary, and ovary and detectable in hypothalamus and mammary gland. Its expression (mRNA and protein) is greater in kidney of 'fertil+' than 'fertil-' cows (P < 0.05). KIRREL (mRNA and protein) is also present in the different ovarian cells with a greater expression in granulosa cells of 'fertil+' than 'fertil-' cows. In cultured granulosa cells, recombinant KIRREL halved steroid secretion in basal state (P < 0.05). It also decreased cell proliferation (P < 0.05) and in vitro oocyte maturation (P < 0.05). These results were associated to a rapid increase in MAPK1/3 and MAPK14 phosphorylation in granulosa cells and to a decrease in MAPK1/3 phosphorylation in oocyte. Thus, KIRREL could be a potential metabolic messenger linking body composition and fertility.

Whole-genome scan to detect quantitative trait loci associated with milk protein composition in 3 French dairy cattle breeds.

In the context of the PhénoFinLait project, a genome-wide analysis was performed to detect quantitative trait loci (QTL) that affect milk protein composition estimated using mid-infrared spectrometry in the Montbéliarde (MO), Normande (NO), and Holstein (HO) French dairy cattle breeds. The 6 main milk proteins (α-lactalbumin, ß-lactoglobulin, and αS1-, αS2-, ß-, and κ-caseins) expressed as grams per 100g of milk (% of milk) or as grams per 100g of protein (% of protein) were estimated in 848,068 test-day milk samples from 156,660 cows. Genotyping was performed for 2,773 MO, 2,673 NO, and 2,208 HO cows using the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA). Individual test-day records were adjusted for environmental effects and then averaged per cow to define the phenotypes analyzed. Quantitative trait loci detection was performed within each breed using a linkage disequilibrium and linkage analysis approach. A total of 39 genomic regions distributed on 20 of the 29 Bos taurus autosomes (BTA) were significantly associated with milk protein composition at a genome-wide level of significance in at least 1 of the 3 breeds. The 9 most significant QTL were located on BTA2 (133 Mbp), BTA6 (38, 47, and 87 Mbp), BTA11 (103 Mbp), BTA14 (1.8 Mbp), BTA20 (32 and 58 Mbp), and BTA29 (8 Mbp). The BTA6 (87 Mbp), BTA11, and BTA20 (58 Mbp) QTL were found in all 3 breeds, and they had highly significant effects on κ-casein, ß-lactoglobulin, and α-lactalbumin, expressed as a percentage of protein, respectively. Each of these QTL explained between 13% (BTA14) and 51% (BTA11) of the genetic variance of the trait. Many other QTL regions were also identified in at least one breed. They were located on 14 additional chromosomes (1, 3, 4, 5, 7, 15, 17, 19, 21, 22, 24, 25, 26, and 27), and they explained 2 to 8% of the genetic variance of 1 or more protein composition traits. Concordance analyses, performed between QTL status and sequence-derived polymorphisms from 13 bulls, revealed previously known causal polymorphisms in LGB (BTA11) and GHR (BTA20 at 32 Mbp) and excluded some other previously described mutations. These results constitute a first step in identifying causal mutations and using routinely collected mid-infrared predictions in future genomic selection programs to improve bovine milk protein composition.

Identification of the putative ancestral allele of bovine single-nucleotide polymorphisms.

Identifying the action of natural selection from patterns of standing genetic variation has long been of interest to the population genetic community. Thanks to the availability of large single-nucleotide polymorphism (SNP) data sets for many species and of high-throughput SNP genotyping methods, whole-genomic surveys to detect selective sweeps are now possible. Knowing the ancestral allele increases the power to detect selection. We present here a comparative genomic approach to determine the putative ancestral allele of bovine SNPs deposited in public databases. We analysed 19,551,488 SNPs and identified the putative ancestral allele for 14,339,107 SNPs. Our predicted ancestral alleles were in agreement with ancestral alleles detected by genotyping outgroup species for 97% SNPs from the BovineSNP50 BeadChip. This comparison indicates that our comparative genomic-based approach to identify putative ancestral alleles is reliable.

Genome-wide linkage disequilibrium in the Blonde d'Aquitaine cattle breed.

We present here the first genome-wide characterization of linkage disequilibrium (LD) in the French Blonde d'Aquitaine (BLA) breed, a well-muscled breed renowned for producing high-yielding beef carcasses. To assess the pattern and extent of LD, we used a sample of 30 unrelated bulls and 36 923 single nucleotide polymorphisms (SNPs) covering all cattle autosomes. The squared correlation of the alleles at two loci (r(2) ) was used as a measure of LD. The analysis of adjacent marker pairs revealed that the level of LD decreases rapidly with physical distance between SNPs. Overall mean r(2) was 0.205 (±0.262). Strong LD (r(2)  > 0.8) and useful LD (measured as r(2 ) > 0.2) were observed within genomic regions of up to 720 and 724 kb, respectively. We analysed the genetic structure of the BLA population and found stratification. The observed genetic sub-structuring is consistent with the known recent demographic history that occurred during BLA breed formation. Our results indicate that LD mapping of phenotypic traits in the BLA population is feasible; however, because of this sub-structuring, special care is needed to reduce the likelihood of false-positive associations between marker loci and traits of interest.

Comparative mapping and genomic annotation of the bovine oncosuppressor gene WWOX.

WWOX (WW domain-containing oxidoreductase) is the gene mapping at FRA16D HSA16q23.1, the second most active common fragile site in the human genome. In this study we characterized at a detailed molecular level WWOX in the bovine genome. First, we sequenced cDNA from various tissues and obtained evidence in support of a 9-exon structure for the gene, similar to the human gene. Then, we recovered BACs using exon tags and annotated the gene to a >1-Mb genomic region of BTA18 using the Btau 4.0 genome assembly as a reference, thus resolving an issue related to exon 9, which is not included in the genomic annotation of the gene in the Entrez database. Finally, BACs spanning WWOX were used as FISH probes to obtain comparative mapping of the gene in Bos taurus, Bubalus bubalis, Ovis aries and Capra hircus to BTA18q12.1, BBU18q13, OAR14q12.1 and CHI18q12.1, respectively. Our data show that the chromosomal location of WWOX is conserved between man and 4 major domesticated species. Moreover, the annotation of the bovine gene also suggests a highly conserved genomic arrangement, including number and size of introns.

Genome screen for asthma and related phenotypes in the French EGEA study.

A genome-wide search was conducted in 107 nuclear families with at least two siblings with asthma, as part of the French EGEA study. A two-stage analysis strategy was applied to the 107 families divided into two independent subsets of 46 and 61 families, where all regions detected in the first set of families were tested for replication in the second set. In addition, all regions reported by published genome scans in different populations were examined in the total sample. A total of 254 markers were typed in the first set of families and 70% of them in the second set. Linkage was investigated by model-free methods for asthma and four asthma-related phenotypes: bronchial responsiveness (BR), skin test response, total immunoglobulin E (IgE) levels, and eosinophil count. The two-stage analysis led to the detection of three regions: 11p13 for IgE, 12q24 for eosinophils, and 17q12-21 for asthma and skin tests. Among the regions reported by published genome screens, seven were found in the 107 French EGEA families: three being already detected by the two-stage analysis, 11p13 (p = 0.005), 12q24 (p = 0.0008), and 17q12-21 (p = 0.001), and four additional ones, 1p31 (p = 0.005) for asthma, 11q13 (p = 0.006) for IgE, 13q31 (p = 0.001) for eosinophils, and 19q13 (p = 0.02) for BR.

Cloning and sequence analysis of a Marek's disease virus origin binding protein (OBP) reveals strict conservation of structural motifs among OBPs of divergent alphaherpesviruses.

Marek's disease virus (MDV) is a highly cell-associated avian herpesvirus. In its natural host, MDV induces Marek's disease (MD), a lethal condition characterized by malignant lymphoma of T cells. Although symptoms of MD may be prevented by vaccination, no practical pharmacological method of control has been widely accepted. Viral replication represents a point at which pharmacological control of herpesvirus infection may be most successful. However, this requires detailed knowledge of viral replication proteins. Studies in HSV-1 DNA replication implicate the UL9 protein as a key initiator of replication. For example, binding of UL9 to HSV-1 origins is a prerequisite for assembly of additional replication proteins. In this study, a protein, whose apparent molecular size is similar to that of HSV-1 UL9, was identified in extracts of MDV infected cells by western blot analysis with anti-HSV-1 UL9 antibody. A putative MDV UL9 gene was subsequently identified through sequencing of MDV genome fragments (BamHI G and C). Extended DNA sequence analysis revealed an open reading frame (ORF) which could encode a protein homologous to HSV-1 UL9. The MDV UL9 ORF encodes 841 amino acids, producing a sequence 49% identical to HSV-1 UL9 and 46% identical to VZV gene 51 product (VZV UL9). MDV UL9 shares numerous structural motifs with HSV-1 and VZV UL9 proteins, including six conserved N-terminal helicase motifs, an N-terminal leucine zipper motif, a C-terminal pseudo-leucine zipper sequence, and a putative helix-turn-helix structure.