Transcriptional Regulatory Pathways Research Articles

Transcription factor HNF4A activity is suppressed by microbiota and inhibits IFN beta production

Intestinal epithelial cells (IECs) separate the intestinal lumen from lamina propria resident immune cells, positioning IECs as mediators of microbial signals to the immune system. Understanding these microbe-IEC-immune interactions may represent potential therapeutic targets for immunological diseases. We previously discovered that microbiota regulate the activity of the nuclear receptor transcription factor Hepatocyte nuclear factor 4 alpha (HNF4A). Microbiota suppress HNF4A-activated genes in zebrafish and lead to genome-wide reduction of HNF4A DNA occupancy in mouse IECs. HNF4A activates lipid metabolism and epithelial differentiation gene programs, while also suppressing inflammatory and immune response gene programs in the intestine. We therefore hypothesized that microbial suppression of HNF4A contributes to activation of immunological pathways. We conducted an unbiased coimmunoprecipitation-mass-spectrometry analysis comparing HNF4A interacting partners in IECs from germ-free and conventionalized (germ-free mice colonized with SPF feces for 2 weeks) mice. HNF4A interacted with significantly more proteins in germ-free IECs, including several proteins involved in promotion of transcription and chromatin regulation, consistent with HNF4A being more active in germ-free mice. Our HNF4A and H3k27ac ChIP-seq data from IECs derived from germ-free and conventionalized mice identified that HNF4A binds ubiquitously to IEC enhancers, but also revealed a small number of enhancers that were activated by microbiota with no HNF4A binding. Motif enrichment analysis showed that these enhancers were enriched for binding sites for Interferon regulatory factors (IRFs). This raises the interesting possibility that there is a selective disadvantage to colocalizing IRF and HNF4A occupancy at a given enhancer, and suggests that these transcriptional regulatory pathways may otherwise communicate. To test this, we used a luciferase reporter assay in HEK293T cells to probe the impact of HNF4A on the IFNβ promoter, which is controlled by IRF3/7. HNF4A potently suppressed IFNβ promoter activity following viral stimulation. HNF4A DNA-binding activity is necessary for this suppression; however, there is no evidence HNF4A binds the IFNβ promoter, suggesting HNF4A likely induces a repressor of the IFNβ promoter. Together these results indicate that HNF4A is both a target of microbiota suppression and a potent suppressor of IFNβ production. Though IRF bound enhancers appear to be insulated against HNF4A binding in IECs, our results show that these two transcriptional pathways communicate, potentially representing a novel mechanism through which microbiota influence innate immunity. NIH P01-DK094779 This is the full abstract presented at the American Physiology Summit 2023 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Genomic and transcriptomic analyses support a silk gland origin of spider venom glands

BackgroundSpiders comprise a hyperdiverse lineage of predators with venom systems, yet the origin of functionally novel spider venom glands remains unclear. Previous studies have hypothesized that spider venom glands originated from salivary glands or evolved from silk-producing glands present in early chelicerates. However, there is insufficient molecular evidence to indicate similarity among them. Here, we provide comparative analyses of genome and transcriptome data from various lineages of spiders and other arthropods to advance our understanding of spider venom gland evolution.ResultsWe generated a chromosome-level genome assembly of a model spider species, the common house spider (Parasteatoda tepidariorum). Module preservation, GO semantic similarity, and differentially upregulated gene similarity analyses demonstrated a lower similarity in gene expressions between the venom glands and salivary glands compared to the silk glands, which questions the validity of the salivary gland origin hypothesis but unexpectedly prefers to support the ancestral silk gland origin hypothesis. The conserved core network in the venom and silk glands was mainly correlated with transcription regulation, protein modification, transport, and signal transduction pathways. At the genetic level, we found that many genes in the venom gland-specific transcription modules show positive selection and upregulated expressions, suggesting that genetic variation plays an important role in the evolution of venom glands.ConclusionsThis research implies the unique origin and evolutionary path of spider venom glands and provides a basis for understanding the diverse molecular characteristics of venom systems.

The LibHLH22 and LibHLH63 from Lilium ‘Siberia’ Can Positively Regulate Volatile Terpenoid Biosynthesis

