Genome-wide RNA-seq Analysis Research Articles

Genome-wide RNA-seq, DNA methylation and small RNA-seq analysis unraveled complex gene regulatory networks in psoriasis pathogenesis

Psoriasis is a complex inflammatory skin disease characterized by reversible albeit relapsing red scaly plaques in the skin of a patient. In addition to the genetic predisposition, involvement of epigenetic and non-coding RNAs have also been liked with the disease. Nevertheless, any comprehensive study involving transcriptomic, small-RNA and DNA methylation at the genomic level from same patients is lacking. To investigate the complex regulation of molecular pathways in psoriasis, we carried out multi-omics integrative analysis of RNA-sequencing, small RNA-sequencing and DNA methylation profiling from the psoriatic and adjacent normal skin tissues. Our multi-omics analysis identified the genes and biological processes regulated either independently or in combination by DNA methylation and microRNAs. We identified miRNAs that specifically regulated keratinocyte hyper-proliferation, and cell cycle progression and checkpoint signaling in psoriasis. On contrary, DNA methylation was found to be more predominant in regulating immune and inflammatory responses, another causative factor in psoriasis pathogenesis. Many characteristic pathways in psoriasis e.g., Th17 cell differentiation and JAK-STAT signaling, were found to be regulated by both miRNAs and DNA methylation. We carried out functional characterization of a downregulated miRNA hsa-let-7c-5p, predicted to target upregulated genes in psoriasis involved in cell cycle processes, Th17 cell differentiation and JAK-STAT signaling pathways. Overexpression of hsa-let-7c-5p in keratinocytes caused the downregulation of its target genes, resulting in reduced cell proliferation and migration rates, demonstrating potential of miRNAs in regulating psoriasis pathogenesis. In conclusion, our findings identified distinct and shared gene-networks regulated by DNA methylation and miRNAs of a complex disease with reversible phenotype.

Genome-Wide Association and RNA-Seq Analyses Reveal a Potential Candidate Gene Related to Oil Content in Soybean Seeds.

Soybean is a crucial crop globally, serving as a significant source of unsaturated fatty acids and protein in the human diet. However, further enhancements are required for the related genes that regulate soybean oil synthesis. In this study, 155 soybean germplasms were cultivated under three different environmental conditions, followed by phenotypic identification and genome-wide association analysis using simplified sequencing data. Genome-wide association analysis was performed using SLAF-seq data. A total of 36 QTLs were significantly associated with oil content (-log10(p) > 3). Out of the 36 QTLs associated with oil content, 27 exhibited genetic overlap with previously reported QTLs related to oil traits. Further transcriptome sequencing was performed on extreme high-low oil soybean varieties. Combined with transcriptome expression data, 22 candidate genes were identified (|log2FC| ≥ 3). Further haplotype analysis of the potential candidate genes showed that three potential candidate genes had excellent haplotypes, including Glyma.03G186200, Glyma.09G099500, and Glyma.18G248900. The identified loci harboring beneficial alleles and candidate genes likely contribute significantly to the molecular network's underlying marker-assisted selection (MAS) and oil content.

Epigenetically active chromatin in neonatal iWAT reveals GABPα as a potential regulator of beige adipogenesis.

Thermogenic beige adipocytes, which dissipate energy as heat, are found in neonates and adults. Recent studies show that neonatal beige adipocytes are highly plastic and contribute to >50% of beige adipocytes in adults. Neonatal beige adipocytes are distinct from recruited beige adipocytes in that they develop independently of temperature and sympathetic innervation through poorly defined mechanisms. We characterized the neonatal beige adipocytes in the inguinal white adipose tissue (iWAT) of C57BL6 postnatal day 3 and 20 mice (P3 and P20) by imaging, genome-wide RNA-seq analysis, ChIP-seq analysis, qRT-PCR validation, and biochemical assays. We found an increase in acetylated histone 3 lysine 27 (H3K27ac) on the promoter and enhancer regions of beige-specific gene UCP1 in iWAT of P20 mice. Furthermore, H3K27ac ChIP-seq analysis in the iWAT of P3 and P20 mice revealed strong H3K27ac signals at beige adipocyte-associated genes in the iWAT of P20 mice. The integration of H3K27ac ChIP-seq and RNA-seq analysis in the iWAT of P20 mice reveal epigenetically active signatures of beige adipocytes, including oxidative phosphorylation and mitochondrial metabolism. We identify the enrichment of GA-binding protein alpha (GABPα) binding regions in the epigenetically active chromatin regions of the P20 iWAT, particularly on beige genes, and demonstrate that GABPα is required for beige adipocyte differentiation. Moreover, transcriptomic analysis and glucose oxidation assays revealed increased glycolytic activity in the neonatal iWAT from P20. Our findings demonstrate that epigenetic mechanisms regulate the development of peri-weaning beige adipocytes via GABPα. Further studies to better understand the upstream mechanisms that regulate epigenetic activation of GABPα and characterization of the metabolic identity of neonatal beige adipocytes will help us harness their therapeutic potential in metabolic diseases.

