Human genomics. The human transcriptome across tissues and individuals

Whole-genome sequencing analysis reveals new susceptibility loci and structural variants associated with progressive supranuclear palsy

BACKGROUND

Progressive supranuclear palsy (PSP) is a rare neurodegenerative disease characterized by the accumulation of aggregated tau proteins in astrocytes, neurons, and oligodendrocytes. Previous genome-wide association studies for PSP were based on genotype array, therefore, were inadequate for the analysis of rare variants as well as larger mutations, such as small insertions/deletions (indels) and structural variants (SVs).

METHOD

In this study, we performed whole genome sequencing (WGS) and conducted association analysis for single nucleotide variants (SNVs), indels, and SVs, in a cohort of 1,718 cases and 2,944 controls of European ancestry. Of the 1,718 PSP individuals, 1,441 were autopsy-confirmed and 277 were clinically diagnosed.

RESULTS

Our analysis of common SNVs and indels confirmed known genetic loci at MAPT, MOBP, STX6, SLCO1A2, DUSP10, and SP1, and further uncovered novel signals in APOE, FCHO1/MAP1S, KIF13A, TRIM24, TNXB, and ELOVL1. Notably, in contrast to Alzheimer's disease (AD), we observed the APOE ε2 allele to be the risk allele in PSP. Analysis of rare SNVs and indels identified significant association in ZNF592 and further gene network analysis identified a module of neuronal genes dysregulated in PSP. Moreover, seven common SVs associated with PSP were observed in the H1/H2 haplotype region (17q21.31) and other loci, including IGH, PCMT1, CYP2A13, and SMCP. In the H1/H2 haplotype region, there is a burden of rare deletions and duplications (P = 6.73 × 10-3) in PSP.

CONCLUSIONS

Through WGS, we significantly enhanced our understanding of the genetic basis of PSP, providing new targets for exploring disease mechanisms and therapeutic interventions.

Phytocannabinoids in neuromodulation: From omics to epigenetics

BACKGROUND

Recent developments in metabolomics, transcriptomic and epigenetics open up new horizons regarding the pharmacological understanding of phytocannabinoids as neuromodulators in treating anxiety, depression, epilepsy, Alzheimer's, Parkinson's disease and autism.

METHODS

The present review is an extensive search in public databases, such as Google Scholar, Scopus, the Web of Science, and PubMed, to collect all the literature about the neurobiological roles of cannabis extract, cannabidiol, 9-tetrahydrocannabinol specially focused on metabolomics, transcriptomic, epigenetic, mechanism of action, in different cell lines, induced animal models and clinical trials. We used bioinformatics, network pharmacology and enrichment analysis to understand the effect of phytocannabinoids in neuromodulation.

RESULTS

Cannabidomics studies show wide variability of metabolites across different strains and varieties, which determine their medicinal and abusive usage, which is very important for its quality control and regulation. CB receptors interact with other compounds besides cannabidiol and Δ9-tetrahydrocannabinol, like cannabinol and Δ8-tetrahydrocannabinol. Phytocannabinoids interact with cannabinoid and non-cannabinoid receptors (GPCR, ion channels, and PPAR) to improve various neurodegenerative diseases. However, its abuse because of THC is also a problem found across different epigenetic and transcriptomic studies. Network enrichment analysis shows CNR1 expression in the brain and its interacting genes involve different pathways such as Rap1 signalling, dopaminergic synapse, and relaxin signalling. CBD protects against diseases like epilepsy, depression, and Parkinson's by modifying DNA and mitochondrial DNA in the hippocampus. Network pharmacology analysis of 8 phytocannabinoids revealed an interaction with 10 (out of 60) targets related to neurodegenerative diseases, with enrichment of ErbB and PI3K-Akt signalling pathways which helps in ameliorating neuro-inflammation in various neurodegenerative diseases. The effects of phytocannabinoids vary across sex, disease state, and age which suggests the importance of a personalized medicine approach for better success.

CONCLUSIONS

Phytocannabinoids present a range of promising neuromodulatory effects. It holds promise if utilized in a strategic way towards personalized neuropsychiatric treatment. However, just like any drug irrational usage may lead to unforeseen negative effects. Exploring neuro-epigenetics and systems pharmacology of major and minor phytocannabinoid combinations can lead to success.

Integrative multi-omic analysis reveals conserved cell-projection deficits in human Down syndrome brains

Down syndrome (DS) is the most common genetic cause of cognitive disability. However, it is largely unclear how triplication of a small gene subset may impinge on diverse aspects of DS brain physiopathology. Here, we took a multi-omic approach and simultaneously analyzed by RNA-seq and proteomics the expression signatures of two diverse regions of human postmortem DS brains. We found that the overexpression of triplicated genes triggered global expression dysregulation, differentially affecting transcripts, miRNAs, and proteins involved in both known and novel biological candidate pathways. Among the latter, we observed an alteration in RNA splicing, specifically modulating the expression of genes involved in cytoskeleton and axonal dynamics in DS brains. Accordingly, we found an alteration in axonal polarization in neurons from DS human iPSCs and mice. Thus, our study provides an integrated multilayer expression database capable of identifying new potential targets to aid in designing future clinical interventions for DS.

Haplotype-resolved assembly of the mule duck genome using high-fidelity sequencing technology

Mule duck is vitally important to the production of global duck meat. Here, we present two high-quality haplotypes of a female mule duck (haplotype 1 (H1):1.28 Gb, haplotype 2 (H2): 1.40 Gb). The continuity (H1: contig N50 = 14.90 Mb, H2: contig N50 = 15.70 Mb) and completeness (BUSCO: H1 = 96.9%, H2 = 97.3%) are substantially better than those of other duck genomes. We detected the structural variations (SVs) in H1 and H2. We observed a positive correlation between autosome length and the number of SVs. Z chromosome was some deficient in deletions and insertions, but W chromosome was some excessive. A total of 1,451 genes were haplotype specific expression (HSEs). Among them, 737 specifically expressed in H1, and 714 specifically expressed in H2. We found that H1 and H2 HSEs tended to be involved in similar biological processes, such as myometrial relaxation and contraction pathways, muscle structure development and phosphorylation. Our haplotype-resolved genome assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in mule duck.

Loss of Function of Vasoactive-intestinal Peptide Alters Sex Ratio and Reduces Male Reproductive Fitness in Zebrafish

Vasoactive-intestinal peptide (Vip) is a pleiotropic peptide with a wide range of distribution and functions. Zebrafish possess 2 isoforms of Vip (a and b), in which Vipa is most homologous to the mammalian form. In female zebrafish, Vipa can stimulate LH secretion from the pituitary but is not essential for female reproduction, as vipa-/- females display normal reproduction. In contrast, we have found that vipa-/- males are severely subfertile and sex ratio of offspring is female-biased. By analyzing all aspects of male reproduction with wild-type (WT) males, we show that the testes of vipa-/- are underdeveloped and contain ∼70% less spermatids compared to WT counterparts. The sperm of vipa-/- males displayed reduced potency in terms of fertilization (by ∼80%) and motility span and duration (by ∼50%). In addition, vipa-/- male attraction to WT females was largely nonexistent, indicating decreased sexual motivation. We show that vipa mRNA and protein is present in Leydig cells and in developing germ cells in the testis of WT, raising the possibility that endogenous Vipa contributes to testicular function. Absence of Vipa in vipa-/- males resulted in downregulation of 3 key genes in the androgen synthesis chain in the testis, 3β-hsd, 17β-hsd1, and cyp11c1 (11β-hydrogenase), associated with a pronounced decrease in 11-ketotestosterone production and, in turn, compromised reproductive fitness. Altogether, this study establishes a crucial role for Vipa in the regulation of male reproduction in zebrafish, like in mammals, with the exception that Vipa is also expressed in zebrafish testis.

