Human transposon tectonics

Evolution of Virus-like Features and Intrinsically Disordered Regions in Retrotransposon-derived Mammalian Genes

Several mammalian genes have originated from the domestication of retrotransposons, selfish mobile elements related to retroviruses. Some of the proteins encoded by these genes have maintained virus-like features; including self-processing, capsid structure formation, and the generation of different isoforms through -1 programmed ribosomal frameshifting. Using quantitative approaches in molecular evolution and biophysical analyses, we studied 28 retrotransposon-derived genes, with a focus on the evolution of virus-like features. By analyzing the rate of synonymous substitutions, we show that the -1 programmed ribosomal frameshifting mechanism in three of these genes (PEG10, PNMA3, and PNMA5) is conserved across mammals and originates alternative proteins. These genes were targets of positive selection in primates, and one of the positively selected sites affects a B-cell epitope on the spike domain of the PNMA5 capsid, a finding reminiscent of observations in infectious viruses. More generally, we found that retrotransposon-derived proteins vary in their intrinsically disordered region content and this is directly associated with their evolutionary rates. Most positively selected sites in these proteins are located in intrinsically disordered regions and some of them impact protein posttranslational modifications, such as autocleavage and phosphorylation. Detailed analyses of the biophysical properties of intrinsically disordered regions showed that positive selection preferentially targeted regions with lower conformational entropy. Furthermore, positive selection introduces variation in binary sequence patterns across orthologues, as well as in chain compaction. Our results shed light on the evolutionary trajectories of a unique class of mammalian genes and suggest a novel approach to study how intrinsically disordered region biophysical characteristics are affected by evolution.

The Diverse Evolutionary Histories of Domesticated Metaviral Capsid Genes in Mammals

Selfish genetic elements comprise significant fractions of mammalian genomes. In rare instances, host genomes domesticate segments of these elements for function. Using a complete human genome assembly and 25 additional vertebrate genomes, we re-analyzed the evolutionary trajectories and functional potential of capsid (CA) genes domesticated from Metaviridae, a lineage of retrovirus-like retrotransposons. Our study expands on previous analyses to unearth several new insights about the evolutionary histories of these ancient genes. We find that at least five independent domestication events occurred from diverse Metaviridae, giving rise to three universally retained single-copy genes evolving under purifying selection and two gene families unique to placental mammals, with multiple members showing evidence of rapid evolution. In the SIRH/RTL family, we find diverse amino-terminal domains, widespread loss of protein-coding capacity in RTL10 despite its retention in several mammalian lineages, and differential utilization of an ancient programmed ribosomal frameshift in RTL3 between the domesticated CA and protease domains. Our analyses also reveal that most members of the PNMA family in mammalian genomes encode a conserved putative amino-terminal RNA-binding domain (RBD) both adjoining and independent from domesticated CA domains. Our analyses lead to a significant correction of previous annotations of the essential CCDC8 gene. We show that this putative RBD is also present in several extant Metaviridae, revealing a novel protein domain configuration in retrotransposons. Collectively, our study reveals the divergent outcomes of multiple domestication events from diverse Metaviridae in the common ancestor of placental mammals.

LINE-1 hypomethylation, increased retrotransposition and tumor-specific insertion in upper gastrointestinal cancer

The long interspersed nuclear element-1 (LINE-1) retrotransposons are a major family of mobile genetic elements, comprising approximately 17% of the human genome. The methylation state of LINE-1 is often used as an indicator of global DNA methylation levels and it regulates the retrotransposition and somatic insertion of the genetic element. We have previously reported the significant relationship between LINE-1 hypomethylation and poor prognosis in upper gastrointestinal (GI) cancers. However, the causal relationships between LINE-1 hypomethylation, retrotransposition, and tumor-specific insertion in upper GI cancers remain unknown. We used bisulfite-pyrosequencing and quantitative real-time PCR to verify LINE-1 methylation and copy number in tissue samples of 101 patients with esophageal and 103 patients with gastric cancer. Furthermore, we analyzed the LINE-1 retrotransposition profile with an originally developed L1Hs-seq. In tumor samples, LINE-1 methylation levels were significantly lower than non-tumor controls, while LINE-1 copy numbers were markedly increased. As such, there was a significant inverse correlation between the LINE-1 methylation level and copy number in tumor tissues, with lower LINE-1 methylation levels corresponding to higher LINE-1 copy numbers. Of particular importance is that somatic LINE-1 insertions were more numerous in tumor than normal tissues. Furthermore, we observed that LINE-1 was inserted evenly across all chromosomes, and most often within genomic regions associated with tumor-suppressive genes. LINE-1 hypomethylation in upper GI cancers is related to increased LINE-1 retrotransposition and tumor-specific insertion events, which may collectively contribute to the acquisition of aggressive tumor features through the inactivation of tumor-suppressive genes.

Pan cancer characterization of genes whose expression has been associated with LINE-1 antisense promoter activity

BACKGROUND

Long interspersed nuclear element-1 (LINE-1 or L1) comprises 17% of the human genome. As the only autonomous and active retrotransposons, L1 may take part in cancer initiation and progression in some ways. The studies of L1 in cancer mainly focus on the impact of L1 insertion into the new genome locus. The L1 5´ untranslated region (UTR) also contains antisense promoter (ASP) activity, generating L1-gene chimeric transcripts to a neighbor exon. Some of these ASP-associated genes have been reported to be overexpressed in cancer and promote cancer cell growth. However, little is known about overall expression patterns and the roles of L1 ASP-associated genes in human cancers.

RESULTS

L1 ASP-associated genes were frequently dysregulated in cancer and associated with the cell cycle, the PI3K/AKT pathway, and the GTPase signaling pathway. The expression of L1 ASP-associated genes was correlated with tumor patient prognosis. Hub L1 ASP-associated genes CENPU and MCM2 showed a correlation with immune infiltration, clinical T stage, and cancer stemness in pan-cancer. Knockdown of L1 ASP-associated gene LINC00491 resulted in a significant decrease in tumor growth and migration ability.

CONCLUSIONS

The expression of L1 ASP-associated genes is significantly dysregulated at the pan-cancer level, which is closely related to the tumor microenvironment, progression, and patient prognosis. Hub genes CENPU and MCM2 are expected to be new tumor diagnostic markers and therapeutic targets.

The diverse evolutionary histories of domesticated metaviral capsid genes in mammals

Selfish genetic elements and their remnants comprise at least half of the human genome. Active transposons duplicate by inserting copies at new sites in a host genome. Following insertion, transposons can acquire mutations that render them inactive; the accrual of additional mutations can render them unrecognizable over time. However, in rare instances, segments of transposons become useful for the host, in a process called gene domestication. Using the first complete human genome assembly and 25 additional vertebrate genomes, we analyzed the evolutionary trajectories and functional potential of genes domesticated from the capsid genes of Metaviridae, a retroviral-like retrotransposon family. Our analysis reveals four families of domesticated capsid genes in placental mammals with varied evolutionary outcomes, ranging from universal retention to lineage-specific duplications or losses and from purifying selection to lineage-specific rapid evolution. The four families of domesticated capsid genes have divergent amino-terminal domains, inherited from four distinct ancestral metaviruses. Structural predictions reveal that many domesticated genes encode a previously unrecognized RNA-binding domain retained in multiple paralogs in mammalian genomes both adjacent to and independent from the capsid domain. Collectively, our study reveals diverse outcomes of domestication of diverse metaviruses, which led to structurally and evolutionarily diverse genes that encode important, but still largely-unknown functions in placental mammals. (207).

Themes and variations on piRNA-guided transposon control

