Post-transcriptional Regulatory Networks Research Articles

Multifaceted regulation of siderophore synthesis by multiple regulatory systems in Shewanella oneidensis

Siderophore-dependent iron uptake is a mechanism by which microorganisms scavenge and utilize iron for their survival, growth, and many specialized activities, such as pathogenicity. The siderophore biosynthetic system PubABC in Shewanella can synthesize a series of distinct siderophores, yet how it is regulated in response to iron availability remains largely unexplored. Here, by whole genome screening we identify TCS components histidine kinase (HK) BarA and response regulator (RR) SsoR as positive regulators of siderophore biosynthesis. While BarA partners with UvrY to mediate expression of pubABC post-transcriptionally via the Csr regulatory cascade, SsoR is an atypical orphan RR of the OmpR/PhoB subfamily that activates transcription in a phosphorylation-independent manner. By combining structural analysis and molecular dynamics simulations, we observe conformational changes in OmpR/PhoB-like RRs that illustrate the impact of phosphorylation on dynamic properties, and that SsoR is locked in the ‘phosphorylated’ state found in phosphorylation-dependent counterparts of the same subfamily. Furthermore, we show that iron homeostasis global regulator Fur, in addition to mediating transcription of its own regulon, acts as the sensor of iron starvation to increase SsoR production when needed. Overall, this study delineates an intricate, multi-tiered transcriptional and post-transcriptional regulatory network that governs siderophore biosynthesis.

Regulation and tumor-suppressive function of the miR-379/miR-656 (C14MC) cluster in cervical cancer.

Cervical cancer (CC) is a key contributor to cancer-related mortality in several countries. The identification of molecular markers and the underlying mechanism may help improve CC management. We studied the regulation and biological function of the chromosome 14 microRNA cluster (C14MC; miR-379/miR-656) in CC. Most C14MC members exhibited considerably lower expression in CC tissues and cell lines in The Cancer Genome Atlas (TCGA) cervical squamous cell carcinoma and endocervical adenocarcinoma patient cohorts. Bisulfite Sanger sequencing revealed hypermethylation of the C14MC promoter in CC tissues and cell lines. 5-aza-2 deoxy cytidine treatment reactivated expression of the C14MC members. We demonstrated that C14MC is a methylation-regulated miRNA cluster via artificial methylation and luciferase reporter assays. C14MC downregulation correlated with poor overall survival and may promote metastasis. C14MC activation via the lentiviral-based CRISPRa approach inhibited growth, proliferation, migration, and invasion; enhanced G2/M arrest; and induced senescence. Post-transcriptional regulatory network analysis of C14MC transcriptomic data revealed enrichment of key cancer-related pathways, such as metabolism, the cell cycle, and phosphatidylinositol 3-kinase (PI3K)-AKT signaling. Reduced cell proliferation, growth, migration, invasion, and senescence correlated with the downregulation of active AKT, MYC, and cyclin E1 (CCNE1) and the overexpression of p16, p21, and p27. We showed that C14MC miRNA activation increases reactive oxygen species (ROS) levels, intracellular Ca2+ levels, and lipid peroxidation rates, and inhibits epithelial-mesenchymal transition (EMT). C14MC targets pyruvate dehydrogenase kinase-3 (PDK3) according to the luciferase reporter assay. PDK3 is overexpressed in CC and is inversely correlated with C14MC. Both miR-494-mimic transfection and C14MC activation inhibited PDK3 expression. Reduced glucose uptake and lactate production, and upregulation of PDK3 upon C14MC activation suggest the potential role of these proteins in metabolic reprogramming. Finally, we showed that C14MC activation may inhibit EMT signaling. Thus, C14MC is a tumor-suppressive and methylation-regulated miRNA cluster in CC. Reactivation of C14MC can be useful in the management of CC.

Computational study of the motor neuron protein KIF5A to identify nsSNPs, bioactive compounds, and its key regulators.

