Clust: automatic extraction of optimal co-expressed gene clusters from gene expression data

3D clustering of gene expression data from systemic autoinflammatory diseases using self-organizing maps (Clust3D)

Background and objective

Systemic autoinflammatory diseases (SAIDs) are characterized by widespread inflammation, but for most of them there is a lack of specific biomarkers for accurate diagnosis. Although a number of machine learning algorithms have been used to analyze SAID datasets, aiding in the discovery of novel biomarkers, there is a growing recognition of the importance of SAID timeseries clustering, as it can capture the temporal dynamics of gene expression patterns.

Methodology

This paper proposes a novel clustering methodology to efficiently associate three-dimensional data. The algorithm utilizes competitive learning to create a self-organizing neural network and adjust neuron positions in time-dependent and high dimensional feature space in order to assign them as clustering centers. The quantitative evaluation of the clustering was based on well-known clustering indices. Furthermore, a differential expression analysis and classification pipeline was employed to assess the capability of the proposed methodology to extract more accurate pathway-specific genes from its clusters. For that, a comparative analysis was also conducted against a heuristic timeseries clustering method.

Results

The proposed methodology achieved better overall clustering indices scores and classification metrics using genes derived from its clusters. Notable cases include a threefold increase in the Calinski-Harabasz clustering index, a twofold improvement in the Davies-Bouldin clustering index and a ∼ 60 % increase in the classification specificity score.

Conclusion

A novel clustering methodology was developed and applied on several gene expression timeseries datasets from systemic autoinflammatory diseases, and its ability to efficiently produce well separated clusters compared to existing heuristic methods was demonstrated.

Photoperiod-1 regulates the wheat inflorescence transcriptome to influence spikelet architecture and flowering time

Photoperiod insensitivity has been selected by breeders to help adapt crops to diverse environments and farming practices. In wheat, insensitive alleles of Photoperiod-1 (Ppd-1) relieve the requirement of long daylengths to flower by promoting expression of floral promoting genes early in the season; however, these alleles also limit yield by reducing the number and fertility of grain-producing florets through processes that are poorly understood. Here, we performed transcriptome analysis of the developing inflorescence using near-isogenic lines that contain either photoperiod-insensitive or null alleles of Ppd-1, during stages when spikelet number is determined and floret development initiates. We report that Ppd-1 influences the stage-specific expression of genes with roles in auxin signaling, meristem identity, and protein turnover, and analysis of differentially expressed transcripts identified bZIP and ALOG transcription factors, namely PDB1 and ALOG1, which regulate flowering time and spikelet architecture. These findings enhance our understanding of genes that regulate inflorescence development and introduce new targets for improving yield potential.

Exploring casual effects and shared molecular mechanism between psoriasis and liver cancer through Mendelian randomization and comprehensive bioinformatic analyses

Psoriasis (Ps), a chronic inflammatory disease affecting approximately 2 % of the global population, has been associated with an increased risk of liver cancer in observational studies. However, their causal relationships as well as underlying shared molecular mechanisms between Ps and liver cancer remain unclear. Using bidirectional Mendelian randomization analysis, we revealed that a genetic predisposition to liver cancer increased the risk of Ps in European and East Asian populations but not the other way around. Moreover, we analyzed three transcriptomic datasets of patients with Ps and liver cancer from open-source databases. Differentially expressed genes (DEGs) and disease-specific gene co-expression module analyses revealed that cell-cycle dysregulation was the shared mechanism of Ps and liver cancer. Moreover, we identified a rank-conservative gene signature shared between these two diseases, which demonstrated significance in diagnostic and prognostic predictions. These findings provided valuable insights into the interconnections between Ps and liver cancer, which may be helpful to guide therapeutic management.

Valproic Acid Causes Redox-Regulated Post-Translational Protein Modifications That Are Dependent upon P19 Cellular Differentiation States

Valproic acid (VPA) is a common anti-epileptic drug and known neurodevelopmental toxicant. Although the exact mechanism of VPA toxicity remains unknown, recent findings show that VPA disrupts redox signaling in undifferentiated cells but has little effect on fully differentiated neurons. Redox imbalances often alter oxidative post-translational protein modifications and could affect embryogenesis if developmentally critical proteins are targeted. We hypothesize that VPA causes redox-sensitive post-translational protein modifications that are dependent upon cellular differentiation states. Undifferentiated P19 cells and P19-derived neurons were treated with VPA alone or pretreated with D3T, an inducer of the nuclear factor erythroid 2-related factor 2 (NRF2) antioxidant pathway, prior to VPA exposure. Undifferentiated cells treated with VPA alone exhibited an oxidized glutathione redox couple and increased overall protein oxidation, whereas differentiated neurons were protected from protein oxidation via increased S-glutathionylation. Pretreatment with D3T prevented the effects of VPA exposure in undifferentiated cells. Taken together, our findings support redox-sensitive post-translational protein alterations in undifferentiated cells as a mechanism of VPA-induced developmental toxicity and propose NRF2 activation as a means to preserve proper neurogenesis.

Investigating the potential of Juglans regia phytoconstituents for the treatment of cervical cancer utilizing network biology and molecular docking approach

The fourth most frequent type of cancer in women and the leading cause of mortality for females worldwide is cervical cancer. Traditionally, medicinal plants have been utilized to treat various illnesses and ailments. The molecular docking method is used in the current study to look into the phytoconstituents of Juglans regia's possible anticancer effects on cervical cancer target proteins. This work uses the microarray dataset analysis of GSE63678 from the NCBI Gene Expression Omnibus database to find differentially expressed genes. Furthermore, protein-protein interactions of differentially expressed genes were constructed using network biology techniques. The top five hub genes (IGF1, FGF2, ESR1, MYL9, and MYH11) are then determined by computing topological parameters with Cytohubba. In addition, molecular docking research was performed on Juglans regia phytocompounds that were extracted from the IMPPAT database versus hub genes that had been identified. Utilizing molecular dynamics, simulation confirmed that prioritized docked complexes with low binding energies were stable.

Nonlinear DNA methylation trajectories in aging male mice

Although DNA methylation data yields highly accurate age predictors, little is known about the dynamics of this quintessential epigenomic biomarker during lifespan. To narrow the gap, we investigate the methylation trajectories of male mouse colon at five different time points of aging. Our study indicates the existence of sudden hypermethylation events at specific stages of life. Precisely, we identify two epigenomic switches during early-to-midlife (3-9 months) and mid-to-late-life (15-24 months) transitions, separating the rodents' life into three stages. These nonlinear methylation dynamics predominantly affect genes associated with the nervous system and enrich in bivalently marked chromatin regions. Based on groups of nonlinearly modified loci, we construct a clock-like classifier STageR (STage of aging estimatoR) that accurately predicts murine epigenetic stage. We demonstrate the universality of our clock in an independent mouse cohort and with publicly available datasets.

Integrated proteogenomic characterization of glioblastoma evolution

The evolutionary trajectory of glioblastoma (GBM) is a multifaceted biological process that extends beyond genetic alterations alone. Here, we perform an integrative proteogenomic analysis of 123 longitudinal glioblastoma pairs and identify a highly proliferative cellular state at diagnosis and replacement by activation of neuronal transition and synaptogenic pathways in recurrent tumors. Proteomic and phosphoproteomic analyses reveal that the molecular transition to neuronal state at recurrence is marked by post-translational activation of the wingless-related integration site (WNT)/ planar cell polarity (PCP) signaling pathway and BRAF protein kinase. Consistently, multi-omic analysis of patient-derived xenograft (PDX) models mirror similar patterns of evolutionary trajectory. Inhibition of B-raf proto-oncogene (BRAF) kinase impairs both neuronal transition and migration capability of recurrent tumor cells, phenotypic hallmarks of post-therapy progression. Combinatorial treatment of temozolomide (TMZ) with BRAF inhibitor, vemurafenib, significantly extends the survival of PDX models. This study provides comprehensive insights into the biological mechanisms of glioblastoma evolution and treatment resistance, highlighting promising therapeutic strategies for clinical intervention.

CAP-RNAseq: an integrated pipeline for functional annotation and prioritization of co-expression clusters

Cluster analysis is one of the most widely used exploratory methods for visualization and grouping of gene expression patterns across multiple samples or treatment groups. Although several existing online tools can annotate clusters with functional terms, there is no all-in-one webserver to effectively prioritize genes/clusters using gene essentiality as well as congruency of mRNA-protein expression. Hence, we developed CAP-RNAseq that makes possible (1) upload and clustering of bulk RNA-seq data followed by identification, annotation and network visualization of all or selected clusters; and (2) prioritization using DepMap gene essentiality and/or dependency scores as well as the degree of correlation between mRNA and protein levels of genes within an expression cluster. In addition, CAP-RNAseq has an integrated primer design tool for the prioritized genes. Herein, we showed using comparisons with the existing tools and multiple case studies that CAP-RNAseq can uniquely aid in the discovery of co-expression clusters enriched with essential genes and prioritization of novel biomarker genes that exhibit high correlations between their mRNA and protein expression levels. CAP-RNAseq is applicable to RNA-seq data from different contexts including cancer and available at http://konulabapps.bilkent.edu.tr:3838/CAPRNAseq/ and the docker image is downloadable from https://hub.docker.com/r/konulab/caprnaseq.

