Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification

Single-cell profiling of lncRNA expression during Ebola virus infection in rhesus macaques

Long non-coding RNAs (lncRNAs) are involved in numerous biological processes and are pivotal mediators of the immune response, yet little is known about their properties at the single-cell level. Here, we generate a multi-tissue bulk RNAseq dataset from Ebola virus (EBOV) infected and not-infected rhesus macaques and identified 3979 novel lncRNAs. To profile lncRNA expression dynamics in immune circulating single-cells during EBOV infection, we design a metric, Upsilon, to estimate cell-type specificity. Our analysis reveals that lncRNAs are expressed in fewer cells than protein-coding genes, but they are not expressed at lower levels nor are they more cell-type specific when expressed in the same number of cells. In addition, we observe that lncRNAs exhibit similar changes in expression patterns to those of protein-coding genes during EBOV infection, and are often co-expressed with known immune regulators. A few lncRNAs change expression specifically upon EBOV entry in the cell. This study sheds light on the differential features of lncRNAs and protein-coding genes and paves the way for future single-cell lncRNA studies.

Noncoding transcripts are linked to brain resting-state activity in non-human primates

Brain-derived transcriptomes are known to correlate with resting-state brain activity in humans. Whether this association holds in nonhuman primates remains uncertain. Here, we search for such molecular correlates by integrating 757 transcriptomes derived from 100 macaque cortical regions with resting-state activity in separate conspecifics. We observe that 150 noncoding genes explain variations in resting-state activity at a comparable level with protein-coding genes. In-depth analysis of these noncoding genes reveals that they are connected to the function of nonneuronal cells such as oligodendrocytes. Co-expression network analysis finds that the modules of noncoding genes are linked to both autism and schizophrenia risk genes. Moreover, genes associated with resting-state noncoding genes are highly enriched in human resting-state functional genes and memory-effect genes, and their links with resting-state functional magnetic resonance imaging (fMRI) signals are altered in the brains of patients with autism. Our results highlight the potential for noncoding RNAs to explain resting-state activity in the nonhuman primate brain.

The origin and structural evolution of de novo genes in Drosophila

Although previously thought to be unlikely, recent studies have shown that de novo gene origination from previously non-genic sequences is a relatively common mechanism for gene innovation in many species and taxa. These young genes provide a unique set of candidates to study the structural and functional origination of proteins. However, our understanding of their protein structures and how these structures originate and evolve are still limited, due to a lack of systematic studies. Here, we combined high-quality base-level whole genome alignments, bioinformatic analysis, and computational structure modeling to study the origination, evolution, and protein structure of lineage-specific de novo genes. We identified 555 de novo gene candidates in D. melanogaster that originated within the Drosophilinae lineage. We found a gradual shift in sequence composition, evolutionary rates, and expression patterns with their gene ages, which indicates possible gradual shifts or adaptations of their functions. Surprisingly, we found little overall protein structural changes for de novo genes in the Drosophilinae lineage. Using Alphafold2, ESMFold, and molecular dynamics, we identified a number of de novo gene candidates with protein products that are potentially well-folded, many of which are more likely to contain transmembrane and signal proteins compared to other annotated protein-coding genes. Using ancestral sequence reconstruction, we found that most potentially well-folded proteins are often born folded. Interestingly, we observed one case where disordered ancestral proteins become ordered within a relatively short evolutionary time. Single-cell RNA-seq analysis in testis showed that although most de novo genes are enriched in spermatocytes, several young de novo genes are biased in the early spermatogenesis stage, indicating potentially important but less emphasized roles of early germline cells in the de novo gene origination in testis. This study provides a systematic overview of the origin, evolution, and structural changes of Drosophilinae-specific de novo genes.

Polycyclic aromatic hydrocarbons exposure and plasma lncRNA signature: A profile and functional analysis

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental pollutants that pose detrimental effects on human health, and the exploration of the associations of PAHs exposure with long non-coding RNA (lncRNA) may provide novel clues to the underlying mechanisms. In the present study, we detected 10 urinary PAHs metabolites by GC-MS and plasma lncRNAs levels by Human LncRNA Array v4 among 230 participants from two panels (160 in the Shiyan panel and 70 in the Wuhan-Zhuhai panel). We applied linear regression models to assess the associations between PAHs metabolites and lncRNAs separately in each panel and combined the results using fixed-effect meta-analysis. To explore the potential origin of PAHs-related lncRNAs in plasma, we estimated their tissue-specificity and associations between lncRNAs levels in plasma and leukocytes. Leukocytes mRNA sequencing data and RNA binding proteins were utilized to explore implicated pathways of identified lncRNAs. We found that urinary 1-hydroxyphenanthrene (1-OH-Phe) was inversely associated with 8 lncRNAs and positively associated with 1 lncRNA, as well as 9-hydroxyphenanthrene (9-OH-Phe) was inversely associated with 11 lncRNAs (FDR < 0.1). Tissue specificity analysis using Genome Tissue Expression database suggested that several identified lncRNAs might specifically express in organs targeted by PAHs exposure (lung, liver, heart, kidney, and brain). Besides, plasma levels of 1-OH-Phe related ENSG00000260616 and 9-OH-Phe related STARD4-AS1 were inversely associated with their intra-leukocytes levels (P value < 0.05). Notably, STARD4-AS1 was positively associated with the expression levels of its neighboring protein-coding gene (CAMK4 and STARD4) in leukocytes and were involved in pathways related to cellular response to DNA damage, which we further confirmed using DNA damage biomarker, 8-hydroxydeoxyguanosine. Functional analysis also revealed vital pathways related to cytokine-mediated signaling and glucose homeostasis. Our findings provided novel insights into plausible biological mechanisms underlying the adverse effects of PAHs exposure.

Tumour mutations in long noncoding RNAs enhance cell fitness

Long noncoding RNAs (lncRNAs) are linked to cancer via pathogenic changes in their expression levels. Yet, it remains unclear whether lncRNAs can also impact tumour cell fitness via function-altering somatic "driver" mutations. To search for such driver-lncRNAs, we here perform a genome-wide analysis of fitness-altering single nucleotide variants (SNVs) across a cohort of 2583 primary and 3527 metastatic tumours. The resulting 54 mutated and positively-selected lncRNAs are significantly enriched for previously-reported cancer genes and a range of clinical and genomic features. A number of these lncRNAs promote tumour cell proliferation when overexpressed in in vitro models. Our results also highlight a dense SNV hotspot in the widely-studied NEAT1 oncogene. To directly evaluate the functional significance of NEAT1 SNVs, we use in cellulo mutagenesis to introduce tumour-like mutations in the gene and observe a significant and reproducible increase in cell fitness, both in vitro and in a mouse model. Mechanistic studies reveal that SNVs remodel the NEAT1 ribonucleoprotein and boost subnuclear paraspeckles. In summary, this work demonstrates the utility of driver analysis for mapping cancer-promoting lncRNAs, and provides experimental evidence that somatic mutations can act through lncRNAs to enhance pathological cancer cell fitness.