Introduction: Kinesin family member 5A (KIF5A) is a motor neuron protein expressed in neurons and involved in anterograde transportation of organelles, proteins, and RNA. Variations in the KIF5A gene that interfere with axonal transport have emerged as a distinguishing feature in several neurodegenerative disorders, including hereditary spastic paraplegia (HSP10), Charcot-Marie-Tooth disease type 2 (CMT2), and Amyotrophic Lateral Sclerosis (ALS). Methods: In this study, we implemented a computational structural and systems biology approach to uncover the role of KIF5A in ALS. Using the computational structural biology method, we explored the role of non-synonymous Single Nucleotide Polymorphism (nsSNPs) in KIF5A. Further, to identify the potential inhibitory molecule against the highly destabilizing structure variant, we docked 24 plant-derived phytochemicals involved in ALS. Results: We found KIF5AS291F variant showed the most structure destabilizing behavior and the phytocompound "epigallocatechin gallate" showed the highest binding affinity (-9.0Kcal/mol) as compared to wild KIF5A (-8.4Kcal/mol). Further, with the systems biology approach, we constructed the KIF5A protein-protein interaction (PPI) network to identify the associated Kinesin Families (KIFs) proteins, modules, and their function. We also constructed a transcriptional and post-transcriptional regulatory network of KIF5A. With the network topological parameters of PPIN (Degree, Bottleneck, Closeness, and MNC) using CytoHubba and computational knock-out experiment using Network Analyzer, we found KIF1A, 5B, and 5C were the significant proteins. The functional modules were highly enriched with microtubule motor activity, chemical synaptic transmission in neurons, GTP binding, and GABA receptor activity. In regulatory network analysis, we found KIF5A post-transcriptionally down-regulated by miR-107 which is further transcriptionally up-regulated by four TFs (HIF1A, PPARA, SREBF1, and TP53) and down-regulated by three TFs (ZEB1, ZEB2, and LIN28A). Discussion: We concluded our study by finding a crucial variant of KIF5A and its potential therapeutic target (epigallocatechin gallate) and KIF5A associated significant genes with important regulators which could decrypt the novel therapeutics in ALS and other neurodegenerative diseases.

Systematic analysis of various RNA transcripts and construction of biological regulatory networks at the post-transcriptional level for chronic obstructive pulmonary disease

BackgroundAlthough chronic inflammation, oxidative stress, airway remodeling, and protease-antiprotease imbalance have been implicated in chronic obstructive pulmonary disease (COPD), the exact pathogenesis is still obscure. Gene transcription and post-transcriptional regulation have been taken into account as key regulators of COPD occurrence and development. Identifying the hub genes and constructing biological regulatory networks at the post-transcriptional level will help extend current knowledge on COPD pathogenesis and develop potential drugs.MethodsAll lung tissues from non-smokers (n = 6), smokers without COPD (smokers, n = 7), and smokers with COPD (COPD, n = 7) were collected to detect messenger RNA (mRNA), microRNA (miRNA), circular RNA (circRNA), and long non-coding RNA (lncRNA) expression and identify the hub genes. Biological regulatory networks were constructed at the post-transcriptional level, including the RNA-binding protein (RBP)-hub gene interaction network and the competitive endogenous RNA (ceRNA) network. In addition, we assessed the composition and abundance of immune cells in COPD lung tissue and predicted potential therapeutic drugs for COPD. Finally, the hub genes were confirmed at both the RNA and protein levels.ResultsAmong the 20 participants, a total of 121169 mRNA transcripts, 1871 miRNA transcripts, 4244 circRNA transcripts, and 122130 lncRNA transcripts were detected. There were differences in the expression of 1561 mRNAs, 48 miRNAs, 33 circRNAs, and 545 lncRNAs between smokers and non-smokers, as well as 1289 mRNAs, 69 miRNAs, 32 circRNAs, and 433 lncRNAs between smokers and COPD patients. 18 hub genes were identified in COPD. TGF-β signaling and Wnt/β-catenin signaling may be involved in the development of COPD. Furthermore, the circRNA/lncRNA-miRNA-mRNA ceRNA networks and the RBP-hub gene interaction network were also constructed. Analysis of the immune cell infiltration level revealed that M2 macrophages and activated NK cells were increased in COPD lung tissues. Finally, we identified that the ITK inhibitor and oxybutynin chloride may be effective in treating COPD.ConclusionsWe identified several novel hub genes involved in COPD pathogenesis. TGF-β signaling and Wnt/β-catenin signaling were the most dysregulated pathways in COPD patients. Our study constructed post-transcriptional biological regulatory networks and predicted small-molecule drugs for the treatment of COPD, which enhanced the existing understanding of COPD pathogenesis and suggested an innovative direction for the therapeutic intervention of the disease.

