Enriched Gene Ontology Categories Research Articles

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein

BackgroundGenetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood.ResultsTo investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes.ConclusionsThis study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

Phenotypic features, transcriptomic signatures and genomic regions associated with body ellipticity in the flatfish Senegalese sole

Body shape is a complex and plastic trait with a high impact on fish performance and commercialization. Shape is particularly relevant in flatfish with highly asymmetric bodies. In this study we investigated the skeleton features, transcriptomic profiles and genomic regions in two Senegalese sole groups from a full-sib family with positive (named as HE) or negative (LE) breeding values for body ellipticity (ELL). Soles from HE group were heavier and displayed more elongated bodies than LE. A skeleton analysis showed that HE had a higher number of vertebrae, longer vertebral bodies and intervertebral spaces and a lower incidence of skeletal fusions than LE. With respect to the gender effect, females showed a higher height of ventral body and some differences in the ratio urostyle to skeleton length than males for similar ELL values. The PCA analysis using 28 dimensionless morphological traits indicated that PC1 and PC2 explained 51.3% of the total variance and they were mainly associated with ELL and gender variation, respectively. Expression analysis using two factors (ELL and gender) revealed a major effect of gender in muscle, pterygiophores and liver with 6764, 6506 and 7559 differentially expressed transcripts (DETs) between sexes, respectively. In brain, only 16 DETs were identified. When HE and LE were compared controlling by gender, we identified the subset of DETs with conserved responses in both sexes or those regulated in a sex-specific manner in each tissue. A comparison across tissues identified 31 ubiquitous transcripts as differentially expressed in at least two of them. Functional analysis showed that liver and pterygiophores had the highest number of enriched gene ontology categories (24 and 14, respectively). Some of these enriched pathways were related to regulation of cell cycle, cell structure and shape, locomotory behavior, immune system, hormonal responses and homeostasis. A BSR-Seq analysis confirmed five significant genomic windows associated with ELL by using the tricube-smoothed G' and Fst statistics. Windows located in linkage groups (LG) 5 and LG14 included the malt1 and vtg3 genes, respectively, that were also detected within the subset of ubiquitous DETs. Moreover, malt1 appeared located very close to the bmpr1b gene within the significance hotspot of LG5 indicating that mediators of the inflammatory responses and BMP signaling pathway could be responsible for the differences found in the skeleton features and body ellipticity.

Transcriptomic landscape of tomato traditional long shelf-life landraces under low water regimes

‘Corbarino’ (COR) and ‘Lucariello’ (LUC) belong to the family of Mediterranean long shelf-life tomato landraces, producing high quality fruits under low water input cultivation regime in their traditional cultivation area. Understanding the morpho-physiological and molecular details of the peculiar drought stress tolerance of these two genotypes may be key to their valorization as breeding material. RNA sequencing of leaf samples of COR and LUC subjected to drought stress by water withholding in a semi-controlled greenhouse identified 3089 and 2135 differentially expressed genes respectively. These included COR- and LUC-specific annotated genes, as well as genes containing single nucleotide polymorphisms as compared to reference genome. Enriched Gene Ontology categories showed that categories such as response to water, oxidoreductase activity, nucleotide salvation and lipid biosynthesis-related processes were enriched among up-regulated DEGs. By contrast, growth and photosynthesis related genes were down-regulated after drought stress, consistent with leaf gas exchange and biomass accumulation measurements. Genes encoding cell wall degrading enzymes of the pectinase family were also down-regulated in drought stress conditions and upregulated in rewatering, indicating that cell wall composition/hardness is important for drought stress responses. Globally our results contribute to understanding the transcriptomic and physiological responses of representative tomato genotypes from Southern Italy, highlighting a promising set of genes to be investigated to improve tomato tolerance to drought.

IRAK2, an Immune and Radiation-Response Gene, Correlates with Advanced Disease Features but Predicts Higher Post-Irradiation Local Control in Non-Metastatic and Resected Oral Cancer Patients.

