An evolutionary classification of genomic function

Pseudogenes in plasmid genomes reveal past transitions in plasmid mobility

Evidence for gene non-functionalization due to mutational processes is found in genomes in the form of pseudogenes. Pseudogenes are known to be rare in prokaryote chromosomes, with the exception of lineages that underwent an extreme genome reduction (e.g. obligatory symbionts). Much less is known about the frequency of pseudogenes in prokaryotic plasmids; those are genetic elements that can transfer between cells and may encode beneficial traits for their host. Non-functionalization of plasmid-encoded genes may alter the plasmid characteristics, e.g. mobility, or their effect on the host. Analyzing 10 832 prokaryotic genomes, we find that plasmid genomes are characterized by threefold-higher pseudogene density compared to chromosomes. The majority of plasmid pseudogenes correspond to deteriorated transposable elements. A detailed analysis of enterobacterial plasmids furthermore reveals frequent gene non-functionalization events associated with the loss of plasmid self-transmissibility. Reconstructing the evolution of closely related plasmids reveals that non-functionalization of the conjugation machinery led to the emergence of non-mobilizable plasmid types. Examples are virulence plasmids in Escherichia and Salmonella. Our study highlights non-functionalization of core plasmid mobility functions as one route for the evolution of domesticated plasmids. Pseudogenes in plasmids supply insights into past transitions in plasmid mobility that are akin to transitions in bacterial lifestyle.

Bridging multiple dimensions: roles of transposable elements in higher-order genome regulation

Transposable elements (TEs) such as endogenous retroviruses (ERVs), long interspersed nuclear elements (LINEs), and short interspersed nuclear elements (SINEs) occupy nearly half of typical mammalian genomes. Previous studies show that these parasitic elements, especially LINEs and ERVs, provide important activities promoting host germ cell and placental development, preimplantation embryogenesis, and maintenance of pluripotent stem cells. Despite being the most numerically abundant type of TEs in the genome, the consequences of SINEs on host genome regulation are less well characterized than those of ERVs and LINEs. Interestingly, recent findings reveal that SINEs recruit the key architectural protein CTCF (CCCTC-binding factor), indicating a role of these elements for 3D genome regulation. Higher-order nuclear structures are linked with important cellular functions such as gene regulation and DNA replication. SINEs and other TEs, therefore, may mediate distinct physiological processes with benefits to the host by modulating the 3D genome.

Moderation of Structural DNA Properties by Coupled Dinucleotide Contents in Eukaryotes

Dinucleotides are known as determinants for various structural and physiochemical properties of DNA and for binding affinities of proteins to DNA. These properties (e.g., stiffness) and bound proteins (e.g., transcription factors) are known to influence important biological functions, such as transcription regulation and 3D chromatin organization. Accordingly, the question arises of how the considerable variations in dinucleotide contents of eukaryotic chromosomes could still provide consistent DNA properties resulting in similar functions and 3D conformations. In this work, we investigate the hypothesis that coupled dinucleotide contents influence DNA properties in opposite directions to moderate each other's influences. Analyzing all 2478 chromosomes of 155 eukaryotic species, considering bias from coding sequences and enhancers, we found sets of correlated and anti-correlated dinucleotide contents. Using computational models, we estimated changes of DNA properties resulting from this coupling. We found that especially pure A/T dinucleotides (AA, TT, AT, TA), known to influence histone positioning and AC/GT contents, are relevant moderators and that, e.g., the Roll property, which is known to influence histone affinity of DNA, is preferably moderated. We conclude that dinucleotide contents might indirectly influence transcription and chromatin 3D conformation, via regulation of histone occupancy and/or other mechanisms.

Non-coding RNAs: the silent regulators of health and diseases

Non-coding RNAs (ncRNAs) like miRNAs, siRNA, lncRNAs, circRNAs, piRNAs, snoRNAs, snRNAs etc. form a collective group of RNAs that is instrumental to the various functions of the genome. With the advent of cutting-edge molecular biology tools and techniques, scientists have unearthed several mechanisms through which these ncRNAs act. Although our understanding may still be limited, yet scientists have been able to establish ncRNAs as major regulators of genetic inter-plays that dictate various pathophysiological conditions. This special issue of Molecular Biology Reports features a collection of research and review articles on ncRNAs and their involvement in different pathophysiological conditions that include different types of cancers. It is expected that this special issue will motivate researchers in the field to delve deeper into the world of ncRNAs and attempt to develop new diagnostic and therapeutic interventions for challenging clinical conditions.

The state of Medusozoa genomics: current evidence and future challenges

Medusozoa is a widely distributed ancient lineage that harbors one-third of Cnidaria diversity divided into 4 classes. This clade is characterized by the succession of stages and modes of reproduction during metagenic lifecycles, and includes some of the most plastic body plans and life cycles among animals. The characterization of traditional genomic features, such as chromosome numbers and genome sizes, was rather overlooked in Medusozoa and many evolutionary questions still remain unanswered. Modern genomic DNA sequencing in this group started in 2010 with the publication of the Hydra vulgaris genome and has experienced an exponential increase in the past 3 years. Therefore, an update of the state of Medusozoa genomics is warranted. We reviewed different sources of evidence, including cytogenetic records and high-throughput sequencing projects. We focused on 4 main topics that would be relevant for the broad Cnidaria research community: (i) taxonomic coverage of genomic information; (ii) continuity, quality, and completeness of high-throughput sequencing datasets; (iii) overview of the Medusozoa specific research questions approached with genomics; and (iv) the accessibility of data and metadata. We highlight a lack of standardization in genomic projects and their reports, and reinforce a series of recommendations to enhance future collaborative research.

Taming, Domestication and Exaptation: Trajectories of Transposable Elements in Genomes

