Retroelements and their impact on genome evolution and functioning

Chromosome-contiguous genome for the Haecon-5 strain of Haemonchus contortus reveals marked genetic variability and enables the discovery of essential gene candidates

Millions of livestock animals worldwide are infected with the haematophagous barber's pole worm, Haemonchus contortus, the aetiological agent of haemonchosis. Despite the major significance of this parasite worldwide and its widespread resistance to current treatments, the lack of a high-quality genome for the well-defined strain of this parasite from Australia, called Haecon-5, has constrained research in a number of areas including host-parasite interactions, drug discovery and population genetics. To enable research in these areas, we report here a chromosome-contiguous genome (∼280 Mb) for Haecon-5 with high-quality models for 19,234 protein-coding genes. Comparative genomic analyses show significant genomic similarity (synteny) with a UK strain of H. contortus, called MHco3(ISE).N1 (abbreviated as "ISE"), but we also discover marked differences in genomic structure/gene arrangements, distribution of nucleotide variability (single nucleotide polymorphisms (SNPs) and indels) and orthology between Haecon-5 and ISE. We used the genome and extensive transcriptomic resources for Haecon-5 to predict a subset of essential single-copy genes employing a "cross-species" machine learning (ML) approach using a range of features from nucleotide/protein sequences, protein orthology, subcellular localisation, single-cell RNA-seq and/or histone methylation data available for the model organisms Caenorhabditis elegans and Drosophila melanogaster. From a set of 1,464 conserved single copy genes, transcribed in key life-cycle stages of H. contortus, we identified 232 genes whose homologs have critical functions in C. elegans and/or D. melanogaster, and prioritised 10 of them for further characterisation; nine of the 10 genes likely play roles in neurophysiological processes, germline, hypodermis and/or respiration, and one is an unknown (orphan) gene for which no detailed functional information exists. Future studies of these genes/gene products are warranted to elucidate their roles in parasite biology, host-parasite interplay and/or disease. Clearly, the present Haecon-5 reference genome and associated resources now underpin a broad range of fundamental investigations of H. contortus and could assist in accelerating the discovery of novel intervention targets and drug candidates to combat haemonchosis.

Number of human protein interactions correlates with structural, but not regulatory conservation of the respective genes

Introduction

The differential ratio of nonsynonymous to synonymous nucleotide substitutions (dN/dS) is a common measure of the rate of structural evolution in proteincoding genes. In addition, we recently suggested that the proportion of transposable elements in gene promoters that host functional genomic sites serves as a marker of the rate of regulatory evolution of genes. Such functional genomic regions may include transcription factor binding sites and modified histone binding loci.

Methods

Here, we constructed a model of the human interactome based on 600,136 documented molecular interactions and investigated the overall relationship between the number of interactions of each protein and the rate of structural and regulatory evolution of the corresponding genes.

Results

By evaluating a total of 4,505 human genes and 1,936 molecular pathways we found a general correlation between structural and regulatory evolution rate metrics (Spearman 0.08-0.16 and 0.25-0.37 for gene and pathway levels, respectively, p < 0.01). Further exploration revealed in the established human interactome model lack of correlation between the rate of gene regulatory evolution and the number of protein interactions on gene level, and weak negative correlation (∼0.15) on pathway level. We also found a statistically significant negative correlation between the rate of gene structural evolution and the number of protein interactions (Spearman -0.11 and -0.3 for gene and pathway levels, respectively, p < 0.01).

Discussion

Our result suggests stronger structural rather than regulatory conservation of genes whose protein products have multiple interaction partners.

Cancer fusion transcripts with human non-coding RNAs

Cancer chimeric, or fusion, transcripts are thought to most frequently appear due to chromosomal aberrations that combine moieties of unrelated normal genes. When being expressed, this results in chimeric RNAs having upstream and downstream parts relatively to the breakpoint position for the 5'- and 3'-fusion components, respectively. As many other types of cancer mutations, fusion genes can be of either driver or passenger type. The driver fusions may have pivotal roles in malignisation by regulating survival, growth, and proliferation of tumor cells, whereas the passenger fusions most likely have no specific function in cancer. The majority of research on fusion gene formation events is concentrated on identifying fusion proteins through chimeric transcripts. However, contemporary studies evidence that fusion events involving non-coding RNA (ncRNA) genes may also have strong oncogenic potential. In this review we highlight most frequent classes of ncRNAs fusions and summarize current understanding of their functional roles. In many cases, cancer ncRNA fusion can result in altered concentration of the non-coding RNA itself, or it can promote protein expression from the protein-coding fusion moiety. Differential splicing, in turn, can enrich the repertoire of cancer chimeric transcripts, e.g. as observed for the fusions of circular RNAs and long non-coding RNAs. These and other ncRNA fusions are being increasingly recognized as cancer biomarkers and even potential therapeutic targets. Finally, we discuss the use of ncRNA fusion genes in the context of cancer detection and therapy.

Indicaxanthin Induces Autophagy in Intestinal Epithelial Cancer Cells by Epigenetic Mechanisms Involving DNA Methylation

Autophagy is an evolutionarily conserved process critical in maintaining cellular homeostasis. Recently, the anticancer potential of autophagy inducers, including phytochemicals, was suggested. Indicaxanthin is a betalain pigment found in prickly pear fruit with antiproliferative and pro-apoptotic activities in colorectal cancer cells associated with epigenetic changes in selected methylation-silenced oncosuppressor genes. Here, we demonstrate that indicaxanthin induces the up-regulation of the autophagic markers LC3-II and Beclin1, and increases autophagolysosome production in Caco-2 cells. Methylomic studies showed that the indicaxanthin-induced pro-autophagic activity was associated with epigenetic changes. In addition to acting as a hypermethylating agent at the genomic level, indicaxanthin also induced significant differential methylation in 39 out of 47 autophagy-related genes, particularly those involved in the late stages of autophagy. Furthermore, in silico molecular modelling studies suggested a direct interaction of indicaxanthin with Bcl-2, which, in turn, influenced the function of Beclin1, a key autophagy regulator. External effectors, including food components, may modulate the epigenetic signature of cancer cells. This study demonstrates, for the first time, the pro-autophagic potential of indicaxanthin in human colorectal cancer cells associated with epigenetic changes and contributes to outlining its potential healthy effect in the pathophysiology of the gastrointestinal tract.

Distinct Traits of Structural and Regulatory Evolutional Conservation of Human Genes with Specific Focus on Major Cancer Molecular Pathways

The evolution of protein-coding genes has both structural and regulatory components. The first can be assessed by measuring the ratio of non-synonymous to synonymous nucleotide substitutions. The second component can be measured as the normalized proportion of transposable elements that are used as regulatory elements. For the first time, we characterized in parallel the regulatory and structural evolutionary profiles for 10,890 human genes and 2972 molecular pathways. We observed a ~0.1 correlation between the structural and regulatory metrics at the gene level, which appeared much higher (~0.4) at the pathway level. We deposited the data in the publicly available database RetroSpect. We also analyzed the evolutionary dynamics of six cancer pathways of two major axes: Notch/WNT/Hedgehog and AKT/mTOR/EGFR. The Hedgehog pathway had both components slower, whereas the Akt pathway had clearly accelerated structural evolution. In particular, the major hub nodes Akt and beta-catenin showed both components strongly decreased, whereas two major regulators of Akt TCL1 and CTMP had outstandingly high evolutionary rates. We also noticed structural conservation of serine/threonine kinases and the genes related to guanosine metabolism in cancer signaling: GPCRs, G proteins, and small regulatory GTPases (Src, Rac, Ras); however, this was compensated by the accelerated regulatory evolution.

Identification and Characterization of the HERV-K (HML-8) Group of Human Endogenous Retroviruses in the Genome