Weighted gene co-expression network analysis identified hub genes critical to fatty acid composition in Gushi chicken breast muscle

BackgroundThe composition and content of fatty acids in the breast muscle are important factors influencing meat quality. In this study, we investigated the fatty acid composition and content in the breast muscle of Gushi chickens at different developmental stages (14 weeks, 22 weeks, and 30 weeks). Additionally, we utilized transcriptomic data from the same tissue and employed WGCNA and module identification methods to identify key genes associated with the fatty acid composition in Gushi chicken breast muscle and elucidate their regulatory networks.ResultsAmong them, six modules (blue, brown, green, light yellow, purple, and red modules) showed significant correlations with fatty acid content and metabolic characteristics. Enrichment analysis revealed that these modules were involved in multiple signaling pathways related to fatty acid metabolism, including fatty acid metabolism, PPAR signaling pathway, and fatty acid biosynthesis. Through analysis of key genes, we identified 136 genes significantly associated with fatty acid phenotypic traits. Protein–protein interaction network analysis revealed that nine of these genes were closely related to fatty acid metabolism. Additionally, through correlation analysis of transcriptome data, we identified 51 key ceRNA regulatory networks, including six central genes, 7 miRNAs, and 28 lncRNAs.ConclusionThis study successfully identified key genes closely associated with the fatty acid composition in Gushi chicken breast muscle, as well as their post-transcriptional regulatory networks. These findings provide new insights into the molecular regulatory mechanisms underlying the flavor characteristics of chicken meat and the composition of fatty acids in the breast muscle.

Key Risk Genes Identified From the Postmortem Brain of Patients With Major Depressive Disorder and Their Potential Clinical Applications.

Major depressive disorder (MDD) is a type of emotional dysfunction, and its pathogenesis has not been fully elucidated. Specifically, the key molecules in depression-related brain regions involved in this disease and their contributions to this disease are currently unclear. GSE53987 and GSE54568 were selected from the Gene Expression Omnibus database. The data were standardized to identify the common differentially expressed genes (DEGs) in the cortex of MDD patients in the 2 datasets. The DEGs were subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses. The STRING database was used to build protein-protein interaction networks, and the cytoHubba plugin was used to identify hub genes. Furthermore, we selected another blood transcriptome dataset that included 161 MDD and 169 control samples to explore the changes in the screened hub genes. Mice were subjected to 4 weeks of chronic unpredictable mild stress to establish an animal model of depression, and the expression of these hub genes in tissues of the prefrontal cortex was then detected by quantitative real time polymerase chain reaction (qRT-PCR). We subsequently predicted the possible posttranscriptional regulatory networks and traditional Chinese medicine according to the hub genes using a few online databases. The analysis identified 147 upregulated genes and 402 downregulated genes were identified in the cortex of MDD patients compared with that of the controls. Enrichment analyses revealed that DEGs were predominantly enriched in synapse-related cell functions, linoleic acid metabolism, and other pathways. Protein-protein interaction analysis identified 20 hub genes based on the total score. The changes in KDM6B, CUX2, NAAA, PHKB, NFYA, GTF2H1, CRK, CCNG2, ACER3, and SLC4A2 in the peripheral blood of MDD patients were consistent with those in the brain. Furthermore, the prefrontal cortex of mice with depressive-like behaviors showed significantly increased Kdm6b, Aridb1, Scaf11, and Thoc2 expression and decreased Ccng2 expression compared with that of normal mice, which was consistent with the results found for the human brain. Potential therapeutic candidates, such as citron, fructus citri, leaves of Panax Notoginseng, sanchi flower, pseudoginseng, and dan-shen root, were selected via traditional Chinese medicine screening. This study identified several novel hub genes in specific brain regions involved in the pathogenesis of MDD, which may not only deepen our understanding of depression but may also provide new ideas for its diagnosis and treatment.

