Origin of functional de novo genes in humans from "hopeful monsters".

For a long time, it was believed that new genes arise only from modifications of preexisting genes, but the discovery of de novo protein-coding genes that originated from noncoding DNA regions demonstrates the existence of a "motherless" origination process for new genes. However, the features, distributions, expression profiles, and origin modes of these genes in humans seem to support the notion that their origin is not a purely "motherless" process; rather, these genes arise preferentially from genomic regions encoding preexisting precursors with gene-like features. In such a case, the gene loci are typically not brand new. In this short review, we will summarize the definition and features of human de novo genes and clarify their process of origination from ancestral non-coding genomic regions. In addition, we define the favored precursors, or "hopeful monsters," for the origin of de novo genes and present a discussion of the functional significance of these young genes in brain development and tumorigenesis in humans. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Adaptive functions of structural variants in human brain development.

Quantifying the structural variants (SVs) in nonhuman primates could provide a niche to clarify the genetic backgrounds underlying human-specific traits, but such resource is largely lacking. Here, we report an accurate SV map in a population of 562 rhesus macaques, verified by in-house benchmarks of eight macaque genomes with long-read sequencing and another one with genome assembly. This map indicates stronger selective constrains on inversions at regulatory regions, suggesting a strategy for prioritizing them with the most important functions. Accordingly, we identified 75 human-specific inversions and prioritized them. The top-ranked inversions have substantially shaped the human transcriptome, through their dual effects of reconfiguring the ancestral genomic architecture and introducing regional mutation hotspots at the inverted regions. As a proof of concept, we linked APCDD1, located on one of these inversions and down-regulated specifically in humans, to neuronal maturation and cognitive ability. We thus highlight inversions in shaping the human uniqueness in brain development.

6mA-Sniper: Quantifying 6mA sites in eukaryotes at single-nucleotide resolution.

While N6-methyldeoxyadenine (6mA) modification is a fundamental regulation in prokaryotes, its prevalence and functions in eukaryotes are controversial. Here, we report 6mA-Sniper to quantify 6mA sites in eukaryotes at single-nucleotide resolution, and delineate a 6mA profile in Caenorhabditis elegans with 2034 sites. Twenty-six of 39 events with Mnl I restriction endonuclease sites were verified, demonstrating the feasibility of this method. The levels of 6mA sites pinpointed by 6mA-Sniper are generally increased after Pseudomonas aeruginosa infection, but decreased in strains with the removal of METL-9, the dominant 6mA methyltransferase. The enrichment of these sites on specific motif of [GC]GAG, the selective constrains on them, and their coordinated changes with METL-9 levels thus support an active shaping of the 6mA profile by methyltransferase. Moreover, for regions marked by 6mA sites that emerged after infection, an enrichment of up-regulated genes was detected, possibly mediated through a mutual exclusive cross-talk between 6mA and H3K27me3 modification. We thus highlight 6mA regulation as a previously neglected regulator in eukaryotes.

A novel N6-Deoxyadenine methyltransferase METL-9 modulates C. elegans immunity via dichotomous mechanisms.

N6-Methyldeoxyadenine (6mA) has been rediscovered as a DNA modification with potential biological function in metazoans. However, the physiological function and regulatory mechanisms regarding the establishment, maintenance and removal of 6mA in eukaryotes are still poorly understood. Here we show that genomic 6mA levels change in response to pathogenic infection in Caenorhabditis elegans (C. elegans). We further identify METL-9 as the methyltransferase that catalyzes DNA 6mA modifications upon pathogen infection. Deficiency of METL-9 impairs the induction of innate immune response genes and renders the animals more susceptible to pathogen infection. Interestingly, METL-9 functions through both 6mA-dependent and -independent mechanisms to transcriptionally regulate innate immunity. Our findings reveal that 6mA is a functional DNA modification in immunomodulation in C. elegans.

STellaris: a web server for accurate spatial mapping of single cells based on spatial transcriptomics data.

