Whole-genome Gene Expression Research Articles

EQTL mapping in transgenic alpha-synuclein carrying Caenorhabditis elegans recombinant inbred lines.

Protein aggregation of α-synuclein (αS) is a genetic and neuropathological hallmark of Parkinson's disease (PD). Studies in the model nematode Caenorhabditis elegans suggested that variation of αS aggregation depends on the genetic background. However, which genes and genetic modifiers underlie individual differences in αS pathology remains unknown. To study the genotypic-phenotypic relationship of αS aggregation, we constructed a Recombinant Inbred Line (RIL) panel derived from a cross between genetically divergent strains C. elegans NL5901 and SCH4856, both harboring the human αS gene. As a first step to discover genetic modifiers 70 αS-RILs were measured for whole-genome gene expression and expression quantitative locus analysis (eQTL) were mapped. We detected multiple eQTL hot-spots, many of which were located on Chromosome V. To confirm a causal locus, we developed Introgression Lines (ILs) that contain SCH4856 introgressions on Chromosome V in an NL5901 background. We detected 74 genes with an interactive effect between αS and the genetic background, including the human p38 MAPK homologue pmk-1 that has previously been associated with PD. Together, we present a unique αS-RIL panel for defining effects of natural genetic variation on αS pathology, which contributes to finding genetic modifiers of PD.

Melanoma-specific mutation hotspots in distal, non-coding, promoter-interacting regions implicate novel candidate driver genes.

To develop targeted treatments, it is crucial to identify the full spectrum of genetic drivers in melanoma, including those in non-coding regions. However, recent efforts to explore non-coding regions have primarily focused on gene-adjacent elements such as promoters and non-coding RNAs, leaving intergenic distal regulatory elements largely unexplored. We used Hi-C chromatin contact data from melanoma cells to map distal, non-coding, promoter-interacting regulatory elements genome-wide in melanoma. Using this "promoter-interaction network", alongside whole-genome sequence and gene expression data from the Pan Cancer Analysis of Whole Genomes, we developed multivariate linear regression models to identify distal somatic mutation hotspots that affect promoter activity. We identified eight recurrently mutated hotspots that are novel, melanoma-specific, located in promoter-interacting distal regulatory elements, alter transcription factor binding motifs, and affect the expression of genes (e.g., HSPB7, CLDN1, ADCY9 and FDXR) previously implicated as tumour suppressors/oncogenes in various cancers. Our study suggests additional non-coding drivers beyond the well-characterised TERT promoter in melanoma, offering new insights into the disruption of complex regulatory networks by non-coding mutations that may contribute to melanoma development. Furthermore, our study provides a framework for integrating multiple levels of biological data to uncover cancer-specific non-coding drivers.

Exploring the dynamic three-dimensional chromatin architecture and transcriptional landscape in goose liver tissues underlying metabolic adaptations induced by a high-fat diet

BackgroundGoose, descendants of migratory ancestors, have undergone extensive selective breeding, resulting in their remarkable ability to accumulate fat in the liver and exhibit a high tolerance for significant energy intake. As a result, goose offers an excellent model for studying obesity, metabolic disorders, and liver diseases in mammals. Although the impact of the three-dimensional arrangement of chromatin within the cell nucleus on gene expression and transcriptional regulation is widely acknowledged, the precise functions of chromatin architecture reorganization during fat deposition in goose liver tissues still need to be fully comprehended.ResultsIn this study, geese exhibited more pronounced changes in the liver index and triglyceride (TG) content following the consumption of the high-fat diet (HFD) than mice without significant signs of inflammation. Additionally, we performed comprehensive analyses on 10 goose liver tissues (5 HFD, 5 normal), including generating high-resolution maps of chromatin architecture, conducting whole-genome gene expression profiling, and identifying H3K27ac peaks in the livers of geese and mice subjected to the HFD. Our results unveiled a multiscale restructuring of chromatin architecture, encompassing Compartment A/B, topologically associated domains, and interactions between promoters and enhancers. The dynamism of the three-dimensional genome architecture, prompted by the HFD, assumed a pivotal role in the transcriptional regulation of crucial genes. Furthermore, we identified genes that regulate chromatin conformation changes, contributing to the metabolic adaptation process of lipid deposition and hepatic fat changes in geese in response to excessive energy intake. Moreover, we conducted a cross-species analysis comparing geese and mice exposed to the HFD, revealing unique characteristics specific to the goose liver compared to a mouse. These chromatin conformation changes help elucidate the observed characteristics of fat deposition and hepatic fat regulation in geese under conditions of excessive energy intake.ConclusionsWe examined the dynamic modifications in three-dimensional chromatin architecture and gene expression induced by an HFD in goose liver tissues. We conducted a cross-species analysis comparing that of mice. Our results contribute significant insights into the chromatin architecture of goose liver tissues, offering a novel perspective for investigating mammal liver diseases.