Identification of a key signaling network regulating perennating bud dormancy in Panax ginseng

BackgroundThe cycle of seasonal dormancy of perennating buds is an essential adaptation of perennial plants to unfavorable winter conditions. Plant hormones are key regulators of this critical biological process, which is intricately connected with diverse internal and external factors. Recently, global warming has increased the frequency of aberrant temperature events that negatively affect the dormancy cycle of perennials. Although many studies have been conducted on the perennating organs of Panax ginseng, the molecular aspects of bud dormancy in this species remain largely unknown. MethodsIn this study, the molecular physiological responses of three P. ginseng cultivars with different dormancy break phenotypes in the spring were dissected using comparative genome-wide RNA-seq and network analyses. These analyses identified a key role for abscisic acid (ABA) activity in the regulation of bud dormancy. Gene set enrichment analysis revealed that a transcriptional network comprising stress-related hormone responses made a major contribution to the maintenance of dormancy. ResultsIncreased expression levels of cold response and photosynthesis-related genes were associated with the transition from dormancy to active growth in perennating buds. Finally, the expression patterns of genes encoding ABA transporters, receptors (PYRs/PYLs), PROTEIN PHOSPHATASE 2Cs (PP2Cs), and DELLAs were highly correlated with different dormancy states in three P. ginseng cultivars. ConclusionThis study provides evidence that ABA and stress signaling outputs are intricately connected with a key signaling network to regulate bud dormancy under seasonal conditions in the perennial plant P. ginseng.

Genome-wide DNA methylation dynamics and RNA-seq analysis during grape (cv. ‘Cabernet Franc’) skin coloration

This study generated whole genome DNA methylation maps to characterize DNA methylomes of grape (cv. ‘Cabernet Franc’) skins and examine their functional significance during grape skin coloration. We sampled grape skin tissues at three key stages (the early stage of grape berry swelling, the late stage of grape berry swelling and the veraison) during which the color of grape berries changed from green to red. DNA methylation levels of grape skins at the three stages were higher in transposable element regions than in the genic regions, and the CG and CHG DNA methylation levels of the genic region were higher than the CHH DNA methylation levels. We identified differentially methylated regions (DMRs) in S2_vs_S1 and S3_vs_S1. The results indicated that DMRs predominantly occurred within the CHH context during grape skin coloration. Many gene ontology (GO)-enriched DMR-related genes were involved in “nucleotide binding,” “catalytic activity” and “ribonucleotide binding” terms; however, many KEGG-enriched DMR-related genes were involved in the “flavonoid biosynthesis” pathway. Our results could provide an important foundation for future research on the development mechanism of grape berries.

LncRNA SOX9-AS1 triggers a transcriptional program involved in lipid metabolic reprogramming, cell migration and invasion in triple-negative breast cancer