PIWI-interacting RNAs (piRNAs) are responsible for preventing the movement of transposable elements in germ cells and protect the integrity of germline genomes. In this review, we examine the common elements of piRNA-guided silencing as well as the differences observed between species. We have categorized the mechanisms of piRNA biogenesis and function into modules. Individual PIWI proteins combine these modules in various ways to produce unique PIWI-piRNA pathways, which nevertheless possess the ability to perform conserved functions. This modular model incorporates conserved core mechanisms and accommodates variable co-factors. Adaptability is a hallmark of this RNA-based immune system. We believe that considering the differences in germ cell biology and resident transposons in different organisms is essential for placing the variations observed in piRNA biology into context, while still highlighting the conserved themes that underpin this process.

Abundance and Diversification of Repetitive Elements in Decapoda Genomes

Repetitive elements are a major component of DNA sequences due to their ability to propagate through the genome. Characterization of Metazoan repetitive profiles is improving; however, current pipelines fail to identify a significant proportion of divergent repeats in non-model organisms. The Decapoda order, for which repeat content analyses are largely lacking, is characterized by extremely variable genome sizes that suggest an important presence of repetitive elements. Here, we developed a new standardized pipeline to annotate repetitive elements in non-model organisms, which we applied to twenty Decapoda and six other Crustacea genomes. Using this new tool, we identified 10% more repetitive elements than standard pipelines. Repetitive elements were more abundant in Decapoda species than in other Crustacea, with a very large number of highly repeated satellite DNA families. Moreover, we demonstrated a high correlation between assembly size and transposable elements and different repeat dynamics between Dendrobranchiata and Reptantia. The patterns of repetitive elements largely reflect the phylogenetic relationships of Decapoda and the distinct evolutionary trajectories within Crustacea. In summary, our results highlight the impact of repetitive elements on genome evolution in Decapoda and the value of our novel annotation pipeline, which will provide a baseline for future comparative analyses.

A pilot study of LINE-1 copy number and telomere length with aging in human sperm

PURPOSE

Unlike other cells in the body, in sperm, telomere length (TL) increases with age. TL can regulate nearby genes, and the subtelomeric region is rich in retrotransposons. We hypothesized that age-related telomere lengthening in sperm might suppress Long Interspersed Element 1 (LINE-1/L1), the only competent retrotransposon in humans.

METHODS

We measured L1 copy number (L1-CN) and sperm telomere length (STL) from young and older men to evaluate the relationship between age, TL and L1-CN. We also evaluated L1-CN and TL in individual sperm to determine whether these variables influence sperm morphology. STL was assayed by Multiplex quantitative polymerase chain reaction method (mmqPCR) and L1-CN by Quantitative polymerase chain reaction (qPCR).

RESULTS

We found that STL increased, and L1-CN decreased significantly with paternal age. STL in normal single sperm was significantly higher than in abnormal sperm. L1-CN did not differ between normal and abnormal sperm. Furthermore, morphologically normal sperm have longer telomeres than abnormal sperm.

CONCLUSIONS

Elongation of telomeres in the male germline could repress retrotransposition, which tends to increase with cellular aging. More studies in larger cohorts across a wide age span are needed to confirm our conclusions and explore their biological and clinical significance.

Transposon insertion profiling by sequencing (TIPseq) identifies novel LINE-1 insertions in human sperm

PURPOSE

Long interspersed nuclear element-1 (LINE-1 or L1) comprises 17% of the human genome. Retrotransposons may perturb gene integrity or alter gene expression by altering regulatory regions in the genome. The germline employs a number of mechanisms, including cytosine methylation, to repress retrotransposon transcription throughout most of life. Demethylation during germ cell and early embryo development de-represses retrotransposons. Intriguingly, de novo genetic variation appearing in sperm has been implicated in a number of disorders in offspring, including autism spectrum disorder, schizophrenia, and bipolar disorder. We hypothesize that human sperm exhibit de novo retrotransposition and employ a new sequencing method, single cell transposon insertion profiling by sequencing (scTIPseq) to map them in small amounts of human sperm.

METHODS

Cross-sectional case-control study of sperm samples (n=10 men; ages 32-55 years old) from consenting men undergoing IVF at NYU Langone Fertility Center. scTIPseq identified novel LINE-1 insertions in individual sperm and TIPseqHunter, a custom bioinformatics pipeline, compared the architecture of sperm LINE-1 to known LINE-1 insertions from the European database of Human specific LINE-1 (L1Hs) retrotransposon insertions (euL1db).

RESULTS

scTIPseq identified 17 novel insertions in sperm. New insertions were mainly intergenic or intronic. Only one sample did not exhibit new insertions. The location or number of novel insertions did not differ by paternal age.

CONCLUSION

This study for the first time reports novel LINE-1 insertions in human sperm, demonstrating the feasibility of scTIPseq, and identifies new contributors to genetic diversity in the human germ line.

Loss of epigenetic suppression of retrotransposons with oncogenic potential in aging mammary luminal epithelial cells

A primary function of DNA methylation in mammalian genomes is to repress transposable elements (TEs). The widespread methylation loss that is commonly observed in cancer cells results in the loss of epigenetic repression of TEs. The aging process is similarly characterized by changes to the methylome. However, the impact of these epigenomic alterations on TE silencing and the functional consequences of this have remained unclear. To assess the epigenetic regulation of TEs in aging, we profiled DNA methylation in human mammary luminal epithelial cells (LEps)-a key cell lineage implicated in age-related breast cancers-from younger and older women. We report here that several TE subfamilies function as regulatory elements in normal LEps, and a subset of these display consistent methylation changes with age. Methylation changes at these TEs occurred at lineage-specific transcription factor binding sites, consistent with loss of lineage specificity. Whereas TEs mainly showed methylation loss, CpG islands (CGIs) that are targets of the Polycomb repressive complex 2 (PRC2) show a gain of methylation in aging cells. Many TEs with methylation loss in aging LEps have evidence of regulatory activity in breast cancer samples. We furthermore show that methylation changes at TEs impact the regulation of genes associated with luminal breast cancers. These results indicate that aging leads to DNA methylation changes at TEs that undermine the maintenance of lineage specificity, potentially increasing susceptibility to breast cancer.

Insertion of 643bp Retrotransposon Upstream of PPARγ CDS Is Associated with Backfat of Large White Pigs

PPARs are essential regulators of mammalian fatty acid and lipid metabolism. Although the effects of genetic variations, including single nucleotide polymorphisms (SNPs) in PPARs genes on the phenotype of domestic animals have been investigated, there is limited information on the impact of retrotransposon insertion polymorphisms (RIPs). In this study, a combined comparative genome and polymerase chain reaction (PCR) was used to excavate the RIPs in porcine PPARs. We also investigated the potential effects of retrotransposon insertion on phenotype and expression patterns. This study identified the two RIPs in PPARs genes, namely an ERV in intron 1 of PPARα and a combined retrotransposon in intron 2 of PPARγ, designated as PPARα-ERV-RIP and PPARγ-COM-RIP, respectively. These RIPs exhibited different distribution patterns among Chinese indigenous breeds and Western commercial breeds. Individuals with the PPARα-ERV-RIP^+/+ genotype (+/+ indicated homozygous with insertion) among Large White pigs had significantly higher (p < 0.05) corrected backfat thickness compared to those with the other two genotypes. Similarly, those with the PPARγ-COM-RIP^-/- genotype had significantly higher (p < 0.05) corrected backfat thickness than those with the other two genotypes in Large White pigs. Moreover, in 30-day-old Sujiang piglets, the PPARγ gene expression in the backfat of those with the PPARγ-COM-RIP^-/- genotype (-/- indicated homozygous without insertion) was significantly greater (p < 0.01) than those with other genotypes. The dual luciferase reporter gene assay demonstrated that the combined retrotransposon insertion significantly reduced the activity of the MYC promoter in both C2C12 and 3T3-L1 cells (p < 0.01). Therefore, the combined retrotransposon insertion could function as a repressor to decrease the expression of PPARγ, making PPARγ-COM-RIP a valuable molecular marker for assisted selection of backfat thickness in pig breeding.