Multi-omic approaches for host-microbiome data integration

The gut microbiome interacts with the host through complex networks that affect physiology and health outcomes. It is becoming clear that these interactions can be measured across many different omics layers, including the genome, transcriptome, epigenome, metabolome, and proteome, among others. Multi-omic studies of the microbiome can provide insight into the mechanisms underlying host-microbe interactions. As more omics layers are considered, increasingly sophisticated statistical methods are required to integrate them. In this review, we provide an overview of approaches currently used to characterize multi-omic interactions between host and microbiome data. While a large number of studies have generated a deeper understanding of host-microbiome interactions, there is still a need for standardization across approaches. Furthermore, microbiome studies would also benefit from the collection and curation of large, publicly available multi-omics datasets.

Identification of driving genes of familial adenomatous polyposis by differential gene expression analysis and weighted gene co-expression network analysis

BACKGROUND

Despite the advancement of new screening strategies and the advances in pharmacological therapies, the cancerization rates of familial adenomatous polyposis (FAP) are stable and even increased in the last years. Therefore, it necessitates additional research to characterize and understand the underlying mechanisms of FAP.

OBJECTIVE

To determine the genes that drive the pathogenesis of familial adenomatous polyposis (FAP).

METHODS

We performed on a cohort (GSE111156) gene profile, which consist of four group of gene expressions (the gene expressions of cancer, adenoma and normal tissue of duodenal cancer from patients with FAP were defined as Case N, Case A and Case C respectively, while that of adenoma tissue from patients with FAP who did not have duodenal cancer was Ctrl A). Tracking Tumor Immunophenotype (TIP) website was applied to reveal immune infiltration profile and signature genes of FAP. We merged the genes of key module (pink and midnight module) with signature genes to obtained the biomarkers related with FAP pathogenesis. The expression of these five biomarkers in FAP intratumoral region (IT) and tumor rim (TR) was detected with Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR).

RESULTS

In total, 220, 23 and 63 DEGs were determined in Cases C, A and N, in comparison to Ctrl A. In total, 196 and 10 DEGs were determined in Cases C and A, separately, as compared to Case N. A total of four biomarkers including CCL5, CD3G, CD2 and TLR3 were finally identified associated with pink module, while only one biomarker (KLF2) associated with midnight module was identified. All biomarkers were evidently raised in FAP IT tissues utilizing qRT-PCR.

CONCLUSION

We identified five potential biomarkers for pathogenesis of FAP to understand the fundamental mechanisms of FAP progression and revealed some probable targets for the diagnosis or treatment of FAP.

Genome-wide analysis of long noncoding RNAs in response to salt stress in Nicotiana tabacum

BACKGROUND

Long noncoding RNAs (lncRNAs) have been shown to play important roles in the response of plants to various abiotic stresses, including drought, heat and salt stress. However, the identification and characterization of genome-wide salt-responsive lncRNAs in tobacco (Nicotiana tabacum L.) have been limited. Therefore, this study aimed to identify tobacco lncRNAs in roots and leaves in response to different durations of salt stress treatment.

RESULTS

A total of 5,831 lncRNAs were discovered, with 2,428 classified as differentially expressed lncRNAs (DElncRNAs) in response to salt stress. Among these, only 214 DElncRNAs were shared between the 2,147 DElncRNAs in roots and the 495 DElncRNAs in leaves. KEGG pathway enrichment analysis revealed that these DElncRNAs were primarily associated with pathways involved in starch and sucrose metabolism in roots and cysteine and methionine metabolism pathway in leaves. Furthermore, weighted gene co-expression network analysis (WGCNA) identified 15 co-expression modules, with four modules strongly linked to salt stress across different treatment durations (MEsalmon, MElightgreen, MEgreenyellow and MEdarkred). Additionally, an lncRNA-miRNA-mRNA network was constructed, incorporating several known salt-associated miRNAs such as miR156, miR169 and miR396.

CONCLUSIONS

This study enhances our understanding of the role of lncRNAs in the response of tobacco to salt stress. It provides valuable information on co-expression networks of lncRNA and mRNAs, as well as networks of lncRNAs-miRNAs-mRNAs. These findings identify important candidate lncRNAs that warrant further investigation in the study of plant-environment interactions.

Transcriptome analysis reveals temporally regulated genetic networks during Drosophila border cell collective migration

BACKGROUND

Collective cell migration underlies many essential processes, including sculpting organs during embryogenesis, wound healing in the adult, and metastasis of cancer cells. At mid-oogenesis, Drosophila border cells undergo collective migration. Border cells round up into a small group at the pre-migration stage, detach from the epithelium and undergo a dynamic and highly regulated migration at the mid-migration stage, and stop at the oocyte, their final destination, at the post-migration stage. While specific genes that promote cell signaling, polarization of the cluster, formation of protrusions, and cell-cell adhesion are known to regulate border cell migration, there may be additional genes that promote these distinct active phases of border cell migration. Therefore, we sought to identify genes whose expression patterns changed during border cell migration.

RESULTS

We performed RNA-sequencing on border cells isolated at pre-, mid-, and post-migration stages. We report that 1,729 transcripts, in nine co-expression gene clusters, are temporally and differentially expressed across the three migration stages. Gene ontology analyses and constructed protein-protein interaction networks identified genes expected to function in collective migration, such as regulators of the cytoskeleton, adhesion, and tissue morphogenesis, but also uncovered a notable enrichment of genes involved in immune signaling, ribosome biogenesis, and stress responses. Finally, we validated the in vivo expression and function of a subset of identified genes in border cells.

CONCLUSIONS

Overall, our results identified differentially and temporally expressed genetic networks that may facilitate the efficient development and migration of border cells. The genes identified here represent a wealth of new candidates to investigate the molecular nature of dynamic collective cell migrations in developing tissues.

Maternal blood transcriptome as a sensor of fetal organ maturation at the end of organogenesis in cattle†

Harnessing information from the maternal blood to predict fetal growth is attractive yet scarcely explored in livestock. The objectives were to determine the transcriptomic modifications in maternal blood and fetal liver, gonads, and heart according to fetal weight and to model a molecular signature based on the fetal organs allowing the prediction of fetal weight from the maternal blood transcriptome in cattle. In addition to a contemporaneous maternal blood sample, organ samples were collected from 10 male fetuses at 42 days of gestation for RNA-sequencing. Fetal weight ranged from 1.25 to 1.69 g (mean = 1.44 ± 0.15 g). Clustering data analysis revealed clusters of co-expressed genes positively correlated with fetal weight and enriching ontological terms biologically relevant for the organ. For the heart, the 1346 co-expressed genes were involved in energy generation and protein synthesis. For the gonads, the 1042 co-expressed genes enriched seminiferous tubule development. The 459 co-expressed genes identified in the liver were associated with lipid synthesis and metabolism. Finally, the cluster of 571 co-expressed genes determined in maternal blood enriched oxidative phosphorylation and thermogenesis. Next, data from the fetal organs were used to train a regression model of fetal weight, which was predicted with the maternal blood data. The best prediction was achieved when the model was trained with 35 co-expressed genes overlapping between heart and maternal blood (root-mean-square error = 0.04, R2 = 0.93). In conclusion, linking transcriptomic information from maternal blood with that from the fetal heart unveiled maternal blood as a predictor of fetal development.

Key processes required for the different stages of fungal carnivory by a nematode-trapping fungus

Nutritional deprivation triggers a switch from a saprotrophic to predatory lifestyle in soil-dwelling nematode-trapping fungi (NTF). In particular, the NTF Arthrobotrys oligospora secretes food and sex cues to lure nematodes to its mycelium and is triggered to develop specialized trapping devices. Captured nematodes are then invaded and digested by the fungus, thus serving as a food source. In this study, we examined the transcriptomic response of A. oligospora across the stages of sensing, trap development, and digestion upon exposure to the model nematode Caenorhabditis elegans. A. oligospora enacts a dynamic transcriptomic response, especially of protein secretion-related genes, in the presence of prey. Two-thirds of the predicted secretome of A. oligospora was up-regulated in the presence of C. elegans at all time points examined, and among these secreted proteins, 38.5% are predicted to be effector proteins. Furthermore, functional studies disrupting the t-SNARE protein Sso2 resulted in impaired ability to capture nematodes. Additionally, genes of the DUF3129 family, which are expanded in the genomes of several NTF, were highly up-regulated upon nematode exposure. We observed the accumulation of highly expressed DUF3129 proteins in trap cells, leading us to name members of this gene family as Trap Enriched Proteins (TEPs). Gene deletion of the most highly expressed TEP gene, TEP1, impairs the function of traps and prevents the fungus from capturing prey efficiently. In late stages of predation, we observed up-regulation of a variety of proteases, including metalloproteases. Following penetration of nematodes, these metalloproteases facilitate hyphal growth required for colonization of prey. These findings provide insights into the biology of the predatory lifestyle switch in a carnivorous fungus and provide frameworks for other fungal-nematode predator-prey systems.