Epigenetic control on transcription of vernalization genes and whole-genome gene expression profile induced by vernalization in common wheat

Vernalization is necessary for winter wheat to flower. However, it is unclear whether vernalization is also required for spring wheat, which is frequently sown in fall, and what molecular mechanisms underlie the vernalization response in wheat varieties. In this study, we examined the molecular mechanisms that regulate vernalization response in winter and spring wheat varieties. For this purpose, we determined how major vernalization genes (VRN1, VRN2, and VRN3) respond to vernalization in these varieties and whether modifications to histones play a role in changes in gene expression. We also identified genes that are differentially regulated in response to vernalization in winter and spring wheat varieties. We found that in winter wheat, but not in spring wheat, VRN1 expression decreases when returned to warm temperature following vernalization. This finding may be associated with differences between spring and winter wheat in the levels of tri-methylation of lysine 27 on histone H3 (H3K27me3) and tri-methylation of lysine 4 on histone H3 (H3K4me3) at the VRN1 gene. Analysis of winter wheat transcriptomes before and after vernalization revealed that vernalization influences the expression of several genes, including those involved in leucine catabolism, cysteine biosynthesis, and flavonoid biosynthesis. These findings provide new candidates for further study on the mechanism of vernalization regulation in wheat.

Preliminary study on whole genome methylation and transcriptomics in age-related cataracts

DNA methylation plays an important role in the occurrence and development of age-related cataracts (ARC). This study aims to reveal potential epigenetic biomarkers of ARC by detecting modifications to the DNA methylation patterns of genes shown to be related to ARC by transcriptomics. The MethylationEPIC BeadChip (850 K) was used to analyze the DNA methylation levels in ARC patients and unaffected controls, and the Pearson correlation test was used to perform genome-wide integration analysis of DNA methylation and transcriptome data. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were used to perform functional analysis of the whole genome, promoter regions (TSS1500/TSS200), and the associated differentially methylated genes (DMG). Pyrosequencing was used to verify the methylation levels of the selected genes. The results showed that, compared with the control group, a total of 52,705 differentially methylated sites were detected in the ARC group, of which 13,858 were hypermethylated and 38,847 were hypomethylated. GO and KEGG analyses identified functions related to the cell membrane, the calcium signaling pathway, and their possible molecular mechanisms. Then, 57 DMGs with negative promoter methylation correlations were screened by association analysis. Pyrosequencing verified that the ARC group had higher methylation levels of C3 and CCKAR and lower methylation levels of NLRP3, LEFTY1, and GPR35 compared with the control group. In summary, our study reveals the whole-genome DNA methylation patterns and gene expression profiles in ARC, and the molecular markers of methylation identified herein may aid in the prevention, diagnosis, treatment, and prognosis of ARC.

Music compensates for altered gene expression in age-related cognitive disorders