Human endogenous retrovirus onco-exaptation counters cancer cell senescence through calbindin

Increased levels and diversity of human endogenous retrovirus (HERV) transcription characterize most cancer types and are linked with disease outcomes. However, the underlying processes are incompletely understood. Here, we show that elevated transcription of HERVH proviruses predicted survival of lung squamous cell carcinoma (LUSC) and identified an isoform of CALB1, encoding calbindin, ectopically driven by an upstream HERVH provirus under the control of KLF5, as the mediator of this effect. HERVH-CALB1 expression was initiated in preinvasive lesions and associated with their progression. Calbindin loss in LUSC cell lines impaired in vitro and in vivo growth and triggered senescence, consistent with a protumor effect. However, calbindin also directly controlled the senescence-associated secretory phenotype (SASP), marked by secretion of CXCL8 and other neutrophil chemoattractants. In established carcinomas, CALB1-negative cancer cells became the dominant source of CXCL8, correlating with neutrophil infiltration and worse prognosis. Thus, HERVH-CALB1 expression in LUSC may display antagonistic pleiotropy, whereby the benefits of escaping senescence early during cancer initiation and clonal competition were offset by the prevention of SASP and protumor inflammation at later stages.

CancerHERVdb: Human Endogenous Retrovirus (HERV) Expression Database for Human Cancer Accelerates Studies of the Retrovirome and Predictions for HERV-Based Therapies

In this study, we sought to create a database summarizing the expression of human endogenous retroviruses (HERVs) in various human cancers. HERVs are suitable therapeutic targets due to their abundance in the human genome, overexpression in various malignancies, and involvement in various cancer pathways. We identified articles on HERVs from PubMed and then prescreened and automatically categorized them using the portable document format (PDF) data extractor (PDE) R package. We discovered 196 primary research articles with HERV expression data from cancer tissues or cancer cell lines. HERV RNA and protein expression was reported in brain, breast, cervical, colorectal, endocrine, gastrointestinal, kidney/renal/pelvis, liver, lung, genital, oral cavity, pharynx, ovary, pancreas, prostate, skin, testicular, urinary/bladder, and uterus cancers, leukemias, lymphomas, and myelomas. Additionally, we discovered reports of HERV RNA-only overexpression in soft tissue cancers including heart, thyroid, bone, and joint cancers. The CancerHERVdb database is hosted in the form of interactive visualizations of the expression data and a summary data table at https://erikstricker.shinyapps.io/cancerHERVdb/. The user can filter the findings according to cancer type, HERV family, HERV gene, or a combination thereof and easily export the results with the corresponding reference list. In our report, we provide examples of potential uses of the CancerHERVdb, such as identification of cancers suitable for off-target treatment with the multiple sclerosis-associated retrovirus (MSRV)-Env-targeting antibody GNbAC1 (now named temelimab) currently in phase 2b clinical trials for multiple sclerosis or the discovery of cancers overexpressing HERV-H long terminal repeat-associating 2 (HHLA2), a newly emerging immune checkpoint. In summary, the CancerHERVdb allows cross-study comparisons, encourages data exploration, and informs about potential off-target effects of HERV-targeting treatments. IMPORTANCE Human endogenous retroviruses (HERVs), which in the past have inserted themselves in various regions of the human genome, are to various degrees activated in virtually every cancer type. While a centralized naming system and resources summarizing HERV levels in cancers are lacking, the CancerHERVdb database provides a consolidated resource for cross-study comparisons, data exploration, and targeted searches of HERV activation. The user can access data extracted from hundreds of articles spanning 25 human cancer categories. Therefore, the CancerHERVdb database can aid in the identification of prognostic and risk markers, drivers of cancer, tumor-specific targets, multicancer spanning signals, and targets for immune therapies. Consequently, the CancerHERVdb database is of direct relevance for clinical as well as basic research.

LINE-1 retrotransposon expression in cancerous, epithelial and neuronal cells revealed by 5' single-cell RNA-Seq

LINE-1 retrotransposons are sequences capable of copying themselves to new genomic loci via an RNA intermediate. New studies implicate LINE-1 in a range of diseases, especially in the context of aging, but without an accurate understanding of where and when LINE-1 is expressed, a full accounting of its role in health and disease is not possible. We therefore developed a method-5' scL1seq-that makes use of a widely available library preparation method (10x Genomics 5' single cell RNA-seq) to measure LINE-1 expression in tens of thousands of single cells. We recapitulated the known pattern of LINE-1 expression in tumors-present in cancer cells, absent from immune cells-and identified hitherto undescribed LINE-1 expression in human epithelial cells and mouse hippocampal neurons. In both cases, we saw a modest increase with age, supporting recent research connecting LINE-1 to age related diseases.

HERVs and Cancer-A Comprehensive Review of the Relationship of Human Endogenous Retroviruses and Human Cancers

Genomic instability and genetic mutations can lead to exhibition of several cancer hallmarks in affected cells such as sustained proliferative signaling, evasion of growth suppression, activated invasion, deregulation of cellular energetics, and avoidance of immune destruction. Similar biological changes have been observed to be a result of pathogenic viruses and, in some cases, have been linked to virus-induced cancers. Human endogenous retroviruses (HERVs), once external pathogens, now occupy more than 8% of the human genome, representing the merge of genomic and external factors. In this review, we outline all reported effects of HERVs on cancer development and discuss the HERV targets most suitable for cancer treatments as well as ongoing clinical trials for HERV-targeting drugs. We reviewed all currently available reports of the effects of HERVs on human cancers including solid tumors, lymphomas, and leukemias. Our review highlights the central roles of HERV genes, such as gag, env, pol, np9, and rec in immune regulation, checkpoint blockade, cell differentiation, cell fusion, proliferation, metastasis, and cell transformation. In addition, we summarize the involvement of HERV long terminal repeat (LTR) regions in transcriptional regulation, creation of fusion proteins, expression of long non-coding RNAs (lncRNAs), and promotion of genome instability through recombination.

Dynamics of histone acetylation during human early embryogenesis

It remains poorly understood about the regulation of gene and transposon transcription during human early embryogenesis. Here, we report that broad H3K27ac domains are genome-widely distributed in human 2-cell and 4-cell embryos and transit into typical peaks in the 8-cell embryos. The broad H3K27ac domains in early embryos before zygotic genome activation (ZGA) are also observed in mouse. It suggests that broad H3K27ac domains play conserved functions before ZGA in mammals. Intriguingly, a large portion of broad H3K27ac domains overlap with broad H3K4me3 domains. Further investigation reveals that histone deacetylases are required for the removal or transition of broad H3K27ac domains and ZGA. After ZGA, the number of typical H3K27ac peaks is dynamic, which is associated with the stage-specific gene expression. Furthermore, P300 is important for the establishment of H3K27ac peaks and the expression of associated genes in early embryos after ZGA. Our data also indicate that H3K27ac marks active transposons in early embryos. Interestingly, H3K27ac and H3K18ac signals rather than H3K9ac signals are enriched at ERVK elements in mouse embryos after ZGA. It suggests that different types of histone acetylations exert distinct roles in the activation of transposons. In summary, H3K27ac modification undergoes extensive reprogramming during early embryo development in mammals, which is associated with the expression of genes and transposons.

Quantification of LINE-1 RNA Expression from Bulk RNA-seq Using L1EM

