Are non-protein coding RNAs junk or treasure?: An attempt to explain and reconcile opposing viewpoints of whether the human genome is mostly transcribed into non-functional or functional RNAs

The human genome project's lasting legacies are the emerging insights into human physiology and disease, and the ascendance of biology as the dominant science of the 21st century. Sequencing revealed that >90% of the human genome is not coding for proteins, as originally thought, but rather is overwhelmingly transcribed into non-protein coding, or non-coding, RNAs (ncRNAs). This discovery initially led to the hypothesis that most genomic DNA is "junk", a term still championed by some geneticists and evolutionary biologists. In contrast, molecular biologists and biochemists studying the vast number of transcripts produced from most of this genome "junk" often surmise that these ncRNAs have biological significance. What gives? This essay contrasts the two opposing, extant viewpoints, aiming to explain their bases, which arise from distinct reference frames of the underlying scientific disciplines. Finally, it aims to reconcile these divergent mindsets in hopes of stimulating synergy between scientific fields.

Mutational profiling of mitochondrial DNA reveals an epithelial ovarian cancer-specific evolutionary pattern contributing to high oxidative metabolism

BACKGROUND

Epithelial ovarian cancer (EOC) heavily relies on oxidative phosphorylation (OXPHOS) and exhibits distinct mitochondrial metabolic reprogramming. Up to now, the evolutionary pattern of somatic mitochondrial DNA (mtDNA) mutations in EOC tissues and their potential roles in metabolic remodelling have not been systematically elucidated.

METHODS

Based on a large somatic mtDNA mutation dataset from private and public EOC cohorts (239 and 118 patients, respectively), we most comprehensively characterised the EOC-specific evolutionary pattern of mtDNA mutations and investigated its biological implication.

RESULTS

Mutational profiling revealed that the mitochondrial genome of EOC tissues was highly unstable compared with non-cancerous ovary tissues. Furthermore, our data indicated the delayed heteroplasmy accumulation of mtDNA control region (mtCTR) mutations and near-complete absence of mtCTR non-hypervariable segment (non-HVS) mutations in EOC tissues, which is consistent with stringent negative selection against mtCTR mutation. Additionally, we observed a bidirectional and region-specific evolutionary pattern of mtDNA coding region mutations, manifested as significant negative selection against mutations in complex V (ATP6/ATP8) and tRNA loop regions, and potential positive selection on mutations in complex III (MT-CYB). Meanwhile, EOC tissues showed higher mitochondrial biogenesis compared with non-cancerous ovary tissues. Further analysis revealed the significant association between mtDNA mutations and both mitochondrial biogenesis and overall survival of EOC patients.

CONCLUSIONS

Our study presents a comprehensive delineation of EOC-specific evolutionary patterns of mtDNA mutations that aligned well with the specific mitochondrial metabolic remodelling, conferring novel insights into the functional roles of mtDNA mutations in EOC tumourigenesis and progression.

Characterization of the complete chloroplast genome of Gypsophila huashanensis Y. W. Tsui & D. Q. Lu, an endemic herb species in China

Gypsophila huashanensis Y. W. Tsui & D. Q. Lu (Caryophyllaceae) is an endemic herb species to the Qinling Mountains in China. In this study, we characterized its whole plastid genome using the Illumina sequencing platform. The complete plastid genome of G. huashanensis is 152,457 bp in length, including a large single-copy DNA region of 83,476 bp, a small single-copy DNA region of 17,345 bp, and a pair of inverted repeat DNA sequences of 25,818 bp. The genome contains 130 genes comprising 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Evolutionary analysis showed that the non-coding regions of Caryophyllaceae exhibit a higher level of divergence than the exon regions. Gene site selection analysis suggested that 11 coding protein genes (accD, atpF, ndhA, ndhB, petB, petD, rpoCl, rpoC2, rps16, ycfl, and ycf2) have some sites under protein sequence evolution. Phylogenetic analysis showed that G. huashanensis is most closely related to the congeneric species G. oldhamiana. These results are very useful for studying phylogenetic evolution and species divergence in the family Caryophyllaceae.

