An endogenous retrovirus insertion disrupting bovine ALKBH8 causes a failure-to-thrive syndrome with immunodeficiency associated with juvenile mortality in Brown Swiss cattle
This study demonstrates that exploring population-scale genomic data and mining thousands of life-history records, followed by veterinary follow-up evaluations and molecular genetic analyses, provides an effective strategy for identifying cryptic recessive disorders that shorten the lifespan of cattle.
Abstract
The Brown Swiss (BS) cattle breed is one of the major Swiss dairy breeds. Intensive selection and the widespread use of few elite sires in artificial insemination have increased inbreeding and the occurrence of deleterious recessive alleles in the homozygous state. Analyzing life trajectories in large, genotyped cohorts can identify hidden recessive disorders that are difficult to detect using traditional case-control association testing. Long-read DNA sequencing enables precise detection of causal alleles, including structural variants. This study aimed to (1) identify cryptic recessive loci affecting rearing performance in Swiss BS cattle, (2) evaluate their impact on survival, (3) characterize the associated phenotype, (4) identify the causal variant using long-read whole-genome sequencing, and (5) assess its functional impact. Using Homozygous Haplotype Enrichment/Depletion (HHED) mapping, we identified a risk haplotype (BH39) on chromosome 15 spanning from 16,276,819 bp to 16,446,984 bp that was associated with increased juvenile mortality within the first 180 days of life when present in the homozygous state. The BH39 occurred at a frequency of approximately 4.5% in Swiss BS cattle and 5.3% in German and Austrian BS cattle, and homozygous carriers exhibited a significantly reduced first-year survival rate. Five females homozygous for BH39 underwent clinical examination. They all showed recurrent respiratory disease, impaired growth, poor body condition, rough hair coat, and brown-discolored teeth. Pathological examination revealed bronchopneumonia and eosinophilic enteritis. Clinicopathological findings indicated failure to thrive and immunodeficiency. Long-read WGS of two BH39 homozygous calves revealed a private homozygous coding variant that was in high linkage disequilibrium with BH39. The identified structural variant was an insertion of a large transposable element (10.4 kb ERVK[2-1-LTR]) into the third exon of ALKBH8 (NM_001080341.2 c.267_268indel). Full-length RNA sequencing of cerebellum and liver from a homozygous calf revealed that the endogenous retrovirus (ERV) insertion introduces a cryptic transcription termination signal, truncating ALKBH8 mRNA. This study demonstrates that exploring population-scale genomic data and mining thousands of life-history records, followed by veterinary follow-up evaluations and molecular genetic analyses, provides an effective strategy for identifying cryptic recessive disorders that shorten the lifespan of cattle. The findings provide strong evidence that the ERV insertion into the coding sequence of ALKBH8 represents a loss-of-function variant that causes a previously undescribed recessive disorder that results in increased rearing loss. Interpretive summary We identified a recessive disorder in Brown Swiss cattle that causes retarded growth, recurrent infections, immunodeficiency, and increased mortality during the first year of life. Using population-scale genomic data, clinical investigations, and long-read sequencing, we linked the disorder to an exonic transposable element insertion disrupting ALKBH8. The identification of the causal variant now enables direct genetic testing and the implementation of genome-based mating strategies to avoid carrier-by-carrier matings and, consequently, prevent the birth of affected homozygous offspring. We demonstrate the utility of integrating large-scale breeding records, veterinary phenotyping, and advanced genomics to identify hidden defects affecting livestock health and productivity.
Genomic signatures of selection can reveal loci underlying adaptation and disease resistance in livestock populations, but such analyses in water buffalo (Bubalus bubalis) have historically been constrained by the absence of a chromosome-level, species-native reference genome for SNP array data. We re-analyzed genotype data from 85 Nili-Ravi buffalo (Axiom Buffalo Genotyping 90K array, originally positioned using bovine (Bos taurus, UMD3.1) proxy coordinates, by performing a full coordinate liftover to the buffalo-native UOA_WB_1 assembly using an independently published SNP remapping resource. Following quality control (51,209 markers retained), haplotype phasing, and genome-wide integrated haplotype score (iHS) and Wright’s Fst (case/control) selection scans, we evaluated 14 classical bovine-tuberculosis (bTB) candidate genes and identified six additional genes with putative immune function through an unbiased genome-wide screen. None of the 14 classical candidates (including SLC11A1, the Toll-like receptors, and IFNG) reached genome-wide significance in either scan. In contrast, six novel loci TNFSF18, IL2RB, TNFRSF19, IRF2, IL15, and CD28 showed significant iHS or Fst signals, four of which (TNFSF18, IL2RB, IL15, CD28) converge functionally on T-cell costimulation and cytokine receptor signaling (KEGG pathways map04660 and map04060, Bos taurus proxy annotation). Using extended haplotype homozygosity (EHH) decay, haplotype furcation structure, and per-marker haplotype counts as three independent lines of corroborating evidence, we classified these six genes into confidence tiers: TNFSF18 and IL2RB showed the strongest, most balanced support, while CD28 and IL15 signals were driven by very few haplotypes (3 and 5 of 30, respectively) and should be interpreted cautiously pending replication. These findings suggest that adaptive, cell-mediated immune signaling rather than the innate/macrophage-centred mechanisms emphasized by existing bTB candidate gene panels may be a more productive avenue for future selection studies in Nili-Ravi buffalo, while underscoring the value of buffalo-native coordinate systems for accurate genomic inference in this species.
Atiq Ahmad, Abu Bakar, Sheikh Muhammad Laeeque et al.· bioRxiv· 0 citations
Bovine alphaherpesvirus 1 (BoAHV1; family Orthoherpesviridae, taxon species Varicellovirus bovinealpha1)-a major cause of reproductive losses in cattle-establishes lifelong latency with periodic reactivation. Although vaccination reduces disease severity and abortion risk, BoAHV1 continues to circulate in cattle populations. Here, we report the whole-genome sequencing (WGS) of a BoAHV1 isolate recovered from fetal tissues associated with a last-trimester abortion in a vaccinated heifer herd. BoAHV1 was confirmed by quantitative PCR (qPCR) and virus isolation in Madin-Darby bovine kidney cells. The cell culture-isolated virus was subsequently sequenced using an Oxford Nanopore platform. Comparative genomic analysis using reference strain MH791339 identified 73 nucleotide variants (~0.12% of the genome), mostly located in the unique short (US) region. Multiple high-impact frameshift mutations were detected in several viral genes, including independent frameshifts in the infective cell protein 4 gene (bICP4) and moderate-impact variants in US3, a gene involved in immune modulation. Overall, genomic features were consistent with a wild-type BoAHV1.1 field strain and not compatible with known modified-live virus vaccine genomes. Our findings highlight the extent of genomic variability of BoAHV1 circulating in vaccinated cattle populations and demonstrate the value of WGS for resolving diagnostic questions related to abortion events and vaccine strain exclusion.
Karen Salas-Briceno, Lester E Buckner, Manoj K. Sekhwal et al.· Journal of Veterinary Diagno...· 0 citations
NHLRC2-associated FINCA disease is an ultra-rare autosomal recessive multisystem disorder caused by biallelic pathogenic variants in NHLRC2. Its mutational spectrum and genotype–phenotype correlations remain incompletely defined, and the contribution of non-coding variants is poorly understood. Here, we report a male infant with a severe FINCA-like phenotype, including early-onset hemolytic anemia, pulmonary involvement, neurodevelopmental impairment, growth failure, recurrent infections, and fatal progression at 8.5 months. Whole-genome sequencing identified a compound heterozygous NHLRC2 genotype comprising the previously reported pathogenic missense variant c.442G>T (p.Asp148Tyr) and a novel deep intronic variant, c.331+6863A>G. Segregation analysis confirmed inheritance from different parents. Integrated genomic and splicing analysis predicted that c.331+6863A>G creates a strong cryptic donor splice site and supports pseudoexon inclusion. Reconstruction of the predicted aberrant transcript indicated premature termination and potential susceptibility to nonsense-mediated mRNA decay. To our knowledge, this is the first reported deep intronic NHLRC2 variant predicted to activate pseudoexon inclusion. Although experimental validation was unavailable, convergent clinical, segregation, population, and computational evidence supports c.331+6863A>G as the most plausible second disease-associated allele. This case expands the genomic spectrum of NHLRC2-associated FINCA disease and highlights the diagnostic value of phenotype-driven whole-genome sequencing.
A. Rozhkova, Anton Esibov, A. Borkovskaia et al.· International Journal of Mol...· 0 citations
Recommendations are formulated for building a genetic health management system in beef cattle breeding, based on mandatory DNA testing, inbreeding control, accounting for pleiotropic effects in selection, and the integration of advanced biotechnologies.
N. Bezborodova, O. Sokolova, O. S. Zaitseva et al.· International Journal of Vet...· 0 citations
Hirschsprung disease (HSCR) is a complex developmental disorder of the enteric nervous system, primarily driven by regulatory variants within enhancer elements of the RET gene. To investigate how these variants lead to aganglionosis, we developed a humanized mouse model by inserting an intact 77kb human RET genomic locus into the Rosa26 safe-harbor locus. Utilizing “big DNA” synthetic biology and Bxb1-mediated recombination, we integrated the complete human locus including all exons, introns, and upstream regulatory elements which we validated via nanopore and short-read sequencing. Functional analysis confirmed in vivo human RET expression; however, our initial HSCR-associated “sensitive” haplotype expressed at only 21% of wild-type levels. This significant reduction proved insufficient to rescue the viability when endogenous mouse Ret was deleted. We identified that this deficiency is partially driven by five risk SNPs within established enhancers. Specifically, using CRISPR/Cas9 to restore a conserved So×10 binding site (converting a sensitive SNP to a protective one) increased RET expression by 1.9-fold and restored transcription factor binding. This study provides a robust framework for modeling human-specific regulatory disorders and demonstrates the critical impact of non-coding variation on disease pathogenesis.
Ryan D. Fine, B. Low, Jarod A. Rollins et al.· bioRxiv· 0 citations
Mastitis remains one of the most prevalent and economically significant diseases in dairy cattle worldwide. Although somatic cell count (SCC) is widely used as an indicator for mastitis diagnosis, its physiological variability limits its utility for predicting individual susceptibility. In this study, we aimed to identify genomic markers associated with recurrent clinical mastitis by defining mastitis-susceptible cows as those experiencing three or more episodes within a single lactation. Whole-genome resequencing was conducted on 50 Holstein cows (25 mastitis-susceptible and 25 healthy controls), yielding 536,184 high-quality SNPs after stringent quality control. A genome-wide association analysis identified 86 SNPs surpassing the significance threshold (–log₁₀P > 5.0), and seven candidate SNPs were evaluated in an independent cohort of 100 cows. Notably, the majority of candidate SNPs were localized on the X chromosome, suggesting a potential role for X-linked variation in mastitis immune response and disease resistance. While the candidate marker panel combining SNP1, SNP2, and SNP7 demonstrated moderate sensitivity (47%), its high specificity (98%) highlights its potential utility as a preliminary screening tool for identifying individuals at increased risk of recurrent mastitis. These findings provide a foundation for further functional validation and large-scale replication studies, which are essential for implementing effective genomic selection strategies to mitigate mastitis incidence in dairy herds.
S. Miyata, Lijie Fan, M. Ito et al.· PLoS ONE· 0 citations