Dynamic Spatial-temporal Expression Ratio of X Chromosome to Autosomes but Stable Dosage Compensation in Mammals

In the evolutionary model of dosage compensation, per-allele expression level of the X chromosome has been proposed to have twofold up-regulation to compensate its dose reduction in males (XY) compared to females (XX). However, the expression regulation of X-linked genes is still controversial, and comprehensive evaluations are still lacking. By integrating multi-omics datasets in mammals, we investigated the expression ratios including X to autosomes (X:AA ratio) and X to orthologs (X:XX ratio) at the transcriptome, translatome, and proteome levels. We revealed a dynamic spatial-temporal X:AA ratio during development in humans and mice. Meanwhile, by tracing the evolution of orthologous gene expression in chickens, platypuses, and opossums, we found a stable expression ratio of X-linked genes in humans to their autosomal orthologs in other species (X:XX ≈ 1) across tissues and developmental stages, demonstrating stable dosage compensation in mammals. We also found that different epigenetic regulations contributed to the high tissue specificity and stage specificity of X-linked gene expression, thus affecting X:AA ratios. It could be concluded that the dynamics of X:AA ratios were attributed to the different gene contents and expression preferences of the X chromosome, rather than the stable dosage compensation.

Weaker selection on genes with treatment-specific expression consistent with a limit on plasticity evolution in Arabidopsis thaliana

Differential gene expression between environments often underlies phenotypic plasticity. However, environment-specific expression patterns are hypothesized to relax selection on genes, and thus limit plasticity evolution. We collated over 27 terabases of RNA-sequencing data on Arabidopsis thaliana from over 300 peer-reviewed studies and 200 treatment conditions to investigate this hypothesis. Consistent with relaxed selection, genes with more treatment-specific expression have higher levels of nucleotide diversity and divergence at nonsynonymous sites but lack stronger signals of positive selection. This result persisted even after controlling for expression level, gene length, GC content, the tissue specificity of expression, and technical variation between studies. Overall, our investigation supports the existence of a hypothesized trade-off between the environment specificity of a gene's expression and the strength of selection on said gene in A. thaliana. Future studies should leverage multiple genome-scale datasets to tease apart the contributions of many variables in limiting plasticity evolution.

Long non-coding RNA and microRNA landscape of two major domesticated cotton species

Allotetraploid cotton plants Gossypium hirsutum and Gossypium barbadense have been widely cultivated for their natural, renewable textile fibres. Even though ncRNAs in domesticated cotton species have been extensively studied, systematic identification and annotation of lncRNAs and miRNAs expressed in various tissues and developmental stages under various biological contexts are limited. This influences the comprehension of their functions and future research on these cotton species. Here, we report high confidence lncRNAs and miRNA collection from G. hirsutum accession and G. barbadense accession using large-scale RNA-seq and small RNA-seq datasets incorporated into a user-friendly database, CoNCRAtlas. This database provides a wide range and depth of lncRNA and miRNA annotation based on the systematic integration of extensive annotations such as expression patterns derived from transcriptome data analysis in thousands of samples, as well as multi-omics annotations. We assume this comprehensive resource will accelerate evolutionary and functional studies in ncRNAs and inform future breeding programs for cotton improvement. CoNCRAtlas is accessible at http://www.nipgr.ac.in/CoNCRAtlas/.

The landscape of the long non-coding RNAs in developing mouse retinas

BACKGROUND

The long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes. Nevertheless, a global view of its expression and function in the mouse retina, a crucial model for neurogenesis study, still needs to be made available.

RESULTS

Herein, by integrating the established gene models and the result from ab initio prediction using short- and long-read sequencing, we characterized 4,523 lncRNA genes (MRLGs) in developing mouse retinas (from the embryonic day of 12.5 to the neonatal day of P28), which was so far the most comprehensive collection of retinal lncRNAs. Next, derived from transcriptomics analyses of different tissues and developing retinas, we found that the MRLGs were highly spatiotemporal specific in expression and played essential roles in regulating the genesis and function of mouse retinas. In addition, we investigated the expression of MRLGs in some mouse mutants and revealed that 97 intergenic MRLGs might be involved in regulating differentiation and development of retinal neurons through Math5, Isl1, Brn3b, NRL, Onecut1, or Onecut2 mediated pathways.

CONCLUSIONS

In summary, this work significantly enhanced our knowledge of lncRNA genes in mouse retina development and provided valuable clues for future exploration of their biological roles.

Evolution of High-Molecular-Mass Hyaluronic Acid is Associated with Subterranean Lifestyle

Hyaluronic acid (HA) is a major component of extracellular matrix (ECM) which plays an important role in development, cellular response to injury and inflammation, cell migration, and cancer. The naked mole-rat (NMR, Heterocephalus glaber ) contains abundant high-molecular-mass HA (HMM-HA) in its tissues, which contributes to this species' cancer resistance and possibly longevity. Here we report that abundant HMM-HA is found in a wide range of subterranean mammalian species, but not in phylogenetically related aboveground species. These species accumulate abundant HMM-HA by regulating the expression of genes involved in HA degradation and synthesis and contain unique mutations in these genes. The abundant high molecular weight HA may benefit the adaptation to subterranean environment by increasing skin elasticity and protecting from oxidative stress due to hypoxic subterranean environment. HMM-HA may also be coopted to confer cancer resistance and longevity to subterranean mammals. Our work suggests that HMM-HA has evolved with subterranean lifestyle.

Selection and geography shape male reproductive tract transcriptomes in Drosophila melanogaster

Transcriptome analysis of several animal clades suggests that male reproductive tract gene expression evolves quickly. However, the factors influencing the abundance and distribution of within-species variation, the ultimate source of interspecific divergence, are poorly known. Drosophila melanogaster, an ancestrally African species that has recently spread throughout the world and colonized the Americas in the last roughly 100 years, exhibits phenotypic and genetic latitudinal clines on multiple continents, consistent with a role for spatially varying selection in shaping its biology. Nevertheless, geographic expression variation in the Americas is poorly described, as is its relationship to African expression variation. Here, we investigate these issues through the analysis of two male reproductive tissue transcriptomes [testis and accessory gland (AG)] in samples from Maine (USA), Panama, and Zambia. We find dramatic differences between these tissues in differential expression between Maine and Panama, with the accessory glands exhibiting abundant expression differentiation and the testis exhibiting very little. Latitudinal expression differentiation appears to be influenced by the selection of Panama expression phenotypes. While the testis shows little latitudinal expression differentiation, it exhibits much greater differentiation than the accessory gland in Zambia vs American population comparisons. Expression differentiation for both tissues is non-randomly distributed across the genome on a chromosome arm scale. Interspecific expression divergence between D. melanogaster and D. simulans is discordant with rates of differentiation between D. melanogaster populations. Strongly heterogeneous expression differentiation across tissues and timescales suggests a complex evolutionary process involving major temporal changes in the way selection influences expression evolution in these organs.