The SARS-CoV-2 UTR's Intrudes Host RBP's and Modulates Cellular Splicing.

SARS-CoV-2 is a novel coronavirus that causes a potentially fatal respiratory disease known as coronavirus disease (COVID-19) and is responsible for the ongoing pandemic with increasing mortality. Understanding the host-virus interaction involved in SARS-CoV-2 pathophysiology will enhance our understanding of the mechanistic basis of COVID-19 infection. The characterization of post-transcriptional gene regulatory networks, particularly pre-mRNA splicing, and the identification and characterization of host proteins interacting with the 5' and 3'UTRs of SARS-CoV-2 will improve our understanding of post-transcriptional gene regulation during SARS-CoV-2 pathogenesis. Here, we demonstrate that either SARS-CoV-2 infection or exogenous overexpression of the 5' and 3'UTRs of the viral genomic RNAs, results in reduced mRNA levels possibly due to modulation of host cell pre-mRNA splicing. Further, we have investigated the potential RNA-binding proteins interacting with the 5' and 3'UTRs, using in-silico approaches. Our results suggest that 5' and 3'UTRs indeed interact with many RNA-binding proteins. Our results provide a primer for further investigations into the UTR-mediated regulation of splicing and related molecular mechanisms in host cells.

A Data-Driven Signaling Network Inference Approach for Phosphoproteomics.

Proteins are rapidly and dynamically post-transcriptionally modified as cells respond to changes in their environment. For example, protein phosphorylation is mediated by kinases while dephosphorylation is mediated by phosphatases. Quantifying and predicting interactions between kinases, phosphatases, and target proteins over time will aid the study of signaling cascades under a variety of environmental conditions. Here, we describe methods to statistically analyze label-free phosphoproteomic data and infer posttranscriptional regulatory networks over time. We provide an R-based method that can be used to normalize and analyze label-free phosphoproteomic data using variance stabilizing normalization and a linear mixed model across multiple time points and conditions. We also provide a method to infer regulator-target interactions over time using a discretization scheme followed by dynamic Bayesian modeling computations to validate our conclusions. Overall, this pipeline is designed to perform functional analyses and predictions of phosphoproteomic signaling cascades.

Species-specific regulatory pathways of small RNAs play sophisticated roles in flower development in Dimocarpus longan Lour.

Flower development plays vital role in horticultural plants. Post-transcriptional regulation via small RNAs is important for plant flower development. To uncover post-transcriptional regulatory networks during the flower development in Dimocarpus longan Lour. ‘Shixia’, an economically important fruit crop in subtropical regions, we collected and analyzed sRNA deep-sequencing datasets and degradome libraries Apart from identifying miRNAs and phased siRNA generating loci (PHAS loci), 120 hairpin loci, producing abundant sRNAs, were identified by in-house protocols. Our results suggested that 56 miRNA-target pairs, 22 21-nt-PHAS loci, and 111 hairpin loci are involved in post-transcriptional gene silencing during longan reproductive development. Lineage-specific or species-specific post-transcriptional regulatory modules have been unveiled, including miR482-PHAS and miRN15. miR482-PHAS might be involved in longan flower development beyond their conserved roles in plant defense, and miRN15 is a novel miRNA likely associated with a hairpin locus (HPL-056) to regulate strigolactone receptor gene DWARF14 (D14) and the biogenesis of phasiRNAs from D14. These small RNAs are enriched in flower buds, suggesting they are likely involved in post-transcriptional regulatory networks essential for longan flower development via the strigolactone signaling pathway.