Early-life stress and ovarian hormones alter transcriptional regulation in the nucleus accumbens resulting in sex-specific responses to cocaine

Early-life stress and ovarian hormones contribute to increased female vulnerability to cocaine addiction. Here, we reveal molecular substrates in the reward area, the nucleus accumbens, through which these female-specific factors affect immediate and conditioning responses to cocaine. We find shared involvement of X chromosome inactivation-related and estrogen signaling-related gene regulation in enhanced conditioning responses following early-life stress and during the low-estrogenic state in females. Low-estrogenic females respond to acute cocaine by opening neuronal chromatin enriched for the sites of ΔFosB, a transcription factor implicated in chronic cocaine response and addiction. Conversely, high-estrogenic females respond to cocaine by preferential chromatin closing, providing a mechanism for limiting cocaine-driven chromatin and synaptic plasticity. We find that physiological estrogen withdrawal, early-life stress, and absence of one X chromosome all nullify the protective effect of a high-estrogenic state on cocaine conditioning in females. Our findings offer a molecular framework to enable understanding of sex-specific neuronal mechanisms underlying cocaine use disorder.

Unravelling molecular mechanisms involved in resistance priming against downy mildew (Plasmopara viticola) in grapevine (Vitis vinifera L.)

Downy mildew (DM; Plasmopara viticola) is amongst the most severe fungal diseases in viticulture and the reason for the majority of fungicide applications. To reduce synthetic and copper-based fungicides, there is an urgent need for natural alternatives, which are being increasingly tested by the industry and the research community. However, their mode of action remains unclear. Therefore, our study aimed to investigate the transcriptomic changes induced by oregano essential oil vapour (OEOV) in DM-infected grapevines. OEOV was applied at different time points before and after DM infection to differentiate between a priming effect and a direct effect. Both pre-DM treatment with OEOV and post-infection treatment resulted in a significant reduction in DM sporulation. RNA-seq, followed by differential gene expression and weighted gene co-expression network analysis, identified co-expressed gene modules associated with secondary metabolism, pathogen recognition and response. Surprisingly, the molecular mechanisms underlying the efficiency of OEOV against DM appear to be independent of stilbene synthesis, and instead involve genes from a putative signalling pathway that has yet to be characterized. This study enhances our understanding of the molecular regulation of innate plant immunity and provides new insights into the mode of action of alternative natural antifungal agents.

Spatiotemporal protein dynamics during early organogenesis in mouse conceptuses treated with valproic acid

Valproic acid (VPA) is an anti-epileptic medication that increases the risk of neural tube defect (NTD) outcomes in infants exposed during gestation. Previous studies into VPA's mechanism of action have focused on alterations in gene expression and metabolism but have failed to consider how exposure changes the abundance of critical developmental proteins over time. This study evaluates the effects of VPA on protein abundance in the developmentally distinct tissues of the mouse visceral yolk sac (VYS) and embryo proper (EMB) using mouse whole embryo culture. Embryos were exposed to 600 μM VPA at 2 h intervals over 10 h during early organogenesis with the aim of identifying protein pathways relevant to VPA's mechanism of action in failed NTC. Protein abundance was measured through tandem mass tag (TMT) labeling followed by liquid chromatography and mass spectrometry. Overall, there were over 1500 proteins with altered abundance after VPA exposure in the EMB or VYS with 428 of these proteins showing previous gene expression associations with VPA exposure. Limited overlap of significant proteins between tissues supported the conclusion of independent roles for the VYS and EMB in response to VPA. Pathway analysis of proteins with increased or decreased abundance identified multiple pathways with mechanistic relevance to NTC and embryonic development including convergent extension, Wnt Signaling/planar cell polarity, cellular migration, cellular proliferation, cell death, and cytoskeletal organization processes as targets of VPA. Clustering of co-regulated proteins to identify shared patterns of protein abundance over time highlighted 4 h and 6/10 h as periods of divergent protein abundance between control and VPA-treated samples in the VYS and EMB, respectively. Overall, this study demonstrated that VPA temporally alters protein content in critical developmental pathways in the VYS and the EMB during early organogenesis in mice.

Spatiotemporal patterns of wheat response to Pyrenophora tritici-repentis in asymptomatic regions revealed by transcriptomic and X-ray fluorescence microscopy analyses

Pathogen attacks elicit dynamic and widespread molecular responses in plants. While our understanding of plant responses has advanced considerably, little is known of the molecular responses in the asymptomatic 'green' regions adjoining lesions. Here, we explore gene expression data and high-resolution elemental imaging to report the spatiotemporal changes in the asymptomatic green region of susceptible and moderately resistant wheat cultivars infected with a necrotrophic fungal pathogen, Pyrenophora tritici-repentis. We show, with improved spatiotemporal resolution, that calcium oscillations are modified in the susceptible cultivar, resulting in 'frozen' host defence signals at the mature disease stage, and silencing of the host's recognition and defence mechanisms that would otherwise protect it from further attacks. In contrast, calcium accumulation and a heightened defence response were observed in the moderately resistant cultivar in the later stage of disease development. Furthermore, in the susceptible interaction, the asymptomatic green region was unable to recover after disease disruption. Our targeted sampling technique also enabled detection of eight previously predicted proteinaceous effectors in addition to the known ToxA effector. Collectively, our results highlight the benefits of spatially resolved molecular analysis and nutrient mapping to provide high-resolution spatiotemporal snapshots of host-pathogen interactions, paving the way for disentangling complex disease interactions in plants.

Gene expression profile of the developing endosperm in durum wheat provides insight into starch biosynthesis

BACKGROUND

Durum wheat (Triticum turgidum subsp. durum) is widely grown for pasta production, and more recently, is gaining additional interest due to its resilience to warm, dry climates and its use as an experimental model for wheat research. Like in bread wheat, the starch and protein accumulated in the endosperm during grain development are the primary contributors to the calorific value of durum grains.

RESULTS

To enable further research into endosperm development and storage reserve synthesis, we generated a high-quality transcriptomics dataset from developing endosperms of variety Kronos, to complement the extensive mutant resources available for this variety. Endosperms were dissected from grains harvested at eight timepoints during grain development (6 to 30 days post anthesis (dpa)), then RNA sequencing was used to profile the transcriptome at each stage. The largest changes in gene expression profile were observed between the earlier timepoints, prior to 15 dpa. We detected a total of 29,925 genes that were significantly differentially expressed between at least two timepoints, and clustering analysis revealed nine distinct expression patterns. We demonstrate the potential of our dataset to provide new insights into key processes that occur during endosperm development, using starch metabolism as an example.

CONCLUSION

We provide a valuable resource for studying endosperm development in this increasingly important crop species.

MOBILE pipeline enables identification of context-specific networks and regulatory mechanisms

Robust identification of context-specific network features that control cellular phenotypes remains a challenge. We here introduce MOBILE (Multi-Omics Binary Integration via Lasso Ensembles) to nominate molecular features associated with cellular phenotypes and pathways. First, we use MOBILE to nominate mechanisms of interferon-γ (IFNγ) regulated PD-L1 expression. Our analyses suggest that IFNγ-controlled PD-L1 expression involves BST2, CLIC2, FAM83D, ACSL5, and HIST2H2AA3 genes, which were supported by prior literature. We also compare networks activated by related family members transforming growth factor-beta 1 (TGFβ1) and bone morphogenetic protein 2 (BMP2) and find that differences in ligand-induced changes in cell size and clustering properties are related to differences in laminin/collagen pathway activity. Finally, we demonstrate the broad applicability and adaptability of MOBILE by analyzing publicly available molecular datasets to investigate breast cancer subtype specific networks. Given the ever-growing availability of multi-omics datasets, we envision that MOBILE will be broadly useful for identification of context-specific molecular features and pathways.

PHOSPHORUS-STARVATION TOLERANCE 1 (OsPSTOL1) is prevalent in upland rice and enhances root growth and hastens low phosphate signaling in wheat