An atlas of regulatory elements in chicken: A resource for chicken genetics and genomics

A comprehensive characterization of regulatory elements in the chicken genome across tissues will have substantial impacts on both fundamental and applied research. Here, we systematically identified and characterized regulatory elements in the chicken genome by integrating 377 genome-wide sequencing datasets from 23 adult tissues. In total, we annotated 1.57 million regulatory elements, representing 15 distinct chromatin states, and predicted about 1.2 million enhancer-gene pairs and 7662 super-enhancers. This functional annotation of the chicken genome should have wide utility on identifying regulatory elements accounting for gene regulation underlying domestication, selection, and complex trait regulation, which we explored. In short, this comprehensive atlas of regulatory elements provides the scientific community with a valuable resource for chicken genetics and genomics.

A pan-grass transcriptome reveals patterns of cellular divergence in crops

Different plant species within the grasses were parallel targets of domestication, giving rise to crops with distinct evolutionary histories and traits1. Key traits that distinguish these species are mediated by specialized cell types2. Here we compare the transcriptomes of root cells in three grass species-Zea mays, Sorghum bicolor and Setaria viridis. We show that single-cell and single-nucleus RNA sequencing provide complementary readouts of cell identity in dicots and monocots, warranting a combined analysis. Cell types were mapped across species to identify robust, orthologous marker genes. The comparative cellular analysis shows that the transcriptomes of some cell types diverged more rapidly than those of others-driven, in part, by recruitment of gene modules from other cell types. The data also show that a recent whole-genome duplication provides a rich source of new, highly localized gene expression domains that favour fast-evolving cell types. Together, the cell-by-cell comparative analysis shows how fine-scale cellular profiling can extract conserved modules from a pan transcriptome and provide insight on the evolution of cells that mediate key functions in crops.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

Splicing complexity as a pivotal feature of alternative exons in mammalian species

BACKGROUND

As a significant process of post-transcriptional gene expression regulation in eukaryotic cells, alternative splicing (AS) of exons greatly contributes to the complexity of the transcriptome and indirectly enriches the protein repertoires. A large number of studies have focused on the splicing inclusion of alternative exons and have revealed the roles of AS in organ development and maturation. Notably, AS takes place through a change in the relative abundance of the transcript isoforms produced by a single gene, meaning that exons can have complex splicing patterns. However, the commonly used percent spliced-in (Ψ) values only define the usage rate of exons, but lose information about the complexity of exons' linkage pattern. To date, the extent and functional consequence of splicing complexity of alternative exons in development and evolution is poorly understood.

RESULTS

By comparing splicing complexity of exons in six tissues (brain, cerebellum, heart, liver, kidney, and testis) from six mammalian species (human, chimpanzee, gorilla, macaque, mouse, opossum) and an outgroup species (chicken), we revealed that exons with high splicing complexity are prevalent in mammals and are closely related to features of genes. Using traditional machine learning and deep learning methods, we found that the splicing complexity of exons can be moderately predicted with features derived from exons, among which length of flanking exons and splicing strength of downstream/upstream splice sites are top predictors. Comparative analysis among human, chimpanzee, gorilla, macaque, and mouse revealed that, alternative exons tend to evolve to an increased level of splicing complexity and higher tissue specificity in splicing complexity. During organ development, not only developmentally regulated exons, but also 10-15% of non-developmentally regulated exons show dynamic splicing complexity.

CONCLUSIONS

Our analysis revealed that splicing complexity is an important metric to characterize the splicing dynamics of alternative exons during the development and evolution of mammals.

Convergent and complementary selection shaped gains and losses of eusociality in sweat bees

Sweat bees have repeatedly gained and lost eusociality, a transition from individual to group reproduction. Here we generate chromosome-length genome assemblies for 17 species and identify genomic signatures of evolutionary trade-offs associated with transitions between social and solitary living. Both young genes and regulatory regions show enrichment for these molecular patterns. We also identify loci that show evidence of complementary signals of positive and relaxed selection linked specifically to the convergent gains and losses of eusociality in sweat bees. This includes two pleiotropic proteins that bind and transport juvenile hormone (JH)-a key regulator of insect development and reproduction. We find that one of these proteins is primarily expressed in subperineurial glial cells that form the insect blood-brain barrier and that brain levels of JH vary by sociality. Our findings are consistent with a role of JH in modulating social behaviour and suggest that eusocial evolution was facilitated by alteration of the proteins that bind and transport JH, revealing how an ancestral developmental hormone may have been co-opted during one of life's major transitions. More broadly, our results highlight how evolutionary trade-offs have structured the molecular basis of eusociality in these bees and demonstrate how both directional selection and release from constraint can shape trait evolution.

A comparison of anatomic and cellular transcriptome structures across 40 human brain diseases

Genes associated with risk for brain disease exhibit characteristic expression patterns that reflect both anatomical and cell type relationships. Brain-wide transcriptomic patterns of disease risk genes provide a molecular-based signature, based on differential co-expression, that is often unique to that disease. Brain diseases can be compared and aggregated based on the similarity of their signatures which often associates diseases from diverse phenotypic classes. Analysis of 40 common human brain diseases identifies 5 major transcriptional patterns, representing tumor-related, neurodegenerative, psychiatric and substance abuse, and 2 mixed groups of diseases affecting basal ganglia and hypothalamus. Further, for diseases with enriched expression in cortex, single-nucleus data in the middle temporal gyrus (MTG) exhibits a cell type expression gradient separating neurodegenerative, psychiatric, and substance abuse diseases, with unique excitatory cell type expression differentiating psychiatric diseases. Through mapping of homologous cell types between mouse and human, most disease risk genes are found to act in common cell types, while having species-specific expression in those types and preserving similar phenotypic classification within species. These results describe structural and cellular transcriptomic relationships of disease risk genes in the adult brain and provide a molecular-based strategy for classifying and comparing diseases, potentially identifying novel disease relationships.

An artificial intelligence network-guided signature for predicting outcome and immunotherapy response in lung adenocarcinoma patients based on 26 machine learning algorithms