Alterations in hepatic mitotic and cell cycle transcriptional networks during the metabolic switch in broiler chicks.

During embryonic life, chicks mainly derive energy from hepatic oxidation of yolk lipids. After hatch, chicks must rely on carbohydrate-rich feed to obtain energy. This requires an abrupt and intensive switch of metabolic processes, particularly in the liver. We recently identified a number of transcriptional and post-transcriptional regulatory networks that work concordantly to tune metabolic processes during the metabolic switch. Here, we used delayed feeding post-hatch (48h) to impede the metabolic switch in broilers. We used RNA-seq to identify hepatic transcriptome differences between late stage embryos (E18) and two-day-old chicks (D2), which were either fed-from-hatch (FED) or not fed (DLY). Between FED and E18, 2,430 genes were differentially expressed (fold-change≥ 2; FDR p-value 0.05), of these 1,237 were downregulated in FED birds and 1,193 were upregulated. Between DLY and E18, 1979 genes were differentially expressed, of these 1,043 were downregulated and 936 were upregulated in DLY birds. Between DLY and FED, 880 genes were differentially expressed, of these 543 were downregulated and 337 were upregulated in DLY birds. We found that in addition to disturbances in a number of metabolic pathways, unfed chicks had a widespread suppression of gene networks associated with cell proliferation, cell cycle progression and mitosis. Expression patterns suggest that hepatocytes of delayed-fed birds have abnormal mitosis and increased polyploidization. This suggests that post-hatch feed consumption maintains the rate and integrity of liver growth immediately, which in turn, likely helps facilitate the appropriate programming of hepatic metabolic networks.

Identification, Phylogeny, Divergence, Structure, and Expression Analysis of A20/AN1 Zinc Finger Domain Containing Stress-Associated Proteins (SAPs) Genes in Jatropha curcas L.

Jatropha is a small woody perennial biofuel-producing shrub. Stress-associated proteins (SAPs) are novel stress regulatory zinc-finger proteins and are mainly associated with tolerance against various environmental abiotic stresses in Jatropha. In the present study, the JcSAP gene family were analyzed comprehensively in Jatropha curcas and 11 JcSAP genes were identified. Phylogenetic analysis classified the JcSAP genes into four groups based on sequence similarity, similar gene structure features, conserved A20 and/or AN1 domains, and their responsive motifs. Moreover, the divergence analysis further evaluated the evolutionary aspects of the JcSAP genes with the predicted time of divergence from 9.1 to 40 MYA. Furthermore, a diverse range of cis-elements including light-responsive elements, hormone-responsive elements, and stress-responsive elements were detected in the promoter region of JcSAP genes while the miRNA target sites predicted the regulation of JcSAP genes via a candid miRNA mediated post-transcriptional regulatory network. In addition, the expression profiles of JcSAP genes in different tissues under stress treatment indicated that many JcSAP genes play functional developmental roles in different tissues, and exhibit significant differential expression under stress treatment. These results collectively laid a foundation for the functional diversification of JcSAP genes.

Genomic Variants Disrupt miRNA-mRNA Regulation.