PHOSPHORUS-STARVATION TOLERANCE 1 (OsPSTOL1) is a variably present gene that benefits crown root growth and phosphorus (P) sufficiency in rice (Oryza sativa). To explore the ecophysiological importance of this gene, we performed a biogeographic survey of landraces and cultivars, confirming that functional OsPSTOL1 alleles prevail in low nutrient and drought-prone rainfed ecosystems, whereas loss-of-function and absence haplotypes predominate in control-irrigated paddy varieties of east Asia. An evolutionary history analysis of OsPSTOL1 and related genes in cereal, determined it and other genes are kinase-only domain derivatives of membrane-associated receptor like kinases. Finally, to evaluate the potential value of this kinase of unknown function in another Gramineae, wheat (Triticum aestivum) lines overexpressing OsPSTOL1 were evaluated under field and controlled low P conditions. OsPSTOL1 enhances growth, crown root number, and overall root plasticity under low P in wheat. Survey of root and shoot crown transcriptomes at two developmental stages identifies transcription factors that are differentially regulated in OsPSTOL1 wheat that are similarly controlled by the gene in rice. In wheat, OsPSTOL1 alters the timing and amplitude of regulators of root development in dry soils and hastens induction of the core P-starvation response. OsPSTOL1 and related genes may aid more sustainable cultivation of cereal crops.

Metamorphic gene regulation programs in Xenopus tropicalis tadpole brain

Amphibian metamorphosis is controlled by thyroid hormone (TH), which binds TH receptors (TRs) to regulate gene expression programs that underlie morphogenesis. Gene expression screens using tissues from premetamorphic tadpoles treated with TH identified some TH target genes, but few studies have analyzed genome-wide changes in gene regulation during spontaneous metamorphosis. We analyzed RNA sequencing data at four developmental stages from the beginning to the end of spontaneous metamorphosis, conducted on the neuroendocrine centers of Xenopus tropicalis tadpole brain. We also conducted chromatin immunoprecipitation sequencing (ChIP-seq) for TRs, and we compared gene expression changes during metamorphosis with those induced by exogenous TH. The mRNA levels of 26% of protein coding genes changed during metamorphosis; about half were upregulated and half downregulated. Twenty four percent of genes whose mRNA levels changed during metamorphosis had TR ChIP-seq peaks. Genes involved with neural cell differentiation, cell physiology, synaptogenesis and cell-cell signaling were upregulated, while genes involved with cell cycle, protein synthesis, and neural stem/progenitor cell homeostasis were downregulated. There is a shift from building neural structures early in the metamorphic process, to the differentiation and maturation of neural cells and neural signaling pathways characteristic of the adult frog brain. Only half of the genes modulated by treatment of premetamorphic tadpoles with TH for 16 h changed expression during metamorphosis; these represented 33% of the genes whose mRNA levels changed during metamorphosis. Taken together, our results provide a foundation for understanding the molecular basis for metamorphosis of tadpole brain, and they highlight potential caveats for interpreting gene regulation changes in premetamorphic tadpoles induced by exogenous TH.

Effective root responses to salinity stress include maintained cell expansion and carbon allocation

Acclimation of root growth is vital for plants to survive salt stress. Halophytes are great examples of plants that thrive even under severe salinity, but their salt tolerance mechanisms, especially those mediated by root responses, are still largely unknown. We compared root growth responses of the halophyte Schrenkiella parvula with its glycophytic relative species Arabidopsis thaliana under salt stress and performed transcriptomic analysis of S. parvula roots to identify possible gene regulatory networks underlying their physiological responses. Schrenkiella parvula roots do not avoid salt and experience less growth inhibition under salt stress. Salt-induced abscisic acid levels were higher in S. parvula roots compared with Arabidopsis. Root transcriptomic analysis of S. parvula revealed the induction of sugar transporters and genes regulating cell expansion and suberization under salt stress. ¹⁴ C-labeled carbon partitioning analyses showed that S. parvula continued allocating carbon to roots from shoots under salt stress while carbon barely allocated to Arabidopsis roots. Further physiological investigation revealed that S. parvula roots maintained root cell expansion and enhanced suberization under severe salt stress. In summary, roots of S. parvula deploy multiple physiological and developmental adjustments under salt stress to maintain growth, providing new avenues to improve salt tolerance of plants using root-specific strategies.

A transcription factor of the NAC family regulates nitrate-induced legume nodule senescence

Legumes establish symbioses with rhizobia by forming nitrogen-fixing nodules. Nitrate is a major environmental factor that affects symbiotic functioning. However, the molecular mechanism of nitrate-induced nodule senescence is poorly understood. Comparative transcriptomic analysis reveals an NAC-type transcription factor in Lotus japonicus, LjNAC094, that acts as a positive regulator in nitrate-induced nodule senescence. Stable overexpression and mutant lines of NAC094 were constructed and used for phenotypic characterization. DNA-affinity purification sequencing was performed to identify NAC094 targeting genes and results were confirmed by electrophoretic mobility shift and transactivation assays. Overexpression of NAC094 induces premature nodule senescence. Knocking out NAC094 partially relieves nitrate-induced degradation of leghemoglobins and abolishes nodule expression of senescence-associated genes (SAGs) that contain a conserved binding motif for NAC094. Nitrate-triggered metabolic changes in wild-type nodules are largely affected in nac094 mutant nodules. Induction of NAC094 and its targeting SAGs was almost blocked in the nitrate-insensitive nlp1, nlp4, and nlp1 nlp4 mutants. We conclude that NAC094 functions downstream of NLP1 and NLP4 by regulating nitrate-induced expression of SAGs. Our study fills in a key gap between nitrate and the execution of nodule senescence, and provides a potential strategy to improve nitrogen fixation and stress tolerance of legumes.

Exploring the shared gene signatures of smoking-related osteoporosis and chronic obstructive pulmonary disease using machine learning algorithms

Objectives: Cigarette smoking has been recognized as a predisposing factor for both osteoporosis (OP) and chronic obstructive pulmonary disease (COPD). This study aimed to investigate the shared gene signatures affected by cigarette smoking in OP and COPD through gene expression profiling. Materials and methods: Microarray datasets (GSE11784, GSE13850, GSE10006, and GSE103174) were obtained from Gene Expression Omnibus (GEO) and analyzed for differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA). Least absolute shrinkage and selection operator (LASSO) regression method and a random forest (RF) machine learning algorithm were used to identify candidate biomarkers. The diagnostic value of the method was assessed using logistic regression and receiver operating characteristic (ROC) curve analysis. Finally, immune cell infiltration was analyzed to identify dysregulated immune cells in cigarette smoking-induced COPD. Results: In the smoking-related OP and COPD datasets, 2858 and 280 DEGs were identified, respectively. WGCNA revealed 982 genes strongly correlated with smoking-related OP, of which 32 overlapped with the hub genes of COPD. Gene Ontology (GO) enrichment analysis showed that the overlapping genes were enriched in the immune system category. Using LASSO regression and RF machine learning, six candidate genes were identified, and a logistic regression model was constructed, which had high diagnostic values for both the training set and external validation datasets. The area under the curves (AUCs) were 0.83 and 0.99, respectively. Immune cell infiltration analysis revealed dysregulation in several immune cells, and six immune-associated genes were identified for smoking-related OP and COPD, namely, mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1), tissue-type plasminogen activator (PLAT), sodium channel 1 subunit alpha (SCNN1A), sine oculis homeobox 3 (SIX3), sperm-associated antigen 9 (SPAG9), and vacuolar protein sorting 35 (VPS35). Conclusion: The findings suggest that immune cell infiltration profiles play a significant role in the shared pathogenesis of smoking-related OP and COPD. The results could provide valuable insights for developing novel therapeutic strategies for managing these disorders, as well as shedding light on their pathogenesis.

Temporal Exposure to Bt Insecticide Causes Oxidative Stress in Larval Midgut Tissue

Bacillus thuringiensis (Bt) three-domain Cry toxins are highly successful biological pesticides; however, the mechanism through which they cause death to targeted larval midgut cells is not fully understood. Herein, we challenged transgenic Bt-susceptible Drosophila melanogaster larvae with moderate doses of activated Cry1Ac toxin and assessed the midgut tissues after one, three, and five hours using transmission electron microscopy and transcriptome sequencing. Larvae treated with Cry1Ac showed dramatic changes to their midgut morphology, including shortened microvilli, enlarged vacuoles, thickened peritrophic membranes, and swelling of the basal labyrinth, suggesting water influx. Transcriptome analysis showed that innate immune responses were repressed, genes involved with cell death pathways were largely unchanged, and mitochondria-related genes were strongly upregulated following toxin exposure. Defective mitochondria produced after toxin exposure were likely to contribute to significant levels of oxidative stress, which represent a common physiological response to a range of toxic chemicals. Significant reductions in both mitochondrial aconitase activity and ATP levels in the midgut tissue supported a rapid increase in reactive oxygen species (ROS) following exposure to Cry1Ac. Overall, these findings support the role of water influx, midgut cell swelling, and ROS activity in response to moderate concentrations of Cry1Ac.

Core cellular and tissue-specific mechanisms enable desiccation tolerance in Craterostigma