Human endogenous retroviruses (HERVs) can be vertically transmitted in a Mendelian fashion, are stably maintained in the human genome, and are estimated to constitute ∼8% of the genome. HERVs affect human physiology and pathology through their provirus-encoded protein or long terminal repeat (LTR) element effect. Characterization of the genomic distribution is an essential step to understanding the relationships between endogenous retrovirus expression and diseases. However, the poor characterization of human MMTV-like (HML)-8 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, the definition of an accurate and updated HERV-K HML-8 genomic map is urgently needed. In this study, we report the results of a comprehensive analysis of HERV-K HML-8 sequence presence and distribution within the human genome and hominoids, with a detailed description of the different structural and phylogenetic aspects characterizing the group. A total of 40 proviruses and 5 solo LTR elements for human were characterized, which included a detailed description of provirus structure, integration time, potentially regulated genes, transcription factor-binding sites, and primer-binding site features. Besides, 9 chimpanzee sequences, 8 gorilla sequences, and 10 orangutan sequences belonging to the HML-8 subgroup were identified. The integration time results showed that the HML-8 elements were integrated into the primate lineage around 35 and 42 million years ago (mya), during primates evolutionary speciation. Overall, the results clarified the composition of the HML-8 groups, providing an exhaustive background for subsequent functional studies.

Endogenous Retrovirus Elements Are Co-Expressed with IFN Stimulation Genes in the JAK-STAT Pathway

Background: Endogenous retrovirus (ERV) elements can act as proximal regulatory elements in promoting interferon (IFN) responses. Previous relevant studies have mainly focused on IFN-stimulated genes (ISGs). However, the role of ERV elements as cis-regulatory motifs in regulating genes of the JAK-STAT pathway remains poorly understood. In our study, we analyzed the changes in ERV elements and genes under both IFN stimulation and blockade of the signaling pathway. Methods: The effects of interferon on cells under normal conditions and knockout of the receptor were compared based on the THP1_IFNAR1_KO and THP1_IFNAR2_mutant cell lines. The correlation between differentially expressed ERVs (DHERVs) and differentially expressed genes (DEGs) as DEHERV-G pairs was explored with construction of gene regulatory networks related to ERV and induced by proinflammatory cytokines. Results: A total of 430 DEHERV loci and 190 DEGs were identified in 842 DEHERV-G pairs that are common to the three groups. More than 87% of DEHERV-G pairs demonstrated a consistent expression pattern. ISGs such as AIM2, IFIT1, IFIT2, IFIT3, STAT1, and IRF were activated via the JAK-STAT pathway in response to interferon stimulation. Thus, STAT1, STAT2, and IRF1 appear to play core roles in regulatory networks and are closely associated with ERVs. Conclusions: The RNA expression of ISGs and ERV elements is correlated, indicating that ERV elements are closely linked to host innate immune responses.

From Nucleus to Organs: Insights of Aryl Hydrocarbon Receptor Molecular Mechanisms

The aryl hydrocarbon receptor (AHR) is a markedly established regulator of a plethora of cellular and molecular processes. Its initial role in the detoxification of xenobiotic compounds has been partially overshadowed by its involvement in homeostatic and organ physiology processes. In fact, the discovery of its ability to bind specific target regulatory sequences has allowed for the understanding of how AHR modulates such processes. Thereby, AHR presents functions in transcriptional regulation, chromatin architecture modifications and participation in different key signaling pathways. Interestingly, such fields of influence end up affecting organ and tissue homeostasis, including regenerative response both to endogenous and exogenous stimuli. Therefore, from classical spheres such as canonical transcriptional regulation in embryonic development, cell migration, differentiation or tumor progression to modern approaches in epigenetics, senescence, immune system or microbiome, this review covers all aspects derived from the balance between regulation/deregulation of AHR and its physio-pathological consequences.

Population genomics of the neotropical palm Copernicia prunifera (Miller) H. E. Moore: Implications for conservation

Copernicia prunifera (Miller) H. E. Moore is a palm tree native to Brazil. The products obtained from its leaf extracts are a source of income for local families and the agroindustry. Owing to the reduction of natural habitats and the absence of a sustainable management plan, the maintenance of the natural populations of this palm tree has been compromised. Therefore, this study aimed to evaluate the diversity and genetic structure of 14 C. prunifera populations using single nucleotide polymorphisms (SNPs) identified through genotyping-by-sequencing (GBS) to provide information that contributes to the conservation of this species. A total of 1,013 SNP markers were identified, of which 84 loci showed outlier behavior and may reflect responses to natural selection. Overall, the level of genomic diversity was compatible with the biological aspects of this species. The inbreeding coefficient (f) was negative for all populations, indicating excess heterozygotes. Most genetic variations occurred within populations (77.26%), and a positive correlation existed between genetic and geographic distances. The population structure evaluated through discriminant analysis of principal components (DAPC) revealed low genetic differentiation between populations. The results highlight the need for efforts to conserve C. prunifera as well as its distribution range to preserve its global genetic diversity and evolutionary potential.

Intragenic L1 Insertion: One Possibility of Brain Disorder

Long interspersed nuclear element 1 (LINE1, L1) is a retrotransposon comprising ~17% of the human genome. A subset of L1s maintains the potential to mobilize and alter the genomic landscape, consequently contributing to the change in genome integrity and gene expression. L1 retrotransposition occurs in the human brain regardless of disease status. However, in the brain of patients with various brain diseases, the expression level and copy number of L1 are significantly increased. In this review, we briefly introduce the methodologies applied to measure L1 mobility and identify genomic loci where new insertion of L1 occurs in the brain. Then, we present a list of genes disrupted by L1 transposition in the genome of patients with brain disorders. Finally, we discuss the association between genes disrupted by L1 and relative brain disorders.

Comprehensive identification and characterization of the HERV-K (HML-9) group in the human genome

BACKGROUND

Human endogenous retroviruses (HERVs) result from ancestral infections caused by exogenous retroviruses that became incorporated into the germline DNA and evolutionarily fixed in the human genome. HERVs can be transmitted vertically in a Mendelian fashion and be stably maintained in the human genome, of which they are estimated to comprise approximately 8%. HERV-K (HML1-10) transcription has been confirmed to be associated with a variety of diseases, such as breast cancer, lung cancer, prostate cancer, melanoma, rheumatoid arthritis, and amyotrophic lateral sclerosis. However, the poor characterization of HML-9 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, a precise and updated HERV-K HML-9 genomic map is urgently needed to better evaluate the role of these elements in human health.

RESULTS

We report a comprehensive analysis of the presence and distribution of HERV-K HML-9 elements within the human genome, with a detailed characterization of the structural and phylogenetic properties of the group. A total of 23 proviruses and 47 solo LTR elements were characterized, with a detailed description of the provirus structure, integration time, potential regulated genes, transcription factor binding sites (TFBS), and primer binding site (PBS) features. The integration time results showed that the HML-9 elements found in the human genome integrated into the primate lineage between 17.5 and 48.5 million years ago (mya).

CONCLUSION

The results provide a clear characterization of HML-9 and a comprehensive background for subsequent functional studies.

Aryl Hydrocarbon Receptor: From Homeostasis to Tumor Progression