Gene Ontology (GO) analysis can provide a comprehensive function analysis for investigating genes, allowing us to identify the potential biological roles of genes. The present study conducted GO analysis to explore the biological function of IRAK2 and performed a case analysis to define its clinical role in disease progression and mediating tumor response to RT. Methods: We performed a GO enrichment analysis on the RNA-seq data to validate radiation-induced gene expression. A total of 172 I-IVB specimens from oral squamous cell carcinoma patients were collected for clinical analysis, from which IRAK2 expression was analyzed by immunohistochemistry. This was a retrospective study conducted between IRAK2 expression and the outcomes of oral squamous cell carcinoma patients after radiotherapy treatment. We conducted Gene Ontology (GO) analysis to explore the biological function of IRAK2 and performed a case analysis to define its clinical role in mediating tumor response to radiotherapy. GO enrichment analysis to validate radiation-induced gene expression was performed. Clinically, 172 stage I-IVB resected oral cancer patients were used to validate IRAK2 expression in predicting clinical outcomes. GO enrichment analysis showed that IRAK2 is involved in 10 of the 14 most enriched GO categories for post-irradiation biological processes, focusing on stress response and immune modulation. Clinically, high IRAK2 expression was correlated with adverse disease features, including pT3-4 status (p = 0.01), advanced overall stage (p = 0.02), and positive bone invasion (p = 0.01). In patients who underwent radiotherapy, the IRAK2-high group was associated with reduced post-irradiation local recurrence (p = 0.025) compared to the IRAK2-low group. IRAK2 plays a crucial role in the radiation-induced response. Patients with high IRAK2 expression demonstrated more advanced disease features but predicted higher post-irradiation local control in a clinical setting. These findings support IRAK2 as a potential predictive biomarker for radiotherapy response in non-metastatic and resected oral cancer patients.

Reconstruction and Differential Expression Profiling Core Target Analyses of the circRNA-miRNA-mRNA Network Based on Competitive Endogenous RNAs in Ulcerative Colitis

Ulcerative colitis (UC) is a common autoimmune disease worldwide. Circular RNA (circRNA) is a type of noncoding ribonucleic acids (ncRNAs). In addition to their roles in numerous biological processes, circRNAs are also linked to a vast range of diseases including UC. Although previous studies have examined many circRNAs, the physiological and pathological roles of the circRNA-associated competing endogenous RNA (ceRNA) network in UC remain unclear. Thus, we constructed a circRNA-miRNA-mRNA network based on the ceRNA hypothesis by analyzing data from the National Center for Biotechnology Information Gene Expression Omnibus (NCBI-GEO) database. Genes with higher degree values than others in the ceRNA network were selected as central nodes when constructing the corresponding core subnetworks. To fully understand the biological function of the ceRNA network, we entered all differentially expressed mRNAs (DEmRNAs) from the ceRNA network into the Database for Annotation and Integrated Discovery (DAVID), which was used to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. We further entered DEmRNAs into the STRING database for protein-protein interaction (PPI) network analysis. The results elucidated that the ceRNA network comprised 403 circRNA nodes, 5 miRNA nodes, 138 mRNA nodes, and 559 edges. Three core ceRNA subnetworks centered on hsa-miR-342-3p, hsa-miR-199a-5p, and hsa-miR-142-3p were reconstructed in this study. GO and KEGG enrichment analyses identified 167 enriched GO categories and 14 enriched KEGG pathway terms. The core PPI network was composed of 15 core targets, of which CD44, HIF1A, and MMP2 were the most significant. In summary, 3 hub miRNAs (hsa-miR-342-3p, hsa-miR-199a-5p, hsa-miR-142-3p) and 3 hub genes (CD44, HIF1A, and MMP2) might play an important role in the development of UC. These hub nodes, first proposed here, might also be used as potential diagnostic markers and therapeutic targets.

Epigenetics: new insights into postoperative adhesion development.

The link between post-operative adhesion development and epigenetic modifications is important in understanding the mechanism behind their formation. The purpose of this study was to determine whether epigenetic differences exist between primary fibroblasts of normal peritoneum and adhesion tissues isolated from the same patient(s). DNA from fibroblasts isolated from normal peritoneum and adhesion tissues was isolated using Qiagen's EZ1 Advanced Kit. Methylation patterns of genes were quantified and compared in both cell lines using the Infinium Human Methylation 27 Beadchip system. A total of 7364 genes had been found to manifest significantly different DNA methylation levels in adhesion fibroblasts as compared to normal peritoneal fibroblasts (p<0.01). A total of 1685 genes were found to have increased DNA methylation by 50% in adhesion compared to peritoneal fibroblasts, and were enriched in Gene Ontology categories, Glycoprotein, and Defense Response. Furthermore, 1287 genes were found to have decreased DNA methylation patterns with enriched Gene Ontology categories, "Homeobox", and Transcription Factor Activity in adhesion fibroblasts. Epigenetic differences in fibroblasts isolated from normal peritoneum and adhesion tissues were observed. Future studies focusing on the precise role of these genes in the development of post operative adhesions will allow us to more fully appreciate regulatory mechanisms leading to adhesion development, thereby establishing targets for therapeutic interventions to prevent or limit adhesion development.

