Ultraconserved elements in the human genome

Lycibarbarspermidine L from Goji Berry Promotes Intestinal Restoration in an Antibiotic-Induced Rat Model through Targeting uc.141A

Oral antibiotics damage gut epithelium while treating infections. Goji berry derives from Lycium barbarum L., which is employed as a dietary and tonic herb with the effect of promoting intestinal health. The spermidine compound lycibarbarspermidine L (LBS L) in it can promote intestinal epithelial growth. The purpose of the research is to examine the reparative effects of LBS L against antibiotic-stimulated gut epithelial damage, as well as the processes associated with T-UCR. Intestinal epithelial cells were stimulated with antibiotics to induce injury. The impact of LBS L and T-UCR-related mechanisms was evaluated through cell assays. Ultimately, a rat model of intestinal injury was developed by administering antibiotics. We found that LBS L reversed the intestinal damage caused by antibiotics through modulating uc.141A/miR195-3p. A reduction in intestinal damage was observed in rat models. In conclusion, LBS L ameliorates antibiotic-induced intestinal epithelial cell damage in rats by modulating uc.141A.

Massively parallel reporter assays and mouse transgenic assays provide correlated and complementary information about neuronal enhancer activity

High-throughput massively parallel reporter assays (MPRAs) and phenotype-rich in vivo transgenic mouse assays are two potentially complementary ways to study the impact of noncoding variants associated with psychiatric diseases. Here, we investigate the utility of combining these assays. Specifically, we carry out an MPRA in induced human neurons on over 50,000 sequences derived from fetal neuronal ATAC-seq datasets and enhancers validated in mouse assays. We also test the impact of over 20,000 variants, including synthetic mutations and 167 common variants associated with psychiatric disorders. We find a strong and specific correlation between MPRA and mouse neuronal enhancer activity. Four out of five tested variants with significant MPRA effects affected neuronal enhancer activity in mouse embryos. Mouse assays also reveal pleiotropic variant effects that could not be observed in MPRA. Our work provides a catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

Deciphering enhancers of hearing loss genes for efficient and targeted gene therapy of hereditary deafness

Hereditary hearing loss accounts for about 60% of congenital deafness. Although adeno-associated virus (AAV)-mediated gene therapy shows substantial potential for treating genetic hearing impairments, there remain significant concerns regarding the specificity and safety of AAV vectors. The sophisticated nature of the cochlea further complicates the challenge of precisely targeting gene delivery. Here, we introduced an AAV-reporter-based in vivo transcriptional enhancer reconstruction (ARBITER) workflow, enabling efficient and reliable dissection of enhancers. With ARBITER, we successfully demonstrated that the conserved non-coding elements (CNEs) within the gene locus collaboratively regulate the expression of Slc26a5, which was further validated using knockout mouse models. We also assessed the potential of identified enhancers to treat hereditary hearing loss by conducting gene therapy in Slc26a5 mutant mice. Based on the original Slc26a5 enhancer with limited efficiency, we engineered a highly efficient and outer hair cell (OHC)-specific enhancer, B8, which successfully restored hearing of Slc26a5 knockout mice.

A deletion at the X-linked ARHGAP36 gene locus is associated with the orange coloration of tortoiseshell and calico cats

The X-linked orange (O) locus in domestic cats controls an unknown molecular mechanism that causes the suppression of black-brownish pigmentation in favor of orange coloration. The alternating black-brownish and orange patches seen in tortoiseshell and calico cats are considered classic examples of the phenotypic expression of random X chromosome inactivation (XCI) occurring in female mammals. However, the O gene in the cat genome has not been identified, and the genetic variation responsible for the orange coloration remains unknown. We report here that a 5.1-kilobase (kb) deletion within an intron of the X-linked ARHGAP36 gene, encoding a Rho GTPase-activating protein, is closely and exclusively associated with orange coloration. The deleted region contains a highly conserved putative regulatory element, whose removal is presumed to alter ARHGAP36 expression. Notably, ARHGAP36 expression in cat skin tissues is linked to the suppression of many melanogenesis genes, potentially shifting pigment synthesis from eumelanin to pheomelanin. Furthermore, we find evidence that the gene undergoes XCI in female human and mouse cells and XCI-dependent CpG island methylation consistent with random XCI in female domestic cats. The 5.1-kb deletion seems widespread in domestic cats with orange coat coloration, suggesting a single origin of this coat color phenotype.

Molecular and genetic characterization of sex-linked orange coat color in the domestic cat

The Sex-linked orange mutation in domestic cats causes variegated patches of reddish/yellow hair and is a defining signature of random X inactivation in female tortoiseshell and calico cats. Unlike the situation for most coat color genes, there is no apparent homolog for Sex-linked orange in other mammals. We show that Sex-linked orange is caused by a 5-kb deletion that leads to ectopic and melanocyte-specific expression of the Rho GTPase Activating Protein 36 (Arhgap36) gene. Single-cell RNA sequencing (RNA-seq) studies from fetal cat skin reveal that red/yellow hair color is caused by reduced expression of melanogenic genes that are normally activated by the melanocortin 1 receptor (Mc1r)-cyclic adenosine monophosphate (cAMP)-protein kinase A (PKA) pathway, but Mc1r and its ability to stimulate cAMP accumulation is intact. Instead, we show that expression of Arhgap36 in melanocytes leads to reduced levels of the PKA catalytic subunit (PKAC); thus, Sex-linked orange is genetically and biochemically downstream of Mc1r. Our findings resolve a longstanding comparative genetic puzzle, provide in vivo evidence for the ability of Arhgap36 to inhibit PKA, and reveal a molecular explanation for a charismatic color pattern with a rich genetic history.

MAFin: motif detection in multiple alignment files

MOTIVATION

Whole Genome and Proteome Alignments, represented by the multiple alignment file format, have become a standard approach in comparative genomics and proteomics. These often require identifying conserved motifs, which is crucial for understanding functional and evolutionary relationships. However, current approaches lack a direct method for motif detection within MAF files. We present MAFin, a novel tool that enables efficient motif detection and conservation analysis in MAF files to address this gap, streamlining genomic and proteomic research.

RESULTS

We developed MAFin, the first motif detection tool for Multiple Alignment Format files. MAFin enables the multithreaded search of conserved motifs using three approaches: (i) using user-specified k-mers to search the sequences. (ii) with regular expressions, in which case one or more patterns are searched, and (iii) with predefined Position Weight Matrices. Once the motif has been found, MAFin detects the motif instances and calculates the conservation across the aligned sequences. MAFin also calculates a conservation percentage, which provides information about the conservation levels of each motif across the aligned sequences, based on the number of matches relative to the length of the motif. A set of statistics enables the interpretation of each motif's conservation level, and the detected motifs are exported in JSON and CSV files for downstream analyses.

AVAILABILITY AND IMPLEMENTATION

