Using GOstats to test gene lists for GO term association

Genome-centric investigation of bile acid metabolizing microbiota of dairy cows and associated diet-induced functional implications

Although the importance of bile acid (BA)-related microbial strains and enzymes is increasingly recognized for monogastric animals, a lack of knowledge about BA metabolism in dairy cows limits functional applications aimed at the targeted modulation of microbe-host interactions for animal production and health. In the present study, 108 content samples from six intestinal regions of dairy cows were used for shotgun metagenomic sequencing. Overall, 372 high-quality metagenome-assembled genomes (MAGs) were involved in BA deconjugation, oxidation, and dehydroxylation pathways. Furthermore, the BA-metabolizing microbiome predominately occurred in the large intestine, resulting in the accumulation of secondary unconjugated BAs. Comparative genomic analysis revealed that the bile salt hydrolase (BSH)-carrying microbial populations managed with the selective environment of the dairy cow intestine by adopting numerous host mucin glycan-degrading abilities. A sequence similarity network analysis classified 439 BSH homologs into 12 clusters and identified different clusters with diverse evolution, taxonomy, signal peptides, and ecological niches. Our omics data further revealed that the strains of Firmicutes bacterium CAG-110 processed the increased abundance of BSHs from Cluster 1, coinciding with the changes in the colon cholic acid concentration after grain introduction, and were intricately related to intestinal inflammation. This study is the first to use a genome-centric approach and whole intestine-targeted metabolomics to reveal microbial BA metabolism and its diet-induced functional implications in dairy cows. These findings provide insight into the manipulation of intestinal microorganisms for improving host health.

Dynamic regulation of gonadal transposon control across the lifespan of the naturally short-lived African turquoise killifish

Although germline cells are considered to be functionally "immortal," both the germline and supporting somatic cells in the gonad within an organism experience aging. With increased age at parenthood, the age-related decline in reproductive success has become an important biological issue for an aging population. However, molecular mechanisms underlying reproductive aging across sexes in vertebrates remain poorly understood. To decipher molecular drivers of vertebrate gonadal aging across sexes, we perform longitudinal characterization of the gonadal transcriptome throughout the lifespan in the naturally short-lived African turquoise killifish (Nothobranchius furzeri). By combining mRNA-seq and small RNA-seq from 26 individuals, we characterize the aging gonads of young-adult, middle-aged, and old female and male fish. We analyze changes in transcriptional patterns of genes, transposable elements (TEs), and piRNAs. We find that testes seem to undergo only marginal changes during aging. In contrast, in middle-aged ovaries, the time point associated with peak female fertility in this strain, PIWI pathway components are transiently down-regulated, TE transcription is elevated, and piRNA levels generally decrease, suggesting that egg quality may already be declining at middle-age. Furthermore, we show that piRNA ping-pong biogenesis declines steadily with age in ovaries, whereas it is maintained in aging testes. To our knowledge, this data set represents the most comprehensive transcriptomic data set for vertebrate gonadal aging. This resource also highlights important pathways that are regulated during reproductive aging in either ovaries or testes, which could ultimately be leveraged to help restore aspects of youthful reproductive function.

Transcription factor PagLBD21 functions as a repressor of secondary xylem development in Populus

During secondary growth in trees, xylem cells differentiated from cambium cells begin to synthesize secondary cell walls that are primarily composed of cellulose, hemicellulose and lignin, and are deposited between primary cell walls and plasma membranes, leading to wood formation. Identification of regulatory genes functioning during this developmental process is valuable for increasing wood production. In this study, we functionally characterized an LBD (LATERAL ORGAN BOUNDARIES DOMAIN) transcription factor PagLBD21 as a repressor of secondary xylem development in Populus. Compared to wild type plants, transgenic plants overexpressing PagLBD21 (PagLBD21OE) exhibited decreased xylem widths in cross-sections. Consistent with the functional analysis, RNA sequencing (RNA-seq) analysis revealed that genes functioning in xylem development and secondary cell wall biosynthesis pathways were significantly down-regulated in PagLBD21OE plants. We also performed DNA affinity purification followed by sequencing (DAP-seq) to identify genome-wide target genes of PagLBD21. Furthermore, we compared the RNA-seq and DAP-seq datasets of PagLBD21 and PagLBD3, and the results showed that there was a significant overlap between their target genes, suggesting these two LBD transcription factors are functionally redundant during secondary growth.

Experimental Life History Evolution Results in Sex-specific Evolution of Gene Expression in Seed Beetles

The patterns of reproductive timing and senescence vary within and across species owing to differences in reproductive strategies, but our understanding of the molecular underpinnings of such variation is incomplete. This is perhaps particularly true for sex differences. We investigated the evolution of sex-specific gene expression associated with life history divergence in replicated populations of the seed beetle Acanthoscelides obtectus, experimentally evolving under (E)arly or (L)ate life reproduction for >200 generations which has resulted in strongly divergent life histories. We detected 1,646 genes that were differentially expressed in E and L lines, consistent with a highly polygenic basis of life history evolution. Only 30% of differentially expressed genes were similarly affected in males and females. The evolution of long life was associated with significantly reduced sex differences in expression, especially in non-reproductive tissues. The expression differences were overall more pronounced in females, in accordance with their greater phenotypic divergence in lifespan. Functional enrichment analysis revealed differences between E and L beetles in gene categories previously implicated in aging, such as mitochondrial function and defense response. The results show that divergent life history evolution can be associated with profound changes in gene expression that alter the transcriptome in a sex-specific way, highlighting the importance of understanding the mechanisms of aging in each sex.

Vector-virus interaction affects viral loads and co-occurrence

BACKGROUND

Vector-borne viral diseases threaten human and wildlife worldwide. Vectors are often viewed as a passive syringe injecting the virus. However, to survive, replicate and spread, viruses must manipulate vector biology. While most vector-borne viral research focuses on vectors transmitting a single virus, in reality, vectors often carry diverse viruses. Yet how viruses affect the vectors remains poorly understood. Here, we focused on the varroa mite (Varroa destructor), an emergent parasite that can carry over 20 honey bee viruses, and has been responsible for colony collapses worldwide, as well as changes in global viral populations. Co-evolution of the varroa and the viral community makes it possible to investigate whether viruses affect vector gene expression and whether these interactions affect viral epidemiology.

RESULTS

Using a large set of available varroa transcriptomes, we identified how abundances of individual viruses affect the vector's transcriptional network. We found no evidence of competition between viruses, but rather that some virus abundances are positively correlated. Furthermore, viruses that are found together interact with the vector's gene co-expression modules in similar ways, suggesting that interactions with the vector affect viral epidemiology. We experimentally validated this observation by silencing candidate genes using RNAi and found that the reduction in varroa gene expression was accompanied by a change in viral load.

CONCLUSIONS

Combined, the meta-transcriptomic analysis and experimental results shed light on the mechanism by which viruses interact with each other and with their vector to shape the disease course.

Mitochondrial dysfunction induces ALK5-SMAD2-mediated hypovascularization and arteriovenous malformations in mouse retinas

Although mitochondrial activity is critical for angiogenesis, its mechanism is not entirely clear. Here we show that mice with endothelial deficiency of any one of the three nuclear genes encoding for mitochondrial proteins, transcriptional factor (TFAM), respiratory complex IV component (COX10), or redox protein thioredoxin 2 (TRX2), exhibit retarded retinal vessel growth and arteriovenous malformations (AVM). Single-cell RNA-seq analyses indicate that retinal ECs from the three mutant mice have increased TGFβ signaling and altered gene expressions associated with vascular maturation and extracellular matrix, correlating with vascular malformation and increased basement membrane thickening in microvesels of mutant retinas. Mechanistic studies suggest that mitochondrial dysfunction from Tfam, Cox10, or Trx2 depletion induces a mitochondrial localization and MAPKs-mediated phosphorylation of SMAD2, leading to enhanced ALK5-SMAD2 signaling. Importantly, pharmacological blockade of ALK5 signaling or genetic deficiency of SMAD2 prevented retinal vessel growth retardation and AVM in all three mutant mice. Our studies uncover a novel mechanism whereby mitochondrial dysfunction via the ALK5-SMAD2 signaling induces retinal vascular malformations, and have therapeutic values for the alleviation of angiogenesis-associated human retinal diseases.

Establishment and validation of a predictive model of preeclampsia based on transcriptional signatures of 43 genes in decidua basalis and peripheral blood

