Impacts of Neanderthal-Introgressed Sequences on the Landscape of Human Gene Expression.

Regulatory variation influencing gene expression is a key contributor to phenotypic diversity, both within and between species. Unfortunately, RNA degrades too rapidly to be recovered from fossil remains, limiting functional genomic insights about our extinct hominin relatives. Many Neanderthal sequences survive in modern humans due to ancient hybridization, providing an opportunity to assess their contributions to transcriptional variation and to test hypotheses about regulatory evolution. We developed a flexible Bayesian statistical approach to quantify allele-specific expression (ASE) in complex RNA-seq datasets. We identified widespread expression differences between Neanderthal and modern human alleles, indicating pervasive cis-regulatory impacts of introgression. Brain regions and testes exhibited significant downregulation of Neanderthal alleles relative to other tissues, consistent with natural selection influencing the tissue-specific regulatory landscape. Our study demonstrates that Neanderthal-inherited sequences are not silent remnants of ancient interbreeding but have measurable impacts on gene expression that contribute to variation in modern human phenotypes.

Opportunities and tradeoffs in single-cell transcriptomic technologies.

Recent technological and algorithmic advances enable single-cell transcriptomic analysis with remarkable depth and breadth. Nonetheless, a persistent challenge is the compromise between the ability to profile high numbers of cells and the achievement of full-length transcript coverage. Currently, the field is progressing and developing new and creative solutions that improve cellular throughput, gene detection sensitivity and full-length transcript capture. Furthermore, long-read sequencing approaches for single-cell transcripts are breaking frontiers that have previously blocked full transcriptome characterization. We here present a comprehensive overview of available options for single-cell transcriptome profiling, highlighting the key advantages and disadvantages of each approach.

Ornaments for efficient allele-specific expression estimation with bias correction.

Allele-specific expression plays a crucial role in unraveling various biological mechanisms, including genomic imprinting and gene expression controlled by cis-regulatory variants. However, existing methods for quantification from RNA-sequencing (RNA-seq) reads do not adequately and efficiently remove various allele-specific read mapping biases, such as reference bias arising from reads containing the alternative allele that do not map to the reference transcriptome or ambiguous mapping bias caused by reads containing the reference allele that map differently from reads containing the alternative allele. We present Ornaments, a computational tool for rapid and accurate estimation of allele-specific transcript expression at unphased heterozygous loci from RNA-seq reads while correcting for allele-specific read mapping biases. Ornaments removes reference bias by mapping reads to a personalized transcriptome and ambiguous mapping bias by probabilistically assigning reads to multiple transcripts and variant loci they map to. Ornaments is a lightweight extension of kallisto, a popular tool for fast RNA-seq quantification, that improves the efficiency and accuracy of WASP, a popular tool for bias correction in allele-specific read mapping. In experiments with simulated and human lymphoblastoid cell-line RNA-seq reads with the genomes of the 1000 Genomes Project, we demonstrate that Ornaments improves the accuracy of WASP and kallisto, is nearly as efficient as kallisto, and is an order of magnitude faster than WASP per sample, with the additional cost of constructing a personalized index for multiple samples. Additionally, we show that Ornaments finds imprinted transcripts with higher sensitivity than WASP, which detects imprinted signals only at gene level.

Long-read transcriptome sequencing reveals allele-specific variants at high resolution.

Disturbance in the regulation of transcript structure plays a crucial role in human disease. In a recent study, Glinos et al. characterized allele-specific transcript alterations in long-read RNA sequencing (RNA-seq) data derived from multiple human tissues and provide a high-resolution view of how disease-associated genetic variants affect transcript structure.

Genome-wide chromatin accessibility reveals transcriptional regulation of heterosis in inter-subspecific hybrid rice.