During evolution, several types of sequences pass through genomes. Along with mutations and internal genetic tinkering, they are a useful source of genetic variability for adaptation and evolution. Most of these sequences are acquired by horizontal transfers (HT), but some of them may come from the genomes themselves. If they are not lost or eliminated quickly, they can be tamed, domesticated, or even exapted. Each of these processes results from a series of events, depending on the interactions between these sequences and the host genomes, but also on environmental constraints, through their impact on individuals or population fitness. After a brief reminder of the characteristics of each of these states (taming, domestication, exaptation), the evolutionary trajectories of these new or acquired sequences will be presented and discussed, emphasizing that they are not totally independent insofar as the first can constitute a step towards the second, and the second is another step towards the third.

Major Evolutionary Transitions and the Roles of Facilitation and Information in Ecosystem Transformations

A small number of extraordinary “Major Evolutionary Transitions” (METs) have attracted attention among biologists. They comprise novel forms of individuality and information, and are defined in relation to organismal complexity, irrespective of broader ecosystem-level effects. This divorce between evolutionary and ecological consequences qualifies unicellular eukaryotes, for example, as a MET although they alone failed to significantly alter ecosystems. Additionally, this definition excludes revolutionary innovations not fitting into either MET type (e.g., photosynthesis). We recombine evolution with ecology to explore how and why entire ecosystems were newly created or radically altered – as Major System Transitions (MSTs). In doing so, we highlight important morphological adaptations that spread through populations because of their immediate, direct-fitness advantages for individuals. These are Major Competitive Transitions, or MCTs. We argue that often multiple METs and MCTs must be present to produce MSTs. For example, sexually-reproducing, multicellular eukaryotes (METs) with anisogamy and exoskeletons (MCTs) significantly altered ecosystems during the Cambrian. Therefore, we introduce the concepts of Facilitating Evolutionary Transitions (FETs) and Catalysts as key events or agents that are insufficient themselves to set a MST into motion, but are essential parts of synergies that do. We further elucidate the role of information in MSTs as transitions across five levels: (I) Encoded; (II) Epigenomic; (III) Learned; (IV) Inscribed; and (V) Dark Information. The latter is ‘authored’ by abiotic entities rather than biological organisms. Level IV has arguably allowed humans to produce a MST, and V perhaps makes us a FET for a future transition that melds biotic and abiotic life into one entity. Understanding the interactive processes involved in past major transitions will illuminate both current events and the surprising possibilities that abiotically-created information may produce.

The Genomics of Plant Satellite DNA

The twenty-first century began with a certain indifference to the research of satellite DNA (satDNA). Neither genome sequencing projects were able to accurately encompass the study of satDNA nor classic methodologies were able to go further in undertaking a better comprehensive study of the whole set of satDNA sequences of a genome. Nonetheless, knowledge of satDNA has progressively advanced during this century with the advent of new analytical techniques. The enormous advantages that genome-wide approaches have brought to its analysis have now stimulated a renewed interest in the study of satDNA. At this point, we can look back and try to assess more accurately many of the key questions that were left unsolved in the past about this enigmatic and important component of the genome. I review here the understanding gathered on plant satDNAs over the last few decades with an eye on the near future.

The role of purifying selection in the origin and maintenance of complex function

Fitness contribution alone should not be the criterion of 'function' in molecular biology and genomics. Disagreement over the use of 'function' in molecular biology and genomics is still with us, almost eight years after publicity surrounding the Encyclopedia of DNA Elements project claimed that 80.4% of the human genome comprises "functional elements". Recent approaches attempt to resolve or reformulate this debate by redefining genomic 'function' in terms of current fitness contribution. In its favour, this redefinition for the genomic context is in apparent conformity with predominant experimental practices, especially in biomedical research, and with ascription of function by selective maintenance. We argue against approaches of this kind, however, on the grounds that they could be seen as non-Darwinian, and fail to properly account for the diversity of non-adaptive processes involved in the origin and maintenance of genomic complexity. We examine cases of molecular and organismal complexity that arise neutrally, showing how purifying selection maintains non-adaptive genomic complexity. Rather than lumping different sorts of genomic complexity together by defining 'function' as fitness contribution, we argue that it is best to separate the heterogeneous contributions of preaptation, exaptation and adaptation to the historical processes of origin and maintenance for complex features.

ConsRM: collection and large-scale prediction of the evolutionarily conserved RNA methylation sites, with implications for the functional epitranscriptome

Motivation N6-methyladenosine (m6A) is the most prevalent RNA modification on mRNAs and lncRNAs. Evidence increasingly demonstrates its crucial importance in essential molecular mechanisms and various diseases. With recent advances in sequencing techniques, tens of thousands of m6A sites are identified in a typical high-throughput experiment, posing a key challenge to distinguish the functional m6A sites from the remaining 'passenger' (or 'silent') sites. Results: We performed a comparative conservation analysis of the human and mouse m6A epitranscriptomes at single site resolution. A novel scoring framework, ConsRM, was devised to quantitatively measure the degree of conservation of individual m6A sites. ConsRM integrates multiple information sources and a positive-unlabeled learning framework, which integrated genomic and sequence features to trace subtle hints of epitranscriptome layer conservation. With a series validation experiments in mouse, fly and zebrafish, we showed that ConsRM outperformed well-adopted conservation scores (phastCons and phyloP) in distinguishing the conserved and unconserved m6A sites. Additionally, the m6A sites with a higher ConsRM score are more likely to be functionally important. An online database was developed containing the conservation metrics of 177 998 distinct human m6A sites to support conservation analysis and functional prioritization of individual m6A sites. And it is freely accessible at: https://www.xjtlu.edu.cn/biologicalsciences/con.

Plant Mitochondria are a Riddle Wrapped in a Mystery Inside an Enigma