Transcription factor aryl hydrocarbon receptor (AHR) has emerged as one of the main regulators involved both in different homeostatic cell functions and tumor progression. Being a member of the family of basic-helix-loop-helix (bHLH) transcriptional regulators, this intracellular receptor has become a key member in differentiation, pluripotency, chromatin dynamics and cell reprogramming processes, with plenty of new targets identified in the last decade. Besides this role in tissue homeostasis, one enthralling feature of AHR is its capacity of acting as an oncogene or tumor suppressor depending on the specific organ, tissue and cell type. Together with its well-known modulation of cell adhesion and migration in a cell-type specific manner in epithelial-mesenchymal transition (EMT), this duality has also contributed to the arise of its clinical interest, highlighting a new potential as therapeutic tool, diagnosis and prognosis marker. Therefore, a deregulation of AHR-controlled pathways may have a causal role in contributing to physiological and homeostatic failures, tumor progression and dissemination. With that firmly in mind, this review will address the remarkable capability of AHR to exert a different function influenced by the phenotype of the target cell and its potential consequences.

Viruses Defined by the Position of the Virosphere within the Replicator Space

Originally, viruses were defined as miniscule infectious agents that passed through filters that retain even the smallest cells. Subsequently, viruses were considered obligate intracellular parasites whose reproduction depends on their cellular hosts for energy supply and molecular building blocks. However, these features are insufficient to unambiguously define viruses as they are broadly understood today. We outline possible approaches to define viruses and explore the boundaries of the virosphere within the virtual space of replicators and the relationships between viruses and other types of replicators. Regardless of how, exactly, viruses are defined, viruses clearly have evolved on many occasions from nonviral replicators, such as plasmids, by recruiting host proteins to become virion components. Conversely, other types of replicators have repeatedly evolved from viruses. Thus, the virosphere is a dynamic entity with extensive evolutionary traffic across its boundaries. We argue that the virosphere proper, here termed orthovirosphere, consists of a distinct variety of replicators that encode structural proteins encasing the replicators' genomes, thereby providing protection and facilitating transmission among hosts. Numerous and diverse replicators, such as virus-derived but capsidless RNA and DNA elements, or defective viruses occupy the zone surrounding the orthovirosphere in the virtual replicator space. We define this zone as the perivirosphere. Although intense debates on the nature of certain replicators that adorn the internal and external boundaries of the virosphere will likely continue, we present an operational definition of virus that recently has been accepted by the International Committee on Taxonomy of Viruses.

Transposons-Based Clonal Diversity in Trematode Involves Parts of CR1 (LINE) in Eu- and Heterochromatin

Trematode parthenitae have long been believed to form clonal populations, but clonal diversity has been discovered in this asexual stage of the lifecycle. Clonal polymorphism in the model species Himasthla elongata has been previously described, but the source of this phenomenon remains unknown. In this work, we traced cercarial clonal diversity using a simplified amplified fragment length polymorphism (SAFLP) method and characterised the nature of fragments in diverse electrophoretic bands. The repetitive elements were identified in both the primary sequence of the H. elongata genome and in the transcriptome data. Long-interspersed nuclear elements (LINEs) and long terminal repeat retrotransposons (LTRs) were found to represent an overwhelming majority of the genome and the transposon transcripts. Most sequenced fragments from SAFLP pattern contained the reverse transcriptase (RT, ORF2) domains of LINEs, and only a few sequences belonged to ORFs of LTRs and ORF1 of LINEs. A fragment corresponding to a CR1-like (LINE) spacer region was discovered and named CR1-renegade (CR1-rng). In addition to RT-containing CR1 transcripts, we found short CR1-rng transcripts in the redia transcriptome and short contigs in the mobilome. Probes against CR1-RT and CR1-rng presented strikingly different pictures in FISH mapping, despite both being fragments of CR1. In silico data and Southern blotting indicated that CR1-rng is not tandemly organised. CR1 involvement in clonal diversity is discussed.

Aux/IAA14 Regulates microRNA-Mediated Cold Stress Response in Arabidopsis Roots

The phytohormone auxin and microRNA-mediated regulation of gene expressions are key regulators of plant growth and development at both optimal and under low-temperature stress conditions. However, the mechanistic link between microRNA and auxin in regulating plant cold stress response remains elusive. To better understand the role of microRNA (miR) in the crosstalk between auxin and cold stress responses, we took advantage of the mutants of Arabidopsis thaliana with altered response to auxin transport and signal. Screening of the mutants for root growth recovery after cold stress at 4 °C revealed that the auxin signaling mutant, solitary root 1 (slr1; mutation in Aux/IAA14), shows a hypersensitive response to cold stress. Genome-wide expression analysis of miRs in the wild-type and slr1 mutant roots using next-generation sequencing revealed 180 known and 71 novel cold-responsive microRNAs. Cold stress also increased the abundance of 26-31 nt small RNA population in slr1 compared with wild type. Comparative analysis of microRNA expression shows significant differential expression of 13 known and 7 novel miRs in slr1 at 4 °C compared with wild type. Target gene expression analysis of the members from one potential candidate miR, miR169, revealed the possible involvement of miR169/NF-YA module in the Aux/IAA14-mediated cold stress response. Taken together, these results indicate that SLR/IAA14, a transcriptional repressor of auxin signaling, plays a crucial role in integrating miRs in auxin and cold responses.

Endogenous Retrovirus Transcript Levels Are Associated with Immunogenic Signatures in Multiple Metastatic Cancer Types

Next-generation sequencing of solid tumors has revealed variable signatures of immunogenicity across tumors, but underlying molecular characteristics driving such variation are not fully understood. Although expression of endogenous retrovirus (ERV)-containing transcripts can provide a source of tumor-specific neoantigen in some cancer models, associations between ERV levels and immunogenicity across different types of metastatic cancer are not well established. We performed bioinformatics analysis of genomic, transcriptomic, and clinical data across an integrated cohort of 199 patients with metastatic breast, colorectal, and pancreatic ductal adenocarcinoma tumors. Within each cancer type, we identified a subgroup of viral mimicry tumors in which increased ERV levels were coupled with transcriptional signatures of autonomous antiviral response and immunogenicity. In addition, viral mimicry colorectal and pancreatic tumors showed increased expression of DNA demethylation gene TET2 Taken together, these data demonstrate the existence of an ERV-associated viral mimicry phenotype across three distinct metastatic cancer types, while indicating links between ERV abundance, epigenetic dysregulation, and immunogenicity.

New Insights on the Effects of Methylphenidate in Attention Deficit Hyperactivity Disorder

This narrative review describes an overview of the multiple effects of methylphenidate (MPH) in attention-deficit/hyperactivity disorder (ADHD) and its potential neurobiological targets. It addressed the following aspects: 1) MPH effects on attention and executive functions in ADHD; 2) the relation between MPH efficacy and dopamine transporter gene (DAT) polymorphism; and 3) the role of MPH as an epigenetic modulator in ADHD. Literature analysis showed that MPH, the most commonly used psychostimulant in the therapy of ADHD, acts on multiple components of the disorder. Marked improvements in attentional and executive dysfunction have been observed in children with ADHD during treatment with MPH, as well as reductions in neurological soft signs. MPH efficacy may be influenced by polymorphisms in the DAT, and better responses to treatment were associated with the 10/10 genotype. Innovative lines of research have suggested that ADHD etiopathogenesis and its neuropsychological phenotypes also depend on the expression levels of human endogenous retrovirus (HERV). In particular, several studies have revealed that ADHD is associated with HERV-H over-expression and that MPH administration results in decreased expression levels of this retroviral family and a reduction in the main symptoms of the disorder. In conclusion, there is a confirmed role for MPH as an elective drug in the therapy of ADHD alone or in association with behavioral therapy. Its effectiveness can vary based on DAT polymorphisms and can act as a modulator of HERV-H gene expression, pointing to targets for a precision medicine approach.

Retroelement-Linked H3K4me1 Histone Tags Uncover Regulatory Evolution Trends of Gene Enhancers and Feature Quickly Evolving Molecular Processes in Human Physiology