The immune cells play an increasingly vital role in influencing the proliferation, progression, and metastasis of lung adenocarcinoma (LUAD) cells. However, the potential of immune cells' specific genes-based model remains largely unknown. In the current study, by analysing single-cell RNA sequencing (scRNA-seq) data and bulk RNA sequencing data, the tumour-infiltrating immune cell (TIIC) associated signature was developed based on a total of 26 machine learning (ML) algorithms. As a result, the TIIC signature score could predict survival outcomes of LUAD patients across five independent datasets. The TIIC signature score showed superior performance to 168 previously established signatures in LUAD. Moreover, the TIIC signature score developed by the immunofluorescence staining of the tissue array of LUAD patients showed a prognostic value. Our research revealed a solid connection between TIIC signature score and tumour immunity as well as metabolism. Additionally, it has been discovered that the TIIC signature score can forecast genomic change, chemotherapeutic drug susceptibility, and-most significantly-immunotherapeutic response. As a newly demonstrated biomarker, the TIIC signature score facilitated the selection of the LUAD population who would benefit from future clinical stratification.

Lycium chinense Miller fruit extract lowers liver enzyme levels in subjects with mild hepatic dysfunction: a randomized, double-blind, placebo-controlled clinical trial

PURPOSE

In our previous study, we showed that Lycium chinense Miller fruit extract (LFE) exerted hepatoprotective effects in mice. In the current study, we examined the effect of LFE on liver enzyme levels in subjects with mild hepatic dysfunction.

METHODS

A total of 90 subjects, aged 19 to 70 years old, with abnormal alanine aminotransferase (ALT) levels, were randomly placed into either an LFE (n = 45) treatment group or a placebo group (n = 45). During the 12-week clinical trial, subjects in each group received either LFE or placebo capsules, and were instructed to take four tablets per day (1760 mg/day). The primary outcome of the study was the changes of ALT and γ-glutamyltransferase (GGT) levels in each subject. The safety of LFE supplementation was assessed and adverse events were recorded.

RESULTS

LFE supplementation for 12 weeks resulted in a significant reduction of ALT (P = 0.0498) and GGT (P = 0.0368) levels in comparison to the placebo. No clinically significant changes were observed in any safety parameters.

CONCLUSION

These results suggest that LFE can be applied to subjects with mild hepatic dysfunction with no possible side effects.

TRIAL REGISTRATION

This study was registered at the Clinical Research Information Service (CRIS) as no. KCT0003985.

SMRT and Illumina sequencing provide insights into mechanisms of lignin and terpenoids biosynthesis in Pinus massoniana Lamb

Wood and oleoresin are important industrial raw materials with high economic value; however, their molecular formation and biosynthesis mechanisms in different tissues of Pinus massoniana remain unexplored. Therefore, we used single-molecule real-time sequencing technology (SMRT) and Illumina RNA sequencing to establish a transcriptome dataset and explore the expression pattern of genes related to secondary metabolites involved in wood formation and oleoresin biosynthesis in six different P. massoniana tissues. In total, 63.58 Gb of polymerase reads were obtained, including 41,407 isoforms with an average length of 1822 bp. We identified 3939 and 8785 isoforms and 161 and 481 transcription factors with tissue expression specificity and in the reproductive and vegetative organs, respectively. Eighty isoforms were annotated as cellulose synthases and 224 isoforms involved in lignin biosynthesis were enriched. Additionally, we identified 217 isoforms involved in the terpenoid biosynthesis pathway, with needles having the most tissue-specific genes for terpenoid biosynthesis. Some isoforms related to lignin biosynthesis were highly expressed in the xylem, according to the results of transcriptome sequencing and real-time quantitative reverse-transcription polymerase chain reaction. Our research confirmed the advantages of SMRT sequencing and provided valuable information for the transcriptional annotation of P. massoniana, which will be beneficial for producing better raw wood and oleoresin materials.

Crosstalk between SOX Genes and Long Non-Coding RNAs in Glioblastoma

Glioblastoma (GBM) continues to be the most devastating primary brain malignancy. Despite significant advancements in understanding basic GBM biology and enormous efforts in developing new therapeutic approaches, the prognosis for most GBM patients remains poor with a median survival time of 15 months. Recently, the interplay between the SOX (SRY-related HMG-box) genes and lncRNAs (long non-coding RNAs) has become the focus of GBM research. Both classes of molecules have an aberrant expression in GBM and play essential roles in tumor initiation, progression, therapy resistance, and recurrence. In GBM, SOX and lncRNAs crosstalk through numerous functional axes, some of which are part of the complex transcriptional and epigenetic regulatory mechanisms. This review provides a systematic summary of current literature data on the complex interplay between SOX genes and lncRNAs and represents an effort to underscore the effects of SOX/lncRNA crosstalk on the malignant properties of GBM cells. Furthermore, we highlight the significance of this crosstalk in searching for new biomarkers and therapeutic approaches in GBM treatment.

A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain

The mammalian brain is composed of millions to billions of cells that are organized into numerous cell types with specific spatial distribution patterns and structural and functional properties. An essential step towards understanding brain function is to obtain a parts list, i.e., a catalog of cell types, of the brain. Here, we report a comprehensive and high-resolution transcriptomic and spatial cell type atlas for the whole adult mouse brain. The cell type atlas was created based on the combination of two single-cell-level, whole-brain-scale datasets: a single-cell RNA-sequencing (scRNA-seq) dataset of ~7 million cells profiled, and a spatially resolved transcriptomic dataset of ~4.3 million cells using MERFISH. The atlas is hierarchically organized into five nested levels of classification: 7 divisions, 32 classes, 306 subclasses, 1,045 supertypes and 5,200 clusters. We systematically analyzed the neuronal, non-neuronal, and immature neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell type organization in different brain regions, in particular, a dichotomy between the dorsal and ventral parts of the brain: the dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. We also systematically characterized cell-type specific expression of neurotransmitters, neuropeptides, and transcription factors. The study uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types across the brain, suggesting they mediate a myriad of modes of intercellular communications. Finally, we found that transcription factors are major determinants of cell type classification in the adult mouse brain and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole-mouse-brain transcriptomic and spatial cell type atlas establishes a benchmark reference atlas and a foundational resource for deep and integrative investigations of cell type and circuit function, development, and evolution of the mammalian brain.

The nearly complete assembly of the Cercis chinensis genome and Fabaceae phylogenomic studies provide insights into new gene evolution