Preeclampsia (PE) has an increasing incidence worldwide, and there is no gold standard for prediction. Recent progress has shown that abnormal decidualization and impaired vascular remodeling are essential to PE pathogenesis. Therefore, it is of great significance to analyze the decidua basalis and blood changes of PE to explore new methods. Here, we performed weighted gene co-expression network analysis based on 9553 differentially expressed genes of decidua basalis data (GSE60438 includes 25 cases of PE and 23 non-cases) from Gene Expression Omnibus to screen relevant module-eigengenes (MEs). Among them, MEblue and MEgrey are the most correlated with PE, which contains 371 core genes. Subsequently, we applied the logistic least absolute shrinkage and selection operator regression, screened 43 genes most relevant to prediction from the intersections of the 371 genes and training set (GSE48424 includes 18 cases of PE and 18 non-cases) genes, and built a predictive model. The specificity and sensitivity are illustrated by receiver operating characteristic curves, and the stability was verified by two validation sets (GSE86200 includes 12 cases of PE and 48 non-cases, and GSE85307 includes 47 cases of PE and 110 non-cases). The results demonstrated that our predictive model shows good predictions, with an area under the curve of 0.991 for the training set, 0.874 and 0.986 for the validation sets. Finally, we found the 43 key marker genes in the model are closely associated with the clinically accepted predictive molecules, including FLT1, PIGF, ENG and VEGF. Therefore, this predictive model provides a potential approach for PE diagnosis and treatment.

Experimental inheritance of antibiotic acquired dysbiosis affects host phenotypes across generations

Microbiomes can enhance the health, fitness and even evolutionary potential of their hosts. Many organisms propagate favorable microbiomes fully or partially via vertical transmission. In the long term, such co-propagation can lead to the evolution of specialized microbiomes and functional interdependencies with the host. However, microbiomes are vulnerable to environmental stressors, particularly anthropogenic disturbance such as antibiotics, resulting in dysbiosis. In cases where microbiome transmission occurs, a disrupted microbiome may then become a contagious pathology causing harm to the host across generations. We tested this hypothesis using the specialized socially transmitted gut microbiome of honey bees as a model system. By experimentally passaging tetracycline-treated microbiomes across worker 'generations' we found that an environmentally acquired dysbiotic phenotype is heritable. As expected, the antibiotic treatment disrupted the microbiome, eliminating several common and functionally important taxa and strains. When transmitted, the dysbiotic microbiome harmed the host in subsequent generations. Particularly, naïve bees receiving antibiotic-altered microbiomes died at higher rates when challenged with further antibiotic stress. Bees with inherited dysbiotic microbiomes showed alterations in gene expression linked to metabolism and immunity, among other pathways, suggesting effects on host physiology. These results indicate that there is a possibility that sublethal exposure to chemical stressors, such as antibiotics, may cause long-lasting changes to functional host-microbiome relationships, possibly weakening the host's progeny in the face of future ecological challenges. Future studies under natural conditions would be important to examine the extent to which negative microbiome-mediated phenotypes could indeed be heritable and what role this may play in the ongoing loss of biodiversity.

Sequestration and efflux largely account for cadmium and copper resistance in the deep-sea Nitratiruptor sp. SB155-2 (phylum Campylobacterota)

In deep-sea hydrothermal vent environments, metal-enriched fluids and sediments abound, making these habitats ideal to study metal resistance in prokaryotes. In this investigation, we employed transcriptomics and shotgun proteomics with scanning transmission electron microscopy and energy-dispersive x-ray spectroscopy (STEM-EDX) to better understand mechanisms of tolerance for cadmium (Cd) and copper (Cu) at stress-inducing concentrations in Nitratiruptor sp. SB155-2 (phylum Campylobacterota). Transcriptomic profiles were remarkably different in the presence of these two metals, displaying 385 (19%) and 629 (31%) differentially transcribed genes (DTG) in the presence of Cd(II) and Cu(II), respectively, while only 7% of differentially transcribed (DT) genes were shared, with genes for non-specific metal transporters and genes involved in oxidative stress-response predominating. Transcriptomic and proteomic analyses confirmed that metal-specific DT pathways under Cu(II) stress, including those involving sulfur, cysteine, and methionine, are likely required for high-affinity efflux systems, while flagella formation and chemotaxis were over-represented under Cd(II) stress. Consistent with these differences, STEM-EDX analysis revealed that polyphosphate-like granules (pPLG), the formation of CdS particles, and the periplasmic space are crucial for Cd(II) sequestration. Overall, this study provides new insights regarding metal-specific adaptations of Campylobacterota to deep-sea hydrothermal vent environments.

Transcriptomic response to ISAV infection in the gills, head kidney and spleen of resistant and susceptible Atlantic salmon

BACKGROUND

Infectious Salmon Anaemia virus (ISAV) is an orthomyxovirus responsible for large losses in Atlantic salmon (Salmo salar) aquaculture. Current available treatments and vaccines are not fully effective, and therefore selective breeding to produce ISAV-resistant strains of Atlantic salmon is a high priority for the industry. Genomic selection and potentially genome editing can be applied to enhance the disease resistance of aquaculture stocks, and both approaches can benefit from increased knowledge on the genomic mechanisms of resistance to ISAV. To improve our understanding of the mechanisms underlying resistance to ISAV in Atlantic salmon we performed a transcriptomic study in ISAV-infected salmon with contrasting levels of resistance to this virus.

RESULTS

Three different tissues (gills, head kidney and spleen) were collected on 12 resistant and 12 susceptible fish at three timepoints (pre-challenge, 7 and 14 days post challenge) and RNA sequenced. The transcriptomes of infected and non-infected fish and of resistant and susceptible fish were compared at each timepoint. The results show that the responses to ISAV are organ-specific; an important response to the infection was observed in the head kidney, with up-regulation of immune processes such as interferon and NLR pathways, while in gills and spleen the response was more moderate. In addition to immune related genes, our results suggest that other processes such as ubiquitination and ribosomal processing are important during early infection with ISAV. Moreover, the comparison between resistant and susceptible fish has also highlighted some interesting genes related to ubiquitination, intracellular transport and the inflammasome.

CONCLUSIONS

Atlantic salmon infection by ISAV revealed an organ-specific response, implying differential function during the infection. An immune response was observed in the head kidney in these early timepoints, while gills and spleen showed modest responses in comparison. Comparison between resistance and susceptible samples have highlighted genes of interest for further studies, for instance those related to ubiquitination or the inflammasome.

DIO3 protects against thyrotoxicosis-derived cranio-encephalic and cardiac congenital abnormalities

Maternal hyperthyroidism is associated with an increased incidence of congenital abnormalities at birth, but it is not clear which of these defects arise from a transient developmental excess of thyroid hormone and which depend on pregnancy stage, antithyroid drug choice, or unwanted subsequent fetal hypothyroidism. To address this issue, we studied a mouse model of comprehensive developmental thyrotoxicosis secondary to a lack of type 3 deiodinase (DIO3). Dio3-/- mice exhibited reduced neonatal viability on most genetic backgrounds and perinatal lethality on a C57BL/6 background. Dio3-/- mice exhibited severe growth retardation during the neonatal period and cartilage loss. Mice surviving after birth manifested brain and cranial dysmorphisms, severe hydrocephalus, choanal atresia, and cleft palate. These abnormalities were noticeable in C57BL/6J Dio3-/- mice at fetal stages, in addition to a thyrotoxic heart with septal defects and thin ventricular walls. Our findings stress the protecting role of DIO3 during development and support the hypothesis that human congenital abnormalities associated with hyperthyroidism during pregnancy are caused by transient thyrotoxicosis before clinical intervention. Our results also suggest thyroid hormone involvement in the etiology of idiopathic pathologies including cleft palate, choanal atresia, Chiari malformations, Kaschin-Beck disease, and Temple and other cranio-encephalic and heart syndromes.

Buffy Coat Transcriptomic Analysis Reveals Alterations in Host Cell Protein Synthesis and Cell Cycle in Severe COVID-19 Patients

Transcriptome studies have reported the dysregulation of cell cycle-related genes and the global inhibition of host mRNA translation in COVID-19 cases. However, the key genes and cellular mechanisms that are most affected by the severe outcome of this disease remain unclear. For this work, the RNA-seq approach was used to study the differential expression in buffy coat cells of two groups of people infected with SARS-CoV-2: (a) Mild, with mild symptoms; and (b) SARS (Severe Acute Respiratory Syndrome), who were admitted to the intensive care unit with the severe COVID-19 outcome. Transcriptomic analysis revealed 1009 up-regulated and 501 down-regulated genes in the SARS group, with 10% of both being composed of long non-coding RNA. Ribosome and cell cycle pathways were enriched among down-regulated genes. The most connected proteins among the differentially expressed genes involved transport dysregulation, proteasome degradation, interferon response, cytokinesis failure, and host translation inhibition. Furthermore, interactome analysis showed Fibrillarin to be one of the key genes affected by SARS-CoV-2. This protein interacts directly with the N protein and long non-coding RNAs affecting transcription, translation, and ribosomal processes. This work reveals a group of dysregulated processes, including translation and cell cycle, as key pathways altered in severe COVID-19 outcomes.