A fundamental paradox motivates the study of plant mitochondrial genomics: the mutation rate is very low (lower than in the nucleus) but the rearrangement rate is high. A landmark paper published in Journal of Molecular Evolution in 1988 established these facts and revealed the paradox. Jeffrey Palmer and Laura Herbon did a prodigious amount of work in the pre-genome sequencing era to identify both the high frequency of rearrangements between closely related species, and the low frequency of mutations, observations that have now been confirmed many times by sequencing. This paper was also the first to use molecular data on rearrangements as a phylogenetic trait to build a parsimonious tree. The work was a technical tour-de-force, its findings are still at the heart of plant mitochondrial genomics, and the underlying molecular mechanisms that produce this paradox are still not completely understood.

Dark Matter of Primate Genomes: Satellite DNA Repeats and Their Evolutionary Dynamics

A substantial portion of the primate genome is composed of non-coding regions, so-called "dark matter", which includes an abundance of tandemly repeated sequences called satellite DNA. Collectively known as the satellitome, this genomic component offers exciting evolutionary insights into aspects of primate genome biology that raise new questions and challenge existing paradigms. A complete human reference genome was recently reported with telomere-to-telomere human X chromosome assembly that resolved hundreds of dark regions, encompassing a 3.1 Mb centromeric satellite array that had not been identified previously. With the recent exponential increase in the availability of primate genomes, and the development of modern genomic and bioinformatics tools, extensive growth in our knowledge concerning the structure, function, and evolution of satellite elements is expected. The current state of knowledge on this topic is summarized, highlighting various types of primate-specific satellite repeats to compare their proportions across diverse lineages. Inter- and intraspecific variation of satellite repeats in the primate genome are reviewed. The functional significance of these sequences is discussed by describing how the transcriptional activity of satellite repeats can affect gene expression during different cellular processes. Sex-linked satellites are outlined, together with their respective genomic organization. Mechanisms are proposed whereby satellite repeats might have emerged as novel sequences during different evolutionary phases. Finally, the main challenges that hinder the detection of satellite DNA are outlined and an overview of the latest methodologies to address technological limitations is presented.

Functional Long Non-coding RNAs Evolve from Junk Transcripts

Transcriptome studies reveal pervasive transcription of complex genomes, such as those of mammals. Despite popular arguments for functionality of most, if not all, of these transcripts, genome-wide analysis of selective constraints indicates that most of the produced RNA are junk. However, junk is not garbage. On the contrary, junk transcripts provide the raw material for the evolution of diverse long non-coding (lnc) RNAs by non-adaptive mechanisms, such as constructive neutral evolution. The generation of many novel functional entities, such as lncRNAs, that fuels organismal complexity does not seem to be driven by strong positive selection. Rather, the weak selection regime that dominates the evolution of most multicellular eukaryotes provides ample material for functional innovation with relatively little adaptation involved.

Genome size evolution: towards new model systems for old questions

Genome size (GS) variation is a fundamental biological characteristic; however, its evolutionary causes and consequences are the topic of ongoing debate. Whether GS is a neutral trait or one subject to selective pressures, and how strong these selective pressures are, may remain open questions. Fundamentally, the genomic sequences responsible for this variation directly impact the potential evolutionary outcomes and, equally, are the targets of different evolutionary pressures. For example, duplications and deletions of genic regions (large or small) can have immediate and drastic phenotypic effects, while an expansion or contraction of non-coding DNA is less likely to cause catastrophic phenotypic effects. However, in the long term, the accumulation or deletion of ncDNA is likely to have larger effects. Modern sequencing technologies are allowing for the dissection of these proximate causes, but a combination of these new technologies with more traditional evolutionary experiments and approaches could revolutionize this debate and potentially resolve many of these arguments. Here, I discuss an ambitious way forward for GS research, putting it in context of historical debates, theories and sometimes contradictory evidence, and highlighting the promise of combining new sequencing technologies and analytical developments with more traditional experimental evolution approaches.

Are Antisense Proteins in Prokaryotes Functional?

Many prokaryotic RNAs are transcribed from loci outside of annotated protein coding genes. Across bacterial species hundreds of short open reading frames antisense to annotated genes show evidence of both transcription and translation, for instance in ribosome profiling data. Determining the functional fraction of these protein products awaits further research, including insights from studies of molecular interactions and detailed evolutionary analysis. There are multiple lines of evidence, however, that many of these newly discovered proteins are of use to the organism. Condition-specific phenotypes have been characterized for a few. These proteins should be added to genome annotations, and the methods for predicting them standardized. Evolutionary analysis of these typically young sequences also may provide important insights into gene evolution. This research should be prioritized for its exciting potential to uncover large numbers of novel proteins with extremely diverse potential practical uses, including applications in synthetic biology and responding to pathogens.

Postgenomics function monism

The ENCODE project has made important new estimates of human genome functionality, now revising the percentage considered functional to more than 80%, which is in stark contrast to the received view, which estimated that less than 10% of the conserved parts of the human genome are functional. ENCODE's unorthodox use of the notion of biological function has stirred the so-called ENCODE controversy, involving conflicting views about the correct notion of function in postgenomics. The debate hinges on the traditional philosophical contrast between the causal role (CR) and selected effects (SE) approaches. In this paper, we examine the ENCODE controversy in terms of the distinction between function monism and pluralism. We propose to apply a weak etiological account to genomic function ascriptions. In this approach, we can ascribe a function to a genomic structure of an organism if and only if performing the function persists in causally contributing to the organism's and its ancestors' fitness. In comparison to the strong etiological (i.e., the selected effects) approach, the present account does not require there to be selection for the structure in question. This is a monistic approach that enables us to avoid the main difficulties of CR, as well as SE's overdependence on natural selection, while still preserving an evolutionary-constrained notion of biological functions. Our proposal is much more moderate in accommodating the estimates of the functionality of the human genome than both ENCODE's proposal itself and the views of the critics relying on a version of the SE account of functions.