Basic helix-loop-helix (bHLH) transcription factors (TFs) play irreplaceable roles in plant growth and development, especially in plant secondary metabolism. However, the functions of most bHLH TFs in Lilium ‘Siberia’ are still unknown, especially their roles in regulating floral fragrance. In this study, two bHLH TFs in lily, i.e., LibHLH22 and LibHLH63, were identified and functionally characterized. A bioinformatics analysis revealed that LibHLH22 and LibHLH63 were unstable proteins. Subcellular localization demonstrated that LibHLH22 and LibHLH63 proteins were in the cell nucleus. Quantitative real-time PCR showed that the highest expression level of LibHLH22 was at the initial flowering stage and in the stigma, and the highest expression level of LibHLH63 was at the budding stage and in the filaments. The results of transient overexpression and virus-induced gene silencing (VIGS) of LibHLH22 and LibHLH63 in lily petals showed that these two transcription factors significantly promoted the expression of LiDXR and LiTPS2, and thus, markedly enhanced the release of floral fragrance. Our results indicated that LibHLH22 and LibHLH63 could effectively regulate the fragrance of Lilium ‘Siberia’, laying the foundation for fragrance breeding and improving the terpenoid transcriptional regulatory pathway.

Computational approaches for detecting disease-associated alternative splicing events.

Alternative splicing (AS) is a key transcriptional regulation pathway. Recent studies have shown that AS events are associated with the occurrence of complex diseases. Various computational approaches have been developed for the detection of disease-associated AS events. In this review, we first describe the metrics used for quantitative characterization of AS events. Second, we review and discuss the three types of methods for detecting disease-associated splicing events, which are differential splicing analysis, aberrant splicing detection and splicing-related network analysis. Third, to further exploit the genetic mechanism of disease-associated AS events, we describe the methods for detecting genetic variants that potentially regulate splicing. For each type of methods, we conducted experimental comparison to illustrate their performance. Finally, we discuss the limitations of these methods and point out potential ways to address them. We anticipate that this review provides a systematic understanding of computational approaches for the analysis of disease-associated splicing.

Degradation mechanism of AtrA mediated by ClpXP and its application in daptomycin production in Streptomyces roseosporus.

The efficiency of drug biosynthesis depends on different transcriptional regulatory pathways in Streptomyces, and the protein degradation system adds another layer of complexity to the regulatory processes. AtrA, a transcriptional regulator in the A-factor regulatory cascade, stimulates the production of daptomycin by binding to the dptE promoter in Streptomyces roseosporus. Using pull-down assays, bacterial two-hybrid system and knockout verification, we demonstrated that AtrA is a substrate for ClpP protease. Furthermore, we showed that ClpX is necessary for AtrA recognition and subsequent degradation. Bioinformatics analysis, truncating mutation, and overexpression proved that the AAA motifs of AtrA were essential for initial recognition in the degradation process. Finally, overexpression of mutated atrA (AAA-QQQ) in S. roseosporus increased the yield of daptomycin by 225% in shake flask and by 164% in the 15 L bioreactor. Thus, improving the stability of key regulators is an effective method to promote the ability of antibiotic synthesis.

Dynamic activity in cis-regulatory elements of leukocytes identifies transcription factor activation and stratifies COVID-19 severity in ICU patients

Transcription factor programs mediating the immune response to coronavirus disease 2019 (COVID-19) are not fully understood. Capturing active transcription initiation from cis-regulatory elements such as enhancers and promoters by capped small RNA sequencing (csRNA-seq), in contrast to capturing steady-state transcripts by conventional RNA-seq, allows unbiased identification of the underlying transcription factor activity and regulatory pathways. Here, we profile transcription initiation in critically ill COVID-19 patients, identifying transcription factor motifs that correlate with clinical lung injury and disease severity. Unbiased clustering reveals distinct subsets of cis-regulatory elements that delineate the cell type, pathway-specific, and combinatorial transcription factor activity. We find evidence of critical roles of regulatory networks, showing that STAT/BCL6 and E2F/MYB regulatory programs from myeloid cell populations are activated in patients with poor disease outcomes and associated with COVID-19 susceptibility genetic variants. More broadly, we demonstrate how capturing acute, disease-mediated changes in transcription initiation can provide insight into the underlying molecular mechanisms and stratify patient disease severity.

Comprehensive characteristics of pathological subtypes in testicular germ cell tumor: Gene expression, mutation and alternative splicing.