MAFin is offered as a Python package under the GPL license as a multi-platform application and is available at: https://github.com/Georgakopoulos-Soares-lab/MAFin.

Machine learning-based prediction reveals kinase MAP4K4 regulates neutrophil differentiation through phosphorylating apoptosis-related proteins

Neutrophils, an essential innate immune cell type with a short lifespan, rely on continuous replenishment from bone marrow (BM) precursors. Although it is established that neutrophils are derived from the granulocyte-macrophage progenitor (GMP), the molecular regulators involved in the differentiation process remain poorly understood. Here we developed a random forest-based machine-learning pipeline, NeuRGI (Neutrophil Regulatory Gene Identifier), which utilized Positive-Unlabeled Learning (PU-learning) and neural network-based in silico gene knockout to identify neutrophil regulators. We interrogated features including gene expression dynamics, physiological characteristics, pathological relatedness, and gene conservation for the model training. Our identified pipeline leads to identifying Mitogen-Activated Protein Kinase-4 (MAP4K4) as a novel neutrophil differentiation regulator. The loss of MAP4K4 in hematopoietic stem cells and progenitors in mice induced neutropenia and impeded the differentiation of neutrophils in the bone marrow. By modulating the phosphorylation level of proteins involved in cell apoptosis, such as STAT5A, MAP4K4 delicately regulates cell apoptosis during the process of neutrophil differentiation. Our work presents a novel regulatory mechanism in neutrophil differentiation and provides a robust prediction model that can be applied to other cellular differentiation processes.

The translatome of glioblastoma

Glioblastoma (GB), the most common and aggressive brain tumor, demonstrates intrinsic resistance to current therapies, resulting in poor clinical outcomes. Cancer progression can be partially attributed to the deregulation of protein translation mechanisms that drive cancer cell growth. In this study, we present the translatome landscape of GB as a valuable data resource. Eight patient-derived GB sphere cultures (GSCs) were analyzed using ribosome profiling and messenger RNA (mRNA) sequencing. We investigated inter-cell-line differences through differential expression analysis at both the translatome and transcriptome levels. Translational changes post-radiotherapy were assessed at 30 and 60 min. The translation of non-coding RNAs (ncRNAs) was validated using in-house and public mass spectrometry (MS) data, whereas RNA expression was confirmed by quantitative PCR (qPCR). Our findings demonstrate that ribosome sequencing provides more detailed information than MS or transcriptional analyses. Transcriptional similarities among GSCs correlate with translational similarities, aligning with previously defined subtypes such as proneural and mesenchymal. Additionally, we identified a broad spectrum of open reading frame types in both coding and non-coding mRNA regions, including long non-coding RNAs (lncRNAs) and pseudogenes undergoing active translation. Translation of ncRNAs into peptides was independently confirmed by in-house data and external MS data. We also observed that translational regulation of histones (downregulated) and splicing factors (upregulated) occurs in response to radiotherapy. These data offer new insights into genome-wide protein synthesis, identifying translationally regulated genes and alternative translation initiation sites in GB under normal and radiotherapeutic conditions, providing a rich resource for GB research. Further functional validation of differentially expressed genes after radiotherapy is needed. Understanding translational control in GB can reveal mechanistic insights and identify currently unknown biomarkers, ultimately enhancing the diagnosis and treatment of this aggressive brain cancer.

The mutational landscape and functional effects of noncoding ultraconserved elements in human cancers

The mutational landscape of phylogenetically ultraconserved elements (UCEs), especially those in noncoding DNAs (ncUCEs), and their functional relevance in cancers remain poorly characterized. Here, we perform a systematic analysis of whole-genome and in-house targeted UCE sequencing datasets from more than 3000 patients with cancer of 13,736 UCEs and demonstrate that ncUCE somatic alterations are common. Using a multiplexed CRISPR knockout screen in colorectal cancer cells, we show that the loss of several altered ncUCEs significantly affects cell proliferation. In-depth functional studies in vitro and in vivo further reveal that specific ncUCEs can be enhancers of tumor suppressors (such as ARID1B) and silencers of oncogenic proteins (such as RPS13). Moreover, several miRNAs located in ncUCEs are recurrently mutated. Mutations in miR-142 locus can affect the Drosha-mediated processing of precursor miRNAs, resulting in the down-regulation of the mature transcript. These results provide systematic evidence that specific ncUCEs play diverse regulatory roles in cancer.

Mapping the regulatory effects of common and rare non-coding variants across cellular and developmental contexts in the brain and heart

Whole genome sequencing has identified over a billion non-coding variants in humans, while GWAS has revealed the non-coding genome as a significant contributor to disease. However, prioritizing causal common and rare non-coding variants in human disease, and understanding how selective pressures have shaped the non-coding genome, remains a significant challenge. Here, we predicted the effects of 15 million variants with deep learning models trained on single-cell ATAC-seq across 132 cellular contexts in adult and fetal brain and heart, producing nearly two billion context-specific predictions. Using these predictions, we distinguish candidate causal variants underlying human traits and diseases and their context-specific effects. While common variant effects are more cell-type-specific, rare variants exert more cell-type-shared regulatory effects, with selective pressures particularly targeting variants affecting fetal brain neurons. To prioritize de novo mutations with extreme regulatory effects, we developed FLARE, a context-specific functional genomic model of constraint. FLARE outperformed other methods in prioritizing case mutations from autism-affected families near syndromic autism-associated genes; for example, identifying mutation outliers near CNTNAP2 that would be missed by alternative approaches. Overall, our findings demonstrate the potential of integrating single-cell maps with population genetics and deep learning-based variant effect prediction to elucidate mechanisms of development and disease-ultimately, supporting the notion that genetic contributions to neurodevelopmental disorders are predominantly rare.

A humanized NOVA1 splicing factor alters mouse vocal communications

NOVA1, a neuronal RNA-binding protein expressed in the central nervous system, is essential for survival in mice and normal development in humans. A single amino acid change (I197V) in NOVA1's second RNA binding domain is unique to modern humans. To study its physiological effects, we generated mice carrying the human-specific I197V variant (Nova1hu/hu) and analyzed the molecular and behavioral consequences. While the I197V substitution had minimal impact on NOVA1's RNA binding capacity, it led to specific effects on alternative splicing, and CLIP revealed multiple binding peaks in mouse brain transcripts involved in vocalization. These molecular findings were associated with behavioral differences in vocalization patterns in Nova1hu/hu mice as pups and adults. Our findings suggest that this human-specific NOVA1 substitution may have been part of an ancient evolutionary selective sweep in a common ancestral population of Homo sapiens, possibly contributing to the development of spoken language through differential RNA regulation during brain development.

MixtureFinder: Estimating DNA Mixture Models for Phylogenetic Analyses