At the molecular level, triple-negative breast cancer (TNBC) is frequently categorized as PAM50 basal-like subtype, but despite the advances in molecular analyses, the clinical outcome for these subtypes is uncertain. Long non-coding RNAs (lncRNAs) are master regulators of genes involved in hallmarks of cancer, which makes them suitable biomarkers for breast cancer (BRCA) diagnosis and prognosis. Here, we evaluated the regulatory role of lncRNA SOX9-AS1 in these subtypes. Using the BRCA-TCGA cohort, we observed that SOX9-AS1 was significantly overexpressed in basal-like and TNBC in comparison with other BRCA subtypes. Survival analyzes showed that SOX9-AS1 overexpression was associated with a favorable prognosis in TNBC and basal-like patients. To study the functions of SOX9-AS1, we determined the expression levels in a panel of nine BRCA cell lines finding increased levels in MDA-MB-468 and HCC1187 TNBC. Using subcellular fractionation in these cell lines, we ascertained that SOX9-AS1 was located in the cytoplasmic compartment. In addition, we performed SOX9-AS1 gene silencing using two short-harping constructs, which were transfected in both cell models and performed a genome-wide RNA-seq analysis. Data showed that 351 lncRNAs and 740 mRNAs were differentially expressed in MDA-MB-468 while 56 lncRNAs and 100 mRNAs were modulated in HCC1187 cells (Log2FC < - 1.5 and > 1.5, p.adj value < 0.05). Pathway analysis revealed that the protein-encoding genes potentially regulate lipid metabolic reprogramming, and epithelial–mesenchymal transition (EMT). Expression of lipid metabolic-related genes LIPE, REEP6, GABRE, FBP1, SCD1, UGT2B11, APOC1 was confirmed by RT-qPCR. Functional analysis demonstrated that the knockdown of SOX9-AS1 increases the triglyceride synthesis, cell migration and invasion in both two TNBC cell lines. In conclusion, high SOX9-AS1 expression predicts an improved clinical course in patients, while the loss of SOX9-AS1 expression enhances the aggressiveness of TNBC cells.

Targeted modulation of intestinal epithelial regeneration and immune response in ulcerative colitis using dual-targeting bilirubin nanoparticles.

Rationale: The therapeutic benefits of bilirubin in the treatment of ulcerative colitis (UC) are considerable, whereas the underlying mechanism of bilirubin on UC remains unclear remains unexplored. In addition, the weak hydrophilicity and toxicity have limited its translational applications. Methods: We have developed a colon dual-targeting nanoparticle, for orally delivering bilirubin through hydrogel encapsulation of hyaluronic acid (HA)-modified poly (lactic-co-glycolic acid) (PLGA) nanoparticles (HA-PLGABilirubin). Confocal microscopy and in vivo imaging were used to evaluate the uptake and the targeted property of HA-PLGABilirubin in UC. Immunohistochemistry, immunofluorescence, and transcriptomic analyses were applied to examine the therapeutic effect and potential mechanism of HA-PLGABilirubin in UC. Results: Our results indicated that HA-PLGAbilirubin can significantly enhance the release of bilirubin at simulated intestinal pH and demonstrate higher cellular uptake in inflammatory macrophages. Moreover, in vivo biodistribution studies revealed high uptake and retention of HA-PLGAbilirubin in inflamed colon tissue of UC mouse model, resulting in effective recovery of intestinal morphology and barrier function. Importantly, HA-PLGAbilirubin exerted potent therapeutic efficacy against ulcerative colitis through modulating the intestinal epithelial/stem cells regeneration, and the improvement of angiogenesis and inflammation. Furthermore, genome-wide RNA-seq analysis revealed transcriptional reprogramming of immune response genes in colon tissue upon HA-PLGAbilirubin treatment in UC mouse model. Conclusion: Overall, our work provides an efficient colon targeted drug delivery system to potentiate the treatment of ulcerative colitis via modulating intestinal epithelium regeneration and immune response in ulcerative colitis.

Novel hippocampal genes involved in enhanced susceptibility to chronic pain-induced behavioral emotionality

Altered mood and psychiatric disorders are commonly associated with chronic pain conditions; however, brain mechanisms linking pain and comorbid clinical depression are still largely unknown. In this study, we aimed to identify whether key genes/cellular mechanisms underlie susceptibility/resiliency to development of depressive-like behaviors during chronic pain state. Genome-wide RNA-seq analysis was used to examine the transcriptomic profile of the hippocampus, a limbic brain region that regulates mood and stress responses, from male rats exposed to chronic inflammatory pain. Pain-exposed animals were separated into either ‘resilient’ or ‘susceptible’ to development of enhanced behavioral emotionality based on behavioral testing. RNA-seq bioinformatic analysis, followed by validation using qPCR, revealed dysregulation of hippocampal genes involved in neuroinflammation, cell cycle/neurogenesis and blood-brain barrier integrity. Specifically, ADAM Metallopeptidase Domain 8 (Adam8) and Aurora Kinase B (Aurkb), genes with functional roles in activation of the NLRP3 inflammasome and microgliosis, respectively, were significantly upregulated in the hippocampus of ‘susceptible’ animals expressing increased behavioral emotionality. In addition, genes associated with blood-brain barrier integrity, such as the Claudin 4 (Cldn4), a tight junction protein and a known marker of astrocyte activation, were also significantly dysregulated between ‘resilient’ or ‘susceptible’ pain groups. Furthermore, differentially expressed genes (DEGs) were further characterized in rodents stress models to determine whether their hippocampal dysregulation is driven by common stress responses vs. affective pain processing. Altogether these results continue to strengthen the connection between dysregulation of hippocampal genes involved in neuroinflammatory and neurodegenerative processes with increased behavioral emotionality often expressed in chronic pain state