The utilization of rice heterosis is essential for ensuring global food security; however, its molecular mechanism remains unclear. In this study, comprehensive analyses of accessible chromatin regions (ACRs), DNA methylation, and gene expression in inter-subspecific hybrid and its parents were performed to determine the potential role of chromatin accessibility in rice heterosis. The hybrid exhibited abundant ACRs, in which the gene ACRs and proximal ACRs were directly related to transcriptional activation rather than the distal ACRs. Regarding the dynamic accessibility contribution of the parents, paternal ZHF1015 transmitted a greater number of ACRs to the hybrid. Accessible genotype-specific target genes were enriched with overrepresented transcription factors, indicating a unique regulatory network of genes in the hybrid. Compared with its parents, the differentially accessible chromatin regions with upregulated chromatin accessibility were much greater than those with downregulated chromatin accessibility, reflecting a stronger regulation in the hybrid. Furthermore, DNA methylation levels were negatively correlated with ACR intensity, and genes were strongly affected by CHH methylation in the hybrid. Chromatin accessibility positively regulated the overall expression level of each genotype. ACR-related genes with maternal Z04A-bias allele-specific expression tended to be enriched during carotenoid biosynthesis, whereas paternal ZHF1015-bias genes were more active in carbohydrate metabolism. Our findings provide a new perspective on the mechanism of heterosis based on chromatin accessibility in inter-subspecific hybrid rice.

Gene and Allele-Specific Expression Underlying the Electric Signal Divergence in African Weakly Electric Fish.

In the African weakly electric fish genus Campylomormyrus, electric organ discharge signals are strikingly different in shape and duration among closely related species, contribute to prezygotic isolation, and may have triggered an adaptive radiation. We performed mRNA sequencing on electric organs and skeletal muscles (from which the electric organs derive) from 3 species with short (0.4 ms), medium (5 ms), and long (40 ms) electric organ discharges and 2 different cross-species hybrids. We identified 1,444 upregulated genes in electric organ shared by all 5 species/hybrid cohorts, rendering them candidate genes for electric organ-specific properties in Campylomormyrus. We further identified several candidate genes, including KCNJ2 and KLF5, and their upregulation may contribute to increased electric organ discharge duration. Hybrids between a short (Campylomormyrus compressirostris) and a long (Campylomormyrus rhynchophorus) discharging species exhibit electric organ discharges of intermediate duration and showed imbalanced expression of KCNJ2 alleles, pointing toward a cis-regulatory difference at this locus, relative to electric organ discharge duration. KLF5 is a transcription factor potentially balancing potassium channel gene expression, a crucial process for the formation of an electric organ discharge. Unraveling the genetic basis of the species-specific modulation of the electric organ discharge in Campylomormyrus is crucial for understanding the adaptive radiation of this emerging model taxon of ecological (perhaps even sympatric) speciation.

Machine learning on alignment features for parent-of-origin classification of simulated hybrid RNA-seq.

BACKGROUND: Parent-of-origin allele-specific gene expression (ASE) can be detected in interspecies hybrids by virtue of RNA sequence variants between the parental haplotypes. ASE is detectable by differential expression analysis (DEA) applied to the counts of RNA-seq read pairs aligned to parental references, but aligners do not always choose the correct parental reference. RESULTS: We used public data for species that are known to hybridize. We measured our ability to assign RNA-seq read pairs to their proper transcriptome or genome references. We tested software packages that assign each read pair to a reference position and found that they often favored the incorrect species reference. To address this problem, we introduce a post process that extracts alignment features and trains a random forest classifier to choose the better alignment. On each simulated hybrid dataset tested, our machine-learning post-processor achieved higher accuracy than the aligner by itself at choosing the correct parent-of-origin per RNA-seq read pair. CONCLUSIONS: For the parent-of-origin classification of RNA-seq, machine learning can improve the accuracy of alignment-based methods. This approach could be useful for enhancing ASE detection in interspecies hybrids, though RNA-seq from real hybrids may present challenges not captured by our simulations. We believe this is the first application of machine learning to this problem domain.

Natural variations of heterosis-related allele-specific expression genes in promoter regions lead to allele-specific expression in maize.