LINE-1 retrotransposons have the potential to cause DNA damage, contribute to genome instability, and induce an interferon response. Thus, accurate measurements of their expression, especially in disease contexts where genome instability and the interferon response are relevant, are of particular importance. Illumina-based bulk RNA sequencing remains the most abundant datatype for measuring gene expression. However, "active" expression from its own internal promoter is only one source of LINE-1 aligning reads in an RNA-seq experiment. With about half a million LINE-1 sequences scattered throughout the genome, many are incorporated into other transcripts that have nothing to do with LINE-1 activity. We call this "passive" co-transcription. Here we will describe how to use L1EM, a computational method that separates active from passive LINE-1 expression at the locus-specific level.

Affinity-Based Interactome Analysis of Endogenous LINE-1 Macromolecules

During their proliferation and the host's concomitant attempts to suppress it, LINE-1 (L1) retrotransposons give rise to a collection of heterogeneous ribonucleoproteins (RNPs); their protein and RNA compositions remain poorly defined. The constituents of L1-associated macromolecules can differ depending on numerous factors, including, for example, position within the L1 life cycle, whether the macromolecule is productive or under suppression, and the cell type within which the proliferation is occurring. This chapter describes techniques that aid the capture and characterization of protein and RNA components of L1 macromolecules from tissues that natively express them. The protocols described have been applied to embryonal carcinoma cell lines that are popular model systems for L1 molecular biology (e.g., N2102Ep, NTERA-2, and PA-1 cells), as well as colorectal cancer tissues. N2102Ep cells are given as the use case for this chapter; the protocols should be applicable to essentially any tissue exhibiting endogenous L1 expression with minor modifications.

As3MT-mediated SAM consumption, which inhibits the methylation of histones and LINE1, is involved in arsenic-induced male reproductive damage

Studies have demonstrated that arsenic (As) induces male reproductive injury, however, the mechanism remains unknown. The high levels of arsenic (3) methyltransferase (As3MT) promote As-induced male reproductive toxicity. For As-exposed mice, the germ cells in seminiferous tubules and sperm quality were reduced. Exposure to As caused lower S-adenosylmethionine (SAM) and 5-methylcytosine (5 mC) levels, histone and DNA hypomethylation, upregulation of long interspersed element class 1 (LINE1, or L1), defective repair of double-strand breaks (DSBs), and the arrest of meiosis, resulting in apoptosis of germ cells and lower litter size. For GC-2spd (GC-2) cells, As induced apoptosis, which was prevented by adding SAM or by reducing the expression of As3MT. The levels of LINE1, affected by SAM content, were involved in As-induced apoptosis. Furthermore, folic acid (FA) and vitamin B12 (VB₁₂) supplements restored SAM, 5 mC, and LINE1 levels and blocked impairment of spermatogenesis and testes and lower litter size. Exposed to As, mice with As3MT knockdown showed less impairment of spermatogenesis and testes and greater litter size compared to As-exposed wild-type (WT) mice. Thus, the high As3MT levels induced by As consume SAM and block histone and LINE1 DNA methylation, elevating LINE1 expression and evoking impairment of spermatogenesis, which causes male reproductive damage. Overall, we have found a mechanism for As-induced male reproductive damage, which provides biological insights into the alleviation of reproductive injury induced by environmental factors.

The balance between NANOG and SOX17 mediated by TET proteins regulates specification of human primordial germ cell fate

Background

Human primordial germ cells (hPGCs) initiate from the early post-implantation embryo at week 2–3 and undergo epigenetic reprogramming during development. However, the regulatory mechanism of DNA methylation during hPGC specification is still largely unknown due to the difficulties in analyzing early human embryos. Using an in vitro model of hPGC induction, we found a novel function of TET proteins and NANOG in the hPGC specification which was different from that discovered in mice.

Methods

Using the CRISPR–Cas9 system, we generated a set of TET1, TET2 and TET3 knockout H1 human embryonic stem cell (hESC) lines bearing a BLIMP1-2A-mKate2 reporter. We determined the global mRNA transcription and DNA methylation profiles of pluripotent cells and induced hPGC-like cells (hPGCLCs) by RNA-seq and whole-genome bisulfite sequencing (WGBS) to reveal the involved signaling pathways after TET proteins knockout. ChIP-qPCR was performed to verify the binding of TET and NANOG proteins in the SOX17 promoter. Real-time quantitative PCR, western blot and immunofluorescence were performed to measure gene expression at mRNA and protein levels. The efficiency of hPGC induction was evaluated by FACS.

Results

In humans, TET1, TET2 and TET3 triple-knockout (TKO) human embryonic stem cells (hESCs) impaired the NODAL signaling pathway and impeded hPGC specification in vitro, while the hyperactivated NODAL signaling pathway led to gastrulation failure when Tet proteins were inactivated in mouse. Specifically, TET proteins stimulated SOX17 through the NODAL signaling pathway and directly regulates NANOG expression at the onset of hPGCLCs induction. Notably, NANOG could bind to SOX17 promoter to regulate its expression in hPGCLCs specification. Furthermore, in TKO hESCs, DNMT3B-mediated hypermethylation of the NODAL signaling-related genes and NANOG/SOX17 promoters repressed their activation and inhibited hPGCLC induction. Knockout of DNMT3B in TKO hESCs partially restored NODAL signaling and NANOG/SOX17 expression, and rescued hPGCLC induction.

Conclusion

Our results show that TETs-mediated oxidation of 5-methylcytosine modulates the NODAL signaling pathway and its downstream genes, NANOG and SOX17, by promoting demethylation in opposition to DNMT3B-mediated methylation, suggesting that the epigenetic balance of DNA methylation and demethylation in key genes plays a fundamental role in early hPGC specification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00917-0.

New Perspectives on the Importance of Cell-Free DNA Biology

Body fluids are constantly replenished with a population of genetically diverse cell-free DNA (cfDNA) fragments, representing a vast reservoir of information reflecting real-time changes in the host and metagenome. As many body fluids can be collected non-invasively in a one-off and serial fashion, this reservoir can be tapped to develop assays for the diagnosis, prognosis, and monitoring of wide-ranging pathologies, such as solid tumors, fetal genetic abnormalities, rejected organ transplants, infections, and potentially many others. The translation of cfDNA research into useful clinical tests is gaining momentum, with recent progress being driven by rapidly evolving preanalytical and analytical procedures, integrated bioinformatics, and machine learning algorithms. Yet, despite these spectacular advances, cfDNA remains a very challenging analyte due to its immense heterogeneity and fluctuation in vivo. It is increasingly recognized that high-fidelity reconstruction of the information stored in cfDNA, and in turn the development of tests that are fit for clinical roll-out, requires a much deeper understanding of both the physico-chemical features of cfDNA and the biological, physiological, lifestyle, and environmental factors that modulate it. This is a daunting task, but with significant upsides. In this review we showed how expanded knowledge on cfDNA biology and faithful reverse-engineering of cfDNA samples promises to (i) augment the sensitivity and specificity of existing cfDNA assays; (ii) expand the repertoire of disease-specific cfDNA markers, thereby leading to the development of increasingly powerful assays; (iii) reshape personal molecular medicine; and (iv) have an unprecedented impact on genetics research.

Isoform-specific and ubiquitination dependent recruitment of Tet1 to replicating heterochromatin modulates methylcytosine oxidation

Oxidation of the epigenetic DNA mark 5-methylcytosine by Tet dioxygenases is an established route to diversify the epigenetic information, modulate gene expression and overall cellular (patho-)physiology. Here, we demonstrate that Tet1 and its short isoform Tet1s exhibit distinct nuclear localization during DNA replication resulting in aberrant cytosine modification levels in human and mouse cells. We show that Tet1 is tethered away from heterochromatin via its zinc finger domain, which is missing in Tet1s allowing its targeting to these regions. We find that Tet1s interacts with and is ubiquitinated by CRL4(VprBP). The ubiquitinated Tet1s is then recognized by Uhrf1 and recruited to late replicating heterochromatin. This leads to spreading of 5-methylcytosine oxidation to heterochromatin regions, LINE 1 activation and chromatin decondensation. In summary, we elucidate a dual regulation mechanism of Tet1, contributing to the understanding of how epigenetic information can be diversified by spatio-temporal directed Tet1 catalytic activity.