Single-cell RNA sequencing (scRNA-seq) provides insights into gene expression heterogeneities in diverse cell types underlying homeostasis, development and pathological states. However, the loss of spatial information hinders its applications in deciphering spatially related features, such as cell-cell interactions in a spatial context. Here, we present STellaris (https://spatial.rhesusbase.com), a web server aimed to rapidly assign spatial information to scRNA-seq data based on their transcriptomic similarity with public spatial transcriptomics (ST) data. STellaris is founded on 101 manually curated ST datasets comprising 823 sections across different organs, developmental stages and pathological states from humans and mice. STellaris accepts raw count matrix and cell type annotation of scRNA-seq data as the input, and maps single cells to spatial locations in the tissue architecture of properly matched ST section. Spatially resolved information for intercellular communications, such as spatial distance and ligand-receptor interactions (LRIs), are further characterized between annotated cell types. Moreover, we also expanded the application of STellaris in spatial annotation of multiple regulatory levels with single-cell multiomics data, using the transcriptome as a bridge. STellaris was applied to several case studies to showcase its utility of adding value to the ever-growing scRNA-seq data from a spatial perspective.

Host-mediated RNA editing in viruses.

Viruses rely on hosts for life and reproduction, cause a variety of symptoms from common cold to AIDS to COVID-19 and provoke public health threats claiming millions of lives around the globe. RNA editing, as a crucial co-/post-transcriptional modification inducing nucleotide alterations on both endogenous and exogenous RNA sequences, exerts significant influences on virus replication, protein synthesis, infectivity and toxicity. Hitherto, a number of host-mediated RNA editing sites have been identified in diverse viruses, yet lacking a full picture of RNA editing-associated mechanisms and effects in different classes of viruses. Here we synthesize the current knowledge of host-mediated RNA editing in a variety of viruses by considering two enzyme families, viz., ADARs and APOBECs, thereby presenting a landscape of diverse editing mechanisms and effects between viruses and hosts. In the ongoing pandemic, our study promises to provide potentially valuable insights for better understanding host-mediated RNA editing on ever-reported and newly-emerging viruses.

A Human-Specific De Novo Gene Promotes Cortical Expansion and Folding.

Newly originated de novo genes have been linked to the formation and function of the human brain. However, how a specific gene originates from ancestral noncoding DNAs and becomes involved in the preexisting network for functional outcomes remains elusive. Here, a human-specific de novo gene, SP0535, is identified that is preferentially expressed in the ventricular zone of the human fetal brain and plays an important role in cortical development and function. In human embryonic stem cell-derived cortical organoids, knockout of SP0535 compromises their growth and neurogenesis. In SP0535 transgenic (TG) mice, expression of SP0535 induces fetal cortex expansion and sulci and gyri-like structure formation. The progenitors and neurons in the SP0535 TG mouse cortex tend to proliferate and differentiate in ways that are unique to humans. SP0535 TG adult mice also exhibit improved cognitive ability and working memory. Mechanistically, SP0535 interacts with the membrane protein Na+ /K+ ATPase subunit alpha-1 (ATP1A1) and releases Src from the ATP1A1-Src complex, allowing increased level of Src phosphorylation that promotes cell proliferation. Thus, SP0535 is the first proven human-specific de novo gene that promotes cortical expansion and folding, and can function through incorporating into an existing conserved molecular network.

De novo genes with an lncRNA origin encode unique human brain developmental functionality.

Human de novo genes can originate from neutral long non-coding RNA (lncRNA) loci and are evolutionarily significant in general, yet how and why this all-or-nothing transition to functionality happens remains unclear. Here, in 74 human/hominoid-specific de novo genes, we identified distinctive U1 elements and RNA splice-related sequences accounting for RNA nuclear export, differentiating mRNAs from lncRNAs, and driving the origin of de novo genes from lncRNA loci. The polymorphic sites facilitating the lncRNA-mRNA conversion through regulating nuclear export are selectively constrained, maintaining a boundary that differentiates mRNAs from lncRNAs. The functional new genes actively passing through it thus showed a mode of pre-adaptive origin, in that they acquire functions along with the achievement of their coding potential. As a proof of concept, we verified the regulations of splicing and U1 recognition on the nuclear export efficiency of one of these genes, the ENSG00000205704, in human neural progenitor cells. Notably, knock-out or over-expression of this gene in human embryonic stem cells accelerates or delays the neuronal maturation of cortical organoids, respectively. The transgenic mice with ectopically expressed ENSG00000205704 showed enlarged brains with cortical expansion. We thus demonstrate the key roles of nuclear export in de novo gene origin. These newly originated genes should reflect the novel uniqueness of human brain development.

CHDbase: A Comprehensive Knowledgebase for Congenital Heart Disease-related Genes and Clinical Manifestations.

Congenital heart disease (CHD) is one of themost common causes of major birth defects, with a prevalence of 1%. Although an increasing number of studies have reported the etiology of CHD, the findings scattered throughout the literature are difficult to retrieve and utilize in research and clinical practice. We therefore developed CHDbase, an evidence-based knowledgebase of CHD-related genes and clinical manifestations manually curated from 1114 publications, linking 1124susceptibility genes and 3591 variations to more than 300 CHD types and related syndromes. Metadata such as the information of each publication and the selected population and samples, the strategy of studies, and the major findings of studies were integrated with each item of the research record. We also integrated functional annotations through parsing ∼ 50 databases/tools to facilitate the interpretation of these genes and variations in disease pathogenicity. We further prioritized the significance of these CHD-related genes with a gene interaction network approach and extracted a core CHD sub-network with 163 genes. The clear genetic landscape of CHD enables the phenotype classification based on the shared genetic origin. Overall, CHDbase provides a comprehensive and freely available resource to study CHD susceptibilities, supporting a wide range of users in the scientific and medical communities. CHDbase is accessible at http://chddb.fwgenetics.org.

ICAM-1-related noncoding RNA accelerates atherosclerosis by amplifying NF-κB signaling.

Long noncoding RNAs (lncRNAs) are critical regulators of inflammation with great potential as new therapeutic targets. However, the role of lncRNAs in early atherosclerosis remains poorly characterized. This study aimed to identify the key lncRNA players in activated endothelial cells (ECs). The lncRNAs in response to pro-inflammatory factors in ECs were screened through RNA sequencing. ICAM-1-related non-coding RNA (ICR) was identified as the most potential candidate for early atherosclerosis. ICR is essential for intercellular adhesion molecule-1 (ICAM1) expression, EC adhesion and migration. In a high fat diet-induced atherosclerosis model in mice, ICR is upregulated in the development of atherosclerosis. After intravenous injection of adenovirus carrying shRNA for mouse ICR, the atherosclerotic plaque area was markedly reduced with the declined expression of ICR and ICAM1. Mechanistically, ICR stabilized the mRNA of ICAM1 in quiescent ECs; while under inflammatory stress, ICR upregulated ICAM1 in a nuclear factor kappa B (NF-κB) dependent manner. RNA-seq analysis showed pro-inflammatory targets of NF-κB were regulated by ICR. Furthermore, the chromatin immunoprecipitation assays showed that p65 binds to ICR promoter and facilitates its transcription. Interestingly, ICR, in turn, promotes p65 accumulation and activity, forming a positive feedback loop to amplify NF-κB signaling. Preventing the degradation of p65 using proteasome inhibitors rescued the expression of NF-κB targets suppressed by ICR. Taken together, ICR acts as an accelerator to amplify NF-κB signaling in activated ECs and suppressing ICR is a promising early intervention for atherosclerosis through ICR/p65 loop blockade.

A comprehensive RNA editome reveals that edited Azin1 partners with DDX1 to enable hematopoietic stem cell differentiation.

Adenosine-to-inosine RNA editing and the catalyzing enzyme adenosine deaminase are both essential for hematopoietic development and differentiation. However, the RNA editome during hematopoiesis and the underlying mechanisms are poorly defined. Here, we sorted 12 murine adult hematopoietic cell populations at different stages and identified 30 796 editing sites through RNA sequencing. The dynamic landscape of the RNA editome comprises stage- and group-specific and stable editing patterns, but undergoes significant changes during lineage commitment. Notably, we found that antizyme inhibitor 1 (Azin1) was highly edited in hematopoietic stem and progenitor cells (HSPCs). Azin1 editing results in an amino acid change to induce Azin1 protein (AZI) translocation to the nucleus, enhanced AZI binding affinity for DEAD box polypeptide 1 to alter the chromatin distribution of the latter, and altered expression of multiple hematopoietic regulators that ultimately promote HSPC differentiation. Our findings have delineated an essential role for Azin1 RNA editing in hematopoietic cells, and our data set is a valuable resource for studying RNA editing on a more general basis.

Polyadenylation-related isoform switching in human evolution revealed by full-length transcript structure.

Rhesus macaque is a unique nonhuman primate model for human evolutionary and translational study, but the error-prone gene models critically limit its applications. Here, we de novo defined full-length macaque gene models based on single molecule, long-read transcriptome sequencing in four macaque tissues (frontal cortex, cerebellum, heart and testis). Overall, 8 588 227 poly(A)-bearing complementary DNA reads with a mean length of 14 106 nt were generated to compile the backbone of macaque transcripts, with the fine-scale structures further refined by RNA sequencing and cap analysis gene expression sequencing data. In total, 51 605 macaque gene models were accurately defined, covering 89.7% of macaque or 75.7% of human orthologous genes. Based on the full-length gene models, we performed a human-macaque comparative analysis on polyadenylation (PA) regulation. Using macaque and mouse as outgroup species, we identified 79 distal PA events newly originated in humans and found that the strengthening of the distal PA sites, rather than the weakening of the proximal sites, predominantly contributes to the origination of these human-specific isoforms. Notably, these isoforms are selectively constrained in general and contribute to the temporospatially specific reduction of gene expression, through the tinkering of previously existed mechanisms of nuclear retention and microRNA (miRNA) regulation. Overall, the protocol and resource highlight the application of bioinformatics in integrating multilayer genomics data to provide an intact reference for model animal studies, and the isoform switching detected may constitute a hitherto underestimated regulatory layer in shaping the human-specific transcriptome and phenotypic changes.

Functional and Adaptive Significance of Promoter Mutations That Affect Divergent Myocardial Expressions of TRIM72 in Primates.

Cis-regulatory elements play important roles in tissue-specific gene expression and in the evolution of various phenotypes, and mutations in promoters and enhancers may be responsible for adaptations of species to environments. TRIM72 is a highly conserved protein that is involved in energy metabolism. Its expression in the heart varies considerably in primates, with high levels of expression in Old World monkeys and near absence in hominids. Here, we combine phylogenetic hypothesis testing and experimentation to demonstrate that mutations in promoter are responsible for the differences among primate species in the heart-specific expression of TRIM72. Maximum likelihood estimates of lineage-specific substitution rates under local-clock models show that relative to the evolutionary rate of introns, the rate of promoter was accelerated by 78% in the common ancestor of Old World monkeys, suggesting a role for positive selection in the evolution of the TRIM72 promoter, possibly driven by selective pressure due to changes in cardiac physiology after species divergence. We demonstrate that mutations in the TRIM72 promoter account for the differential myocardial TRIM72 expression of the human and the rhesus macaque. Furthermore, changes in TRIM72 expression alter the expression of genes involved in oxidative phosphorylation, which in turn affects mitochondrial respiration and cardiac energy capacity. On a broader timescale, phylogenetic regression analyses of data from 29 mammalian species show that mammals with high cardiac expression of TRIM72 have high heart rate, suggesting that the expression changes of TRIM72 may be related to differences in the heart physiology of those species.

HLH-11 modulates lipid metabolism in response to nutrient availability.

The ability of organisms to sense nutrient availability and tailor their metabolic states to withstand nutrient deficiency is critical for survival. To identify previously unknown regulators that couple nutrient deficiency to body fat utilization, we performed a cherry-picked RNAi screen in C. elegans and found that the transcription factor HLH-11 regulates lipid metabolism in response to food availability. In well-fed worms, HLH-11 suppresses transcription of lipid catabolism genes. Upon fasting, the HLH-11 protein level is reduced through lysosome- and proteasome-mediated degradation, thus alleviating the inhibitory effect of HLH-11, activating the transcription of lipid catabolism genes, and utilizing fat. Additionally, lipid profiling revealed that reduction in the HLH-11 protein level remodels the lipid landscape in C. elegans. Moreover, TFAP4, the mammalian homolog of HLH-11, plays an evolutionarily conserved role in regulating lipid metabolism in response to starvation. Thus, TFAP4 may represent a potential therapeutic target for lipid storage disorders.

Histone deacetylase HDA-1 modulates mitochondrial stress response and longevity.

The ability to detect, respond and adapt to mitochondrial stress ensures the development and survival of organisms. Caenorhabditis elegans responds to mitochondrial stress by activating the mitochondrial unfolded protein response (UPRmt) to buffer the mitochondrial folding environment, rewire the metabolic state, and promote innate immunity and lifespan extension. Here we show that HDA-1, the C. elegans ortholog of mammalian histone deacetylase (HDAC) is required for mitochondrial stress-mediated activation of UPRmt. HDA-1 interacts and coordinates with the genome organizer DVE-1 to induce the transcription of a broad spectrum of UPRmt, innate immune response and metabolic reprogramming genes. In rhesus monkey and human tissues, HDAC1/2 transcript levels correlate with the expression of UPRmt genes. Knocking down or pharmacological inhibition of HDAC1/2 disrupts the activation of the UPRmt and the mitochondrial network in mammalian cells. Our results underscore an evolutionarily conserved mechanism of HDAC1/2 in modulating mitochondrial homeostasis and regulating longevity.

Prohibitin 2 deficiency impairs cardiac fatty acid oxidation and causes heart failure.

Fatty acids are the most major substrate source for adult cardiac energy generation. Prohibitin 2 (PHB2), a highly conserved protein located in mitochondrial inner membrane, plays key roles in cellular energy metabolic homeostasis. However, its functions in regulating cardiac fatty acid metabolism have remained largely unknown. Our study demonstrates that cardiac-specific knockout of Phb2 leads to accumulation of lipid droplets and causes heart failure. Mechanistically, ablation of PHB2 impairs cardiac fatty acid oxidation (FAO) through downregulating carnitine palmitoyltransferase1b (CPT1b), a rate-limiting enzyme of cardiac mitochondrial FAO. Moreover, overexpression of CPT1b alleviates impaired FAO in PHB2-deficient cardiomyocytes. Thus, our study provides direct evidence for the link between PHB2 and cardiac fatty acid metabolism. Our study points out that PHB2 is a potential FAO regulator in cardiac mitochondrial inner membrane, as well as the connection between PHB2 and CPT1b and their relationships to cardiac pathology especially to cardiac fatty acid metabolic disorder.

CaMKII-δ9 promotes cardiomyopathy through disrupting UBE2T-dependent DNA repair.

Ca2+/calmodulin-dependent kinase II (CaMKII) is a multifunctional serine/threonine kinase family, and its δ isoform is predominant in the heart. Excessive CaMKII activation plays a pivotal role in the pathogenesis of severe heart conditions, including myocardial infarction, cardiomyopathy and heart failure. However, the identity of CaMKII splice variants and the mechanism(s) underlying CaMKII-mediated cardiac pathology remain elusive. Here, we show that CaMKII-δ9, the most abundant CaMKII-δ splice variant in human heart, potently promotes cardiomyocyte death, cardiomyopathy and heart failure by disrupting cardiomyocyte genome stability. Mechanistically, CaMKII-δ9, but not the previously well-studied CaMKII-δ2 and CaMKII-δ3, targets the ubiquitin-conjugating enzyme E2T (UBE2T) for phosphorylation and degradation, disrupting UBE2T-dependent DNA repair and leading to the accumulation of DNA damage and genome instability. These findings not only reveal a crucial role of CaMKII in the regulation of DNA repair, but also mark the CaMKII-δ9-UBE2T-DNA damage pathway as an important therapeutic target for cardiomyopathy and heart failure.

Replication-Independent Histone Turnover Underlines the Epigenetic Homeostasis in Adult Heart.

RATIONALE: Replication-independent histone turnover has been linked to cis-regulatory chromatin domains in cultured cell lines, but its molecular underpinnings and functional relevance in adult mammalian tissues remain yet to be defined. OBJECTIVE: We investigated regulatory functions of replication-independent histone turnover in chromatin states of postmitotic cardiomyocytes from adult mouse heart. METHODS AND RESULTS: We used H2B-GFP (histone 2B-green fluorescent protein) fusion protein pulse-and-chase approaches to measure histone turnover rate in mouse cardiomyocytes. Surprisingly, we found that the short histone half-life (≈2 weeks) contrasted dramatically with the long lifetime of cardiomyocytes, and rapid histone turnover regions corresponded to cis-regulatory domains of heart genes. Interestingly, recruitment of chromatin modifiers, including Polycomb EED (embryonic ectoderm development), was positively correlated with histone turnover rate at enhancers. Mechanistically, through directly interacting with and engaging the BAF (BRG1 [Brahma-related gene-1]-associated factor) complex for nucleosome exchange for stereotyped histone modifications from the free histone pool, EED augmented histone turnover to restrain enhancer overactivation. CONCLUSIONS: We propose a model in which replication-independent histone turnover reinforces robustness of local chromatin states for adult tissue homeostasis.

ADAR1 promotes robust hypoxia signaling via distinct regulation of multiple HIF-1α-inhibiting factors.

Adenosine deaminase acting on RNA (ADAR)-catalyzed adenosine-to-inosine RNA editing is potentially dysregulated in neoplastic progression. However, how this transcriptome recoding process is functionally correlated with tumorigenesis remains largely elusive. Our analyses of RNA editome datasets identify hypoxia-related genes as A-to-I editing targets. In particular, two negative regulators of HIF-1A-the natural antisense transcript HIF1A-AS2 and the ubiquitin ligase scaffold LIMD1-are directly but differentially modulated by ADAR1. We show that HIF1A-AS2 antagonizes the expression of HIF-1A in the immediate-early phase of hypoxic challenge, likely through a convergent transcription competition in cis ADAR1 in turn suppresses transcriptional progression of the antisense gene. In contrast, ADAR1 affects LIMD1 expression post-transcriptionally, by interfering with the cytoplasmic translocation of LIMD1 mRNA and thus protein translation. This multi-tier regulation coordinated by ADAR1 promotes robust and timely accumulation of HIF-1α upon oxygen depletion and reinforces target gene induction and downstream angiogenesis. Our results pinpoint ADAR1-HIF-1α axis as a hitherto unrecognized key regulator in hypoxia.

Evolutionarily significant A-to-I RNA editing events originated through G-to-A mutations in primates.

BACKGROUND: Recent studies have revealed thousands of A-to-I RNA editing events in primates, but the origination and general functions of these events are not well addressed. RESULTS: Here, we perform a comparative editome study in human and rhesus macaque and uncover a substantial proportion of macaque A-to-I editing sites that are genomically polymorphic in some animals or encoded as non-editable nucleotides in human. The occurrence of these recent gain and loss of RNA editing through DNA point mutation is significantly more prevalent than that expected for the nearby regions. Ancestral state analyses further demonstrate that an increase in recent gain of editing events contribute to the over-representation, with G-to-A mutation site as a favorable location for the origination of robust A-to-I editing events. Population genetics analyses of the focal editing sites further reveal that a portion of these young editing events are evolutionarily significant, indicating general functional relevance for at least a fraction of these sites. CONCLUSIONS: Overall, we report a list of A-to-I editing events that recently originated through G-to-A mutations in primates, representing a valuable resource to investigate the features and evolutionary significance of A-to-I editing events at the population and species levels. The unique subset of primate editome also illuminates the general functions of RNA editing by connecting it to particular gene regulatory processes, based on the characterized outcome of a gene regulatory level in different individuals or primate species with or without these editing events.