Human Plasma Proteome During Normal Pregnancy

The human plasma proteome is underexplored despite its potential value for monitoring health and disease. Herein, using a recently developed aptamer-based platform, we profiled 7288 proteins in 528 plasma samples from 91 normal pregnancies (Gene Expression Omnibus identifier GSE206454). The coefficient of variation was <20% for 93% of analytes (median 7%), and a cross-platform correlation for selected key angiogenic and anti-angiogenic proteins was significant. Gestational age was associated with changes in 953 proteins, including highly modulated placenta- and decidua-specific proteins, and they were enriched in biological processes including regulation of growth, angiogenesis, immunity, and inflammation. The abundance of proteins corresponding to RNAs specific to populations of cells previously described by single-cell RNA-Seq analysis of the placenta was highly modulated throughout gestation. Furthermore, machine learning-based prediction of gestational age and of time from sampling to term delivery compared favorably with transcriptomic models (mean absolute error of 2 weeks). These results suggested that the plasma proteome may provide a non-invasive readout of placental cellular dynamics and serve as a blueprint for investigating obstetrical disease.

Aortic Cellular Diversity and Quantitative Genome-Wide Association Study Trait Prioritization Through Single-Nuclear RNA Sequencing of the Aneurysmal Human Aorta

BACKGROUND

Mural cells in ascending aortic aneurysms undergo phenotypic changes that promote extracellular matrix destruction and structural weakening. To explore this biology, we analyzed the transcriptional features of thoracic aortic tissue.

METHODS

Single-nuclear RNA sequencing was performed on 13 samples from human donors, 6 with thoracic aortic aneurysm, and 7 without aneurysm. Individual transcriptomes were then clustered based on transcriptional profiles. Clusters were used for between-disease differential gene expression analyses, subcluster analysis, and analyzed for intersection with genetic aortic trait data.

RESULTS

We sequenced 71 689 nuclei from human thoracic aortas and identified 14 clusters, aligning with 11 cell types, predominantly vascular smooth muscle cells (VSMCs) consistent with aortic histology. With unbiased methodology, we found 7 vascular smooth muscle cell and 6 fibroblast subclusters. Differentially expressed genes analysis revealed a vascular smooth muscle cell group accounting for the majority of differential gene expression. Fibroblast populations in aneurysm exhibit distinct behavior with almost complete disappearance of quiescent fibroblasts. Differentially expressed genes were used to prioritize genes at aortic diameter and distensibility genome-wide association study loci highlighting the genes JUN, LTBP4 (latent transforming growth factor beta-binding protein 1), and IL34 (interleukin 34) in fibroblasts, ENTPD1, PDLIM5 (PDZ and LIM domain 5), ACTN4 (alpha-actinin-4), and GLRX in vascular smooth muscle cells, as well as LRP1 in macrophage populations.

CONCLUSIONS

Using nuclear RNA sequencing, we describe the cellular diversity of healthy and aneurysmal human ascending aorta. Sporadic aortic aneurysm is characterized by differential gene expression within known cellular classes rather than by the appearance of novel cellular forms. Single-nuclear RNA sequencing of aortic tissue can be used to prioritize genes at aortic trait loci.

Adapting to climate with limited genetic diversity: Nucleotide, DNA methylation and microbiome variation among populations of the social spider Stegodyphus dumicola

Understanding the role of genetic and nongenetic variants in modulating phenotypes is central to our knowledge of adaptive responses to local conditions and environmental change, particularly in species with such low population genetic diversity that it is likely to limit their evolutionary potential. A first step towards uncovering the molecular mechanisms underlying population-specific responses to the environment is to carry out environmental association studies. We associated climatic variation with genetic, epigenetic and microbiome variation in populations of a social spider with extremely low standing genetic diversity. We identified genetic variants that are associated strongly with environmental variation, particularly with average temperature, a pattern consistent with local adaptation. Variation in DNA methylation in many genes was strongly correlated with a wide set of climate parameters, thereby revealing a different pattern of associations than that of genetic variants, which show strong correlations to a more restricted range of climate parameters. DNA methylation levels were largely independent of cis-genetic variation and of overall genetic population structure, suggesting that DNA methylation can work as an independent mechanism. Microbiome composition also correlated with environmental variation, but most strong associations were with precipitation-related climatic factors. Our results suggest a role for both genetic and nongenetic mechanisms in shaping phenotypic responses to local environments.

A possible genomic footprint of polygenic adaptation on population divergence in seed beetles?