Resurrection plants can survive prolonged life without water (anhydrobiosis) in regions with seasonal drying. This desiccation tolerance requires the coordination of numerous cellular processes across space and time, and individual plant tissues face unique constraints related to their function. Here, we analyzed the complex, octoploid genome of the model resurrection plant Craterostigma (C. plantagineum), and surveyed spatial and temporal expression dynamics to identify genetic elements underlying desiccation tolerance. Homeologous genes within the Craterostigma genome have divergent expression profiles, suggesting the subgenomes contribute differently to desiccation tolerance traits. The Craterostigma genome contains almost 200 tandemly duplicated early light-induced proteins, a hallmark trait of desiccation tolerance, with massive upregulation under water deficit. We identified a core network of desiccation-responsive genes across all tissues, but observed almost entirely unique expression dynamics in each tissue during recovery. Roots and leaves have differential responses related to light and photoprotection, autophagy and nutrient transport, reflecting their divergent functions. Our findings highlight a universal set of likely ancestral desiccation tolerance mechanisms to protect cellular macromolecules under anhydrobiosis, with secondary adaptations related to tissue function.

A-MYB and BRDT-dependent RNA Polymerase II pause release orchestrates transcriptional regulation in mammalian meiosis

During meiotic prophase I, spermatocytes must balance transcriptional activation with homologous recombination and chromosome synapsis, biological processes requiring extensive changes to chromatin state. We explored the interplay between chromatin accessibility and transcription through prophase I of mammalian meiosis by measuring genome-wide patterns of chromatin accessibility, nascent transcription, and processed mRNA. We find that Pol II is loaded on chromatin and maintained in a paused state early during prophase I. In later stages, paused Pol II is released in a coordinated transcriptional burst mediated by the transcription factors A-MYB and BRDT, resulting in ~3-fold increase in transcription. Transcriptional activity is temporally and spatially segregated from key steps of meiotic recombination: double strand breaks show evidence of chromatin accessibility earlier during prophase I and at distinct loci from those undergoing transcriptional activation, despite shared chromatin marks. Our findings reveal mechanisms underlying chromatin specialization in either transcription or recombination in meiotic cells.

Gene Coexpression Analysis Identifies Genes Associated with Chlorophyll Content and Relative Water Content in Pearl Millet

Pearl millet is a significant crop that is tolerant to abiotic stresses and is a staple food of arid regions. However, its underlying mechanisms of stress tolerance are not fully understood. Plant survival is regulated by the ability to perceive a stress signal and induce appropriate physiological changes. Here, we screened for genes regulating physiological changes such as chlorophyll content (CC) and relative water content (RWC) in response to abiotic stress by using "weighted gene coexpression network analysis" (WGCNA) and clustering changes in physiological traits, i.e., CC and RWC associated with gene expression. Genes' correlations with traits were defined in the form of modules, and different color names were used to denote a particular module. Modules are groups of genes with similar patterns of expression, which also tend to be functionally related and co-regulated. In WGCNA, the dark green module (7082 genes) showed a significant positive correlation with CC, and the black (1393 genes) module was negatively correlated with CC and RWC. Analysis of the module positively correlated with CC highlighted ribosome synthesis and plant hormone signaling as the most significant pathways. Potassium transporter 8 and monothiol glutaredoxin were reported as the topmost hub genes in the dark green module. In Clust analysis, 2987 genes were found to display a correlation with increasing CC and RWC. Furthermore, the pathway analysis of these clusters identified the ribosome and thermogenesis as positive regulators of RWC and CC, respectively. Our study provides novel insights into the molecular mechanisms regulating CC and RWC in pearl millet.

Transcription Factor Dynamics in Cross-Regulation of Plant Hormone Signaling Pathways

Cross-regulation between hormone signaling pathways is indispensable for plant growth and development. However, the molecular mechanisms by which multiple hormones interact and co-ordinate activity need to be understood. Here, we generated a cross-regulation network explaining how hormone signals are integrated from multiple pathways in etiolated Arabidopsis (Arabidopsis thaliana) seedlings. To do so we comprehensively characterized transcription factor activity during plant hormone responses and reconstructed dynamic transcriptional regulatory models for six hormones; abscisic acid, brassinosteroid, ethylene, jasmonic acid, salicylic acid and strigolactone/karrikin. These models incorporated target data for hundreds of transcription factors and thousands of protein-protein interactions. Each hormone recruited different combinations of transcription factors, a subset of which were shared between hormones. Hub target genes existed within hormone transcriptional networks, exhibiting transcription factor activity themselves. In addition, a group of MITOGEN-ACTIVATED PROTEIN KINASES (MPKs) were identified as potential key points of cross-regulation between multiple hormones. Accordingly, the loss of function of one of these (MPK6) disrupted the global proteome, phosphoproteome and transcriptome during hormone responses. Lastly, we determined that all hormones drive substantial alternative splicing that has distinct effects on the transcriptome compared with differential gene expression, acting in early hormone responses. These results provide a comprehensive understanding of the common features of plant transcriptional regulatory pathways and how cross-regulation between hormones acts upon gene expression.

Time-Course Transcriptome Profiling Reveals Differential Resistance Responses of Tomato to a Phytotoxic Effector of the Pathogenic Oomycete Phytophthora cactorum

Blight caused by Phytophthora pathogens has a devastating impact on crop production. Phytophthora species secrete an array of effectors, such as Phytophthora cactorum-Fragaria (PcF)/small cysteine-rich (SCR) phytotoxic proteins, to facilitate their infections. Understanding host responses to such proteins is essential to developing next-generation crop resistance. Our previous work identified a small, 8.1 kDa protein, SCR96, as an important virulence factor in Phytophthora cactorum. Host responses to SCR96 remain obscure. Here, we analyzed the effect of SCR96 on the resistance of tomato treated with this recombinant protein purified from yeast cells. A temporal transcriptome analysis of tomato leaves infiltrated with 500 nM SCR96 for 0, 3, 6, and 12 h was performed using RNA-Seq. In total, 36,779 genes, including 2704 novel ones, were detected, of which 32,640 (88.7%) were annotated. As a whole, 5929 non-redundant genes were found to be significantly co-upregulated in SCR96-treated leaves (3, 6, 12 h) compared to the control (0 h). The combination of annotation, enrichment, and clustering analyses showed significant changes in expression beginning at 3 h after treatment in genes associated with defense and metabolism pathways, as well as temporal transcriptional accumulation patterns. Noticeably, the expression levels of resistance-related genes encoding receptor-like kinases/proteins, resistance proteins, mitogen-activated protein kinases (MAPKs), transcription factors, pathogenesis-related proteins, and transport proteins were significantly affected by SCR96. Quantitative reverse transcription PCR (qRT-PCR) validated the transcript changes in the 12 selected genes. Our analysis provides novel information that can help delineate the molecular mechanism and components of plant responses to effectors, which will be useful for the development of resistant crops.

Wound healing and regeneration in the reef building coral Acropora millepora

Branching scleractinian corals are niche-constructing organisms, providing continuously-growing, structural foundation for spectacularly biodiverse coral reef ecosystems. A large part of their success lies in the ability to quickly regenerate following mechanical damage. Even now, when the corals undergo great decline due to anthropogenic weather and storm extremes, it is surprising how little is known about molecular mechanisms governing regeneration in these iconic organisms. In this study, we used RNA-seq to identify genes involved in the regeneration of Acropora millepora, starting with the initial wound closure up to complete rebuilding of lost structures. Many of the differentially expressed genes we found in the wound healing steps are homologues of genes known to be involved in wound healing and regeneration of bilaterian and other cnidarian species, prominently including multiple components of FGF and Wnt signalling pathways. Comparison between genes involved in wound healing and continuous growth of the colony demonstrates both similarity and distinctiveness of the genetic programmes controlling these processes. A striking example is specific expression of c-Fos, a transcription factor with conserved role in early injury response, during the earliest stages of wound healing of A. millepora. By comparing results obtained in diverse experimental conditions including a closed-loop, recirculating aquarium and a flow-through system of marine station, we have demonstrated feasibility of using zooxanthellate scleractinian corals as experimental models in fundamental biology research, including studies of regeneration.

Global transcriptome and gene co-expression network analyses reveal regulatory and non-additive effects of drought and heat stress in grapevine

Despite frequent co-occurrence of drought and heat stress, the molecular mechanisms governing plant responses to these stresses in combination have not often been studied. This is particularly evident in non-model, perennial plants. We conducted large scale physiological and transcriptome analyses to identify genes and pathways associated with grapevine response to drought and/or heat stress during stress progression and recovery. We identified gene clusters with expression correlated to leaf temperature and water stress and five hub genes for the combined stress co-expression network. Several differentially expressed genes were common to the individual and combined stresses, but the majority were unique to the individual or combined stress treatments. These included heat-shock proteins, mitogen-activated kinases, sugar metabolizing enzymes, and transcription factors, while phenylpropanoid biosynthesis and histone modifying genes were unique to the combined stress treatment. Following physiological recovery, differentially expressed genes were found only in plants under heat stress, both alone and combined with drought. Taken collectively, our results suggest that the effect of the combined stress on physiology and gene expression is more severe than that of individual stresses, but not simply additive, and that epigenetic chromatin modifications may play an important role in grapevine responses to combined drought and heat stress.

