Differences in breeds and sire lines suggest the presence of a genetic component for heifer livability (HLIV). Genomic evaluation for this trait can increase profitability and improve animal health and welfare. Evaluations for HLIV were examined from 3,362,499 calf data records from heifers of all breeds born from 2009 to 2016. Data were obtained from the national cooperator database maintained by the Council on Dairy Cattle Breeding (https://www.uscdcb.com/). The total number of deaths reported was 134,753 (4.01%), which included herds with death loss between 1.5 and 25.5%. Age at death was evaluated and ranged from >2 d of age until the heifer left the herd, with a maximum of 18 mo of age. Records were not included until 3 yr after the birthdate so that live status of contemporaries could be confirmed by a calving date for those animals. Deaths observed until 2 d after birth were considered to be a stillbirth rather than a failure of HLIV. The scale used for analysis of HLIV was 0 (died) or 100 (live), and the heritability estimate was 0.7% based on sire model with restricted maximum likelihood estimation. Genomic predicted transmitting abilities for Holstein ranged from -1.6% to +1.6% with a standard deviation of 0.5%, and genomic predicted transmitting abilities for Jersey ranged from -0.5% to +0.5% with a standard deviation of 0.2%. The mean overall death loss was about 4%. Reliabilities of genomic predictions for young animals averaged 46% for Holsteins and 30% for Jerseys, and corresponding traditional parent average reliabilities averaged 16% and 12%, respectively. Correlations of HLIV were 0.44 with productive life, 0.18 to 0.22 with yield traits, and 0.29 with early first calving on proven Holstein bulls. The HLIV trait had a favorable genetic trend in recent years, likely because of the indirect selection associated with the correlated traits. The trait HLIV should receive 1% of emphasis on the Lifetime Net Merit index, resulting in economic progress worth $50,000/yr. By encouraging more comprehensive recording on calf mortality, the reliabilities of genetic predictions could increase significantly.

A recessive haplotype resulting in elevated calf mortality but with apparent incomplete penetrance was previously linked to the end of chromosome 16 (78.7-80.7 Mbp). Genotype analysis of 5.6 million Holsteins indicated that the haplotype was common and traced back to 1952, with a key ancestor born in 1984 (HOUSA1964484, Southwind) identified from chip genotypes as homozygous for the suspect haplotype. Sequence data from Southwind (an affected calf) and the sire of the affected calf were scanned for candidate mutations. A missense mutation with a deleterious projected impact at 79,613,592 bp was homozygous in the affected calf and heterozygous in the calf's sire and Southwind. Sequence data available from the Cooperative Dairy DNA Repository for 299 other Holsteins indicated a 97% concordance with the haplotype and an 89% call rate. The exon amino acid sequence appears to be broadly conserved in the CACNA1S gene, and mutations in humans and mice can cause phenotypes of temporary or permanent paralysis analogous to those in calves with the haplotype causing muscle weakness (HMW). Improved methods for using pedigree to track new mutations within existing haplotypes were developed and applied to the haplotypes for both muscle weakness and Holstein cholesterol deficiency (HCD). For HCD, concordance of the gene test with its haplotype status was greatly improved. For both defects, haplotype status was matched to heifer livability records for 558,000 calves. For HMW, only 46 heifers with livability records were homozygous and traced only to Southwind on both sides. Of those, 52% died before 18 mo at an average age of 1.7 ± 1.6 mo, but that death rate may be underestimated if only healthier calves were genotyped. The death rate was 2.4% for noncarriers. Different reporting methods or dominance effects may be needed to include HMW and other partially lethal effects in selection and mating. Direct tests are needed for new mutations within existing common haplotypes because tracking can be difficult even with accurate pedigrees when the original haplotype has a high frequency.

Background: India harbors the world’s largest cattle population, encompassing over 50 distinct Bos indicus breeds. This rich genetic diversity underscores the inadequacy of a single reference genome to fully capture the genomic landscape of Indian cattle. To comprehensively characterize the genomic variation within Bos indicus and, specifically, dairy breeds, we aim to identify non-reference sequences and construct a comprehensive pangenome. Results: Five representative genomes of prominent dairy breeds, including Gir, Kankrej, Tharparkar, Sahiwal, and Red Sindhi, were sequenced using 10X Genomics ‘linked-read’ technology. Assemblies generated from these linked-reads ranged from 2.70 Gb to 2.77 Gb, comparable to the Bos indicus Brahman reference genome. A pangenome of Bos indicus cattle was constructed by comparing the newly assembled genomes with the reference using alignment and graph-based methods, revealing 8 Mb and 17.7 Mb of novel sequence respectively. A confident set of 6,844 Non-reference Unique Insertions (NUIs) spanning 7.57 Mb was identified through both methods, representing the pangenome of Indian Bos indicus breeds. Comparative analysis with previously published pangenomes unveiled 2.8 Mb (37%) commonality with the Chinese indicine pangenome and only 1% commonality with the Bos taurus pangenome. Among these, 2,312 NUIs encompassing ~ 2 Mb, were commonly found in 98 samples of the 5 breeds and designated as Bos indicus Common Insertions (BICIs) in the population. Furthermore, 926 BICIs were identified within 682 protein-coding genes, 54 long non-coding RNAs (lncRNA), and 18 pseudogenes. These protein-coding genes were enriched for functions such as chemical synaptic transmission, cell junction organization, cell-cell adhesion, and cell morphogenesis. The protein-coding genes were found in various prominent quantitative trait locus (QTL) regions, suggesting potential roles of BICIs in traits related to milk production, reproduction, exterior, health, meat, and carcass. Notably, 63.21% of the bases within the BICIs call set contained interspersed repeats, predominantly Long Interspersed Nuclear Elements (LINEs). Additionally, 70.28% of BICIs are shared with other domesticated and wild species, highlighting their evolutionary significance. Conclusions: This is the first report unveiling a robust set of NUIs defining the pangenome of Bos indicus breeds of India. The analyses contribute valuable insights into the genomic landscape of desi cattle breeds.

Presence-absence variations (PAVs) are important structural variations, wherein a genomic segment containing one or more genes is present in some individuals but absent in others. While PAVs have been extensively studied in plants, research in cattle remains limited. This study identified PAVs in 173 Holstein bulls using whole-genome sequencing data and assessed their associations with 46 economically important traits. Out of 28,772 cattle genes (from the longest transcripts), a total of 26,979 (93.77%) core genes were identified (present in all individuals), while variable genes included 928 softcore (present in 95–99% of individuals), 494 shell (present in 5–94%), and 371 cloud genes (present in <5%). Cloud genes were enriched in functions associated with hormonal and antimicrobial activities, while shell genes were enriched in immune functions. PAV-based genome-wide association studies identified associations between gene PAVs and 16 traits including milk, fat, and protein yields, as well as traits related to health and reproduction. Associations were found on multiple chromosomes, illustrating important associations on cattle chromosomes 7 and 15, involving olfactory receptor and immune-related genes, respectively. By examining the PAVs at the population level, the results of this research provided crucial insights into the genetic structures underlying the complex traits of Holstein cattle.

Summary: Late-term abortions >152 days and <251 days of gestation that terminate a lactation or initiate a new lactation have long been recorded by Dairy Herd Improvement (DHI). For 22.7 million DHI lactations, the average recorded incidence of late-term abortions across all years (2001–2018) was 1.2%. However, the 1.3% incidence of late-term abortions reported in 2012 has declined to 75%). Thus, PTA for late-term abortions should not be used as a separate fertility trait because national evaluations for current fertility traits already account for those economic losses.

Background: Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms. Results: We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull’s family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats. Conclusion: Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.