Sphingolipids in Childhood Asthma and Obesity (SOAP Study): A Protocol of a Cross-Sectional Study.

Asthma and obesity are two of the most common chronic conditions in children and adolescents. There is increasing evidence that sphingolipid metabolism is altered in childhood asthma and is linked to airway hyperreactivity. Dysregulated sphingolipid metabolism is also reported in obesity. However, the functional link between sphingolipid metabolism, asthma, and obesity is not completely understood. This paper describes the protocol of an ongoing study on sphingolipids that aims to examine the pathophysiology of sphingolipids in childhood asthma and obesity. In addition, this study aims to explore the novel biomarkers through a comprehensive multi-omics approach including genomics, genome-wide DNA methylation, RNA-Seq, microRNA (miRNA) profiling, lipidomics, metabolomics, and cytokine profiling. This is a cross-sectional study aiming to recruit 440 children from different groups: children with asthma and normal weight (n = 100), asthma with overweight or obesity (n = 100), overweight or obesity (n = 100), normal weight (n = 70), and siblings of asthmatic children with normal weight, overweight, or obesity (n = 70). These participants will be recruited from the pediatric pulmonology, pediatric endocrinology, and general pediatric outpatient clinics at Sidra Medicine, Doha, Qatar. Information will be obtained from self-reported questionnaires on asthma, quality of life, food frequency (FFQ), and a 3-day food diary that are completed by the children and their parents. Clinical measurements will include anthropometry, blood pressure, biochemistry, bioelectrical impedance, and pulmonary function tests. Blood samples will be obtained for sphingolipid analysis, serine palmitoyltransferase (SPT) assay, whole-genome sequencing (WGS), genome-wide DNA methylation study, RNA-Seq, miRNA profiling, metabolomics, lipidomics, and cytokine analysis. Group comparisons of continuous outcome variables will be carried out by a one-way analysis of variance or the Kruskal-Wallis test using an appropriate pairwise multiple comparison test. The chi-squared test or a Fisher's exact test will be used to test the associations between categorical variables. Finally, multivariate analysis will be carried out to integrate the clinical data with multi-omics data. This study will help us to understand the role of dysregulated sphingolipid metabolism in obesity and asthma. In addition, the multi-omics data from the study will help to identify novel genetic and epigenetic signatures, inflammatory markers, and mechanistic pathways that link asthma and obesity in children. Furthermore, the integration of clinical and multi-omics data will help us to uncover the potential interactions between these diseases and to offer a new paradigm for the treatment of pediatric obesity-associated asthma.

Genome-wide association and RNA-seq analyses reveal a potential gene related to linolenic acid in soybean seeds

Linolenic acid (LA) has poor oxidative stability since it is a polyunsaturated fatty acid. Soybean oil has a high LA content and thus has poor oxidative stability. To identify candidate genes that affect the linolenic acid (LA) content in soybean seeds, a genome-wide association study (GWAS) was performed with 1,060 soybean cultivars collected in China between 2019–2021 and which LA content was measured using matrix-assisted laser desorption/ionization time-of-flight imaging mass spectrometry (MALDI-TOF IMS). A candidate gene, GmWRI14, encoding an APETALA2 (AP2)-type transcription factor, was detected by GWAS in cultivars from all three study years. Multiple sequence alignments showed that GmWRI14 belongs to the plant WRI1 family. The fatty acid contents of different soybean lines were evaluated in transgenic lines with a copy of GmWRI14, control lines without GmWRI14, and the gmwri14 mutant. MALDI-TOF IMS revealed that GmWRI14 transgenic soybeans had a lower LA content with a significant effect on seed size and shape, whereas gmwri14 mutants had a higher LA content. compared to control. The RNA-seq results showed that GmWRI14 suppresses GmFAD3s (GmFAD3B and GmFAD3C) and GmbZIP54 expression in soybean seeds, leading to decreased LA content. Based on the RNA-seq data, yeast one-hybrid (Y1H) and qRT-PCR were performed to confirm the transcriptional regulation of FAD3s by GmWRI14. Our results suggest that FAD3 is indirectly regulated by GmWRI14, representing a new molecular mechanism of fatty acid biosynthesis, in which GmWRI14 regulates LA content in soybean seeds.