Testicular germ cell tumor (TGCT) is the most common tumor in young men, but molecular signatures, especially the alternative splicing (AS) between its subtypes have not yet been explored. To investigate the differences between TGCT subtypes, we comprehensively analyzed the data of gene expression, alternative splicing (AS), and somatic mutation in TGCT patients from the TCGA database. The gene ontology (GO) enrichment analyses were used to explore the function of differentially expressed genes and spliced genes respectively, and Spearman correlation analysis was performed to explore the correlation between differential genes and AS events. In addition, the possible patterns in which AS regulates gene expression were elaborated by the ensemble database transcript atlas. And, we identified important transcription factors that regulate gene expression and AS and functionally validated them in TGCT cell lines. We found significant differences between expression and AS in embryonal carcinoma and seminoma, while mixed cell tumors were in between. GO enrichment analyses revealed that both differentially expressed and spliced genes were enriched in transcriptional regulatory pathways, and obvious correlation between expression and AS events was determined. By analyzing the transcript map and the sites where splicing occurs, we have demonstrated that AS regulates gene expression in a variety of ways. We further identified two pivot AS-related molecules (SOX2 and HDAC9) involved in AS regulation, which were validated in embryonal carcinoma and seminoma cell lines. Differences in somatic mutations between subtypes are also of concern, with our results suggesting that mutations in some genes (B3GNT8, CAPN7, FAT4, GRK1, TACC2, and TRAM1L1) occur only in embryonal carcinoma, while mutations in KIT, KARS, and NRAS are observed only in seminoma. In conclusion, our analysis revealed the differences in gene expression, AS and somatic mutation among TGCT subtypes, providing a molecular basis for clinical diagnosis and precise therapy of TGCT patients.

PiRNAs regulate a Hedgehog germline-to-soma pro-aging signal

The reproductive system regulates somatic aging through competing anti- and pro-aging signals. Germline removal extends somatic lifespan through conserved pathways including insulin and mammalian target-of-rapamycin signaling, while germline hyperactivity shortens lifespan through unknown mechanisms. Here we show that mating-induced germline hyperactivity downregulates piRNAs, in turn desilencing their targets, including the Hedgehog-like ligand-encoding genes wrt-1 and wrt-10, ultimately causing somatic collapse and death. Germline-produced Hedgehog signals require PTR-6 and PTR-16 receptors for mating-induced shrinking and death. Our results reveal an unconventional role of the piRNA pathway in transcriptional regulation of Hedgehog signaling and a new role of Hedgehog signaling in the regulation of longevity and somatic maintenance: Hedgehog signaling is controlled by the tunable piRNA pathway to encode the previously unknown germline-to-soma pro-aging signal. Mating-induced piRNA downregulation in the germline and subsequent Hedgehog signaling to the soma enable the animal to tune somatic resource allocation to germline needs, optimizing reproductive timing and survival.

A signalling pathway for transcriptional regulation of sleep amount in mice.

In mice and humans, sleep quantity is governed by genetic factors and exhibits age-dependent variation1-3. However, the core molecular pathways and effector mechanisms that regulate sleep duration in mammals remain unclear. Here, we characterize a major signalling pathway for the transcriptional regulation of sleep in mice using adeno-associated virus-mediated somatic genetics analysis4. Chimeric knockout of LKB1 kinase-an activator of AMPK-related protein kinase SIK35-7-in adult mouse brain markedly reduces the amountand delta power-a measure of sleep depth-of non-rapid eye movement sleep (NREMS). Downstream of the LKB1-SIK3 pathway, gain or loss-of-function of the histone deacetylases HDAC4 and HDAC5 in adult brain neurons causes bidirectional changes of NREMS amount and delta power. Moreover, phosphorylation of HDAC4 and HDAC5 is associated with increasedsleep need, and HDAC4 specifically regulates NREMS amount in posterior hypothalamus. Genetic and transcriptomic studies reveal that HDAC4 cooperates with CREB in both transcriptional and sleep regulation. These findings introduce the concept of signalling pathways targeting transcription modulators to regulate daily sleep amount and demonstrate the power of somatic genetics in mouse sleep research.

Quantitative Sequencing Analysis of the Striatal Transcriptome in a Mouse Model of Alzheimer Disease.