Histamine H3 Receptor Isoforms: Insights from Alternative Splicing to Functional Complexity

Alternative splicing significantly enhances the diversity of the G protein-coupled receptor (GPCR) family, including the histamine H3 receptor (H3R). This post-transcriptional modification generates multiple H3R isoforms with potentially distinct pharmacological and physiological profiles. H3R is primarily involved in the presynaptic inhibition of neurotransmitter release in the central nervous system. Despite the approval of pitolisant for narcolepsy (Wakix®) and daytime sleepiness in adults with obstructive sleep apnea (Ozawade®) and ongoing clinical trials for other H3R antagonists/inverse agonists, the functional significance of the numerous H3R isoforms remains largely enigmatic. Recent publicly available RNA sequencing data have confirmed the expression of multiple H3R isoforms in the brain, with some isoforms exhibiting unique tissue-specific distribution patterns hinting at isoform-specific functions and interactions within neural circuits. In this review, we discuss the complexity of H3R isoforms with a focus on their potential roles in central nervous system (CNS) function. Comparative analysis across species highlights evolutionary conservation and divergence in H3R splicing, suggesting species-specific regulatory mechanisms. Understanding the functionality of H3R isoforms is crucial for the development of targeted therapeutics. This knowledge will inform the design of more precise pharmacological interventions, potentially enhancing therapeutic efficacy and reducing adverse effects in the treatment of neurological and psychiatric disorders.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From > 85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

Biallelic variants in Plexin B2 (PLXNB2) cause amelogenesis imperfecta, hearing loss and intellectual disability

BACKGROUND

Plexins are large transmembrane receptors for the semaphorin family of signalling proteins. Semaphorin-plexin signalling controls cellular interactions that are critical during development as well as in adult life stages. Nine plexin genes have been identified in humans, but despite the apparent importance of plexins in development, only biallelic PLXND1 and PLXNA1 variants have so far been associated with Mendelian genetic disease.

METHODS

Eight individuals from six families presented with a recessively inherited variable clinical condition, with core features of amelogenesis imperfecta (AI) and sensorineural hearing loss (SNHL), with variable intellectual disability. Probands were investigated by exome or genome sequencing. Common variants and those unlikely to affect function were excluded. Variants consistent with autosomal recessive inheritance were prioritised. Variant segregation analysis was performed by Sanger sequencing. RNA expression analysis was conducted in C57Bl6 mice.

RESULTS

Rare biallelic pathogenic variants in plexin B2 (PLXNB2), a large transmembrane semaphorin receptor protein, were found to segregate with disease in all six families. The variants identified include missense, nonsense, splicing changes and a multiexon deletion. Plxnb2 expression was detected in differentiating ameloblasts.

CONCLUSION

We identify rare biallelic pathogenic variants in PLXNB2 as a cause of a new autosomal recessive, phenotypically diverse syndrome with AI and SNHL as core features. Intellectual disability, ocular disease, ear developmental abnormalities and lymphoedema were also present in multiple cases. The variable syndromic human phenotype overlaps with that seen in Plxnb2 knockout mice, and, together with the rarity of human PLXNB2 variants, may explain why pathogenic variants in PLXNB2 have not been reported previously.

Comparison of RNA-Sequencing Methods for Degraded RNA

RNA sequencing (RNA-Seq) is a powerful technique and is increasingly being used in clinical research and drug development. Currently, several RNA-Seq methods have been developed. However, the relative advantage of each method for degraded RNA and low-input RNA, such as RNA samples collected in the field of clinical setting, has remained unknown. The Standard method of RNA-Seq captures mRNA by poly(A) capturing using Oligo dT beads, which is not suitable for degraded RNA. Here, we used three commercially available RNA-Seq library preparation kits (SMART-Seq, xGen Broad-range, and RamDA-Seq) using random primer instead of Oligo dT beads. To evaluate the performance of these methods, we compared the correlation, the number of detected expressing genes, and the expression levels with the Standard RNA-Seq method. Although the performance of RamDA-Seq was similar to that of Standard RNA-Seq, the performance for low-input RNA and degraded RNA has decreased. The performance of SMART-Seq was better than xGen and RamDA-Seq in low-input RNA and degraded RNA. Furthermore, the depletion of ribosomal RNA (rRNA) improved the performance of SMART-Seq and xGen due to increased expression levels. SMART-Seq with rRNA depletion has relative advantages for RNA-Seq using low-input and degraded RNA.

Housekeeping protein-coding genes interrogated with tissue and individual variations

Housekeeping protein-coding genes are stably expressed genes in cells and tissues that are thought to be engaged in fundamental cellular biological functions. They are often utilized as normalization references in molecular biology research and are especially important in integrated bioinformatic investigations. Prior studies have examined human housekeeping protein-coding genes by analyzing various gene expression datasets. The inclusion of different tissue types significantly impacted the discovery of housekeeping genes. In this report, we investigated particularly individual human subject expression differences in protein-coding genes across different tissue types. We used GTEx V8 gene expression datasets obtained from more than 16,000 human normal tissue samples. Furthermore, the Gini index is utilized to investigate the expression variations of protein-coding genes between tissue and individual donor subjects. Housekeeping protein-coding genes found using Gini index profiles may vary depending on the tissue subtypes investigated, particularly given the diverse sample size collections across the GTEx tissue subtypes. We subsequently selected major tissues and identified subsets of housekeeping genes with stable expression levels among human donors within those tissues. In this work, we provide alternative sets of housekeeping protein-coding genes that show more consistent expression patterns in human subjects across major solid organs. Weblink: https://hpsv.ibms.sinica.edu.tw .

Long noncoding RNA GATA2AS influences human erythropoiesis by transcription factor and chromatin landscape modulation

ABSTRACT

Long noncoding RNAs (lncRNAs) are extensively expressed in eukaryotic cells and have been revealed to be important for regulating cell differentiation. Many lncRNAs have been found to regulate erythroid differentiation in the mouse. However, given the low sequence conservation of lncRNAs between mouse and human, our understanding of lncRNAs in human erythroid differentiation remains incomplete. lncRNAs are often transcribed opposite to protein coding genes and regulate their expression. Here, we characterized a human erythrocyte-expressed lncRNA, GATA2AS, which is transcribed opposite to erythroid transcription regulator GATA2. GATA2AS is a 2080-bp long, primarily nucleus-localized noncoding RNA that is expressed in erythroid progenitor cells and decreases during differentiation. Knockout of GATA2AS in human HUDEP2 erythroid progenitor cells using CRISPR-Cas9 genome editing to remove the transcription start site accelerated erythroid differentiation and dysregulated erythroblast gene expression. We identified GATA2AS as a novel GATA2 and HBG activator. Chromatin isolation by RNA purification showed that GATA2AS binds to thousands of genomic sites and colocalizes at a subset of sites with erythroid transcription factors including LRF and KLF1. RNA pulldown and RNA immunoprecipitation confirmed interaction between GATA2AS and LRF and KLF1. Chromatin immunoprecipitation sequencing (ChIP-seq) showed that knockout of GATA2AS reduces binding of these transcription factors genome wide. Assay for transposase-accessible chromatin sequencing (ATAC-seq) and H3K27ac ChIP-seq showed that GATA2AS is essential to maintain the chromatin regulatory landscape during erythroid differentiation. Knockdown of GATA2AS in human primary CD34+ cells mimicked results in HUDEP2 cells. Overall, our results implicate human-specific lncRNA GATA2AS as a regulator of erythroid differentiation by influencing erythroid transcription factor binding and the chromatin regulatory landscape.