Micro RNA (miRNA) and its regulatory effect on messenger RNA (mRNA) gene expression are a major focus in cancer research. Disruption in the normal miRNA-mRNA regulation network can result in serious cascading biological repercussions. In this study, we curated miRNA-related variants from major genomic consortiums and thoroughly evaluated how these variants could exert their effects by cross-validating with independent functional knowledge bases. Nearly all known variants (more than 664 million) categorized by type (germline, somatic, epigenetic) were mapped to the genomic regions involved in miRNA-mRNA binding (miRNA seeds and miRNA-mRNA 3'-UTR binding sequence). Subsets of miRNA-related variants supported by additional functional evidence, such as expression Quantitative Trait Loci (eQTL) and Genome-Wide Association Study (GWAS), were identified and scrutinized. Our results show that variants in miRNA seeds can substantially alter the composition of a miRNA's target mRNA set. Various functional analyses converged to reveal a post-transcriptional complex regulatory network where miRNA, eQTL, and RNA-binding protein intertwined to disseminate the impact of genomic variants. These results may potentially explain how certain variants affect disease/trait risks in genome wide association studies.

Esrrb Regulates Specific Feed-Forward Loops to Transit From Pluripotency Into Early Stages of Differentiation.

Characterization of pluripotent states, in which cells can both self-renew or differentiate, with the irreversible loss of pluripotency, are important research areas in developmental biology. Although microRNAs (miRNAs) have been shown to play a relevant role in cellular differentiation, the role of miRNAs integrated into gene regulatory networks and its dynamic changes during these early stages of embryonic stem cell (ESC) differentiation remain elusive. Here we describe the dynamic transcriptional regulatory circuitry of stem cells that incorporate protein-coding and miRNA genes based on miRNA array expression and quantitative sequencing of short transcripts upon the downregulation of the Estrogen Related Receptor Beta (Esrrb). The data reveals how Esrrb, a key stem cell transcription factor, regulates a specific stem cell miRNA expression program and integrates dynamic changes of feed-forward loops contributing to the early stages of cell differentiation upon its downregulation. Together these findings provide new insights on the architecture of the combined transcriptional post-transcriptional regulatory network in embryonic stem cells.

Identification of an interactome network between lncRNAs and miRNAs in thyroid cancer reveals SPTY2D1-AS1 as a new tumor suppressor

Thyroid cancer is the most common primary endocrine malignancy in adults and its incidence is rapidly increasing. Long non-coding RNAs (lncRNAs), generally defined as RNA molecules longer than 200 nucleotides with no protein-encoding capacity, are highly tissue-specific molecules that serve important roles in gene regulation through a variety of different mechanisms, including acting as competing endogenous RNAs (ceRNAs) that ‘sponge’ microRNAs (miRNAs). In the present study, using an integrated approach through RNA-sequencing of paired thyroid tumor and non-tumor samples, we have identified an interactome network between lncRNAs and miRNAs and examined the functional consequences in vitro and in vivo of one of such interactions. We have identified a likely operative post-transcriptional regulatory network in which the downregulated lncRNA, SPTY2D1-AS1, is predicted to target the most abundant and upregulated miRNAs in thyroid cancer, particularly miR-221, a well-known oncomiRNA in cancer. Indeed, SPTY2D1-AS1 functions as a potent tumor suppressor in vitro and in vivo, it is downregulated in the most advanced stages of human thyroid cancer, and it seems to block the processing of the primary form of miR-221. Overall, our results link SPTY2D1-AS1 to thyroid cancer progression and highlight the potential use of this lncRNA as a therapeutic target of thyroid cancer.

Genome-wide identification, classification and expression analysis of the JmjC domain-containing histone demethylase gene family in Jatropha curcas L.

JmjC domain-containing proteins, an important family of histone lysine demethylase, play significant roles in maintaining the homeostasis of histone methylation. In this study, we comprehensively analyzed the JmjC domain-containing gene family in Jatropha curcas and found 20 JmjC domain-containing genes (JcJMJ genes). Phylogenetic analysis revealed that these JcJMJ genes can be classified into five major subgroups, and genes in each subgroup had similar motif and domain composition. Cis-regulatory element analysis showed that the number and types of cis-regulatory elements owned by the promoter of JcJMJ genes in different subgroup were significantly different. Moreover, miRNA target prediction result revealed a complicated miRNA-mediated post-transcriptional regulatory network, in which JcJMJ genes were regulated by different numbers and types of miRNAs. Further analysis of the tissue and stress expression profiles showed that many JcJMJ genes had tissue and stress expression specificity. All these results provided valuable information for understanding the evolution of JcJMJ genes and the complex transcriptional and post transcriptional regulation involved, and laid the foundation for further functional analysis of JcJMJ genes.

Regulatory Role of Circadian Clocks on ABA Production and Signaling, Stomatal Responses, and Water-Use Efficiency under Water-Deficit Conditions.

Plants deploy molecular, physiological, and anatomical adaptations to cope with long-term water-deficit exposure, and some of these processes are controlled by circadian clocks. Circadian clocks are endogenous timekeepers that autonomously modulate biological systems over the course of the day–night cycle. Plants’ responses to water deficiency vary with the time of the day. Opening and closing of stomata, which control water loss from plants, have diurnal responses based on the humidity level in the rhizosphere and the air surrounding the leaves. Abscisic acid (ABA), the main phytohormone modulating the stomatal response to water availability, is regulated by circadian clocks. The molecular mechanism of the plant’s circadian clock for regulating stress responses is composed not only of transcriptional but also posttranscriptional regulatory networks. Despite the importance of regulatory impact of circadian clock systems on ABA production and signaling, which is reflected in stomatal responses and as a consequence influences the drought tolerance response of the plants, the interrelationship between circadian clock, ABA homeostasis, and signaling and water-deficit responses has to date not been clearly described. In this review, we hypothesized that the circadian clock through ABA directs plants to modulate their responses and feedback mechanisms to ensure survival and to enhance their fitness under drought conditions. Different regulatory pathways and challenges in circadian-based rhythms and the possible adaptive advantage through them are also discussed.

The core and accessory Hfq interactomes across Pseudomonas aeruginosa lineages

The major RNA-binding protein Hfq interacts with mRNAs, either alone or together with regulatory small noncoding RNAs (sRNAs), affecting mRNA translation and degradation in bacteria. However, studies tend to focus on single reference strains and assume that the findings may apply to the entire species, despite the important intra-species genetic diversity known to exist. Here, we use RIP-seq to identify Hfq-interacting RNAs in three strains representing the major phylogenetic lineages of Pseudomonas aeruginosa. We find that most interactions are in fact not conserved among the different strains. We identify growth phase-specific and strain-specific Hfq targets, including previously undescribed sRNAs. Strain-specific interactions are due to different accessory gene sets, RNA abundances, or potential context- or sequence- dependent regulatory mechanisms. The accessory Hfq interactome includes most mRNAs encoding Type III Secretion System (T3SS) components and secreted toxins in two strains, as well as a cluster of CRISPR guide RNAs in one strain. Conserved Hfq targets include the global virulence regulator Vfr and metabolic pathways involved in the transition from fast to slow growth. Furthermore, we use rGRIL-seq to show that RhlS, a quorum sensing sRNA, activates Vfr translation, thus revealing a link between quorum sensing and virulence regulation. Overall, our work highlights the important intra-species diversity in post-transcriptional regulatory networks in Pseudomonas aeruginosa.

Circ_0004354 might compete with circ_0040039 to induce NPCs death and inflammatory response by targeting miR-345-3p-FAF1/TP73 axis in intervertebral disc degeneration.

The abnormal function of nucleus pulposus cells (NPCs) plays a crucial role in the pathogenesis of intervertebral disc degeneration (IVDD). Recent studies have demonstrated that circular RNAs (circRNAs) are involved in the pathological process of IVDD by regulating NPCs' function. Nevertheless, the investigation on circRNA-circRNA interaction has not yet been reported. Here, we identified the top upregulated circ_0040039 and circ_0004354 in IVDD, derived from the syntrophin beta 2 gene but had different degrees of biological functions. Accumulating studies have reported PANoptosis is composed of apoptosis, pyroptosis, and necroptosis. Based on this, we think there should be a new pro-inflammatory cell death PAoptosis in the form of apoptosis and pyroptosis. Circ_0004354 might compete with circ_0040039 to induce the development of IVDD by modulating miR-345-3p-FAF1/TP73 axis-mediated PAoptosis, inflammatory response, growth inhibition, and ECM degradation of NPCs. Thus, these findings offer a novel insight into the circRNAs-mediated posttranscriptional regulatory network in IVDD, contributing to further clarification of the pathological mechanism of IVDD to develop a promising therapeutic target for IVDD diseases.

Abstract PR-02: Characterization of pseudogene expression and immune infiltration in breast cancer subtypes for African American and Caucasian patients

Abstract Introduction: Breast Cancer (BC) occurrence, especially triple-negative BC (TNBC), is relatively high among African American women. In addition, Caucasian BC mortality has decreased much faster than African American BC mortality even with advancing medicines over the last few years. Researchers have reported socioeconomic and biological factors that may contribute to African American BC disparity, but concrete genetic biomarkers have not been elucidated. In this study, we conducted a systematic bioinformatics analysis to compare the gene expression pattern and examine the immunological differences in Black and White patients for different breast cancer subtypes. Methods: Using the public portal of The Cancer Genome Atlas (TCGA), Level 3 RNA-Seq expression data and corresponding clinical data for 183 African American and 764 Caucasian BC patients were downloaded. These patients were then categorized by four subtypes (luminal A, luminal B, HER2-enriched, and TNBC) and subsequently analyzed for gene expression and immune infiltration association. The average expression in RPKM (reads per kilobase of transcript per million mapped reads) of transcribed genes for each race and subtype was calculated using a custom R script. For each subtype, the top 30 genes with the highest expression in each race were selected and cross-checked for any overlap. Association between their expression and tumor infiltrates for six immune cell types (B cell, Dendritic cell, Macrophage, Neutrophil, T-cell CD4, T-cell CD8) was assessed using TIMER2.0. Results: For all subtypes, the majority of the top 30 genes overlapped between both races. In particular, 24 out of the top 30 genes (80%) were common between both races for TNBC. Five genes (ADAM6, CD74, LOC96610, B2M, HLA-B) were reported to show strong positive correlation between their expression and immune infiltration. Interestingly, out of the unique gene set for both races in the TNBC subtype, none of them showed strong positive correlation. ADAM6, a transcribed pseudogene, was identified to be positively correlated to immune infiltration in African American BC patients uniquely across all subtypes excluding TNBC in our study. Conclusion: Our computation result provides a candidate list of unique and common genes having high expression for four BC subtypes in African American and Caucasian patients. Our work also suggests that expressed pseudogenes could serve as potential biomarkers for BC patients. We plan to continue to investigate biomarkers in BC subtypes and immune response over post-transcriptional regulatory network in the context of breast cancer etiology. Citation Format: Eric W. Li, Yongsheng Bai. Characterization of pseudogene expression and immune infiltration in breast cancer subtypes for African American and Caucasian patients [abstract]. In: Proceedings of the AACR Virtual Conference: 14th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2021 Oct 6-8. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2022;31(1 Suppl):Abstract nr PR-02.

Comprehensive transcriptome characterization of Grus japonensis using PacBio SMRT and Illumina sequencing

The red-crowned crane (Grus japonensis) is an endangered species distributed across southeast Russia, northeast China, Korea, and Japan. Here, we sequenced for the first time the full-length unreferenced transcriptome of red-crowned crane mixed samples using a PacBio Sequel platform. A total of 359,136 circular consensus sequences (CCS) were obtained via clustering to remove redundancy. A total of 303,544 full-length non-chimeric sequences were identified by judging whether CCS contained 5′ and 3′ adapters, and the poly(A) tail. Eight samples were sequenced using Illumina, and PacBio sequencing data were corrected according to the collected Illumina data to obtain more accurate full-length transcripts. A total of 4,100 long non-coding RNAs, 13,115 simple sequences repeat loci and 29 transcription factor families were identified. The expression of lncRNAs and TFs in pancreas was lowest comparing with other tissues. Many enriched immune-related transmission pathways (MHC and IL receptors) were identified in the spleen. This study will contribute to a better understanding of the gene structure and post-transcriptional regulatory network, and provide references for future studies on red-crowned cranes.