Chromosomal-level assembly of Bactericera cockerelli reveals rampant gene family expansions impacting genome structure, function and insect-microbe-plant-interactions

Lineage specific expansions and gene duplications are some of the most important sources of evolutionary novelty in eukaryotes. Although not as prevalent in eukaryotes compared to bacteria, horizontal gene transfer events can also result in key adaptations for insects, especially for those involved in insect-microbe interactions. In this study we assemble the first chromosomal assembly of the psyllid Bactericera cockerelli and reveal that the B. cockerelli genome has experienced significantly more gene expansion events compared to other Hemipteran representatives with fully sequenced genomes. We also reveal that B. cockerelli's genome is the largest psyllid genome (567 Mb) sequenced to date and is ~15% larger than the other two psyllid species genomes sequenced (Pachypsylla venusta and Diaphorina citri). Structurally, B. cockerelli appears to have an additional chromosome compared to the distantly related psyllid species P. venusta due to a previous chromosomal fission or fusion event. The increase in genome size and dynamic nature of the B. cockerelli genome may largely be contributed to the widespread expansion of type I and type II repeat elements that are rampant across all of B. cockerelli's. chromosomes. These repeat elements are distributed near equally in both euchromatic and heterochromatic regions. Furthermore, significant gene family expansions and gene duplications were uncovered for genes that are expected to be important in its adaptation to insect-plant and microbe interactions, which include transcription factors, proteases, odorant receptors, and horizontally transferred genes that are involved in the nutritional symbioses with their long-term nutritional endosymbiont Carsonella.

3D genome, on repeat: Higher-order folding principles of the heterochromatinized repetitive genome

Nearly half of the human genome is comprised of diverse repetitive sequences ranging from satellite repeats to retrotransposable elements. Such sequences are susceptible to stepwise expansions, duplications, inversions, and recombination events which can compromise genome function. In this review, we discuss the higher-order folding mechanisms of compartmentalization and loop extrusion and how they shape, and are shaped by, heterochromatin. Using primarily mammalian model systems, we contrast mechanisms governing H3K9me3-mediated heterochromatinization of the repetitive genome and highlight emerging links between repetitive elements and chromatin folding.

The Happy Hopping of Transposons: The Origins of V(D)J Recombination in Adaptive Immunity

Nearly 50% of the human genome is derived from transposable elements (TEs). Though dysregulated transposons are deleterious to humans and can lead to diseases, co-opted transposons play an important role in generating alternative or new DNA sequence combinations to perform novel cellular functions. The appearance of an adaptive immune system in jawed vertebrates, wherein the somatic rearrangement of T and B cells generates a repertoire of antibodies and receptors, is underpinned by Class II TEs. This review follows the evolution of recombination activation genes (RAGs), components of adaptive immunity, from TEs, focusing on the structural and mechanistic similarities between RAG recombinases and DNA transposases. As evolution occurred from a transposon precursor, DNA transposases developed a more targeted and constrained mechanism of mobilization. As DNA repair is integral to transposition and recombination, we note key similarities and differences in the choice of DNA repair pathways following these processes. Understanding the regulation of V(D)J recombination from its evolutionary origins may help future research to specifically target RAG proteins to rectify diseases associated with immune dysregulation.

A retrotransposon storm marks clinical phenoconversion to late-onset Alzheimer's disease

Recent reports have suggested that the reactivation of otherwise transcriptionally silent transposable elements (TEs) might induce brain degeneration, either by dysregulating the expression of genes and pathways implicated in cognitive decline and dementia or through the induction of immune-mediated neuroinflammation resulting in the elimination of neural and glial cells. In the work we present here, we test the hypothesis that differentially expressed TEs in blood could be used as biomarkers of cognitive decline and development of AD. To this aim, we used a sample of aging subjects (age > 70) that developed late-onset Alzheimer’s disease (LOAD) over a relatively short period of time (12–48 months), for which blood was available before and after their phenoconversion, and a group of cognitive stable subjects as controls. We applied our developed and validated customized pipeline that allows the identification, characterization, and quantification of the differentially expressed (DE) TEs before and after the onset of manifest LOAD, through analyses of RNA-Seq data. We compared the level of DE TEs within more than 600,000 TE-mapping RNA transcripts from 25 individuals, whose specimens we obtained before and after their phenotypic conversion (phenoconversion) to LOAD, and discovered that 1790 TE transcripts showed significant expression differences between these two timepoints (logFC ± 1.5, logCMP > 5.3, nominal p value < 0.01). These DE transcripts mapped both over- and under-expressed TE elements. Occurring before the clinical phenoconversion, this TE storm features significant increases in DE transcripts of LINEs, LTRs, and SVAs, while those for SINEs are significantly depleted. These dysregulations end with signs of manifest LOAD. This set of highly DE transcripts generates a TE transcriptional profile that accurately discriminates the before and after phenoconversion states of these subjects. Our findings suggest that a storm of DE TEs occurs before phenoconversion from normal cognition to manifest LOAD in risk individuals compared to controls, and may provide useful blood-based biomarkers for heralding such a clinical transition, also suggesting that TEs can indeed participate in the complex process of neurodegeneration.

RNA sensor MDA5 suppresses LINE-1 retrotransposition by regulating the promoter activity of LINE-1 5'-UTR

Background

Type 1 long interspersed elements, or LINE-1, are the only retroelements that replicate autonomously in human cells. The retrotransposition process of LINE-1 can trigger the activation of the innate immune system and has been proposed to play a role in the development of several autoimmune diseases, including Aicardi-Goutières syndrome (AGS). In contrast, all known AGS-associated proteins, except MDA5, have been reported to affect LINE-1 activity. Thus, MDA5 is likely to also function as a LINE-1 suppressor.

Results

MDA5 was found to potently suppress LINE-1 activity in a reporter-based LINE-1 retrotransposition assay. Although MDA5 is an endogenous RNA sensor able to activate the innate immune system, increased interferon (IFN) expression only contributed in part to MDA5-mediated LINE-1 suppression. Instead, MDA5 potently regulated the promoter activity of LINE-1 5′-UTR, as confirmed by transiently expressed myc-tagged MDA5 or knockdown of endogenous MDA5 expression. Consequently, MDA5 effectively reduced the generation of LINE-1 RNA and the subsequent expression of LINE-1 ORF1p and ORF2p. Interestingly, despite MDA5 being a multi-domain protein, the N-terminal 2CARD domain alone is sufficient to interact with LINE-1 5′-UTR and inhibit LINE-1 promoter activity.

Conclusion

Our data reveal that MDA5 functions as a promoter regulator; it directly binds to the LINE-1 5′-UTR and suppresses its promoter activity. Consequently, MDA5 reduces LINE-1 RNA and protein levels, and ultimately inhibits LINE-1 retrotransposition. In contrast, MDA5-induced IFN expression only plays a mild role in MDA5-mediated LINE-1 suppression. In addition, the N-terminal 2CARD domain was found to be a functional region for MDA5 upon inhibition of LINE-1 replication. Thus, our data suggest that besides being an initiator of the innate immune system, MDA5 is also an effector against LINE-1 activity, potentially forming a feedback loop by suppressing LINE-1-induced innate immune activation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00268-0.

Identification and Characterization of Extrachromosomal Circular DNA in Human Placentas With Fetal Growth Restriction