Background: Retroelements (REs) are mobile genetic elements comprising ~40% of human DNA. They can reshape expression patterns of nearby genes by providing various regulatory sequences. The proportion of regulatory sequences held by REs can serve a measure of regulatory evolution rate of the respective genes and molecular pathways. Methods: We calculated RE-linked enrichment scores for individual genes and molecular pathways based on ENCODE project epigenome data for enhancer-specific histone modification H3K4me1 in five human cell lines. We identified consensus groups of molecular processes that are enriched and deficient in RE-linked H3K4me1 regulation. Results: We calculated H3K4me1 RE-linked enrichment scores for 24,070 human genes and 3095 molecular pathways. We ranked genes and pathways and identified those statistically significantly enriched and deficient in H3K4me1 RE-linked regulation. Conclusion: Non-coding RNA genes were statistically significantly enriched by RE-linked H3K4me1 regulatory modules, thus suggesting their high regulatory evolution rate. The processes of gene silencing by small RNAs, DNA metabolism/chromatin structure, sensory perception/neurotransmission and lipids metabolism showed signs of the fastest regulatory evolution, while the slowest processes were connected with immunity, protein ubiquitination/degradation, cell adhesion, migration and interaction, metals metabolism/ion transport, cell death, intracellular signaling pathways.

H3K4me3, H3K9ac, H3K27ac, H3K27me3 and H3K9me3 Histone Tags Suggest Distinct Regulatory Evolution of Open and Condensed Chromatin Landmarks

BACKGROUND

Transposons are selfish genetic elements that self-reproduce in host DNA. They were active during evolutionary history and now occupy almost half of mammalian genomes. Close insertions of transposons reshaped structure and regulation of many genes considerably. Co-evolution of transposons and host DNA frequently results in the formation of new regulatory regions. Previously we published a concept that the proportion of functional features held by transposons positively correlates with the rate of regulatory evolution of the respective genes.

METHODS

We ranked human genes and molecular pathways according to their regulatory evolution rates based on high throughput genome-wide data on five histone modifications (H3K4me3, H3K9ac, H3K27ac, H3K27me3, H3K9me3) linked with transposons for five human cell lines.

RESULTS

Based on the total of approximately 1.5 million histone tags, we ranked regulatory evolution rates for 25075 human genes and 3121 molecular pathways and identified groups of molecular processes that showed signs of either fast or slow regulatory evolution. However, histone tags showed different regulatory patterns and formed two distinct clusters: promoter/active chromatin tags (H3K4me3, H3K9ac, H3K27ac) vs. heterochromatin tags (H3K27me3, H3K9me3).

CONCLUSION

In humans, transposon-linked histone marks evolved in a coordinated way depending on their functional roles.

Evolution of Protein Domain Architectures

This chapter reviews current research on how protein domain architectures evolve. We begin by summarizing work on the phylogenetic distribution of proteins, as this will directly impact which domain architectures can be formed in different species. Studies relating domain family size to occurrence have shown that they generally follow power law distributions, both within genomes and larger evolutionary groups. These findings were subsequently extended to multi-domain architectures. Genome evolution models that have been suggested to explain the shape of these distributions are reviewed, as well as evidence for selective pressure to expand certain domain families more than others. Each domain has an intrinsic combinatorial propensity, and the effects of this have been studied using measures of domain versatility or promiscuity. Next, we study the principles of protein domain architecture evolution and how these have been inferred from distributions of extant domain arrangements. Following this, we review inferences of ancestral domain architecture and the conclusions concerning domain architecture evolution mechanisms that can be drawn from these. Finally, we examine whether all known cases of a given domain architecture can be assumed to have a single common origin (monophyly) or have evolved convergently (polyphyly). We end by a discussion of some available tools for computational analysis or exploitation of protein domain architectures and their evolution.

Stable coevolutionary regimes for genetic parasites and their hosts: you must differ to coevolve

Background

Genetic parasites are ubiquitous satellites of cellular life forms most of which host a variety of mobile genetic elements including transposons, plasmids and viruses. Theoretical considerations and computer simulations suggest that emergence of genetic parasites is intrinsic to evolving replicator systems.

Results

Using methods of bifurcation analysis, we investigated the stability of simple models of replicator-parasite coevolution in a well-mixed environment. We first analyze what appears to be the simplest imaginable system of this type, one in which the parasite evolves during the replication of the host genome through a minimal mutation that renders the genome of the emerging parasite incapable of producing the replicase but able to recognize and recruit it for its own replication. This model has only trivial or “semi-trivial”, parasite-free equilibria: an inefficient parasite is outcompeted by the host and dies off, whereas an efficient one pushes the host out of existence, leading to the collapse of the entire system. We show that stable host-parasite coevolution (a non-trivial equilibrium) is possible in a modified model where the parasite is qualitatively distinct from the host replicator in that the replication of the parasite depends solely on the availability of the host but not on the carrying capacity of the environment.

Conclusions

We analytically determine the conditions for stable coevolution of genetic parasites and their hosts coevolution in simple mathematical models. It is shown that the evolutionary dynamics of a parasite that initially evolves from the host through the loss of the ability to replicate autonomously must substantially differ from that of the host, for a stable host-parasite coevolution regime to be established.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0230-9) contains supplementary material, which is available to authorized users.

LINE- and Alu-containing genomic instability hotspot at 16q24.1 associated with recurrent and nonrecurrent CNV deletions causative for ACDMPV

Transposable elements modify human genome by inserting into new loci or by mediating homology-, microhomology-, or homeology-driven DNA recombination or repair, resulting in genomic structural variation. Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal neonatal developmental lung disorder caused by point mutations or copy-number variant (CNV) deletions of FOXF1 or its distant tissue-specific enhancer. Eighty-five percent of 45 ACDMPV-causative CNV deletions, of which junctions have been sequenced, had at least one of their two breakpoints located in a retrotransposon, with more than half of them being Alu elements. We describe a novel ∼35 kb-large genomic instability hotspot at 16q24.1, involving two evolutionarily young LINE-1 (L1) elements, L1PA2 and L1PA3, flanking AluY, two AluSx, AluSx1, and AluJr elements. The occurrence of L1s at this location coincided with the branching out of the Homo-Pan-Gorilla clade, and was preceded by the insertion of AluSx, AluSx1, and AluJr. Our data show that, in addition to mediating recurrent CNVs, L1 and Alu retrotransposons can predispose the human genome to formation of variably sized CNVs, both of clinical and evolutionary relevance. Nonetheless, epigenetic or other genomic features of this locus might also contribute to its increased instability.

The Mediator co-activator complex regulates Ty1 retromobility by controlling the balance between Ty1i and Ty1 promoters

The Ty1 retrotransposons present in the genome of Saccharomyces cerevisiae belong to the large class of mobile genetic elements that replicate via an RNA intermediary and constitute a significant portion of most eukaryotic genomes. The retromobility of Ty1 is regulated by numerous host factors, including several subunits of the Mediator transcriptional co-activator complex. In spite of its known function in the nucleus, previous studies have implicated Mediator in the regulation of post-translational steps in Ty1 retromobility. To resolve this paradox, we systematically examined the effects of deleting non-essential Mediator subunits on the frequency of Ty1 retromobility and levels of retromobility intermediates. Our findings reveal that loss of distinct Mediator subunits alters Ty1 retromobility positively or negatively over a >10,000-fold range by regulating the ratio of an internal transcript, Ty1i, to the genomic Ty1 transcript. Ty1i RNA encodes a dominant negative inhibitor of Ty1 retromobility that blocks virus-like particle maturation and cDNA synthesis. These results resolve the conundrum of Mediator exerting sweeping control of Ty1 retromobility with only minor effects on the levels of Ty1 genomic RNA and the capsid protein, Gag. Since the majority of characterized intrinsic and extrinsic regulators of Ty1 retromobility do not appear to effect genomic Ty1 RNA levels, Mediator could play a central role in integrating signals that influence Ty1i expression to modulate retromobility.