The purpose of this study was to analyze the transcriptomic changes in the striatum of amyloid precursor protein/presenilin 1 (APP/PS1) transgenic mice and uncover its association with the methyl-CpG binding protein 2 (MeCP2) mediated-changes in striatal epigenetic signature during Alzheimer disease (AD) pathological progression. To observe transcriptomic alterations in the striatum before the onset of cognitive impairment in APP/PS1 mice, quantitative 3'mRNA sequencing was performed with RNA extracted from the striatum of 6-month-old and 12-month-old wildtype and APP/PS1 mice. In addition, chromatin immunoprecipitation sequencing was conducted with the DNA from wildtype and APP/PS1 mice of the same age as aforementioned. For transcriptomic analysis, comparison terms were constructed based on aging and transgene expression-normal-aging (12-month-old wildtype/6-month-old wildtype), early-AD (6-month-old APP/PS1/6-month-old wildtype), and late-AD (12-month-old APP/PS1/6-month-old wildtype). To compare the changes in biological pathways and networks, we analyzed gene lists from each comparison term via bioinformatics tools including DAVID (Database for Annotation, Visualization, and Integrated Discovery), STRING (Search Tool for the Retrieval of Interacting Genes/Proteins), and SynGO (Synaptic Gene Ontologies). Furthermore, to assume the effect MeCP2 in AD pathological conditions may have on the transcriptome regulation, analysis of the common genes from Quant-Seq and MeCP2-ChIP-Seq was performed. Enriched pathways including immune system and inflammatory response were confirmed in normal- aging and lateAD, respectively. In particular, enriched pathways of gene expression regulation, transcriptional regulation, and protein catabolic pathways were found to be significantly altered in early-AD. MeCP2-bound genes that were significantly altered in the transcriptome were suggested to be target genes that have a role in the striatum of the early-stage AD model. This study confirmed that the alteration of the striatal transcriptomic profile in APP/PS1 mice was involved with several biological pathways. Additionally, comparative analysis of the transcriptomic changes and the MeCP2 bound regions found that a group of differentially expressed genes may be regulated under the epigenetic control of MeCP2.

Global DNA and protein interactomes of FLT1P1 (Fms-related tyrosine kinase 1 pseudogene 1) revealed its molecular regulatory functions associated with preeclampsia.

Preeclampsia (PE) is one of the most serious pregnancy complications with unknown pathogenesis. Emerging evidence has demonstrated that Fms-related tyrosine kinase 1 (FLT1) is highly involved in PE development. As a pseudogene of FLT1, FLT1P1 increased in PE samples. However, its functions remain largely unknown. In this study, co-expression analysis was performed to identify the potential target genes of FTL1P1. Then chromatin isolation using RNA purification (ChIRP) method was employed to explore the interactomes of FLT1P1, including interacting with DNA fragments and proteins. We found that in PE samples, both FLT1P1 and FLT1 were highly expressed and closely correlated. ChIRP-protein data revealed that FLT1P1 interacts with translation- and transcription-related proteins, including 4 transcription factors (TFs). ChIRP-DNA analysis revealed that FLT1P1 preferentially interacted with DNA fragments downstream of transcription start sites (TSSs). Functional analysis of its interacting genes revealed that they were enriched in transcriptional regulation and apoptosis-related pathways. Twenty-six TFs, including CREB1 and SRF, were extracted from the potential FLT1P1-interacting gene sets and were potential targets of FLT1P1. CREB1 could bind to FLT1 promoter, and was negatively correlated with FLT1 at the expression level, making it a potential regulator of FLT1. Our study extensively investigated the interactome profiles of FLT1P1, especially the prompter region of TF gene CREB1, and revealed the potential molecular regulatory mechanisms of FLT1 expression in PE samples. Our results provide a novel view of PE pathogenesis, and suggest that FLT1P1 could serve as a potential therapeutic target in PE diagnosis and treatment.

Inherited Germline DNA Repair Gene Defects Are Prevalent Among Thai Patients with Myeloid Neoplasms