Many studies have confirmed that extrachromosomal circular DNAs (eccDNAs/ecDNAs) exist in tumor and normal cells independently of the chromosome and are essential for oncogene plasticity and drug resistance. Studies have confirmed that there are many eccDNAs/ecDNAs in maternal plasma derived from the fetus. Fetal growth restriction (FGR) is a pregnancy-related disease associated with high newborn morbidity and mortality. However, the characteristics and nature of eccDNAs/ecDNAs in FGR are poorly understood. This study aims to deconstruct the properties and potential functions of eccDNAs/ecDNAs in FGR. We performed circle-seq to identify the expression profile of eccDNAs/ecDNAs, analyzed by bioinformatics, and verified by real-time Polymerase Chain Reaction (PCR) combined with southern blot in FGR compared with the normal groups. A total of 45,131 eccDNAs/ecDNAs (including 2,118 unique ones) were identified, which had significantly higher abundance in FRG group than in normal group, and was bimodal in length, peaking at ~146bp and ~340bp, respectively. Gestational age may be one independent factor affecting the production of eccDNAs/ecDNAs, most of which come from genomic regions with high gene density, with a 4~12bp repeat around the junction, and their origin had a certain genetic preference. In addition, some of the host-genes overlapped with non-coding RNAs (ncRNAs) partially or even completely. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed that host-genes on the differentially expressed eccDNAs/ecDNAs (DEEECs/DEECs) were mainly enriched in immune-related functions and pathways. The presence of some ecDNAs were verified, and whose variability were consistent with the circle-seq results. We identified and characterized eccDNAs/ecDNAs in placentas with FGR, and elucidated the formation mechanisms and the networks with ncRNAs, which provide a new vision for the screening of new biomarkers and therapeutic targets for FGR.

Gene drives gaining speed

Gene drives are selfish genetic elements that are transmitted to progeny at super-Mendelian (>50%) frequencies. Recently developed CRISPR-Cas9-based gene-drive systems are highly efficient in laboratory settings, offering the potential to reduce the prevalence of vector-borne diseases, crop pests and non-native invasive species. However, concerns have been raised regarding the potential unintended impacts of gene-drive systems. This Review summarizes the phenomenal progress in this field, focusing on optimal design features for full-drive elements (drives with linked Cas9 and guide RNA components) that either suppress target mosquito populations or modify them to prevent pathogen transmission, allelic drives for updating genetic elements, mitigating strategies including trans-complementing split-drives and genetic neutralizing elements, and the adaptation of drive technology to other organisms. These scientific advances, combined with ethical and social considerations, will facilitate the transparent and responsible advancement of these technologies towards field implementation.

Alu insertion variants alter gene transcript levels

Alu are high copy number interspersed repeats that have accumulated near genes during primate and human evolution. They are a pervasive source of structural variation in modern humans. Impacts that Alu insertions may have on gene expression are not well understood, although some have been associated with expression quantitative trait loci (eQTLs). Here, we directly test regulatory effects of polymorphic Alu insertions in isolation of other variants on the same haplotype. To screen insertion variants for those with such effects, we used ectopic luciferase reporter assays and evaluated 110 Alu insertion variants, including more than 40 with a potential role in disease risk. We observed a continuum of effects with significant outliers that up- or down-regulate luciferase activity. Using a series of reporter constructs, which included genomic context surrounding the Alu, we can distinguish between instances in which the Alu disrupts another regulator and those in which the Alu introduces new regulatory sequence. We next focused on three polymorphic Alu loci associated with breast cancer that display significant effects in the reporter assay. We used CRISPR to modify the endogenous sequences, establishing cell lines varying in the Alu genotype. Our findings indicate that Alu genotype can alter expression of genes implicated in cancer risk, including PTHLH, RANBP9, and MYC These data show that commonly occurring polymorphic Alu elements can alter transcript levels and potentially contribute to disease risk.

Long Terminal Repeat Retrotransposon Afut4 Promotes Azole Resistance of Aspergillus fumigatus by Enhancing the Expression of sac1 Gene

Aspergillus fumigatus causes a series of invasive diseases, including the high-mortality invasive aspergillosis, and has been a serious global health threat because of its increased resistance to the first-line clinical triazoles. We analyzed the whole-genome sequence of 15 A. fumigatus strains from China and found that long terminal repeat retrotransposons (LTR-RTs), including Afut1, Afut2, Afut3, and Afut4, are most common and have the largest total nucleotide length among all transposable elements in A. fumigatus. Deleting one of the most enriched Afut4_977-sac1 in azole-resistant strains decreased azole resistance and downregulated its nearby gene, sac1, but it did not significantly affect the expression of genes of the ergosterol synthesis pathway. We then discovered that 5'LTR of Afut4_977-sac1 had promoter activity and enhanced the adjacent sac1 gene expression. We found that sac1 is important to A. fumigatus, and the upregulated sac1 caused elevated resistance of A. fumigatus to azoles. Finally, we showed that Afut4_977-sac1 has an evolution pattern similar to that of the whole genome of azole-resistant strains due to azoles; phylogenetic analysis of both the whole genome and Afut4_977-sac1 suggests that the insertion of Afut4_977-sac1 might have preceded the emergence of azole-resistant strains. Taking these data together, we found that the Afut4_977-sac1 LTR-RT might be involved in the regulation of azole resistance of A. fumigatus by upregulating its nearby sac1 gene.

A comparative analysis of L1 retrotransposition activities in human genomes suggests an ongoing increase in L1 number despite an evolutionary trend towards lower activity

Background

LINE-1 (Long Interspersed Nuclear Elements, L1) retrotransposons are the only autonomously active transposable elements in the human genome. The evolution of L1 retrotransposition rates and its implications for L1 dynamics are poorly understood. Retrotransposition rates are commonly measured in cell culture-based assays, but it is unclear how well these measurements provide insight into L1 population dynamics. This study applied comparative methods to estimate parameters for the evolution of retrotransposition rates, and infer L1 dynamics from these estimates.

Results

Our results show that the rates at which new L1s emerge in the human population correlate positively to cell-culture based retrotransposition activities, that there is an evolutionary trend towards lower retrotransposition activity, and that this evolutionary trend is not sufficient to counter-balance the increase in active L1s resulting from continuing retrotransposition.

Conclusions

Together, these findings support a model of the population-level L1 retrotransposition dynamics that is consistent with prior expectations and indicate the remaining gaps in the understanding of L1 dynamics in human genomes.

A Genomic Survey of Mayetiola destructor Mobilome Provides New Insights into the Evolutionary History of Transposable Elements in the Cecidomyiid Midges

The availability of the Whole-Genome Sequence of the wheat pest Mayetiola destructor offers the opportunity to investigate the Transposable Elements (TEs) content and their relationship with the genes involved in the insect virulence. In this study, de novo annotation carried out using REPET pipeline showed that TEs occupy approximately 16% of the genome and are represented by 1038 lineages. Class II elements were the most frequent and most TEs were inactive due to the deletions they have accumulated. The analyses of TEs ages revealed a first burst at 20% of divergence from present that mobilized many TE families including mostly Tc1/mariner and Gypsy superfamilies and a second burst at 2% of divergence, which involved mainly the class II elements suggesting new TEs invasions. Additionally, 86 TEs insertions involving recently transposed elements were identified. Among them, several MITEs and Gypsy retrotransposons were inserted in the vicinity of SSGP and chemosensory genes. The findings represent a valuable resource for more in-depth investigation of the TE impact onto M. destructor genome and their possible influence on the expression of the virulence and chemosensory genes and consequently the behavior of this pest towards its host plants.

Transposable elements as new players in neurodegenerative diseases

