RNA-seq analysis identifies age-dependent changes in expression of mRNAs - encoding N-glycosylation pathway enzymes in mouse gonadotropes.

Follicle-stimulating hormone (FSH) is a glycoprotein that is assembled as a heterodimer of α/ß subunits in gonadotropes. Each subunit contains two N-glycan chains. Our previous in vivo genetic studies identified that at least one N-glycan chain must be present on the FSHß subunit for efficient FSH dimer assembly and secretion. Moreover, macroheterogeneity observed uniquely on human FSHß results in ratiometric changes in age-specific FSH glycoforms, particularly during menopausal transition. Despite the recognition of many prominent roles of sugars on FSH including dimer assembly and secretion, serum half-life, receptor binding and signal transduction, the N-glycosylation machinery in gonadotropes has never been defined. Here, we used a mouse model in which gonadotropes are GFP-labeled in vivo and achieved rapid purification of GFP+ gonadotropes from pituitaries of female mice at reproductively young, middle, and old ages. We identified by RNA-seq analysis 52 mRNAs encoding N-glycosylation pathway enzymes expressed in 3- and 8-10-month-old mouse gonadotropes. We hierarchically mapped and localized the enzymes to distinct subcellular organelles within the N-glycosylation biosynthetic pathway. Of the 52 mRNAs, we found 27 mRNAs are differentially expressed between the 3- and 8-10-month old mice. We subsequently selected 8 mRNAs which showed varying changes in expression for confirmation of abundance in vivo via qPCR analysis, using more expanded aging time points with distinct 8-month and 14-month age groups. Real time qPCR analysis indicated dynamic changes in expression of N-glycosylation pathway enzyme-encoding mRNAs across the life span. Notably, computational analysis predicted the promoters of genes encoding these 8 mRNAs contain multiple high probability binding sites for estrogen receptor-1 and progesterone receptor. Collectively, our studies define the N-glycome and identify age-specific dynamic changes in mRNAs encoding N-glycosylation pathway enzymes in mouse gonadotropes. Our studies suggest the age-related decline in ovarian steroids may regulate expression of N-glycosylation enzymes in mouse gonadotropes and explain the age-related N-glycosylation shift previously observed on human FSHß subunit in pituitaries of women.

Molecular Subtyping and Survival Analysis of Osteosarcoma Reveals Prognostic Biomarkers and Key Canonical Pathways.

Osteosarcoma (OS) is a common bone malignancy in children and adolescents. Although histological subtyping followed by improved OS treatment regimens have helped achieve favorable outcomes, a lack of understanding of the molecular subtypes remains a challenge to characterize its genetic heterogeneity and subsequently to identify diagnostic and prognostic biomarkers for developing effective treatments. In the present study, global analysis of DNA methylation, and mRNA and miRNA gene expression in OS patient samples were correlated with their clinical characteristics. The mucin family of genes, MUC6, MUC12, and MUC4, were found to be highly mutated in the OS patients. Results revealed the enrichment of molecular pathways including Wnt signaling, Calcium signaling, and PI3K-Akt signaling in the OS tumors. Survival analyses showed that the expression levels of several genes such as RAMP1, CRIP1, CORT, CHST13, and DDX60L, miRNAs and lncRNAs were associated with survival of OS patients. Molecular subtyping using Cluster-Of-Clusters Analysis (COCA) for mRNA, lncRNA, and miRNA expression; DNA methylation; and mutation data from the TARGET dataset revealed two distinct molecular subtypes, each with a distinctive gene expression profile. Between the two subtypes, three upregulated genes, POP4, HEY1, CERKL, and seven downregulated genes, CEACAM1, ABLIM1, LTBP2, ISLR, LRRC32, PTPRF, and GPX3, associated with OS metastasis were found to be differentially regulated. Thus, the molecular subtyping results provide a strong basis for classification of OS patients that could be used to develop better prognostic treatment strategies.

Cell-free DNA in maternal blood and artificial intelligence: accurate prenatal detection of fetal congenital heart defects.