Amazing symmetrical clustering in chloroplast genomes

Background

Previously, a seven-cluster pattern claiming to be a universal one in bacterial genomes has been reported. Keeping in mind the most popular theory of chloroplast origin, we checked whether a similar pattern is observed in chloroplast genomes.

Results

Surprisingly, eight cluster structure has been found, for chloroplasts. The pattern observed for chloroplasts differs rather significantly, from bacterial one, and from that latter observed for cyanobacteria. The structure is provided by clustering of the fragments of equal length isolated within a genome so that each fragment is converted in triplet frequency dictionary with non-overlapping triplets with no gaps in frame tiling. The points in 63-dimensional space were clustered due to elastic map technique. The eight cluster found in chloroplasts comprises the fragments of a genome bearing tRNA genes and exhibiting excessively high GC-content, in comparison to the entire genome.

Conclusion

Chloroplasts exhibit very specific symmetry type in distribution of coding and non-coding fragments of a genome in the space of triplet frequencies: this is mirror symmetry. Cyanobacteria may have both mirror symmetry, and the rotational symmetry typical for other bacteria.

The meanings of 'function' in biology and the problematic case of de novo gene emergence

The word function has many different meanings in molecular biology. Here we explore the use of this word (and derivatives like functional) in research papers about de novo gene birth. Based on an analysis of 20 abstracts we propose a simple lexicon that, we believe, will help scientists and philosophers discuss the meaning of function more clearly.

Deep Conservation of cis-Element Variants Regulating Plant Hormonal Responses

Phytohormones regulate many aspects of plant life by activating transcription factors (TFs) that bind sequence-specific response elements (REs) in regulatory regions of target genes. Despite their short length, REs are degenerate, with a core of just 3 to 4 bp. This degeneracy is paradoxical, as it reduces specificity and REs are extremely common in the genome. To study whether RE degeneracy might serve a biological function, we developed an algorithm for the detection of regulatory sequence conservation and applied it to phytohormone REs in 45 angiosperms. Surprisingly, we found that specific RE variants are highly conserved in core hormone response genes. Experimental evidence showed that specific variants act to regulate the magnitude and spatial profile of hormonal response in Arabidopsis (Arabidopsis thaliana) and tomato (Solanum lycopersicum). Our results suggest that hormone-regulated TFs bind a spectrum of REs, each coding for a distinct transcriptional response profile. Our approach has implications for precise genome editing and for rational promoter design.

The Protein-Coding Human Genome: Annotating High-Hanging Fruits

The major transcript variants of human protein-coding genes are annotated to a certain degree of accuracy combining manual curation, transcript data, and proteomics evidence. However, there is considerable disagreement on the annotation of about 2000 genes-they can be protein-coding, noncoding, or pseudogenes-and on the annotation of most of the predicted alternative transcripts. Pure transcriptome mapping approaches seem to be limited in discriminating functional expression from noise. These limitations have partially been overcome by dedicated algorithms to detect alternative spliced micro-exons and wobble splice variants. Recently, knowledge about splice mechanism and protein structure are incorporated into an algorithm to predict neighboring homologous exons, often spliced in a mutually exclusive manner. Predicted exons are evaluated by transcript data, structural compatibility, and evolutionary conservation, revealing hundreds of novel coding exons and splice mechanism re-assignments. The emerging human pan-genome is necessitating distinctive annotations incorporating differences between individuals and between populations.

The various and shared roles of lncRNAs during development

lncRNAs, genes transcribed but not translated, longer than 200 nucleotides, are classified as a separate class of nonprotein coding genes. Since their discovery, largely from RNAseq data, a number of pioneer studies have begun to unravel its myriad functions, centered on gene expression regulation, suggesting developmental and evolutionary conservation. Since they do not code for proteins and have no open reading frames, their functional constraints likely differ from that of protein coding genes, or of genes where the majority of the nucleotide sequence is required for function, like tRNAs. This has complicated assessment of both developmental and evolutionary conservation, and the identification of homologs in different species. Here we argue that other characteristics: general synteny and particular chromosomal placement regardless of sequence, sequence micro-motifs, and secondary structure allow for "homologs" to be identified and compared, confirming developmental and evolutionary conservation of lncRNAs. We conclude exemplifying a case in point: that of the evolutionarily conserved lncRNA acal, characterized and required for embryogenesis in Drosophila.

An Evaluation of Function of Multicopy Noncoding RNAs in Mammals Using ENCODE/FANTOM Data and Comparative Genomics

Mammalian diversification has coincided with a rapid proliferation of various types of noncoding RNAs, including members of both snRNAs and snoRNAs. The significance of this expansion however remains obscure. While some ncRNA copy-number expansions have been linked to functionally tractable effects, such events may equally likely be neutral, perhaps as a result of random retrotransposition. Hindering progress in our understanding of such observations is the difficulty in establishing function for the diverse features that have been identified in our own genome. Projects such as ENCODE and FANTOM have revealed a hidden world of genomic expression patterns, as well as a host of other potential indicators of biological function. However, such projects have been criticized, particularly from practitioners in the field of molecular evolution, where many suspect these data provide limited insight into biological function. The molecular evolution community has largely taken a skeptical view, thus it is important to establish tests of function. We use a range of data, including data drawn from ENCODE and FANTOM, to examine the case for function for the recent copy number expansion in mammals of six evolutionarily ancient RNA families involved in splicing and rRNA maturation. We use several criteria to assess evidence for function: conservation of sequence and structure, genomic synteny, evidence for transposition, and evidence for species-specific expression. Applying these criteria, we find that only a minority of loci show strong evidence for function and that, for the majority, we cannot reject the null hypothesis of no function.