Retrotransposons are mobile genetic elements that copy their RNA genomes into DNA and insert the DNA copies into the host genome. These elements contribute to genome instability, control of host gene expression and adaptation to changing environments. Retrotransposons depend on numerous host factors for their own propagation and control. The retrovirus-like retrotransposon, Ty1, in the yeast Saccharomyces cerevisiae has been an invaluable model for retrotransposon research, and hundreds of host factors that regulate Ty1 retrotransposition have been identified. Non-essential subunits of the Mediator transcriptional co-activator complex have been identified as one set of host factors implicated in Ty1 regulation. Here, we report a systematic investigation of the effects of loss of these non-essential subunits of Mediator on Ty1 retrotransposition. Our findings reveal a heretofore unknown mechanism by which Mediator influences the balance between transcription from two promoters in Ty1 to modulate expression of an autoinhibitory transcript known as Ty1i RNA. Our results provide new insights into host control of retrotransposon activity via promoter choice and elucidate a novel mechanism by which the Mediator co-activator governs this choice.

The epigenetic alterations of endogenous retroelements in aging

Endogenous retroelements, transposons that mobilize through RNA intermediates, include some of the most abundant repetitive sequences of the human genome, such as Alu and LINE-1 sequences, and human endogenous retroviruses. Recent discoveries demonstrate that these mobile genetic elements not only act as intragenomic parasites, but also exert regulatory roles in living cells. The risk of genomic instability represented by endogenous retroelements is normally counteracted by a series of epigenetic control mechanisms which include, among the most important, CpG DNA methylation. Indeed, most of the genomic CpG sites subjected to DNA methylation in the nuclear DNA are carried by these repetitive elements. As other parts of the genome, endogenous retroelements and other transposable elements are subjected to deep epigenetic alterations during aging, repeatedly observed in the context of organismal and cellular senescence, in human and other species. This review summarizes the current status of knowledge about the epigenetic alterations occurring in this large, non-genic portion of the genome in aging and age-related conditions, with a focus on the causes and the possible functional consequences of these alterations.

Transposable elements shape the human proteome landscape via formation of cis-acting upstream open reading frames

Transposons are major drivers of mammalian genome evolution. To obtain new insights into the contribution of transposons to the regulation of protein translation, we here examined how transposons affected the genesis and function of upstream open reading frames (uORFs), which serve as cis-acting elements to regulate translation from annotated ORFs (anORFs) located downstream of the uORFs in eukaryotic mRNAs. Among 39,786 human uORFs, 3,992 had ATG trinucleotides of a transposon origin, termed "transposon-derived upstream ATGs" or TuATGs. Luciferase reporter assays suggested that many TuATGs modulate translation from anORFs. Comparisons with transposon consensus sequences revealed that most TuATGs were generated by nucleotide substitutions in non-ATG trinucleotides of integrated transposons. Among these non-ATG trinucleotides, GTG and ACG were converted into TuATGs more frequently, indicating a CpG methylation-mediated process of TuATG formation. Interestingly, it is likely that this process accelerated human-specific upstream ATG formation within transposon sequences in 5' untranslated regions after divergence between human and nonhuman primates. Methylation-mediated TuATG formation seems to be ongoing in the modern human population and could alter the expression of disease-related proteins. This study shows that transposons have potentially been shaping the human proteome landscape via cis-acting uORF creation.

Profiling of Human Molecular Pathways Affected by Retrotransposons at the Level of Regulation by Transcription Factor Proteins

Endogenous retroviruses and retrotransposons also termed retroelements (REs) are mobile genetic elements that were active until recently in human genome evolution. REs regulate gene expression by actively reshaping chromatin structure or by directly providing transcription factor binding sites (TFBSs). We aimed to identify molecular processes most deeply impacted by the REs in human cells at the level of TFBS regulation. By using ENCODE data, we identified ~2 million TFBS overlapping with putatively regulation-competent human REs located in 5-kb gene promoter neighborhood (~17% of all TFBS in promoter neighborhoods; ~9% of all RE-linked TFBS). Most of REs hosting TFBS were highly diverged repeats, and for the evolutionary young (0–8% diverged) elements we identified only ~7% of all RE-linked TFBS. The gene-specific distributions of RE-linked TFBS generally correlated with the distributions for all TFBS. However, several groups of molecular processes were highly enriched in the RE-linked TFBS regulation. They were strongly connected with the immunity and response to pathogens, with the negative regulation of gene transcription, ubiquitination, and protein degradation, extracellular matrix organization, regulation of STAT signaling, fatty acids metabolism, regulation of GTPase activity, protein targeting to Golgi, regulation of cell division and differentiation, development and functioning of perception organs and reproductive system. By contrast, the processes most weakly affected by the REs were linked with the conservative aspects of embryo development. We also identified differences in the regulation features by the younger and older fractions of the REs. The regulation by the older fraction of the REs was linked mainly with the immunity, cell adhesion, cAMP, IGF1R, Notch, Wnt, and integrin signaling, neuronal development, chondroitin sulfate and heparin metabolism, and endocytosis. The younger REs regulate other aspects of immunity, cell cycle progression and apoptosis, PDGF, TGF beta, EGFR, and p38 signaling, transcriptional repression, structure of nuclear lumen, catabolism of phospholipids, and heterocyclic molecules, insulin and AMPK signaling, retrograde Golgi-ER transport, and estrogen signaling. The immunity-linked pathways were highly represented in both categories, but their functional roles were different and did not overlap. Our results point to the most quickly evolving molecular pathways in the recent and ancient evolution of human genome.

Activity analysis of LTR12C as an effective regulatory element of the RAE1 gene

Ribonucleic acid export 1 (RAE1) plays an important role in the export of mature mRNAs from the nucleus to the cytoplasm. Long terminal repeats (LTRs) became integrated into the human genome during primate evolution. One such repeat element, LTR12C, lies within a predicted regulatory region located upstream of the RAE1 gene. We examined the transcriptional activity of LTR12C by using the luciferase assay, and showed that the tandem repeat region (TRR) located within LTR12C was required for its regulatory function. A bioinformatics analysis revealed that the LTR12C element had multiple transcription factor binding sites specific for nuclear transcription factor Y (NF-Y), and the promoter activity of LTR12C was significantly decreased after NF-Y knockdown. Additionally, we discovered novel data indicating that LTR12C was initially inserted into the gorilla genome. Taken together, our results reveal that the TRR of LTR12C has powerful regulatory activity due to its NF-Y binding sites, and the integration of the LTR12C element into the primate genome during evolution may have affected RAE1 transcription.

piRNAs and Their Functions in the Brain

Piwi-interacting RNAs (piRNAs) are the non-coding RNAs with 24-32 nucleotides (nt). They exhibit stark differences in length, expression pattern, abundance, and genomic organization when compared to micro-RNAs (miRNAs). There are hundreds of thousands unique piRNA sequences in each species. Numerous piRNAs have been identified and deposited in public databases. Since the piRNAs were originally discovered and well-studied in the germline, a few other studies have reported the presence of piRNAs in somatic cells including neurons. This paper reviewed the common features, biogenesis, functions, and distributions of piRNAs and summarized their specific functions in the brain. This review may provide new insights and research direction for brain disorders.

Small RNA Pathways That Protect the Somatic Genome