BACKGROUND: DNA cytosine nucleotide methylation (epigenomics and epigenetics) is an important mechanism for controlling gene expression in cardiac development. Combined artificial intelligence and whole-genome epigenomic analysis of circulating cell-free DNA in maternal blood has the potential for the detection of fetal congenital heart defects. OBJECTIVE: This study aimed to use genome-wide DNA cytosine methylation and artificial intelligence analyses of circulating cell-free DNA for the minimally invasive detection of fetal congenital heart defects. STUDY DESIGN: In this prospective study, whole-genome cytosine nucleotide methylation analysis was performed on circulating cell-free DNA using the Illumina Infinium MethylationEPIC BeadChip array. Multiple artificial intelligence approaches were evaluated for the detection of congenital hearts. The Ingenuity Pathway Analysis program was used to identify gene pathways that were epigenetically altered and important in congenital heart defect pathogenesis to further elucidate the pathogenesis of isolated congenital heart defects. RESULTS: There were 12 cases of isolated nonsyndromic congenital heart defects and 26 matched controls. A total of 5918 cytosine nucleotides involving 4976 genes had significantly altered methylation, that is, a P value of <.05 along with ≥5% whole-genome cytosine nucleotide methylation difference, in congenital heart defect cases vs controls. Artificial intelligence analysis of the methylation data achieved excellent congenital heart defect predictive accuracy (areas under the receiver operating characteristic curve, ≥0.92). For example, an artificial intelligence model using a combination of 5 whole-genome cytosine nucleotide markers achieved an area under the receiver operating characteristic curve of 0.97 (95% confidence interval, 0.87-1.0) with 98% sensitivity and 94% specificity. We found epigenetic changes in genes and gene pathways involved in the following important cardiac developmental processes: "cardiovascular system development and function," "cardiac hypertrophy," "congenital heart anomaly," and "cardiovascular disease." This lends biologic plausibility to our findings. CONCLUSION: This study reported the feasibility of minimally invasive detection of fetal congenital heart defect using artificial intelligence and DNA methylation analysis of circulating cell-free DNA for the prediction of fetal congenital heart defect. Furthermore, the findings supported an important role of epigenetic changes in congenital heart defect development.

mintRULS: Prediction of miRNA-mRNA Target Site Interactions Using Regularized Least Square Method.

Identification of miRNA-mRNA interactions is critical to understand the new paradigms in gene regulation. Existing methods show suboptimal performance owing to inappropriate feature selection and limited integration of intuitive biological features of both miRNAs and mRNAs. The present regularized least square-based method, mintRULS, employs features of miRNAs and their target sites using pairwise similarity metrics based on free energy, sequence and repeat identities, and target site accessibility to predict miRNA-target site interactions. We hypothesized that miRNAs sharing similar structural and functional features are more likely to target the same mRNA, and conversely, mRNAs with similar features can be targeted by the same miRNA. Our prediction model achieved an impressive AUC of 0.93 and 0.92 in LOOCV and LmiTOCV settings, respectively. In comparison, other popular tools such as miRDB, TargetScan, MBSTAR, RPmirDIP, and STarMir scored AUCs at 0.73, 0.77, 0.55, 0.84, and 0.67, respectively, in LOOCV setting. Similarly, mintRULS outperformed other methods using metrics such as accuracy, sensitivity, specificity, and MCC. Our method also demonstrated high accuracy when validated against experimentally derived data from condition- and cell-specific studies and expression studies of miRNAs and target genes, both in human and mouse.

Ecdysoneless Overexpression Drives Mammary Tumorigenesis through Upregulation of C-MYC and Glucose Metabolism.