Prediction of tumor-specific splicing from somatic mutations as a source of neoantigen candidates

Motivation

Neoantigens are promising targets for cancer immunotherapies and might arise from alternative splicing. However, detecting tumor-specific splicing is challenging because many non-canonical splice junctions identified in tumors also appear in healthy tissues. To increase tumor-specificity, we focused on splicing caused by somatic mutations as a source for neoantigen candidates in individual patients.

Results

We developed the tool splice2neo with multiple functionalities to integrate predicted splice effects from somatic mutations with splice junctions detected in tumor RNA-seq and to annotate the resulting transcript and peptide sequences. Additionally, we provide the tool EasyQuant for targeted RNA-seq read mapping to candidate splice junctions. Using a stringent detection rule, we predicted 1.7 splice junctions per patient as splice targets with a false discovery rate below 5% in a melanoma cohort. We confirmed tumor-specificity using independent, healthy tissue samples. Furthermore, using tumor-derived RNA, we confirmed individual exon-skipping events experimentally. Most target splice junctions encoded neoepitope candidates with predicted major histocompatibility complex (MHC)-I or MHC-II binding. Compared to neoepitope candidates from non-synonymous point mutations, the splicing-derived MHC-I neoepitope candidates had lower self-similarity to corresponding wild-type peptides. In conclusion, we demonstrate that identifying mutation-derived, tumor-specific splice junctions can lead to additional neoantigen candidates to expand the target repertoire for cancer immunotherapies.

Availability and implementation

The R package splice2neo and the python package EasyQuant are available at https://github.com/TRON-Bioinformatics/splice2neo and https://github.com/TRON-Bioinformatics/easyquant, respectively.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

CircCamsap1 is dispensable for male fertility in mice

Background

Circular RNAs (circRNAs) are a large class of RNAs present in mammals. Among these, circCamsap1 is a well-acknowledged circRNA with significant implications, particularly in the development and progression of diverse tumors. However, the potential consequences of circCamsap1 depletion in vivo on male reproduction are yet to be thoroughly investigated.

Methods

The presence of circCamsap1 in the mouse testes was confirmed, and gene expression analysis was performed using reverse transcription quantitative polymerase chain reaction. CircCamsap1 knockout mice were generated utilizing the CRISPR/Cas9 system. Phenotypic analysis of both the testes and epididymis was conducted using histological and immunofluorescence staining. Additionally, fertility and sperm motility were assessed.

Results

Here, we successfully established a circCamsap1 knockout mouse model without affecting the expression of parental gene. Surprisingly, male mice lacking circCamsap1 (circCamsap1-/-) exhibited normal fertility, with no discernible differences in testicular and epididymal histology, spermatogenesis, sperm counts or sperm motility compared to circCamsap1+/+ mice. These findings suggest that circCamsap1 may not play an essential role in physiological spermatogenesis. Nonetheless, this result also underscores the complexity of circRNA function in male reproductive biology. Therefore, further research is necessary to elucidate the precise roles of other circRNAs in regulating male fertility.

Sex affects transcriptional associations with schizophrenia across the dorsolateral prefrontal cortex, hippocampus, and caudate nucleus

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

The genetic regulation of the gastric transcriptome is associated with metabolic and obesity-related traits and diseases

Tissue-specific gene expression and gene regulation lead to a better understanding of tissue-specific physiology and pathophysiology. We analyzed the transcriptome and genetic regulatory profiles of two distinct gastric sites, corpus and antrum, to identify tissue-specific gene expression and its regulation. Gastric corpus and antrum mucosa biopsies were collected during routine gastroscopies from up to 431 healthy individuals. We obtained genotype and transcriptome data and performed transcriptome profiling and expression quantitative trait locus (eQTL) studies. We further used data from genome-wide association studies (GWAS) of various diseases and traits to partition their heritability and to perform transcriptome-wide association studies (TWAS). The transcriptome data from corpus and antral mucosa highlights the heterogeneity of gene expression in the stomach. We identified enriched pathways revealing distinct and common physiological processes in gastric corpus and antrum. Furthermore, we found an enrichment of the single nucleotide polymorphism (SNP)-based heritability of metabolic, obesity-related, and cardiovascular traits and diseases by considering corpus- and antrum-specifically expressed genes. Particularly, we could prioritize gastric-specific candidate genes for multiple metabolic traits, like NQO1 which is involved in glucose metabolism, MUC1 which contributes to purine and protein metabolism or RAB27B being a regulator of weight and body composition. Our findings show that gastric corpus and antrum vary in their transcriptome and genetic regulatory profiles indicating physiological differences which are mostly related to digestion and epithelial protection. Moreover, our findings demonstrate that the genetic regulation of the gastric transcriptome is linked to biological mechanisms associated with metabolic, obesity-related, and cardiovascular traits and diseases. NEW & NOTEWORTHY We analyzed the transcriptomes and genetic regulatory profiles of gastric corpus and for the first time also of antrum mucosa in 431 healthy individuals. Through tissue-specific gene expression and eQTL analyses, we uncovered unique and common physiological processes across both primary gastric sites. Notably, our findings reveal that stomach-specific eQTLs are enriched in loci associated with metabolic traits and diseases, highlighting the pivotal role of gene expression regulation in gastric physiology and potential pathophysiology.

Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits

BACKGROUND

Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits.

RESULTS

In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits.

CONCLUSIONS

Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production traits.

Comprehensive transcriptome analysis reveals altered mRNA splicing and post-transcriptional changes in the aged mouse brain

A comprehensive understanding of molecular changes during brain aging is essential to mitigate cognitive decline and delay neurodegenerative diseases. The interpretation of mRNA alterations during brain aging is influenced by the health and age of the animal cohorts studied. Here, we carefully consider these factors and provide an in-depth investigation of mRNA splicing and dynamics in the aging mouse brain, combining short- and long-read sequencing technologies with extensive bioinformatic analyses. Our findings encompass a spectrum of age-related changes, including differences in isoform usage, decreased mRNA dynamics and a module showing increased expression of neuronal genes. Notably, our results indicate a reduced abundance of mRNA isoforms leading to nonsense-mediated RNA decay and suggest a regulatory role for RNA-binding proteins, indicating that their regulation may be altered leading to the reshaping of the aged brain transcriptome. Collectively, our study highlights the importance of studying mRNA splicing events during brain aging.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Time is ticking faster for long genes in aging

Recent studies of aging organisms have identified a systematic phenomenon, characterized by a negative correlation between gene length and their expression in various cell types, species, and diseases. We term this phenomenon gene-length-dependent transcription decline (GLTD) and suggest that it may represent a bottleneck in the transcription machinery and thereby significantly contribute to aging as an etiological factor. We review potential links between GLTD and key aging processes such as DNA damage and explore their potential in identifying disease modification targets. Notably, in Alzheimer's disease, GLTD spotlights extremely long synaptic genes at chromosomal fragile sites (CFSs) and their vulnerability to postmitotic DNA damage. We suggest that GLTD is an integral element of biological aging.

Sex-biased gene regulation varies across human populations as a result of adaptive evolution

OBJECTIVES

Studies of human sexual dimorphism and gender disparities in health focus on ostensibly universal molecular sex differences, such as sex chromosomes and circulating hormone levels, while ignoring the extraordinary diversity in biology, behavior, and culture acquired by different human populations over their unique evolutionary histories.