In phylogenetic studies, both partitioned models and mixture models are used to account for heterogeneity in molecular evolution among the sites of DNA sequence alignments. Partitioned models require the user to specify the grouping of sites into subsets, and then assume that each subset of sites can be modeled by a single common process. Mixture models do not require users to prespecify subsets of sites, and instead calculate the likelihood of every site under every model, while co-estimating the model weights and parameters. While much research has gone into the optimization of partitioned models by merging user-specified subsets, there has been less attention paid to the optimization of mixture models for DNA sequence alignments. In this study, we first ask whether a key assumption of partitioned models-that each user-specified subset can be modeled by a single common process-is supported by the data. Having shown that this is not the case, we then design, implement, test, and apply an algorithm, MixtureFinder, to select the optimum number of classes for a mixture model of Q-matrices for the standard models of DNA sequence evolution. We show this algorithm performs well on simulated and empirical datasets and suggest that it may be useful for future empirical studies. MixtureFinder is available in IQ-TREE2, and a tutorial for using MixtureFinder can be found here: http://www.iqtree.org/doc/Complex-Models#mixture-models.

The lawc gene emerged de novo from conserved genomic elements and acquired a broad expression pattern in Drosophila

It has recently become evident that the de novo emergence of genes is widespread and documented for a variety of organisms. De novo genes frequently emerge in proximity to existing genes, forming gene overlaps. Here, we present an analysis of the evolutionary history of a putative de novo gene, lawc, which overlaps with the conserved Trf2 gene, which encodes a general transcription factor in Drosophila melanogaster. We demonstrate that lawc emerged approximately 68 million years ago in the 5'-untranslated region (UTR) of Trf2 and displays an extensive spatiotemporal expression pattern. One of the most remarkable features of the lawc evolutionary history is that its emergence was facilitated by the engagement of Drosophilidae-specific short, highly conserved regions located in Trf2 introns. This represents a unique example of putative de novo gene birth involving conserved DNA regions localized in introns of conserved genes. The observed lawc expression pattern may be due to the overlap of lawc with the 5'-UTR of Trf2. This study not only enriches our understanding of gene evolution but also highlights the complex interplay between genetic conservation and innovation.

A comprehensive analysis of Transcribed Ultra Conserved Regions uncovers important regulatory functions of novel non-coding transcripts in gliomas

Transcribed Ultra-Conserved Regions (TUCRs) represent a severely understudied class of putative non-coding RNAs (ncRNAs) that are 100% conserved across multiple species. We performed the first-ever analysis of TUCRs in glioblastoma (GBM) and low-grade gliomas (LGG). We leveraged large human datasets to identify the genomic locations, chromatin accessibility, transcription, differential expression, correlation with survival, and predicted functions of all 481 TUCRs, and identified TUCRs that are relevant to glioma biology. Of these, we investigated the expression, function, and mechanism of action of the most highly upregulated intergenic TUCR, uc.110, identifying it as a new tumor enhancer. Uc.110 was highly overexpressed in GBM and LGG, where it promoted malignancy and tumor growth. Uc.110 activated the WNT pathway by upregulating the expression of membrane frizzled-related protein (MFRP), by sponging the tumor suppressor microRNA miR-544. This pioneering study shows important roles for TUCRs in gliomas and provides an extensive database and novel methods for future TUCR research.

New Challenges in Bladder Cancer Diagnosis: How Biosensing Tools Can Lead to Population Screening Opportunities

Bladder cancer is one of the most common cancers worldwide. Despite its high incidence, cystoscopy remains the currently used diagnostic gold standard, although it is invasive, expensive and has low sensitivity. As a result, the cancer diagnosis is mostly late, as it occurs following the presence of hematuria in urine, and population screening is not allowed. It would therefore be desirable to be able to act promptly in the early stage of the disease with the aid of biosensing. The use of devices/tools based on genetic assessments would be of great help in this field. However, the genetic differences between populations do not allow accurate analysis in the context of population screening. Current research is directed towards the discovery of universal biomarkers present in urine with the aim of providing an approach based on a non-invasive, easy-to-perform, rapid, and accurate test that can be widely used in clinical practice for the early diagnosis and follow-up of bladder cancer. An efficient biosensing device may have a disruptive impact in terms of patient health and disease management, contributing to a decrease in mortality rate, as well as easing the social and economic burden on the national healthcare system. Considering the advantage of accessing population screening for early diagnosis of cancer, the main challenges and future perspectives are critically discussed to address the research towards the selection of suitable biomarkers for the development of a very sensitive biosensor for bladder cancer.

Comparative Genomics Provides Insights into Adaptive Evolution and Demographics of Bats

Bats possess a range of distinctive characteristics, including flight, echolocation, impressive longevity, and the ability to harbor various zoonotic pathogens. Additionally, they account for the second-highest species diversity among mammalian orders, yet their phylogenetic relationships and demographic history remain underexplored. Here, we generated de novo assembled genomes for 17 bat species and 2 of their mammalian relatives (the Amur hedgehog and Chinese mole shrew), with 12 genomes reaching chromosome-level assembly. Comparative genomics and ChIP-seq assays identified newly gained genomic regions in bats potentially linked to the regulation of gene activity and expression. Notably, some antiviral infection-related gene under positive selection exhibited the activity of suppressing cancer, evidencing the linkage between virus tolerance and cancer resistance in bats. By integrating published bat genome assemblies, phylogenetic reconstruction established the proximity of noctilionoid bats to vesper bats. Interestingly, we found 2 distinct patterns of ancient population dynamics in bats and population changes since the last glacial maximum does not reflect species phylogenetic relationships. These findings enriched our understanding of adaptive mechanisms and demographic history of bats.

Posttranscriptional Regulation of Intestinal Mucosal Growth and Adaptation by Noncoding RNAs in Critical Surgical Disorders

The human intestinal epithelium has an impressive ability to respond to insults and its homeostasis is maintained by well-regulated mechanisms under various pathophysiological conditions. Nonetheless, acute injury and inhibited regeneration of the intestinal epithelium occur commonly in critically ill surgical patients, leading to the translocation of luminal toxic substances and bacteria to the bloodstream. Effective therapies for the preservation of intestinal epithelial integrity and for the prevention of mucosal hemorrhage and gut barrier dysfunction are limited, primarily because of a poor understanding of the mechanisms underlying mucosal disruption. Noncoding RNAs (ncRNAs), which include microRNAs (miRNAs), long ncRNAs (lncRNAs), circular RNAs (circRNAs), and small vault RNAs (vtRNAs), modulate a wide array of biological functions and have been identified as orchestrators of intestinal epithelial homeostasis. Here, we feature the roles of many important ncRNAs in controlling intestinal mucosal growth, barrier function, and repair after injury-particularly in the context of postoperative recovery from bowel surgery. We review recent literature surrounding the relationships between lncRNAs, microRNAs, and RNA-binding proteins and how their interactions impact cell survival, proliferation, migration, and cell-to-cell interactions in the intestinal epithelium. With advancing knowledge of ncRNA biology and growing recognition of the importance of ncRNAs in maintaining the intestinal epithelial integrity, ncRNAs provide novel therapeutic targets for treatments to preserve the gut epithelium in individuals suffering from critical surgical disorders.