Extensive literature has explored the beneficial effects of music in age-related cognitive disorders (ACD), but limited knowledge exists regarding its impact on gene expression. We analyzed transcriptomes of ACD patients and healthy controls, pre-post a music session (n = 60), and main genes/pathways were compared to those dysregulated in mild cognitive impairment (MCI) and Alzheimer’s disease (AD) as revealed by a multi-cohort study (n = 1269 MCI/AD and controls). Music was associated with 2.3 times more whole-genome gene expression, particularly on neurodegeneration-related genes, in ACD than in controls. Co-expressed gene-modules and pathways analysis demonstrated that music impacted autophagy, vesicle and endosome organization, biological processes commonly dysregulated in MCI/AD. Notably, the data indicated a strong negative correlation between musically-modified genes/pathways in ACD and those dysregulated in MCI/AD. These findings highlight the compensatory effect of music on genes/biological processes affected in MCI/AD, providing insights into the molecular mechanisms underlying the benefits of music on these disorders.

Mechanism of Androgen-Independent Stromal Proliferation in Benign Prostatic Hyperplasia.

Benign prostatic hyperplasia (BPH) is a chronic proliferative disease showing stromal-dominant proliferation. However, the detailed proliferation mechanism has remained unclear. Although aging and androgen have been reported as definitive risk factors for BPH, recent studies have focused on the involvement of androgen-independent factors. Androgen-independent factors include ischemia, oxidative stress, metabolic syndrome, infection, autoimmune reactions, and inflammation, with inflammation in BPH tissues playing a central role in the BPH proliferative process. Inflammation in BPH tissues by various factors finally leads to tissue remodeling and stromal proliferation through the wound healing process of the prostate. To elucidate the proliferative mechanism of BPH, a study using whole-genome gene expression analysis in a stromal-dominant BPH rat model was performed and showed that immune response-related pathways and complement classical pathways are activated. Furthermore, expression analysis using this BPH rat model showed that the autoimmune reaction triggered complement pathway activation in the proliferative process of BPH. BPH is a multifactorial disease, and understanding the role of androgen-independent factors including immune responses contributes to elucidating the pathogenesis of BPH. Androgen-independent factors may lead to new therapeutic targets for BPH, and further development of this research is expected.

Whole-genome DNA methylation and gene expression profiling in the livers of mice with nonalcoholic steatohepatitis

AimsNon-alcoholic fatty liver disease (NAFLD) has emerged as one of the major causes of liver-related morbidity and mortality. It ranges simple steatosis to non-alcoholic steatohepatitis (NASH). Previous studies have shown that epigenetic factors, such as DNA methylation, can contribute to the development and progression of simple steatosis. However, the profiling of whole-genome DNA methylation remains poorly characterized in NASH. Main methodsIn this study, we established a mouse model of diet-induced NASH, by maintaining male mice on a high-fructose-high-cholesterol diet (HFHC), to generate hepatic steatosis, inflammation and injury. We profiled hepatic gene expression by RNA-Sequencing and locus-specific 5-methylcytosine level, using Whole Genome Bisulfite Sequencing (WGBS). Key findingsWe identified >1000 differentially methylated regions in NASH versus control group, indicating that NASH diet could modulate the liver methylome. Furthermore, integrated analysis of methylome and transcriptome identified certain key methylated genes and pathways, which may be involved in steroid metabolism and inflammation response. The liver methylation levels of key genes especially Tgfb, Msn, Iqgap1, Cyba, Fcgr1 decreased, and their consequent increased expression may lead to NASH development. SignificanceWe found that HFHC diet-induced NASH could induces genome-wide differential DNA methylation changes. Thus, we proposed that DNA methylation profiles of genomes may be a useful signature of gene transcription and may play an important role in the development of NASH. We also screened and validated the changes of key genes, which may provide new perspectives for the mechanistic study of NASH in future.

Abstract 4901: The role of protein disulfide isomerase and copper in the paraptotic cell death of clinically investigated anticancer thiosemicarbazones