Neurodegenerative diseases (NDs), including the most prevalent Alzheimer's disease and Parkinson disease, share common pathological features. Despite decades of gene-centric approaches, the molecular mechanisms underlying these diseases remain widely elusive. In recent years, transposable elements (TEs), long considered 'junk' DNA, have gained growing interest as pathogenic players in NDs. Age is the major risk factor for most NDs, and several repressive mechanisms of TEs, such as heterochromatinization, fail with age. Indeed, heterochromatin relaxation leading to TE derepression has been reported in various models of neurodegeneration and NDs. There is also evidence that certain pathogenic proteins involved in NDs (e.g., tau, TDP-43) may control the expression of TEs. The deleterious consequences of TE activation are not well known but they could include DNA damage and genomic instability, altered host gene expression, and/or neuroinflammation, which are common hallmarks of neurodegeneration and aging. TEs might thus represent an overlooked pathogenic culprit for both brain aging and neurodegeneration. Certain pathological effects of TEs might be prevented by inhibiting their activity, pointing to TEs as novel targets for neuroprotection.

ADAR1 interaction with Z-RNA promotes editing of endogenous double-stranded RNA and prevents MDA5-dependent immune activation

Loss of function of adenosine deaminase acting on double-stranded RNA (dsRNA)-1 (ADAR1) causes the severe autoinflammatory disease Aicardi-Goutières syndrome (AGS). ADAR1 converts adenosines into inosines within dsRNA. This process called A-to-I editing masks self-dsRNA from detection by the antiviral dsRNA sensor MDA5. ADAR1 binds to dsRNA in both the canonical A-form and the poorly defined Z conformation (Z-RNA). Mutations in the Z-RNA-binding Zα domain of ADAR1 are common in patients with AGS. How loss of ADAR1/Z-RNA interaction contributes to disease development is unknown. We demonstrate that abrogated binding of ADAR1 to Z-RNA leads to reduced A-to-I editing of dsRNA structures formed by base pairing of inversely oriented short interspersed nuclear elements. Preventing ADAR1 binding to Z-RNA triggers an MDA5/MAVS-mediated type I interferon response and leads to the development of lethal autoinflammation in mice. This shows that the interaction between ADAR1 and Z-RNA restricts sensing of self-dsRNA and prevents AGS development.

Driving to Safety: CRISPR-Based Genetic Approaches to Reducing Antibiotic Resistance

Bacterial resistance to antibiotics has reached critical levels, skyrocketing in hospitals and the environment and posing a major threat to global public health. The complex and challenging problem of reducing antibiotic resistance (AR) requires a network of both societal and science-based solutions to preserve the most lifesaving pharmaceutical intervention known to medicine. In addition to developing new classes of antibiotics, it is essential to safeguard the clinical efficacy of existing drugs. In this review, we examine the potential application of novel CRISPR-based genetic approaches to reducing AR in both environmental and clinical settings and prolonging the utility of vital antibiotics.

Population genomics of transposable element activation in the highly repressive genome of an agricultural pathogen

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

Background

Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation.

Results

We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating.

Conclusions

Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01081-6.

Role of long interspersed nuclear element-1 in the regulation of chromatin landscapes and genome dynamics

LINE-1 retrotransposon, the most active mobile element of the human genome, is subject to tight regulatory control. Stressful environments and disease modify the recruitment of regulatory proteins leading to unregulated activation of LINE-1. The activation of LINE-1 influences genome dynamics through altered chromatin landscapes, insertion mutations, deletions, and modulation of cellular plasticity. To date, LINE-1 retrotransposition has been linked to various cancer types and may in fact underwrite the genetic basis of various other forms of chronic human illness. The occurrence of LINE-1 polymorphisms in the human population may define inter-individual differences in susceptibility to disease. This review is written in honor of Dr Peter Stambrook, a friend and colleague who carried out highly impactful cancer research over many years of professional practice. Dr Stambrook devoted considerable energy to helping others live up to their full potential and to navigate the complexities of professional life. He was an inspirational leader, a strong advocate, a kind mentor, a vocal supporter and cheerleader, and yes, a hard critic and tough friend when needed. His passionate stand on issues, his witty sense of humor, and his love for humanity have left a huge mark in our lives. We hope that that the knowledge summarized here will advance our understanding of the role of LINE-1 in cancer biology and expedite the development of innovative cancer diagnostics and treatments in the ways that Dr Stambrook himself had so passionately envisioned.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

Elevated retrotransposon activity and genomic instability in primed pluripotent stem cells

BACKGROUND

Naïve and primed pluripotent stem cells (PSCs) represent two different pluripotent states. Primed PSCs following in vitro culture exhibit lower developmental potency as evidenced by failure in germline chimera assays, unlike mouse naïve PSCs. However, the molecular mechanisms underlying the lower developmental competency of primed PSCs remain elusive.

RESULTS

We examine the regulation of telomere maintenance, retrotransposon activity, and genomic stability of primed PSCs and compare them with naïve PSCs. Surprisingly, primed PSCs only minimally maintain telomeres and show fragile telomeres, associated with declined DNA recombination and repair activity, in contrast to naïve PSCs that robustly elongate telomeres. Also, we identify LINE1 family integrant L1Md_T as naïve-specific retrotransposon and ERVK family integrant IAPEz to define primed PSCs, and their transcription is differentially regulated by heterochromatic histones and Dnmt3b. Notably, genomic instability of primed PSCs is increased, in association with aberrant retrotransposon activity.

CONCLUSIONS

Our data suggest that fragile telomere, retrotransposon-associated genomic instability, and declined DNA recombination repair, together with reduced function of cell cycle and mitochondria, increased apoptosis, and differentiation properties may link to compromised developmental potency of primed PSCs, noticeably distinguishable from naïve PSCs.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

SINE Insertion in the Intron of Pig GHR May Decrease Its Expression by Acting as a Repressor

Simple Summary

GH/IGF axis genes play a central role in the regulation of skeletal accretion during development and growth, and thus represent candidate genes for growth traits. Retrotransposon insertion polymorphisms are major contributors to structural variations. They tend to generate large effect mutations resulting in variations in target gene activity and phenotype due to the fact that they carry functional elements, such as enhancers, insulators, or promoters. In the present study, RIPs in four GH/IGF axis genes (GH, GHR, IGF1, and IGF1R) were investigated by comparative genomics and PCR. Four RIPs in the GHR gene and one RIP in the IGF1 gene were identified. Further analysis revealed that one RIP in the first intron of GHR might play a role in the regulation of GHR expression by acting as a repressor. These findings contribute to the understanding of the role of RIPs in the genetic variation of GH/IGF axis genes and phenotypic variation in pigs.

Abstract

The genetic diversity of the GH/IGF axis genes and their association with the variation of gene expression and phenotypic traits, principally represented by SNPs, have been extensively reported. Nevertheless, the impact of retrotransposon insertion polymorphisms (RIPs) on the GH/IGF axis gene activity has not been reported. In the present study, bioinformatic prediction and PCR verification were performed to screen RIPs in four GH/IGF axis genes (GH, GHR, IGF1 and IGF1R). In total, five RIPs, including one SINE RIP in intron 3 of IGF1, one L1 RIP in intron 7 of GHR, and three SINE RIPs in intron 1, intron 5 and intron 9 of GHR, were confirmed by PCR, displaying polymorphisms in diverse breeds. Dual luciferase reporter assay revealed that the SINE insertion in intron 1 of GHR significantly repressed the GHR promoter activity in PK15, Hela, C2C12 and 3T3-L1 cells. Furthermore, qPCR results confirmed that this SINE insertion was associated with a decreased expression of GHR in the leg muscle and longissimus dorsi, indicating that it may act as a repressor involved in the regulation of GHR expression. In summary, our data revealed that RIPs contribute to the genetic variation of GH/IGF axis genes, whereby one SINE RIP in the intron 1 of GHR may decrease the expression of GHR by acting as a repressor.

