GWAS reveals determinants of mobilization rate and dynamics of an active endogenous retrovirus of cattle.

Five to ten percent of mammalian genomes is occupied by multiple clades of endogenous retroviruses (ERVs), that may count thousands of members. New ERV clades arise by retroviral infection of the germline followed by expansion by reinfection and/or retrotransposition. ERV mobilization is a source of deleterious variation, driving the emergence of ERV silencing mechanisms, leaving "DNA fossils". Here we show that the ERVK[2-1-LTR] clade is still active in the bovine and a source of disease-causing alleles. We develop a method to measure the rate of ERVK[2-1-LTR] mobilization, finding an average of 1 per ~150 sperm cells, with >10-fold difference between animals. We perform a genome-wide association study and identify eight loci affecting ERVK[2-1-LTR] mobilization. We provide evidence that polymorphic ERVK[2-1-LTR] elements in four of these loci cause the association. We generate a catalogue of full length ERVK[2-1-LTR] elements, and show that it comprises 15% of C-type autonomous elements, and 85% of D-type non-autonomous elements lacking functional genes. We show that >25% of the variance of mobilization rate is determined by the number of C-type elements, yet that de novo insertions are dominated by D-type elements. We propose that D-type elements act as parasite-of-parasite gene drives that may contribute to the observed demise of ERV elements.

PCIP-seq: simultaneous sequencing of integrated viral genomes and their insertion sites with long reads.

The integration of a viral genome into the host genome has a major impact on the trajectory of the infected cell. Integration location and variation within the associated viral genome can influence both clonal expansion and persistence of infected cells. Methods based on short-read sequencing can identify viral insertion sites, but the sequence of the viral genomes within remains unobserved. We develop PCIP-seq, a method that leverages long reads to identify insertion sites and sequence their associated viral genome. We apply the technique to exogenous retroviruses HTLV-1, BLV, and HIV-1, endogenous retroviruses, and human papillomavirus.

A de novo missense mutation of FGFR2 causes facial dysplasia syndrome in Holstein cattle.

BACKGROUND: Surveillance for bovine genetic diseases in Denmark identified a hitherto unreported congenital syndrome occurring among progeny of a Holstein sire used for artificial breeding. A genetic aetiology due to a dominant inheritance with incomplete penetrance or a mosaic germline mutation was suspected as all recorded cases were progeny of the same sire. Detailed investigations were performed to characterize the syndrome and to reveal its cause. RESULTS: Seven malformed calves were submitted examination. All cases shared a common morphology with the most striking lesions being severe facial dysplasia and complete prolapse of the eyes. Consequently the syndrome was named facial dysplasia syndrome (FDS). Furthermore, extensive brain malformations, including microencephaly, hydrocephalus, lobation of the cerebral hemispheres and compression of the brain were present. Subsequent data analysis of progeny of the sire revealed that around 0.5% of his offspring suffered from FDS. High density single nucleotide polymorphism (SNP) genotyping data of the seven cases and their parents were used to map the defect in the bovine genome. Significant genetic linkage was obtained for three regions, including chromosome 26 where whole genome sequencing of a case-parent trio revealed two de novo variants perfectly associated with the disease: an intronic SNP in the DMBT1 gene and a single non-synonymous variant in the FGFR2 gene. This FGFR2 missense variant (c.927G>T) affects a gene encoding a member of the fibroblast growth factor receptor family, where amino acid sequence is highly conserved between members and across species. It is predicted to change an evolutionary conserved tryptophan into a cysteine residue (p.Trp309Cys). Both variant alleles were proven to result from de novo mutation events in the germline of the sire. CONCLUSIONS: FDS is a novel genetic disorder of Holstein cattle. Mutations in the human FGFR2 gene are associated with various dominant inherited craniofacial dysostosis syndromes. Given the phenotypic similarities in FDS affected calves, the genetic mapping and absence of further high impact variants in the critical genome regions, it is highly likely that the missense mutation in the FGFR2 gene caused the FDS phenotype in a dominant mode of inheritance.

Cis-perturbation of cancer drivers by the HTLV-1/BLV proviruses is an early determinant of leukemogenesis.