Efforts to unravel the genomic basis of incipient speciation are hampered by a mismatch between our toolkit and our understanding of the ecology and genetics of adaptation. While the former is focused on detecting selective sweeps involving few independently acting or linked speciation genes, the latter states that divergence typically occurs in polygenic traits under stabilizing selection. Here, we ask whether a role of stabilizing selection on polygenic traits in population divergence may be unveiled by using a phenotypically informed integrative approach, based on genome‐wide variation segregating in divergent populations. We compare three divergent populations of seed beetles (Callosobruchus maculatus) where previous work has demonstrated a prominent role for stabilizing selection on, and population divergence in, key life history traits that reflect rate‐dependent metabolic processes. We derive and assess predictions regarding the expected pattern of covariation between genetic variation segregating within populations and genetic differentiation between populations. Population differentiation was considerable (mean F_ST = 0.23–0.26) and was primarily built by genes showing high selective constraints and an imbalance in inferred selection in different populations (positive Tajima's D_NS in one and negative in one), and this set of genes was enriched with genes with a metabolic function. Repeatability of relative population differentiation was low at the level of individual genes but higher at the level of broad functional classes, again spotlighting metabolic genes. Absolute differentiation (d_XY) showed a very different general pattern at this scale of divergence, more consistent with an important role for genetic drift. Although our exploration is consistent with stabilizing selection on polygenic metabolic phenotypes as an important engine of genome‐wide relative population divergence and incipient speciation in our study system, we note that it is exceedingly difficult to firmly exclude other scenarios.

Effects of urinary extracellular vesicles from prostate cancer patients on the transcriptomes of cancer-associated and normal fibroblasts

BACKGROUND

Increasing evidence suggests that cancer-derived extracellular vesicles (EVs) alter the phenotype and functions of fibroblasts and trigger the reprogramming of normal fibroblasts into cancer-associated fibroblasts (CAFs). Here, we for the first time studied the effects of urinary EVs from PC patients and healthy males on the transcriptional landscape of prostate CAFs and normal foreskin fibroblasts.

METHODS

Patient-derived prostate fibroblast primary cultures PCF-54 and PCF-55 were established from two specimens of PC tissues. EVs were isolated from urine samples of 3 patients with PC and 2 healthy males and used for the treatment of prostate fibroblast primary cultures and normal foreskin fibroblasts. The EV-treated fibroblasts were subjected to RNA sequencing analysis.

RESULTS

RNA sequencing analysis showed that the fibroblast cultures differed significantly in their response to urinary EVs. The transcriptional response of foreskin fibroblasts to the urinary EVs isolated from PC patients and healthy controls was very similar and mostly related to the normal functions of fibroblasts. On the contrary, PCF-54 cells responded very differently - EVs from PC patients elicited transcriptional changes related to the regulation of the cell division and chromosome segregation, whereas EVs from healthy males affected mitochondrial respiration. In PCF-55 cells, EVs from both, PC-patients and controls induced the expression of a number of chemokines such as CCL2, CCL13, CXCL1, CXCL8, whereas pathways related to regulation of apoptotic signaling and production of cell adhesion molecules were triggered specifically by EVs from PC patients.

CONCLUSION

This study demonstrates that urinary EVs from PC patients and healthy controls elicit distinct transcriptional responses in prostate CAFs and supports the idea that EVs contribute to the generation of functional heterogeneity of CAFs. Moreover, this study suggests that the changes in the gene expression pattern in EV recipient cells might serve as a novel type of functional cancer biomarkers.

Transcriptomic and Metabolomic Analysis of a Pseudomonas-Resistant versus a Susceptible Arabidopsis Accession

Accessions of one plant species may show significantly different levels of susceptibility to stresses. The Arabidopsis thaliana accessions Col-0 and C24 differ significantly in their resistance to the pathogen Pseudomonas syringae pv. tomato (Pst). To help unravel the underlying mechanisms contributing to this naturally occurring variance in resistance to Pst, we analyzed changes in transcripts and compounds from primary and secondary metabolism of Col-0 and C24 at different time points after infection with Pst. Our results show that the differences in the resistance of Col-0 and C24 mainly involve mechanisms of salicylic-acid-dependent systemic acquired resistance, while responses of jasmonic-acid-dependent mechanisms are shared between the two accessions. In addition, arginine metabolism and differential activity of the biosynthesis pathways of aliphatic glucosinolates and indole glucosinolates may also contribute to the resistance. Thus, this study highlights the difference in the defense response strategies utilized by different genotypes.

Profiling of epigenetic marker regions in murine ILCs under homeostatic and inflammatory conditions

Epigenetic modifications such as DNA methylation play an essential role in imprinting specific transcriptional patterns in cells. We performed genome-wide DNA methylation profiling of murine lymph node–derived ILCs, which led to the identification of differentially methylated regions (DMRs) and the definition of epigenetic marker regions in ILCs. Marker regions were located in genes with a described function for ILCs, such as Tbx21, Gata3, or Il23r, but also in genes that have not been related to ILC biology. Methylation levels of the marker regions and expression of the associated genes were strongly correlated, indicating their functional relevance. Comparison with T helper cell methylomes revealed clear lineage differences, despite partial similarities in the methylation of specific ILC marker regions. IL-33–mediated challenge affected methylation of ILC2 epigenetic marker regions in the liver, while remaining relatively stable in the lung. In our study, we identified a set of epigenetic markers that can serve as a tool to study phenotypic and functional properties of ILCs.

The gut microbiota affects the social network of honeybees

The gut microbiota influences animal neurodevelopment and behaviour but has not previously been documented to affect group-level properties of social organisms. Here, we use honeybees to probe the effect of the gut microbiota on host social behaviour. We found that the microbiota increased the rate and specialization of head-to-head interactions between bees. Microbiota colonization was associated with higher abundances of one-third of the metabolites detected in the brain, including amino acids with roles in synaptic transmission and brain energetic function. Some of these metabolites were significant predictors of the number of social interactions. Microbiota colonization also affected brain transcriptional processes related to amino acid metabolism and epigenetic modifications in a brain region involved in sensory perception. These results demonstrate that the gut microbiota modulates the emergent colony social network of honeybees and suggest changes in chromatin accessibility and amino acid biosynthesis as underlying processes.

Claudin-6 increases SNAI1, NANOG and SOX2 gene expression in human gastric adenocarcinoma AGS cells

BACKGROUND

Gastric cancer is a heterogeneous disease associated to deregulated gastric epithelia tight junction barrier function and di novo expression of claudin-6; these changes are associated with epithelial-mesenchymal transition, enhanced invasiveness, metastatic progression, resistance to chemotherapy, and poor prognosis. Gastric cancer stem cells represent a rare population of cells within the tumor implicated in tumor growth and higher tumorigenic capacity. The possible relation between claudin-6 expression and the expression of some markers associated to epithelial mesenchymal transition and cancer stem cells in gastric cancer cells have never been explored.

METHODS AND RESULTS

CD44, CD24, Twist, Villin, DCLK1, claudin-6, NANOG, E-Cadherin, SOX2, and SNAI1 expression was evaluated by immunofluorescence and cytofluorometry in wild type and Claudin-6 transfected AGS cells. Cell migration assays were also performed. Differentially expressed genes and biological processes analysis was performed to determine gene preponderance. The results showed that claudin-6 overexpression enriched the CD44 + /CD24- subpopulation with an overall increase in the expression and the number of CD44 + cells. A significant increase in NANOG, SOX2 and SNAI1 expression and enhanced cell migration was observed in claudin-6 transfected cells. Transcriptome analysis revealed 271 genes involved in enhanced biological processes with only 31 with a significantly p value; thirteen of those genes are closely associated to epithelial mesenchymal transition processes and folding and unfolding processes of proteins in the endoplasmic reticulum.

CONCLUSIONS

The pro-tumorigenic effect of claudin-6 in gastric cancer could be associated to dedifferentiation of epithelial cells and an increase in di novo cancer stem cell genesis.

Class I TCP transcription factor AtTCP8 modulates key brassinosteroid-responsive genes

The plant-specific TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) transcription factor family is most closely associated with regulating plant developmental programs. Recently, TCPs were also shown to mediate host immune signaling, both as targets of pathogen virulence factors and as regulators of plant defense genes. However, comprehensive characterization of TCP gene targets is still lacking. Loss of function of the class I TCP gene AtTCP8 attenuates early immune signaling and, when combined with mutations in AtTCP14 and AtTCP15, additional layers of defense signaling in Arabidopsis (Arabidopsis thaliana). Here, we focus on TCP8, the most poorly characterized of the three to date. We used chromatin immunoprecipitation and RNA sequencing to identify TCP8-bound gene promoters and differentially regulated genes in the tcp8 mutant; these datasets were heavily enriched in signaling components for multiple phytohormone pathways, including brassinosteroids (BRs), auxin, and jasmonic acid. Using BR signaling as a representative example, we showed that TCP8 directly binds and activates the promoters of the key BR transcriptional regulatory genes BRASSINAZOLE-RESISTANT1 (BZR1) and BRASSINAZOLE-RESISTANT2 (BZR2/BES1). Furthermore, tcp8 mutant seedlings exhibited altered BR-responsive growth patterns and complementary reductions in BZR2 transcript levels, while TCP8 protein demonstrated BR-responsive changes in subnuclear localization and transcriptional activity. We conclude that one explanation for the substantial targeting of TCP8 alongside other TCP family members by pathogen effectors may lie in its role as a modulator of BR and other plant hormone signaling pathways.

Erythropoietin-driven dynamic proteome adaptations during erythropoiesis prevent iron overload in the developing embryo

Erythropoietin (Epo) ensures survival and proliferation of colony-forming unit erythroid (CFU-E) progenitor cells and their differentiation to hemoglobin-containing mature erythrocytes. A lack of Epo-induced responses causes embryonic lethality, but mechanisms regulating the dynamic communication of cellular alterations to the organismal level remain unresolved. By time-resolved transcriptomics and proteomics, we show that Epo induces in CFU-E cells a gradual transition from proliferation signature proteins to proteins indicative for differentiation, including heme-synthesis enzymes. In the absence of the Epo receptor (EpoR) in embryos, we observe a lack of hemoglobin in CFU-E cells and massive iron overload of the fetal liver pointing to a miscommunication between liver and placenta. A reduction of iron-sulfur cluster-containing proteins involved in oxidative phosphorylation in these embryos leads to a metabolic shift toward glycolysis. This link connecting erythropoiesis with the regulation of iron homeostasis and metabolic reprogramming suggests that balancing these interactions is crucial for protection from iron intoxication and for survival.

Global proteomic profiling of multiple organs of cats (Felis catus) and proteome-transcriptome correlation during acute Toxoplasma gondii infection

Background

Toxoplasma gondii is a protozoan parasite which can infect almost all warm-blooded animals and humans. Understanding the differential expression of proteins and transcripts associated with T. gondii infection in its definitive host (cat) may improve our knowledge of how the parasite manipulates the molecular microenvironment of its definitive host. The aim of this study was to explore the global proteomic alterations in the major organs of cats during acute T. gondii infection.

Methods

iTRAQ-based quantitative proteomic profiling was performed on six organs (brain, liver, lung, spleen, heart and small intestine) of cats on day 7 post-infection by cysts of T. gondii PRU strain (Genotype II). Mascot software was used to conduct the student’s t-test. Proteins with P values < 0.05 and fold change > 1.2 or < 0.83 were considered as differentially expressed proteins (DEPs).

Results

A total of 32,657 proteins were identified in the six organs, including 2556 DEPs; of which 1325 were up-regulated and 1231 were down-regulated. The brain, liver, lung, spleen, heart and small intestine exhibited 125 DEPs, 463 DEPs, 255 DEPs, 283 DEPs, 855 DEPs and 575 DEPs, respectively. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of all proteins and DEPs in all organs showed that many proteins were enriched in binding, cell part, cell growth and death, signal transduction, translation, sorting and degradation, extracellular matrix remodeling, tryptophan catabolism, and immune system. Correlations between differentially expressed proteins and transcripts were detected in the liver (n = 19), small intestine (n = 17), heart (n = 9), lung (n = 9) and spleen (n = 3).