Abstract α-N-Heterocyclic thiosemicarbazones (TSCs) have long been investigated as anticancer compounds. Triapine is one of the best-known TSCs for anticancer therapy and is currently tested in a clinical phase III trial. To further improve the anticancer activity of TSCs, novel derivatives (such as DpC and COTI-2) have been developed and clinically investigated for their activity against solid tumors. These novel TSCs belong to a subclass (incl. the tetra-methylated Triapine derivative, Me2NNMe2) that are terminally disubstituted and demonstrated enhanced anticancer activity in vitro and in vivo. The improved activity of these TSCs correlated with their ability to form stable copper complexes as well as induce paraptosis, a form of programmed but caspase-independent cell death, that is characterized by the formation of vesicles originating from the endoplasmic reticulum (ER). However, the molecular events behind paraptosis induction are still not fully elucidated. Consequently, the aim of this study was to investigate the molecular signaling behind paraptosis as well as the reasons behind the induction of this quite unknown form of cell death by TSCs. A whole-genome gene expression analysis of cells treated with a terminally disubstituted TSC revealed an upregulation of thiol homeostasis-related genes, copper metabolism, and metallothionein genes as well as ER stress response genes. In addition, we found an upregulation of the copper-sensitive and thiol-containing protein disulfide isomerase (PDI), which was confirmed on protein level. In cell-free activity studies, PDI was indeed inhibited by terminally disubstituted TSCs, although only in form of their copper complexes. In agreement, removal of copper ions by specific chelators greatly reduced their anticancer activity as well as paraptotic potential in cell culture. This inhibition was hypothesized to be induced by a general disruption of the thiol redox homeostasis affecting also other responsive proteins and resulting in a more oxidative environment, especially in the ER. In confirmation, thiol-containing small molecules could reduce the compound’s anticancer activity in viability assays. Furthermore, glutathione as well as PDI were detected predominantly in their oxidized form in cells treated with terminally disubstituted TSCs. Consequently, terminally disubstituted TSCs, which are characterized by higher anticancer activity, induce paraptosis due to their high copper complex stability. The formation of the complex results in the interaction with thiol-containing proteins and subsequently in the disruption of the thiol-redox homeostasis, which leads to paraptotic cell death. Overall this work shed light on the paraptotic cell death and will be of interest for future clinical development of anticancer TSCs, as paraptosis is hypothesized to overcome apoptosis-resistance of cancer cells. Citation Format: Sonja Hager, Walter Berger, Christian R. Kowol, Eva A. Enyedy, Petra Heffeter. The role of protein disulfide isomerase and copper in the paraptotic cell death of clinically investigated anticancer thiosemicarbazones. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4901.

Tidal emersion effects on universal metrics, elemental contents, and health conditions of Pacific oyster (Crassostrea gigas) on artificial reefs in Yellow River Delta, China

Oysters form natural reefs and are important coastal foundation species that provide a wide range of ecosystem services. Clarification of the mechanisms by which tidal emersion affects oyster reef performance will help guide and direct future oyster reef restoration efforts. In the present study, we constructed artificial Pacific oyster (Crassostrea gigas) reefs in the coastal zones of the Yellow River Delta and quantified the effects of tidal emersion on the growth, density, biomass, element (carbon, nitrogen, phosphorus) content within tissues, and health conditions of oysters. The results showed that temperature change driven by the tidal emersion effect was one of the major forces affecting reef performance. In the experimental field, fully immersed oysters performed better than partially exposed ones in terms of growth and density. The mean biomass of the fully immersed oysters was 2.88-fold higher than that of the partially exposed oysters (p < 0.05). However, we found that tidal emersion had no significant effect on the carbon, nitrogen, and phosphorus content within oyster tissues. Oyster health condition revealed by the biochemical marker and whole-genome gene expression analyses demonstrated that the oysters partially exposed to air presented enhanced oxidative and lipid metabolism and were forced to allocate metabolic energy towards defense against stress rather than growth, which suggested sub-optimal health conditions. Overall, the study indicated that full immersion is optimal for oyster reef restoration in the experimental field. The results of this research will provide future conservationists, aquaculturists, and other specialists with a better understanding of tidal emersion effects on the success of oyster reef restoration to achieve region-specific restoration goals and improve protection of the coastal environment.