Certain germline mutations increase the risk of developing myeloid neoplasms (MNs) accompanied by second hit somatic events, such as therapy-related or oxidation-induced DNA damages. Aberrancy in DNA damage response (DDR) pathways leads to genomic instability giving rise to cancer. Previous data showed that the prevalence of pathogenic germline mutations in young adult Caucasian MN patients was 13-19%.1, 2 The most commonly affected pathways were involved in telomere biology and transcription factors, while the defects in DDRs were found in 15-25%. However, data in Asians are limited. In this study, we investigated the prevalence and landscape of germline variants in cancer predisposition genes in a cohort of Thai patients diagnosed with MNs. Targeted next generation sequencing of 23 germline predisposition genes was performed using the SureSelect custom kit and the Illumina HiSeq 2500 system. The total of 127 unselected adult (> 15 years old) patients from 3 medical centers comprising MDS (66%), MDS/MPN (14%) secondary AML (sAML, 20%) were enrolled. The median age was 66 years old (range 16-93) and 50% were male. Variants reported in COSMIC, the frequency of higher than 0.01% in the Genome Aggregation Database (gnomAD) or Thai exome database (consisted of 1084 Thai exome controls) and variant allelic frequency of less than 30% were excluded. New or exceeding rare nonsense variants or indels, as well as missense variants predicted to be deleterious by at least 4/6 in silico software, were included. The variants were classified according to American College of Medical Genetics (ACMG) classification. 14pathogenic and likely pathogenic germline variants were identified in 15 patients (15/127; 12%). One patient had 2 variants including heterozygous CEBPAP75fs and MSH2A604G. Of the 16 variants, there were 8 (50%) missense, 7 (44%) frameshift, and 1 (6%) nonsense, with a median VAF of 51%, ranging 36-83%. The variants affected the DDR pathway (8/16; 50%), transcriptional regulation pathway (4/16; 25%) and signal transduction pathway (2/16; 13%). Two variants were recurrent i.e.,FANCAAla1219GlyfsTer59 and BRCA2Glu626delinsAspPro and the others were identified in BRCA1/2, MSH2, MSH6, GATA2, CEBPA, DDX41, CBL and NF1. The median age of patients harboring germline variants was younger (58 vs. 66 years old) than those who did not. Half of the cases (6/12) were characterized by normal cytogenetics without enrichment of -7 or complex karyotype, comparable to patients with no germline variants. Of note, 28% (4/14) patients had previous history of gastrointestinal stromal tumor, Graves’ disease, Immune thrombocytopenia, and aplastic anemia (FANCA mutation) for years before MNs were diagnosed. None of the patients reported family history of hematological malignancies. In conclusion, the prevalence of germline predisposition to MNs in Southeast Asian populations is as high as in Caucasians and the identified variants are more commonly involved in the DDR pathway. The defects in the DDR pathway suggested an important role of genomic instability in leukemogenesis. This subset of patients should preferably receive non-cytotoxic treatments to avoid oxidative damages that possibly lead to further clonal evolution. ReferencesKeel SB, Scott A, Sanchez-Bonilla M, Ho PA, Gulsuner S, Pritchard CC, et al. Genetic features of myelodysplastic syndrome and aplastic anemia in pediatric and young adult patients. Haematologica 2016 Nov; 101(11): 1343-1350. Feurstein S, Churpek JE, Walsh T, Keel S, Hakkarainen M, Schroeder T, et al. Germline variants drive myelodysplastic syndrome in young adults. Leukemia 2021 2021/08/01; 35(8): 2439-2444.

Abstract 15391: Transcriptomic Expression Changes in Embryonic Mouse Hearts Following Single-Dose Maternal Binge Alcohol Consumption: Implications for Congenital Heart Defects

Introduction: Prenatal alcohol exposure is a risk factor for congenital heart defects, but it is unclear how alcohol exposure in utero leads to disruption of normal cardiovascular development. To gain insight into the transcriptomic pathways involved, we performed RNA-sequencing (RNA-seq) of embryonic mouse hearts 24h (E10.5) and 48h (E11.5) after a single dose of maternal alcohol. Hypothesis: Maternal binge alcohol consumption at E9.5 disrupts expression of specific embryonic genes important for heart development. Methods: Alcohol treatment (2.5g/Kg ethanol in saline) or control solution (0.9% saline) was administered to pregnant C57Bl/6 mice via single oral gavage at E9.5. Embryonic hearts were collected 24h and 48h later for RNA isolation and analysis by RNA-Seq and quantitative PCR (qPCR) methods. Reactome® pathway analysis was performed to determine which were most affected. Results: Of the nearly 20,000 genes expressed in the embryonic heart, we found remarkably few (less than 0.1%) showing significant alterations in expression. Most of the gene expression changes induced by alcohol occurred within 24h and were largely genes involved in transcriptional regulation or signaling pathways. In particular, Wnt7a showed significant ( p=0.013 ) and substantial change (nearly 5-fold decrease relative to control expression, n=11 per group) in response to maternal binge alcohol 24h post-gavage. Similar results were observed for Wnt7a at 48h post-binge (p=0.043, n=8/group). Further analyses showed that Wnt -related genes, Dkk2, Rtl1,Wnt3, Wnt6, and Wnt10a showed significantly decreased expression and Wnt11 showed significantly increased expression in the embryonic heart 24h post-binge alcohol. Conclusions: These results suggest that Wnt7a and related signaling pathway genes are disrupted by alcohol action in the embryonic heart, thereby opening new avenues for investigation into molecular underpinnings of alcohol-induced congenital heart defects.

A predictive ensemble classifier for the gene expression diagnosis of ASD at ages 1 to 4 years