Human T-cell leukaemia virus type-1 (HTLV-1) and bovine leukaemia virus (BLV) infect T- and B-lymphocytes, respectively, provoking a polyclonal expansion that will evolve into an aggressive monoclonal leukaemia in ∼5% of individuals following a protracted latency period. It is generally assumed that early oncogenic changes are largely dependent on virus-encoded products, especially TAX and HBZ, while progression to acute leukaemia/lymphoma involves somatic mutations, yet that both are independent of proviral integration site that has been found to be very variable between tumours. Here, we show that HTLV-1/BLV proviruses are integrated near cancer drivers which they affect either by provirus-dependent transcription termination or as a result of viral antisense RNA-dependent cis-perturbation. The same pattern is observed at polyclonal non-malignant stages, indicating that provirus-dependent host gene perturbation contributes to the initial selection of the multiple clones characterizing the asymptomatic stage, requiring additional alterations in the clone that will evolve into full-blown leukaemia/lymphoma.

A missense mutation accelerating the gating of the lysosomal Cl-/H+-exchanger ClC-7/Ostm1 causes osteopetrosis with gingival hamartomas in cattle.

Chloride-proton exchange by the lysosomal anion transporter ClC-7/Ostm1 is of pivotal importance for the physiology of lysosomes and bone resorption. Mice lacking either ClC-7 or Ostm1 develop a lysosomal storage disease and mutations in either protein have been found to underlie osteopetrosis in mice and humans. Some human disease-causing CLCN7 mutations accelerate the usually slow voltage-dependent gating of ClC-7/Ostm1. However, it has remained unclear whether the fastened kinetics is indeed causative for the disease. Here we identified and characterized a new deleterious ClC-7 mutation in Belgian Blue cattle with a severe symptomatology including perinatal lethality and in most cases gingival hamartomas. By autozygosity mapping and genome-wide sequencing we found a handful of candidate variants, including a cluster of three private SNPs causing the substitution of a conserved tyrosine in the CBS2 domain of ClC-7 by glutamine. The case for ClC-7 was strengthened by subsequent examination of affected calves that revealed severe osteopetrosis. The Y750Q mutation largely preserved the lysosomal localization and assembly of ClC-7/Ostm1, but drastically accelerated its activation by membrane depolarization. These data provide first evidence that accelerated ClC-7/Ostm1 gating per se is deleterious, highlighting a physiological importance of the slow voltage-activation of ClC-7/Ostm1 in lysosomal function and bone resorption.

Methylation analysis of the imprinted DLK1-GTL2 domain supports the random parental origin of the IGH-involving del(14q) in B-cell malignancies.

Leukemias/lymphomas with IGH-involving del(14q)(1) commonly lose the DLK1-GTL2 imprinted domain that comprises several paternally and maternally expressed genes, including a cluster of microRNAs. Given that deletion of this region could lead to inactivation of a monoallelically expressed tumor suppressor gene, our study aimed at determination of the parental origin of del(14q/IGH). The designed allele-specific methylation study of the DLK1/GTL2 intergenic differentially methylated region allowed us to determine the parental origin of del(14q/IGH) in 9/20 analyzed cases. In six cases del(14q/IGH) was of the paternal origin and in three cases of the maternal origin. These findings argue against the concept that a TSG/anti-oncomir located in the imprinted region is systematically inactivated by a targeted deletion of its functional allele.

RNAi-mediated allelic trans-interaction at the imprinted Rtl1/Peg11 locus.

The Dlk1-Gtl2 imprinted domain, encompassing the callipyge (CLPG) locus in sheep, has recently been shown to harbor a large number of maternally expressed miRNA genes [1, 2]. Two of these (mir127 and mir136) are processed from a transcript (antiPeg11) that is antisense to Rtl1/Peg11, a paternally expressed intronless gene with homology to the gag and pol polyproteins of Sushi-like retroelements [3]. We herein demonstrate that several additional miRNAs are processed from antiPeg11 and that these regulate Rtl1/Peg11 in trans by guiding RISC-mediated cleavage of its mRNA. This is the first demonstration of miRNA-mediated RNAi involving imprinted genes in mammals.