Integrated whole-genome gene expression analysis reveals an atlas of dynamic immune landscapes after myocardial infarction.

Myocardial infarction (MI) is a deadly medical condition leading to irreversible damage to the inflicted cardiac tissue. Elevated inflammatory response marks the severity of MI and is associated with the development of heart failure (HF), a long-term adverse outcome of MI. However, the efficacy of anti-inflammatory therapies for MI remains controversial. Deciphering the dynamic transcriptional signatures in peripheral blood mononuclear cells (PBMCs) is a viable and translatable route to better understand post-MI inflammation, which may help guide post-MI anti-inflammatory treatments. In this work, integrated whole-genome gene expression analysis was performed to explore dynamic immune landscapes associated with MI. GSEA and GSVA showed that pathways involved in the inflammatory response and metabolic reprogramming were significantly enriched in PBMCs from MI patients. Based on leukocyte profiles generated by xCell algorithm, the relative abundance of monocytes and neutrophils was significantly increased in PBMCs from MI patients and had positive correlations with typical inflammation-associated transcripts. Mfuzz clustering revealed temporal gene expression profiles of PBMCs during the 6-month post-MI follow-up. Analysis of DEGs and gene sets indicated that PBMCs from HF group were characterized by elevated and lasting expression of genes implicated in inflammation and coagulation. Consensus clustering generated 4 metabolic subtypes of PBMCs with molecular heterogeneity in HF patients. In summary, integrated whole-genome gene expression analysis here outlines a transcriptomic framework that may improve the understanding of dynamic signatures present in PBMCs, as well as the heterogeneity of PBMCs in MI patients with or without long-term clinical outcome of HF. Moreover, the work here uncovers the diversity and heterogeneity of PBMCs from HF patients, providing novel bioinformatic evidence supporting the mechanistic implications of metabolic reprogramming and mitochondrial dysfunction in the post-MI inflammation and HF. Therefore, our work here supports the notion that individualized anti-inflammatory therapies are needed to improve the clinical management of post-MI patients.

Gene co-expression network analysis for the selection of candidate early warning indicators of heat and nutrient stress in Posidonia oceanica

The continuous worldwide seagrasses decline calls for immediate actions in order to preserve this precious marine ecosystem. The main stressors that have been linked with decline in seagrasses are 1) the increasing ocean temperature due to climate change and 2) the continuous inputs of nutrients (eutrophication) associated with coastal human activities. To avoid the loss of seagrass populations, an “early warning” system is needed. We used Weighed Gene Co-expression Network Analysis (WGCNA), a systems biology approach, to identify potential candidate genes that can provide an early warning signal of stress in the Mediterranean iconic seagrass Posidonia oceanica, anticipating plant mortality. Plants were collected from both eutrophic (EU) and oligotrophic (OL) environments and were exposed to thermal and nutrient stress in a dedicated mesocosm. By correlating the whole-genome gene expression after 2-weeks exposure with the shoot survival percentage after 5-weeks exposure to stressors, we were able to identify several transcripts that indicated an early activation of several biological processes (BP) including: protein metabolic process, RNA metabolic process, organonitrogen compound biosynthetic process, catabolic process and response to stimulus, which were shared among OL and EU plants and among leaf and shoot apical meristem (SAM), in response to excessive heat and nutrients. Our results suggest a more dynamic and specific response of the SAM compared to the leaf, especially the SAM from plants coming from a stressful environment appeared more dynamic than the SAM from a pristine environment. A vast list of potential molecular markers is also provided that can be used as targets to assess field samples.

Transcriptomic study of gastrointestinal stromal tumors with liver metastasis.