Phylogenomic insights into the worldwide evolutionary relationships of the stingless bees (Apidae, Meliponini)

Stingless bees (tribe Meliponini) are remarkable for their characteristically large social colonies, their capacity to produce honey and other useful products, and their morphological and behavioral diversity. They have a disjunct pan-tropical distribution, primarily occurring in warm and humid environments in the Neotropical, Afrotropical, and Indo-Australasian regions. Even though phylogenetic hypotheses have been proposed for Meliponini based on morphology and molecular data, many questions are still unsolved regarding the evolutionary relationships and systematics of the tribe. In this contribution, we present a large phylogenomic dataset comprising over 2500 ultra-conserved element (UCE) loci sequenced for 153 species of Meliponini, representing all known genera of stingless bees. The genera Camargoia, Paratrigonoides, Plectoplebeia, Cleptotrigona, Ebaiotrigona, Papuatrigona, Pariotrigona, Platytrigona, and Sahulotrigona were included in molecular phylogenetic analyses for the first time. Concatenated and species-tree analyses were performed using different partitioning strategies and summary methods. We performed gene-genealogy interrogation (GGI) on several recalcitrant nodes to resolve discordances among recovered tree topologies. Results were mostly consistent among analyses, recovering three main lineages of Meliponini congruent with the biogeographic domains to which they are associated. Within major clades, discordances were found in relation to previous works. The genus Frieseomelitta was recovered as paraphyletic in relation to Trichotrigona, and the genus Lepidotrigona was revealed to be composed of two independent lineages. Even though concatenated and weighted ASTRAL analyses were mostly effective in recovering the relationships favored by GGI, they retrieved different results in relation to the phylogenetic placements of Oxytrigona and Cephalotrigona. The most favored hypothesis in GGI analyses was not found in any other analyses, being more congruent with morphological evidence and highlighting the relevance of exploring the support given to alternative hypotheses through topological tests. Recent advances in our capacity to generate molecular sequences from old specimens using modern sequencing methods allowed for unparalleled sampling across genera, yielding a backbone for the phylogenetic relationships of stingless bees, which will further investigations into their systematics and evolution.

Historical biogeography and the evolution of habitat preference in the North American camel spider family, Eremobatidae (Arachnida:Solifugae)

Abiotic variables can influence species distributions, often restricting taxa to an acquired climatic signature or conversely, related species are conserved in the same ecological space over millions of years. An investigation into how abiotic change has shaped geographic distributions of taxa may be key to understanding diversification of lineages, and in the absence of reliable morphological characteristics, such information may support taxonomic units at multiple scales. Here, we examine the historical biogeography and patterns of habitat preference within the North American solifuge family, Eremobatidae. A previous study demonstrated that a major taxonomic revision of Eremobatidae is warranted, however recent studies demonstrate high levels of morphological convergence within the group, thus a re-classification of generic boundaries using additional information must be prioritized before we can formally begin solid revisionary efforts. In this study, we aimed to reconstruct a well-resolved phylogenetic hypothesis of Eremobatidae by filtering UCE loci based on informativeness, by mitigating the effect of cogenic UCE on phylogenetic estimation, and by supplementing our curated UCE loci with mitochondrial information. Using our preferred topology, in conjunction with published estimated divergence dates for Eremobatidae, we inferred a time-calibrated phylogenetic hypothesis to inform the historical biogeography and patterns of habitat preference. The two major habitat types that were observed for Eremobatidae were warm deserts for early diverging taxa and a subsequent evolution to cold deserts and Mediterranean California ecoregions for later diverging taxa. Eremobatid niche space, determined by temperature and precipitation, has been conserved for at least 25 million years in North America, supporting a warm desert origin, and thus supporting high species richness in the Sonoran and Mexican Plateau. Overall, our study provides support for new generic level designations within Eremobatidae.

The role of polypeptides encoded by ncRNAs in cancer

It was previously thought that ncRNA could not encode polypeptides, but recent reports have challenged this notion. As research into ncRNA progresses, it is increasingly clear that it serves roles beyond traditional mechanisms, playing significant regulatory roles in various diseases, notably cancer, which is responsible for 70% of human deaths. Numerous studies have highlighted the diverse regulatory mechanisms of ncRNA that are pivotal in cancer initiation and progression. The role of ncRNA-encoded polypeptides in cancer regulation has gained prominence. This article explores the newly identified regulatory functions of these polypeptides in three types of ncRNA-lncRNA, pri-miRNA, and circRNA. These polypeptides can interact with proteins, influence signaling pathways, enhance miRNA stability, and regulate cancer progression, malignancy, resistance, and other clinical challenges. Furthermore, we discuss the evolutionary significance of these polypeptides in the transition from RNA to protein, examining their emergence and conservation throughout evolution.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. Here, we examine the positioning and functional significance of human genes, grouped by their lineage restriction level, within the 3D organization of the genome. We reveal that genes of different lineage restriction levels have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young lineage restricted genes to ancient genes present in most species. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Molecular and genetic characterization of sex-linked orange coat color in the domestic cat

The Sex-linked orange mutation in domestic cats causes variegated patches of reddish/yellow hair and is a defining signature of random X-inactivation in female tortoiseshell and calico cats. Unlike the situation for most coat color genes, there is no apparent homolog for Sex-linked orange in other mammals. We show that the Sex-linked orange is caused by a 5 kb deletion that leads to ectopic and melanocyte-specific expression of the Rho GTPase Activating Protein 36 ( Arhgap36 ) gene. Single cell RNA-seq studies from fetal cat skin reveal that red/yellow hair color is caused by reduced expression of melanogenic genes that are normally activated by the Melanocortin 1 receptor (Mc1r)-cyclic adenosine monophosphate (cAMP)-protein kinase A (PKA) pathway, but the Mc1r gene and its ability to stimulate cAMP accumulation is intact. Instead, we show that increased expression of Arhgap36 in melanocytes leads to reduced levels of the PKA catalytic subunit (PKA C ); thus, Sex-linked orange is genetically and biochemically downstream of Mc1r . Our findings solve a comparative genomic conundrum, provide in vivo evidence for the ability of Arhgap36 to inhibit PKA, and reveal a molecular explanation for a charismatic color pattern with a rich genetic history.

Non-coding RNA notations, regulations and interactive resources