Ecdysoneless (ECD) protein is essential for embryogenesis, cell-cycle progression, and cellular stress mitigation with an emerging role in mRNA biogenesis. We have previously shown that ECD protein as well as its mRNA are overexpressed in breast cancer and ECD overexpression predicts shorter survival in patients with breast cancer. However, the genetic evidence for an oncogenic role of ECD has not been established. Here, we generated transgenic mice with mammary epithelium-targeted overexpression of an inducible human ECD transgene (ECDTg). Significantly, ECDTg mice develop mammary hyperplasia, preneoplastic lesions, and heterogeneous tumors with occasional lung metastasis. ECDTg tumors exhibit epithelial to mesenchymal transition and cancer stem cell characteristics. Organoid cultures of ECDTg tumors showed ECD dependency for in vitro oncogenic phenotype and in vivo growth when implanted in mice. RNA sequencing (RNA-seq) analysis of ECDTg tumors showed a c-MYC signature, and alterations in ECD levels regulated c-MYC mRNA and protein levels as well as glucose metabolism. ECD knockdown-induced decrease in glucose uptake was rescued by overexpression of mouse ECD as well as c-MYC. Publicly available expression data analyses showed a significant correlation of ECD and c-MYC overexpression in breast cancer, and ECD and c-MYC coexpression exhibits worse survival in patients with breast cancer. Taken together, we establish a novel role of overexpressed ECD as an oncogenesis driver in the mouse mammary gland through upregulation of c-MYC-mediated glucose metabolism. IMPLICATIONS: We demonstrate ECD overexpression in the mammary gland of mice led to the development of a tumor progression model through upregulation of c-MYC signaling and glucose metabolism.

Maternal opioid use disorder: Placental transcriptome analysis for neonatal opioid withdrawal syndrome.

Excessive prenatal opioid exposure may lead to the development of Neonatal Opioid Withdrawal Syndrome (NOWS). RNA-seq was done on 64 formalin-fixed paraffin-embedded placental tissue samples from 32 mothers with opioid use disorder, with newborns with NOWS that required treatment, and 32 prenatally unexposed controls. We identified 93 differentially expressed genes in the placentas of infants with NOWS compared to unexposed controls. There were 4 up- and 89 downregulated genes. Among these, 7 genes CYP1A1, APOB, RPH3A, NRXN1, LINC01206, AL157396.1, UNC80 achieved an FDR p-value of <0.01. The remaining 87 genes were significant with FDR p-value <0.05. The 4 upregulated, CYP1A1, FP671120.3, RAD1, RN7SL856P, and the 10 most significantly downregulated genes were RNA5SP364, GRIN2A, UNC5D, DMBT1P1, MIR3976HG, LINC02199, LINC02822, PANTR1, AC012178.1, CTNNA2. Ingenuity Pathway Analysis identified the 7 most likely to play an important role in the etiology of NOWS. Our study expands insights into the genetic mechanisms of NOWS development.

Hypo-glycosylated hFSH drives ovarian follicular development more efficiently than fully-glycosylated hFSH: enhanced transcription and PI3K and MAPK signaling.