Transposable elements (TEs) are DNA elements that can change their position within the genome, with the potential to create mutations and destabilize the genome. As such, special molecular systems have been adopted in animals to control TE activity in order to protect the genome. PIWI proteins, in collaboration with PIWI-interacting RNAs (piRNAs), are well known to play a critical role in silencing germline TEs. Although initially thought to be germline-specific, the role of PIWI–piRNA pathways in controlling TEs in somatic cells has recently begun to be explored in various organisms, together with the role of endogenous small interfering RNAs (endo-siRNAs). This review summarizes recent results suggesting that these small RNA pathways have been critically implicated in the silencing of somatic TEs underlying various physiological traits, with a special focus on the Drosophila model organism.

LINE-1 distribution in six rodent genomes follow a species-specific pattern

L1 distribution in mammal's genomes is yet a huge riddle. However, these repetitive sequences were already found in all chromosomic regions, and in general, they seem to be nonrandomly distributed in the genome. It also seems that after insertion and when they are not deleterious, they are always involved in dynamic processes occurring on that particular chromosomic region. Furthermore, it seems that large-scale genome rearrangements and L1 activity and accumulation are somehow interconnected. In the present study, we analysed L1 genomic distribution in Tatera gambiana (Muridae, Gerbillinae), Acomys sp. (Muridae, Deomyinae), Cricetomys sp. (Nesomyidae, Cricetomyinae), Microtus arvalis (Cricetidae, Arvicolinae), Phodopus roborovskii and P. sungorus (Cricetidae, Cricetinae). All the species studied here seems to exhibit a species-specific pattern.Possible mechanisms, and processes involved in L1 distribution and preferential accumulation in certain regions are di scussed.

The intron-enriched HERV-K(HML-10) family suppresses apoptosis, an indicator of malignant transformation

Background

Human endogenous retroviruses (HERVs) constitute 8% of the human genome and contribute substantially to the transcriptome. HERVs have been shown to generate RNAs that modulate host gene expression. However, experimental evidence for an impact of these regulatory transcripts on the cellular phenotype has been lacking.

Results

We characterized the previously little described HERV-K(HML-10) endogenous retrovirus family on a genome-wide scale. HML-10 invaded the ancestral genome of Old World monkeys about 35 Million years ago and is enriched within introns of human genes when compared to other HERV families. We show that long terminal repeats (LTRs) of HML-10 exhibit variable promoter activity in human cancer cell lines. One identified HML-10 LTR-primed RNA was in opposite orientation to the pro-apoptotic Death-associated protein 3 (DAP3). In HeLa cells, experimental inactivation of HML-10 LTR-primed transcripts induced DAP3 expression levels, which led to apoptosis.

Conclusions

Its enrichment within introns suggests that HML-10 may have been evolutionary co-opted for gene regulation more than other HERV families. We demonstrated such a regulatory activity for an HML-10 RNA that suppressed DAP3-mediated apoptosis in HeLa cells. Since HML-10 RNA appears to be upregulated in various tumor cell lines and primary tumor samples, it may contribute to evasion of apoptosis in malignant cells. However, the overall weak expression of HML-10 transcripts described here raises the question whether our result described for HeLa represent a rare event in cancer. A possible function in other cells or tissues requires further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0081-9) contains supplementary material, which is available to authorized users.

Profuse evolutionary diversification and speciation on volcanic islands: transposon instability and amplification bursts explain the genetic paradox

BACKGROUND

Species-rich adaptive radiations arising from rare plant and animal colonizers are common on remote volcanic archipelagoes. However, they present a paradox. The severe genetic bottleneck of founder events and effects of inbreeding depression, coupled with the inherently stressful volcanic environment, would seem to predict reduced evolutionary potential and increased risk of extinction, rather than rapid adaptive divergence and speciation. Significantly, eukaryotic genomes harbor many families of transposable elements (TEs) that are mobilized by genome shock; these elements may be the primary drivers of genetic reorganization and speciation on volcanic islands.

PRESENTATION OF THE HYPOTHESIS

Here I propose that a central factor in the spectacular radiation and diversification of the endemic Hawaiian Drosophila and other terrestrial lineages on the Hawaiian and other oceanic islands has been repeated bursts of transposition of multiple TEs induced by the unique ecological features of volcanic habitats. Founder individuals and populations on remote volcanic islands experience significant levels of physiological and genomic stress as a consequence of both biotic and abiotic factors. This results in disruption of the usual epigenetic suppression of TEs, unleashing them to proliferate and spread, which in turn gives rise to novel genetic variation and remodels genomic regulatory circuits, facilitating rapid morphological, ecological and behavioral change, and adaptive radiation.

TESTING THE HYPOTHESIS

To obtain empirical support for the hypothesis, test organisms should be exposed to prolonged heat stress, high levels of carbon dioxide and other volcanic gases, along with inbreeding. Data from subsequent whole genome sequencing and bioinformatics screening for TE numbers and locations would then be compared with initial pre-exposure TE information for the test strains, a labor-intensive project. Several predicted outcomes arising from the hypothesis are discussed. Currently available data are consistent with the proposed concept of stress-induced TE mobilization as a trigger of evolutionary diversification and speciation on volcanic islands.

IMPLICATIONS OF THE HYPOTHESIS

The main implication is that both TEs and species should proliferate at a much higher rate on volcanic islands than elsewhere. Second, the evolvability of a lineage may correlate with the abundance and distribution of TEs in the genome. Successful colonizers of volcanic habitats with high genomic proportions of TEs may be best poised to found a speciose lineage that gives rise to a dramatic adaptive radiation. Colonizers that are depauperate in TEs are likely to be evolutionarily constrained and diversify little, if at all.

REVIEWERS

This article was reviewed by Dr. James Shapiro and Dr. Wolfgang Miller (nominated by Editorial Board member Dr. I. King Jordan).

Evidence of extensive non-allelic gene conversion among LTR elements in the human genome

Long Terminal Repeats (LTRs) are nearly identical DNA sequences found at either end of Human Endogenous Retroviruses (HERVs). The high sequence similarity that exists among different LTRs suggests they could be substrate of ectopic gene conversion events. To understand the extent to which gene conversion occurs and to gain new insights into the evolutionary history of these elements in humans, we performed an intra-species phylogenetic study of 52 LTRs on different unrelated Y chromosomes. From this analysis, we obtained direct evidence that demonstrates the occurrence of ectopic gene conversion in several LTRs, with donor sequences located on both sex chromosomes and autosomes. We also found that some of these elements are characterized by an extremely high density of polymorphisms, showing one of the highest nucleotide diversities in the human genome, as well as a complex patchwork of sequences derived from different LTRs. Finally, we highlighted the limits of current short-read NGS studies in the analysis of genetic diversity of the LTRs in the human genome. In conclusion, our comparative re-sequencing analysis revealed that ectopic gene conversion is a common event in the evolution of LTR elements, suggesting complex genetic links among LTRs from different chromosomes.

[STRUCTURE AND DISTRIBUTION OF THE RETROTRANSPOSON BOV-B LINE]

The classification of mobile elements was discussed. Special attention was devoted to the retroelement of the LINE group: retrotransposon Bov-B LINE. The history of its origin and distribution in the nature was considered. The results of the phenomenon of horizontal transition of the retrotransposon Bov-B LINE between evolutionally distant classes were discussed.

New bioinformatic tool for quick identification of functionally relevant endogenous retroviral inserts in human genome

Endogenous retroviruses (ERVs) and LTR retrotransposons (LRs) occupy ∼8% of human genome. Deep sequencing technologies provide clues to understanding of functional relevance of individual ERVs/LRs by enabling direct identification of transcription factor binding sites (TFBS) and other landmarks of functional genomic elements. Here, we performed the genome-wide identification of human ERVs/LRs containing TFBS according to the ENCODE project. We created the first interactive ERV/LRs database that groups the individual inserts according to their familial nomenclature, number of mapped TFBS and divergence from their consensus sequence. Information on any particular element can be easily extracted by the user. We also created a genome browser tool, which enables quick mapping of any ERV/LR insert according to genomic coordinates, known human genes and TFBS. These tools can be used to easily explore functionally relevant individual ERV/LRs, and for studying their impact on the regulation of human genes. Overall, we identified ∼110,000 ERV/LR genomic elements having TFBS. We propose a hypothesis of “domestication” of ERV/LR TFBS by the genome milieu including subsequent stages of initial epigenetic repression, partial functional release, and further mutation-driven reshaping of TFBS in tight coevolution with the enclosing genomic loci.