Comprehensive Deciphering the Complexity of the Deep Bite: Insight from Animal Model to Human Subjects.

Deep bite is a malocclusion phenotype, defined as the misalignment in the vertical dimension of teeth and jaws and characterized by excessive overlap of the upper front teeth over the lower front teeth. Numerous factors, including genetics, environmental factors, and behavioral ones, might contribute to deep bite. In this study, we discuss the current clinical treatment strategies for deep bite, summarize the already published findings of genetic analysis associated with this complex phenotype, and their constraints. Finally, we propose a comprehensive roadmap to facilitate investigations for determining the genetic bases of this complex phenotype development. Initially, human deep bite phenotype, genetics of human deep bite, the prevalence of human deep bite, diagnosis, and treatment of human deep bite were the search terms for published publications. Here, we discuss these findings and their limitations and our view on future strategies for studying the genetic bases of this complex phenotype. New preventative and treatment methods for this widespread dental issue can be developed with the help of an understanding of the genetic and epigenetic variables that influence malocclusion. Additionally, malocclusion treatment may benefit from technological developments like 3D printing and computer-aided design and manufacture (CAD/CAM). These technologies enable the development of personalized surgical and orthodontic guidelines, enhancing the accuracy and effectiveness of treatment. Overall, the most significant results for the patient can only be achieved with a customized treatment plan created by an experienced orthodontic professional. To design a plan that meets the patient's specific requirements and expectations, open communication between the patient and the orthodontist is essential. Here, we propose to conduct a genome-wide association study (GWAS), RNAseq analysis, integrating GWAS and expression quantitative trait loci (eQTL), micro and small RNA, and long noncoding RNA analysis in tissues associated with deep bite malocclusion in human, and complement it by the same approaches in the collaborative cross (CC) mouse model which offer a novel platform for identifying genetic factors as a cause of deep bite in mice, and subsequently can then be translated to humans. An additional direct outcome of this study is discovering novel genetic elements to advance our knowledge of how this malocclusion phenotype develops and open the venue for early identification of patients carrying the susceptible genetic factors so that we can offer early prevention and treatment strategies, a step towards applying a personalized medicine approach.

Deciphering the genetic architecture of resistance to Corynespora cassiicola in soybean (Glycine max L.) by integrating genome-wide association mapping and RNA-Seq analysis.

Target spot caused by Corynespora cassiicola is a problematic disease in tropical and subtropical soybean (Glycine max) growing regions. Although resistant soybean genotypes have been identified, the genetic mechanisms underlying target spot resistance has not yet been studied. To address this knowledge gap, this is the first genome-wide association study (GWAS) conducted using the SoySNP50K array on a panel of 246 soybean accessions, aiming to unravel the genetic architecture of resistance. The results revealed significant associations of 14 and 33 loci with resistance to LIM01 and SSTA C. cassiicola isolates, respectively, with six loci demonstrating consistent associations across both isolates. To identify potential candidate genes within GWAS-identified loci, dynamic transcriptome profiling was conducted through RNA-Seq analysis. The analysis involved comparing gene expression patterns between resistant and susceptible genotypes, utilizing leaf tissue collected at different time points after inoculation. Integrating results of GWAS and RNA-Seq analyses identified 238 differentially expressed genes within a 200 kb region encompassing significant quantitative trait loci (QTLs) for disease severity ratings. These genes were involved in defense response to pathogen, innate immune response, chitinase activity, histone H3-K9 methylation, salicylic acid mediated signaling pathway, kinase activity, and biosynthesis of flavonoid, jasmonic acid, phenylpropanoid, and wax. In addition, when combining results from this study with previous GWAS research, 11 colocalized regions associated with disease resistance were identified for biotic and abiotic stress. This finding provides valuable insight into the genetic resources that can be harnessed for future breeding programs aiming to enhance soybean resistance against target spot and other diseases simultaneously.

A Parallel Computing Approach to Gene Expression and Phenotype Correlation for Identifying Retinitis Pigmentosa Modifiers in Drosophila