Fabaceae is a large family of angiosperms with high biodiversity that contains a variety of economically important crops and model plants for the study of biological nitrogen fixation. Polyploidization events have been extensively studied in some Fabaceae plants, but the occurrence of new genes is still concealed, owing to a lack of genomic information on certain species of the basal clade of Fabaceae. Cercis chinensis (Cercidoideae) is one such species; it diverged earliest from Fabaceae and is essential for phylogenomic studies and new gene predictions in Fabaceae. To facilitate genomic studies on Fabaceae, we performed genome sequencing of C. chinensis and obtained a 352.84 Mb genome, which was further assembled into seven pseudochromosomes with 30 612 predicted protein-coding genes. Compared with other legume genomes, that of C. chinensis exhibits no lineage-specific polyploidization event. Further phylogenomic analyses of 22 legumes and 11 other angiosperms revealed that many gene families are lineage specific before and after the diversification of Fabaceae. Among them, dozens of genes are candidates for new genes that have evolved from intergenic regions and are thus regarded as de novo-originated genes. They differ significantly from established genes in coding sequence length, exon number, guanine-cytosine content, and expression patterns among tissues. Functional analysis revealed that many new genes are related to asparagine metabolism. This study represents an important advance in understanding the evolutionary pattern of new genes in legumes and provides a valuable resource for plant phylogenomic studies.

Age-dependent expression of ion channel genes in rat

Ion channels regulate a large number of cellular functions and their functional role in many diseases makes them potential therapeutic targets. Given their diverse distribution across multiple organs, the roles of ion channels, particularly in age-associated transcriptomic changes in specific organs, are yet to be fully revealed. Using RNA-seq data, we investigated the rat transcriptomic profiles of ion channel genes across 11 organs/tissues and 4 developmental stages in both sexes of Fischer 344 rats and identify tissue-specific and age-dependent changes in ion channel gene expression. Organ-enriched ion channel genes were identified. In particular, the brain showed higher tissue-specificity of ion channel genes, including Gabrd, Gabra6, Gabrg2, Grin2a, and Grin2b. Notably, age-dependent changes in ion channel gene expression were prominently observed in the thymus, including in Aqp1, Clcn4, Hvcn1, Itpr1, Kcng2, Kcnj11, Kcnn3, and Trpm2. Our comprehensive study of ion channel gene expression will serve as a primary resource for biological studies of aging-related diseases caused by abnormal ion channel functions.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

Detecting macroevolutionary genotype-phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype-phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω_C, which measures the error-corrected convergence rate of protein evolution. While ω_C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype-phenotype associations, even in lineages that diverged for hundreds of millions of years.

Independent Evolution of Sex Chromosomes and Male Pregnancy-Related Genes in Two Seahorse Species

Unlike birds and mammals, many teleosts have homomorphic sex chromosomes, and changes in the chromosome carrying the sex-determining locus, termed "turnovers", are common. Recent turnovers allow studies of several interesting questions. One question is whether the new sex-determining regions evolve to become completely non-recombining, and if so, how and why. Another is whether (as predicted) evolutionary changes that benefit one sex accumulate in the newly sex-linked region. To study these questions, we analyzed the genome sequences of two seahorse species of the Syngnathidae, a fish group in which many species evolved a unique structure, the male brood pouch. We find that both seahorse species have XY sex chromosome systems, but their sex chromosome pairs are not homologs, implying that at least one turnover event has occurred. The Y-linked regions occupy 63.9% and 95.1% of the entire sex chromosome of the two species and do not exhibit extensive sequence divergence with their X-linked homologs. We find evidence for occasional recombination between the extant sex chromosomes that may account for their homomorphism. We argue that these Y-linked regions did not evolve by recombination suppression after the turnover, but by the ancestral nature of the low crossover rates in these chromosome regions. With such an ancestral crossover landscape, a turnover can instantly create an extensive Y-linked region. Finally, we test for adaptive evolution of male pouch-related genes after they became Y-linked in the seahorse.

Linolenic Acid-Derived Oxylipins Inhibit Aflatoxin Biosynthesis in Aspergillus flavus through Activation of Imizoquin Biosynthesis

Oxylipins play important signaling roles in aflatoxin (AF) biosynthesis in Aspergillus flavus. We previously showed that exogenous supply of autoxidated linolenic acid (AL) inhibited AF biosynthesis in A. flavus via oxylipins, but the molecular mechanism is still unknown. Here, we performed multiomics analyses of A. flavus grown in media with or without AL. Targeted metabolite analyses and quantitative reverse transcription (qRT)-polymerase chain reaction (PCR) showed that the imizoquin (IMQ) biosynthetic pathway was distinctly upregulated in the presence of AL. 13C-glucose labeling confirmed in parallel that the tricarboxylic acid cycle was also enhanced by AL, consistent with observed increases in mycelial growth. Moreover, we integrated thermal proteome profiling and molecular dynamics simulations to identify a potential receptor of AL; AL was found to interact with a transporter (ImqJ) located in the IMQ gene cluster, primarily through hydrophobic interactions. Further analyses of strains with an IMQ pathway transcription factor overexpressed or knocked out confirmed that this pathway was critical for AL-mediated inhibition of AF biosynthesis. Comparison of 22 assembled A. flavus and Aspergillus oryzae genomes showed that genes involved in the IMQ pathway were positively selected in A. oryzae. Taken together, the results of our study provide novel insights into oxylipin-mediated regulation of AF biosynthesis and suggest potential methods for preventing AF contamination of crops.

Diversity in Expression Biases of Lineage-Specific Genes During Development and Anhydrobiosis Among Tardigrade Species

Lineage-specific genes can contribute to the emergence and evolution of novel traits and adaptations. Tardigrades are animals that have adapted to tolerate extreme conditions by undergoing a form of cryptobiosis called anhydrobiosis, a physical transformation to an inactive desiccated state. While studies to understand the genetics underlying the interspecies diversity in anhydrobiotic transitions have identified tardigrade-specific genes and family expansions involved in this process, the contributions of species-specific genes to the variation in tardigrade development and cryptobiosis are less clear. We used previously published transcriptomes throughout development and anhydrobiosis (5 embryonic stages, 7 juvenile stages, active adults, and tun adults) to assess the transcriptional biases of different classes of genes between 2 tardigrade species, Hypsibius exemplaris and Ramazzottius varieornatus. We also used the transcriptomes of 2 other tardigrades, Echiniscoides sigismundi and Richtersius coronifer, and data from 3 non-tardigrade species (Adenita vaga, Drosophila melanogaster, and Caenorhabditis elegans) to help identify lineage-specific genes. We found that lineage-specific genes have generally low and narrow expression but are enriched among biased genes in different stages of development depending on the species. Biased genes tend to be specific to early and late development, but there is little overlap in functional enrichment of biased genes between species. Gene expansions in the 2 tardigrades also involve families with different functions despite homologous genes being expressed during anhydrobiosis in both species. Our results demonstrate the interspecific variation in transcriptional contributions and biases of lineage-specific genes during development and anhydrobiosis in 2 tardigrades.

Genome-wide identification and comprehensive analysis reveal potential roles of long non-coding RNAs in fruit development of southern highbush blueberry (Vaccinium corymbosum L.)

Introduction