LncRNA DUXAP10 modulates cell proliferation in esophageal squamous cell carcinoma through epigenetically silencing p21

Esophageal Squamous Cell Carcinoma (ESCC) belongs to malignant tumor of human digestive system. It has greatly threatened human health both mentally and physically. Long non-coding RNAs have been discovered to be special molecular regulators in various cancers, including ESCC. LncRNA DUXAP10 is a newfound RNA, which is able to improve the progression of cancers ^1-3 . In this study, DUXAP10 was certified to be upregulated in ESCC tissues and cells. Besides, it was positively correlated with short survival time. Moreover, down-expression of DUXAP10 contributed to decreased cell proliferation and metastasis. Silenced DUXAP10 led to increased apoptosis rate and stagnation of cell cycle. Results of mechanism experiments suggested that DUXAP10 motivated ESCC progression through recruiting enhancer of zeste homolog 2 (EZH2) to the promoter of p21. Our findings suggested that the pseudogene-derived from lncRNA DUXAP10 drove the biological progression of ESCC. DUXAP10 was likely to be a potential therapeutic target for ESCC.

Correlated Selection on Amino Acid Deletion and Replacement in Mammalian Protein Sequences

A low ratio of nonsynonymous and synonymous substitution rates (dN/dS) at a codon is an indicator of functional constraint caused by purifying selection. Intuitively, the functional constraint would also be expected to prevent such a codon from being deleted. However, to the best of our knowledge, the correlation between the rates of deletion and substitution has never actually been estimated. Here, we use 8595 protein-coding region sequences from nine mammalian species to examine the relationship between deletion rate and dN/dS. We find significant positive correlations at the levels of both sites and genes. We compared our data against controls consisting of simulated coding sequences evolving along identical phylogenetic trees, where deletions occur independently of substitutions. A much weaker correlation was found in the corresponding simulated sequences, probably caused by alignment errors. In the real data, the correlations cannot be explained by alignment errors. Separate investigations on nonsynonymous (dN) and synonymous (dS) substitution rates indicate that the correlation is most likely due to a similarity in patterns of selection rather than in mutation rates.

Unsolvable Problems of Biology: It Is Impossible to Create Two Identical Organisms, to Defeat Cancer, or to Map Organisms onto Their Genomes

The review is devoted to unsolvable problems of biology. 1) Problems unsolvable due to stochastic mutations occurring during DNA replication that make it impossible to create two identical organisms or even two identical complex cells (Sverdlov, E. D. (2009) Biochemistry (Moscow), 74, 939-944) and to "defeat" cancer. 2) Problems unsolvable due to multiple interactions in complex systems leading to the appearance of unpredictable emergent properties that prevent establishment of unambiguous relationships between the genetic architecture and phenotypic manifestation of the genome and make impossible to predict with certainty responses of the organism, its parts, or pathological processes to external factors. 3) Problems unsolvable because of the uncertainty principle and observer effect in biology, due to which it is impossible to obtain adequate information about cells in their tissue microenvironment by isolating and analyzing individual cells. In particular, we cannot draw conclusions on the properties of stem cells in their niches based on the properties of stem cell cultures. A strategy is proposed for constructing the pattern most closely approximated to the relationship of genotypes with their phenotypes by designing networks of intermediate phenotypes (endophenotypes).

Ribosomal RNA Genes Contribute to the Formation of Pseudogenes and Junk DNA in the Human Genome

Approximately 35% of the human genome can be identified as sequence devoid of a selected-effect function, and not derived from transposable elements or repeated sequences. We provide evidence supporting a known origin for a fraction of this sequence. We show that: 1) highly degraded, but near full length, ribosomal DNA (rDNA) units, including both 45S and Intergenic Spacer (IGS), can be found at multiple sites in the human genome on chromosomes without rDNA arrays, 2) that these rDNA sequences have a propensity for being centromere proximal, and 3) that sequence at all human functional rDNA array ends is divergent from canonical rDNA to the point that it is pseudogenic. We also show that small sequence strings of rDNA (from 45S + IGS) can be found distributed throughout the genome and are identifiable as an “rDNA-like signal”, representing 0.26% of the q-arm of HSA21 and ∼2% of the total sequence of other regions tested. The size of sequence strings found in the rDNA-like signal intergrade into the size of sequence strings that make up the full-length degrading rDNA units found scattered throughout the genome. We conclude that the displaced and degrading rDNA sequences are likely of a similar origin but represent different stages in their evolution towards random sequence. Collectively, our data suggests that over vast evolutionary time, rDNA arrays contribute to the production of junk DNA. The concept that the production of rDNA pseudogenes is a by-product of concerted evolution represents a previously under-appreciated process; we demonstrate here its importance.

Efficient termination of nuclear lncRNA transcription promotes mitochondrial genome maintenance

Most DNA in the genomes of higher organisms does not code for proteins. RNA Polymerase II (Pol II) transcribes non-coding DNA into long non-coding RNAs (lncRNAs), but biological roles of lncRNA are unclear. We find that mutations in the yeast lncRNA CUT60 result in poor growth. Defective termination of CUT60 transcription causes read-through transcription across the ATP16 gene promoter. Read-through transcription localizes chromatin signatures associated with Pol II elongation to the ATP16 promoter. The act of Pol II elongation across this promoter represses functional ATP16 expression by a Transcriptional Interference (TI) mechanism. Atp16p function in the mitochondrial ATP-synthase complex promotes mitochondrial DNA stability. ATP16 repression by TI through inefficient termination of CUT60 therefore triggers mitochondrial genome loss. Our results expand the functional and mechanistic implications of non-coding DNA in eukaryotes by highlighting termination of nuclear lncRNA transcription as mechanism to stabilize an organellar genome.

