The impact of retrotransposons on human genome evolution

Mechanisms of disease-associated SINE-VNTR-Alus

SINE-VNTR-Alus (SVAs) are the youngest retrotransposon family in the human genome. Their ongoing mobilization has generated genetic variation within the human population. At least 24 insertions to date, detailed in this review, have been associated with disease. The predominant mechanisms through which this occurs are alterations to normal splicing patterns, exonic insertions causing loss-of-function mutations, and large genomic deletions. Dissecting the functional impact of these SVAs and the mechanism through which they cause disease provides insight into the consequences of their presence in the genome and how these elements could influence phenotypes. Many of these disease-associated SVAs have been difficult to characterize and would not have been identified through routine analyses. However, the number identified has increased in recent years as DNA and RNA sequencing data became more widely available. Therefore, as the search for complex structural variation in disease continues, it is likely to yield further disease-causing SVA insertions.

A hypothesis: Retrotransposons as a relay of epigenetic marks in intergenerational epigenetic inheritance

Epigenetic marks in gametes, which both respond to the parental environmental factors and shape offspring phenotypes, are usually positioned to mediate intergenerational or transgenerational epigenetic inheritance. Nonetheless, the mechanisms through which gametic epigenetic signatures encode parental acquired phenotypes, and further initiate a cascade of molecular events to affect offspring phenotypes during early embryonic development, remain unclear. Retrotransposons are mobile DNA elements that could resist to genomic epigenetic reprogramming at specific loci and rewire the core regulatory networks of embryogenesis. Increasing evidences show that retrotransposons in the embryonic genome could interact with gametic epigenetic marks, which provides a tentative possibility that retrotransposons may serve as a relay of gametic epigenetic marks to transmit parental acquired traits. Here, we summarize the recent progress in exploring the crosstalk between gametic epigenetic marks and retrotransposons, and the regulation of gene expression and early embryonic development by retrotransposons. Accordingly, deciphering the mystery of interactions between gametic epigenetic marks and retrotransposons during early embryonic development will provide valuable insights into the intergenerational or transgenerational transmission of acquired traits.

The genomic architecture of the passerine MHC region: high repeat content and contrasting evolutionary histories of single copy and tandemly duplicated MHC genes

The Major Histocompatibility Complex (MHC) is of central importance to the immune system, and an optimal MHC diversity is believed to maximize pathogen elimination. Birds show substantial variation in MHC diversity, ranging from few genes in most bird orders to very many genes in passerines. Our understanding of the evolutionary trajectories of the MHC in passerines is hampered by lack of data on genomic organization. Therefore, we assemble and annotate the MHC genomic region of the great reed warbler (Acrocephalus arundinaceus), using long-read sequencing and optical mapping. The MHC region is large (>5.5Mb), characterized by structural changes compared to hitherto investigated bird orders and shows higher repeat content than the genome average. These features were supported by analyses in three additional passerines. MHC genes in passerines are found in two different chromosomal arrangements, either as single copy MHC genes located among non-MHC genes, or as tandemly duplicated tightly linked MHC genes. Some single copy MHC genes are old and putative orthologs among species. In contrast tandemly duplicated MHC genes are monophyletic within species and have evolved by simultaneous gene duplication of several MHC genes. Structural differences in the MHC genomic region among bird orders seem substantial compared to mammals and have possibly been fuelled by clade-specific immune system adaptations. Our study provides methodological guidance in characterizing complex genomic regions, constitutes a resource for MHC research in birds, and calls for a revision of the general belief that avian MHC has a conserved gene order and small size compared to mammals.

Cu(I) Binding to Designed Proteins Reveals a Putative Copper Binding Site of the Human Line1 Retrotransposon Protein ORF1p

Long interspersed nuclear elements-1 (L1) are autonomous retrotransposons that encode two proteins in different open reading frames (ORF1 and ORF2). The ORF1p, which may be an RNA binding and chaperone protein, contains a three-stranded coiled coil (3SCC) domain that facilitates the formation of the biologically active homotrimer. This 3SCC domain is composed of seven amino acid (heptad) repeats as found in native and designed peptides and a stammer that modifies the helical structure. Cysteine residues occur at three hydrophobic positions (2 a and 1 d sites) within this domain. We recently showed that the cysteine layers in ORF1p and model de novo designed peptides bind the toxic metalloid lead(II) with high affinities, a feature that had not been previously recognized. However, there is little understanding of how essential metal ions might interact with this metal binding domain. We have, therefore, investigated the copper(I) binding properties of analogous de novo designed 3SCCs that contain cysteine layers within the hydrophobic core. The results from UV-visible and X-ray absorption spectroscopy show that these designed peptides bind Cu(I) with high affinity in a pH-dependent manner. At pH 9, monomeric trigonal planar Cu(I)S3 centers are formed with 1 equiv of metal, while dinuclear centers form with a second equivalent of metal. At physiologic pH conditions, the dinuclear center forms cooperatively. These data suggest that ORF1p is capable of binding two copper ions to its tris(cysteine) layers. This has major implications for ORF1p coiled coil domain stability and dynamics, ultimately potentially impacting the resulting biological activity.

A hidden layer of structural variation in transposable elements reveals potential genetic modifiers in human disease-risk loci

Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease-associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease-associated risk loci and in the BCKDK Parkinson's disease-associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

Chromatin modifier HUSH co-operates with RNA decay factor NEXT to restrict transposable element expression

Transposable elements (TEs) are widespread genetic parasites known to be kept under tight transcriptional control. Here, we describe a functional connection between the mouse-orthologous “nuclear exosome targeting” (NEXT) and “human silencing hub” (HUSH) complexes, involved in nuclear RNA decay and the epigenetic silencing of TEs, respectively. Knocking out the NEXT component ZCCHC8 in embryonic stem cells results in elevated TE RNA levels. We identify a physical interaction between ZCCHC8 and the MPP8 protein of HUSH and establish that HUSH recruits NEXT to chromatin at MPP8-bound TE loci. However, while NEXT and HUSH both dampen TE RNA expression, their activities predominantly affect shorter non-polyadenylated and full-length polyadenylated transcripts, respectively. Indeed, our data suggest that the repressive action of HUSH promotes a condition favoring NEXT RNA decay activity. In this way, transcriptional and post-transcriptional machineries synergize to suppress the genotoxic potential of TE RNAs.

Garland et al. report a physical and functional connection between the NEXT complex, involved in RNA decay, and the HUSH complex, involved in chromatin regulation. Together, NEXT and HUSH cooperate to control transposable element (TE) RNA expression in embryonic stem cells, suppressing pA⁻ and pA⁺ transcripts, respectively.

Transposable Elements: Major Players in Shaping Genomic and Evolutionary Patterns