MATERIALS AND METHODS

Using RNA-Seq data and whole genome sequences from 1000G and HGDP, we investigate variation in sex-biased gene expression across 11 human populations and test whether population-level variation in sex-biased expression may have resulted from adaptive evolution in regions containing sex-specific regulatory variants.

RESULTS

We find that sex-biased gene expression in humans is highly variable, mostly population-specific, and demonstrates between population reversals. Expression quantitative trait locus mapping reveals sex-specific regulatory regions with evidence of recent positive natural selection, suggesting that variation in sex-biased expression may have evolved as an adaptive response to ancestral environments experienced by human populations.

DISCUSSION

These results indicate that sex-biased gene expression is more flexible than previously thought and is not generally shared among human populations. Instead, molecular phenotypes associated with sex depend on complex interactions between population-specific molecular evolution and physiological responses to contemporary socioecologies.

Identification of Interpretable Clusters and Associated Signatures in Breast Cancer Single-Cell Data: A Topic Modeling Approach

Topic modeling is a popular technique in machine learning and natural language processing, where a corpus of text documents is classified into themes or topics using word frequency analysis. This approach has proven successful in various biological data analysis applications, such as predicting cancer subtypes with high accuracy and identifying genes, enhancers, and stable cell types simultaneously from sparse single-cell epigenomics data. The advantage of using a topic model is that it not only serves as a clustering algorithm, but it can also explain clustering results by providing word probability distributions over topics. Our study proposes a novel topic modeling approach for clustering single cells and detecting topics (gene signatures) in single-cell datasets that measure multiple omics simultaneously. We applied this approach to examine the transcriptional heterogeneity of luminal and triple-negative breast cancer cells using patient-derived xenograft models with acquired resistance to chemotherapy and targeted therapy. Through this approach, we identified protein-coding genes and long non-coding RNAs (lncRNAs) that group thousands of cells into biologically similar clusters, accurately distinguishing drug-sensitive and -resistant breast cancer types. In comparison to standard state-of-the-art clustering analyses, our approach offers an optimal partitioning of genes into topics and cells into clusters simultaneously, producing easily interpretable clustering outcomes. Additionally, we demonstrate that an integrative clustering approach, which combines the information from mRNAs and lncRNAs treated as disjoint omics layers, enhances the accuracy of cell classification.

Genome Sequencing of Idiopathic Speech Delay

Genetic investigations of people with speech and language disorders can provide windows into key aspects of human biology. Most genomic research into impaired speech development has so far focused on childhood apraxia of speech (CAS), a rare neurodevelopmental disorder characterized by difficulties with coordinating rapid fine motor sequences that underlie proficient speech. In 2001, pathogenic variants of FOXP2 provided the first molecular genetic accounts of CAS aetiology. Since then, disruptions in several other genes have been implicated in CAS, with a substantial proportion of cases being explained by high-penetrance variants. However, the genetic architecture underlying other speech-related disorders remains less well understood. Thus, in the present study, we used systematic DNA sequencing methods to investigate idiopathic speech delay, as characterized by delayed speech development in the absence of a motor speech diagnosis (such as CAS), a language/reading disorder, or intellectual disability. We performed genome sequencing in a cohort of 23 children with a rigorous diagnosis of idiopathic speech delay. For roughly half of the sample (ten probands), sufficient DNA was also available for genome sequencing in both parents, allowing discovery of de novo variants. In the thirteen singleton probands, we focused on identifying loss-of-function and likely damaging missense variants in genes intolerant to such mutations. We found that one speech delay proband carried a pathogenic frameshift deletion in SETD1A, a gene previously implicated in a broader variable monogenic syndrome characterized by global developmental problems including delayed speech and/or language development, mild intellectual disability, facial dysmorphisms, and behavioural and psychiatric symptoms. Of note, pathogenic SETD1A variants have been independently reported in children with CAS in two separate studies. In other probands in our speech delay cohort, likely pathogenic missense variants were identified affecting highly conserved amino acids in key functional domains of SPTBN1 and ARF3. Overall, this study expands the phenotype spectrum associated with pathogenic SETD1A variants, to also include idiopathic speech delay without CAS or intellectual disability, and suggests additional novel potential candidate genes that may harbour high-penetrance variants that can disrupt speech development.

Using rare genetic mutations to revisit structural brain asymmetry

Asymmetry between the left and right hemisphere is a key feature of brain organization. Hemispheric functional specialization underlies some of the most advanced human-defining cognitive operations, such as articulated language, perspective taking, or rapid detection of facial cues. Yet, genetic investigations into brain asymmetry have mostly relied on common variants, which typically exert small effects on brain-related phenotypes. Here, we leverage rare genomic deletions and duplications to study how genetic alterations reverberate in human brain and behavior. We designed a pattern-learning approach to dissect the impact of eight high-effect-size copy number variations (CNVs) on brain asymmetry in a multi-site cohort of 552 CNV carriers and 290 non-carriers. Isolated multivariate brain asymmetry patterns spotlighted regions typically thought to subserve lateralized functions, including language, hearing, as well as visual, face and word recognition. Planum temporale asymmetry emerged as especially susceptible to deletions and duplications of specific gene sets. Targeted analysis of common variants through genome-wide association study (GWAS) consolidated partly diverging genetic influences on the right versus left planum temporale structure. In conclusion, our gene-brain-behavior data fusion highlights the consequences of genetically controlled brain lateralization on uniquely human cognitive capacities.

Quantifying transcriptome diversity: a review

Following the central dogma of molecular biology, gene expression heterogeneity can aid in predicting and explaining the wide variety of protein products, functions and, ultimately, heterogeneity in phenotypes. There is currently overlapping terminology used to describe the types of diversity in gene expression profiles, and overlooking these nuances can misrepresent important biological information. Here, we describe transcriptome diversity as a measure of the heterogeneity in (1) the expression of all genes within a sample or a single gene across samples in a population (gene-level diversity) or (2) the isoform-specific expression of a given gene (isoform-level diversity). We first overview modulators and quantification of transcriptome diversity at the gene level. Then, we discuss the role alternative splicing plays in driving transcript isoform-level diversity and how it can be quantified. Additionally, we overview computational resources for calculating gene-level and isoform-level diversity for high-throughput sequencing data. Finally, we discuss future applications of transcriptome diversity. This review provides a comprehensive overview of how gene expression diversity arises, and how measuring it determines a more complete picture of heterogeneity across proteins, cells, tissues, organisms and species.

Uncovering a multitude of stage-specific splice variants and putative protein isoforms generated along mouse spermatogenesis

BACKGROUND

Mammalian testis is a highly complex and heterogeneous tissue. This complexity, which mostly derives from spermatogenic cells, is reflected at the transcriptional level, with the largest number of tissue-specific genes and long noncoding RNAs (lncRNAs) compared to other tissues, and one of the highest rates of alternative splicing. Although it is known that adequate alternative-splicing patterns and stage-specific isoforms are critical for successful spermatogenesis, so far only a very limited number of reports have addressed a detailed study of alternative splicing and isoforms along the different spermatogenic stages.

RESULTS