An Upper Limit on the Functional Fraction of the Human Genome

For the human population to maintain a constant size from generation to generation, an increase in fertility must compensate for the reduction in the mean fitness of the population caused, among others, by deleterious mutations. The required increase in fertility due to this mutational load depends on the number of sites in the genome that are functional, the mutation rate, and the fraction of deleterious mutations among all mutations in functional regions. These dependencies and the fact that there exists a maximum tolerable replacement level fertility can be used to put an upper limit on the fraction of the human genome that can be functional. Mutational load considerations lead to the conclusion that the functional fraction within the human genome cannot exceed 15%.

Most m6A RNA Modifications in Protein-Coding Regions Are Evolutionarily Unconserved and Likely Nonfunctional

Methylation of the adenosine base at the nitrogen-6 position (m6A) is the most prevalent internal posttranscriptional modification of mRNAs in many eukaryotes. Despite the rapid progress in the transcriptome-wide mapping of m6As, identification of proteins responsible for writing, reading, and erasing m6As, and elucidation of m6A functions in splicing, RNA stability, translation, and other processes, it is unknown whether most observed m6A modifications are functional. To address this question, we respectively analyze the evolutionary conservation of yeast and human m6As in protein-coding regions. Relative to comparable unmethylated As, m6As are overall no more conserved in yeasts and only slightly more conserved in mammals. Furthermore, yeast m6As and comparable unmethylated As have no significant difference in single nucleotide polymorphism (SNP) density or SNP site frequency spectrum. The same is true in human. The methylation status of a gene, not necessarily the specific sites methylated in the gene, is subject to purifying selection for no more than ∼20% of m6A-modified genes. These observations suggest that most m6A modifications in protein-coding regions are nonfunctional and nonadaptive, probably resulting from off-target activities of m6A methyltransferases. In addition, our reanalysis invalidates the recent claim of positive selection for newly acquired m6A modifications in human evolution. Regarding the small number of evolutionarily conserved m6As, evidence suggests that a large proportion of them are likely functional; they should be prioritized in future functional characterizations of m6As. Together, these findings have important implications for understanding the biological significance of m6A and other posttranscriptional modifications.

On causal roles and selected effects: our genome is mostly junk

The idea that much of our genome is irrelevant to fitness-is not the product of positive natural selection at the organismal level-remains viable. Claims to the contrary, and specifically that the notion of "junk DNA" should be abandoned, are based on conflating meanings of the word "function". Recent estimates suggest that perhaps 90% of our DNA, though biochemically active, does not contribute to fitness in any sequence-dependent way, and possibly in no way at all. Comparisons to vertebrates with much larger and smaller genomes (the lungfish and the pufferfish) strongly align with such a conclusion, as they have done for the last half-century.

Hierarchy, determinism, and specificity in theories of development and evolution

The concepts of hierarchical organization, genetic determinism and biological specificity (for example of species, biologically relevant macromolecules, or genes) have played a crucial role in biology as a modern experimental science since its beginnings in the nineteenth century. The idea of genetic information (specificity) and genetic determination was at the basis of molecular biology that developed in the 1940s with macromolecules, viruses and prokaryotes as major objects of research often labelled "reductionist". However, the concepts have been marginalized or rejected in some of the research that in the late 1960s began to focus additionally on the molecularization of complex biological structures and functions using systems approaches. This paper challenges the view that 'molecular reductionism' has been successfully replaced by holism and a focus on the collective behaviour of cellular entities. It argues instead that there are more fertile replacements for molecular 'reductionism', in which genomics, embryology, biochemistry, and computer science intertwine and result in research that is as exact and causally predictive as earlier molecular biology.

The Evolution of Dark Matter in the Mitogenome of Seed Beetles

Animal mitogenomes are generally thought of as being economic and optimized for rapid replication and transcription. We use long-read sequencing technology to assemble the remarkable mitogenomes of four species of seed beetles. These are the largest circular mitogenomes ever assembled in insects, ranging from 24,496 to 26,613 bp in total length, and are exceptional in that some 40% consists of non-coding DNA. The size expansion is due to two very long intergenic spacers (LIGSs), rich in tandem repeats. The two LIGSs are present in all species but vary greatly in length (114-10,408 bp), show very low sequence similarity, divergent tandem repeat motifs, a very high AT content and concerted length evolution. The LIGSs have been retained for at least some 45 my but must have undergone repeated reductions and expansions, despite strong purifying selection on protein coding mtDNA genes. The LIGSs are located in two intergenic sites where a few recent studies of insects have also reported shorter LIGSs (>200 bp). These sites may represent spaces that tolerate neutral repeat array expansions or, alternatively, the LIGSs may function to allow a more economic translational machinery. Mitochondrial respiration in adult seed beetles is based almost exclusively on fatty acids, which reduces the need for building complex I of the oxidative phosphorylation pathway (NADH dehydrogenase). One possibility is thus that the LIGSs may allow depressed transcription of NAD genes. RNA sequencing showed that LIGSs are partly transcribed and transcriptional profiling suggested that all seven mtDNA NAD genes indeed show low levels of transcription and co-regulation of transcription across sexes and tissues.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

The pseudogene derived from long non-coding RNA DUXAP10 promotes colorectal cancer cell growth through epigenetically silencing of p21 and PTEN