Blueberries have a high antioxidant content and are produced as healthy food worldwide. Long non-coding RNAs (lncRNAs) are a type of regulatory RNAs that play a variety of roles in plants. Nonetheless, information on lncRNAs and their functions during blueberry fruit development is scarce in public databases.

Methods

In the present study, we performed genome-wide identification of lncRNAs in a southern highbush blueberry using strand-specific RNA sequencing (ssRNA-Seq). Differentially expressed lncRNAs (DE-lncRNAs) and their potential target genes were analyzed at four stages of fruit development. Cis-regulatory DE-lncRNAs were predicted using co-localization analysis.

Results

These findings included a total of 25,036 lncRNAs from 17,801 loci. Blueberry lncRNAs had shorter transcript lengths, smaller open reading frame (ORF) sizes, fewer exons, and fewer isoforms than protein-coding RNAs, as well as lower expression levels and higher stage-specificity during fruit development. A total of 105 DE-lncRNAs were identified among the comparison group of PAD vs. CUP, 443 DE-lncRNAs were detected when comparing CUP with PINK fruits, and 285 DE-lncRNAs were revealed when comparing PINK and BLUE fruits. According to Kyoto Encyclopedia of Genes and Genomes annotation, target genes of DE-lncRNAs were primarily enriched in the "Autophagy-other", "DNA replication", "Endocytosis", 'photosynthesis' and 'chlorophyll metabolism' pathways, suggesting that lncRNAs may pay potential roles in fruit expansion and ripening. Moreover, several lncRNAs have been proposed as cis-regulators of the key genes involved in flavonoid biosynthesis. MSTRG.107242.6, and its putative target gene, BTB/POZ and TAZ domain-containing protein, might play critical roles in anthocyanin accumulation in blueberries.

Discussion

These findings highlight the regulatory function of lncRNAs and aid in elucidating the molecular mechanism underlying blueberry fruit growth.

Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

Robust and rigorous identification of tissue-specific genes by statistically extending tau score

OBJECTIVES

In this study, we aimed to identify tissue-specific genes for various human tissues/organs more robustly and rigorously by extending the tau score algorithm.

INTRODUCTION

Tissue-specific genes are a class of genes whose functions and expressions are preferred in one or several tissues restrictedly. Identification of tissue-specific genes is essential for discovering multi-cellular biological processes such as tissue-specific molecular regulations, tissue development, physiology, and the pathogenesis of tissue-associated diseases.

MATERIALS AND METHODS

Gene expression data derived from five large RNA sequencing (RNA-seq) projects, spanning 96 different human tissues, were retrieved from ArrayExpress and ExpressionAtlas. The first step is categorizing genes using significant filters and tau score as a specificity index. After calculating tau for each gene in all datasets separately, statistical distance from the maximum expression level was estimated using a new meaningful procedure. Specific expression of a gene in one or several tissues was calculated after the integration of tau and statistical distance estimation, which is called as extended tau approach. Obtained tissue-specific genes for 96 different human tissues were functionally annotated, and some comparisons were carried out to show the effectiveness of the extended tau method.

RESULTS AND DISCUSSION

Categorization of genes based on expression level and identification of tissue-specific genes for a large number of tissues/organs were executed. Genes were successfully assigned to multiple tissues by generating the extended tau approach as opposed to the original tau score, which can assign tissue specificity to single tissue only.

Pan-cancer surveys indicate cell cycle-related roles of primate-specific genes in tumors and embryonic cerebrum

BACKGROUND

Despite having been extensively studied, it remains largely unclear why humans bear a particularly high risk of cancer. The antagonistic pleiotropy hypothesis predicts that primate-specific genes (PSGs) tend to promote tumorigenesis, while the molecular atavism hypothesis predicts that PSGs involved in tumors may represent recently derived duplicates of unicellular genes. However, these predictions have not been tested.

RESULTS

By taking advantage of pan-cancer genomic data, we find the upregulation of PSGs across 13 cancer types, which is facilitated by copy-number gain and promoter hypomethylation. Meta-analyses indicate that upregulated PSGs (uPSGs) tend to promote tumorigenesis and to play cell cycle-related roles. The cell cycle-related uPSGs predominantly represent derived duplicates of unicellular genes. We prioritize 15 uPSGs and perform an in-depth analysis of one unicellular gene-derived duplicate involved in the cell cycle, DDX11. Genome-wide screening data and knockdown experiments demonstrate that DDX11 is broadly essential across cancer cell lines. Importantly, non-neutral amino acid substitution patterns and increased expression indicate that DDX11 has been under positive selection. Finally, we find that cell cycle-related uPSGs are also preferentially upregulated in the highly proliferative embryonic cerebrum.

CONCLUSIONS

Consistent with the predictions of the atavism and antagonistic pleiotropy hypotheses, primate-specific genes, especially those PSGs derived from cell cycle-related genes that emerged in unicellular ancestors, contribute to the early proliferation of the human cerebrum at the cost of hitchhiking by similarly highly proliferative cancer cells.

GeneFriends: gene co-expression databases and tools for humans and model organisms

Gene co-expression analysis has emerged as a powerful method to provide insights into gene function and regulation. The rapid growth of publicly available RNA-sequencing (RNA-seq) data has created opportunities for researchers to employ this abundant data to help decipher the complexity and biology of genomes. Co-expression networks have proven effective for inferring the relationship between the genes, for gene prioritization and for assigning function to poorly annotated genes based on their co-expressed partners. To facilitate such analyses we created previously an online co-expression tool for humans and mice entitled GeneFriends. To continue providing a valuable tool to the scientific community, we have now updated the GeneFriends database and website. Here, we present the new version of GeneFriends, which includes gene and transcript co-expression networks based on RNA-seq data from 46 475 human and 34 322 mouse samples. The new database also encompasses tissue-specific gene co-expression networks for 20 human and 21 mouse tissues, dataset-specific gene co-expression maps based on TCGA and GTEx projects and gene co-expression networks for additional seven model organisms (fruit fly, zebrafish, worm, rat, yeast, cow and chicken). GeneFriends is freely available at http://www.genefriends.org/.

Comparative transcriptomics reveals commonalities and differences in the genetic underpinnings of a floral dimorphism

Distyly, a floral dimorphism associated with heteromorphic self-incompatibility and controlled by the S-locus supergene, evolved independently multiple times. Comparative analyses of the first transcriptome atlas for the main distyly model, Primula veris, with other distylous species produced the following findings. A set of 53 constitutively expressed genes in P. veris did not include any of the housekeeping genes commonly used to normalize gene expression in qPCR experiments. The S-locus gene CYP^T acquired its role in controlling style elongation via a change in expression profile. Comparison of genes differentially expressed between floral morphs revealed that brassinosteroids and auxin are the main hormones controlling style elongation in P. veris and Fagopyrum esculentum, respectively. Furthermore, shared biochemical pathways might underlie the expression of distyly in the distantly related P. veris, F. esculentum and Turnera subulata, suggesting a degree of correspondence between evolutionary convergence at phenotypic and molecular levels. Finally, we provide the first evidence supporting the previously proposed hypothesis that distyly supergenes of distantly related species evolved via the recruitment of genes related to the phytochrome-interacting factor (PIF) signaling network. To conclude, this is the first study that discovered homologous genes involved in the control of distyly in distantly related taxa.