In the present work, using highly purified stage-specific testicular cell populations, we detected 33,002 transcripts expressed throughout mouse spermatogenesis not annotated so far. These include both splice variants of already annotated genes, and of hitherto unannotated genes. Using conservative criteria, we uncovered 13,471 spermatogenic lncRNAs, which reflects the still incomplete annotation of lncRNAs. A distinctive feature of lncRNAs was their lower number of splice variants compared to protein-coding ones, adding to the conclusion that lncRNAs are, in general, less complex than mRNAs. Besides, we identified 2,794 unannotated transcripts with high coding potential (including some arising from yet unannotated genes), many of which encode unnoticed putative testis-specific proteins. Some of the most interesting coding splice variants were chosen, and validated through RT-PCR. Remarkably, the largest number of stage-specific unannotated transcripts are expressed during early meiotic prophase stages, whose study has been scarcely addressed in former transcriptomic analyses.

CONCLUSIONS

We detected a high number of yet unannotated genes and alternatively spliced transcripts along mouse spermatogenesis, hence showing that the transcriptomic diversity of the testis is considerably higher than previously reported. This is especially prominent for specific, underrepresented stages such as those of early meiotic prophase, and its unveiling may constitute a step towards the understanding of their key events.

Escape from X inactivation is directly modulated by levels of Xist non-coding RNA

In placental females, one copy of the two X chromosomes is largely silenced during a narrow developmental time window, in a process mediated by the non-coding RNA Xist1. Here, we demonstrate that Xist can initiate X-chromosome inactivation (XCI) well beyond early embryogenesis. By modifying its endogenous level, we show that Xist has the capacity to actively silence genes that escape XCI both in neuronal progenitor cells (NPCs) and in vivo, in mouse embryos. We also show that Xist plays a direct role in eliminating TAD-like structures associated with clusters of escapee genes on the inactive X chromosome, and that this is dependent on Xist's XCI initiation partner, SPEN2. We further demonstrate that Xist's function in suppressing gene expression of escapees and topological domain formation is reversible for up to seven days post-induction, but that sustained Xist up-regulation leads to progressively irreversible silencing and CpG island DNA methylation of facultative escapees. Thus, the distinctive transcriptional and regulatory topologies of the silenced X chromosome is actively, directly - and reversibly - controlled by Xist RNA throughout life.

Biological basis of extensive pleiotropy between blood traits and cancer risk

BACKGROUND

The immune system has a central role in preventing carcinogenesis. Alteration of systemic immune cell levels may increase cancer risk. However, the extent to which common genetic variation influences blood traits and cancer risk remains largely undetermined. Here, we identify pleiotropic variants and predict their underlying molecular and cellular alterations.

METHODS

Multivariate Cox regression was used to evaluate associations between blood traits and cancer diagnosis in cases in the UK Biobank. Shared genetic variants were identified from the summary statistics of the genome-wide association studies of 27 blood traits and 27 cancer types and subtypes, applying the conditional/conjunctional false-discovery rate approach. Analysis of genomic positions, expression quantitative trait loci, enhancers, regulatory marks, functionally defined gene sets, and bulk- and single-cell expression profiles predicted the biological impact of pleiotropic variants. Plasma small RNAs were sequenced to assess association with cancer diagnosis.

RESULTS

The study identified 4093 common genetic variants, involving 1248 gene loci, that contributed to blood-cancer pleiotropism. Genomic hotspots of pleiotropism include chromosomal regions 5p15-TERT and 6p21-HLA. Genes whose products are involved in regulating telomere length are found to be enriched in pleiotropic variants. Pleiotropic gene candidates are frequently linked to transcriptional programs that regulate hematopoiesis and define progenitor cell states of immune system development. Perturbation of the myeloid lineage is indicated by pleiotropic associations with defined master regulators and cell alterations. Eosinophil count is inversely associated with cancer risk. A high frequency of pleiotropic associations is also centered on the regulation of small noncoding Y-RNAs. Predicted pleiotropic Y-RNAs show specific regulatory marks and are overabundant in the normal tissue and blood of cancer patients. Analysis of plasma small RNAs in women who developed breast cancer indicates there is an overabundance of Y-RNA preceding neoplasm diagnosis.

CONCLUSIONS

This study reveals extensive pleiotropism between blood traits and cancer risk. Pleiotropism is linked to factors and processes involved in hematopoietic development and immune system function, including components of the major histocompatibility complexes, and regulators of telomere length and myeloid lineage. Deregulation of Y-RNAs is also associated with pleiotropism. Overexpression of these elements might indicate increased cancer risk.

Increased burden of rare protein-truncating variants in constrained, brain-specific and synaptic genes in extremely impulsively violent males with antisocial personality disorder

The genetic correlates of extreme impulsive violence are poorly understood, and there have been few studies that have characterized a large group of affected individuals both clinically and genetically. We performed whole exome sequencing (WES) in 290 males with the life-course-persistent, extremely impulsively violent form of antisocial personality disorder (APD) and analyzed the spectrum of rare protein-truncating variants (rPTVs). Comparisons were made with 314 male controls and publicly available genotype data. Functional annotation tools were used for biological interpretation. Participants were significantly more likely to harbor rPTVs in genes that are intolerant to loss-of-function variants (odds ratio [OR] 2.06; p < 0.001), specifically expressed in brain (OR 2.80; p = 0.036) and enriched for those involved in neurotransmitter transport and synaptic processes. In 60 individuals (20%), we identified rPTVs that we classified as clinically relevant based on their clinical associations, biological function and gene expression patterns. Of these, 37 individuals harbored rPTVs in 23 genes that are associated with a monogenic neurological disorder, and 23 individuals harbored rPTVs in 20 genes reportedly intolerant to loss-of-function variants. The analysis presents evidence in support of a model where presence of either one or several private, functionally relevant mutations contribute significantly to individual risk of life-course-persistent APD and reveals multiple individuals who could be affected by clinically unrecognized neuropsychiatric Mendelian disease. Thus, Mendelian diseases and increased rPTV burden may represent important factors for the development of extremely impulsive violent life-course-persistent forms of APD irrespective of their clinical presentation.

BORIS/CTCFL epigenetically reprograms clustered CTCF binding sites into alternative transcriptional start sites

BACKGROUND

Pervasive usage of alternative promoters leads to the deregulation of gene expression in carcinogenesis and may drive the emergence of new genes in spermatogenesis. However, little is known regarding the mechanisms underpinning the activation of alternative promoters.

RESULTS

Here we describe how alternative cancer-testis-specific transcription is activated. We show that intergenic and intronic CTCF binding sites, which are transcriptionally inert in normal somatic cells, could be epigenetically reprogrammed into active de novo promoters in germ and cancer cells. BORIS/CTCFL, the testis-specific paralog of the ubiquitously expressed CTCF, triggers the epigenetic reprogramming of CTCF sites into units of active transcription. BORIS binding initiates the recruitment of the chromatin remodeling factor, SRCAP, followed by the replacement of H2A histone with H2A.Z, resulting in a more relaxed chromatin state in the nucleosomes flanking the CTCF binding sites. The relaxation of chromatin around CTCF binding sites facilitates the recruitment of multiple additional transcription factors, thereby activating transcription from a given binding site. We demonstrate that the epigenetically reprogrammed CTCF binding sites can drive the expression of cancer-testis genes, long noncoding RNAs, retro-pseudogenes, and dormant transposable elements.

CONCLUSIONS

Thus, BORIS functions as a transcription factor that epigenetically reprograms clustered CTCF binding sites into transcriptional start sites, promoting transcription from alternative promoters in both germ cells and cancer cells.

Exploring Multi-Tissue Alternative Splicing and Skeletal Muscle Metabolism Regulation in Obese- and Lean-Type Pigs