Recently, substantial evidence has demonstrated that pseudogene derived lncRNAs are crucial regulators of cancer development and progression. DUXAP10,a pseudogene derived long non-coding RNA(lncRNA), is overexpression in colorectal cancer (CRC), but its expression pattern, biological function and underlying mechanism in CRC is still undetermined. In this study, we observed that DUXAP10 was up-regulated in CRC tissues which was positively correlated with advanced pathological stages, larger tumor sizes and lymph node metastasis. Additionally, knockdown of DUXAP10 inhibited cell proliferation, induced cell apoptosis and increase the number of G0/G1 cells significantly in the HCT116 and SW480 cell lines. Moreover, DUXAP10 silencing inhibited tumor growth in vivo. Further mechanism study showed that, by binding to histone demethylase lysine-specific demethylase 1 (LSD1), DUXAP10 promote CRC cell growth and reduced cell apoptosis through silencing the expression of p21 and phosphatase and tensin homolog (PTEN) tumor suppressor. Our findings suggested that the pseudogene-derived from lncRNA DUXAP10 promotes the biological progression of CRC and is likely to be a potential therapeutic target for CRC intervention.

Identification of introns harboring functional sequence elements through positional conservation

Many human introns carry out a function, in the sense that they are critical to maintain normal cellular activity. Their identification is fundamental to understanding cellular processes and disease. However, being noncoding elements, such functional introns are poorly predicted based on traditional approaches of sequence and structure conservation. Here, we generated a dataset of human functional introns that carry out different types of functions. We showed that functional introns share common characteristics, such as higher positional conservation along the coding sequence and reduced loss rates, regardless of their specific function. A unique property of the data is that if an intron is unknown to be functional, it still does not mean that it is indeed non-functional. We developed a probabilistic framework that explicitly accounts for this unique property, and predicts which specific human introns are functional. We show that we successfully predict function even when the algorithm is trained on introns with a different type of function. This ability has many implications in studying regulatory networks, gene regulation, the effect of mutations outside exons on human disease, and on our general understanding of intron evolution and their functional exaptation in mammals.

Estimating the prevalence of functional exonic splice regulatory information

In addition to coding information, human exons contain sequences necessary for correct splicing. These elements are known to be under purifying selection and their disruption can cause disease. However, the density of functional exonic splicing information remains profoundly uncertain. Several groups have experimentally investigated how mutations at different exonic positions affect splicing. They have found splice information to be distributed widely in exons, with one estimate putting the proportion of splicing-relevant nucleotides at >90%. These results suggest that splicing could place a major pressure on exon evolution. However, analyses of sequence conservation have concluded that the need to preserve splice regulatory signals only slightly constrains exon evolution, with a resulting decrease in the average human rate of synonymous evolution of only 1–4%. Why do these two lines of research come to such different conclusions? Among other reasons, we suggest that the methods are measuring different things: one assays the density of sites that affect splicing, the other the density of sites whose effects on splicing are visible to selection. In addition, the experimental methods typically consider short exons, thereby enriching for nucleotides close to the splice junction, such sites being enriched for splice-control elements. By contrast, in part owing to correction for nucleotide composition biases and to the assumption that constraint only operates on exon ends, the conservation-based methods can be overly conservative.

When TADs go bad: chromatin structure and nuclear organisation in human disease

Chromatin in the interphase nucleus is organised as a hierarchical series of structural domains, including self-interacting domains called topologically associating domains (TADs). This arrangement is thought to bring enhancers into closer physical proximity with their target genes, which often are located hundreds of kilobases away in linear genomic distance. TADs are demarcated by boundary regions bound by architectural proteins, such as CTCF and cohesin, although much remains to be discovered about the structure and function of these domains. Recent studies of TAD boundaries disrupted in engineered mouse models show that boundary mutations can recapitulate human developmental disorders as a result of aberrant promoter-enhancer interactions in the affected TADs. Similar boundary disruptions in certain cancers can result in oncogene overexpression, and CTCF binding sites at boundaries appear to be hyper-mutated across cancers. Further insights into chromatin organisation, in parallel with accumulating whole genome sequence data for disease cohorts, are likely to yield additional valuable insights into the roles of noncoding sequence variation in human disease.

Unveiling the Impact of the Genomic Architecture on the Evolution of Vertebrate microRNAs

Eukaryotic genomes frequently exhibit interdependency between transcriptional units, as evidenced by regions of high gene density. It is well recognized that vertebrate microRNAs (miRNAs) are usually embedded in those regions. Recent work has shown that the genomic context is of utmost importance to determine miRNA expression in time and space, thus affecting their evolutionary fates over long and short terms. Consequently, understanding the inter- and intraspecific changes on miRNA genomic architecture may bring novel insights on the basic cellular processes regulated by miRNAs, as well as phenotypic evolution and disease-related mechanisms.

Splendor and misery of adaptation, or the importance of neutral null for understanding evolution

The study of any biological features, including genomic sequences, typically revolves around the question: what is this for? However, population genetic theory, combined with the data of comparative genomics, clearly indicates that such a "pan-adaptationist" approach is a fallacy. The proper question is: how has this sequence evolved? And the proper null hypothesis posits that it is a result of neutral evolution: that is, it survives by sheer chance provided that it is not deleterious enough to be efficiently purged by purifying selection. To claim adaptation, the neutral null has to be falsified. The adaptationist fallacy can be costly, inducing biologists to relentlessly seek function where there is none.

The mutational hazard hypothesis of organelle genome evolution: 10 years on

Why is there such a large variation in size and noncoding DNA content among organelle genomes? One explanation is that this genomic variation results from differences in the rates of organelle mutation and random genetic drift, as opposed to being the direct product of natural selection. Along these lines, the mutational hazard hypothesis (MHH) holds that 'excess' DNA is a mutational liability (because it increases the potential for harmful mutations) and, thus, has a greater tendency to accumulate in an organelle system with a low mutation rate as opposed to one with a high rate of mutation. Various studies have explored this hypothesis and, more generally, the relationship between organelle genome architecture and the mode and efficiency of organelle DNA repair. Although some of these investigations are in agreement with the MHH, others have contradicted it; nevertheless, they support a central role of mutation, DNA maintenance pathways and random genetic drift in fashioning organelle chromosomes. Arguably, one of the most important contributions of the MHH is that it has sparked crucial, widespread discussions about the importance of nonadaptive processes in genome evolution.