Mutational signature of mtDNA confers mechanistic insight into oxidative metabolism remodeling in colorectal cancer

Rationale: Mitochondrial dysfunction caused by mitochondrial DNA (mtDNA) mutations and subsequent metabolic defects are closely involved in tumorigenesis and progression in a cancer-type specific manner. To date, the mutational pattern of mtDNA somatic mutations in colorectal cancer (CRC) tissues and its clinical implication are still not completely clear. Methods: In the present study, we generated a large mtDNA somatic mutation dataset from three CRC cohorts (432, 1,015, and 845 patients, respectively) and then most comprehensively characterized the CRC-specific evolutionary pattern and its clinical implication. Results: Our results showed that the mtDNA control region (mtCTR) with a high mutation density exhibited a distinct mutation spectrum characterizing a high enrichment of L-strand C > T mutations, which was contrary to the H-strand C > T mutational bias observed in the mtDNA coding region (mtCDR) (P < 0.001). Further analysis clearly confirmed the relaxed evolutionary selection of mtCTR mutations, which was mainly characterized by the similar distribution of hypervariable region (HVS) and non-HVS mutation density. Moreover, significant negative selection was identified in mutations of mtDNA complex V (ATP6/ATP8) and tRNA loop regions. Although our data showed that oxidative metabolism was commonly increased in CRC cells, mtDNA somatic mutations in CRC tissues were not closely associated with mitochondrial biogenesis, oxidative metabolism, and clinical progression, suggesting a cancer-type specific relationship between mtDNA mutations and mitochondrial metabolic functions in CRC cells. Conclusion: Our study identified the CRC-specific evolutionary mode of mtDNA mutations, which is possibly matched to specific mitochondrial metabolic remodeling and confers new mechanic insight into CRC tumorigenesis.

The Evolution of G-quadruplex Structure in mRNA Untranslated Region

The RNA G-quadruplex (rG4) is a kind of non-canonical high-order secondary structure with important biological functions and is enriched in untranslated regions (UTRs) of protein-coding genes. However, how rG4 structures evolve is largely unknown. Here, we systematically investigated the evolution of RNA sequences around UTR rG4 structures in 5 eukaryotic organisms. We found universal selection on UTR sequences, which facilitated rG4 formation in all the organisms that we analyzed. While G-rich sequences were preferred in the rG4 structural region, C-rich sequences were selectively not preferred. The selective pressure acting on rG4 structures in the UTRs of genes with higher G content was significantly smaller. Furthermore, we found that rG4 structures experienced smaller evolutionary selection near the translation initiation region in the 5' UTR, near the polyadenylation signals in the 3' UTR, and in regions flanking the miRNA targets in the 3' UTR. These results suggest universal selection for rG4 formation in the UTRs of eukaryotic genomes and the selection may be related to the biological functions of rG4s.

Phylogenetic relationships and molecular evolution of woody forest tree family Aceraceae based on plastid phylogenomics and nuclear gene variations

The forest tree family Aceraceae is widespread in the northern hemisphere and it has ecological and economic importance. However, the phylogenetic relationships and classifications within the family are still controversial due to transitional intraspecific morphological characteristics and introgression hybridization among species. In this study, we determined the evolutionary relationships and molecular evolution of Aceraceae based on plastid phylogenomics and two nuclear gene variations. Phylogenetic analysis based on the plastid genomes suggested that Aceraceae species can be divided into two larger sub-clades corresponding to the two genera Acer and Dipteronia. Conjoint analysis of the plastid and nuclear gene sequences supported the classification with two genera in the family. Molecular dating showed that the two genera diverged 60.2 million years ago, which is generally consistently with previously reported results. Divergence hotspots and positively selected genes identified in the plastid genomes could be useful genetic resources in Aceraceae.

Coevolutionary and Phylogenetic Analysis of Mimiviral Replication Machinery Suggest the Cellular Origin of Mimiviruses