Antisense Transcription of Retrotransposons in Drosophila: An Origin of Endogenous Small Interfering RNA Precursors

Movement of transposons causes insertions, deletions, and chromosomal rearrangements potentially leading to premature lethality in Drosophila melanogaster. To repress these elements and combat genomic instability, eukaryotes have evolved several small RNA-mediated defense mechanisms. Specifically, in Drosophila somatic cells, endogenous small interfering (esi)RNAs suppress retrotransposon mobility. EsiRNAs are produced by Dicer-2 processing of double-stranded RNA precursors, yet the origins of these precursors are unknown. We show that most transposon families are transcribed in both the sense (S) and antisense (AS) direction in Dmel-2 cells. LTR retrotransposons Dm297, mdg1, and blood, and non-LTR retrotransposons juan and jockey transcripts, are generated from intraelement transcription start sites with canonical RNA polymerase II promoters. We also determined that retrotransposon antisense transcripts are less polyadenylated than sense. RNA-seq and small RNA-seq revealed that Dicer-2 RNA interference (RNAi) depletion causes a decrease in the number of esiRNAs mapping to retrotransposons and an increase in expression of both S and AS retrotransposon transcripts. These data support a model in which double-stranded RNA precursors are derived from convergent transcription and processed by Dicer-2 into esiRNAs that silence both sense and antisense retrotransposon transcripts. Reduction of sense retrotransposon transcripts potentially lowers element-specific protein levels to prevent transposition. This mechanism preserves genomic integrity and is especially important for Drosophila fitness because mobile genetic elements are highly active.

Gain of a New Exon by a Lineage-Specific Alu Element-Integration Event in the BCS1L Gene during Primate Evolution

BCS1L gene encodes mitochondrial protein and is a member of conserved AAA protein family. This gene is involved in the incorporation of Rieske FeS and Qcr10p into complex III of respiratory chain. In our previous study, AluYRa2-derived alternative transcript in rhesus monkey genome was identified. However, this transcript has not been reported in human genome. In present study, we conducted evolutionary analysis of AluYRa2-exonized transcript with various primate genomic DNAs and cDNAs from humans, rhesus monkeys, and crab-eating monkeys. Remarkably, our results show that AluYRa2 element has only been integrated into genomes of Macaca species. This Macaca lineage-specific integration of AluYRa2 element led to exonization event in the first intron region of BCS1L gene by producing a conserved 3' splice site. Intriguingly, in rhesus and crab-eating monkeys, more diverse transcript variants by alternative splicing (AS) events, including exon skipping and different 5' splice sites from humans, were identified. Alignment of amino acid sequences revealed that AluYRa2-exonized transcript has short N-terminal peptides. Therefore, AS events play a major role in the generation of various transcripts and proteins during primate evolution. In particular, lineage-specific integration of Alu elements and species-specific Alu-derived exonization events could be important sources of gene diversification in primates.

Molecular functions of human endogenous retroviruses in health and disease

Human endogenous retroviruses (HERVs) and related genetic elements form 504 distinct families and occupy ~8% of human genome. Recent success of high-throughput experimental technologies facilitated understanding functional impact of HERVs for molecular machinery of human cells. HERVs encode active retroviral proteins, which may exert important physiological functions in the body, but also may be involved in the progression of cancer and numerous human autoimmune, neurological and infectious diseases. The spectrum of related malignancies includes, but not limits to, multiple sclerosis, psoriasis, lupus, schizophrenia, multiple cancer types and HIV. In addition, HERVs regulate expression of the neighboring host genes and modify genomic regulatory landscape, e.g., by providing regulatory modules like transcription factor binding sites (TFBS). Indeed, recent bioinformatic profiling identified ~110,000 regulatory active HERV elements, which formed at least ~320,000 human TFBS. These and other peculiarities of HERVs might have played an important role in human evolution and speciation. In this paper, we focus on the current progress in understanding of normal and pathological molecular niches of HERVs, on their implications in human evolution, normal physiology and disease. We also review the available databases dealing with various aspects of HERV genetics.

Regulation of gene expression via retrotransposon insertions and the noncoding RNA 4.5S RNAH

Short interspersed elements (SINEs) comprise a significant portion of mammalian genomes and regulate gene expression through a variety of mechanisms. Here, we show that Myodonta clade-specific 4.5S RNAH (4.5SH), an abundant nuclear noncoding RNA that is highly homologous to the retrotransposon SINE B1, controls the expression of reporter gene that contains the antisense insertion of SINE B1 via nuclear retention. The depletion of endogenous 4.5SH with antisense oligonucleotides neutralizes the nuclear retention and changes the subcellular distribution of the reporter transcripts containing the antisense SINE B1 insertion. Importantly, endogenous transcripts with antisense SINE B1 were increased in the cytoplasm after knockdown of 4.5SH, leading to a decrease in cellular growth. We propose a tentative hypothesis that the amplification of the 4.5SH cluster in specific rodent species might delineate their evolutionary direction via the regulation of genes containing the antisense insertion of SINE B1.

A Deluge of Complex Repeats: The Solanum Genome

Repetitive elements have lately emerged as key components of genome, performing varieties of roles. It has now become necessary to have an account of repeats for every genome to understand its dynamics and state. Recently, genomes of two major Solanaceae species, Solanum tuberosum and Solanum lycopersicum, were sequenced. These species are important crops having high commercial significance as well as value as model species. However, there is a reasonable gap in information about repetitive elements and their possible roles in genome regulation for these species. The present study was aimed at detailed identification and characterization of complex repetitive elements in these genomes, along with study of their possible functional associations as well as to assess possible transcriptionally active repetitive elements. In this study, it was found that ~50-60% of genomes of S. tuberosum and S. lycopersicum were composed of repetitive elements. It was also found that complex repetitive elements were associated with >95% of genes in both species. These two genomes are mostly composed of LTR retrotransposons. Two novel repeat families very similar to LTR/ERV1 and LINE/RTE-BovB have been reported for the first time. Active existence of complex repeats was estimated by measuring their transcriptional abundance using Next Generation Sequencing read data and Microarray platforms. A reasonable amount of regulatory components like transcription factor binding sites and miRNAs appear to be under the influence of these complex repetitive elements in these species, while several genes appeared to possess exonized repeats.

Transposable element-driven transcript diversification and its relevance to genetic disorders

The human genome project and subsequent gene annotation projects have shown that the human genome contains 22,000-25,000 functional genes. Therefore, it is believed that the diversity of protein repertoire is achieved by the alternative splicing (AS) mechanism. Transposable elements (TEs) are mobile in nature and can therefore alter their position in the genome. The insertion of TEs into a new gene region can result in AS of a particular transcript through various mechanisms, including intron retention, and alternative donor or acceptor splice sites. TE-derived AS is thought to have played a part in primate evolution and in hominid radiation. However, TE-derived AS or genetic instability may sometimes result in genetic disorders. For the past two decades, numerous studies have been performed on TEs and their role in genomes. Accumulating evidence shows that the term 'junk DNA', previously used for TEs is a misnomer. Recent research has indicated that TEs may have clinical potential. However, to explore the feasibility of using TEs in clinical practice, additional studies are required. This review summarizes the available literature on TE-derived AS, alternative promoter, and alternative polyadenylation. The review covers the effects of TEs on coding genes and their clinical implications, and provides our perspectives and directions for future research.