Introduction: GIST (gastrointestinal stromal tumor) is the most prominent mesenchymal neoplasms of the gastrointestinal tract, and liver is the most common metastasis site for GIST. The molecular mechanism leading to liver metastasis of GIST is currently unclear. Methods: With the goal of revealing the underlying mechanism, we performed whole-genome gene expression profiling on 18 pairs of RNA samples comprised of GIST tissues (with liver metastasis) and corresponding non-tumor tissues. After identifying differentially expressed gene, functional annotation and signal pathway analyses were conducted. GSE13861, datasets that compare GIST (without liver metastasis) with adjacent tissues, served as a comparison. Results: A total of 492 up-regulated genes and 629 down-regulated genes were identified as differentially expressed genes between liver metastasis tissues and non-tumor tissues. We characterized expression patterns of DEGs identified from our cohort and GSE13861 that show signatures of enrichment for functionality. In subsequent gene set enrichment analysis, differentially expressed genes were mainly enriched in Epithelial Mesenchymal Transition in both datasets. 493 genes were overlapped among our whole-genome gene expression profiling results and GSE13861, consisting 188 up-regulated genes and 305 down-regulated genes. By using CytoHubba plugin of Cytoscape, CDH1, CD34, KIT, PROM1, SOX9, FGF2, CD24, ALDH1A1, JAG1 and NES were identified as top ten hub genes in tumorigenesis and liver metastasis of GIST. higher expression levels of FGF2, JAG1, CD34, ALDH1A1 and the lower expression level of CDH1 were respectively associated with unfavorable overall survival. Meanwhile higher expression levels of CD34, FGF2, KIT, JAG1, ALDH1A were correlated with worse disease-free survival. Discussion: The present study may help to provide candidate pathways and targets for treatment of GIST and prevention methods to liver metastasis.

Host Transcriptome Analysis of Ferret Tissues Following Henipavirus Infection.

Ferrets are commonly used as experimental models of infection for a variety of viruses due to their susceptibility to human respiratory viruses and the close resemblance of pathological outcomes found in human infections. Even though ferret-specific reagents are limited, the use of ferrets as a preclinical experimental model of infection has gained considerable interest since the publication of the ferret transcriptome and draft ferret genome. These advances have made it feasible to easily perform whole-genome gene expression analysis in the ferret infection model. Here, we describe methods for genome-wide gene expression analysis using RNA sequence (RNAseq) data obtained from the lung and brain tissues obtained from experimental infections of Hendra (HeV) and Nipah (NiV) viruses in ferrets. We provide detailed methods for RNAseq and representative data for host gene expression profiles of the lung tissues that show early activation of interferon pathways and later activation of inflammation-related pathways.

Integrative Methylome and Transcriptome Characterization Identifies SERINC2 as a Tumor-Driven Gene for Papillary Thyroid Carcinoma

Most papillary thyroid carcinomas (PTCs) can be diagnosed preoperatively by routine evaluation, such as thyroid ultrasonography and fine-needle aspiration biopsy. Nevertheless, understanding how to differentiate indolent thyroid tumors from aggressive thyroid cancers remains a challenge, which may cause overtreatment. This study aimed to identify papillary thyroid cancer-specific indicators with whole-genome DNA methylation and gene expression profiles utilizing Infinium Methylation EPIC BeadChip (850k) and RNA arrays. In this paper, we report SERINC2 as a potential tumor-driven indicator in PTC. The up-regulated expression levels of SERINC2 were verified in PTC cell lines via qPCR. Then, cell counting kit 8 (CCK-8), wound healing, and flow cytometric assays were performed to confirm the influence of SERINC2 on proliferation and apoptosis in PTC cell lines after intervention or overexpression. Moreover, the investigation of data from the Cancer Dependency Map (DepMap) provided a potential pathway targeted by SERINC2. The activation of the tryptophan metabolic pathway may reduce the dependency of SERINC2 in thyroid cancers. In conclusion, our results demonstrate the whole-genome DNA methylation and gene expression profiles of papillary thyroid carcinoma, identify SERINC2 as a potential tumor-driven biomarker, and preliminarily verify its function in PTC.