Targeted and explorative profiling of kallikrein proteases and global proteome biology of pancreatic ductal adenocarcinoma, chronic pancreatitis, and normal pancreas highlights disease-specific proteome remodelling

Pancreatic ductal adenocarcinoma (PDAC) represents one of the most aggressive and lethal malignancies worldwide with an urgent need for new diagnostic and therapeutic strategies. One major risk factor for PDAC is the pre-indication of chronic pancreatitis (CP), which represents highly inflammatory pancreatic tissue. Kallikreins (KLKs) are secreted serine proteases that play an important role in various cancers as components of the tumor microenvironment. Previous studies of KLKs in solid tumors largely relied on either transcriptomics or immunodetection. We present one of the first targeted mass spectrometry profiling of kallikrein proteases in PDAC, CP, and normal pancreas. We show that KLK6 and KLK10 are significantly upregulated in PDAC (n=14) but not in CP (n=7) when compared to normal pancreas (n=16), highlighting their specific intertwining with malignancy. Additional explorative proteome profiling identified 5936 proteins in our pancreatic cohort and observed disease-specific proteome rearrangements in PDAC and CP. As such, PDAC features an enriched proteome motif for extracellular matrix (ECM) and cell adhesion while there is depletion of mitochondrial energy metabolism proteins, reminiscent of the Warburg effect. Although often regarded as a PDAC hallmark, the ECM fingerprint was also observed in CP, alongside with a prototypical inflammatory proteome motif as well as with an increased wound healing process and proteolytic activity, thereby possibly illustrating tissue autolysis. Proteogenomic analysis based on publicly accessible data sources identified 112 PDAC-specific and 32 CP-specific single amino acid variants, which among others affect KRAS and ANKHD1. Our study emphasizes the diagnostic potential of kallikreins and provides novel insights into proteomic characteristics of PDAC and CP.

Clustering and Annotation of T Cell Receptor Repertoires

Immunological protection against a wide variety of pathogens is largely mediated by the diverse and dynamic T cell receptor (TCR) repertoire, a crucial component of the adaptive immune system. An encounter with infectious agents stimulates specific T cells to initiate a direct immune response to combat intruders. Hence, the TCR repertoire may conceal crucial information regarding current and past infections and might assist in the development and monitoring of vaccines. To unlock its knowledge, we describe a computational workflow involving both supervised and unsupervised machine learning techniques to analyze and annotate full TCR repertoire data. The method is explained using data from a published yellow fever virus (YFV) vaccination study in healthy individuals. The TCR repertoire of one individual is studied before and 2 weeks after vaccination, using an efficient clustering method and identification of YFV-specific TCRs.

Targeting integrated stress response with ISRIB combined with imatinib treatment attenuates RAS/RAF/MAPK and STAT5 signaling and eradicates chronic myeloid leukemia cells

The integrated stress response (ISR) facilitates cellular adaptation to unfavorable conditions by reprogramming the cellular response. ISR activation was reported in neurological disorders and solid tumors; however, the function of ISR and its role as a possible therapeutic target in hematological malignancies still remain largely unexplored. Previously, we showed that the ISR is activated in chronic myeloid leukemia (CML) cells and correlates with blastic transformation and tyrosine kinase inhibitor (TKI) resistance. Moreover, the ISR was additionally activated in response to imatinib as a type of protective internal signaling. Here, we show that ISR inhibition combined with imatinib treatment sensitized and more effectively eradicated leukemic cells both in vitro and in vivo compared to treatment with single agents. The combined treatment specifically inhibited the STAT5 and RAS/RAF/MEK/ERK pathways, which are recognized as drivers of resistance. Mechanistically, this drug combination attenuated both interacting signaling networks, leading to BCR-ABL1- and ISR-dependent STAT5 activation. Consequently, leukemia engraftment in patient-derived xenograft mice bearing CD34+ TKI-resistant CML blasts carrying PTPN11 mutation responsible for hyperactivation of the RAS/RAF/MAPK and JAK/STAT5 pathways was decreased upon double treatment. This correlated with the downregulation of genes related to the RAS/RAF/MAPK, JAK/STAT5 and stress response pathways and was associated with lower expression of STAT5-target genes regulating proliferation, viability and the stress response. Collectively, these findings highlight the effect of imatinib plus ISRIB in the eradication of leukemic cells resistant to TKIs and suggest potential clinical benefits for leukemia patients with TKI resistance related to RAS/RAF/MAPK or STAT5 signaling. We propose that personalized treatment based on the genetic selection of patients carrying mutations that cause overactivation of the targeted pathways and therefore make their sensitivity to such treatment probable should be considered as a possible future direction in leukemia treatment.

Identification of a naturally-occurring canine model for early detection and intervention research in high grade urothelial carcinoma

BACKGROUND

Early detection and intervention research is expected to improve the outcomes for patients with high grade muscle invasive urothelial carcinoma (InvUC). With limited patients in suitable high-risk study cohorts, relevant animal model research is critical. Experimental animal models often fail to adequately represent human cancer. The purpose of this study was to determine the suitability of dogs with high breed-associated risk for naturally-occurring InvUC to serve as relevant models for early detection and intervention research. The feasibility of screening and early intervention, and similarities and differences between canine and human tumors, and early and later canine tumors were determined.

METHODS

STs (n=120) ≥ 6 years old with no outward evidence of urinary disease were screened at 6-month intervals for 3 years with physical exam, ultrasonography, and urinalysis with sediment exam. Cystoscopic biopsy was performed in dogs with positive screening tests. The pathological, clinical, and molecular characteristics of the "early" cancer detected by screening were determined. Transcriptomic signatures were compared between the early tumors and published findings in human InvUC, and to more advanced "later" canine tumors from STs who had the typical presentation of hematuria and urinary dysfunction. An early intervention trial of an oral cyclooxygenase inhibitor, deracoxib, was conducted in dogs with cancer detected through screening.

RESULTS

Biopsy-confirmed bladder cancer was detected in 32 (27%) of 120 STs including InvUC (n=29, three starting as dysplasia), grade 1 noninvasive cancer (n=2), and carcinoma in situ (n=1). Transcriptomic signatures including druggable targets such as EGFR and the PI3K-AKT-mTOR pathway, were very similar between canine and human InvUC, especially within luminal and basal molecular subtypes. Marked transcriptomic differences were noted between early and later canine tumors, particularly within luminal subtype tumors. The deracoxib remission rate (42% CR+PR) compared very favorably to that with single-agent cyclooxygenase inhibitors in more advanced canine InvUC (17-25%), supporting the value of early intervention.

CONCLUSIONS

The study defined a novel naturally-occurring animal model to complement experimental models for early detection and intervention research in InvUC. Research incorporating the canine model is expected to lead to improved outcomes for humans, as well as pet dogs, facing bladder cancer.

Salt Stress Tolerance in Casuarina glauca: Insights from the Branchlets Transcriptome

Climate change and the accelerated rate of population growth are imposing a progressive degradation of natural ecosystems worldwide. In this context, the use of pioneer trees represents a powerful approach to reverse the situation. Among others, N₂-fixing actinorhizal trees constitute important elements of plant communities and have been successfully used in land reclamation at a global scale. In this study, we have analyzed the transcriptome of the photosynthetic organs of Casuarina glauca (branchlets) to unravel the molecular mechanisms underlying salt stress tolerance. For that, C. glauca plants supplied either with chemical nitrogen (KNO₃⁺) or nodulated by Frankia (NOD⁺) were exposed to a gradient of salt concentrations (200, 400, and 600 mM NaCl) and RNA-Seq was performed. An average of ca. 25 million clean reads was obtained for each group of plants, corresponding to 86,202 unigenes. The patterns of differentially expressed genes (DEGs) clearly separate two groups: (i) control- and 200 mM NaCl-treated plants, and (ii) 400 and 600 mM NaCl-treated plants. Additionally, although the number of total transcripts was relatively high in both plant groups, the percentage of significant DEGs was very low, ranging from 6 (200 mM NaCl/NOD⁺) to 314 (600 mM NaCl/KNO₃⁺), mostly involving down-regulation. The vast majority of up-regulated genes was related to regulatory processes, reinforcing the hypothesis that some ecotypes of C. glauca have a strong stress-responsive system with an extensive set of constitutive defense mechanisms, complemented by a tight mechanism of transcriptional and post-transcriptional regulation. The results suggest that the robustness of the stress response system in C. glauca is regulated by a limited number of genes that tightly regulate detoxification and protein/enzyme stability, highlighting the complexity of the molecular interactions leading to salinity tolerance in this species.

Proteome alterations during clonal isolation of established human pancreatic cancer cell lines