Machine learning-based tumor-infiltrating immune cell-associated lncRNAs for predicting prognosis and immunotherapy response in patients with glioblastoma

Long noncoding ribonucleic acids (RNAs; lncRNAs) have been associated with cancer immunity regulation. However, the roles of immune cell-specific lncRNAs in glioblastoma (GBM) remain largely unknown. In this study, a novel computational framework was constructed to screen the tumor-infiltrating immune cell-associated lncRNAs (TIIClnc) for developing TIIClnc signature by integratively analyzing the transcriptome data of purified immune cells, GBM cell lines and bulk GBM tissues using six machine learning algorithms. As a result, TIIClnc signature could distinguish survival outcomes of GBM patients across four independent datasets, including the Xiangya in-house dataset, and more importantly, showed superior performance than 95 previously established signatures in gliomas. TIIClnc signature was revealed to be an indicator of the infiltration level of immune cells and predicted the response outcomes of immunotherapy. The positive correlation between TIIClnc signature and CD8, PD-1 and PD-L1 was verified in the Xiangya in-house dataset. As a newly demonstrated predictive biomarker, the TIIClnc signature enabled a more precise selection of the GBM population who would benefit from immunotherapy and should be validated and applied in the near future.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications

BACKGROUND

Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates.

RESULTS

We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements.

CONCLUSIONS

In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

Systematic analysis and prediction of genes associated with monogenic disorders on human chromosome X

Disease gene discovery on chromosome (chr) X is challenging owing to its unique modes of inheritance. We undertook a systematic analysis of human chrX genes. We observe a higher proportion of disorder-associated genes and an enrichment of genes involved in cognition, language, and seizures on chrX compared to autosomes. We analyze gene constraints, exon and promoter conservation, expression, and paralogues, and report 127 genes sharing one or more attributes with known chrX disorder genes. Using machine learning classifiers trained to distinguish disease-associated from dispensable genes, we classify 247 genes, including 115 of the 127, as having high probability of being disease-associated. We provide evidence of an excess of variants in predicted genes in existing databases. Finally, we report damaging variants in CDK16 and TRPC5 in patients with intellectual disability or autism spectrum disorders. This study predicts large-scale gene-disease associations that could be used for prioritization of X-linked pathogenic variants.

TMEM158 expression is negatively regulated by AR signaling and associated with favorite survival outcomes in prostate cancers

BACKGROUND

Membrane protein TMEM158 was initially reported as a Ras-induced gene during senescence and has been implicated as either an oncogenic factor or tumor suppressor, depending on tumor types. It is unknown if TMEM158 expression is altered in prostate cancers.

METHODS

Multiple public gene expression datasets from RNA-seq and cDNA microarray assays were utilized to analyze candidate gene expression profiles. TMEM158 protein expression was assessed using an immunohistochemistry approach on a tissue section array from benign and malignant prostate tissues. Comparisons of gene expression profiles were conducted using the bioinformatics software R package.

RESULTS

COX regression-based screening identified the membrane protein TMEM158 gene as negatively associated with disease-specific and progression-free survival in prostate cancer patients. Gene expression at the mRNA and protein levels revealed that TMEM158 expression was significantly reduced in malignant tissues compared to benign compartments. Meanwhile, TMEM158 downregulation was strongly correlated with advanced clinicopathological features, including late-stage diseases, lymph node invasion, higher PSA levels, residual tumors after surgery, and adverse Gleason scores. In castration-resistant prostate cancers, TMEM158 expression was negatively correlated with AR signaling activity but positively correlated with neuroendocrinal progression index. Consistently, in cell culture models, androgen treatment reduced TMEM158 expression, while androgen deprivation led to upregulation of TMEM158 expression. Correlation analysis showed a tight correlation of TMEM158 expression with the level of R-Ras gene expression, which was also significantly downregulated in prostate cancers. Tumor immune infiltration profiling analysis discovered a strong association of TMEM158 expression with NK cell and Mast cell enrichment.

CONCLUSION

The membrane protein TMEM158 is significantly downregulated in prostate cancer and is tightly associated with disease progression, anti-tumor immune infiltration, and patient survival outcome.

A network view of human immune system and virus-human interaction

The immune system is highly networked and complex, which is continuously changing as encountering old and new pathogens. However, reductionism-based researches do not give a systematic understanding of the molecular mechanism of the immune response and viral pathogenesis. Here, we present HUMPPI-2022, a high-quality human protein-protein interaction (PPI) network, containing > 11,000 protein-coding genes with > 78,000 interactions. The network topology and functional characteristics analyses of the immune-related genes (IRGs) reveal that IRGs are mostly located in the center of the network and link genes of diverse biological processes, which may reflect the gene pleiotropy phenomenon. Moreover, the virus-human interactions reveal that pan-viral targets are mostly hubs, located in the center of the network and enriched in fundamental biological processes, but not for coronavirus. Finally, gene age effect was analyzed from the view of the host network for IRGs and virally-targeted genes (VTGs) during evolution, with IRGs gradually became hubs and integrated into host network through bridging functionally differentiated modules. Briefly, HUMPPI-2022 serves as a valuable resource for gaining a better understanding of the composition and evolution of human immune system, as well as the pathogenesis of viruses.

RiceNCexp: a rice non-coding RNA co-expression atlas based on massive RNA-seq and small-RNA seq data

Non-coding RNAs (ncRNAs) play important roles in regulating expression of protein-coding genes. Although gene expression databases have emerged in a timely manner, a comprehensive expression database for ncRNAs is still lacking. Herein, we constructed a rice ncRNA co-expression atlas (RiceNCexp), based on 491 RNA-seq and 274 small RNA (sRNA)-seq datasets. RiceNCexp hosts four types of ncRNAs, namely lncRNAs, PHAS genes, miRNAs, and phasiRNAs. RiceNCexp provides comprehensive expression information for rice ncRNAs in 22 tissues/organs, an efficient tau-based mining tool for tissue-specific ncRNAs, and the robust co-expression analysis among ncRNAs or between ncRNAs and protein-coding genes, based on 116 pairs of RNA-seq and sRNA-seq libraries from the same experiments. In summary, RiceNCexp is a user-friendly and comprehensive rice ncRNA co-expression atlas and can be freely accessed at https://cbi.njau.edu.cn/RiceNCexp/.