RNA-Seq of amniotic fluid cell-free RNA: a discovery phase study of the pathophysiology of congenital cytomegalovirus infection

Congenital cytomegalovirus infection is the most common perinatal infection and a significant cause of sensorineural hearing loss, cerebral palsy, and neurodevelopmental disability. There is a paucity of human gene expression studies examining the pathophysiology of cytomegalovirus infection. This study aimed to perform a whole transcriptomic assessment of amniotic fluid from pregnancies with live fetuses to identify differentially expressed genes and enriched Gene Ontology categories associated with congenital cytomegalovirus infection. Amniotic fluid supernatant was prospectively collected from pregnant women undergoing amniocentesis for suspected congenital cytomegalovirus infection because of first-trimester maternal primary infection or ultrasound features suggestive of fetal infection. Women who had received therapy to prevent fetal infection were excluded. Congenital cytomegalovirus infection was diagnosed via viral polymerase chain reaction of amniotic fluid; cytomegalovirus-infected fetuses were paired with noninfected controls, matched for gestational age and fetal sex. Paired-end RNA sequencing was performed on amniotic fluid cell-free RNA with the Novaseq 6000 at a depth of 30 million reads per sample. Following quality control and filtering, reads were mapped to the human genome and counts summarized across genes. Differentially expressed genes were identified using 2 approaches: voomWithQualityWeights in conjunction with limma and RUVSeq with edgeR. Genes with a false discovery rate <0.05 were considered statistically significant. Differential exon use was analyzed using DEXSeq. Functional analysis was performed using gene set enrichment analysis and Ingenuity Pathway Analysis. Manual curation of differentially regulated genes was also performed. Amniotic fluid samples were collected from 50 women; 16 (32%) had congenital cytomegalovirus infection confirmed by polymerase chain reaction. After excluding 3 samples without matched controls, 13 cytomegalovirus-infected samples collected at 18 to 23 weeks and 13 cytomegalovirus-negative gestation-matched controls were submitted for RNA sequencing and analysis (N=26). Ten of the 13 pregnancies with cytomegalovirus-infected fetuses had amniocentesis because of serologic evidence of maternal primary infection with normal fetal ultrasound, and 3 had amniocentesis because of ultrasound abnormality suggestive of cytomegalovirus infection. Four cytomegalovirus-infected pregnancies ended in termination (n=3) or fetal death (n=1), and 9 resulted in live births. Pregnancy outcomes were available for 11 of the 13 cytomegalovirus-negative controls; all resulted in live births of clinically-well infants. Differential gene expression analysis revealed 309 up-regulated and 32 down-regulated genes in the cytomegalovirus-infected group compared with the cytomegalovirus-negative group. Gene set enrichment analysis showed significant enrichment of multiple Gene Ontology categories involving the innate immune response to viral infection and interferon signaling. Of the 32 significantly down-regulated genes, 8 were known to be involved in neurodevelopment and preferentially expressed by the brain. Six specific cellular restriction factors involved in host defense to cytomegalovirus infection were up-regulated in the cytomegalovirus-infected group. Ingenuity Pathway Analysis predicted the activation of pathways involved in progressive neurologic disease and inflammatory neurologic disease. In this next-generation sequencing study, we revealed new insights into the pathophysiology of congenital cytomegalovirus infection. These data on the up-regulation of the intraamniotic innate immune response to cytomegalovirus infection and the dysregulation of neurodevelopmental genes may inform future approaches to developing prognostic markers and assessing fetal responses to in utero therapy.

Expression of lncRNAs in children with pancreaticobiliary maljunction: functional analysis and potential biomarkers.