Clonal isolation is an integral step of numerous workflows in genome editing and cell engineering. It comprises the isolation of a single progenitor cell from a defined cell line population with subsequent expansion to obtain a monoclonal cell population. This process is associated with transient loss of cell–cell contacts and absence of a multicellular microenvironment. Previous studies have revealed transcriptomic changes upon clonal isolation with cell line specific extent. Since transcriptome alterations are only partially reflected on the proteome level, we sought to investigate the impact of clonal isolation on the cellular proteome to a depth of > 6000 proteins in three established pancreatic cancer cell lines. We show that clonal isolation does have an impact on the cellular proteome, however, with cell line specific extent, affecting different biological processes, and also depending on the isolation method. We demonstrate a different impact of clonal isolation on mesenchymal- and epithelial-derived cell lines mainly affecting cell proliferation, metabolism, cell adhesion and cellular stress. The results bear relevance to the field of genomic editing and cell engineering and highlight the need to consider the impact of clonal isolation when interpreting data stemming from experiments that include this step.

SMGR: a joint statistical method for integrative analysis of single-cell multi-omics data

Unravelling the regulatory programs from single-cell multi-omics data has long been one of the major challenges in genomics, especially in the current emerging single-cell field. Currently there is a huge gap between fast-growing single-cell multi-omics data and effective methods for the integrative analysis of these inherent sparse and heterogeneous data. In this study, we have developed a novel method, Single-cell Multi-omics Gene co-Regulatory algorithm (SMGR), to detect coherent functional regulatory signals and target genes from the joint single-cell RNA-sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) data obtained from different samples. Given that scRNA-seq and scATAC-seq data can be captured by zero-inflated Negative Binomial distribution, we utilize a generalized linear regression model to identify the latent representation of consistently expressed genes and peaks, thus enables the identification of co-regulatory programs and the elucidation of regulating mechanisms. Results from both simulation and experimental data demonstrate that SMGR outperforms the existing methods with considerably improved accuracy. To illustrate the biological insights of SMGR, we apply SMGR to mixed-phenotype acute leukemia (MPAL) and identify the MPAL-specific regulatory program with significant peak-gene links, which greatly enhance our understanding of the regulatory mechanisms and potential targets of this complex tumor.

Acute, Chronic, and Treated Aortic Diseases Present Distinguishable Serum Proteome Fingerprints with Protein Profiles That Correlate with Disease Severity

Aortic diseases are a rare but potentially life-threatening condition. We present a serum proteomic study for a spectrum of aortic diseases including thoracic aortic aneurysms (n = 11), chronic dissections (n = 9), acute aortic dissections (n = 11), and surgically treated dissections (n = 19) as well as healthy controls (n = 10) and patients of coronary heart disease (n = 10) to represent non-aortic cardiovascular disease. In total, we identified and quantified 425 proteins across all 70 samples. The different aortic diseases represented distinguishable proteome profiles. We identified protein clusters that positively or negatively correlate with disease severity, including increase of cytosolic tissue leakage proteins and decrease of components of the coagulation and complement system. Further, we identified a serum proteome fingerprint of acute aortic dissections, consisting, among others, of enriched inflammatory markers such as C-reactive protein and members of the S100 protein family. The study underlines the applicability of serum proteomics for the investigation of aortic diseases and highlights the possibility to establish disease-specific prognostic markers.

Temporally coordinated expression of nuclear genes encoding chloroplast proteins in wheat promotes Puccinia striiformis f. sp. tritici infection

Targeting host processes that allow pathogens to thrive can be invaluable in resistance breeding. Here, we generated a deep-sequencing transcriptome time course for Puccinia striiformis f. sp. tritici (Pst) infection on wheat and compared datasets from three wheat varieties with different levels of susceptibility to two tested pathogen isolates. We sought genes specifically altered in a susceptible host as candidates that might support colonisation. Host responses differed between Pst-varietal pairs most prominently early during infection. Notably, however, nuclear genes encoding chloroplast-localised proteins (NGCPs) exhibited temporal coordination of expression profiles that differed at later time points in relation to Pst susceptibility. Disrupting one such NGCP, encoding the chloroplast-localised RNA binding protein TaCSP41a, led to lower Pst susceptibility. These analyses thus highlight NGCPs as prime targets for Pst manipulation during infection and point to TaCSP41a disruption as a potential source of Pst resistance for breeding programmes.

A transcriptome time course of Puccinia striiformis f. sp. tritici (Pst) infection reveals nuclear genes encoding chloroplast-localized proteins are manipulated during infection and highlights TaCSP41a disruption as a target for resistance breeding.

Differential Expression of Genes between a Tolerant and a Susceptible Maize Line in Response to a Sugarcane Mosaic Virus Infection

To uncover novel genes associated with the Sugarcane mosaic virus (SCMV) response, we used RNA-Seq data to analyze differentially expressed genes (DEGs) and transcript expression pattern clusters between a tolerant/resistant (CI-RL1) and a susceptible (B73) line, in addition to the F1 progeny (CI-RL1xB73). A Gene Ontology (GO) enrichment of DEGs led us to propose three genes possibly associated with the CI-RL1 response: a heat shock 90-2 protein and two ABC transporters. Through a clustering analysis of the transcript expression patterns (CTEPs), we identified two genes putatively involved in viral systemic spread: the maize homologs to the PIEZO channel (ZmPiezo) and to the Potyvirus VPg Interacting Protein 1 (ZmPVIP1). We also observed the complex behavior of the maize eukaryotic factors ZmeIF4E and Zm-elfa (involved in translation), homologs to eIF4E and eEF1α in A. thaliana. Together, the DEG and CTEPs results lead us to suggest that the tolerant/resistant CI-RL1 response to the SCMV encompasses the action of diverse genes and, for the first time, that maize translation factors are associated with viral interaction.

Transforming growth factor β1 impairs the transcriptomic response to contraction in myotubes from women with polycystic ovary syndrome

Abstract

Polycystic ovary syndrome (PCOS) is characterised by a hormonal imbalance affecting the reproductive and metabolic health of reproductive‐aged women. Exercise is recommended as a first‐line therapy for women with PCOS to improve their overall health; however, women with PCOS are resistant to the metabolic benefits of exercise training. Here, we aimed to gain insight into the mechanisms responsible for such resistance to exercise in PCOS. We employed an in vitro approach with electrical pulse stimulation (EPS) of cultured skeletal muscle cells to explore whether myotubes from women with PCOS have an altered gene expression signature in response to contraction. Following EPS, 4719 genes were differentially expressed (false discovery rate <0.05) in myotubes from women with PCOS compared to 173 in healthy women. Both groups included genes involved in skeletal muscle contraction. We also determined the effect of two transforming growth factor β (TGFβ) ligands that are elevated in plasma of women with PCOS, TGFβ1 and anti‐Müllerian hormone (AMH), alone and on the EPS‐induced response. While AMH (30 ng/ml) had no effect, TGFβ1 (5 ng/ml) induced the expression of extracellular matrix genes and impaired the exercise‐like transcriptional signature in myotubes from women with and without PCOS in response to EPS by interfering with key processes related to muscle contraction, calcium transport and actin filament. Our findings suggest that while the fundamental gene expression responses of skeletal muscle to contraction is intact in PCOS, circulating factors like TGFβ1 may be responsible for the impaired adaptation to exercise in women with PCOS.

Key points

Gene expression responses to in vitro contraction (electrical pulse stimulation, EPS) are altered in myotubes from women with polycystic ovary syndrome (PCOS) compared to healthy controls, with an increased expression of genes related to pro‐inflammatory pathways.

Transforming growth factor β1 (TGFβ1) upregulates genes related to extracellular matrix remodelling and reduces the expression of contractile genes in myotubes, regardless of the donor's health status.

TGFβ1 alters the gene expression response to EPS, providing a possible mechanism for the impaired exercise adaptations in women with PCOS.

Serum proteome alterations during conventional and extracorporeal resuscitation in pigs

Background

Only a small number of patients survive an out-of-hospital cardiac arrest (CA) and can be discharged from hospital alive with a large percentage of these patients retaining neurological impairments. In recent years, extracorporeal cardiopulmonary resuscitation (ECPR) has emerged as a beneficial strategy to optimize cardiac arrest treatment. However, ECPR is still associated with various complications. To reduce these problems, a profound understanding of the underlying mechanisms is required. This study aims to investigate the effects of CA, conventional cardiopulmonary resuscitation (CPR) and ECPR using a whole-body reperfusion protocol (controlled and automated reperfusion of the whole body—CARL) on the serum proteome profiles in a pig model of refractory CA.

Methods

N = 7 pigs underwent 5 min of untreated CA followed by 30 min CPR and 120 min perfusion with CARL. Blood samples for proteomic analysis were drawn at baseline, after CPR and at the end of the CARL period. Following albumin-depletion, proteomic analysis was performed using liquid chromatography–tandem mass spectrometry.

Results