STUDY QUESTION: Does hypo-glycosylated human recombinant FSH (hFSH18/21) have greater in vivo bioactivity that drives follicle development in vivo compared to fully-glycosylated human recombinant FSH (hFSH24)? SUMMARY ANSWER: Compared with fully-glycosylated hFSH, hypo-glycosylated hFSH has greater bioactivity, enabling greater follicular health and growth in vivo, with enhanced transcriptional activity, greater activation of receptor tyrosine kinases (RTKs) and elevated phosphatidylinositol 3-kinase (PI3K)/protein kinase B (AKT) and Mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK) signaling. WHAT IS KNOWN ALREADY: Glycosylation of FSH is necessary for FSH to effectively activate the FSH receptor (FSHR) and promote preantral follicular growth and formation of antral follicles. In vitro studies demonstrate that compared to fully-glycosylated recombinant human FSH, hypo-glycosylated FSH has greater activity in receptor binding studies, and more effectively stimulates the PKA pathway and steroidogenesis in human granulosa cells. STUDY DESIGN, SIZE, DURATION: This is a cross-sectional study evaluating the actions of purified recombinant human FSH glycoforms on parameters of follicular development, gene expression and cell signaling in immature postnatal day (PND) 17 female CD-1 mice. To stimulate follicle development in vivo, PND 17 female CD-1 mice (n = 8-10/group) were treated with PBS (150 µl), hFSH18/21 (1 µg/150 µl PBS) or hFSH24 (1 µg/150 µl PBS) by intraperitoneal injection (i.p.) twice daily (8:00 a.m. and 6:00 p.m.) for 2 days. Follicle numbers, serum anti-Müllerian hormone (AMH) and estradiol levels, and follicle health were quantified. PND 17 female CD-1 mice were also treated acutely (2 h) in vivo with PBS, hFSH18/21 (1 µg) or hFSH24 (1 µg) (n = 3-4/group). One ovary from each mouse was processed for RNA sequencing analysis and the other ovary processed for signal transduction analysis. An in vitro ovary culture system was used to confirm the relative signaling pathways. PARTICIPANTS/MATERIALS, SETTING, METHODS: The purity of different recombinant hFSH glycoforms was analyzed using an automated western blot system. Follicle numbers were determined by counting serial sections of the mouse ovary. Real-time quantitative RT-PCR, western blot and immunofluorescence staining were used to determine growth and apoptosis markers related with follicle health. RNA sequencing and bioinformatics were used to identify pathways and processes associated with gene expression profiles induced by acute FSH glycoform treatment. Analysis of RTKs was used to determine potential FSH downstream signaling pathways in vivo. Western blot and in vitro ovarian culture system were used to validate the relative signaling pathways. MAIN RESULTS AND THE ROLE OF CHANCE: Our present study shows that both hypo- and fully-glycosylated recombinant human FSH can drive follicular growth in vivo. However, hFSH18/21 promoted development of significantly more large antral follicles compared to hFSH24 (P < 0.01). In addition, compared with hFSH24, hFSH18/21 also promoted greater indices of follicular health, as defined by lower BAX/BCL2 ratios and reduced cleaved Caspase 3. Following acute in vivo treatment with FSH glycoforms RNA-sequencing data revealed that both FSH glycoforms rapidly induced ovarian transcription in vivo, but hypo-glycosylated FSH more robustly stimulated Gαs and cAMP-mediated signaling and members of the AP-1 transcription factor complex. Moreover, hFSH18/21 treatment induced significantly greater activation of RTKs, PI3K/AKT and MAPK/ERK signaling compared to hFSH24. FSH-induced indices of follicle growth in vitro were blocked by inhibition of PI3K and MAPK. LARGE SCALE DATA: RNA sequencing of mouse ovaries. Data will be shared upon reasonable request to the corresponding author. LIMITATIONS, REASONS FOR CAUTION: The observations that hFSH glycoforms have different bioactivities in the present study employing a mouse model of follicle development should be verified in nonhuman primates. The gene expression studies reflect transcriptomes of whole ovaries. WIDER IMPLICATIONS OF THE FINDINGS: Commercially prepared recombinant human FSH used for ovarian stimulation in human ART is fully-glycosylated FSH. Our findings that hypo-glycosylated hFSH has greater bioactivity enabling greater follicular health and growth without exaggerated estradiol production in vivo, demonstrate the potential for its development for application in human ART. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by NIH 1P01 AG029531, NIH 1R01 HD 092263, VA I01 BX004272, and the Olson Center for Women's Health. JSD is the recipient of a VA Senior Research Career Scientist Award (1IK6 BX005797). This work was also partially supported by National Natural Science Foundation of China (No. 31872352). The authors declared there are no conflicts of interest.

KSR1- and ERK-dependent translational regulation of the epithelial-to-mesenchymal transition.

The epithelial-to-mesenchymal transition (EMT) is considered a transcriptional process that induces a switch in cells from a polarized state to a migratory phenotype. Here, we show that KSR1 and ERK promote EMT-like phenotype through the preferential translation of Epithelial-Stromal Interaction 1 (EPSTI1), which is required to induce the switch from E- to N-cadherin and coordinate migratory and invasive behavior. EPSTI1 is overexpressed in human colorectal cancer (CRC) cells. Disruption of KSR1 or EPSTI1 significantly impairs cell migration and invasion in vitro, and reverses EMT-like phenotype, in part, by decreasing the expression of N-cadherin and the transcriptional repressors of E-cadherin expression, ZEB1 and Slug. In CRC cells lacking KSR1, ectopic EPSTI1 expression restored the E- to N-cadherin switch, migration, invasion, and anchorage-independent growth. KSR1-dependent induction of EMT-like phenotype via selective translation of mRNAs reveals its underappreciated role in remodeling the translational landscape of CRC cells to promote their migratory and invasive behavior.