Alternative splicing (AS) is a crucial mechanism in post-transcriptional regulation, contributing significantly to the diversity of the transcriptome and proteome. In this study, we performed a comprehensive AS profile in nine tissues obtained from Duroc (lean-type) and Luchuan (obese-type) pigs. Notably, 94,990 AS events from 14,393 genes were identified. Among these AS events, it was observed that 80% belonged to the skipped exon (SE) type. Functional enrichment analysis showed that genes with more than ten AS events were closely associated with tissue-specific functions. Additionally, the analysis of overlap between differentially alternative splicing genes (DSGs) and differentially expressed genes (DEGs) revealed the highest number of overlapped genes in the heart and skeletal muscle. The novelty of our study is that it identified and validated three genes (PYGM, MAPK11 and CAMK2B) in the glucagon signaling pathway, and their alternative splicing differences were highly significant across two pig breeds. In conclusion, our study offers novel insights into the molecular regulation of diverse tissue physiologies and the phenotypic differences between obese- and lean-type pigs, which are helpful for pig breeding.

Whole-Genome Sequencing Analysis Reveals New Susceptibility Loci and Structural Variants Associated with Progressive Supranuclear Palsy

Background

Progressive supranuclear palsy (PSP) is a rare neurodegenerative disease characterized by the accumulation of aggregated tau proteins in astrocytes, neurons, and oligodendrocytes. Previous genome-wide association studies for PSP were based on genotype array, therefore, were inadequate for the analysis of rare variants as well as larger mutations, such as small insertions/deletions (indels) and structural variants (SVs).

Method

In this study, we performed whole genome sequencing (WGS) and conducted association analysis for single nucleotide variants (SNVs), indels, and SVs, in a cohort of 1,718 cases and 2,944 controls of European ancestry. Of the 1,718 PSP individuals, 1,441 were autopsy-confirmed and 277 were clinically diagnosed.

Results

Our analysis of common SNVs and indels confirmed known genetic loci at MAPT, MOBP, STX6, SLCO1A2, DUSP10, and SP1, and further uncovered novel signals in APOE, FCHO1/MAP1S, KIF13A, TRIM24, TNXB, and ELOVL1. Notably, in contrast to Alzheimer's disease (AD), we observed the APOE ε2 allele to be the risk allele in PSP. Analysis of rare SNVs and indels identified significant association in ZNF592 and further gene network analysis identified a module of neuronal genes dysregulated in PSP. Moreover, seven common SVs associated with PSP were observed in the H1/H2 haplotype region (17q21.31) and other loci, including IGH, PCMT1, CYP2A13, and SMCP. In the H1/H2 haplotype region, there is a burden of rare deletions and duplications (P = 6.73×10-3) in PSP.

Conclusions

Through WGS, we significantly enhanced our understanding of the genetic basis of PSP, providing new targets for exploring disease mechanisms and therapeutic interventions.

Computational inference of eIF4F complex function and structure in human cancers

The canonical eukaryotic initiation factor 4F (eIF4F) complex, composed of eIF4G1, eIF4A1, and the cap-binding protein eIF4E, plays a crucial role in cap-dependent translation initiation in eukaryotic cells. An alternative cap-independent initiation can occur, involving only eIF4G1 and eIF4A1 through internal ribosome entry sites (IRESs). This mechanism is considered complementary to cap-dependent initiation, particularly in tumors under stress conditions. However, the selection and molecular mechanism of specific translation initiation remains poorly understood in human cancers. Thus, we analyzed gene copy number variations (CNVs) in TCGA tumor samples and found frequent amplification of genes involved in translation initiation. Copy number gains in EIF4G1 and EIF3E frequently co-occur across human cancers. Additionally, EIF4G1 expression strongly correlates with genes from cancer cell survival pathways including cell cycle and lipogenesis, in tumors with EIF4G1 amplification or duplication. Furthermore, we revealed that eIF4G1 and eIF4A1 protein levels strongly co-regulate with ribosomal subunits, eIF2, and eIF3 complexes, while eIF4E co-regulates with 4E-BP1, ubiquitination, and ESCRT proteins. Utilizing Alphafold predictions, we modeled the eIF4F structure with and without eIF4E binding. For cap-dependent initiation, our modeling reveals extensive interactions between the N-terminal eIF4E-binding domain of eIF4G1 and eIF4E. Furthermore, the eIF4G1 HEAT-2 domain positions eIF4E near the eIF4A1 N-terminal domain (NTD), resulting in the collaborative enclosure of the RNA binding cavity within eIF4A1. In contrast, during cap-independent initiation, the HEAT-2 domain directly binds the eIF4A1-NTD, leading to a stronger interaction between eIF4G1 and eIF4A1, thus closing the mRNA binding cavity without the involvement of eIF4E.

Should we really use graph neural networks for transcriptomic prediction?

The recent development of deep learning methods have undoubtedly led to great improvement in various machine learning tasks, especially in prediction tasks. This type of methods have also been adapted to answer various problems in bioinformatics, including automatic genome annotation, artificial genome generation or phenotype prediction. In particular, a specific type of deep learning method, called graph neural network (GNN) has repeatedly been reported as a good candidate to predict phenotypes from gene expression because its ability to embed information on gene regulation or co-expression through the use of a gene network. However, up to date, no complete and reproducible benchmark has ever been performed to analyze the trade-off between cost and benefit of this approach compared to more standard (and simpler) machine learning methods. In this article, we provide such a benchmark, based on clear and comparable policies to evaluate the different methods on several datasets. Our conclusion is that GNN rarely provides a real improvement in prediction performance, especially when compared to the computation effort required by the methods. Our findings on a limited but controlled simulated dataset shows that this could be explained by the limited quality or predictive power of the input biological gene network itself.

Complexity of human death: its physiological, transcriptomic, and microbiological implications

Human death is a complex, time-governed phenomenon that leads to the irreversible cessation of all bodily functions. Recent molecular and genetic studies have revealed remarkable experimental evidence of genetically programmed cellular death characterized by several physiological processes; however, the basic physiological function that occurs during the immediate postmortem period remains inadequately described. There is a paucity of knowledge connecting necrotic pathologies occurring in human organ tissues to complete functional loss of the human organism. Cells, tissues, organs, and organ systems show a range of differential resilience and endurance responses that occur during organismal death. Intriguingly, a persistent ambiguity in the study of postmortem physiological systems is the determination of the trajectory of a complex multicellular human body, far from life-sustaining homeostasis, following the gradual or sudden expiry of its regulatory systems. Recent groundbreaking investigations have resulted in a paradigm shift in understanding the cell biology and physiology of death. Two significant findings are that (i) most cells in the human body are microbial, and (ii) microbial cell abundance significantly increases after death. By addressing the physiological as well as the microbiological aspects of death, future investigations are poised to reveal innovative insights into the enigmatic biological activities associated with death and human decomposition. Understanding the elaborate crosstalk of abiotic and biotic factors in the context of death has implications for scientific discoveries important to informing translational knowledge regarding the transition from living to the non-living. There are important and practical needs for a transformative reestablishment of accepted models of biological death (i.e., artificial intelligence, AI) for more precise determinations of when the regulatory mechanisms for homeostasis of a living individual have ceased. In this review, we summarize mechanisms of physiological, genetic, and microbiological processes that define the biological changes and pathways associated with human organismal death and decomposition.

Long-read RNA sequencing identifies region- and sex-specific C57BL/6J mouse brain mRNA isoform expression and usage