An increasing number of non-coding RNAs (ncRNAs) are found to have roles in gene expression and cellular regulations. However, there are still a large number of ncRNAs whose functions remain to be studied. Despite decades of research, the field continues to evolve, with each newly identified ncRNA undergoing processes such as biogenesis, identification, and functional annotation. Bioinformatics methodologies, alongside traditional biochemical experimental methods, have played an important role in advancing ncRNA research across various stages. Presently, over 50 types of ncRNAs have been characterized, each exhibiting diverse functions. However, there remains a need for standardization and integration of these ncRNAs within a unified framework. In response to this gap, this review traces the historical trajectory of ncRNA research and proposes a unified notation system. Additionally, we comprehensively elucidate the ncRNA interactome, detailing its associations with DNAs, RNAs, proteins, complexes, and chromatin. A web portal named ncRNA Hub ( https://bis.zju.edu.cn/nchub/ ) is also constructed to provide detailed notations of ncRNAs and share a collection of bioinformatics resources. This review aims to provide a broader perspective and standardized paradigm for advancing ncRNA research.

Convergent Degenerated Regulatory Elements Associated with Limb Loss in Limbless Amphibians and Reptiles

Limbs are a defining characteristic of tetrapods, yet numerous taxa, primarily among amphibians and reptiles, have independently lost limbs as an adaptation to new ecological niches. To elucidate the genetic factors contributing to this convergent limb loss, we present a 12 Gb chromosome-level assembly of the Banna caecilian (Ichthyophis bannanicus), a limbless amphibian. Our comparative analysis, which includes the reconstruction of amphibian karyotype evolution, reveals constrained gene length evolution in a subset of developmental genes across 3 large genomes. Investigation of limb development genes uncovered the loss of Grem1 in caecilians and Tulp3 in snakes. Interestingly, caecilians and snakes share a significantly larger number of convergent degenerated conserved noncoding elements than limbless lizards, which have a shorter evolutionary history of limb loss. These convergent degenerated conserved noncoding elements overlap significantly with active genomic regions during mouse limb development and are conserved in limbed species, suggesting their essential role in limb patterning in the tetrapod common ancestor. While most convergent degenerated conserved noncoding elements emerged in the jawed vertebrate ancestor, coinciding with the origin of paired appendage, more recent degenerated conserved noncoding elements also contribute to limb development, as demonstrated through functional experiments. Our study provides novel insights into the regulatory elements associated with limb development and loss, offering an evolutionary perspective on the genetic basis of morphological specialization.

The Role of Long Intergenic Noncoding RNA in Fetal Development

The role of long intergenic noncoding RNAs (lincRNAs) in fetal development has emerged as a significant area of study, challenging the traditional protein-centric view of gene expression. While messenger RNAs (mRNAs) have long been recognized for their role in encoding proteins, recent advances have illuminated the critical functions of lincRNAs in various biological processes. Initially identified through high-throughput sequencing technologies, lincRNAs are transcribed from intergenic regions between protein-coding genes and exhibit unique regulatory functions. Unlike mRNAs, lincRNAs are involved in complex interactions with chromatin and chromatin-modifying complexes, influencing gene expression and chromatin structure. LincRNAs are pivotal in regulating tissue-specific development and embryogenesis. For example, they are crucial for proper cardiac, neural, and reproductive system development, with specific lincRNAs being associated with organogenesis and differentiation processes. Their roles in embryonic development include regulating transcription factors and modulating chromatin states, which are essential for maintaining developmental programs and cellular identity. Studies using RNA sequencing and genetic knockout models have highlighted the importance of lincRNAs in processes such as cell differentiation, tissue patterning, and organ development. Despite their functional significance, the comprehensive annotation and understanding of lincRNAs remain limited. Ongoing research aims to elucidate their mechanisms of action and potential applications in disease diagnostics and therapeutics. This review summarizes current knowledge on the functional roles of lincRNAs in fetal development, emphasizing their contributions to tissue-specific gene regulation and developmental processes.

The case against simplistic genetic explanations of evolution

Humans are curious to understand the causes of traits that distinguish us from other animals and that distinguish vastly different species from one another. We also have a proclivity for simple stories and sometimes tend toward seeking and accepting simple genetic explanations for large evolutionary shifts, even to a single gene. Here, I reveal how a biased expectation of mechanistic simplicity threads through the long history of evolutionary and developmental genetics. I argue, however, that expecting a simple mechanism threatens a deeper understanding of evolution, and I define the limitations for interpreting experimental evidence in evolutionary developmental genetics.

Ultraconserved elements from transcriptome and genome data provide insight into the phylogenomics of Sternorrhyncha (Insecta: Hemiptera)

Sternorrhyncha, one of the four major suborders of Hemiptera, is a phytophagous taxon inclusive of nearly 18 000 described species. The phylogenetic relationships within the taxon and the earliest-branching lineage of its infraorders remain incompletely understood. This study attempted to illuminate the phylogenetic relationships within Sternorrhyncha through the use of maximum likelihood, Bayesian inference and maximum parsimony analyses, employing ultraconserved element (UCE) data from 39 genomic and 62 transcriptomic datasets and thereby representing most families within the taxon. The probe set Hemiptera 2.7Kv1 was used to recover a total of 2731 UCE loci: from 547 to 1699 (with an average of 1084) across all genomic datasets and from 108 to 849 (with an average of 329) across all transcriptomic datasets. All three types of phylogenetic analyses employed in this study produced robust statistical support for Sternorrhyncha being a monophyletic group. The different methods of phylogenetic analysis produced inconsistent descriptions of topological structure at the infraorder level: while maximum likelihood and Bayesian inference analyses produced strong statistical evidence (100%) indicating the clade Psylloidea + Aleyrodoidea to be a sister of the clade Aphidoidea (Aphidomorpha) + Coccoidea (Coccomorpha), the maximum parsimony analysis failed to recover a similar result. Our results also provide detail on the phylogenetic relationships within each infraorder. This study presents the first use of UCE data to investigate the phylogeny of Sternorrhyncha. It also shows the viability of amalgamating genomic and transcriptomic data in studies of phylogenetic relationships, potentially highlighting a resource-efficient approach for future inquiries into diverse taxa through the integration of varied data sources.

Non-coding transcribed ultraconserved region uc.290 in colon mucosa promotes intestinal fibrosis in chronic active ulcerative colitis

BACKGROUND AND AIMS

TGF-β1 induces epithelial-mesenchymal transition (EMT) and leads to intestinal fibrosis in ulcerative colitis (UC). We aimed to investigate the expression of transcribed ultraconserved region uc.290 in chronic UC and its role in intestinal fibrosis.

METHODS

Colon specimens were taken from thirty chronic active UC, chronic inactive UC and healthy controls respectively. Modified Mayo score, expressions of uc.290, TGF-β1, EMT biomarkers (Vimentin, α-SMA and E-cadherin) and intestinal fibrosis biomarker (collagen Ⅲ) in colon biopsy specimens were determined in human. Expressions of TGF-β1, EMT biomarkers and collagen Ⅲ were determined in uc.290 overexpressed or silenced epithelial colon cells (HT29).

RESULTS

Uc.290, TGF-β1 and collagen Ⅲ were overexpressed, and EMT was prominent in chronic active UC. Uc.290 level had a positive correlation with modified Mayo score in chronic active UC. TGF-β1 and collagen Ⅲ were overexpressed, and EMT was prominent in uc.290 overexpressed HT29 cells.

CONCLUSIONS

Uc.290 was overexpressed in chronic active UC and might promote intestinal fibrosis by TGF-β1/EMT/collagen Ⅲ pathway.

Genomics-informed captive breeding can reduce inbreeding depression and the genetic load in zoo populations

Zoo populations of threatened species are a valuable resource for the restoration of wild populations. However, their small effective population size poses a risk to long-term viability, especially in species with high genetic load. Recent bioinformatic developments can identify harmful genetic variants in genome data. Here, we advance this approach, analysing the genetic load in the threatened pink pigeon (Nesoenas mayeri). We lifted the mutation-impact scores that had been calculated for the chicken (Gallus gallus) to estimate the genetic load in six pink pigeons. Additionally, we perform in silico crossings to predict the genetic load and realized load of potential offspring. We thus identify the optimal mate pairs that are theoretically expected to produce offspring with the least inbreeding depression. We use computer simulations to show how genomics-informed conservation can reduce the genetic load whilst reducing the loss of genome-wide diversity. Genomics-informed management is likely to become instrumental in maintaining the long-term viability of zoo populations.

Nr2f1 enhancers have distinct functions in controlling Nr2f1 expression during cortical development

There is evidence that transcription factor (TF) encoding genes, which temporally control development in multiple cell types, can have tens of enhancers that regulate their expression. The NR2F1 TF developmentally promotes caudal and ventral cortical regional fates. Here, we epigenomically compared the activity of Nr2f1's enhancers during mouse cortical development with their activity in a transgenic assay. We identified at least six that are likely to be important in prenatal cortical development, with three harboring de novo mutants identified in ASD individuals. We chose to study the function of two of the most robust enhancers by deleting them singly or together. We found that they have distinct and overlapping functions in driving Nr2f1's regional and laminar expression in the developing cortex. Thus, these two enhancers, probably in combination with the others that we defined epigenetically, precisely tune Nr2f1's regional, cell type, and temporal expression during corticogenesis.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Autoregulated splicing of TRA2β programs T cell fate in response to antigen-receptor stimulation

T cell receptor (TCR) sensitivity to peptide-major histocompatibility complex (MHC) dictates T cell fate. Canonical models of TCR sensitivity cannot be fully explained by transcriptional regulation. In this work, we identify a posttranscriptional regulatory mechanism of TCR sensitivity that guides alternative splicing of TCR signaling transcripts through an evolutionarily ultraconserved poison exon (PE) in the RNA-binding protein (RBP) TRA2β in mouse and human. TRA2β-PE splicing, seen during cancer and infection, was required for TCR-induced effector T cell expansion and function. Tra2β-PE skipping enhanced T cell response to antigen by increasing TCR sensitivity. As antigen levels decreased, Tra2β-PE reinclusion allowed T cell survival. Finally, we found that TRA2β-PE was first included in the genome of jawed vertebrates that were capable of TCR gene rearrangements. We propose that TRA2β-PE splicing acts as a gatekeeper of TCR sensitivity to shape T cell fate.

Recombinase-aid amplification combined with lateral flow detection assay for sex identification of the great white pelican (Pelecanus onocrotalus)

Sex identification in avian species is essential for biodiversity conservation and ecological studies. However, the sex of nearly half of the birds could not be identified based on their external appearance. It is difficult to visually identify sex to monitor the ecology and conservation of wild populations. In this study, we designed primer pairs for large white pelican using recombinase-based isothermal amplification combined with a lateral flow dipstick (RAA-LFD) assay for chromo-helicase-DNA binding protein (CHD) genes mapped to W chromosomes and an ultra-conserved element (UCE) located on chromosome 6, respectively. Our result showed that the raaW4-RAA-LFD can detect up to 0.1 ng of genomic DNA (gDNA) templates of female pelicans in 30 min at 39 ℃ and accurately distinguish female from male without any cross reactivity. RaaUCE2-RAA-LFD can amplify both male and female pelicans with a detection limit of 25 pg. To further evaluate the assay, 15 white pelicans of unknown sex were tested using the RAA-LFD assay and conventional polymerase chain reaction (PCR). The results of the raaW4-RAA-LFD assay were consistent with those of the conventional PCR. The developed RAA-LFD assay is equipped with field-deployable instruments and offers a field platform for rapid and reliable sex identification in pelicans.

Computational identification of ultra-conserved elements in the human genome: a hypothesis on homologous DNA pairing

Thousands of prolonged sequences of human ultra-conserved non-coding elements (UCNEs) share only one common feature: peculiarities in the unique composition of their dinucleotides. Here we investigate whether the numerous weak signals emanating from these dinucleotide arrangements can be used for computational identification of UCNEs within the human genome. For this purpose, we analyzed 4272 UCNE sequences, encompassing 1 393 448 nucleotides, alongside equally sized control samples of randomly selected human genomic sequences. Our research identified nine different features of dinucleotide arrangements that enable differentiation of UCNEs from the rest of the genome. We employed these nine features, implementing three Machine Learning techniques - Support Vector Machine, Random Forest, and Artificial Neural Networks - to classify UCNEs, achieving an accuracy rate of 82-84%, with specific conditions allowing for over 90% accuracy. Notably, the strongest feature for UCNE identification was the frequency ratio between GpC dinucleotides and the sum of GpG and CpC dinucleotides. Additionally, we investigated the entire pool of 31 046 SNPs located within UCNEs for their representation in the ClinVar database, which catalogs human SNPs with known phenotypic effects. The presence of UCNE-associated SNPs in ClinVar aligns with the expectation of a random distribution, emphasizing the enigmatic nature of UCNE phenotypic manifestation.

Thermoregulated transcriptomics: the molecular basis and biological significance of temperature-dependent alternative splicing

Temperature-dependent alternative splicing (AS) is a crucial mechanism for organisms to adapt to varying environmental temperatures. In mammals, even slight fluctuations in body temperature are sufficient to drive significant AS changes in a concerted manner. This dynamic regulation allows organisms to finely tune gene expression and protein isoform diversity in response to temperature cues, ensuring proper cellular function and physiological adaptation. Understanding the molecular mechanisms underlying temperature-dependent AS thus provides valuable insights into the intricate interplay between environmental stimuli and gene expression regulation. In this review, we provide an overview of recent advances in understanding temperature-regulated AS across various biological processes and systems. We will discuss the machinery sensing and translating temperature cues into changed AS patterns, the adaptation of the splicing regulatory machinery to extreme temperatures, the role of temperature-dependent AS in shaping the transcriptome, functional implications and the development of potential therapeutics targeting temperature-sensitive AS pathways.

Genomic data provide insights into the classification of extant termites

The higher classification of termites requires substantial revision as the Neoisoptera, the most diverse termite lineage, comprise many paraphyletic and polyphyletic higher taxa. Here, we produce an updated termite classification using genomic-scale analyses. We reconstruct phylogenies under diverse substitution models with ultraconserved elements analyzed as concatenated matrices or within the multi-species coalescence framework. Our classification is further supported by analyses controlling for rogue loci and taxa, and topological tests. We show that the Neoisoptera are composed of seven family-level monophyletic lineages, including the Heterotermitidae Froggatt, Psammotermitidae Holmgren, and Termitogetonidae Holmgren, raised from subfamilial rank. The species-rich Termitidae are composed of 18 subfamily-level monophyletic lineages, including the new subfamilies Crepititermitinae, Cylindrotermitinae, Forficulitermitinae, Neocapritermitinae, Protohamitermitinae, and Promirotermitinae; and the revived Amitermitinae Kemner, Microcerotermitinae Holmgren, and Mirocapritermitinae Kemner. Building an updated taxonomic classification on the foundation of unambiguously supported monophyletic lineages makes it highly resilient to potential destabilization caused by the future availability of novel phylogenetic markers and methods. The taxonomic stability is further guaranteed by the modularity of the new termite classification, designed to accommodate as-yet undescribed species with uncertain affinities to the herein delimited monophyletic lineages in the form of new families or subfamilies.

Illuminating the pathogenic role of SARS-CoV-2: Insights into competing endogenous RNAs (ceRNAs) regulatory networks

The appearance of SARS-CoV-2 in 2019 triggered a significant economic and health crisis worldwide, with heterogeneous molecular mechanisms that contribute to its development are not yet fully understood. Although substantial progress has been made in elucidating the mechanisms behind SARS-CoV-2 infection and therapy, it continues to rank among the top three global causes of mortality due to infectious illnesses. Non-coding RNAs (ncRNAs), being integral components across nearly all biological processes, demonstrate effective importance in viral pathogenesis. Regarding viral infections, ncRNAs have demonstrated their ability to modulate host reactions, viral replication, and host-pathogen interactions. However, the complex interactions of different types of ncRNAs in the progression of COVID-19 remains understudied. In recent years, a novel mechanism of post-transcriptional gene regulation known as "competing endogenous RNA (ceRNA)" has been proposed. Long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and viral ncRNAs function as ceRNAs, influencing the expression of associated genes by sequestering shared microRNAs. Recent research on SARS-CoV-2 has revealed that disruptions in specific ceRNA regulatory networks (ceRNETs) contribute to the abnormal expression of key infection-related genes and the establishment of distinctive infection characteristics. These findings present new opportunities to delve deeper into the underlying mechanisms of SARS-CoV-2 pathogenesis, offering potential biomarkers and therapeutic targets. This progress paves the way for a more comprehensive understanding of ceRNETs, shedding light on the intricate mechanisms involved. Further exploration of these mechanisms holds promise for enhancing our ability to prevent viral infections and develop effective antiviral treatments.

The Evolution of Ultraconserved Elements in Vertebrates

Ultraconserved elements were discovered two decades ago, arbitrarily defined as sequences that are identical over a length ≥ 200 bp in the human, mouse, and rat genomes. The definition was subsequently extended to sequences ≥ 100 bp identical in at least three of five mammalian genomes (including dog and cow), and shown to have undergone rapid expansion from ancestors in fish and strong negative selection in birds and mammals. Since then, many more genomes have become available, allowing better definition and more thorough examination of ultraconserved element distribution and evolutionary history. We developed a fast and flexible analytical pipeline for identifying ultraconserved elements in multiple genomes, dedUCE, which allows manipulation of minimum length, sequence identity, and number of species with a detectable ultraconserved element according to specified parameters. We suggest an updated definition of ultraconserved elements as sequences ≥ 100 bp and ≥97% sequence identity in ≥50% of placental mammal orders (12,813 ultraconserved elements). By mapping ultraconserved elements to ∼200 species, we find that placental ultraconserved elements appeared early in vertebrate evolution, well before land colonization, suggesting that the evolutionary pressures driving ultraconserved element selection were present in aquatic environments in the Cambrian-Devonian periods. Most (>90%) ultraconserved elements likely appeared after the divergence of gnathostomes from jawless predecessors, were largely established in sequence identity by early Sarcopterygii evolution-before the divergence of lobe-finned fishes from tetrapods-and became near fixed in the amniotes. Ultraconserved elements are mainly located in the introns of protein-coding and noncoding genes involved in neurological and skeletomuscular development, enriched in regulatory elements, and dynamically expressed throughout embryonic development.

Patterns of Genetic Diversity and Gene Flow Associated With an Aridity Gradient in Populations of Common Mole-rats, Cryptomys hottentotus hottentotus

Genetic adaptation is the change of a population toward a phenotype that best fits the present ecological conditions of the environment it inhabits. As environmental conditions change, allele frequencies shift, resulting in different populations of the same species possessing genetic variation and divergent phenotypes. Cooperatively breeding common mole-rats (Cryptomys hottentotus hottentotus) inhabit environments along an aridity gradient in South Africa, which provides an opportunity for local genetic adaptations to occur. Using one mitochondrial gene (cytochrome b) and 3,540 SNP loci across the whole genome, we determined the phylogenetic relationship, population structure and genetic diversity of five populations of C. h. hottentotus located along an aridity gradient. Mitochondrial data identified population-specific clades that were less distinct in the two mesic populations, potentially indicating historical or recent gene flow, or the retention of ancestral haplotypes. Arid and semi-arid populations formed a distinct cluster from the non-arid populations. Genetic diversity and gene flow were higher in arid-dwelling individuals, suggesting greater connectivity and interactions between colonies in arid regions in comparison to mesic ones. Using an Aridity Index, we determined that isolation by environment, rather than isolation by geographical distance, best explains the genetic distance between the populations. Further analyses using target loci may determine if there are differing underlying genetic adaptations among populations of C. h. hottentotus. These analyses could help unravel population differences in response to environmental factors within a subspecies of bathyergid mole-rat and determine the adaptive capacity of this small nonmigratory subterranean rodent species in response to aridification in the face of climate change.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Long non-coding RNAs as pathophysiological regulators, therapeutic targets and novel extracellular vesicle biomarkers for the diagnosis of inflammatory bowel disease

Inflammatory bowel disease (IBD) is a chronic relapsing disease of the gastrointestinal (GI) system that includes two groups, Crohn's disease (CD) and ulcerative colitis (UC). To cope with these two classes of IBD, the investigation of pathogenic mechanisms and the discovery of new diagnostic and therapeutic approaches are crucial. Long non-coding RNAs (lncRNAs) which are non-coding RNAs with a length of longer than 200 nucleotides have indicated significant association with the pathology of IBD and strong potential to be used as accurate biomarkers in diagnosing and predicting responses to the IBD treatment. In the current review, we aim to investigate the role of lncRNAs in the pathology and development of IBD. We first describe recent advances in research on dysregulated lncRNAs in the pathogenesis of IBD from the perspective of epithelial barrier function, intestinal immunity, mitochondrial function, and intestinal autophagy. Then, we highlight the possible translational role of lncRNAs as therapeutic targets, diagnostic biomarkers, and predictors of therapeutic response in colon tissues and plasma samples. Finally, we discuss the potential of extracellular vesicles and their lncRNA cargo in the pathophysiology, diagnosis, and treatment of IBD.

DLX genes and proteins in mammalian forebrain development

The vertebrate Dlx gene family encode homeobox transcription factors that are related to the Drosophila Distal-less (Dll) gene and are crucial for development. Over the last ∼35 years detailed information has accrued about the redundant and unique expression and function of the six mammalian Dlx family genes. DLX proteins interact with general transcriptional regulators, and co-bind with other transcription factors to enhancer elements with highly specific activity in the developing forebrain. Integration of the genetic and biochemical data has yielded a foundation for a gene regulatory network governing the differentiation of forebrain GABAergic neurons. In this Primer, we describe the discovery of vertebrate Dlx genes and their crucial roles in embryonic development. We largely focus on the role of Dlx family genes in mammalian forebrain development revealed through studies in mice. Finally, we highlight questions that remain unanswered regarding vertebrate Dlx genes despite over 30 years of research.

The future of rapid and automated single-cell data analysis using reference mapping

As the number of single-cell datasets continues to grow rapidly, workflows that map new data to well-curated reference atlases offer enormous promise for the biological community. In this perspective, we discuss key computational challenges and opportunities for single-cell reference-mapping algorithms. We discuss how mapping algorithms will enable the integration of diverse datasets across disease states, molecular modalities, genetic perturbations, and diverse species and will eventually replace manual and laborious unsupervised clustering pipelines.

Disruption of maternal vascular remodeling by a fetal endoretrovirus-derived gene in preeclampsia

BACKGROUND

Preeclampsia, one of the most lethal pregnancy-related diseases, is associated with the disruption of uterine spiral artery remodeling during placentation. However, the early molecular events leading to preeclampsia remain unknown.

RESULTS

By analyzing placentas from preeclampsia, non-preeclampsia, and twin pregnancies with selective intrauterine growth restriction, we show that the pathogenesis of preeclampsia is attributed to immature trophoblast and maldeveloped endothelial cells. Delayed epigenetic reprogramming during early extraembryonic tissue development leads to generation of excessive immature trophoblast cells. We find reduction of de novo DNA methylation in these trophoblast cells results in selective overexpression of maternally imprinted genes, including the endoretrovirus-derived gene PEG10 (paternally expressed gene 10). PEG10 forms virus-like particles, which are transferred from the trophoblast to the closely proximate endothelial cells. In normal pregnancy, only a low amount of PEG10 is transferred to maternal cells; however, in preeclampsia, excessive PEG10 disrupts maternal vascular development by inhibiting TGF-beta signaling.

CONCLUSIONS

Our study reveals the intricate epigenetic mechanisms that regulate trans-generational genetic conflict and ultimately ensure proper maternal-fetal interface formation.

De novo variants in SP9 cause a novel form of interneuronopathy characterized by intellectual disability, autism spectrum disorder, and epilepsy with variable expressivity

PURPOSE

Interneuronopathies are a group of neurodevelopmental disorders characterized by deficient migration and differentiation of gamma-aminobutyric acidergic interneurons resulting in a broad clinical spectrum, including autism spectrum disorders, early-onset epileptic encephalopathy, intellectual disability, and schizophrenic disorders. SP9 is a transcription factor belonging to the Krüppel-like factor and specificity protein family, the members of which harbor highly conserved DNA-binding domains. SP9 plays a central role in interneuron development and tangential migration, but it has not yet been implicated in a human neurodevelopmental disorder.

METHODS

Cases with SP9 variants were collected through international data-sharing networks. To address the specific impact of SP9 variants, in silico and in vitro assays were carried out.

RESULTS

De novo heterozygous variants in SP9 cause a novel form of interneuronopathy. SP9 missense variants affecting the glutamate 378 amino acid result in severe epileptic encephalopathy because of hypomorphic and neomorphic DNA-binding effects, whereas SP9 loss-of-function variants result in a milder phenotype with epilepsy, developmental delay, and autism spectrum disorder.

CONCLUSION

De novo heterozygous SP9 variants are responsible for a neurodevelopmental disease. Interestingly, variants located in conserved DNA-binding domains of KLF/SP family transcription factors may lead to neomorphic DNA-binding functions resulting in a combination of loss- and gain-of-function effects.

Massively parallel reporter assays and mouse transgenic assays provide complementary information about neuronal enhancer activity

Genetic studies find hundreds of thousands of noncoding variants associated with psychiatric disorders. Massively parallel reporter assays (MPRAs) and in vivo transgenic mouse assays can be used to assay the impact of these variants. However, the relevance of MPRAs to in vivo function is unknown and transgenic assays suffer from low throughput. Here, we studied the utility of combining the two assays to study the impact of non-coding variants. We carried out an MPRA on over 50,000 sequences derived from enhancers validated in transgenic mouse assays and from multiple fetal neuronal ATAC-seq datasets. We also tested over 20,000 variants, including synthetic mutations in highly active neuronal enhancers and 177 common variants associated with psychiatric disorders. Variants with a high impact on MPRA activity were further tested in mice. We found a strong and specific correlation between MPRA and mouse neuronal enhancer activity including changes in neuronal enhancer activity in mouse embryos for variants with strong MPRA effects. Mouse assays also revealed pleiotropic variant effects that could not be observed in MPRA. Our work provides a large catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

Shared evolutionary processes shape landscapes of genomic variation in the great apes

For at least the past 5 decades, population genetics, as a field, has worked to describe the precise balance of forces that shape patterns of variation in genomes. The problem is challenging because modeling the interactions between evolutionary processes is difficult, and different processes can impact genetic variation in similar ways. In this paper, we describe how diversity and divergence between closely related species change with time, using correlations between landscapes of genetic variation as a tool to understand the interplay between evolutionary processes. We find strong correlations between landscapes of diversity and divergence in a well-sampled set of great ape genomes, and explore how various processes such as incomplete lineage sorting, mutation rate variation, GC-biased gene conversion and selection contribute to these correlations. Through highly realistic, chromosome-scale, forward-in-time simulations, we show that the landscapes of diversity and divergence in the great apes are too well correlated to be explained via strictly neutral processes alone. Our best fitting simulation includes both deleterious and beneficial mutations in functional portions of the genome, in which 9% of fixations within those regions is driven by positive selection. This study provides a framework for modeling genetic variation in closely related species, an approach which can shed light on the complex balance of forces that have shaped genetic variation.

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.