As a genetic eye disorder, retinitis pigmentosa (RP) has been a focus of researchers to find a diagnosis through either genome-wide association (GWA) or RNAseq analysis. In fact, GWA and RNAseq are considered two complementary approaches to gaining a more comprehensive understanding of the genetics of different diseases. However, RNAseq analysis can provide information about the specific mechanisms underlying the disease and the potential targets for therapy. This research proposes a new approach to differential gene expression (DGE) analysis, which is the heart of the core-analysis phase in any RNAseq study. Based on the Drosophila Genetic Reference Panel (DGRP), the gene expression dataset is computationally analyzed in light of eye-size phenotypes. We utilized the foreach and the doParallel R packages to run the code on a multicore machine to reduce the running time of the original algorithm, which exhibited an exponential time complexity. Experimental results showed an outstanding performance, reducing the running time by 95% while using 32 processes. In addition, more candidate modifier genes for RP were identified by increasing the scope of the analysis and considering more datasets that represent different phenotype models.

OA07 Clinical and pathogenic significance of S100A4 overexpression in systemic sclerosis

Abstract Background/Aims Preclinical studies suggest that the S100A4, a Damage Associated Molecular Pattern (DAMP) protein upregulated upon stress or injury, may be a target for antifibrotic therapy. We have studied expression and function of S100A4 in systemic sclerosis (SSc) to further explore its relevance as a driver of fibroblast activation. Methods S100A4 protein concentration was measured by ELISA (Cusabio, USA) in serum of SSc (n = 94) and healthy controls (n = 15). Association of serum level with major organ complications and skin score (mRSS) was determined. S100A4 protein expression in supernatant and cell layer of skin fibroblast cultures from diffuse cutaneous SSc (SScF, n = 6) and healthy controls (NF, n = 6) was assessed. Recombinant S100A4 and a high affinity anti-S100A4 neutralising monoclonal antibody (AX-202) were used to investigate activation of profibrotic molecular and functional phenotype in SSc and control fibroblasts. Results Median [range] S100A4 (ng/ml) was higher in serum of SSc (89.9 [15.0-240.0]) than healthy controls (71.4 [7.9-131.8]; p = 0.027). There was significant association with SSc-ILD (p = 0.025, n = 55), scleroderma renal crisis (p = 0.026, n = 4) and trend of increase in pulmonary hypertension (p = 0.075, n = 7). There was a trend for S100A4 to fall in longer disease duration (r=-0.18, p = 0.09) and a weak (r = 0.28, p = 0.07) correlation with current mRSS but not peak mRSS in dcSSc. Median [range] S100A4 (ng/ml) was higher in culture supernatants of SScF (4.19 [0.52-8.42]) than NF controls (0.28 [0.02-3.29]; p &amp;lt; 0.0001) and SSc cell layer had increased levels by Western blot analysis (mean [sd] RDU NF = 0.50 [0.91]; SScF 8.16 [2.49]; p &amp;lt; 0.0001). S100A4 induced a fibrotic phenotype in NF with concentration dependent increased in Col1, SMA and CTGF expression, and promotion of migration (scratch wound assay) and 3D type I collagen gel contraction that was attenuated by AX-202. AX-202 also reduced the constitutive profibrotic gene and protein expression phenotype of SScF. The effect of S100A4 was significantly attenuated by SB431524, a specific inhibitor of ALK5 (p &amp;lt; 0.0001 ANOVA) using qPCR assays across treatment groups for Col1, SMA and CTGF. Gene expression data (qPCR) confirmed protein expression. Genome-wide RNAseq analysis identified an S100A4 activated signature in NF overlapping the hallmark gene expression signature of SScF. Thus, 464 differentially expressed genes (FDR&amp;lt;0.001 and FC &amp;gt; 1.5) induced in NF by S100A4 were also constitutively overexpressed in SScF, and significantly downregulated by AX-202. Pathway mapping of these S100A4 dependent transcripts in SSc showed the most significant enriched Kegg® pathways (FDR&amp;lt;0.001) were regulation of stem cell pluripotency (4.6-fold) and metabolic pathways (1.9-fold). Conclusion Our findings provide compelling evidence supporting a profibrotic role for S100A4 in fibroblast activation in SSc and suggest that serum level may be a biomarker of major organ manifestations and disease severity. Studies to elucidate the therapeutic potential of targeting S100A4 in SSc are warranted. Disclosure C.P. Denton: None. S. Xu: None. R.H. Maclean: None. K.E. Clark: None. S. Borchert: Corporate appointments; Employee of ARXX Therapeutics. R. Hussain: Corporate appointments; Employee of ARXX Therapeutics. J. Klingelhöfer: Corporate appointments; Employee of ARXX Therapeutics. J. Hallén: Corporate appointments; Employee of ARXX Therapeutics. V.H. Ong: None.