In-vivo fusion of human cancer and hamster stromal cells permanently transduces and transcribes human DNA

After demonstrating, with karyotyping, polymerase chain reaction (PCR) and fluorescence in-situ hybridization, the retention of certain human chromosomes and genes following the spontaneous fusion of human tumor and hamster cells in-vivo, it was postulated that cell fusion causes the horizontal transmission of malignancy and donor genes. Here, we analyzed gene expression profiles of 3 different hybrid tumors first generated in the hamster cheek pouch after human tumor grafting, and then propagated in hamsters and in cell cultures for years: two Hodgkin lymphomas (GW-532, GW-584) and a glioblastoma multiforme (GB-749). Based on the criteria of MAS 5.0 detection P-values ≤0.065 and at least a 2-fold greater signal expression value than a hamster melanoma control, we identified 3,759 probe sets (ranging from 1,040 to 1,303 in each transplant) from formalin-fixed, paraffin-embedded sections of the 3 hybrid tumors, which unambiguously mapped to 3,107 unique Entrez Gene IDs, representative of all human chromosomes; however, by karyology, one of the hybrid tumors (GB-749) had a total of 15 human chromosomes in its cells. Among the genes mapped, 39 probe sets, representing 33 unique Entrez Gene IDs, complied with the detection criteria in all hybrid tumor samples. Five of these 33 genes encode transcription factors that are known to regulate cell growth and differentiation; five encode cell adhesion- and transmigration-associated proteins that participate in oncogenesis and/or metastasis and invasion; and additional genes encode proteins involved in signaling pathways, regulation of apoptosis, DNA repair, and multidrug resistance. These findings were corroborated by PCR and reverse transcription PCR, showing the presence of human alphoid (α)-satellite DNA and the F11R transcripts in additional tumor transplant generations. We posit that in-vivo fusion discloses genes implicated in tumor progression, and gene families coding for the organoid phenotype. Thus, cancer cells can transduce adjacent stromal cells, with the resulting progeny having permanently transcribed genes with malignant and other gene functions of the donor DNA. Using heterospecific in-vivo cell fusion, genes encoding oncogenic and organogenic traits may be identified.

Alu elements: at the crossroads between disease and evolution

The cost of DNA sequencing is decreasing year by year, and the era of personalized medicine and the $1000 genome seems to be just around the corner. In order to link genetic variation to gene function, however, we need to learn more about the function of the non-coding genomic elements. The advance of high-throughput sequencing enabled rapid progress in mapping the functional elements in our genome. In the present article, I discuss how intronic mutations acting at Alu elements enable formation of new exons. I review the mutations that cause disease when promoting a major increase in the inclusion of Alu exon into mature transcripts. Moreover, I present the mechanism that represses such a major inclusion of Alu exons and instead enables a gradual evolution of Alu elements into new exons.

Human endogenous retroviruses and ADHD

OBJECTIVES

Several lines of evidences suggest that human endogenous retroviruses (HERVs) are implicated in the development of many complex diseases with a multifactorial aetiology and a strong heritability, such as neurological and psychiatric diseases. Attention deficit hyperactivity Disorder (ADHD) is a neurodevelopmental disorder that results from a complex interaction of environmental, biological and genetic factors. Our aim was to analyse the expression levels of three HERV families (HERV-H, K and W) in patients with ADHD.

METHODS

The expression of retroviral mRNAs from the three HERV families was evaluated in peripheral blood mononuclear cells (PBMCs) from 30 patients with ADHD and 30 healthy controls by quantitative RT-PCR.

RESULTS

The expression levels of HERV-H are significantly higher in patients with ADHD compared to healthy controls, while there are no differences in the expression levels of HERV-K and W.

CONCLUSIONS

Since the ADHD aetiology is due to a complex interaction of environmental, biological and genetic factors, HERVs may represent one link among these factors and clinical phenotype of ADHD. A future confirmation of HERV-H overexpression in a larger number of ADHD patients will make possible to identify it as a new parameter for this clinical condition, also contributing to deepen the study on the role of HERVs in the neurodevelopment diseases.

Identification of structural variation in mouse genomes

Structural variation is variation in structure of DNA regions affecting DNA sequence length and/or orientation. It generally includes deletions, insertions, copy-number gains, inversions, and transposable elements. Traditionally, the identification of structural variation in genomes has been challenging. However, with the recent advances in high-throughput DNA sequencing and paired-end mapping (PEM) methods, the ability to identify structural variation and their respective association to human diseases has improved considerably. In this review, we describe our current knowledge of structural variation in the mouse, one of the prime model systems for studying human diseases and mammalian biology. We further present the evolutionary implications of structural variation on transposable elements. We conclude with future directions on the study of structural variation in mouse genomes that will increase our understanding of molecular architecture and functional consequences of structural variation.

Physiology of the read-write genome

Discoveries in cytogenetics, molecular biology, and genomics have revealed that genome change is an active cell-mediated physiological process. This is distinctly at variance with the pre-DNA assumption that genetic changes arise accidentally and sporadically. The discovery that DNA changes arise as the result of regulated cell biochemistry means that the genome is best modelled as a read-write (RW) data storage system rather than a read-only memory (ROM). The evidence behind this change in thinking and a consideration of some of its implications are the subjects of this article. Specific points include the following: cells protect themselves from accidental genome change with proofreading and DNA damage repair systems; localized point mutations result from the action of specialized trans-lesion mutator DNA polymerases; cells can join broken chromosomes and generate genome rearrangements by non-homologous end-joining (NHEJ) processes in specialized subnuclear repair centres; cells have a broad variety of natural genetic engineering (NGE) functions for transporting, diversifying and reorganizing DNA sequences in ways that generate many classes of genomic novelties; natural genetic engineering functions are regulated and subject to activation by a range of challenging life history events; cells can target the action of natural genetic engineering functions to particular genome locations by a range of well-established molecular interactions, including protein binding with regulatory factors and linkage to transcription; and genome changes in cancer can usefully be considered as consequences of the loss of homeostatic control over natural genetic engineering functions.

Genomic localization of AtRE1 and AtRE2, copia-type retrotransposons, in natural variants of Arabidopsis thaliana

Retrotransposons are ubiquitous components of plant genomes. They affect genome organization, and can also affect the expression patterns of neighboring genes. Retrotransposons are therefore important elements for changing genomic information. To understand the evolution of the Arabidopsis genome, we examined the distribution of certain retrotransposons, AtRE1s and AtRE2s, in the genomes of 12 natural variants (accessions) of Arabidopsis thaliana. AtRE1 and AtRE2 are copia-type retrotransposons that are potentially active. Their copy numbers are low, and they are absent from the genomes of some accessions. We detected four loci with AtRE1s inserted in six accessions, and one locus with an insertion of a solo-LTR-like sequence derived from AtRE1 in two accessions. Seven loci with AtRE2s inserted were detected on eight accessions. These loci were distributed in euchromatic regions of chromosomes 1, 2, 3, and 4. The AtRE1 and AtRE2 sequences at some loci identified in this study have not been recorded in the database of the 1001 Genome project. The sequences of AtRE1s and those of AtRE2s in different accessions and at different loci were highly conserved. There was a complete or almost complete conservation of sequences of both long terminal repeats in each AtRE1 and in each AtRE2. These results suggest that AtRE1 and AtRE2 appeared quite recently in the Arabidopsis genome. Furthermore, sequence comparisons of AtRE1 and AtRE2 loci among accessions revealed the possibility that large deletions containing entire sequences of AtRE1 and AtRE2 have occurred in some accessions.