Conclusions

The present study identified 2556 DEPs in six cat tissues on day 7 after infection by T. gondii PRU strain, and functional enrichment analyses showed that these DEPs were associated with various cellular and metabolic processes. These findings provide a solid base for further in-depth investigation of the complex proteotranscriptomic reprogramming that mediates the dynamic interplays between T. gondii and the different feline tissues.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40249-022-01022-7.

Single-cell RNA-seq uncovered hemocyte functional subtypes and their differentiational characteristics and connectivity with morphological subpopulations in Litopenaeus vannamei

Hemocytes play central roles in shrimp immune system, whereas whose subclasses have not yet been completely defined. At present, the morphological classification of hemocytes is inadequate to classify the complete hemocyte repertoire and elucidate the functions and differentiation and maturation processes. Based on single-cell RNA sequencing (scRNA-seq) of hemocytes in healthy Litopenaeus vannamei, combined with RNA-FISH and flow cytometric sorting, we identified three hemocyte clusters including TGase+ cells, CTL+ cells and Crustin+ cells, and further determined their functional properties, potential differentiation trajectory and correspondence with morphological subpopulations. The TGase+ cells were mainly responsible for the coagulation, exhibiting distinguishable characteristics of hyalinocyte, and appeared to be developmentally arrested at an early stage of hemocyte differentiation. The CTL+ cells and Crustin+ cells arrested at terminal stages of differentiation mainly participated in recognizing foreign pathogens and initiating immune defense responses, owning distinctive features of granule-containing hemocytes. Furthermore, we have revealed the functional sub-clusters of three hemocyte clusters and their potential differentiation pathways according to the expression of genes involved in cell cycle, cell differentiation and immune response, and the successive differentiation and maturation of hyalinocytes to granule-containing hemocytes have also mapped. The results revealed the diversity of shrimp hemocytes and provide new theoretical rationale for hemocyte classification, which also facilitate systematic research on crustacean immunity.

Purifying selection enduringly acts on the sequence evolution of highly expressed proteins in Escherichia coli

The evolutionary speed of a protein sequence is constrained by its expression level, with highly expressed proteins evolving relatively slowly. This negative correlation between expression levels and evolutionary rates (known as the E–R anticorrelation) has already been widely observed in past macroevolution between species from bacteria to animals. However, it remains unclear whether this seemingly general law also governs recent evolution, including past and de novo, within a species. However, the advent of genomic sequencing and high-throughput phenotyping, particularly for bacteria, has revealed fundamental gaps between the 2 evolutionary processes and has provided empirical data opposing the possible underlying mechanisms which are widely believed. These conflicts raise questions about the generalization of the E–R anticorrelation and the relevance of plausible mechanisms. To explore the ubiquitous impact of expression levels on molecular evolution and test the relevance of the possible underlying mechanisms, we analyzed the genome sequences of 99 strains of Escherichia coli for evolution within species in nature. We also analyzed genomic mutations accumulated under laboratory conditions as a model of de novo evolution within species. Here, we show that E–R anticorrelation is significant in both past and de novo evolution within species in E. coli. Our data also confirmed ongoing purifying selection on highly expressed genes. Ongoing selection included codon-level purifying selection, supporting the relevance of the underlying mechanisms. However, the impact of codon-level purifying selection on the constraints in evolution within species might be smaller than previously expected from evolution between species.

Nitrogen and phosphorus significantly alter growth, nitrogen fixation, anatoxin-a content, and the transcriptome of the bloom-forming cyanobacterium, Dolichospermum

While freshwater cyanobacteria are traditionally thought to be limited by the availability of phosphorus (P), fixed nitrogen (N) supply can promote the growth and/or toxin production of some genera. This study characterizes how growth on N2 (control), nitrate (NO3 -), ammonium (NH4 +), and urea as well as P limitation altered the growth, toxin production, N2 fixation, and gene expression of an anatoxin-a (ATX-A) - producing strain of Dolichospermum sp. 54. The transcriptomes of fixed N and P-limited cultures differed significantly from those of fixed N-deplete, P-replete (control) cultures, while the transcriptomes of P-replete cultures amended with either NH4 + or NO3 - were not significantly different relative to those of the control. Growth rates of Dolichospermum (sp. 54) were significantly higher when grown on fixed N relative to without fixed N; growth on NH4 + was also significantly greater than growth on NO3 -. NH4 + and urea significantly lowered N2 fixation and nifD gene transcript abundance relative to the control while cultures amended with NO3 - exhibited N2 fixation and nifD gene transcript abundance that was not different from the control. Cultures grown on NH4 + exhibited the lowest ATX-A content per cell and lower transcript abundance of genes associated ATX-A synthesis (ana), while the abundance of transcripts of several ana genes were highest under fixed N and P - limited conditions. The significant negative correlation between growth rate and cellular anatoxin quota as well as the significantly higher number of transcripts of ana genes in cultures deprived of fixed N and P relative to P-replete cultures amended with NH4 + suggests ATX-A was being actively synthesized under P limitation. Collectively, these findings indicate that management strategies that do not regulate fixed N loading will leave eutrophic water bodies vulnerable to more intense and toxic (due to increased biomass) blooms of Dolichospermum.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.

Family dinner: Transcriptional plasticity of five Noctuidae (Lepidoptera) feeding on three host plant species

Polyphagous insects often show specialization in feeding on different host plants in terms of survival and growth and, therefore, can be considered minor or major pests of particular hosts. Whether polyphagous insects employ a common transcriptional response to cope with defenses from diverse host plants is under-studied. We focused on patterns of transcriptional plasticity in polyphagous moths (Noctuidae), of which many species are notorious pests, in relation to herbivore performance on different host plants. We compared the transcriptional plasticity of five polyphagous moth species feeding and developing on three different host plant species. Using a comparative phylogenetic framework, we evaluated if successful herbivory, as measured by larval performance, is determined by a shared or lineage-specific transcriptional response. The upregulated transcriptional activity, or gene expression pattern, of larvae feeding on the different host plants and artificial control diet was highly plastic and moth species-specific. Specialization, defined as high herbivore success for specific host plants, was not generally linked to a lower number of induced genes. Moths that were more distantly related and showing high herbivore success for certain host plants showed shared expression of multiple homologous genes, indicating convergence. We further observed specific transcriptional responses within phylogenetic lineages. These expression patterns for specific host plant species are likely caused by shared evolutionary histories, for example, symplesiomorphic patterns, and could therefore not be associated with herbivore success alone. Multiple gene families, with roles in plant digestion and detoxification, were widely expressed in response to host plant feeding but again showed highly moth species-specific. Consequently, high herbivore success for specific host plants is also driven by species-specific transcriptional plasticity. Thus, potential pest moths display a complex and species-specific transcriptional plasticity.

Intron size minimisation in teleosts

Background

Spliceosomal introns are parts of primary transcripts that are removed by RNA splicing. Although introns apparently do not contribute to the function of the mature transcript, in vertebrates they comprise the majority of the transcribed region increasing the metabolic cost of transcription. The persistence of long introns across evolutionary time suggests functional roles that can offset this metabolic cost. The teleosts comprise one of the largest vertebrate clades. They have unusually compact and variable genome sizes and provide a suitable system for analysing intron evolution.

Results

We have analysed intron lengths in 172 vertebrate genomes and show that teleost intron lengths are relatively short, highly variable and bimodally distributed. Introns that were long in teleosts were also found to be long in mammals and were more likely to be found in regulatory genes and to contain conserved sequences. Our results argue that intron length has decreased in parallel in a non-random manner throughout teleost evolution and represent a deviation from the ancestral state.

Conclusion

Our observations indicate an accelerated rate of intron size evolution in the teleosts and that teleost introns can be divided into two classes by their length. Teleost intron sizes have evolved primarily as a side-effect of genome size evolution and small genomes are dominated by short introns (<256 base pairs). However, a non-random subset of introns has resisted this process across the teleosts and these are more likely have functional roles in all vertebrate clades.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08760-w).

Single-nucleus profiling of human dilated and hypertrophic cardiomyopathy