Genome-wide identification and characterization of long non-coding RNA in barley roots in response to Piriformospora indica colonization

Currently, there is very limited information about long noncoding RNAs (lncRNAs) found in barley. It remains unclear whether barley lncRNAs are responsive to Piriformospora indica (P. indica) colonization.We found that barley roots exhibited fast development and that large roots branched after P. indica colonization. Genome-wide high-throughput RNA-seq and bioinformatic analysis showed that 4356 and 5154 differentially expressed LncRNAs (DELs) were found in response to P. indica at 3 and 7 days after colonization (dai), respectively, and 2456 DELs were found at 7 dai compared to 3 dai. Based on the coexpression correlation of lncRNAmRNA, we found that 98.6% of lncRNAs were positively correlated with 3430 mRNAs at 3 dai and 7 dai. Further GO analysis showed that 30 lncRNAs might be involved in the regulation of gene transcription; 23 lncRNAs might participate in cell cycle regulation. Moreover, the metabolite analysis indicated that chlorophyll a, sucrose, protein, gibberellin, and auxin were in accordance with the results of the transcriptome, and the respective lncRNAs were positively correlated with these target RNAs. Gene silencing suggested that lncRNA TCONS_00262342 is probably a key regulator of GA3 synthesis pathway, which participates in P. indica and barley interactions. We concluded that acting as a molecular material basis and resource, lncRNAs respond to P. indica colonization by regulating metabolite content in barley and coordinate the complex regulatory process of higher life by constructing highly positive correlations with their target mRNAs.

Integrating genome-wide association study with RNA-seq revealed DBI as a good candidate gene for intramuscular fat content in Beijing black pigs.

Increasing intramuscular fat (IMF) content can enhance the sensory quality of meat, including tenderness, juiciness, flavor, and color. Genome-wide association study and RNA-sequencing (RNA-seq) analysis were used to identify candidate IMF genes in Beijing Black pigs, a popular species among consumers in northern China. Two and three single nucleotide polymorphisms were significantly associated with IMF in SSC13 and SSC15 respectively. Solute carrier family 4 member 7 (SLC4A7) on SSC13 and insulin induced gene 2 (INSIG2), coiled-coil domain containing 93 (CCDC93), and diazepam binding inhibitor acyl-CoA binding protein (DBI) on SSC15 are good candidate genes in this population. Furthermore, RNA-seq analysis was performed between high and low IMF groups, and 534 differentially expressed genes were identified. In addition, based on differentially expressed genes, Kyoto Encyclopedia of Genes and Genomes analysis revealed that peroxisome proliferator-activated receptors and FoxO signaling pathway pathways might contribute to IMF. Moreover, the DBI gene was identified as a candidate for IMF both by genome-wide association study and RNA-seq analysis, suggesting that it might be a crucial candidate gene for influencing IMF in Beijing Black pigs.

Genes and pathways associated with fear discrimination identified by genome-wide DNA methylation and RNA-seq analyses in nucleus accumbens in mice

Fear memory is critical for individual survival. However, the maladaptive fear response is one of the hallmarks of fear-related disorders, which is characterized by the failure to discriminate threatening signals from neutral or safe cues. The biological mechanisms of fear discrimination remain to be clarified. In this study, we found that the nucleus accumbens (NAc) was indispensable for the formation of cued fear memory in mice, during which the expression of DNA methyltransferase 3a gene (DNMT3a) increased. Injection of Zebularine, a nonspecific DNMT inhibitor, into NAc immediately after conditioning induced a maladaptive fear response to neutral cue (CS−). Using whole-genome bisulfite sequencing (WGBS), differentially methylated sites and methylated regions (DMRs) were investigated. 16,226 DMRs in the genenome were identified, in which, 214 genes with significant differences in their methylation levels and mRNA expression profiles were identified through correlation analysis. Notably, 15 genes were synaptic function-related and 8 genes were enriched in the cGMP-PKG signaling pathway. Moreover, inhibition of PKG impaired fear discrimination. Together, our results revealed the profile and role of genome-wide DNA methylation in NAc in the regulation of fear discrimination.

Modeling HPV-Associated Disease and Cancer Using the Cottontail Rabbit Papillomavirus.