BACKGROUND: Heterosis has successfully enhanced maize productivity and quality. Although significant progress has been made in delineating the genetic basis of heterosis, the molecular mechanisms underlying its genetic components remain less explored. Allele-specific expression (ASE), the imbalanced expression between two parental alleles in hybrids, is increasingly being recognized as a factor contributing to heterosis. ASE is a complex process regulated by both epigenetic and genetic variations in response to developmental and environmental conditions. RESULTS: In this study, we explored the differential characteristics of ASE by analyzing the transcriptome data of two maize hybrids and their parents under four light conditions. On the basis of allele expression patterns in different hybrids under various conditions, ASE genes were divided into three categories: bias-consistent genes involved in basal metabolic processes in a functionally complementary manner, bias-reversal genes adapting to the light environment, and bias-specific genes maintaining cell homeostasis. We observed that 758 ASE genes (ASEGs) were significantly overlapped with heterosis quantitative trait loci (QTLs), and high-frequency variations in the promoter regions of heterosis-related ASEGs were identified between parents. In addition, 10 heterosis-related ASEGs participating in yield heterosis were selected during domestication. CONCLUSIONS: The comprehensive analysis of ASEGs offers a distinctive perspective on how light quality influences gene expression patterns and gene-environment interactions, with implications for the identification of heterosis-related ASEGs to enhance maize yield.

Lineage-specific selection and the evolution of virulence in the Candida clade.

Candida albicans is the most common cause of systemic fungal infections in humans and is considerably more virulent than its closest known relative, Candida dubliniensis. To investigate this difference, we constructed interspecies hybrids and quantified mRNA levels produced from each genome in the hybrid. This approach systematically identified expression differences in orthologous genes arising from cis-regulatory sequence changes that accumulated since the two species last shared a common ancestor, some 10 million y ago. We documented many orthologous gene-expression differences between the two species, and we pursued one striking observation: All 15 genes coding for the enzymes of glycolysis showed higher expression from the C. albicans genome than the C. dubliniensis genome in the interspecies hybrid. This pattern requires evolutionary changes to have occurred at each gene; the fact that they all act in the same direction strongly indicates lineage-specific natural selection as the underlying cause. To test whether these expression differences contribute to virulence, we created a C. dubliniensis strain in which all 15 glycolysis genes were produced at modestly elevated levels and found that this strain had significantly increased virulence in the standard mouse model of systemic infection. These results indicate that small expression differences across a deeply conserved set of metabolism enzymes can play a significant role in the evolution of virulence in fungal pathogens.

Parental bias in expression and interaction of genes in the equine placenta.

Most autosomal genes in the placenta show a biallelic expression pattern. However, some genes exhibit allele-specific transcription depending on the parental origin of the chromosomes on which the copy of the gene resides. Parentally expressed genes are involved in the reciprocal interaction between maternal and paternal genes, coordinating the allocation of resources between fetus and mother. One of the main challenges of studying parental-specific allelic expression (allele-specific expression [ASE]) in the placenta is the maternal cellular remnant at the fetomaternal interface. Horses (Equus caballus) have an epitheliochorial placenta in which both the endometrial epithelium and the epithelium of the chorionic villi are juxtaposed with minimal extension into the uterine mucosa, yet there is no information available on the allelic gene expression of equine chorioallantois (CA). In the current study, we present a dataset of 1,336 genes showing ASE in the equine CA (https://pouya-dini.github.io/equine-gene-db/) along with a workflow for analyzing ASE genes. We further identified 254 potentially imprinted genes among the parentally expressed genes in the equine CA and evaluated the expression pattern of these genes throughout gestation. Our gene ontology analysis implies that maternally expressed genes tend to decrease the length of gestation, while paternally expressed genes extend the length of gestation. This study provides fundamental information regarding parental gene expression during equine pregnancy, a species with a negligible amount of maternal cellular remnant in its placenta. This information will provide the basis for a better understanding of the role of parental gene expression in the placenta during gestation.