Autism Spectrum Disorder (ASD) diagnosis remains behavior-based and the median age of diagnosis is ~52 months, nearly 5 years after its first-trimester origin. Accurate and clinically-translatable early-age diagnostics do not exist due to ASD genetic and clinical heterogeneity. Here we collected clinical, diagnostic, and leukocyte RNA data from 240 ASD and typically developing (TD) toddlers (175 toddlers for training and 65 for test). To identify gene expression ASD diagnostic classifiers, we developed 42,840 models composed of 3570 gene expression feature selection sets and 12 classification methods. We found that 742 models had AUC-ROC ≥ 0.8 on both Training and Test sets. Weighted Bayesian model averaging of these 742 models yielded an ensemble classifier model with accurate performance in Training and Test gene expression datasets with ASD diagnostic classification AUC-ROC scores of 85–89% and AUC-PR scores of 84–92%. ASD toddlers with ensemble scores above and below the overall ASD ensemble mean of 0.723 (on a scale of 0 to 1) had similar diagnostic and psychometric scores, but those below this ASD ensemble mean had more prenatal risk events than TD toddlers. Ensemble model feature genes were involved in cell cycle, inflammation/immune response, transcriptional gene regulation, cytokine response, and PI3K-AKT, RAS and Wnt signaling pathways. We additionally collected targeted DNA sequencing smMIPs data on a subset of ASD risk genes from 217 of the 240 ASD and TD toddlers. This DNA sequencing found about the same percentage of SFARI Level 1 and 2 ASD risk gene mutations in TD (12 of 105) as in ASD (13 of 112) toddlers, and classification based only on the presence of mutation in these risk genes performed at a chance level of 49%. By contrast, the leukocyte ensemble gene expression classifier correctly diagnostically classified 88% of TD and ASD toddlers with ASD risk gene mutations. Our ensemble ASD gene expression classifier is diagnostically predictive and replicable across different toddler ages, races, and ethnicities; out-performs a risk gene mutation classifier; and has potential for clinical translation.

Genome-Wide Identification and Expression Analysis of HSF Transcription Factors in Alfalfa (Medicago sativa) under Abiotic Stress.

Alfalfa (Medicago sativa) is one of the most important legume forage species in the world. It is often affected by several abiotic stressors that result in reduced yields and poor growth. Therefore, it is crucial to study the resistance of M. sativa to abiotic stresses. Heat shock transcription factors (HSF) are key players in a number of transcriptional regulatory pathways. These pathways play an essential role in controlling how plants react to different abiotic stressors. Studies on the HSF gene family have been reported in many species but have not yet undergone a thorough analysis in M. sativa. Therefore, in order to identify a more comprehensive set of HSF genes, from the genomic data, we identified 16 members of the MsHSF gene, which were unevenly distributed over six chromosomes. We also looked at their gene architectures and protein motifs, and phylogenetic analysis allowed us to divide them into 3 groups with a total of 15 subgroups. Along with these aspects, we then examined the physicochemical properties, subcellular localization, synteny analysis, GO annotation and enrichment, and protein interaction networks of amino acids. Finally, the analysis of 16 MsHSF genes’ expression levels across all tissues and under four abiotic stresses using publicly available RNA-Seq data revealed that these genes had significant tissue-specific expression. Moreover, the expression of most MsHSF genes increased dramatically under abiotic stress, further validating the critical function played by the MsHSF gene family in abiotic stress. These results provided basic information about MsHSF gene family and laid a foundation for further study on the biological role of MsHSF gene in response to stress in M. sativa.

Circulating Tumor DNA: Less Invasive, More Representative Method to Unveil the Genomic Landscape of Newly Diagnosed Multiple Myeloma Than Bone Marrow Aspirates.