Approximately 5% of all human cancers are attributable to human papillomavirus (HPV) infections. HPV-associated diseases and cancers remain a substantial public health and economic burden worldwide despite the availability of prophylactic HPV vaccines. Current diagnosis and treatments for HPV-associated diseases and cancers are predominantly based on cell/tissue morphological examination and/or testing for the presence of high-risk HPV types. There is a lack of robust targets/markers to improve the accuracy of diagnosis and treatments. Several naturally occurring animal papillomavirus models have been established as surrogates to study HPV pathogenesis. Among them, the Cottontail rabbit papillomavirus (CRPV) model has become known as the gold standard. This model has played a pivotal role in the successful development of vaccines now available to prevent HPV infections. Over the past eighty years, the CRPV model has been widely applied to study HPV carcinogenesis. Taking advantage of a large panel of functional mutant CRPV genomes with distinct, reproducible, and predictable phenotypes, we have gained a deeper understanding of viral–host interaction during tumor progression. In recent years, the application of genome-wide RNA-seq analysis to the CRPV model has allowed us to learn and validate changes that parallel those reported in HPV-associated cancers. In addition, we have established a selection of gene-modified rabbit lines to facilitate mechanistic studies and the development of novel therapeutic strategies. In the current review, we summarize some significant findings that have advanced our understanding of HPV pathogenesis and highlight the implication of the development of novel gene-modified rabbits to future mechanistic studies.

Nitrate enhances the secondary growth of storage roots in Panax ginseng

BackgroundNitrogen (N) is an essential macronutrient for plant growth and development. To support agricultural production and enhance crop yield, two major N sources, nitrate and ammonium, are applied as fertilizers to the soil. Although many studies have been conducted on N uptake and signal transduction, the molecular genetic mechanisms of N-mediated physiological roles, such as the secondary growth of storage roots, remain largely unknown. MethodsOne-year-old P. ginseng seedlings treated with KNO3 were analyzed for the secondary growth of storage roots. The histological paraffin sections were subjected to bright and polarized light microscopic analysis. Genome-wide RNA-seq and network analysis were carried out to dissect the molecular mechanism of nitrate-mediated promotion of ginseng storage root thickening. ResultsHere, we report the positive effects of nitrate on storage root secondary growth in Panax ginseng. Exogenous nitrate supply to ginseng seedlings significantly increased the root secondary growth. Histological analysis indicated that the enhancement of root secondary growth could be attributed to the increase in cambium stem cell activity and the subsequent differentiation of cambium-derived storage parenchymal cells. RNA-seq and gene set enrichment analysis (GSEA) revealed that the formation of a transcriptional network comprising auxin, brassinosteroid (BR)-, ethylene-, and jasmonic acid (JA)-related genes mainly contributed to the secondary growth of ginseng storage roots. In addition, increased proliferation of cambium stem cells by a N-rich source inhibited the accumulation of starch granules in storage parenchymal cells. ConclusionThus, through the integration of bioinformatic and histological tissue analyses, we demonstrate that nitrate assimilation and signaling pathways are integrated into key biological processes that promote the secondary growth of P. ginseng storage roots.

Early and late genome-wide gastric epithelial transcriptome response during infection with the human carcinogen Helicobacterpylori.

Infection of the stomach by Helicobacter pylori is a major risk factor for the development of gastric cancer. Colonization of the gastric epithelium leads to the activation of multiple disease-related signaling pathways. Serine protease HtrA represents an important secreted virulence factor that mediates cleavage of cellular junctions. However, its potential role in nuclear responses is unknown. Here, we performed a genome-wide RNA-seq analysis of polarized gastric epithelial cells infected by wild-type (wt) and ΔhtrA mutant bacteria. Fluorescence microscopy showed that H.pylori wt, but not ΔhtrA bacteria, preferably localized at cellular junctions. Our results pinpointed early (2h) and late (6h) transcriptional responses, with most differentially expressed genes at 6h post infection. The transcriptomes revealed HtrA-dependent targeting of genes associated with inflammation and apoptosis (e.g. IL8, ZFP36, TNF). Accordingly, infection with the ΔhtrA mutant induced increased apoptosis rates in host cells, which was associated with reduced H.pylori CagA expression. In contrast, transcription of various carcinogenesis-associated genes (e.g. DKK1, DOCK8) was affected by H.pylori independent of HtrA. These findings suggest that H.pylori disturbs previously unknown molecular pathways in an HtrA-dependent and HtrA-independent manner, and provide valuable new insights of this significant pathogen in humans and thus potential targets for better controlling the risk of malignant transformation.