Alternative splicing (AS) contributes to the biological heterogeneity between species, sexes, tissues, and cell types. Many diseases are either caused by alterations in AS or by alterations to AS. Therefore, measuring AS accurately and efficiently is critical for assessing molecular phenotypes, including those associated with disease. Long-read sequencing enables more accurate quantification of differentially spliced isoform expression than short-read sequencing approaches, and third-generation platforms facilitate high-throughput experiments. To assess differences in AS across the cerebellum, cortex, hippocampus, and striatum by sex, we generated and analyzed Oxford Nanopore Technologies (ONT) long-read RNA sequencing (lrRNA-Seq) C57BL/6J mouse brain cDNA libraries. From >85 million reads that passed quality control metrics, we calculated differential gene expression (DGE), differential transcript expression (DTE), and differential transcript usage (DTU) across brain regions and by sex. We found significant DGE, DTE, and DTU across brain regions and that the cerebellum had the most differences compared to the other three regions. Additionally, we found region-specific differential splicing between sexes, with the most sex differences in DTU in the cortex and no DTU in the hippocampus. We also report on two distinct patterns of sex DTU we observed, sex-divergent and sex-specific, that could potentially help explain sex differences in the prevalence and prognosis of various neurological and psychiatric disorders in future studies. Finally, we built a Shiny web application for researchers to explore the data further. Our study provides a resource for the community; it underscores the importance of AS in biological heterogeneity and the utility of long-read sequencing to better understand AS in the brain.

OXR1 maintains the retromer to delay brain aging under dietary restriction

Dietary restriction (DR) delays aging, but the mechanism remains unclear. We identified polymorphisms in mtd, the fly homolog of OXR1, which influenced lifespan and mtd expression in response to DR. Knockdown in adulthood inhibited DR-mediated lifespan extension in female flies. We found that mtd/OXR1 expression declines with age and it interacts with the retromer, which regulates trafficking of proteins and lipids. Loss of mtd/OXR1 destabilized the retromer, causing improper protein trafficking and endolysosomal defects. Overexpression of retromer genes or pharmacological restabilization with R55 rescued lifespan and neurodegeneration in mtd-deficient flies and endolysosomal defects in fibroblasts from patients with lethal loss-of-function of OXR1 variants. Multi-omic analyses in flies and humans showed that decreased Mtd/OXR1 is associated with aging and neurological diseases. mtd/OXR1 overexpression rescued age-related visual decline and tauopathy in a fly model. Hence, OXR1 plays a conserved role in preserving retromer function and is critical for neuronal health and longevity.

The human Y and inactive X chromosomes similarly modulate autosomal gene expression

Somatic cells of human males and females have 45 chromosomes in common, including the "active" X chromosome. In males the 46th chromosome is a Y; in females it is an "inactive" X (Xi). Through linear modeling of autosomal gene expression in cells from individuals with zero to three Xi and zero to four Y chromosomes, we found that Xi and Y impact autosomal expression broadly and with remarkably similar effects. Studying sex chromosome structural anomalies, promoters of Xi- and Y-responsive genes, and CRISPR inhibition, we traced part of this shared effect to homologous transcription factors-ZFX and ZFY-encoded by Chr X and Y. This demonstrates sex-shared mechanisms by which Xi and Y modulate autosomal expression. Combined with earlier analyses of sex-linked gene expression, our studies show that 21% of all genes expressed in lymphoblastoid cells or fibroblasts change expression significantly in response to Xi or Y chromosomes.

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

Identification of pan-cancer/testis genes and validation of therapeutic targeting in triple-negative breast cancer: Lin28a-based and Siglece-based vaccination induces antitumor immunity and inhibits metastasis

BACKGROUND

Cancer-testis (CT) genes are targets for tumor antigen-specific immunotherapy given that their expression is normally restricted to the immune-privileged testis in healthy individuals with aberrant expression in tumor tissues. While they represent targetable germ tissue antigens and play important functional roles in tumorigenesis, there is currently no standardized approach for identifying clinically relevant CT genes. Optimized algorithms and validated methods for accurate prediction of reliable CT antigens (CTAs) with high immunogenicity are also lacking.

METHODS

Sequencing data from the Genotype-Tissue Expression (GTEx) and The Genomic Data Commons (GDC) databases was used for the development of a bioinformatic pipeline to identify CT exclusive genes. A CT germness score was calculated based on the number of CT genes expressed within a tumor type and their degree of expression. The impact of tumor germness on clinical outcome was evaluated using healthy GTEx and GDC tumor samples. We then used a triple-negative breast cancer mouse model to develop and test an algorithm that predicts epitope immunogenicity based on the identification of germline sequences with strong major histocompatibility complex class I (MHCI) and MHCII binding affinities. Germline sequences for CT genes were synthesized as long synthetic peptide vaccines and tested in the 4T1 triple-negative model of invasive breast cancer with Poly(I:C) adjuvant. Vaccine immunogenicity was determined by flow cytometric analysis of in vitro and in vivo T-cell responses. Primary tumor growth and lung metastasis was evaluated by histopathology, flow cytometry and colony formation assay.

RESULTS

We developed a new bioinformatic pipeline to reliably identify CT exclusive genes as immunogenic targets for immunotherapy. We identified CT genes that are exclusively expressed within the testis, lack detectable thymic expression, and are significantly expressed in multiple tumor types. High tumor germness correlated with tumor progression but not with tumor mutation burden, supporting CTAs as appealing targets in low mutation burden tumors. Importantly, tumor germness also correlated with markers of antitumor immunity. Vaccination of 4T1 tumor-bearing mice with Siglece and Lin28a antigens resulted in increased T-cell antitumor immunity and reduced primary tumor growth and lung metastases.

CONCLUSION

Our results present a novel strategy for the identification of highly immunogenic CTAs for the development of targeted vaccines that induce antitumor immunity and inhibit metastasis.

Population-level annotation of lncRNAs in Arabidopsis reveals extensive expression variation associated with transposable element-like silencing

Long noncoding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing data sets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

Evaluation of noninvasive biospecimens for transcriptome studies

Transcriptome studies disentangle functional mechanisms of gene expression regulation and may elucidate the underlying biology of disease processes. However, the types of tissues currently collected typically assay a single post-mortem timepoint or are limited to investigating cell types found in blood. Noninvasive tissues may improve disease-relevant discovery by enabling more complex longitudinal study designs, by capturing different and potentially more applicable cell types, and by increasing sample sizes due to reduced collection costs and possible higher enrollment from vulnerable populations. Here, we develop methods for sampling noninvasive biospecimens, investigate their performance across commercial and in-house library preparations, characterize their biology, and assess the feasibility of using noninvasive tissues in a multitude of transcriptomic applications. We collected buccal swabs, hair follicles, saliva, and urine cell pellets from 19 individuals over three to four timepoints, for a total of 300 unique biological samples, which we then prepared with replicates across three library preparations, for a final tally of 472 transcriptomes. Of the four tissues we studied, we found hair follicles and urine cell pellets to be most promising due to the consistency of sample quality, the cell types and expression profiles we observed, and their performance in disease-relevant applications. This is the first study to thoroughly delineate biological and technical features of noninvasive samples and demonstrate their use in a wide array of transcriptomic and clinical analyses. We anticipate future use of these biospecimens will facilitate discovery and development of clinical applications.

Epitranscriptome analysis of NAD-capped RNA by spike-in-based normalization and prediction of chronological age