QTL-mapping in the obese Berlin Fat Mouse identifies additional candidate genes for obesity and fatty liver disease

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with liver weight, liver triglycerides, and body weight using the obese BFMI sub-line BFMI861-S1. BFMI861-S1 mice are insulin resistant and store ectopic fat in the liver. In generation 10, 58 males and 65 females of the advanced intercross line (AIL) BFMI861-S1xB6N were phenotyped under a standard diet over 20 weeks. QTL analysis was performed after genotyping with the MiniMUGA Genotyping Array. Whole-genome sequencing and gene expression data of the parental lines was used for the prioritization of positional candidate genes. Three QTLs associated with liver weight, body weight, and subcutaneous adipose tissue (scAT) weight were identified. A highly significant QTL on chromosome (Chr) 1 (157–168 Mb) showed an association with liver weight. A QTL for body weight at 20 weeks was found on Chr 3 (34.1–40 Mb) overlapping with a QTL for scAT weight. In a multiple QTL mapping approach, an additional QTL affecting body weight at 16 weeks was identified on Chr 6 (9.5–26.1 Mb). Considering sequence variants and expression differences, Sec16b and Astn1 were prioritized as top positional candidate genes for the liver weight QTL on Chr 1; Met and Ica1 for the body weight QTL on Chr 6. Interestingly, all top candidate genes have previously been linked with metabolic traits. This study shows once more the power of an advanced intercross line for fine mapping. QTL mapping combined with a detailed prioritization approach allowed us to identify additional and plausible candidate genes linked to metabolic traits in the BFMI861-S1xB6N AIL. By reidentifying known candidate genes in a different crossing population the causal link with specific traits is underlined and additional evidence is given for further investigations.

Development and validation of an RNA-seq-based transcriptomic risk score for asthma

Recent progress in RNA sequencing (RNA-seq) allows us to explore whole-genome gene expression profiles and to develop predictive model for disease risk. The objective of this study was to develop and validate an RNA-seq-based transcriptomic risk score (RSRS) for disease risk prediction that can simultaneously accommodate demographic information. We analyzed RNA-seq gene expression data from 441 asthmatic and 254 non-asthmatic samples. Logistic least absolute shrinkage and selection operator (Lasso) regression analysis in the training set identified 73 differentially expressed genes (DEG) to form a weighted RSRS that discriminated asthmatics from healthy subjects with area under the curve (AUC) of 0.80 in the testing set after adjustment for age and gender. The 73-gene RSRS was validated in three independent RNA-seq datasets and achieved AUCs of 0.70, 0.77 and 0.60, respectively. To explore their biological and molecular functions in asthma phenotype, we examined the 73 genes by enrichment pathway analysis and found that these genes were significantly (p < 0.0001) enriched for DNA replication, recombination, and repair, cell-to-cell signaling and interaction, and eumelanin biosynthesis and developmental disorder. Further in-silico analyses of the 73 genes using Connectivity map shows that drugs (mepacrine, dactolisib) and genetic perturbagens (PAK1, GSR, RBM15 and TNFRSF12A) were identified and could potentially be repurposed for treating asthma. These findings show the promise for RNA-seq risk scores to stratify and predict disease risk.

Genome-Wide Identification and Expression Analysis of Chitinase-like Genes in Petunia axillaris.