Mimivirus is one of the most complex and largest viruses known. The origin and evolution of Mimivirus and other giant viruses have been a subject of intense study in the last two decades. The two prevailing hypotheses on the origin of Mimivirus and other viruses are the reduction hypothesis, which posits that viruses emerged from modern unicellular organisms; whereas the virus-first hypothesis proposes viruses as relics of precellular forms of life. In this study, to gain insights into the origin of Mimivirus, we have carried out extensive phylogenetic, correlation, and multidimensional scaling analyses of the putative proteins involved in the replication of its 1.2-Mb large genome. Correlation analysis and multidimensional scaling methods were validated using bacteriophage, bacteria, archaea, and eukaryotic replication proteins before applying to Mimivirus. We show that a large fraction of mimiviral replication proteins, including polymerase B, clamp, and clamp loaders are of eukaryotic origin and are coevolving. Although phylogenetic analysis places some components along the lineages of phage and bacteria, we show that all the replication-related genes have been homogenized and are under purifying selection. Collectively our analysis supports the idea that Mimivirus originated from a complex cellular ancestor. We hypothesize that Mimivirus has largely retained complex replication machinery reminiscent of its progenitor while losing most of the other genes related to processes such as metabolism and translation.

Hypothalamic transcriptome of tame and aggressive silver foxes (Vulpes vulpes) identifies gene expression differences shared across brain regions

The underlying neurological events accompanying dog domestication remain elusive. To reconstruct the domestication process in an experimental setting, silver foxes (Vulpes vulpes) have been deliberately bred for tame vs aggressive behaviors for more than 50 generations at the Institute for Cytology and Genetics in Novosibirsk, Russia. The hypothalamus is an essential part of the hypothalamic-pituitary-adrenal axis and regulates the fight-or-flight response, and thus, we hypothesized that selective breeding for tameness/aggressiveness has shaped the hypothalamic transcriptomic profile. RNA-seq analysis identified 70 differentially expressed genes (DEGs). Seven of these genes, DKKL1, FBLN7, NPL, PRIMPOL, PTGRN, SHCBP1L and SKIV2L, showed the same direction expression differences in the hypothalamus, basal forebrain and prefrontal cortex. The genes differentially expressed across the three tissues are involved in cell division, differentiation, adhesion and carbohydrate processing, suggesting an association of these processes with selective breeding. Additionally, 159 transcripts from the hypothalamus demonstrated differences in the abundance of alternative spliced forms between the tame and aggressive foxes. Weighted gene coexpression network analyses also suggested that gene modules in hypothalamus were significantly associated with tame vs aggressive behavior. Pathways associated with these modules include signal transduction, interleukin signaling, cytokine-cytokine receptor interaction and peptide ligand-binding receptors (eg, G-protein coupled receptor [GPCR] ligand binding). Current studies show the selection for tameness vs aggressiveness in foxes is associated with unique hypothalamic gene profiles partly shared with other brain regions and highlight DEGs involved in biological processes such as development, differentiation and immunological responses. The role of these processes in fox and dog domestication remains to be determined.

Plasmodium Genomics and Genetics: New Insights into Malaria Pathogenesis, Drug Resistance, Epidemiology, and Evolution

Protozoan Plasmodium parasites are the causative agents of malaria, a deadly disease that continues to afflict hundreds of millions of people every year. Infections with malaria parasites can be asymptomatic, with mild or severe symptoms, or fatal, depending on many factors such as parasite virulence and host immune status. Malaria can be treated with various drugs, with artemisinin-based combination therapies (ACTs) being the first-line choice. Recent advances in genetics and genomics of malaria parasites have contributed greatly to our understanding of parasite population dynamics, transmission, drug responses, and pathogenesis. However, knowledge gaps in parasite biology and host-parasite interactions still remain. Parasites resistant to multiple antimalarial drugs have emerged, while advanced clinical trials have shown partial efficacy for one available vaccine. Here we discuss genetic and genomic studies of Plasmodium biology, host-parasite interactions, population structures, mosquito infectivity, antigenic variation, and targets for treatment and immunization. Knowledge from these studies will advance our understanding of malaria pathogenesis, epidemiology, and evolution and will support work to discover and develop new medicines and vaccines.