A benchmark and an algorithm for detecting germline transposon insertions and measuring de novo transposon insertion frequencies

Transposons are genomic parasites, and their new insertions can cause instability and spur the evolution of their host genomes. Rapid accumulation of short-read whole-genome sequencing data provides a great opportunity for studying new transposon insertions and their impacts on the host genome. Although many algorithms are available for detecting transposon insertions, the task remains challenging and existing tools are not designed for identifying de novo insertions. Here, we present a new benchmark fly dataset based on PacBio long-read sequencing and a new method TEMP2 for detecting germline insertions and measuring de novo ‘singleton’ insertion frequencies in eukaryotic genomes. TEMP2 achieves high sensitivity and precision for detecting germline insertions when compared with existing tools using both simulated data in fly and experimental data in fly and human. Furthermore, TEMP2 can accurately assess the frequencies of de novo transposon insertions even with high levels of chimeric reads in simulated datasets; such chimeric reads often occur during the construction of short-read sequencing libraries. By applying TEMP2 to published data on hybrid dysgenic flies inflicted by de-repressed P-elements, we confirmed the continuous new insertions of P-elements in dysgenic offspring before they regain piRNAs for P-element repression. TEMP2 is freely available at Github: https://github.com/weng-lab/TEMP2.

SARS-CoV-2-Host Chimeric RNA-Sequencing Reads Do Not Necessarily Arise From Virus Integration Into the Host DNA

The human genome bears evidence of extensive invasion by retroviruses and other retroelements, as well as by diverse RNA and DNA viruses. High frequency of somatic integration of the RNA virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into the DNA of infected cells was recently suggested, based on a number of observations. One key observation was the presence of chimeric RNA-sequencing (RNA-seq) reads between SARS-CoV-2 RNA and RNA transcribed from human host DNA. Here, we examined the possible origin specifically of human-SARS-CoV-2 chimeric reads in RNA-seq libraries and provide alternative explanations for their origin. Chimeric reads were frequently detected also between SARS-CoV-2 RNA and RNA transcribed from mitochondrial DNA or episomal adenoviral DNA present in transfected cell lines, which was unlikely the result of SARS-CoV-2 integration. Furthermore, chimeric reads between SARS-CoV-2 RNA and RNA transcribed from nuclear DNA were highly enriched for host exonic, rather than intronic or intergenic sequences and often involved the same, highly expressed host genes. Although these findings do not rule out SARS-CoV-2 somatic integration, they nevertheless suggest that human-SARS-CoV-2 chimeric reads found in RNA-seq data may arise during library preparation and do not necessarily signify SARS-CoV-2 reverse transcription, integration in to host DNA and further transcription.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

Retrotransposons as Drivers of Mammalian Brain Evolution

Retrotransposons, a large and diverse class of transposable elements that are still active in humans, represent a remarkable force of genomic innovation underlying mammalian evolution. Among the features distinguishing mammals from all other vertebrates, the presence of a neocortex with a peculiar neuronal organization, composition and connectivity is perhaps the one that, by affecting the cognitive abilities of mammals, contributed mostly to their evolutionary success. Among mammals, hominids and especially humans display an extraordinarily expanded cortical volume, an enrichment of the repertoire of neural cell types and more elaborate patterns of neuronal connectivity. Retrotransposon-derived sequences have recently been implicated in multiple layers of gene regulation in the brain, from transcriptional and post-transcriptional control to both local and large-scale three-dimensional chromatin organization. Accordingly, an increasing variety of neurodevelopmental and neurodegenerative conditions are being recognized to be associated with retrotransposon dysregulation. We review here a large body of recent studies lending support to the idea that retrotransposon-dependent evolutionary novelties were crucial for the emergence of mammalian, primate and human peculiarities of brain morphology and function.

Phase separation of the LINE-1 ORF1 protein is mediated by the N-terminus and coiled-coil domain

Long interspersed nuclear element-1 (L1) is a retrotransposable element that autonomously replicates in the human genome, resulting in DNA damage and genomic instability. Activation of L1 in senescent cells triggers a type I interferon response and age-associated inflammation. Two open reading frames encode an ORF1 protein functioning as messenger RNA chaperone and an ORF2 protein providing catalytic activities necessary for retrotransposition. No function has been identified for the conserved, disordered N-terminal region of ORF1. Using microscopy and NMR spectroscopy, we demonstrate that ORF1 forms liquid droplets in vitro in a salt-dependent manner and that interactions between its N-terminal region and coiled-coil domain are necessary for phase separation. Mutations disrupting blocks of charged residues within the N-terminus impair phase separation, whereas some mutations within the coiled-coil domain enhance phase separation. Demixing of the L1 particle from the cytosol may provide a mechanism to protect the L1 transcript from degradation.

Characterization of an active LINE-1 in the naked mole-rat genome

Naked mole-rats (NMRs, Heterocephalus glaber) are the longest-living rodent species. A reason for their long lifespan is pronounced cancer resistance. Therefore, researchers believe that NMRs have unknown secrets of cancer resistance and seek to find them. Here, to reveal the secrets, we noticed a retrotransposon, long interspersed nuclear element 1 (L1). L1s can amplify themselves and are considered endogenous oncogenic mutagens. Since the NMR genome contains fewer L1-derived sequences than other mammalian genomes, we reasoned that the retrotransposition activity of L1s in the NMR genome is lower than those in other mammalian genomes. In this study, we successfully cloned an intact L1 from the NMR genome and named it NMR-L1. An L1 retrotransposition assay using the NMR-L1 reporter revealed that NMR-L1 was active retrotransposon, but its activity was lower than that of human and mouse L1s. Despite lower retrotrasposition activity, NMR-L1 was still capable of inducing cell senescence, a tumor-protective system. NMR-L1 required the 3′ untranslated region (UTR) for retrotransposition, suggesting that NMR-L1 is a stringent-type of L1. We also confirmed the 5′ UTR promoter activity of NMR-L1. Finally, we identified the G-quadruplex structure of the 3′ UTR, which modulated the retrotransposition activity of NMR-L1. Taken together, the data indicate that NMR-L1 retrotranspose less efficiently, which may contribute to the cancer resistance of NMRs.

Mechanisms of Action of Hypomethylating Agents: Endogenous Retroelements at the Epicenter

Abnormal DNA methylation patterns are thought to drive the pathobiology of high-risk myelodysplastic syndromes (HR-MDS) and acute myeloid leukemia (AML). Sixteen years after their initial approval, the hypomethylating agents (HMAs), 5-azacytidine (AZA) and 5-aza-2′-deoxycytidine, remain the mainstay of treatment for HR-MDS and AML. However, a connection of the hypomethylating or additional effects of HMAs with clinical responses remains yet to be shown, and the mode of action of HMAs remains obscure. Given the relatively short-lived responses and the inevitable development of resistance in HMAs, a thorough understanding of the antineoplastic mechanisms employed by HMAs holds critical importance. Recent data in cancer cell lines demonstrate that reactivation of endogenous retroelements (EREs) and induction of a cell-intrinsic antiviral response triggered by RNA neotranscripts may underlie the antitumor activity of HMAs. However, data on primary CD34⁺ cells derived from patients with HR-MDS failed to confirm a link between HMA-mediated ERE modulation and clinical response. Though difficult to reconcile the apparent discrepancy, it is possible that HMAs mediate their effects in more advanced levels of differentiation where cells become responsive to interferon, whereas, inter-individual variations in the process of RNA editing and, in particular, in the ADAR1/OAS/RNase L pathway may also confound the associations of clinical response with the induction of viral mimicry. Further ex vivo studies along with clinical correlations in well-annotated patient cohorts are warranted to decipher the role of ERE derepression in the antineoplastic mechanisms of HMAs.