Nicotinamide adenine dinucleotide (NAD) can be used as an initiating nucleotide in RNA transcription to produce NAD-capped RNA (NAD-RNA). RNA modification by NAD that links metabolite with expressed transcript is a poorly studied epitranscriptomic modification. Current NAD-RNA profiling methods involve multi-steps of chemo-enzymatic labeling and affinity-based enrichment, thus presenting a critical analytical challenge to remove unwanted variations, particularly batch effects. Here, we propose a computational framework, enONE, to remove unwanted variations. We demonstrate that designed spike-in RNA, together with modular normalization procedures and evaluation metrics, can mitigate technical noise, empowering quantitative and comparative assessment of NAD-RNA across different datasets. Using enONE and a human aging cohort, we reveal age-associated features of NAD-capping and further develop an accurate RNA-based aging clock that combines signatures from both transcriptome and NAD-modified epitranscriptome. enONE facilitates the discovery of NAD-RNA responsive to physiological changes, laying an important foundation for functional investigations into this modification.

Transposon wave remodeled the epigenomic landscape in the rapid evolution of X-Chromosome dosage compensation

Sex chromosome dosage compensation is a model to understand the coordinated evolution of transcription; however, the advanced age of the sex chromosomes in model systems makes it difficult to study how the complex regulatory mechanisms underlying chromosome-wide dosage compensation can evolve. The sex chromosomes of Poecilia picta have undergone recent and rapid divergence, resulting in widespread gene loss on the male Y, coupled with complete X Chromosome dosage compensation, the first case reported in a fish. The recent de novo origin of dosage compensation presents a unique opportunity to understand the genetic and evolutionary basis of coordinated chromosomal gene regulation. By combining a new chromosome-level assembly of P. picta with whole-genome bisulfite sequencing and RNA-seq data, we determine that the YY1 transcription factor (YY1) DNA binding motif is associated with male-specific hypomethylated regions on the X, but not the autosomes. These YY1 motifs are the result of a recent and rapid repetitive element expansion on the P. picta X Chromosome, which is absent in closely related species that lack dosage compensation. Taken together, our results present compelling support that a disruptive wave of repetitive element insertions carrying YY1 motifs resulted in the remodeling of the X Chromosome epigenomic landscape and the rapid de novo origin of a dosage compensation system.

Structural basis for ligand recognition and signaling of the lysophosphatidylserine receptors GPR34 and GPR174

Lysophosphatidylserine (LysoPS) is a naturally occurring lipid mediator involved in various physiological and pathological processes especially those related to the immune system. GPR34, GPR174, and P2Y10 have been identified as the receptors for LysoPS, and its analogues have been developed as agonists or antagonists for these receptors. However, the lack of structural information hinders the drug development with novel characteristics, such as nonlipid ligands and allosteric modulators. Here, we determined the structures of human GPR34 and GPR174 in complex with LysoPS and G protein by cryo-EM. Combined with structural analysis and functional studies, we elucidated the lipid-binding modes of these receptors. By structural comparison, we identified the structural features of GPR34 and GPR174 in active state. Taken together, our findings provide insights into ligand recognition and signaling of LysoPS receptors and will facilitate the development of novel therapeutics for related inflammatory diseases and autoimmune diseases.

Isoform-level transcriptome-wide association uncovers genetic risk mechanisms for neuropsychiatric disorders in the human brain

Methods integrating genetics with transcriptomic reference panels prioritize risk genes and mechanisms at only a fraction of trait-associated genetic loci, due in part to an overreliance on total gene expression as a molecular outcome measure. This challenge is particularly relevant for the brain, in which extensive splicing generates multiple distinct transcript-isoforms per gene. Due to complex correlation structures, isoform-level modeling from cis-window variants requires methodological innovation. Here we introduce isoTWAS, a multivariate, stepwise framework integrating genetics, isoform-level expression and phenotypic associations. Compared to gene-level methods, isoTWAS improves both isoform and gene expression prediction, yielding more testable genes, and increased power for discovery of trait associations within genome-wide association study loci across 15 neuropsychiatric traits. We illustrate multiple isoTWAS associations undetectable at the gene-level, prioritizing isoforms of AKT3, CUL3 and HSPD1 in schizophrenia and PCLO with multiple disorders. Results highlight the importance of incorporating isoform-level resolution within integrative approaches to increase discovery of trait associations, especially for brain-relevant traits.

Insights into the roles of long noncoding RNAs in the communication between plants and the environment

In addition to coding proteins, RNA molecules, especially long noncoding RNAs (lncRNAs), have well-established functions in regulating gene expression. The number of studies focused on the roles played by different types of lncRNAs in a variety of plant biological processes has markedly increased. These lncRNA roles involve plant vegetative and reproductive growth and responses to biotic and abiotic stresses. In this review, we examine the classification, mechanisms, and functions of lncRNAs and then emphasize the roles played by these lncRNAs in the communication between plants and the environment mainly with respect to the following environmental factors: temperature, light, water, salt stress, and nutrient deficiencies. We also discuss the consensus among researchers and the remaining challenges and underscore the exciting ways lncRNAs may affect the biology of plants.

PTBP2 attenuation facilitates fibroblast to neuron conversion by promoting alternative splicing of neuronal genes

The direct conversion of human skin fibroblasts to neurons has a low efficiency and unclear mechanism. Here, we show that the knockdown of PTBP2 significantly enhanced the transdifferentiation induced by ASCL1, MIR9/9∗-124, and p53 shRNA (AMp) to generate mostly GABAergic neurons. Longitudinal RNA sequencing analyses identified the continuous induction of many RNA splicing regulators. Among these, the knockdown of RBFOX3 (NeuN), significantly abrogated the transdifferentiation. Overexpression of RBFOX3 significantly enhanced the conversion induced by AMp; the enhancement was occluded by PTBP2 knockdown. We found that PTBP2 attenuation significantly favored neuron-specific alternative splicing (AS) of many genes involved in synaptic transmission, signal transduction, and axon formation. RBFOX3 knockdown significantly reversed the effect, while RBFOX3 overexpression occluded the enhancement. The study reveals the critical role of neuron-specific AS in the direct conversion of human skin fibroblasts to neurons by showing that PTBP2 attenuation enhances this mechanism in concert with RBFOX3.

DNA damage in testicular germ cells and spermatozoa. When and how is it induced? How should we measure it? What does it mean?

This review surveys the causes and consequences of DNA damage in the male germ line from spermatogonial stem cells to fully differentiated spermatozoa. Within the stem cell population, DNA integrity is well maintained as a result of excellent DNA surveillance and repair; however, a progressive increase in background mutation rates does occur with paternal age possibly as a result of aberrant DNA repair as well as replication error. Once a germ cell has committed to spermatogenesis, it responds to genetic damage via a range of DNA repair pathways or, if this process fails, by the induction of apoptosis. When fully-differentiated spermatozoa are stressed, they also activate a truncated intrinsic apoptotic pathway which results in the activation of nucleases in the mitochondria and cytoplasm; however, the physical architecture of these cells prevents these enzymes from translocating to the nucleus to induce DNA fragmentation. Conversely, hydrogen peroxide released from the sperm midpiece during apoptosis is able to penetrate the nucleus and induce DNA damage. The base excision repair pathway responds to such damage by cleaving oxidized bases from the DNA, leaving abasic sites that are alkali-labile and readily detected with the comet assay. As levels of oxidative stress increase and these cells enter the perimortem, topoisomerase integrated into the sperm chromatin becomes activated by SUMOylation. Such activation may initially facilitate DNA repair by reannealing double strand breaks but ultimately prepares the DNA for destruction by nucleases released from the male reproductive tract. The abasic sites and oxidized base lesions found in live spermatozoa are mutagenic and may increase the mutational load carried by the offspring, particularly in the context of assisted conception. A variety of strategies are described for managing patients expressing high levels of DNA damage in their spermatozoa, to reduce the risks such lesions might pose to offspring health.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.