Pancreaticobiliary maljunction (PBM) leads to higher rates of complications, including cholangitis, pancreatitis, and malignancies. The aim of the present study was to investigate the expression profile of long non-coding RNAs (lncRNAs) and their potential role as biomarkers in children with pancreaticobiliary maljunction. The differential expression of lncRNAs and messenger RNA (mRNAs) from pediatric patients with pancreaticobiliary maljunction and control subjects was analyzed using a commercial microarray and later validated with qRT-PCR. The potential biological functions of differentially expressed genes were explored based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment. The ability of potential lncRNA biomarkers to predict pancreaticobiliary maljunction was assessed based on the area under the receiver operating characteristic curve (AUC). There were 2915 mRNAs and 173 lncRNAs upregulated, and 2121 mRNAs and 316 lncRNAs downregulated in PBM cases compared to controls. The enriched Gene Ontology categories associated with differentially expressed mRNAs were extracellular matrix, extracellular region, and kinetochore. The most enriched Kyoto Encyclopedia pathway was protein digestion and absorption, which was associated with cancer and PI3K-Akt signaling. Analysis of cis- and trans-target genes predicted that a single lncRNA was able to regulate several mRNAs. The qRT-PCR results for NR_110876, NR_132344, XR_946886, and XR_002956345 were consistent with the microarray results, and the difference was statistically significant for NR_132344, XR_946886, and XR_002956345 (p < 0.05). AUC was significant only for XR_946886 (0.837, p < 0.001). Our results implicate lncRNAs in common bile duct pathogenesis in PBM, and they identify XR_946886 as a potential biomarker for the disease.

Coexpression gene network analysis of cold‐tolerant Solanum commersonii reveals new insights in response to low temperatures

AbstractAmong abiotic stressors, cold is one of the most harmful for the cultivated potato (Solanum tuberosum L.), a frost‐sensitive crop. RNA sequencing (RNA‐seq) profiling of two different clones of wild potato (S. commersonii Dun.) contrasting in their capacity to withstand low temperatures revealed a higher number of differentially expressed genes (DEGs) under nonacclimated conditions (NAC) in tolerant clone cmm1T vs. the susceptible cmm6‐6 (1,002 and 8,055 DEGs, respectively). By contrast, the number of DEGs was much more comparable when both genotypes were under acclimated conditions (AC). Indeed, a total of 5,650 and 8,936 DEGs were detected in the tolerant genotype vs. the susceptible. Gene ontology (GO) classification under NAC showed a significant role for transcription regulation, lignin catabolic genes, and regulation of plant type secondary cell wall in the cold‐tolerant genotypes, suggesting an important role in conferring tolerance response. By contrast, response to stress and response to stimuli were enriched GO categories in both clones under AC. Unsigned weighted correlation networks analysis (WGCNA) allowed identification of coexpressed hub genes with possible main regulatory functions and major impacts on the phenotype. Among those identified, we clarified the role of CBF4. This gene showed contrasting expression profiles in the two clones under NAC, being induced in cold‐tolerant cmm1T but suppressed in susceptible cmm6‐6. By contrast, under AC, CBF4 was upregulated in both clones. Our study provides a global understanding of mechanisms involved following exposure to NAC and AC in S. commersonii. The mechanisms described here will inform future investigations for detailed validation in studies regarding cold tolerance in plants.

Abstract 2695: Identification of an age-related breast cancer gene expression signature with strong prognostic value