Simple SummaryThe study of malignant plasma cell DNA genomics in multiple myeloma is a hot topic, which is mainly based on one-site bone marrow aspirates. In this study, we showed that circulating tumor DNA targeted next-generation sequencing analysis revealed a more comprehensive genomic architecture than bone marrow aspirates in newly diagnosed multiple myeloma. Circulating tumor DNA mutation in the transcriptional regulation pathway and DNA repair pathway were independent predictors of progression-free survival. ctDNA alterations correlated with prognosis and therapy response in newly diagnosed multiple myeloma.Multiple myeloma (MM) is highly heterogenous and dynamic in its genomic abnormalities. Capturing a representative image of these alterations is essential in understanding the molecular pathogenesis and progression of the disease but was limited by single-site invasive bone marrow (BM) biopsy-based genomics studies. We compared the mutational landscapes of circulating tumor DNA (ctDNA) and BM in 82 patients with newly diagnosed MM. A 413-gene panel was used in the sequencing. Our results showed that more than 70% of MM patients showed one or more genes with somatic mutations and at least half of the mutated genes were shared between ctDNA and BM samples. Compared to the BM samples, ctDNA exhibited more types of driver mutations in the shared driver genes, higher numbers of uniquely mutated genes and subclonal clusters, more translocation-associated mutations, and higher frequencies of mutated genes enriched in the transcriptional regulation pathway. Multivariate Cox analysis showed that age, ctDNA mutations in the transcriptional regulation pathway and DNA repair pathway were independent predictors of progression-free survival (PFS). Our results demonstrated sequencing of ctDNA provides more thorough information on the genomic instability and is a potential representative biomarker for risk stratification and in newly diagnosed MM than bone marrow.

Transcriptomic profiling reveals connection to lung autonomic nervous system dysfunction and inflammation in pulmonary artery hypertension rats

Abstract Background Pulmonary arterial hypertension (PAH) is a progressive disorder characterized by elevated pulmonary artery pressure, right ventricular failure, and premature death. Previous studies showed autonomic nervous system (ANS) such as sympathetic nerve overactivity or impaired parasympathetic activity was implicated in the pathogenesis of PAH. However, PA remodeling and molecular mechanisms involved in ANS in PAH remain unclear. Purpose To unravel transcriptional regulation, underlying genes and signaling pathways linking to ANS in PAH. Methods Publicly available two RNA sequencing (RNAseq) datasets and one single-cell RNA sequencing (scRNAseq) data of lung from rats with PAH induced by high-dose monocrotaline (MCT) were downloaded and reprocessed. For RNA-seq data, STAR and RSEM were used for read alignment and gene expression qualification. Differential expression genes (DEGs) were identified using DSEseq2, after which enrichment analyses were implemented by clusterProfiler. scRNAseq data were analyzed using Seurat. Results Two RNAseq datasets identified a total of 6,046 and 2,172 DEGs respectively, including 2,945 and 1,228 up-regulated and 3,101 and 944 down-regulated. Overlapped DEGs genes among datasets were 1,123 genes. Functional enrichment analysis of downregulated DEGs in one data pointed out many dysregulated pathways related to nerve system, such as nervous system process, neuron projection, integral component of synaptic membrane in MCT-induced PAH rat (Fig. 1a). Both RNAseq datasets supported upregulated DEGs were also involved in pathways related to ANS including humoral immune response and synapse pruning (P<0.001) (Fig. 1a). Intriguingly, all four genes (Trem2, C1qa, C1qb and C1qc) in synapse pruning were progressively upregulated with severe PAH condition (P<0.05) (Fig. 1b, c). And RT-PCR results confirmed the expression level of these four genes were remarkably increased in PAH rats (P<0.0001) (Fig. 1d). These four genes are known marker genes for macrophages from lung and microglial cell from brain. scRNAseq data of lungs demonstrated these four genes were significantly upregulated in alveolar macrophages from MCT-PAH rats (P<0.001) (Fig. 1e). However, the pathogenic mechanisms of these four genes in nerve fibers from lungs underlying PAH have never been studied and would be further explored using histological analysis in this study. Conclusions Integrated bioinformatic analysis indicated the association of lung nervous system with PAH. Further, we identified four genes implicated in nervous system, previously reported to link to immune system and provide novel insights into the mechanisms underlying relationship between lung ANS and inflammation in the pathogenesis of PAH. The pathogenetic mechanism of these four genes is being further investigated. Funding Acknowledgement Type of funding sources: Public grant(s) – National budget only. Main funding source(s): The National Natural Science Foundation of Chinathe Nature Science Foundation of Hubei Province

The moss-specific transcription factor PpERF24 positively modulates immunity against fungal pathogens in Physcomitrium patens.