Genetic diversity of CHC22 clathrin impacts its function in glucose metabolism

CHC22 clathrin plays a key role in intracellular membrane traffic of the insulin-responsive glucose transporter GLUT4 in humans. We performed population genetic and phylogenetic analyses of the CHC22-encoding CLTCL1 gene, revealing independent gene loss in at least two vertebrate lineages, after arising from gene duplication. All vertebrates retained the paralogous CLTC gene encoding CHC17 clathrin, which mediates endocytosis. For vertebrates retaining CLTCL1, strong evidence for purifying selection supports CHC22 functionality. All human populations maintained two high frequency CLTCL1 allelic variants, encoding either methionine or valine at position 1316. Functional studies indicated that CHC22-V1316, which is more frequent in farming populations than in hunter-gatherers, has different cellular dynamics than M1316-CHC22 and is less effective at controlling GLUT4 membrane traffic, altering its insulin-regulated response. These analyses suggest that ancestral human dietary change influenced selection of allotypes that affect CHC22's role in metabolism and have potential to differentially influence the human insulin response.

Behind the potential evolution towards prion resistant species

Historically, the observation of naturally occurring cases of prion disease led to the classification of different susceptibility grades and to the designation of prion resistant species. However, the development of highly efficient in vitro prion propagation systems and the generation of ad hoc transgenic models allowed determining that leporidae and equidae families have been erroneously considered resistant to prion infection. On the contrary, similar approaches revealed an unexpected high level of resistance of the canidae family. In PLoS Pathogens [ 1 ], we describe experiments directed toward elucidating which are the determinants of the alleged prion resistance of this family. Studies based on the sequence of the canine prion protein coupled with structural in silico analysis identified a key residue probably implicated in this resistance. Cell and brain-based PMCA highlighted that the presence of aspartic or glutamic acid at codon 163 of the canid PrP, strongly inhibits prion replication in vitro. Transgenic animals carrying this substitution in mouse PrP were resistant to prion infection after intracerebral challenge with different mouse prion strains. The confirmation of the importance of this substitution and its exclusivity in this family, suggests it could have been evolutionarily favored, due to their diet based on carrion and small ruminants.

Punishment diminishes the benefits of network reciprocity in social dilemma experiments

Network reciprocity has been widely advertised in theoretical studies as one of the basic cooperation-promoting mechanisms, but experimental evidence favoring this type of reciprocity was published only recently. When organized in an unchanging network of social contacts, human subjects cooperate provided the following strict condition is satisfied: The benefit of cooperation must outweigh the total cost of cooperating with all neighbors. In an attempt to relax this condition, we perform social dilemma experiments wherein network reciprocity is aided with another theoretically hypothesized cooperation-promoting mechanism-costly punishment. The results reveal how networks promote and stabilize cooperation. This stabilizing effect is stronger in a smaller-size neighborhood, as expected from theory and experiments. Contrary to expectations, punishment diminishes the benefits of network reciprocity by lowering assortment, payoff per round, and award for cooperative behavior. This diminishing effect is stronger in a larger-size neighborhood. An immediate implication is that the psychological effects of enduring punishment override the rational response anticipated in quantitative models of cooperation in networks.

Orphan CpG islands as alternative promoters

CpG islands (CGIs) are associated with ∼60% of mammalian promoters. Most unmethylated CGIs exhibit transcriptional activity, which has led to their co-option as promoters by retrogenes. CGIs may also serve as alternative promoters for downstream genes with methylated promoters, with implications on aberrant activation of oncogenes in cancer phenotypes.

Evidence of a Direct Evolutionary Selection for Strong Folding and Mutational Robustness Within HIV Coding Regions