Abstract Background and Objective: Breast cancer (BC) diagnosed at ages &amp;lt;40 years is known to present more aggressive tumor phenotypes and poorer clinical outcome compared to older BC patients. To gain a better understanding of the possible age-related expression difference in BC promoting genes, we aimed to identify transcriptional alterations supporting cancer progression in BC of the young. Methods: We studied mRNA gene expression data from METABRIC discovery (n=997), and validation cohort (n=995), and explored the transcriptional patterns associated with BC of the young. We investigated differentially expressed genes (DEGs) between primary BC for patients below and above 40 years at time of diagnosis and explored gene sets (by Gene Set Enrichment Analysis) enriched in BC of the young. Protein-protein interaction (PPI) networks were explored by the STRING database and visualized in Cytoscape software. Hub genes were identified by MCODE and CytoHubba plugins. Results: Among genes differentially expressed between BC of the young and older, 91 DEGs were commonly (for the two cohorts) upregulated and 99 commonly downregulated in the young (Fold change ≥1.5/≤ -1.5). Oncogenic- and proliferation associated signatures were enriched in the young BC, including signatures reflecting KRAS, MTOR and MYC. Proliferation and cell cycle processes were the dominating enriched gene ontology categories (FDR &amp;lt;10%). Six hub genes presenting highest PPI network connectivity were identified, described to be associated with the cell cycle and cell proliferation. A signature score was made by summarizing the hub genes expression values. High signature score was significantly associated with high tumor size, high histologic grade, lymph node metastasis, ER negativity, basal-like, and HER2 enriched subtypes (P &amp;lt;0.001). The signature score showed strong correlation with the proliferation signatures OncotypeDx and PCNA (ρ=0.90-0.96, P &amp;lt;0.001). Our signature score associated with reduced cancer specific survival (P &amp;lt;0.001), also when adjusted for tumor size, histologic grade, and lymph node status in multivariate analyses in both cohorts (HR: 1.074-1.081, 95% CI: 1.046-1.109, P ≤0.003). The signature score maintained independent, significant association with survival when adjusting for the traditional clinico-pathologic variables in OncotypeDx-low tumors in discovery (HR: 1.077, 95% CI: 1.034-1.122) and validation cohort (HR: 1.154, 95% CI: 1.065-1.251), both P &amp;lt;0.001. Conclusion: By genome-wide gene expression analyses, the current study provides new insights into age-related gene expression alterations in breast cancer. We demonstrate enrichment of oncogenic signaling and evidence of higher tumor cell proliferation in young BC patients, and identify a gene expression signature reflecting tumor proliferation, aggressive tumor features and reduced survival, also in subsets of low OncotypeDx score. Citation Format: Lise Martine Ingebriktsen, Lars Andreas Akslen, Elisabeth Wik. Identification of an age-related breast cancer gene expression signature with strong prognostic value [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2695.

Multi-gene signature of microcalcification and risk prediction among Taiwanese breast cancer

Microcalcification is one of the most common radiological and pathological features of breast ductal carcinoma in situ (DCIS), and to a lesser extent, invasive ductal carcinoma. We evaluated messenger RNA (mRNA) transcriptional profiles associated with ectopic mammary mineralization. A total of 109 breast cancers were assayed with oligonucleotide microarrays. The associations of mRNA abundance with microcalcifications and relevant clinical features were evaluated. Microcalcifications were present in 86 (79%) patients by pathological examination, and 81 (94%) were with coexistent DCIS, while only 13 (57%) of 23 patients without microcalcification, the invasive diseases were accompanied with DCIS (χ2-test, P < 0.001). There were 69 genes with differential mRNA abundance between breast cancers with and without microcalcifications, and 11 were associated with high-grade (comedo) type DCIS. Enriched Gene Ontology categories included glycosaminoglycan and aminoglycan metabolic processes and protein ubiquitination, indicating an active secretory process. The intersection (18 genes) of microcalcificaion-associated and DCIS-associated genes provided the best predictive accuracy of 82% with Bayesian compound covariate predictor. Ten genes were further selected for prognostic index score construction, and five-year relapse free survival was 91% for low-risk and 83% for high-risk group (log-rank test, P = 0.10). Our study suggested that microcalcification is not only the earliest detectable radiological sign for mammography screening but the phenomenon itself may reflect the underling events during mammary carcinogenesis. Future studies to evaluate the prognostic significance of microcalcifications are warranted.

Intratumoral Cancer Chemotherapy with a Carrier-Based Immunogenic Cell-Death Eliciting Platinum (IV) Agent.

A carrier-based, immunogenic cell death (ICD)-eliciting platinum(IV) chemotherapeutic agent was synthesized via complexation between an axially derivatized Pt(IV)-tocopherol and hyaluronan (HA)-tocopherol nanocarrier. The resultant HA-Pt(IV) complex demonstrated antiproliferative activity and induced calreticulin translocation, an indicator of ICD, in murine and human head and neck cancer (HNC) cells. The intratumorally administered HA-Pt(IV) treatments were tolerable and efficacious in both immunocompetent and immunodeficient mice with HNC, partially because of the direct cytotoxicity. Superior efficacy and survival were observed in the immunocompetent group, suggesting a possible Pt(IV)-induced immunological response, which would only manifest in animals with an intact immune system. Subsequent imaging of tumor tissues demonstrated increased macrophage infiltration in the HA-Pt(IV)-treated tumors compared to the nontreated controls and the cisplatin-treated tumors, suggesting favorable inflammatory activation. RNA sequencing of HA-Pt(IV)-treated tumors indicated that carbohydrate and vitamin metabolisms were the most important Kyoto Encyclopedia of Genes and Genomes pathways, and molecular function, biological process, and cellular component were highly enriched gene ontology categories.

Global gene expression profiling of blast lung injury of goats exposed to shock wave

PurposeBlast lung injury (BLI) is the most common damage resulted from explosion-derived shock wave in military, terrorism and industrial accidents. However, the molecular mechanisms underlying BLI induced by shock wave are still unclear. MethodsIn this study, a goat BLI model was established by a fuel air explosive power. The key genes involved in were identified. The goats of the experimental group were fixed on the edge of the explosion cloud, while the goats of the control group were 3 km far away from the explosive environment. After successful modeling for 24 h, all the goats were sacrificed and the lung tissue was harvested for histopathological observation and RNA sequencing. Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) analysis were performed to identify the main enriched biological functions of differentially expressed genes (DEGs). Quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the consistency of gene expression. ResultsOf the sampled goat lungs, 895 genes were identified to be significantly differentially expressed, and they were involved in 52 significantly enriched GO categories. KEGG analysis revealed that DEGs were highly enriched in 26 pathways, such as cytokine-cytokine receptor interaction, antifolate resistance, arachidonic acid metabolism, amoebiasis and bile secretion, JAK-STAT, and IL-17 signaling pathway. Furthermore, 15 key DEGs involved in the biological processes of BLI were confirmed by qRT-PCR, and the results were consistent with RNA sequencing. ConclusionGene expression profiling provide a better understanding of the molecular mechanisms of BLI, which will help to set strategy for treating lung injury and preventing secondary lung injury induced by shock wave.