Parlez-vous VUS?

Human genome sequencing is routine and will soon be a staple in research and clinical genetics. However, the promise of sequencing is often just that, with genome data routinely failing to reveal useful insights about disease in general or a person's health in particular. Nowhere is this chasm between promise and progress more evident than in the designation, "variant of uncertain significance" (VUS). Although it serves an important role, careful consideration of VUS reveals it to be a nebulous description of genomic information and its relationship to disease, symptomatic of our inability to make even crude quantitative assertions about the disease risks conferred by many genetic variants. In this perspective, I discuss the challenge of "variant interpretation" and the value of comparative and functional genomic information in meeting that challenge. Although already essential, genomic annotations will become even more important as our analytical focus widens beyond coding exons. Combined with more genotype and phenotype data, they will help facilitate more quantitative and insightful assessments of the contributions of genetic variants to disease.

The meaning of biological information

Biological information encoded in genomes is fundamentally different from and effectively orthogonal to Shannon entropy. The biologically relevant concept of information has to do with 'meaning', i.e. encoding various biological functions with various degree of evolutionary conservation. Apart from direct experimentation, the meaning, or biological information content, can be extracted and quantified from alignments of homologous nucleotide or amino acid sequences but generally not from a single sequence, using appropriately modified information theoretical formulae. For short, information encoded in genomes is defined vertically but not horizontally. Informally but substantially, biological information density seems to be equivalent to 'meaning' of genomic sequences that spans the entire range from sharply defined, universal meaning to effective meaninglessness. Large fractions of genomes, up to 90% in some plants, belong within the domain of fuzzy meaning. The sequences with fuzzy meaning can be recruited for various functions, with the meaning subsequently fixed, and also could perform generic functional roles that do not require sequence conservation. Biological meaning is continuously transferred between the genomes of selfish elements and hosts in the process of their coevolution. Thus, in order to adequately describe genome function and evolution, the concepts of information theory have to be adapted to incorporate the notion of meaning that is central to biology.

Multilevel Selection Theory and the Evolutionary Functions of Transposable Elements

One of several issues at play in the renewed debate over “junk DNA” is the organizational level at which genomic features might be seen as selected, and thus to exhibit function, as etiologically defined. The intuition frequently expressed by molecular geneticists that junk DNA is functional because it serves to “speed evolution” or as an “evolutionary repository” could be recast as a claim about selection between species (or clades) rather than within them, but this is not often done. Here, we review general arguments for the importance of selection at levels above that of organisms in evolution, and develop them further for a common genomic feature: the carriage of transposable elements (TEs). In many species, not least our own, TEs comprise a large fraction of all nuclear DNA, and whether they individually or collectively contribute to fitness—or are instead junk— is a subject of ongoing contestation. Even if TEs generally owe their origin to selfish selection at the lowest level (that of genomes), their prevalence in extant organisms and the prevalence of extant organisms bearing them must also respond to selection within species (on organismal fitness) and between species (on rates of speciation and extinction). At an even higher level, the persistence of clades may be affected (positively or negatively) by TE carriage. If indeed TEs speed evolution, it is at these higher levels of selection that such a function might best be attributed to them as a class.

Satellite DNA in Plants: More than Just Rubbish

For decades, satellite DNAs have been the hidden part of genomes. Initially considered as junk DNA, there is currently an increasing appreciation of the functional significance of satellite DNA repeats and of their sequences. Satellite DNA families accumulate in the heterochromatin in different parts of the eukaryotic chromosomes, mainly in pericentromeric and subtelomeric regions, but they also span the functional centromere. Tandem repeat sequences may spread from subtelomeric to interstitial loci, leading to the formation of chromosome-specific loci or to the accumulation in equilocal sites in different chromosomes. They also appear as the main components of the heterochromatin in the sex-specific region of sex chromosomes. Satellite DNA, required for chromosome organization, also plays a role in pairing and segregation. Some satellite repeats are transcribed and can participate in the formation and maintenance of heterochromatin structure and in the modulation of gene expression. In addition to the identification of the different satellite DNA families, their characteristics and location, we are interested in determining their impact on the genomes, by identifying the mechanisms leading to their appearance and amplification as well as in understanding how they change over time, the factors affecting these changes, and the influence exerted by the evolutionary history of the organisms. On the other hand, satellite DNA sequences are rapidly evolving sequences that may cause reproductive barriers between organisms and promote speciation. The accumulation of experimental data collected in recent years and the emergence of new approaches based on next-generation sequencing and high-throughput genome analysis are opening new perspectives that are changing our understanding of satellite DNA. This review examines recent data to provide a timely update on the overall information gathered about this part of the genome, focusing on the advances in the knowledge of its origin, its evolution, and its potential functional roles.

The differential view of genotype-phenotype relationships

An integrative view of diversity and singularity in the living world requires a better understanding of the intricate link between genotypes and phenotypes. Here we re-emphasize the old standpoint that the genotype-phenotype (GP) relationship is best viewed as a connection between two differences, one at the genetic level and one at the phenotypic level. As of today, predominant thinking in biology research is that multiple genes interact with multiple environmental variables (such as abiotic factors, culture, or symbionts) to produce the phenotype. Often, the problem of linking genotypes and phenotypes is framed in terms of genotype and phenotype maps, and such graphical representations implicitly bring us away from the differential view of GP relationships. Here we show that the differential view of GP relationships is a useful explanatory framework in the context of pervasive pleiotropy, epistasis, and environmental effects. In such cases, it is relevant to view GP relationships as differences embedded into differences. Thinking in terms of differences clarifies the comparison between environmental and genetic effects on phenotypes and helps to further understand the connection between genotypes and phenotypes.