isomiRdb: microRNA expression at isoform resolution

A significant fraction of mature miRNA transcripts carries sequence and/or length variations, termed isomiRs. IsomiRs are differentially abundant in cell types, tissues, body fluids or patients' samples. Not surprisingly, multiple studies describe a physiological and pathophysiological role. Despite their importance, systematically collected and annotated isomiR information available in databases remains limited. We thus developed isomiRdb, a comprehensive resource that compiles miRNA expression data at isomiR resolution from various sources. We processed 42 499 human miRNA-seq datasets (5.9 × 1011 sequencing reads) and consistently analyzed them using miRMaster and sRNAbench. Our database provides online access to the 90 483 most abundant isomiRs (>1 RPM in at least 1% of the samples) from 52 tissues and 188 cell types. Additionally, the full set of over 3 million detected isomiRs is available for download. Our resource can be queried at the sample, miRNA or isomiR level so users can quickly answer common questions about the presence/absence of a particular miRNA/isomiR in tissues of interest. Further, the database facilitates to identify whether a potentially interesting new isoform has been detected before and its frequency. In addition to expression tables, isomiRdb can generate multiple interactive visualisations including violin plots and heatmaps. isomiRdb is free to use and publicly available at: https://www.ccb.uni-saarland.de/isomirdb.

Transcriptional signatures of wheat inflorescence development

In order to maintain global food security, it will be necessary to increase yields of the cereal crops that provide most of the calories and protein for the world's population, which includes common wheat (Triticum aestivum L.). An important wheat yield component is the number of grain-holding spikelets which form on the spike during inflorescence development. Characterizing the gene regulatory networks controlling the timing and rate of inflorescence development will facilitate the selection of natural and induced gene variants that contribute to increased spikelet number and yield. In the current study, co-expression and gene regulatory networks were assembled from a temporal wheat spike transcriptome dataset, revealing the dynamic expression profiles associated with the progression from vegetative meristem to terminal spikelet formation. Consensus co-expression networks revealed enrichment of several transcription factor families at specific developmental stages including the sequential activation of different classes of MIKC-MADS box genes. This gene regulatory network highlighted interactions among a small number of regulatory hub genes active during terminal spikelet formation. Finally, the CLAVATA and WUSCHEL gene families were investigated, revealing potential roles for TtCLE13, TtWOX2, and TtWOX7 in wheat meristem development. The hypotheses generated from these datasets and networks further our understanding of wheat inflorescence development.

Phylogenetic analysis and expression profiles of jasmonate ZIM-domain gene family provide insight into abiotic stress resistance in sunflower

Jasmonate ZIM-domain (JAZ) proteins act as inhibitory factors of the jasmonic acid (JA) pathway, which is involved in regulating plant development and defense responses. However, there are no extensive studies available on JAZ genes in sunflower (Helianthus annuus L.). In this study, the phylogenetic analysis of 139 putative JAZ genes from eight plants demonstrated that these JAZs could be divided into five groups (Groups I-V), and the 27 sunflower JAZs (HaJAZs) were classified into these five groups. All groups contained genes from both monocotyledons and dicotyledons, indicating that the emergence of JAZ genes predates the differentiation of monocotyledons and dicotyledons. Both segmental and tandem duplications contributed greatly to this gene family's expansion in sunflower, especially in Group II. Moreover, the expression profiles of HaJAZ genes under normal conditions, hormone treatments or abiotic stresses were analyzed based on RNA-seq data. HaJAZ2 may be undergoing pseudogenization as a nonfunctional gene because it was not expressed in any tissue. Many HaJAZ genes in roots upregulated their expression when involved in responding to exogenous hormones, especially methyl-jasmonate. The abiotic stress treatments of sunflower showed that HaJAZ5, HaJAZ15, HaJAZ17, HaJAZ20, and HaJAZ21 tend to be sensitive to certain abiotic stresses. HaJAZs from different groups may share similar functions but also exercise their unique functions when responding to abiotic stresses. We speculated that this gene family was conserved in sequence but varied in its expression among duplicated HaJAZ genes, which implies that they may confer neofunctionalization in the adaptation to abiotic stresses; this work provides insight into the resistance of sunflowers and their adaptation to diverse environmental conditions.

Tissue-specific impacts of aging and genetics on gene expression patterns in humans

Age is the primary risk factor for many common human diseases. Here, we quantify the relative contributions of genetics and aging to gene expression patterns across 27 tissues from 948 humans. We show that the predictive power of expression quantitative trait loci is impacted by age in many tissues. Jointly modelling the contributions of age and genetics to transcript level variation we find expression heritability (h2) is consistent among tissues while the contribution of aging varies by >20-fold with [Formula: see text] in 5 tissues. We find that while the force of purifying selection is stronger on genes expressed early versus late in life (Medawar's hypothesis), several highly proliferative tissues exhibit the opposite pattern. These non-Medawarian tissues exhibit high rates of cancer and age-of-expression-associated somatic mutations. In contrast, genes under genetic control are under relaxed constraint. Together, we demonstrate the distinct roles of aging and genetics on expression phenotypes.

Rhythmicity is linked to expression cost at the protein level but to expression precision at the mRNA level

Many genes have nycthemeral rhythms of expression, i.e. a 24-hours periodic variation, at either mRNA or protein level or both, and most rhythmic genes are tissue-specific. Here, we investigate and discuss the evolutionary origins of rhythms in gene expression. Our results suggest that rhythmicity of protein expression could have been favored by selection to minimize costs. Trends are consistent in bacteria, plants and animals, and are also supported by tissue-specific patterns in mouse. Unlike for protein level, cost cannot explain rhythm at the RNA level. We suggest that instead it allows to periodically reduce expression noise. Noise control had the strongest support in mouse, with limited evidence in other species. We have also found that genes under stronger purifying selection are rhythmically expressed at the mRNA level, and we propose that this is because they are noise sensitive genes. Finally, the adaptive role of rhythmic expression is supported by rhythmic genes being highly expressed yet tissue-specific. This provides a good evolutionary explanation for the observation that nycthemeral rhythms are often tissue-specific.

For many genes, their expression, i.e. the production of RNA and proteins, is rhythmic with a 24-hour period. Here, we study and discuss the evolutionary origins of these rhythms. Our analyses of data from different species suggest that the rhythmicity of protein level may have been favored by selection for cost minimization. Furthermore, we have shown that cost cannot explain the rhythmic variations in RNA levels. Instead, we suggest that it periodically reduces the stochasticity of gene expression. We also found that genes under stronger purifying selection are rhythmically expressed at the mRNA level, and propose that this is because they are noise-sensitive genes. Finally, rhythmic expression involves genes that are often highly expressed and tissue-specific. This provides a good evolutionary explanation for the tissue-specificity of these rhythms.