Abstract 3661: Association of mammographic density with blood DNA methylation

Abstract Background: Mammographic density is a strong risk factor for breast cancer. While numerous genetic variations have been shown to be associated with mammographic density, an epigenetic signature of mammographic density has not yet been identified. The objective of this study was to identify a DNA methylation signature associated with mammographic density. Methods: Blood samples were collected from 197 postmenopausal women 45-65 years old with no personal history of breast cancer. For each participant, mammographic density according to the BI-RADs density classification was obtained from the most recent mammogram report. White blood cell DNA methylation at more than 850,000 CpG sites was measured using the Illumina HumanMethylationEPIC BeadChip. Differential methylation was assessed using genome-wide, probe-level, and regional approaches after filtering and normalization using the recommended steps for DNA methylation microarray data. All analyses were adjusted for appropriate confounding variables. Results: On average, genome-wide DNA methylation was higher in women with elevated mammographic density. Additionally, DNA methylation at 343 individual sites and in 34 regions was significantly associated with mammographic density after false discovery rate correction. Most density-associated sites were hypermethylated. The top differentially methylated sites were located near cell cycle regulatory genes (ICMT, DESI2, TCF21/TARID, ROBO3), as well as those encoding cytoskeletal or plasma membrane proteins (PALLD, ANO7, ATP8B3, MTUS2). We observed hypermethylation near the promoters of genes including RUFY1 and SLFN12 and within the bodies of genes including LCN6, SNED1, and LRFN1. Functional analysis of the differentially methylated sites revealed enrichment of gene ontology categories including angiogenesis and regulation of cell adhesion, as well as apoptotic signaling in response to DNA damage. Conclusions: White blood cell DNA methylation differed according to mammographic density in a sample of postmenopausal women. Probe-based and regional analysis identified differentially methylated regions located near genes implicated in cancer. These differentially methylated regions may serve as a signature of mammographic density, and future studies of the genes implicated may provide a deeper understanding of the connection between mammographic density and breast cancer risk. Citation Format: Rachel McFarland Lucia, Wei-Lin Huang, Andrea Alvarez, Irene Masunaka, Argyrios Ziogas, Deborah Goodman, Andrew O. Odegaard, Trina M. Norden-Krichmar, Hannah Lui Park. Association of mammographic density with blood DNA methylation [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3661.

Pseudophoenix ekmanii (Arecaceae) seeds at suboptimal temperature show reduced imbibition rates and enhanced expression of genes related to germination inhibition.

Pseudophoenix ekmanii is a critically endangered palm species that can be found in the southeast of the Dominican Republic. The temperatures to which P. ekmanii seeds are exposed upon dispersal range from 32 to 23°C (max and min) and can reach a low of approximately 20°C in January. Our aim was to analyse the effect of suboptimal (20°C) and optimal (30°C) germination temperature on seed imbibition, moisture content, embryo growth and gene expression patterns in this tropical palm species. Seed imbibition was tracked using whole seeds, while moisture content was assessed for individual seed sections. Embryo and whole seed size were measured before and after full imbibition. For transcriptome sequencing, mRNA was extracted from embryo tissues only and the resulting reads were mapped against the Elaeis guineensis reference genome. Differentially expressed genes were identified after statistical analysis and subsequently probed for enrichment of Gene Ontology categories 'Biological process' and 'Cellular component'. Our results show that prolonged exposure to 20°C slows whole seed and embryo imbibition and causes germination to be both delayed and inhibited. Embryonic transcriptome patterns associated with the negative regulation of germination by suboptimal temperature include up-regulation of ABA biosynthesis genes, ABA-responsive genes, as well as other genes previously related to physiological dormancy and inhibition of germination. Thus, our manuscript provides the first insights into the gene expression patterns involved in the response to suboptimal temperature during seed imbibition in a tropical palm species.

Expression Profile Analysis of Selenium-Related Genes in Peripheral Blood Mononuclear Cells of Patients with Keshan Disease.

Keshan disease (KD) is an endemic cardiomyopathy, which mainly occurs in China. Selenium deficiency is believed to play an important role in the pathogenesis of KD, but the molecular mechanism of selenium-induced damage remains unclear. To identify the key genes involved in selenium-induced damage, we compared the expression profiles of selenium-related genes between patients with KD and normal controls. Total RNA was isolated, amplified, labeled, and hybridized to Agilent human 4 × 44 K whole genome microarrays. Selenium-related genes were screened using the Comparative Toxicogenomics Database. The microarray data were subjected to single-gene and gene ontology (GO) expression analysis using R Studio and Gene Set Enrichment Analysis (GSEA) software. Quantitative real-time PCR was conducted to validate the microarray results. We identified 16 upregulated and 11 downregulated selenium-related genes in patients. These genes are involved in apoptosis, metabolism, transcription regulation, ion transport, and growth and development. Of the significantly enriched GO categories in KD patients, we identified four apoptosis-related, two metabolism-related, four growth and development-related, and four ion transport-related GOs. Based on our results, we suggest that selenium might contribute to the development of KD through dysfunction of selenium-related genes involved in apoptosis, metabolism, ion transport, and growth and development in the myocardium.

Transcriptome analysis of the curry tree (Bergera koenigii L., Rutaceae) during leaf development

The curry tree (Bergera koenigii L.) is a widely cultivated plant used in South Asian cooking. Next-generation sequencing was used to generate the transcriptome of the curry leaf to detect changes in gene expression during leaf development, such as those genes involved in the production of oils which lend the leaf its characteristic taste, aroma, and medicinal properties. Using abundance estimation (RSEM) and differential expression analysis, genes that were significantly differentially expressed were identified. The transcriptome was annotated with BLASTx using the non-redundant (nr) protein database, and Gene Ontology (GO) terms were assigned based on the top BLAST hit using Blast2GO. Lastly, functional enrichment of the assigned GO terms was analyzed for genes that were significantly differentially expressed. Of the most enriched GO categories, pathways involved in cell wall, membrane, and lignin synthesis were found to be most upregulated in immature leaf tissue, possibly due to the growth and expansion of the leaf tissue. Terpene synthases, which synthesize monoterpenes and sesquiterpenes, which comprise much of the curry essential oil, were found to be significantly upregulated in mature leaf tissue, suggesting that oil production increases later in leaf development. Enzymes involved in pigment production were also significantly upregulated in mature leaves. The findings were based on computational estimates of gene expression from RNA-seq data, and further study is warranted to validate these results using targeted techniques, such as quantitative PCR.

GeneMANIA update 2018.

GeneMANIA (http://genemania.org) is a flexible user-friendly web site for generating hypotheses about gene function, analyzing gene lists and prioritizing genes for functional assays. Given a query gene list, GeneMANIA finds functionally similar genes using a wealth of genomics and proteomics data. In this mode, it weights each functional genomic dataset according to its predictive value for the query. Another use of GeneMANIA is gene function prediction. Given a single query gene, GeneMANIA finds genes likely to share function with it based on their interactions with it. Enriched Gene Ontology categories among this set can point to the function of the gene. Nine organisms are currently supported (Arabidopsis thaliana, Caenorhabditis elegans, Danio rerio, Drosophila melanogaster, Escherichia coli, Homo sapiens, Mus musculus, Rattus norvegicus and Saccharomyces cerevisiae). Hundreds of data sets and hundreds of millions of interactions have been collected from GEO, BioGRID, IRefIndex and I2D, as well as organism-specific functional genomics data sets. Users can customize their search by selecting specific data sets to query and by uploading their own data sets to analyze. We have recently updated the user interface to GeneMANIA to make it more intuitive and make more efficient use of visual space. GeneMANIA can now be used effectively on a variety of devices.

Deep sequencing-based characterization of transcriptome of Pyrus ussuriensis in response to cold stress

Pyrus ussuriensis is extremely cold hardy when fully acclimated, but knowledge relevant to the molecular mechanisms underlying this economically valuable trait is still limited so far. In this study, global transcriptome profiles of Pyrus ussuriensis under cold conditions (4 °C) over a time course were generated by high-throughput sequencing. In total, >57,121,199 high quality clean reads were obtained with approximately 11.0 M raw data for each library. Among them, the values of 66.84%-72.03% of clean reads in the digital transcript abundance measurement could be well mapped to the pear genome database, resulting in the identification of 8544 differentially expressed genes (DEGs) having 43 Gene Ontology (GO) terms and 17 clusters of orthologous groups (COG) involved in 385 Kyoto Encyclopedia of Genes and Genomes (KEGG) defined pathways. These comprised 3124 (1033 up-regulated, 2091 down-regulated), 1243 (729 up-regulated, 514 down-regulated), and 750 (458 up-regulated, 292 down-regulated) genes from the cold-treated samples at 5, 12 and 24 h, respectively. The accuracy of the RNA-Seq derived transcript expression data was validated by analyzing the expression patterns of 16 DGEs by quantitative real-time PCR. Plant-pathogen interaction, plant hormone signal transduction, Photosynthesis, signal transduction, innate immune response and response to biotic stimulus were the most significantly enriched GO categories among in the DEGs. A total of 335 transcription factors were shown to be cold responsive. In addition, a number of genes involved in the catabolism and signaling of hormones were significantly affected by the cold stress. The RNA-Seq and digital expression profiling provides valuable insights into the understanding the molecular events related to cold responses in Pyrus ussuriensis and dataset may help guide future identification and functional analysis of potential genes that are important for enhancing cold hardiness.

Array probe density and pathobiological relevant CpG calling bias in human disease and physiological DNA methylation profiling.

The HumanMethylation450 BeadChip array (450K; Infinium) is a widely used tool in epigenomics. A recognized concern in the 450K platform is the potential effect of the number of probes/gene (PG) on ranking differentially methylated (DM) CpGs (DM-CpGs) before testing for enrichment of gene ontology categories. We previously showed in a fatty acid (FA)-induced DNA methylation profiling study that when DM-CpGs are ranked by the number of called DM-CpGs-to-PG ratio, the 150 top-ranking gene list is enriched in pathways that overlap with the corresponding Affymetrix array-based expression data. In this study, a comparative analysis of thirteen 450K-based studies representing FA-stimulated cellular models, aging, diseased and normal tissues, revealed that the 150 top-ranking DM-CpGs are in high PG genes. This points to a significant false-negative rate in the low PG gene set when delta-beta-based ranking is performed. We show that PG is not related to the density of methylation-prone sites, as it does not follow gene length or GC content. Conversely, ranking genes by the number of DM-CpGs-to-PG ratio and analysing the 150 top-ranking entries yields significantly enriched gene disease- or tissue-specific function categories that are increased both in number and in the degree of overlap with expression data compared with delta-beta-only ranking or to the previously published gometh-based pipeline. The 15 top-ranking loci list is also significantly enriched in non-coding RNAs, a greatly underrepresented transcript type in 450K. In summary, the proposed simple normalization method yields pathobiologically relevant DM-CpGs. This method is relevant for the newly developed MethylationEPIC (Infinium) microarray.