Heart failure encompasses a heterogeneous set of clinical features that converge on impaired cardiac contractile function1,2 and presents a growing public health concern. Previous work has highlighted changes in both transcription and protein expression in failing hearts3,4, but may overlook molecular changes in less prevalent cell types. Here we identify extensive molecular alterations in failing hearts at single-cell resolution by performing single-nucleus RNA sequencing of nearly 600,000 nuclei in left ventricle samples from 11 hearts with dilated cardiomyopathy and 15 hearts with hypertrophic cardiomyopathy as well as 16 non-failing hearts. The transcriptional profiles of dilated or hypertrophic cardiomyopathy hearts broadly converged at the tissue and cell-type level. Further, a subset of hearts from patients with cardiomyopathy harbour a unique population of activated fibroblasts that is almost entirely absent from non-failing samples. We performed a CRISPR-knockout screen in primary human cardiac fibroblasts to evaluate this fibrotic cell state transition; knockout of genes associated with fibroblast transition resulted in a reduction of myofibroblast cell-state transition upon TGFβ1 stimulation for a subset of genes. Our results provide insights into the transcriptional diversity of the human heart in health and disease as well as new potential therapeutic targets and biomarkers for heart failure.

Identification of Chilling-Responsive Genes in Litchi chinensis by Transcriptomic Analysis Underlying Phytohormones and Antioxidant Systems

Litchi (Litchi chinensis Sonn.) is an important subtropical and tropical evergreen fruit tree that is seriously affected by chilling stress. In order to identify genes that may be involved in the response to chilling in litchi, we investigate the physiological and biochemical changes under chilling stress and construct 12 RNA-Seq libraries of leaf samples at 0, 4, 8, and 12 days of chilling. The results show that antioxidant enzymes are activated by chilling treatments. Comparing the transcriptome data of the four time points, we screen 2496 chilling-responsive genes (CRGs), from which we identify 63 genes related to the antioxidant system (AO-CRGs) and 54 ABA, 40 IAA, 37 CTK, 27 ETH, 21 BR, 13 GA, 35 JA, 29 SA, and 4 SL signal transduction-related genes. Expression pattern analysis shows that the expression trends of the 28 candidate genes detected by qRT-PCR are similar to those detected by RNA-Seq, indicating the reliability of our RNA-Seq data. Partial Least Squares Structural Equation Modeling (PLS-SEM) analysis of the RNA-Seq data suggests a model for the litchi plants in response to chilling stress that alters the expression of the plant hormone signaling-related genes, the transcription factor-encoding genes LcICE1, LcCBFs, and LcbZIPs, and the antioxidant system-related genes. This study provides candidate genes for the future breeding of litchi cultivars with high chilling resistance, and elucidates possible pathways for litchi in response to chilling using transcriptomic data.

Blood transcriptome analysis revealing aging gene expression profiles in red panda

The red panda is an endangered forest species distributed on the edge of the Qinghai Tibet Plateau. The species has been conserved in ex-situ in many countries and its survival is threatened by many diseases. Its immune system is vulnerable to age-associated alterations, which accumulate and result in a progressive deterioration that leads to an increased incidence of diseases. We identified 2,219 differentially expressed genes (DEGs) between geriatric (11-16 years) and adult individuals (4-8 years), and 1690 DEGs between adults and juveniles (1 year). The gene expression and functional annotation results showed that the innate immunity of red pandas increases significantly in geriatric individuals, whereas its change remains unclear when comparing adults and juveniles. We found that the adaptive immunity of red pandas first increased and then decreased with age. We identified CXCR3, BLNK, and CCR4 as the hub genes in the age-related protein-protein interaction network, which showed their central role in age-related immune changes. Many DNA repair genes were down-regulated in geriatric red pandas, suggesting that the DNA repair ability of the blood tissue in geriatric red pandas is significantly reduced. The significantly up-regulated TLR5 in geriatric individuals also suggests the possibility of enhancing the vaccination immune response by incorporating flagellin, which could be used to address decreased vaccine responses caused by age-related declines in immune system function. This work provides an insight into gene expression changes associated with aging and paves the way for effective disease prevention and treatment strategies for red pandas in the future.

A study of strong nucleosomes in the human genome

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, nucleosomes are highly dynamic and susceptible to experimental conditions, resulting in extreme variability across nucleosome maps, which complicates the generation of accurate nucleosome organization data. We mapped nucleosomes from different individuals using improved MNase-seq. The improvements included setting different digestion levels (low, medium, high) and naked DNA correction to remove the noise caused by experimental manipulation and comparing maps to obtain the accurate position and occupancy of strong nucleosomes (SNs) in the whole genome. In addition, the characteristics of SNs were further excavated. SNs were enriched in Alu elements and near the centromere of Chr12. SNs contain some specific sequences, and the GC content of SNs is different from that of dynamic nucleosomes. The findings suggest that nucleosome location in the genome and the DNA sequence may affect nucleosome stability.

•

Naked DNA correction improved the accuracy of nucleosome map in partial digestion

•

Level of MNase digestion has effects on nucleosome organization

•

A type of strong nucleosomes (SNs) exist across different nucleosome maps

•

Nucleosome stability may be related to its location and the DNA sequence

The 22q11.2 region regulates presynaptic gene-products linked to schizophrenia

It is unclear how the 22q11.2 deletion predisposes to psychiatric disease. To study this, we generated induced pluripotent stem cells from deletion carriers and controls and utilized CRISPR/Cas9 to introduce the heterozygous deletion into a control cell line. Here, we show that upon differentiation into neural progenitor cells, the deletion acted in trans to alter the abundance of transcripts associated with risk for neurodevelopmental disorders including autism. In excitatory neurons, altered transcripts encoded presynaptic factors and were associated with genetic risk for schizophrenia, including common and rare variants. To understand how the deletion contributed to these changes, we defined the minimal protein-protein interaction network that best explains gene expression alterations. We found that many genes in 22q11.2 interact in presynaptic, proteasome, and JUN/FOS transcriptional pathways. Our findings suggest that the 22q11.2 deletion impacts genes that may converge with psychiatric risk loci to influence disease manifestation in each deletion carrier.

Remission of obesity and insulin resistance is not sufficient to restore mitochondrial homeostasis in visceral adipose tissue

Metabolic plasticity is the ability of a biological system to adapt its metabolic phenotype to different environmental stressors. We used a whole-body and tissue-specific phenotypic, functional, proteomic, metabolomic and transcriptomic approach to systematically assess metabolic plasticity in diet-induced obese mice after a combined nutritional and exercise intervention. Although most obesity and overnutrition-related pathological features were successfully reverted, we observed a high degree of metabolic dysfunction in visceral white adipose tissue, characterized by abnormal mitochondrial morphology and functionality. Despite two sequential therapeutic interventions and an apparent global healthy phenotype, obesity triggered a cascade of events in visceral adipose tissue progressing from mitochondrial metabolic and proteostatic alterations to widespread cellular stress, which compromises its biosynthetic and recycling capacity. In humans, weight loss after bariatric surgery showed a transcriptional signature in visceral adipose tissue similar to our mouse model of obesity reversion. Overall, our data indicate that obesity prompts a lasting metabolic fingerprint that leads to a progressive breakdown of metabolic plasticity in visceral adipose tissue.

Spatiotemporal Regulation of a Single Adaptively Evolving Trans-Regulatory Element Contributes to Spermatogenetic Expression Divergence in Drosophila

Due to extensive pleiotropy, trans-acting elements are often thought to be evolutionarily constrained. While the impact of trans-acting elements on gene expression evolution has been extensively studied, relatively little is understood about the contribution of a single trans regulator to interspecific expression and phenotypic divergence. Here, we disentangle the effects of genomic context and miR-983, an adaptively evolving young microRNA, on expression divergence between Drosophila melanogaster and D. simulans. We show miR-983 effects promote interspecific expression divergence in testis despite its antagonism with the often-predominant context effects. Single-cyst RNA-seq reveals that distinct sets of genes gain and lose miR-983 influence under disruptive or diversifying selection at different stages of spermatogenesis, potentially helping minimize antagonistic pleiotropy. At the round spermatid stage, the effects of miR-983 are weak and distributed, coincident with the transcriptome undergoing drastic expression changes. Knocking out miR-983 causes reduced sperm length with increased within-individual variation in D. melanogaster but not in D. simulans, and the D. melanogaster knockout also exhibits compromised sperm defense ability. Our results provide empirical evidence for the resolution of antagonistic pleiotropy and also have broad implications for the function and evolution of new trans regulators.

Emerging role of a systems biology approach to elucidate factors of reduced penetrance: transcriptional changes in THAP1-linked dystonia as an example