Transposable elements (TEs) are ubiquitous genetic elements, able to jump from one location of the genome to another, in all organisms. For this reason, on the one hand, TEs can induce deleterious mutations, causing dysfunction, disease and even lethality in individuals. On the other hand, TEs can increase genetic variability, making populations better equipped to respond adaptively to environmental change. To counteract the deleterious effects of TEs, organisms have evolved strategies to avoid their activation. However, their mobilization does occur. Usually, TEs are maintained silent through several mechanisms, but they can be reactivated during certain developmental windows. Moreover, TEs can become de-repressed because of drastic changes in the external environment. Here, we describe the ‘double life’ of TEs, being both ‘parasites’ and ‘symbionts’ of the genome. We also argue that the transposition of TEs contributes to two important evolutionary processes: the temporal dynamic of evolution and the induction of genetic variability. Finally, we discuss how the interplay between two TE-dependent phenomena, insertional mutagenesis and epigenetic plasticity, plays a role in the process of evolution.

SEEMLIS: a flexible semi-automated method for enrichment of methylated DNA from low-input samples

Background

DNA methylation alterations have emerged as hallmarks of cancer and have been proposed as screening, prognostic, and predictive biomarkers. Traditional approaches for methylation analysis have relied on bisulfite conversion of DNA, which can damage DNA and is not suitable for targeted gene analysis in low-input samples. Here, we have adapted methyl-CpG-binding domain protein 2 (MBD2)-based DNA enrichment for use on a semi-automated exclusion-based sample preparation (ESP) platform for robust and scalable enrichment of methylated DNA from low-input samples, called SEEMLIS.

Results

We show that combining methylation-sensitive enzyme digestion with ESP-based MBD2 enrichment allows for single gene analysis with high sensitivity for GSTP1 in highly impure, heterogenous samples. We also show that ESP-based MBD2 enrichment coupled with targeted pre-amplification allows for analysis of multiple genes with sensitivities approaching the single cell level in pure samples for GSTP1 and RASSF1 and sensitivity down to 14 cells for these genes in highly impure samples. Finally, we demonstrate the potential clinical utility of SEEMLIS by successful detection of methylated gene signatures in circulating tumor cells (CTCs) from patients with prostate cancer with varying CTC number and sample purity.

Conclusions

SEEMLIS is a robust assay for targeted DNA methylation analysis in low-input samples, with flexibility at multiple steps. We demonstrate the feasibility of this assay to analyze DNA methylation in prostate cancer cells using CTCs from patients with prostate cancer as a real-world example of a low-input analyte of clinical importance. In summary, this novel assay provides a platform for determining methylation signatures in rare cell populations with broad implications for research as well as clinical applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01252-4.

Genetic mutation spectrum of pantothenate kinase-associated neurodegeneration expanded by breakpoint sequencing in pantothenate kinase 2 gene

BACKGROUND

Neurodegeneration with brain iron accumulation describes a group of rare heterogeneous progressive neurodegenerative disorders characterized by excessive iron accumulation in the basal ganglia region. Pantothenate kinase-associated neurodegeneration (PKAN) is a major form of this disease.

RESULTS

A total of 7 unrelated patients were diagnosed with PKAN in a single tertiary center from August 2009 to February 2018. Ten variants in PANK2 including three novel sequence variants and one large exonic deletion were detected. Sequencing of the breakpoint was performed to predict the mechanism of large deletion and AluSx3 and AluSz6 were found with approximately 97.3% sequence homology.

CONCLUSION

The findings support the disease-causing role of PANK2 and indicate the possibility that exonic deletion of PANK2 found in PKAN is mediated through Alu-mediated homologous recombination.

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor–CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell–cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Characterizing mobile element insertions in 5675 genomes

Mobile element insertions (MEIs) are a major class of structural variants (SVs) and have been linked to many human genetic disorders, including hemophilia, neurofibromatosis, and various cancers. However, human MEI resources from large-scale genome sequencing are still lacking compared to those for SNPs and SVs. Here, we report a comprehensive map of 36 699 non-reference MEIs constructed from 5675 genomes, comprising 2998 Chinese samples (∼26.2×, NyuWa) and 2677 samples from the 1000 Genomes Project (∼7.4×, 1KGP). We discovered that LINE-1 insertions were highly enriched in centromere regions, implying the role of chromosome context in retroelement insertion. After functional annotation, we estimated that MEIs are responsible for about 9.3% of all protein-truncating events per genome. Finally, we built a companion database named HMEID for public use. This resource represents the latest and largest genomewide study on MEIs and will have broad utility for exploration of human MEI findings.

An Alu insertion map of the Indian population: identification and analysis in 1021 genomes of the IndiGen project