N = 21 serum samples were measured resulting in the identification and quantification of 308–360 proteins per sample and 388 unique proteins in total. The three serum proteome profiles at the investigated time points clustered individually and segregated almost completely when considering a 90% confidence interval. Differential expression analysis showed significant abundance changes in 27 proteins between baseline and after CPR and in 9 proteins after CARL compared to CPR. Significant findings were further validated through a co-abundance cluster analysis corroborating the observed abundance changes.

Conclusions

The presented data highlight the impact of systemic ischemia and reperfusion on the entire serum proteome during resuscitation with a special focus on changes regarding haemolysis, coagulation, inflammation, and cell-death processes. Generally, the observed changes contribute to post-ischemic complications. Better understanding of the underlying mechanisms during CA and resuscitation may help to limit these complications and improve therapeutic options.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03441-4.

Exploration of induced sputum BIRC3 levels and clinical implications in asthma

BACKGROUND

Baculoviral IAP repeat-containing 3 (BIRC3) which encodes a member of the IAP family of proteins upregulated in the asthma expression profile dataset. However, there was few research on studying the clinical implication of BIRC3 in asthma.

OBJECTIVE

To validate BIRC3 expression and its clinical implications in induced sputum of asthma.

METHODS

Based on the GSE76262 (118 asthma cases and 21 healthy controls) dataset, differentially expressed genes were screened using R software. Subsequently, BIRC3 mRNA and protein were clinically verified in induced sputum samples through quantitative real-time polymerase chain reaction (qRT-PCR) and enzyme-linked immunosorbent assay (ELISA). Besides, the correlations between BIRC3 expression and asthmatic eosinophilic/allergic inflammation indicators (FeNO, IgE, and EOS%), pulmonary function (FEV₁, FEV₁% pred, FVC% pred, and FEV₁/FVC), and inflammatory cytokines (IL-4, IL-5, IL-13, IL-25, IL-10, IL-33, and TSLP) were analyzed. Finally, BIRC3 mRNA was detected in human primary bronchial epithelial cells stimulated by cytokines (IL-4 or IL-13).

RESULTS

BIRC3 was screened as a candidate gene in the GSE76262, which was highly expressed in asthma. Highly expressed BIRC3 was positively correlated with eosinophilic and allergic indicators, including FeNO, blood eosinophil, and serum IgE. Moreover, BIRC3 protein was positively associated with inflammation cytokines, like IL-4, IL-5, IL-13, IL-25, IL-10, IL-33, and TSLP, while negatively correlated with FEV1, FEV1%pred, FVC% pred, and FEV1/FVC. Furthermore, the expression of BIRC3 could be induced in primary bronchial epithelial cells treated by cytokines IL-4 or IL-13.

CONCLUSIONS

BIRC3 significantly increased in induced sputum of asthma and positively correlated with airway eosinophilic and peripheral blood allergic inflammation, type 2 cytokines, and airway obstruction. Increased BIRC3 might be involved in the pathogenesis of asthma by affecting the eosinophilic and allergic inflammation.

Diminished Auxin Signaling Triggers Cellular Reprogramming by Inducing a Regeneration Factor in the Liverwort Marchantia polymorpha

Regeneration in land plants is accompanied by the establishment of new stem cells, which often involves reactivation of the cell division potential in differentiated cells. The phytohormone auxin plays pivotal roles in this process. In bryophytes, regeneration is enhanced by the removal of the apex and repressed by exogenously applied auxin, which has long been proposed as a form of apical dominance. However, the molecular basis behind these observations remains unexplored. Here, we demonstrate that in the liverwort Marchantia polymorpha, the level of endogenous auxin is transiently decreased in the cut surface of decapitated explants, and identify by transcriptome analysis a key transcription factor gene, LOW-AUXIN RESPONSIVE (MpLAXR), which is induced upon auxin reduction. Loss of MpLAXR function resulted in delayed cell cycle reactivation, and transient expression of MpLAXR was sufficient to overcome the inhibition of regeneration by exogenously applied auxin. Furthermore, ectopic expression of MpLAXR caused cell proliferation in normally quiescent tissues. Together, these data indicate that decapitation causes a reduction of auxin level at the cut surface, where, in response, MpLAXR is up-regulated to trigger cellular reprogramming. MpLAXR is an ortholog of Arabidopsis ENHANCER OF SHOOT REGENERATION 1/DORNRÖSCHEN, which has dual functions as a shoot regeneration factor and a regulator of axillary meristem initiation, the latter of which requires a low auxin level. Thus, our findings provide insights into stem cell regulation as well as apical dominance establishment in land plants.

Partially Resistant Avocado Rootstock Dusa® Shows Prolonged Upregulation of Nucleotide Binding-Leucine Rich Repeat Genes in Response to Phytophthora cinnamomi Infection

Avocado is an important agricultural food crop in many countries worldwide. Phytophthora cinnamomi, a hemibiotrophic oomycete, remains one of the most devastating pathogens within the avocado industry, as it is near impossible to eradicate from areas where the pathogen is present. A key aspect to Phytophthora root rot disease management is the use of avocado rootstocks partially resistant to P. cinnamomi, which demonstrates an increased immune response following infection. In plant species, Nucleotide binding-Leucine rich repeat (NLR) proteins form an integral part of pathogen recognition and Effector triggered immune responses (ETI). To date, a comprehensive set of Persea americana NLR genes have yet to be identified, though their discovery is crucial to understanding the molecular mechanisms underlying P. americana-P. cinnamomi interactions. In this study, a total of 161 PaNLR genes were identified in the P. americana West-Indian pure accession genome. These putative resistance genes were characterized using bioinformatic approaches and grouped into 13 distinct PaNLR gene clusters, with phylogenetic analysis revealing high sequence similarity within these clusters. Additionally, PaNLR expression levels were analyzed in both a partially resistant (Dusa®) and a susceptible (R0.12) avocado rootstock infected with P. cinnamomi using an RNA-sequencing approach. The results showed that the partially resistant rootstock has increased expression levels of 84 PaNLRs observed up to 24 h post-inoculation, while the susceptible rootstock only showed increased PaNLR expression during the first 6 h post-inoculation. Results of this study may indicate that the partially resistant avocado rootstock has a stronger, more prolonged ETI response which enables it to suppress P. cinnamomi growth and combat disease caused by this pathogen. Furthermore, the identification of PaNLRs may be used to develop resistant rootstock selection tools, which can be employed in the avocado industry to accelerate rootstock screening programs.

Estimation of the proportion of true null hypotheses under sparse dependence: Adaptive FDR controlling in microarray data

The proportion of non-differentially expressed genes is an important quantity in microarray data analysis and an appropriate estimate of the same is used to construct adaptive multiple testing procedures. Most of the estimators for the proportion of true null hypotheses based on the thresholding, maximum likelihood and density estimation approaches assume independence among the gene expressions. Usually, sparse dependence structure is natural in modelling associations in microarray gene expression data and hence it is necessary to develop methods for accommodating the sparse dependence well within the framework of existing estimators. We propose a clustering based method to put genes in the same group that are not coexpressed using the estimated high dimensional correlation structure under sparse assumption as dissimilarity matrix. This novel method is applied to three existing estimators for the proportion of true null hypotheses. Extensive simulation study shows that the proposed method improves an existing estimator by making it less conservative and the corresponding adaptive Benjamini-Hochberg algorithm more powerful. The proposed method is applied to a microarray gene expression dataset of colorectal cancer patients and the results show gain in terms of number of differentially expressed genes. The R code is available at https://github.com/aniketstat/Proportiontion-of-true-null-under-sparse-dependence-2021.

Non-cell autonomous mechanisms control mitochondrial gene dysregulation in polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is a common endocrine disorder associated with insulin resistance and impaired energy metabolism in skeletal muscle, the aetiology of which is currently unclear. Here, we mapped the gene expression profile of skeletal muscle from women with PCOS and determined if cultured primary myotubes retain the gene expression signature of PCOS in vivo. Transcriptomic analysis of vastus lateralis biopsies collected from PCOS women showed lower expression of genes associated with mitochondrial function, while the expression of genes associated with the extracellular matrix was higher compared to controls. Altered skeletal muscle mRNA expression of mitochondrial-associated genes in PCOS was associated with lower protein expression of mitochondrial complex II-V, but not complex I, with no difference in mitochondrial DNA content. Transcriptomic analysis of primary myotube cultures established from biopsies did not display any differentially expressed genes between controls and PCOS. Comparison of gene expression profiles in skeletal muscle biopsies and primary myotube cultures showed lower expression of mitochondrial and energy metabolism-related genes in vitro, irrespective of the group. Together, our results show that the altered mitochondrial-associated gene expression in skeletal muscle in PCOS is not preserved in cultured myotubes, indicating that the in vivo extracellular milieu, rather than genetic or epigenetic factors, may drive this alteration. Dysregulation of mitochondrial-associated genes in skeletal muscle by extracellular factors may contribute to the impaired energy metabolism associated with PCOS.