Pathogenic variants in THAP1 can cause dystonia with a penetrance of about 50 %. The underlying mechanisms are unknown and can be considered as means of endogenous disease protection. Since THAP1 encodes a transcription factor, drivers of this variability putatively act at the transcriptome level. Several transcriptome studies tried to elucidate THAP1 function in diverse cellular and mouse models, including mutation carrier-derived cells and iPSC-derived neurons, unveiling various differentially expressed genes and affected pathways. These include nervous system development, dopamine signalling, myelination, or cell-cell adhesion. A network diffusion analysis revealed mRNA splicing, mitochondria, DNA repair, and metabolism as significant pathways that may represent potential targets for therapeutic interventions.

Lymphocyte crosstalk is required for monocyte-intrinsic trained immunity to Plasmodium falciparum

Plasmodium falciparum (P. falciparum) induces trained innate immune responses in vitro, where initial stimulation of adherent PBMCs with P. falciparum-infected RBCs (iRBCs) results in hyperresponsiveness to subsequent ligation of TLR2. This response correlates with the presence of T and B lymphocytes in adherent PBMCs, suggesting that innate immune training is partially due to adaptive immunity. We found that T cell-depleted PBMCs and purified monocytes alone did not elicit hyperproduction of IL-6 and TNF-α under training conditions. Analysis of P. falciparum-trained PBMCs showed that DCs did not develop under control conditions, and IL-6 and TNF-α were primarily produced by monocytes and DCs. Transwell experiments isolating purified monocytes from either PBMCs or purified CD4+ T cells, but allowing diffusion of secreted proteins, enabled monocytes trained with iRBCs to hyperproduce IL-6 and TNF-α after TLR restimulation. Purified monocytes stimulated with IFN-γ hyperproduced IL-6 and TNF-α, whereas blockade of IFN-γ in P. falciparum-trained PBMCs inhibited trained responses. Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-Seq) on monocytes from patients with malaria showed persistently open chromatin at genes that appeared to be trained in vitro. Together, these findings indicate that the trained immune response of monocytes to P. falciparum is not completely cell intrinsic but depends on soluble signals from lymphocytes.

Differential cofactor dependencies define distinct types of human enhancers

All multicellular organisms rely on differential gene transcription regulated by genomic enhancers, which function through cofactors that are recruited by transcription factors^1,2. Emerging evidence suggests that not all cofactors are required at all enhancers^3-5, yet whether these observations reflect more general principles or distinct types of enhancers remained unknown. Here we categorized human enhancers by their cofactor dependencies and show that these categories provide a framework to understand the sequence and chromatin diversity of enhancers and their roles in different gene-regulatory programmes. We quantified enhancer activities along the entire human genome using STARR-seq⁶ in HCT116 cells, following the rapid degradation of eight cofactors. This analysis identified different types of enhancers with distinct cofactor requirements, sequences and chromatin properties. Some enhancers were insensitive to the depletion of the core Mediator subunit MED14 or the bromodomain protein BRD4 and regulated distinct transcriptional programmes. In particular, canonical Mediator⁷ seemed dispensable for P53-responsive enhancers, and MED14-depleted cells induced endogenous P53 target genes. Similarly, BRD4 was not required for the transcription of genes that bear CCAAT boxes and a TATA box (including histone genes and LTR12 retrotransposons) or for the induction of heat-shock genes. This categorization of enhancers through cofactor dependencies reveals distinct enhancer types that can bypass broadly utilized cofactors, which illustrates how alternative ways to activate transcription separate gene expression programmes and provide a conceptual framework to understand enhancer function and regulatory specificity.

Time-series transcriptomic screening of factors contributing to the cross-tolerance to UV radiation and anhydrobiosis in tardigrades

BACKGROUND

Tardigrades are microscopic animals that are capable of tolerating extreme environments by entering a desiccated state of suspended animation known as anhydrobiosis. While antioxidative stress proteins, antiapoptotic pathways and tardigrade-specific intrinsically disordered proteins have been implicated in the anhydrobiotic machinery, conservation of these mechanisms is not universal within the phylum Tardigrada, suggesting the existence of overlooked components.

RESULTS

Here, we show that a novel Mn-dependent peroxidase is an important factor in tardigrade anhydrobiosis. Through time-series transcriptome analysis of Ramazzottius varieornatus specimens exposed to ultraviolet light and comparison with anhydrobiosis entry, we first identified several novel gene families without similarity to existing sequences that are induced rapidly after stress exposure. Among these, a single gene family with multiple orthologs that is highly conserved within the phylum Tardigrada and enhances oxidative stress tolerance when expressed in human cells was identified. Crystallographic study of this protein suggested Zn or Mn binding at the active site, and we further confirmed that this protein has Mn-dependent peroxidase activity in vitro.

CONCLUSIONS

Our results demonstrated novel mechanisms for coping with oxidative stress that may be a fundamental mechanism of anhydrobiosis in tardigrades. Furthermore, localization of these sets of proteins mainly in the Golgi apparatus suggests an indispensable role of the Golgi stress response in desiccation tolerance.

Metapone: a Bioconductor package for joint pathway testing for untargeted metabolomics data

MOTIVATION

Testing for pathway enrichment is an important aspect in the analysis of untargeted metabolomics data. Due to the unique characteristics of untargeted metabolomics data, some key issues have not been fully addressed in existing pathway testing algorithms: (1) matching uncertainty between data features and metabolites; (2) lacking of method to analyze positive mode and negative mode LC/MS data simultaneously on the same set of subjects; (3) the incompleteness of pathways in individual software packages.

RESULTS

We developed an innovative R/Bioconductor package: metabolic pathway testing with positive and negative mode data (metapone), which can perform two novel statistical tests that take matching uncertainty into consideration - (1) a weighted GSEA-type test, and (2) a permutation-based weighted hypergeometric test. The package is capable of combining positive and negative ion mode results in a single testing scheme. For comprehensiveness, the built-in pathways were manually curated from three sources: KEGG, Mummichog, and SMPDB.

AVAILABILITY

The package is available at https://bioconductor.org/packages/devel/bioc/html/metapone.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Neuronal Cells Display Distinct Stability Controls of Alternative Polyadenylation mRNA Isoforms, Long Non-Coding RNAs, and Mitochondrial RNAs

RNA stability plays an important role in gene expression. Here, using 3' end sequencing of newly made and pre-existing poly(A)+ RNAs, we compare transcript stability in multiple human cell lines, including HEK293T, HepG2, and SH-SY5Y. We show that while mRNA stability is generally conserved across the cell lines, specific transcripts having a high GC content and possibly more stable secondary RNA structures are relatively more stable in SH-SY5Y cells compared to the other 2 cell lines. These features also differentiate stability levels of alternative polyadenylation (APA) 3'UTR isoforms in a cell type-specific manner. Using differentiation of a neural stem cell line as a model, we show that mRNA stability difference could contribute to gene expression changes in neurogenesis and confirm the neuronal identity of SH-SY5Y cells at both gene expression and APA levels. In addition, compared to transcripts using 3'-most exon cleavage/polyadenylation sites (PASs), those using intronic PASs are generally less stable, especially when the PAS-containing intron is large and has a strong 5' splice site, suggesting that intronic polyadenylation mostly plays a negative role in gene expression. Interestingly, the differential mRNA stability among APA isoforms appears to buffer PAS choice in these cell lines. Moreover, we found that several other poly(A)+ RNA species, including promoter-associated long noncoding RNAs and transcripts encoded by the mitochondrial genome, are more stable in SH-SY5Y cells than the other 2 cell lines, further highlighting distinct RNA metabolism in neuronal cells. Together, our results indicate that distinct RNA stability control in neuronal cells may contribute to the gene expression and APA programs that define their cell identity.

Population-Scale Polymorphic Short Tandem Repeat Provides an Alternative Strategy for Allele Mining in Cotton

Short tandem repeats (STRs), which vary in size due to featuring variable numbers of repeat units, are present throughout most eukaryotic genomes. To date, few population-scale studies identifying STRs have been reported for crops. Here, we constructed a high-density polymorphic STR map by investigating polymorphic STRs from 911 Gossypium hirsutum accessions. In total, we identified 556,426 polymorphic STRs with an average length of 21.1 bp, of which 69.08% were biallelic. Moreover, 7,718 (1.39%) were identified in the exons of 6,021 genes, which were significantly enriched in transcription, ribosome biogenesis, and signal transduction. Only 5.88% of those exonic STRs altered open reading frames, of which 97.16% were trinucleotide. An alternative strategy STR-GWAS analysis revealed that 824 STRs were significantly associated with agronomic traits, including 491 novel alleles that undetectable by previous SNP-GWAS methods. For instance, a novel polymorphic STR consisting of GAACCA repeats was identified in GH_D06G1697, with its (GAACCA)₅ allele increasing fiber length by 1.96-4.83% relative to the (GAACCA)₄ allele. The database CottonSTRDB was further developed to facilitate use of STR datasets in breeding programs. Our study provides functional roles for STRs in influencing complex traits, an alternative strategy STR-GWAS for allele mining, and a database serving the cotton community as a valuable resource.