The majority of cancer deaths result from tumor cells spreading to other parts of the body via a process known as metastasis. 90% of all cancers originate in epithelial cells that line the inner and outer surface of organs in our bodies. Epithelial cells, however, are typically stationary and must undergo various chemical and physical changes to transform in to migratory cells that can invade other tissues. This transformation process alters the amount of protein cells use to interact with one another. For example, epithelial cells from the colon produce less of a protein called E-cadherin as they transition into migrating cancer cells and make another protein called N-cadherin instead. A protein called KSR1 is a key component of a signaling pathway that is responsible for generating the proteins colon cancer cells need to survive. But it is unknown which proteins KSR1 helps synthesize and whether it plays a role in the metastasis of colon cancer cells. To investigate this, Rao et al. studied the proteins generated by cancerous colon cells cultured in the laboratory, in the presence and absence of KSR1. The experiment showed that KSR1 increases the levels of a protein called EPSTI1, which colon cancer cells need to transform into migratory cells. Depleting KSR1 caused cancer cells to generate less EPSTI1 and to share more features with healthy cells, such as higher levels of E-cadherin on their surface and reduced mobility. Adding EPSTI1 to the cancer cells that lacked KSR1 restored the traits associated with metastasis, such as high levels of N-cadherin, and allowed the cells to move more easily. These findings suggest that KSR1 and EPSTI1 could be new drug targets for reducing, or potentially reversing, the invasive behavior of colon cancer cells. However, further investigation is needed to reveal how EPSTI1 is generated and how this protein helps colon cancer cells move and invade other tissues.

Pan-Cancer Analysis of Human Kinome Gene Expression and Promoter DNA Methylation Identifies Dark Kinase Biomarkers in Multiple Cancers.

Kinases are a group of intracellular signaling molecules that play critical roles in various biological processes. Even though kinases comprise one of the most well-known therapeutic targets, many have been understudied and therefore warrant further investigation. DNA methylation is one of the key epigenetic regulators that modulate gene expression. In this study, the human kinome's DNA methylation and gene expression patterns were analyzed using the level-3 TCGA data for 32 cancers. Unsupervised clustering based on kinome data revealed the grouping of cancers based on their organ level and tissue type. We further observed significant differences in overall kinase methylation levels (hyper- and hypomethylation) between the tumor and adjacent normal samples from the same tissue. Methylation expression quantitative trait loci (meQTL) analysis using kinase gene expression with the corresponding methylated probes revealed a highly significant and mostly negative association (~92%) within 1.5 kb from the transcription start site (TSS). Several understudied (dark) kinases (PKMYT1, PNCK, BRSK2, ERN2, STK31, STK32A, and MAPK4) were also identified with a significant role in patient survival. This study leverages results from multi-omics data to identify potential kinase markers of prognostic and diagnostic importance and further our understanding of kinases in cancer.

Identification of Prognostic Markers in Cholangiocarcinoma Using Altered DNA Methylation and Gene Expression Profiles.

BACKGROUND: Cholangiocarcinoma (CCA) is a rare disease, but it is amongst the most lethal cancers with a median survival under 1 year. Variations in DNA methylation and gene expression have been extensively studied in other cancers for their role in pathogenesis and disease prognosis, but these studies are very limited in CCA. This study focusses on the identification of DNA methylation and gene expression prognostic biomarkers using multi-omics data of CCA tumors from The Cancer Genome Atlas (TCGA). METHOD: We have conducted a genome-wide analysis of differential DNA methylation and gene/miRNA expression using data from 36 CCA tumor and 9 normal samples from TCGA. The impact of DNA methylation in promoters and long-range distal enhancers on the regulation and expression of CCA-associated genes was examined using linear regression. Next, we conducted network analyses on genes which are regulated by DNA methylation as well as by miRNA. Finally, we performed Kaplan-Meier and Cox proportional hazards regression analyses in order to identify the role of selected methylation sites and specific genes and miRNAs in patient survival. We also performed real-time quantitative PCR (qPCR) to confirm the change in gene expression in CCA patients' tumor and adjacent normal samples. RESULTS: Altered DNA methylation was observed on 12,259 CpGs across all chromosomes, of which 78% were hypermethylated. We observed a strong negative relationship between promoter hypermethylation and corresponding gene expression in 92% of the CpGs. Differential expression analyses revealed altered expression patterns in 3,305 genes and 101 miRNAs. Finally, we identified 17 differentially methylated promoter CpGs, 72 differentially expressed genes, and two miRNAs that are likely associated with patient survival. Pathway analysis suggested that cell division, bile secretion, amino acid metabolism, PPAR signaling, hippo signaling were highly affected by gene expression and DNA methylation alterations. The qPCR analysis further confirmed that MDK, HNF1B, PACS1, and GLUD1 are differentially expressed in CCA. CONCLUSION: Based on the survival analysis, we conclude that DEPDC1, FUT4, MDK, PACS1, PIWIL4 genes, miR-22, miR-551b microRNAs, and cg27362525 and cg26597242 CpGs can strongly support their use as prognostic markers of CCA.