APETALA2/ethylene response factors (AP2/ERFs) transcription factors (TFs) have greatly expanded in land plants compared to algae. In angiosperms, AP2/ERFs play important regulatory functions in plant defenses against pathogens and abiotic stress by controlling the expression of target genes. In the moss Physcomitrium patens, a high number of members of the ERF family are induced during pathogen infection, suggesting that they are important regulators in bryophyte immunity. In the current study, we analyzed a P. patens pathogen-inducible ERF family member designated as PpERF24. Orthologs of PpERF24 were only found in other mosses, while they were absent in the bryophytes Marchantia polymorpha and Anthoceros agrestis, the vascular plant Selaginella moellendorffii, and angiosperms. We show that PpERF24 belongs to a moss-specific clade with distinctive amino acids features in the AP2 domain that binds to the DNA. Interestingly, all P. patens members of the PpERF24 subclade are induced by fungal pathogens. The function of PpERF24 during plant immunity was assessed by an overexpression approach and transcriptomic analysis. Overexpressing lines showed increased defenses to infection by the fungal pathogens Botrytis cinerea and Colletotrichum gloeosporioides evidenced by reduced cellular damage and fungal biomass compared to wild-type plants. Transcriptomic and RT-qPCR analysis revealed that PpERF24 positively regulates the expression levels of defense genes involved in transcriptional regulation, phenylpropanoid and jasmonate pathways, oxidative burst and pathogenesis-related (PR) genes. These findings give novel insights into potential mechanism by which PpERF24 increases plant defenses against several pathogens by regulating important players in plant immunity.

CRHR1 mediates the transcriptional expression of pituitary hormones and their receptors under hypoxia.

Hypothalamus-pituitary-adrenal (HPA) axis plays critical roles in stress responses under challenging conditions such as hypoxia, via regulating gene expression and integrating activities of hypothalamus-pituitary-targets cells. However, the transcriptional regulatory mechanisms and signaling pathways of hypoxic stress in the pituitary remain to be defined. Here, we report that hypoxia induced dynamic changes in the transcription factors, hormones, and their receptors in the adult rat pituitary. Hypoxia-inducible factors (HIFs), oxidative phosphorylation, and cAMP signaling pathways were all differentially enriched in genes induced by hypoxic stress. In the pituitary gene network, hypoxia activated c-Fos and HIFs with specific pituitary transcription factors (Prop1), targeting the promoters of hormones and their receptors. HIF and its related signaling pathways can be a promising biomarker during acute or constant hypoxia. Hypoxia stimulated the transcription of marker genes for microglia, chemokines, and cytokine receptors of the inflammatory response. Corticotropin-releasing hormone receptor 1 (CRHR1) mediated the transcription of Pomc, Sstr2, and Hif2a, and regulated the function of HPA axis. Together with HIF, c-Fos initiated and modulated dynamic changes in the transcription of hormones and their receptors. The receptors were also implicated in the regulation of functions of target cells in the pituitary network under hypoxic stress. CRHR1 played an integrative role in the hypothalamus-pituitary-target axes. This study provides new evidence for CRHR1 involved changes of hormones, receptors, signaling molecules and pathways in the pituitary induced by hypoxia.

N6-methyladenosine methylation analysis reveals transcriptome-wide expression response to salt stress in rice roots

The N6-methyladenosine (m6A) modifications play crucial regulatory roles in plant yield development, nutritional growth and stress adaption. Rice ( Oryza sativa ) is relatively sensitive to salt stress among cereal crops, while its regulatory mechanism conferring salt tolerance is still unclear. Here, we performed transcriptome-wide m6A analysis in roots of rice cultivar Nipponbare under salt conditions (0 and 75 mM NaCl). After salt treatments for 4 d, a total of 18,115 m6A peak callings of 12,516 genes and 22,052 peaks of 14,235 genes were identified under salt and control conditions, respectively. These m6A peaks were mostly enriched in CDS and 3′UTR regions. Totally, 346 m6A differentially peaks (DPs) and 1522 differentially expressed genes (DEGs) were identified in response to salt stress, which were mainly enriched in "MAPK signaling module" and "metal ion binding". Among them, the m6A related genes including 2 writers ( MTB and TRM4A ), 1 reader ( CPSF30 ) and 1 eraser ( ALKBH9A ). The m6A modified salt-responding genes including HAK4 , CIPK06 , RBOHH , Myb10 , ERF067 and other transcription factors might participate in transcriptional regulation or multiple Ca 2+ , ABA and ROS signal pathways. The m6A level and gene expression of these genes were further confirmed by m6A-IP-qPCR and qRT-PCR, respectively, which showed a high consistency with the m6A-seq and transcriptome data. This study provides epigenetic responses of molecular regulatory mechanisms of salt tolerance in rice. ● Identified 346 salt-induced m6A DPs and 1522 DEGs in rice roots. ● Most m6A DPs were enriched in stop codon and CDS related to gene length. ● Proposed the roles of m6A modified salt-responding genes in transcriptional regulation and multiple signal pathways. ● Performed m6A-IP-qPCR and qRT-PCR experiments to validate the accuracy of data.