The First Genome of the Balearic Shearwater (Puffinus mauretanicus) Provides a Valuable Resource for Conservation Genomics and Sheds Light on Adaptation to a Pelagic lifestyle

The Balearic shearwater (Puffinus mauretanicus) is the most threatened seabird in Europe and a member of the most speciose group of pelagic seabirds, the order Procellariiformes, which exhibit extreme adaptations to a pelagic lifestyle. The fossil record suggests that human colonisation of the Balearic Islands resulted in a sharp decrease of the Balearic shearwater population size. Currently, populations of the species continue to be decimated mainly due to predation by introduced mammals and bycatch in longline fisheries, with some studies predicting its extinction by 2070. Here, using a combination of short and long reads, we generate the first high-quality reference genome for the Balearic shearwater, with a completeness amongst the highest across available avian species. We used this reference genome to study critical aspects relevant to the conservation status of the species and to gain insights into the adaptation to a pelagic lifestyle of the order Procellariiformes. We detected relatively high levels of genome-wide heterozygosity in the Balearic shearwater despite its reduced population size. However, the reconstruction of its historical demography uncovered an abrupt population decline potentially linked to a reduction of the neritic zone during the Penultimate Glacial Period (∼194-135 ka). Comparative genomics analyses uncover a set of candidate genes that may have played an important role into the adaptation to a pelagic lifestyle of Procellariiformes, including those for the enhancement of fishing capabilities, night vision, and the development of natriuresis. The reference genome obtained will be the crucial in the future development of genetic tools in conservation efforts for this Critically Endangered species.

Conservation and architecture of housekeeping genes in the model marine diatom Thalassiosira pseudonana

Housekeeping genes (HKGs) are constitutively expressed with low variation across tissues/conditions. They are thought to be highly conserved and fundamental to cellular maintenance, with distinctive genomic features. Here, we identify 1505 HKGs in the unicellular marine diatom Thalassiosira pseudonana based on an RNA-seq analysis of 232 samples taken under 12 experimental conditions over 0-72 h. We identify promising internal reference genes (IRGs) for T. pseudonana from the most stably expressed HKGs. A comparative analysis indicates < 18% of HKGs in T. pseudonana have orthologs in other eukaryotes, including other diatom species. Contrary to work on human tissues, T. pseudonana HKGs are longer than non-HKGs, due to elongated introns. More ancient HKGs tend to be shorter than more recent HKGs, and expression levels of HKGs decrease more rapidly with gene length relative to non-HKGs. Our results indicate that HKGs are highly variable across the tree of life and thus unlikely to be universally fundamental for cellular maintenance. We hypothesize that the distinct genomic features of HKGs of T. pseudonana may be a consequence of selection pressures associated with high expression and low variance across conditions.

Dissociation of nanosilicates induces downstream endochondral differentiation gene expression program

Bioactive materials harness the body's innate regenerative potential by directing endogenous progenitor cells to facilitate tissue repair. Dissolution products of inorganic biomaterials provide unique biomolecular signaling for tissue-specific differentiation. Inorganic ions (minerals) are vital to biological processes and play crucial roles in regulating gene expression patterns and directing cellular fate. However, mechanisms by which ionic dissolution products affect cellular differentiation are not well characterized. We demonstrate the role of the inorganic biomaterial synthetic two-dimensional nanosilicates and its ionic dissolution products on human mesenchymal stem cell differentiation. We use whole-transcriptome sequencing (RNA-sequencing) to characterize the contribution of nanosilicates and its ionic dissolution products on endochondral differentiation. Our study highlights the modulatory role of ions in stem cell transcriptome dynamics by regulating lineage-specific gene expression patterns. This work paves the way for leveraging biochemical characteristics of inorganic biomaterials to direct cellular processes and promote in situ tissue regeneration.

RdsA Is a Global Regulator That Controls Cell Shape and Division in Rhizobium etli

Unlike other bacteria, cell growth in rhizobiales is unipolar and asymmetric. The regulation of cell division, and its coordination with metabolic processes is an active field of research. In Rhizobium etli, gene RHE_PE00024, located in a secondary chromosome, is essential for growth. This gene encodes a predicted hybrid histidine kinase sensor protein, participating in a, as yet undescribed, two-component signaling system. In this work, we show that a conditional knockdown mutant (cKD24) in RHE_PE00024 (hereby referred as rdsA, after rhizobium division and shape) generates a striking phenotype, where nearly 64% of the cells present a round shape, with stochastic and uncoordinated cell division. For rod-shaped cells, a large fraction (12 to 29%, depending on their origin) present growth from the old pole, a sector that is normally inactive for growth in a wild-type cell. A fraction of the cells (1 to 3%) showed also multiple ectopic polar growths. Homodimerization of RdsA appears to be required for normal function. RNAseq analysis of mutant cKD24 reveals global changes, with downregulated genes in at least five biological processes: cell division, wall biogenesis, respiration, translation, and motility. These modifications may affect proper structuring of the divisome, as well as peptidoglycan synthesis. Together, these results indicate that the hybrid histidine kinase RdsA is an essential global regulator influencing cell division and cell shape in R. etli.

Transcriptomic Analysis of Fumarate Compounds Identifies Unique Effects of Isosorbide Di-(Methyl Fumarate) on NRF2, NF-kappaB and IRF1 Pathway Genes

Dimethyl fumarate (DMF) has emerged as a first-line therapy for relapsing-remitting multiple sclerosis (RRMS). This treatment, however, has been limited by adverse effects, which has prompted development of novel derivatives with improved tolerability. We compared the effects of fumarates on gene expression in astrocytes. Our analysis included diroximel fumarate (DRF) and its metabolite monomethyl fumarate (MMF), along with a novel compound isosorbide di-(methyl fumarate) (IDMF). Treatment with IDMF resulted in the largest number of differentially expressed genes. The effects of DRF and MMF were consistent with NRF2 activation and NF-κB inhibition, respectively. IDMF responses, however, were concordant with both NRF2 activation and NF-κB inhibition, and we confirmed IDMF-mediated NF-κB inhibition using a reporter assay. IDMF also down-regulated IRF1 expression and IDMF-decreased gene promoters were enriched with IRF1 recognition sequences. Genes altered by each fumarate overlapped significantly with those near loci from MS genetic association studies, but IDMF had the strongest overall effect on MS-associated genes. These results show that next-generation fumarates, such as DRF and IDMF, have effects differing from those of the MMF metabolite. Our findings support a model in which IDMF attenuates oxidative stress via NRF2 activation, with suppression of NF-κB and IRF1 contributing to mitigation of inflammation and pyroptosis.

The bud dormancy disconnect: latent buds of grapevine are dormant during summer despite a high metabolic rate

Grapevine (Vitis vinifera L.) displays wide plasticity to climate; however, the physiology of dormancy along a seasonal continuum is poorly understood. Here we investigated the apparent disconnect between dormancy and the underlying respiratory physiology and transcriptome of grapevine buds, from bud set in summer to bud burst in spring. The establishment of dormancy in summer was pronounced and reproducible; however, this was coupled with little or no change in physiology, indicated by respiration, hydration, and tissue oxygen tension. The release of dormancy was biphasic; the depth of dormancy declined substantially by mid-autumn, while the subsequent decline towards spring was moderate. Observed changes in physiology failed to explain the first phase of dormancy decline, in particular. Transcriptome data contrasting development from summer through to spring also indicated that dormancy was poorly reflected by metabolic quiescence during summer and autumn. Gene Ontology and enrichment data revealed the prevailing influence of abscisic acid (ABA)-related gene expression during the transition from summer to autumn, and promoter motif analysis suggested that photoperiod may play an important role in regulating ABA functions during the establishment of dormancy. Transcriptomic data from later transitions reinforced the importance of oxidation and hypoxia as physiological cues to regulate the maintenance of quiescence and resumption of growth. Collectively these data reveal a novel disconnect between growth and metabolic quiescence in grapevine following bud set, which requires further experimentation to explain the phenology and dormancy relationships.