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.

Identification and characterisation of mitochondrial sequences integrated into the ovine nuclear genome.

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, generating nuclear mitochondrial DNA sequences (NUMTs). Here, we analysed, for the first time, NUMTs in the ovine genome. We obtained 760 alignment matches covering 513.8 kbp of the sheep nuclear genome. After a merging step, we identified 390 NUMT regions with a total length of ~720 kbp, representing 0.02% of the nuclear genome. We discovered copies of all mitochondrial regions and found that most NUMT regions are intergenic or intronic. Ovine NUMTs are mostly not transcribed. However, we identified within some of the NUMTs, potential new genes encoding nuclear humanin isoforms. To rule out the possibility that the identified NUMTs could be artifacts of the Oar Rambouillet v1.0 genome assembly, we validated experimentally nine NUMT regions by PCR amplification. As we found several NUMT regions showing high similarity to the mitochondrial genome that potentially could pose a risk to ovine DNA mitochondrial studies, special care must be taken for the selection of primers for PCR amplification of mitochondrial DNA sequences.

Detection of selection signatures in Limousin cattle using whole-genome resequencing.

Limousin, a renowned beef breed originating from central France, has been selectively bred over the last 100 years to improve economically important traits. We used whole-genome sequencing data from 10 unrelated Limousin bull calves to detect polymorphisms and identify regions under selection. A total of 13 943 766 variants were identified. Moreover, 311 852 bi-allelic SNPs and 92 229 indels located on autosomes were fixed for the alternative allele in all sequenced animals, including the previously reported missense deleterious F94L mutation in MSTN. We performed a whole-genome screen to discover genomic regions with excess homozygosity, using the pooled heterozygosity score and identified 171 different candidate selective sweeps. In total, 68 candidate genes were found in only 57 of these regions, indicating that a large fraction of the genome under selection might lie in non-coding regions and suggesting that a majority of adaptive mutations might be regulatory in nature. Many QTL were found within candidate selective sweep regions, including QTL associated with shear force or carcass weight. Among the putative selective sweeps, we located genes (MSTN, NCKAP5, RUNX2) that potentially contribute to important phenotypes in Limousin. Several candidate regions and genes under selection were also found in previous genome-wide selection scans performed in Limousin. In addition, we were able to pinpoint candidate causative regulatory polymorphisms in GRIK3 and RUNX2 that might have been under selection. Our results will contribute to improved understanding of the mechanisms and targets of artificial selection and will facilitate the interpretation of GWASs performed in Limousin.