Actively retrotransposing primate-specific Alu repeats display insertion-deletion (InDel) polymorphism through their insertion at new loci. In the global datasets, Indian populations remain under-represented and so do their Alu InDels. Here, we report the genomic landscape of Alu InDels from the recently released 1021 Indian Genomes (IndiGen) (available at https://clingen.igib.res.in/indigen). We identified 9239 polymorphic Alu insertions that include private (3831), rare (3974) and common (1434) insertions with an average of 770 insertions per individual. We achieved an 89% PCR validation of the predicted genotypes in 94 samples tested. About 60% of identified InDels are unique to IndiGen when compared to other global datasets; 23% of sites were shared with both SGDP and HGSVC; among these, 58% (1289 sites) were common polymorphisms in IndiGen. The insertions not only show a bias for genic regions, with a preference for introns but also for the associated genes showing enrichment for processes like cell morphogenesis and neurogenesis (P-value < 0.05). Approximately, 60% of InDels mapped to genes present in the OMIM database. Finally, we show that 558 InDels can serve as ancestry informative markers to segregate global populations. This study provides a valuable resource for baseline Alu InDels that would be useful in population genomics.

AnnoSINE: a short interspersed nuclear elements annotation tool for plant genomes

Short interspersed nuclear elements (SINEs) are a widespread type of small transposable element (TE). With increasing evidence for their impact on gene function and genome evolution in plants, accurate genome-scale SINE annotation becomes a fundamental step for studying the regulatory roles of SINEs and their relationship with other components in the genomes. Despite the overall promising progress made in TE annotation, SINE annotation remains a major challenge. Unlike some other TEs, SINEs are short and heterogeneous, and they usually lack well-conserved sequence or structural features. Thus, current SINE annotation tools have either low sensitivity or high false discovery rates. Given the demand and challenges, we aimed to provide a more accurate and efficient SINE annotation tool for plant genomes. The pipeline starts with maximizing the pool of SINE candidates via profile hidden Markov model-based homology search and de novo SINE search using structural features. Then, it excludes the false positives by integrating all known features of SINEs and the features of other types of TEs that can often be misannotated as SINEs. As a result, the pipeline substantially improves the tradeoff between sensitivity and accuracy, with both values close to or over 90%. We tested our tool in Arabidopsis thaliana and rice (Oryza sativa), and the results show that our tool competes favorably against existing SINE annotation tools. The simplicity and effectiveness of this tool would potentially be useful for generating more accurate SINE annotations for other plant species. The pipeline is freely available at https://github.com/yangli557/AnnoSINE.

Deletion and replacement of long genomic sequences using prime editing

Genomic insertions, duplications and insertion/deletions (indels), which account for ~14% of human pathogenic mutations, cannot be accurately or efficiently corrected by current gene-editing methods, especially those that involve larger alterations (>100 base pairs (bp)). Here, we optimize prime editing (PE) tools for creating precise genomic deletions and direct the replacement of a genomic fragment ranging from ~1 kilobases (kb) to ~10 kb with a desired sequence (up to 60 bp) in the absence of an exogenous DNA template. By conjugating Cas9 nuclease to reverse transcriptase (PE-Cas9) and combining it with two PE guide RNAs (pegRNAs) targeting complementary DNA strands, we achieve precise and specific deletion and repair of target sequences via using this PE-Cas9-based deletion and repair (PEDAR) method. PEDAR outperformed other genome-editing methods in a reporter system and at endogenous loci, efficiently creating large and precise genomic alterations. In a mouse model of tyrosinemia, PEDAR removed a 1.38-kb pathogenic insertion within the Fah gene and precisely repaired the deletion junction to restore FAH expression in liver.

MxB inhibits long interspersed element type 1 retrotransposition

Long interspersed element type 1 (LINE-1, also L1 for short) is the only autonomously transposable element in the human genome. Its insertion into a new genomic site may disrupt the function of genes, potentially causing genetic diseases. Cells have thus evolved a battery of mechanisms to tightly control LINE-1 activity. Here, we report that a cellular antiviral protein, myxovirus resistance protein B (MxB), restricts the mobilization of LINE-1. This function of MxB requires the nuclear localization signal located at its N-terminus, its GTPase activity and its ability to form oligomers. We further found that MxB associates with LINE-1 protein ORF1p and promotes sequestration of ORF1p to G3BP1-containing cytoplasmic granules. Since knockdown of stress granule marker proteins G3BP1 or TIA1 abolishes MxB inhibition of LINE-1, we conclude that MxB engages stress granule components to effectively sequester LINE-1 proteins within the cytoplasmic granules, thus hindering LINE-1 from accessing the nucleus to complete retrotransposition. Thus, MxB protein provides one mechanism for cells to control the mobility of retroelements.

Retrotransposons occupy more than 40% of human genome, and have co-evolved with humans for millions of years. Long interspersed element type 1 (LINE-1, or L1) is the only retrotransposon that is able to jump to a new locus. LINE-1 retrotransposition causes genome instability, and is associated with genetic diseases including autoimmune diseases and cancer. To suppress this genome toxicity caused by LINE-1, humans have developed multi-layered mechanisms to control LINE-1 activity. MxB has been previously shown to inhibit LINE-1 mobility, thus contributing to host restriction of LINE-1. Here, we further demonstrate that MxB effectively restricts LINE-1 retrotransposition by sequestering LINE-1 ribonucleoprotein (RNP) within the cytoplasmic stress granules, thus guards genome stability. Hence our data attribute the restriction function of MxB to sequestering LINE-1 RNP to stress granules.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of ‘transitional’ L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Transposons: Unexpected players in cancer

Transposons are repetitive DNA sequences encompassing about half of the human genome. They play a vital role in genome stability maintenance and contribute to genomic diversity and evolution. Their activity is regulated by various mechanisms considering the deleterious effects of these mobile elements. Various genetic risk factors and environmental stress conditions affect the regulatory pathways causing alteration of transposon expression. Our knowledge of the biological role of transposons is limited especially in various types of cancers. Retrotransposons of different types (LTR-retrotransposons, LINEs and SINEs) regulate a plethora of genes that have a role in cell reprogramming, tumor suppression, cell cycle, apoptosis, cell adhesion and migration, and DNA repair. The regulatory mechanisms of transposons, their deregulation and different mechanisms underlying transposon-mediated carcinogenesis in humans focusing on the three most prevalent types, lung, breast and colorectal cancers, were reviewed. The modes of regulation employed include alternative splicing, deletion, insertion, duplication in genes and promoters resulting in upregulation, downregulation or silencing of genes.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

The A to I editing landscape in melanoma and its relation to clinical outcome

RNA editing refers to non-transient RNA modifications that occur after transcription and prior to translation by the ribosomes. RNA editing is more widespread in cancer cells than in non-transformed cells and is associated with tumorigenesis of various cancer tissues. However, RNA editing can also generate neo-antigens that expose tumour cells to host immunosurveillance. Global RNA editing in melanoma and its relevance to clinical outcome currently remain poorly characterized. The present study compared RNA editing as well as gene expression in tumour cell lines from melanoma patients of short or long metastasis-free survival, patients relapsing or not after immuno- and targeted therapy and tumours harbouring BRAF or NRAS mutations. Overall, our results showed that NTRK gene expression can be a marker of resistance to BRAF and MEK inhibition and gives some insights of candidate genes as potential biomarkers. In addition, this study revealed an increase in Adenosine-to-Inosine editing in Alu regions and in non-repetitive regions, including the hyperediting of the MOK and DZIP3 genes in relapsed tumour samples during targeted therapy and of the ZBTB11 gene in NRAS mutated melanoma cells. Therefore, RNA editing could be a promising tool for identifying predictive markers, tumour neoantigens and targetable pathways that could help in preventing relapses during immuno- or targeted therapies.

Emerging roles and potential application of PIWI-interacting RNA in urological tumors

The piRNA (PIWI-interacting RNA) is P-Element induced wimpy testis (PIWI)-interacting RNA which is a small molecule, non-coding RNA with a length of 24-32nt. It was originally found in germ cells and is considered a regulator of germ cell function. It can interact with PIWI protein, a member of the Argonaute family, and play a role in the regulation of gene transcription and epigenetic silencing of transposable factors in the nucleus. More and more studies have shown that piRNAs are abnormally expressed in a variety of cancer tissues and patient fluids, and may become diagnostic tools, therapeutic targets, staging markers, and prognostic evaluation tools for cancer. This article reviews the recent research on piRNA and summarizes the structural characteristics, production mechanism, applications, and its role in urological tumors, to provide a reference value for piRNA to regulate urological tumors.

The Emerging Role of PIWI-Interacting RNAs (piRNAs) in Gastrointestinal Cancers: An Updated Perspective

Gastrointestinal (GI) cancers produce ~3.4 million related deaths worldwide, comprising 35% of all cancer-related deaths. The high mortality among GI cancers is due to late diagnosis, the presence of metastasis and drug resistance development. Additionally, current clinical markers do not adequately guide patient management, thereby new and more reliable biomarkers and therapeutic targets are still needed for these diseases. RNA-seq technology has allowed the discovery of new types of RNA transcripts including PIWI-interacting RNAs (piRNAs), which have particular characteristics that enable these molecules to act via diverse molecular mechanisms for regulating gene expression. Cumulative evidence has described the potential role of piRNAs in the development of several tumor types as a likely explanation for certain genomic abnormalities and signaling pathways' deregulations observed in cancer. In addition, these piRNAs might be also proposed as promising diagnostic or prognostic biomarkers or as potential therapeutic targets in malignancies. This review describes important topics about piRNAs including their molecular characteristics, biosynthesis processes, gene expression silencing mechanisms, and the manner in which these transcripts have been studied in samples and cell lines of GI cancers to elucidate their implications in these diseases. Moreover, this article discusses the potential clinical usefulness of piRNAs as biomarkers and therapeutic targets in GI cancers.

PRDM9-directed recombination hotspots depleted near meiotically transcribed genes

PRDM9 drives recombination hotspots in some mammals, including mice and apes. Non-functional orthologs of PRDM9 are present in a wide variety of vertebrates, but why it is functionally maintained in some lineages is not clear. One possible explanation is that PRDM9 plays a role in ensuring that meiosis is successful. During meiosis, available DNA may be a limiting resource given the tight packaging of chromosomes and could lead to competition between two key processes: meiotic transcription and recombination. Here we explore this potential competition and the role that PRDM9 might play in their interaction. Leveraging existing mouse genomic data, we use resampling schemes that simulate shuffled features along the genome and models that account for the rarity of features in the genome, to test if PRDM9 influences interactions between recombination hotspots and meiotic transcription in a whole genome framework. We also explored possible DNA sequence motifs associated to clusters of hotspots not tied to transcription or PRDM9. We find evidence of competition between meiotic transcription and recombination, with PRDM9 appearing to relocate recombination to avoid said conflict. We also find that retrotransposons may be playing a role in directing hotspots in the absence of other factors.

Complex Genetic Interactions between Piwi and HP1a in the Repression of Transposable Elements and Tissue-Specific Genes in the Ovarian Germline

Insertions of transposable elements (TEs) in eukaryotic genomes are usually associated with repressive chromatin, which spreads to neighbouring genomic sequences. In ovaries of Drosophila melanogaster, the Piwi-piRNA pathway plays a key role in the transcriptional silencing of TEs considered to be exerted mostly through the establishment of H3K9me3 histone marks recruiting Heterochromatin Protein 1a (HP1a). Here, using RNA-seq, we investigated the expression of TEs and the adjacent genomic regions upon Piwi and HP1a germline knockdowns sharing a similar genetic background. We found that the depletion of Piwi and HP1a led to the derepression of only partially overlapping TE sets. Several TEs were silenced predominantly by HP1a, whereas the upregulation of some other TEs was more pronounced upon Piwi knockdown and, surprisingly, was diminished upon a Piwi/HP1a double-knockdown. We revealed that HP1a loss influenced the expression of thousands of protein-coding genes mostly not adjacent to TE insertions and, in particular, downregulated a putative transcriptional factor required for TE activation. Nevertheless, our results indicate that Piwi and HP1a cooperatively exert repressive effects on the transcription of euchromatic loci flanking the insertions of some Piwi-regulated TEs. We suggest that this mechanism controls the silencing of a small set of TE-adjacent tissue-specific genes, preventing their inappropriate expression in ovaries.

RNA modifications as emerging therapeutic targets

The field of epitranscriptome, posttranscriptional modifications to RNAs, is still growing up and has presented substantial evidences for the role of RNA modifications in diseases. In terms of new drug development, RNA modifications have a great promise for therapy. For example, more than 170 type of modifications exist in various types of RNAs. Regulatory genes and their roles in critical biological process have been identified and they are associated with several diseases. Current data, for example, identification of inhibitors targeting RNA modifications regulatory genes, strongly support the idea that RNA modifications have potential as emerging therapeutic targets. Therefore, in this review, RNA modifications and regulatory genes were comprehensively documented in terms of drug development by summarizing the findings from previous studies. It was discussed how RNA modifications or regulatory genes can be targeted by altering molecular mechanisms. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > RNA Editing and Modification.

A hypothetical model of trans-acting R-loops-mediated promoter-enhancer interactions by Alu elements

Enhancers modulate gene expression by interacting with promoters. Models of enhancer-promoter interactions (EPIs) in the literature involve the activity of many components, including transcription factors and nucleic acid. However, the role that sequence similarity plays in EPIs remains largely unexplored. Herein, we report that Alu-derived sequences dominate sequence similarity between enhancers and promoters. After rejecting alternative DNA:DNA and DNA:RNA triplex models, we propose that enhancer-associated RNAs (eRNAs) may directly contact their targeted promoters by forming trans-acting R-loops at those Alu sequences. We show how the characteristic distribution of functional genomic data, such as RNA-DNA proximate ligation reads, binding of transcription factors, and RNA-binding proteins, all align with the Alu sequences of EPIs. We also show that these aligned Alu sequences may be subject to the constraint of coevolution, further implying the functional significance of these R-loop hybrids. Finally, our results imply that eRNA and Alu elements associate in a manner previously unrecognized in EPIs and the evolution of gene regulation networks in mammals.

Systems Biology to Understand and Regulate Human Retroviral Proinflammatory Response

The majority of human genome are non-coding genes. Recent research have revealed that about half of these genome sequences make up of transposable elements (TEs). A branch of these belong to the endogenous retroviruses (ERVs), which are germline viral infection that occurred over millions of years ago. They are generally harmless as evolutionary mutations have made them unable to produce viral agents and are mostly epigenetically silenced. Nevertheless, ERVs are able to express by still unknown mechanisms and recent evidences have shown links between ERVs and major proinflammatory diseases and cancers. The major challenge is to elucidate a detailed mechanistic understanding between them, so that novel therapeutic approaches can be explored. Here, we provide a brief overview of TEs, human ERVs and their links to microbiome, innate immune response, proinflammatory diseases and cancer. Finally, we recommend the employment of systems biology approaches for future HERV research.

Stochastic Effects in Retrotransposon Dynamics Revealed by Modeling under Competition for Cellular Resources

Transposons are genomic elements that can relocate within a host genome using a ‘cut’- or ‘copy-and-paste’ mechanism. They make up a significant part of many genomes, serve as a driving force for genome evolution, and are linked with Mendelian diseases and cancers. Interactions between two specific retrotransposon types, autonomous (e.g., LINE1/L1) and nonautonomous (e.g., Alu), may lead to fluctuations in the number of these transposons in the genome over multiple cell generations. We developed and examined a simple model of retrotransposon dynamics under conditions where transposon replication machinery competed for cellular resources: namely, free ribosomes and available energy (i.e., ATP molecules). Such competition is likely to occur in stress conditions that a malfunctioning cell may experience as a result of a malignant transformation. The modeling revealed that the number of actively replicating LINE1 and Alu elements in a cell decreases with the increasing competition for resources; however, stochastic effects interfere with this simple trend. We stochastically simulated the transposon dynamics in a cell population and showed that the population splits into pools with drastically different transposon behaviors. The early extinction of active Alu elements resulted in a larger number of LINE1 copies occurring in the first pool, as there was no competition between the two types of transposons in this pool. In the other pool, the competition process remained and the number of L1 copies was kept small. As the level of available resources reached a critical value, both types of dynamics demonstrated an increase in noise levels, and both the period and the amplitude of predator–prey oscillations rose in one of the cell pools. We hypothesized that the presented dynamical effects associated with the impact of the competition for cellular resources inflicted on the dynamics of retrotransposable elements could be used as a characteristic feature to assess a cell state, or to control the transposon activity.

Estimating DNA methylation potential energy landscapes from nanopore sequencing data

High-throughput third-generation nanopore sequencing devices have enormous potential for simultaneously observing epigenetic modifications in human cells over large regions of the genome. However, signals generated by these devices are subject to considerable noise that can lead to unsatisfactory detection performance and hamper downstream analysis. Here we develop a statistical method, CpelNano, for the quantification and analysis of 5mC methylation landscapes using nanopore data. CpelNano takes into account nanopore noise by means of a hidden Markov model (HMM) in which the true but unknown (“hidden”) methylation state is modeled through an Ising probability distribution that is consistent with methylation means and pairwise correlations, whereas nanopore current signals constitute the observed state. It then estimates the associated methylation potential energy function by employing the expectation-maximization (EM) algorithm and performs differential methylation analysis via permutation-based hypothesis testing. Using simulations and analysis of published data obtained from three human cell lines (GM12878, MCF-10A, and MDA-MB-231), we show that CpelNano can faithfully estimate DNA methylation potential energy landscapes, substantially improving current methods and leading to a powerful tool for the modeling and analysis of epigenetic landscapes using nanopore sequencing data.

DNA methylation profile in patients with indolent systemic mastocytosis

BACKGROUND

Mastocytosis is a clinically heterogeneous, usually acquired disease of the mast cells with a survival time that depends on the onset of the disease and ranges from skin-limited to systemic disease, including indolent and more aggressive variants. The crucial element in pathogenesis is the presence of oncogenic KIT somatic mutation D816V. Further epigenetic alterations are responsible for regulating the expression of genes. It is essential to identify indicators of disease progression, and the specific clinical picture to establish an appropriate therapeutic strategy.

OBJECTIVE

The aim of this study was to analyze the relation of mastocytosis symptoms and epigenetic changes, and to identify epigenetic predictors of the disease.

METHODS

Global DNA methylation profile analysis was performed in peripheral blood collected from 73 patients with indolent systemic mastocytosis (ISM) and 43 healthy adult volunteers. Levels of 5-methylcytosine (5-mC) and 5-hydroxymethylcytosine (5-hmC) were determined using an ELISA-based method, while the methylation of the Alu and LINE-1 repeats were assayed with the quantitative methylation-specific PCR technique. A questionnaire interview was conducted among the study participants to collect data on possible epigenetic modifiers. Additionally, the methylation profile was compared between three human mast cell lines: ROSA KIT D816V, ROSA KIT WT, and HMC-1.1 KIT V560G, in order to assess the association between KIT mutations and methylation profile.

RESULTS

A significantly lower level of DNA hydroxymethylation (5-hmC) in the blood was found in patients with ISM as compared to the controls (0.022% vs. 0.042%, p = 0.0001). Differences in the markers of global DNA methylation (5-mC, Alu, LINE-1) were not statistically significant, although they did indicate generally higher DNA methylation in patients with mastocytosis. The 5-hmC level was significantly associated with allergy (p = 0.011) in patients with ISM, showing a higher level of 5-hmC in patients with allergy as compared to patients without allergy. The in vitro study revealed significant differences between the studied cell lines at the level of 5-mC, Alu, and LINE-1.

CONCLUSIONS

This study confirms that epigenetic changes are involved in mastocytosis, and suggests that allergy may be an important epigenetic modifier of the disease. A possible association between KIT mutations and methylation status observed in human mast cell lines requires further investigation in human studies.

CLINICAL IMPLICATIONS

Epigenetic alterations are involved in mastocytosis pathology. The possible role of allergy as an important epigenetic modifier suggests the more impaired function of mast cells in ISM patients without allergy.

CAPSULE SUMMARY

Decreased DNA demethylation in the blood DNA of patients with ISM confirms that epigenetic alterations are involved in mastocytosis pathology. We observed a possible role of allergy as an important epigenetic modifier. There is a possible association between KIT mutations and the methylation status observed in human mast cell lines.

Programmed DNA elimination: silencing genes and repetitive sequences in somatic cells

In a multicellular organism, the genomes of all cells are in general the same. Programmed DNA elimination is a notable exception to this genome constancy rule. DNA elimination removes genes and repetitive elements in the germline genome to form a reduced somatic genome in various organisms. The process of DNA elimination within an organism is highly accurate and reproducible; it typically occurs during early embryogenesis, coincident with germline-soma differentiation. DNA elimination provides a mechanism to silence selected genes and repeats in somatic cells. Recent studies in nematodes suggest that DNA elimination removes all chromosome ends, resolves sex chromosome fusions, and may also promote the birth of novel genes. Programmed DNA elimination processes are diverse among species, suggesting DNA elimination likely has evolved multiple times in different taxa. The growing list of organisms that undergo DNA elimination indicates that DNA elimination may be more widespread than previously appreciated. These various organisms will serve as complementary and comparative models to study the function, mechanism, and evolution of programmed DNA elimination in metazoans.

The local integration preference of the Tf1 retrotransposon in Schizosaccharomyces pombe

Transposons are mobile DNAs that can move to different locations in host genomes. The integration site selection of transposons is critical for both themselves and host cells. Studies on the integration of retrotransposons and retroviruses have focused more on the global preference than on the local preference. The local preferences of retrotransposons are usually weak and of large diversity. Here, we analyzed hundreds of thousands of independent integration events of the Tf1 retrotransposon in Schizosaccharomyces pombe. The consensus sequence at the Tf1 integration sites shows a palindromic pattern, which can be divided into four sections, each of them contains one or more CGnTA units with a period of 10 base pairs, indicating interaction with subunits of the integrase oligomer in the pre-integration complex. Moreover, the analysis on the nucleosome occupancy flanking Tf1 target sites shows that Tf1 integration favors regions with one entire nucleosome depletion.

Identification and Characterization of New Alu Element Insertion in the BRCA1 Exon 14 Associated with Hereditary Breast and Ovarian Cancer

Hereditary breast and ovarian cancer syndrome (HBOC) is an autosomal dominant cancer predisposition syndrome characterized by an increased risk of breast and ovarian cancers. Germline pathogenic variants in BRCA1 are found in about 7–10% of all familial breast cancers and 10% of ovarian cancers. Alu elements are the most abundant mobile DNA element in the human genome and are known to affect the human genome by different mechanisms leading to human disease. We report here the detection, by next-generation sequencing (NGS) analysis coupled with a suitable bioinformatics pipeline, of an AluYb8 element in exon 14 of the BRCA1 gene in a family with HBOC history first classified as BRCA-negative by Sanger sequencing and first NGS analysis. The c.4475_c.4476insAluYb8 mutation impacts splicing and induces the skipping of exon 14. As a result, the produced mRNA contains a premature stop, leading to the production of a short and likely non-functional protein (pAla1453Glyfs*10). Overall, our study allowed us to identify a novel pathogenic variant in BRCA1 and showed the importance of bioinformatics tool improvement and versioning.

The Evolutionary Volte-Face of Transposable Elements: From Harmful Jumping Genes to Major Drivers of Genetic Innovation

Transposable elements (TEs) are self-replicating DNA elements that constitute major fractions of eukaryote genomes. Their ability to transpose can modify the genome structure with potentially deleterious effects. To repress TE activity, host cells have developed numerous strategies, including epigenetic pathways, such as DNA methylation or histone modifications. Although TE neo-insertions are mostly deleterious or neutral, they can become advantageous for the host under specific circumstances. The phenomenon leading to the appropriation of TE-derived sequences by the host is known as TE exaptation or co-option. TE exaptation can be of different natures, through the production of coding or non-coding DNA sequences with ultimately an adaptive benefit for the host. In this review, we first give new insights into the silencing pathways controlling TE activity. We then discuss a model to explain how, under specific environmental conditions, TEs are unleashed, leading to a TE burst and neo-insertions, with potential benefits for the host. Finally, we review our current knowledge of coding and non-coding TE exaptation by providing several examples in various organisms and describing a method to identify TE co-option events.

SINEUPs: a novel toolbox for RNA therapeutics

RNA molecules have emerged as a new class of promising therapeutics to expand the range of druggable targets in the genome. In addition to ‘canonical’ protein-coding mRNAs, the emerging richness of sense and antisense long non-coding RNAs (lncRNAs) provides a new reservoir of molecular tools for RNA-based drugs. LncRNAs are composed of modular structural domains with specific activities involving the recruitment of protein cofactors or directly interacting with nucleic acids. A single therapeutic RNA transcript can then be assembled combining domains with defined secondary structures and functions, and antisense sequences specific for the RNA/DNA target of interest.

As the first representative molecules of this new pharmacology, we have identified SINEUPs, a new functional class of natural antisense lncRNAs that increase the translation of partially overlapping mRNAs. Their activity is based on the combination of two domains: an embedded mouse inverted SINEB2 element that enhances mRNA translation (effector domain) and an overlapping antisense region that provides specificity for the target sense transcript (binding domain). By genetic engineering, synthetic SINEUPs can potentially target any mRNA of interest increasing translation and therefore the endogenous level of the encoded protein.

In this review, we describe the state-of-the-art knowledge of SINEUPs and discuss recent publications showing their potential application in diseases where a physiological increase of endogenous protein expression can be therapeutic.

Hallmarks of Retroelement Expression in T-Cells Treated With HDAC Inhibitors

A wide spectrum of drugs have been assessed as latency reversal agents (LRA) to reactivate HIV-1 from cellular reservoirs and aid in viral eradication strategies. Histone deacetylase inhibitors (HDACi) have been studied in vitro and in vivo as potential candidates for HIV-1 latency reversion. Suberoylanilide hydroxamic acid (SAHA) and romidepsin (RMD) are two HDACi able to reverse HIV latency, however studies of potential off-target effects on retroelement expression have been limited. Retroelements constitute a large portion of the human genome, and some are considered “fossil viruses” as they constitute remnants of ancient exogenous retroviruses infections. Retroelements are reactivated during certain disease conditions like cancer or during HIV-1 infection. In this study, we analyzed differential expression of retroelements using publicly available RNA-seq datasets (GSE102187 and GSE114883) obtained from uninfected CD4+, and HIV-1 latently infected CD4+ T-cells treated with HDACi (SAHA and RMD). We found a total of 712 and 1,380 differentially expressed retroelements in HIV-1 latently infected cells following a 24-h SAHA and RMD treatment, respectively. Furthermore, we found that 531 retroelement sequences (HERVs and L1) were differentially expressed under both HDACi treatments, while 1,030 HERV/L1 were exclusively regulated by each drug. Despite differences in specific HERV loci expression, the overall pattern at the HERV family level was similar for both treatments. We detected differential expression of full-length HERV families including HERV-K, HERV-W and HERV-H. Furthermore, we analyzed the link between differentially expressed retroelements and nearby immune genes. TRAF2 (TNF receptor) and GBP5 (inflammasome activator) were upregulated in HDACi treated samples and their expression was correlated with nearby HERV (MERV101_9q34.3) and L1 (L1FLnI_1p22.2k, L1FLnI_1p22.2j, L1FLnI_1p22.2i). Our findings suggest that HDACi have an off-target effect on the expression of retroelements and on the expression of immune associated genes in treated CD4+ T-cells. Furthermore, our data highlights the importance of exploring the interaction between HIV-1 and retroelement expression in LRA treated samples to understand their role and impact on “shock and kill” strategies and their potential use as reservoir biomarkers.

Genomic approaches to trace the history of human brain evolution with an emerging opportunity for transposon profiling of ancient humans

Transposable elements (TEs) significantly contribute to shaping the diversity of the human genome, and lines of evidence suggest TEs as one of driving forces of human brain evolution. Existing computational approaches, including cross-species comparative genomics and population genetic modeling, can be adapted for the study of the role of TEs in evolution. In particular, diverse ancient and archaic human genome sequences are increasingly available, allowing reconstruction of past human migration events and holding the promise of identifying and tracking TEs among other evolutionarily important genetic variants at an unprecedented spatiotemporal resolution. However, highly degraded short DNA templates and other unique challenges presented by ancient human DNA call for major changes in current experimental and computational procedures to enable the identification of evolutionarily important TEs. Ancient human genomes are valuable resources for investigating TEs in the evolutionary context, and efforts to explore ancient human genomes will potentially provide a novel perspective on the genetic mechanism of human brain evolution and inspire a variety of technological and methodological advances. In this review, we summarize computational and experimental approaches that can be adapted to identify and validate evolutionarily important TEs, especially for human brain evolution. We also highlight strategies that leverage ancient genomic data and discuss unique challenges in ancient transposon genomics.

Genomic impact of stress-induced transposable element mobility in Arabidopsis

Transposable elements (TEs) have long been known to be major contributors to plant evolution, adaptation and crop domestication. Stress-induced TE mobilization is of particular interest because it may result in novel gene regulatory pathways responding to stresses and thereby contribute to stress adaptation. Here, we investigated the genomic impacts of stress induced TE mobilization in wild type Arabidopsis plants. We find that the heat-stress responsive ONSEN TE displays an insertion site preference that is associated with specific chromatin states, especially those rich in H2A.Z histone variant and H3K27me3 histone mark. In order to better understand how novel ONSEN insertions affect the plant's response to heat stress, we carried out an in-depth transcriptomic analysis. We find that in addition to simple gene knockouts, ONSEN can produce a plethora of gene expression changes such as: constitutive activation of gene expression, alternative splicing, acquisition of heat-responsiveness, exonisation and genesis of novel non-coding and antisense RNAs. This report shows how the mobilization of a single TE-family can lead to a rapid rise of its copy number increasing the host's genome size and contribute to a broad range of transcriptomic novelty on which natural selection can then act.

Alu repetitive sequence CpG methylation changes in burn scars

Alu elements are retrotransposons related to epigenetic modifications. To date, the role of epigenetics in hypertrophic scars from burn remains unknown. Here, our aim was to examine the pathophysiology of hypertrophic scars from an epigenetic perspective. For that, we performed a cross-sectional analytical study using tissue and blood samples from burned and healthy patients (n = 23 each) to detect Alu methylation levels and patterns. The results of the combined bisulfite restriction analysis technique were categorized into four groups based on the methylation status at the CpG dinucleotides from the 5' to the 3' ends of the Alu sequence: hypermethylated (^mC^mC), hypomethylated (^uC^uC), and partially methylated (^uC^mC and ^mC^uC). Alu methylation levels were significantly lower in hypertrophic scar tissues than in normal skin (29.37 ± 2.49% vs. 35.56 ± 3.18%, p = 0.0002). In contrast, the levels were significantly higher in white blood cells from blood samples of burned patients than in those of control blood samples (26.92 ± 4.04% vs. 24.58 ± 3.34%, p = 0.0278). Alu total methylation (^mC) and the ^uC^mC pattern were significantly lower, whereas ^uC^uC was significantly higher, in hypertrophic scar tissues than in normal skin (p < 0.0001). Receiver operating characteristic analysis indicated that the ^uC^mC and ^uC^uC patterns are useful as hypertrophic scar DNA methylation markers after burn, with 91.30% sensitivity and 96.23% specificity and 100% sensitivity and 94.23% specificity, respectively. Our findings suggest that epigenetic modifications play a major role in hypertrophic scar pathogenesis, and may be the starting point for developing a novel technique for burn scar treatment.

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.

The Long Terminal Repeats of ERV6 Are Activated in Pre-Implantation Embryos of Cynomolgus Monkey

Precise gene regulation is critical during embryo development. Long terminal repeat elements (LTRs) of endogenous retroviruses (ERVs) are dynamically expressed in blastocysts of mammalian embryos. However, the expression pattern of LTRs in monkey blastocyst is still unknown. By single-cell RNA-sequencing (seq) data of cynomolgus monkeys, we found that LTRs of several ERV families, including MacERV6, MacERV3, MacERV2, MacERVK1, and MacERVK2, were highly expressed in pre-implantation embryo cells including epiblast (EPI), trophectoderm (TrB), and primitive endoderm (PrE), but were depleted in post-implantation. We knocked down MacERV6-LTR1a in cynomolgus monkeys with a short hairpin RNA (shRNA) strategy to examine the potential function of MacERV6-LTR1a in the early development of monkey embryos. The silence of MacERV6-LTR1a mainly postpones the differentiation of TrB, EPI, and PrE cells in embryos at day 7 compared to control. Moreover, we confirmed MacERV6-LTR1a could recruit Estrogen Related Receptor Beta (ESRRB), which plays an important role in the maintenance of self-renewal and pluripotency of embryonic and trophoblast stem cells through different signaling pathways including FGF and Wnt signaling pathways. In summary, these results suggest that MacERV6-LTR1a is involved in gene regulation of the pre-implantation embryo of the cynomolgus monkeys.

Factors Regulating the Activity of LINE1 Retrotransposons

LINE-1 (L1) is a class of autonomous mobile genetic elements that form somatic mosaicisms in various tissues of the organism. The activity of L1 retrotransposons is strictly controlled by many factors in somatic and germ cells at all stages of ontogenesis. Alteration of L1 activity was noted in a number of diseases: in neuropsychiatric and autoimmune diseases, as well as in various forms of cancer. Altered activity of L1 retrotransposons for some pathologies is associated with epigenetic changes and defects in the genes involved in their repression. This review discusses the molecular genetic mechanisms of the retrotransposition and regulation of the activity of L1 elements. The contribution of various factors controlling the expression and distribution of L1 elements in the genome occurs at all stages of the retrotransposition. The regulation of L1 elements at the transcriptional, post-transcriptional and integration into the genome stages is described in detail. Finally, this review also focuses on the evolutionary aspects of L1 accumulation and their interplay with the host regulation system.

Processed pseudogenes: A substrate for evolutionary innovation: Retrotransposition contributes to genome evolution by propagating pseudogene sequences with rich regulatory potential throughout the genome

Processed pseudogenes may serve as a genetic reservoir for evolutionary innovation. Here, we argue that through the activity of long interspersed element-1 retrotransposons, processed pseudogenes disperse coding and noncoding sequences rich with regulatory potential throughout the human genome. While these sequences may appear to be non-functional, a lack of contemporary function does not prohibit future development of biological activity. Here, we discuss the dynamic evolution of certain processed pseudogenes into coding and noncoding genes and regulatory elements, and their implication in wide-ranging biological and pathological processes. Also see the video abstract here: https://youtu.be/iUY_mteVoPI.

Tinzaparin inhibits VL30 retrotransposition induced by oxidative stress and/or VEGF in HC11 mouse progenitor mammary cells: Association between inhibition of cancer stem cell proliferation and mammosphere disaggregation

Tinzaparin is an anticoagulant and antiangiogenic drug with inhibitory properties against tumor growth. VEGF stimulates angiogenesis, while an association between reactive oxygen species (ROS) and angiogenesis is involved in tumor progression. The present study aimed to investigate the effect of tinzaparin on VL30 retrotransposition-positive mouse HC11 mammary stem-like epithelial cells, previously reported to be associated with induced mammosphere/cancer stem cell (CSC) generation and tumorigenesis. Under 24 h serum starvation, 15.2% nominal retrotransposition frequency was increased to 29%. Additionally, while treatment with 3–12 ng/ml VEGF further induced retrotransposition frequency in a dose-dependent manner (up to 40.3%), pre-incubation with tinzaparin (2 IU/ml) for 0.5–4 h reduced this frequency to 18.3% in a time-dependent manner, confirmed by analogous results in NIH3T3 fibroblasts. Treatment with 10–40 pg/ml glucose oxidase (GO) for 24 h induced HC11 cell retrotransposition in a dose-dependent manner (up to 82.5%), while a 3 h pre-incubation with tinzaparin (1 or 2 IU/ml) elicited a 13.5 or 25.5% reduction in retrotransposition, respectively. Regarding tumorigenic VL30 retrotransposition-positive HC11 cells, treatment with 2 IU/ml tinzaparin for 5 days reduced proliferation rate in a time-dependent manner (up to ~55%), and after 3 weeks, disaggregated soft agar-formed foci, as well as low-adherent mammospheres, producing single mesenchymal-like cells with a ~50% reduced retrotransposition. With respect to the VL30 retrotransposition mechanism: While 12 ng/ml VEGF increased the level of VL30 and endogenous reverse transcriptase (enRT) transcripts ~1.41- and ~1.16-fold, respectively, subsequent tinzaparin treatment reduced both endogenous/ROS- and VEGF-induced levels 1.15- and 0.40-fold (VL30) and 0.60- and 0.52-fold (enRT), respectively. To the best of our knowledge, these data demonstrate for the first time, the novel inhibition activity of tinzaparin against ROS- and VEGF-induced VL30 retrotransposition, and the proliferation and/or aggregation of mouse HC11 mammosphere/tumor-initiating CSCs, thus contributing to the inhibition of VL30 retrotransposition-induced primary tumor growth.

Molecular Features of Polycystic Ovary Syndrome Revealed by Transcriptome Analysis of Oocytes and Cumulus Cells

Polycystic ovary syndrome (PCOS) is typically characterized by a polycystic ovarian morphology, hyperandrogenism, ovulatory dysfunction, and infertility. Furthermore, PCOS patients undergoing ovarian stimulation have more oocytes; however, the poor quality of oocytes leads to lower fertilization and implantation rates, decreased pregnancy rates, and increased miscarriage rates. The complex molecular mechanisms underlying PCOS and the poor quality of oocytes remain to be elucidated. We obtained matched oocytes and cumulus cells (CCs) from PCOS patients, compared them with age-matched controls, and performed RNA sequencing analysis to explore the transcriptional characteristics of their oocytes and CCs. Moreover, we validated our newly confirmed candidate genes for PCOS by immunofluorescence. Unsupervised clustering analysis showed that the overall global gene expression patterns and transposable element (TE) expression profiles of PCOS patients tightly clustered together, clearly distinct from those of controls. Abnormalities in functionally important pathways are found in PCOS oocytes. Notably, genes involved in microtubule processes, TUBB8 and TUBA1C, are overexpressed in PCOS oocytes. The metabolic and oxidative phosphorylation pathways are also dysregulated in both oocytes and CCs from PCOS patients. Moreover, in oocytes, differentially expressed TEs are not uniformly dispersed in human chromosomes. Endogenous retrovirus 1 (ERV1) elements located on chromosomes 2, 3, 4, and 5 are rather highly upregulated. Interestingly, these correlate with the most highly expressed protein-coding genes, including tubulin-associated genes TUBA1C, TUBB8P8, and TUBB8, linking the ERV1 elements to the occurrence of PCOS. Our comprehensive analysis of gene expression in oocytes and CCs, including TE expression, revealed the specific molecular features of PCOS. The aberrantly elevated expression of TUBB8 and TUBA1C and ERV1 provides additional markers for PCOS and may contribute to the compromised oocyte developmental competence in PCOS patients. Our findings may also have implications for treatment strategies to improve oocyte maturation and the pregnancy outcomes for women with PCOS.

Open Reading Frame 1 Protein of the Human Long Interspersed Nuclear Element 1 Retrotransposon Binds Multiple Equivalents of Lead

The human long interspersed nuclear element 1 (LINE1) has been implicated in numerous diseases and has been suggested to play a significant role in genetic evolution. Open reading frame 1 protein (ORF1p) is one of the two proteins encoded in this self-replicating mobile genetic element, both of which are essential for retrotransposition. The structure of the three-stranded coiled-coil domain of ORF1p was recently solved and showed the presence of tris-cysteine layers in the interior of the coiled-coil that could function as metal binding sites. Here, we demonstrate that ORF1p binds Pb(II). We designed a model peptide, GRCSL16CL23C, to mimic two of the ORF1p Cys3 layers and crystallized the peptide both as the apo-form and in the presence of Pb(II). Structural comparison of the ORF1p with apo-(GRCSL16CL23C)3 shows very similar Cys3 layers, preorganized for Pb(II) binding. We propose that exposure to heavy metals, such as lead, could influence directly the structural parameters of ORF1p and thus impact the overall LINE1 retrotransposition frequency, directly relating heavy metal exposure to genetic modification.

A Model-Driven Quantitative Analysis of Retrotransposon Distributions in the Human Genome

Retrotransposons, DNA sequences capable of creating copies of themselves, compose about half of the human genome and played a central role in the evolution of mammals. Their current position in the host genome is the result of the retrotranscription process and of the following host genome evolution. We apply a model from statistical physics to show that the genomic distribution of the two most populated classes of retrotransposons in human deviates from random placement, and that this deviation increases with time. The time dependence suggests a major role of the host genome dynamics in shaping the current retrotransposon distributions. Focusing on a neutral scenario, we show that a simple model based on random placement followed by genome expansion and sequence duplications can reproduce the empirical retrotransposon distributions, even though more complex and possibly selective mechanisms can have contributed. Besides the inherent interest in understanding the origin of current retrotransposon distributions, this work sets a general analytical framework to analyze quantitatively the effects of genome evolutionary dynamics on the distribution of genomic elements.

PRE-1 Revealed Previous Unknown Introgression Events in Eurasian Boars during the Middle Pleistocene

Introgression events and population admixture occurred among Sus species across the Eurasian mainland in the Middle Pleistocene, which reflects the local adaption of different populations and contributes to evolutionary novelty. Previous findings on these population introgressions were largely based on extensive genome-wide single-nucleotide polymorphism information, ignoring structural variants (SVs) as an important alternative resource of genetic variations. Here, we profiled the genome-wide SVs and explored the formation of pattern-related SVs, indicating that PRE1-SS is a recently active subfamily that was strongly associated with introgression events in multiple Asian and European pig populations. As reflected by the three different combination haplotypes from two specific patterns and known phylogenetic relationships in Eurasian boars, we identified the Asian Northern wild pigs as having experienced introgression from European wild boars around 0.5–0.2 Ma and having received latitude-related selection. During further exploration of the influence of pattern-related SVs on gene functions, we found substantial sequence changes in 199 intron regions of 54 genes and 3 exon regions of 3 genes (HDX, TRO, and SMIM1), implying that the pattern-related SVs were highly related to positive selection and adaption of pigs. Our findings revealed novel introgression events in Eurasian wild boars, providing a timeline of population admixture and divergence across the Eurasian mainland in the Middle Pleistocene.

One-carbon metabolism and global DNA methylation in mothers of individuals with Down syndrome

Down syndrome (DS) is the most common chromosomal disorder, resulting from the failure of normal chromosome 21 segregation. Studies have suggested that impairments within the one-carbon metabolic pathway can be of relevance for the global genome instability observed in mothers of individuals with DS. Based on the association between global DNA hypomethylation, genome instability, and impairments within the one-carbon metabolic pathway, the present study aimed to identify possible predictors, within the one-carbon metabolism, of global DNA methylation, measured by methylation patterns of LINE-1 and Alu repetitive sequences, in mothers of individuals with DS and mothers of individuals without the syndrome. In addition, we investigated one-carbon genetic polymorphisms and metabolites as maternal predisposing factors for the occurrence of trisomy 21 in children. Eighty-three samples of mothers of children with DS with karyotypically confirmed free trisomy 21 (case group) and 84 of mothers who had at least one child without DS or any other aneuploidy were included in the study. Pyrosequencing assays were performed to access global methylation. The results showed that group affiliation (case or control), betaine-homocysteine methyltransferase (BHMT) G742A and transcobalamin 2 (TCN2) C776G polymorphisms, and folate concentration were identified as predictors of global Alu DNA methylation values. In addition, thymidylate synthase (TYMS) 28-bp repeats 2R/3R or 3R/3R genotypes are independent maternal predisposing factors for having a child with DS. This study adds evidence that supports the association of impairments in the one-carbon metabolism, global DNA methylation, and the possibility of having a child with DS.