A large number of studies demonstrated the importance of different HIV RNA structural elements at all stages of the viral life cycle. Nevertheless, the significance of many of these structures is unknown, and plausibly new regions containing RNA structure-mediated regulatory signals remain to be identified. An important characteristic of genomic regions carrying functionally significant secondary structures is their mutational robustness, that is, the extent to which a sequence remains constant in spite of despite mutations in terms of its underlying secondary structure. Structural robustness to mutations is expected to be important in the case of functional RNA structures in viruses with high mutation rate; it may prevent fitness loss due to disruption of possibly functional conformations, pointing to the specific significance of the corresponding genomic region. In the current work, we perform a genome-wide computational analysis to detect signals of a direct evolutionary selection for strong folding and RNA structure-based mutational robustness within HIV coding sequences. We provide evidence that specific regions of HIV structural genes undergo an evolutionary selection for strong folding; in addition, we demonstrate that HIV Rev responsive element seems to undergo a direct evolutionary selection for increased secondary structure robustness to point mutations. We believe that our analysis may enable a better understanding of viral evolutionary dynamics at the RNA structural level and may benefit to practical efforts of engineering antiviral vaccines and novel therapeutic approaches.

Evolutionarily conserved pressure for the existence of distinct G2/M cell cycle arrest and A3H inactivation functions in HIV-1 Vif

HIV-1 Vif assembles the Cul5-EloB/C E3 ubiquitin ligase to induce proteasomal degradation of the cellular antiviral APOBEC3 proteins. Detailed structural studies have confirmed critical functional domains in Vif that we have previously identified as important for the interaction of EloB/C, Cul5, and CBFβ. However, the mechanism by which Vif recognizes substrates remains poorly understood. Specific regions of Vif have been identified as being responsible for binding and depleting APOBEC3G and APOBEC3F. Interestingly, we have now identified distinct yet overlapping domains that are required for HIV-1 Vif-mediated G2/M-phase cell cycle arrest and APOBEC3H degradation, but not for the inactivation of APOBEC3G or APOBEC3F. Surprisingly, Vif molecules from primary HIV-1 variants that caused G2/M arrest were unable to inactivate APOBEC3H; on the other hand, HIV-1 Vif variants that could inactivate APOBEC3H were unable to induce G2/M arrest. All of these Vif variants still maintained the ability to inactivate APOBEC3G/F. Thus, primary HIV-1 variants have evolved to possess distinct functional activities that allow them to suppress APOBEC3H or cause G2 cell cycle arrest, using mutually exclusive interface domains. APOBEC3H depletion and G2 arrest are apparently evolutionary selected features that cannot co-exist on a single Vif molecule. The existence and persistence of both types of HIV-1 Vif variant suggests the importance of APOBEC3H suppression and cell cycle regulation for HIV-1's survival in vivo.

Human SNPs resulting in premature stop codons and protein truncation

Single nucleotide polymorphisms (SNPs) constitute the most common type of genetic variation in humans. SNPs introducing premature termination codons (PTCs), herein called X-SNPs, can alter the stability and function of transcripts and proteins and thus are considered to be biologically important. Initial studies suggested a strong selection against such variations/mutations. In this study, we undertook a genome-wide systematic screening to identify human X-SNPs using the dbSNP database. Our results demonstrated the presence of 28 X-SNPs from 28 genes with known minor allele frequencies. Eight X-SNPs (28.6 per cent) were predicted to cause transcript degradation by nonsense-mediated mRNA decay. Seventeen X-SNPs (60.7 per cent) resulted in moderate to severe truncation at the C-terminus of the proteins (deletion of >50 per cent of the amino acids). The majority of the X-SNPs (78.6 per cent) represent commonly occurring SNPs, by contrast with the rarely occurring disease-causing PTC mutations. Interestingly, X-SNPs displayed a non-uniform distribution across human populations: eight X-SNPs were reported to be prevalent across three different human populations, whereas six X-SNPs were found exclusively in one or two population(s). In conclusion, we have systematically investigated human SNPs introducing PTCs with respect to their possible biological consequences, distributions across different human populations and evolutionary aspects. We believe that the SNPs reported here are likely to affect gene/protein function, although their biological and evolutionary roles need to be further investigated.