Proximal tubule cyclophilin D regulates fatty acid oxidation in cisplatin-induced acute kidney injury.

Regardless of the etiology, acute kidney injury involves aspects of mitochondrial dysfunction and ATP depletion. Fatty acid oxidation is the preferred energy source of the kidney and is inhibited during acute kidney injury. A pivotal role for the mitochondrial matrix protein, cyclophilin D in regulating overall cell metabolism is being unraveled. We hypothesize that mitochondrial interaction of proximal tubule cyclophilin D and the transcription factor PPARα modulate fatty acid beta-oxidation in cisplatin-induced acute kidney injury. Cisplatin injury resulted in histological and functional damage in the kidney with downregulation of fatty acid oxidation genes and increase of intrarenal lipid accumulation. However, proximal tubule-specific deletion of cyclophilin D protected the kidneys from the aforementioned effects. Mitochondrial translocation of PPARα, its binding to cyclophilin D, and sequestration led to inhibition of its nuclear translocation and transcription of PPARα-regulated fatty acid oxidation genes during cisplatin-induced acute kidney injury. Genetic or pharmacological inhibition of cyclophilin D preserved nuclear expression and transcriptional activity of PPARα and prevented the impairment of fatty acid oxidation and intracellular lipid accumulation. Docking analysis identified potential binding sites between PPARα and cyclophilin D. Thus, our results indicate that proximal tubule cyclophilin D elicits impaired mitochondrial fatty acid oxidation via mitochondrial interaction between cyclophilin D and PPARα. Hence, targeting their interaction may be a potential therapeutic strategy to prevent energy depletion, lipotoxicity and cell death in cisplatin-induced acute kidney injury.

Survival Analysis of Multi-Omics Data Identifies Potential Prognostic Markers of Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is the most common and among the deadliest of pancreatic cancers. Its 5-year survival is only ∼8%. Pancreatic cancers are a heterogeneous group of diseases, of which PDAC is particularly aggressive. Like many other cancers, PDAC also starts as a pre-invasive precursor lesion (known as pancreatic intraepithelial neoplasia, PanIN), which offers an opportunity for both early detection and early treatment. Even advanced PDAC can benefit from prognostic biomarkers. However, reliable biomarkers for early diagnosis or those for prognosis of therapy remain an unfulfilled goal for PDAC. In this study, we selected 153 PDAC patients from the TCGA database and used their clinical, DNA methylation, gene expression, and micro-RNA (miRNA) and long non-coding RNA (lncRNA) expression data for multi-omics analysis. Differential methylations at about 12,000 CpG sites were observed in PDAC tumor genomes, with about 61% of them hypermethylated, predominantly in the promoter regions and in CpG-islands. We correlated promoter methylation and gene expression for mRNAs and identified 17 genes that were previously recognized as PDAC biomarkers. Similarly, several genes (B3GNT3, DMBT1, DEPDC1B) and lncRNAs (PVT1, and GATA6-AS) are strongly correlated with survival, which have not been reported in PDAC before. Other genes such as EFR3B, whose biological roles are not well known in mammals are also found to strongly associated with survival. We further identified 406 promoter methylation target loci associated with patients survival, including known esophageal squamous cell carcinoma biomarkers, cg03234186 (ZNF154), and cg02587316, cg18630667, and cg05020604 (ZNF382). Overall, this is one of the first studies that identified survival associated genes using multi-omics data from PDAC patients.