Chitinase (EC 3.2.1.14) is a kind of chitin-degrading glycosidase, which plays important roles in the abiotic and biotic defense of plants. In this study, we conducted whole-genome annotation, molecular evolution, and gene expression analyses on the chitinase-like (CTL) gene family members of Petunia axillaris. Thirty-three Petunia axillaris chitinase-like genes (PaCTLs) were identified from the latest Petunia genome database. According to the phylogenetic analyses, these genes were divided into GH18 and GH19 subgroups and further subdivided into five classes (Class I to Class V). Conserved motif arrangements indicated their functional relevance within each group. The expansion and homeology analyses showed that gene replication events played an important role in the evolution of PaCTLs and the increase of the GH18 subgroup members was the main reason for the expansion of the PaCTL gene family in the evolution progress. By qRT-PCR analysis, we found that most of the PaCTLs showed a very low expression level in the normal growing plants. But lots of PaCTLs showed upregulated expression profiles when the plants suffered different abiotic stress conditions. Among them, five PaCTLs responded to high temperature and exhibited significantly upregulate expression level. Correspondingly, many hormone responses, as well as biotic and abiotic stress elements were found in the promoters of PaCTLs by using cis-acting element analysis. These results provide a foundation for the exploration of PaCTLs’ function and enrich the evolutionary process of the CTL gene family.

Preliminary Study on the Sequencing of Whole Genomic Methylation and Transcriptome-Related Genes in Thyroid Carcinoma.

Thyroid carcinoma is the most prevalent endocrine cancer globally and the primary cause of cancer-related mortality. Epigenetic modifications are progressively being linked to metastasis. This study aimed to examine whole-genome DNA methylation patterns and the gene expression profiles in thyroid cancer tissue samples using a MethylationEPIC BeadChip (850K), RNA sequencing, and a targeted bisulfite sequencing assay. The results of the Illumina Infinium human methylation kit (850K) analyses identified differentially methylated CpG locations (DMPs) and differentially methylated CpG regions (DMRs) encompassing nearly the entire genome with high resolution and depth. Gene ontology and KEGG pathway analyses revealed that the genes associated with DMRs belonged to various domain-specific ontologies, including cell adhesion, molecule binding, and proliferation. The RNA-Seq study found 1627 differentially expressed genes, 1174 of which that were up-regulated and 453 of which that were down-regulated. The targeted bisulfite sequencing assay revealed that CHST2, DPP4, DUSP6, ITGA2, SLC1A5, TIAM1, TNIK, and ABTB2 methylation levels were dramatically lowered in thyroid cancer patients when compared to the controls, but GALNTL6, HTR7, SPOCD1, and GRM5 methylation levels were significantly raised. Our study revealed that the whole-genome DNA methylation patterns and gene expression profiles in thyroid cancer shed new light on the tumorigenesis of thyroid cancer.

Meta-analysis of whole-genome gene expression datasets assessing the effects of IDH1 and IDH2 mutations in isogenic disease models

Mutations in isocitrate dehydrogenase 1 (IDH1) and IDH2 are oncogenic drivers to a variable extent in several tumors, including gliomas, acute myeloid leukemia (AML), cholangiocarcinoma, melanoma, and thyroid carcinoma. The pathobiological effects of these mutations vary considerably, impeding the identification of common expression profiles. We performed an expression meta-analysis between IDH-mutant (IDHmut) and IDH-wild-type (IDHwt) conditions in six human and mouse isogenic disease models. The datasets included colon cancer cells, glioma cells, heart tissue, hepatoblasts, and neural stem cells. Among differentially expressed genes (DEGs), serine protease 23 (PRSS23) was upregulated in four datasets, i.e., in human colon carcinoma cells, mouse heart tissue, mouse neural stem cells, and human glioma cells. Carbonic anhydrase 2 (CA2) and prolyl 3-hydroxylase 2 (P3H2) were upregulated in three datasets, and SOX2 overlapping transcript (SOX2-OT) was downregulated in three datasets. The most significantly overrepresented protein class was termed intercellular signal molecules. An additional DEG set contained genes that were both up- and downregulated in different datasets and included oxidases and extracellular matrix structural proteins as the most significantly overrepresented protein classes. In conclusion, this meta-analysis provides a comprehensive overview of the expression effects of IDH mutations shared between different isogenic disease models. The generated dataset includes biomarkers, e.g., PRSS23 that may gain relevance for further research or clinical applications in IDHmut tumors.