UCSC genome browser tutorial

Mitotic Spindle Positioning (MISP) Facilitates Colorectal Cancer Progression by Forming a Complex with Opa Interacting Protein 5 (OIP5) and Activating the JAK2-STAT3 Signaling Pathway

Patients with inflammatory bowel disease (IBD) who experience long-term chronic inflammation of the colon are at an increased risk of developing colorectal cancer (CRC). Mitotic spindle positioning (MISP), an actin-binding protein, plays a role in mitosis and spindle positioning. MISP is found on the apical membrane of the intestinal mucosa and helps stabilize and elongate microvilli, offering protection against colitis. This study explored the role of MISP in colorectal tumorigenesis using a database, human CRC cells, and a mouse model for colitis-induced colorectal tumors triggered by azoxymethane (AOM)/dextran sodium sulfate (DSS) treatment. We found that MISP was highly expressed in colon cancer patient tissues and that reduced MISP expression inhibited cell proliferation. Notably, MISP-deficient mice showed reduced colon tumor formation in the AOM/DSS-induced colitis model. Furthermore, MISP was found to form a complex with Opa interacting protein 5 (OIP5) in the cytoplasm, influencing the expression of OIP5 in a unidirectional manner. We also observed that MISP increased the levels of phosphorylated STAT3 in the JAK2-STAT3 signaling pathway, which is linked to tumorigenesis. These findings indicate that MISP could be a risk factor for CRC, and targeting MISP might provide insights into the mechanisms of colitis-induced colorectal tumorigenesis.

CD2AP is a potential prognostic biomarker of renal clear cell carcinoma

BACKGROUND

CD2-associated protein (CD2AP) is a podocyte-associated gene and its reduced expression is associated with the development of proteinuria and glomerulosclerosis. However, few studies have focused on the correlation between the expression and prognosis of CD2AP in renal clear cell carcinoma (ccRCC). Therefore, we aimed to assess the regulation of CD2AP expression and prognostic value in ccRCC.

METHODS

Multiple databases were employed to examine the expression of CD2AP in ccRCC. RT-qPCR, Western Blot and immunohistochemistry were used to validate CD2AP expression in different cell lines and tissue samples. Kaplan-Meier analysis and ROC curve analysis were performed on the predictive prognostic performance of CD2AP. COX regression was used to construct CD2AP-related prognostic models. The TIMER and TISIDB databases were used to analyze the correlation of tumor-infiltrating immune cells with gene expression, mutations, somatic copy number variation, and immune molecules. Mass spectrometry was used to detect methylation status of the promoter CpG site of CD2AP in multiple cells.

RESULTS

We found that CD2AP expression was downregulated in ccRCC and its lower expression level was correlation with worse patient prognosis, higher tumor stage and grade and distant metastasis through analysis of databases, ccRCC cell lines and clinical tissue samples. Moreover, database and mass spectrometry techniques identified and validated cg12968598 hypermethylation as one of the key reasons for the downregulation of CD2AP expression. CD2AP expression was also associated with macrophage and neutrophil infiltration.

CONCLUSIONS

Taken together, our results suggest that CD2AP can be used as a diagnostic and prognostic biomarker in ccRCC patients and that DNA hypermethylation plays an important role in reducing CD2AP expression.

Defining Candidate Imprinted loci in Bos taurus

Using a whole-genome assembly of Bos taurus, I applied my bioinformatics strategy to locate candidate imprinting control regions (ICRs) genome-wide. In mammals, genomic imprinting plays essential roles in embryogenesis. In my strategy, peaks in plots mark the locations of known, inferred, and candidate ICRs. Genes in the vicinity of candidate ICRs correspond to potential imprinted genes. By displaying my datasets on the UCSC genome browser, one could view peak positions with respect to genomic landmarks. I give two examples of candidate ICRs in loci that influence spermatogenesis in bulls: CNNM1 and CNR1. I also give examples of candidate ICRs in loci that influence muscle development: SIX1 and BCL6. By examining the ENCODE data reported for mice, I deduced regulatory clues about cattle. I focused on DNase I hypersensitive sites (DHSs). Such sites reveal accessibility of chromatin to regulators of gene expression. For inspection, I chose DHSs in chromatin from mouse embryonic stem cells (ESCs) ES-E14, mesoderm, brain, heart, and skeletal muscle. The ENCODE data revealed that the SIX1 promoter was accessible to the transcription initiation apparatus in mouse ESCs, mesoderm, and skeletal muscles. The data also revealed accessibility of BCL6 locus to regulatory proteins in mouse ESCs and examined tissues.

Prognostic value and immune-infiltration pattern of FOXD3-AS1 in patients with glioma

Gliomas are difficult-to-treat brain tumors due to their aggressive nature, rapid proliferation, and high invasiveness (Zhang et al., J Cell Biochem, 2019, 120 (9), 15106–15118; Ge et al., Int J Biochem Cell Biol, 2021, 139, 106054). FOXD3-AS1 has been identified as an emerging potential target for tumor prediction and treatment in many studies (Qin et al., Front Oncol, 2021, 11, 688027). However, the utility of FOXD3-AS1 has not been reported in glioma patients (Li et al., Cancer Manag Res, 2021, 13, 9037–9048). The differential profiles of FOXD3-AS1 in TCGA–GBMLGG database were analyzed across clinical subgroups. The analysis of overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI) revealed that a high level of FOXD3-AS1 was associated with a poor prognosis and survival outcome. Based on the Cox regression analysis, FOXD3-AS1 was found to be a high-risk factor for glioma that affects prognosis outcomes independently. More importantly, because oxidative stress is closely linked to glioma prognosis, we focused on the potential mechanisms of six oxidative stress co-expressed genes with FOXD3-AS1. In addition, the predictive value of FOXD3-AS1 was determined for each clinical subgroup status. The ROC curve results showed that FOXD3-AS1 had a good predictive performance. A stratified clinicopathological subgroup analysis revealed that high expression of FOXD3-AS1 is associated with a poor prognosis. This also indicates a link between FOXD3-AS1 and tumorigenesis and prognosis, which has potential application value. Furthermore, the immune cell infiltration of FOXD3-AS1 and the signal marker correlation suggested that immune cell infiltration differed significantly between immune cell subsets. To the best of our knowledge, this is the first report to investigate FOXD3-AS1 in glioma and how it may modulate GBM and LGG immune microenvironments. Furthermore, FOXD3-AS1 was detected in tumor and paraneoplastic tissues using RT–qPCR. Transwell analysis verified the migration and invasion of the FOXD3-AS1 knockout group in vitro to a certain extent. In conclusion, FOXD3-AS1 can be used as a prognostic indicator for GBM and LGG, and it is closely related to immune infiltration and response to oxidative stress, which may contribute to the advancement of glioma immunotherapy research.

Pathogenicity of PKCγ Genetic Variants-Possible Function as a Non-Invasive Diagnostic Biomarker in Ovarian Cancer

Ovarian cancer has the highest mortality rate among gynecologic malignancies, owing to its misdiagnosis or late diagnosis. Identification of its genetic determinants could improve disease outcomes. Conventional Protein Kinase C-γ (PKCγ) dysregulation is reported in several cancers. Similarly, its variant rs1331262028 is also reported to have an association with hepatocellular carcinoma. Therefore, the aim of the present study was to analyze the variant rs1331262028 association with ovarian cancer and to determine its impact on PKCγ's protein interactions. Association of variation was determined through genotyping PCR (cohort size:100). Protein-protein docking and molecular dynamic simulation were carried out to study the variant impact of PKCγ interactions. The study outcome indicated the positive association of variant rs1331262028 with ovarian cancer and its clinicopathological features. Molecular dynamics simulation depicted the potential influence of variation on PKCγ molecular signaling. Hence, this study provided the foundations for assessing variant rs1331262028 as a potential prognostic marker for ovarian cancer. Through further validation, it can be applied at the clinical level.

Methylation-Mediated Silencing of RBP7 Promotes Breast Cancer Progression through PPAR and PI3K/AKT Pathway

Retinoid-binding protein7 (RBP7) is a member of the cellular retinol-binding protein (CRBP) family, which is involved in the pathogenesis of breast cancer. The study aims to illustrate the prognostic value and the potential regulatory mechanisms of RBP7 expression in breast cancer. Bioinformatics analysis with the TCGA and CPTAC databases revealed that the mRNA and protein expression levels of RBP7 in normal were higher compared to breast cancer tissues. Survival analysis displayed that the lower expression of RBP7, the worse the prognosis in ER-positive (ER+) breast cancer patients. Genomic analysis showed that low expression of RBP7 correlates with its promoter hypermethylation in breast cancer. Functional enrichment analysis demonstrated that downregulation of RBP7 expression may exert its biological influence on breast cancer through the PPAR pathway and the PI3K/AKT pathway. In summary, we identified RBP7 as a novel biomarker that is helpful for the prognosis of ER+ breast cancer patients. Promoter methylation of RBP7 is involved in its gene silencing in breast cancer, thus regulating the occurrence and development of ER+ breast cancer through the PPAR and PI3K/AKT pathways.

Dynamic Expression of Palmitoylation Regulators across Human Organ Development and Cancers Based on Bioinformatics

Protein palmitoylation is a reversible modification process that links palmitate to cysteine residues via a reversible thioester bond. Palmitoylation exerts an important role in human organ development and tumor progression. However, a comprehensive landscape regarding the dynamic expression of palmitoylation regulators in human organ development remains unclear. In this study, we analyzed the dynamic expression of palmitoylation regulators in seven organ development and eight cancer types based on bioinformatics. We found that the expression levels of most palmitoylation regulators were altered after birth. In particular, ZDHHC7/20/21 exhibited converse expression patterns in multiple cancer types. Survival analysis showed that the poor prognosis in patients with kidney renal clear carcinoma (KIRC) is related to low expression of ZDHHC7/20/21, and a high expression of ZDHHC7/20/21 is related to worse survival in patients with liver hepatocellular carcinoma (LIHC). Furthermore, we found that the expression of ZDHHC7 is associated with infiltration levels of some types of immune cells in the tumor microenvironment (TME), and we explored the relationship between ZDHHC7 expression and immune checkpoint (ICP) genes across 33 cancer types. In addition, gene set enrichment analysis (GSEA) results indicated that ZDHHC7 might regulate different genes to mediate the same pathway in different organs. In summary, the comprehensive analysis of palmitoylation regulators reveals their functions in human organ development and cancer, which may provide new insights for developing new tumor markers.

Along the Bos taurus genome, uncover candidate imprinting control regions

Background

In mammals, Imprinting Control Regions (ICRs) regulate a subset of genes in a parent-of-origin-specific manner. In both human and mouse, previous studies identified a set of CpG-rich motifs occurring as clusters in ICRs and germline Differentially Methylated Regions (gDMRs). These motifs consist of the ZFP57 binding site (ZFBS) overlapping a subset of MLL binding units known as MLL morphemes. MLL or MLL1 (Mixed Lineage Leukemia 1) is a relatively large multidomain protein that plays a central role in the regulation of transcription. The structures of both MLL1 and MLL2 include a domain (MT) that binds CpG-rich DNA and a conserved domain (SET) that methylates lysine 4 in histone H3 producing H3K4me3 marks in chromatin.

Results

Since genomic imprinting impacts many developmental and key physiological processes, we followed a previous bioinformatics strategy to pinpoint ICR positions in the Bos taurus genome. Initial genome-wide analyses involved finding the positions of ZFP57 binding sites, and the CpG-rich motifs (ZFBS-morph overlaps) along cattle chromosomal DNA. By creating plots displaying the density of ZFBS-morph overlaps, we removed background noise and thus improved signal detection. With the density-plots, we could view the positions of peaks locating known and candidate ICRs in cattle DNA. Our evaluations revealed the correspondence of peaks in plots to reported known and inferred ICRs/DMRs in cattle. Beside peaks pinpointing such ICRs, the density-plots also revealed additional peaks. Since evaluations validated the robustness of our approach, we inferred that the additional peaks may correspond to candidate ICRs for imprinted gene expression.

Conclusion

Our bioinformatics strategy offers the first genome-wide approach for systematically localizing candidate ICRs. Furthermore, we have tailored our datasets for upload onto the UCSC genome browser so that researchers could find known and candidate ICRs with respect to a wide variety of annotations at all scales: from the positions of Single Nucleotide Polymorphisms (SNPs), to positions of genes, transcripts, and repeated DNA elements. Furthermore, the UCSC genome browser offers tools to produce enlarged views: to uncover the genes in the vicinity of candidate ICRs and thus discover potential imprinted genes for experimental validations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08694-3.

GATA4 and estrogen receptor alpha bind at SNPs rs9921222 and rs10794639 to regulate AXIN1 expression in osteoblasts

Osteoporosis is a serious public health problem that affects 200 million people worldwide. Genome-wide association studies have revealed the association between several single nucleotide polymorphisms (SNPs) near WNT/β-catenin signaling genes and bone mineral density (BMD). The activation of β-catenin by WNT ligands is required for osteoblast differentiation. SNP rs9921222 is an intronic variant of AXIN1 (a scaffold protein in the destruction complex that regulates β-catenin signaling) that correlates with BMD. However, the biological mechanism of SNP rs9921222 has never been reported. Here, we show that the genotype of SNP rs9921222 correlates with the expression of AXIN1 in human osteoblasts. RNA and genomic DNA were analyzed from primary osteoblasts from 111 different individuals. Homozygous TT at rs9921222 correlates with a higher expression of AXIN1 than homozygous CC. Regional association analysis showed that rs9921222 is in high linkage disequilibrium (LD) with SNP rs10794639. In silico transcription factor analysis predicted that rs9921222 is within a GATA4 motif and rs10794639 is adjacent to an estrogen receptor alpha (ERα) motif. Mechanistically, GATA4 and ERα bind at SNPs rs9921222 and rs10794639 as detected by ChIP-qPCR. Luciferase assays demonstrate that rs9921222 is the causal SNP to alter ERα and GATA4 binding. GATA4 promoted the expression, and in contrast, ERα suppressed the expression of AXIN1 via the histone deacetylase complex member SIN3A. Functionally, the level of AXIN1 negatively correlates with the level of transcriptionally active β-catenin. In summary, we have discovered a molecular mechanism of the SNP rs9921222 to regulate AXIN1 through GATA4 and ERα binding in human osteoblasts.

Association Study of TAF1 Variants in Parkinson’s Disease

Increasing evidence reveals sex as an important factor in the development of Parkinson’s disease (PD), but associations between genes on the sex chromosomes and PD remain unknown. TAF1 is a gene located on the X chromosome which is known to cause X-linked syndromic mental retardation-33 (MRXS33) and X-linked Dystonia-Parkinsonism (XDP). In this study, we conducted whole-exome sequencing (WES) among 1,917 patients with early-onset or familial PD and 1,652 controls in a Chinese population. We detected a hemizygous frameshift variant c.29_53dupGGA(CAG)₂CTACCATCA(CTG)₂C (p.A19Dfs*50) in two unrelated male patients. Further segregation analysis showed an unaffected family member carried this variant, which suggested the penetrance of the variant may be age-related and incomplete. To verify the effects of TAF1 on PD, genetic analyses were carried separately by gender. Analysis of rare variants by optimal sequence kernel association (SKAT-O) test showed a nominally significant difference in variant burden between the male PD patients and controls (2.01 vs. 1.38%, p = 0.027). In the female group, none of the variant types showed significant association with PD in this study. In conclusion, we found rare variants in TAF1 may be implicated in PD, but further genetic and functional analyses were needed.

Full-length transcriptome analysis of Adiantum flabellulatum gametophyte

Ferns are important components of plant communities on earth, but their genomes are generally very large, with many redundant genes, making whole genome sequencing of ferns prohibitively expensive and time-consuming. This means there is a significant lack of fern reference genomes, making molecular biology research difficult. The gametophytes of ferns can survive independently, are responsible for sexual reproduction and the feeding of young sporophytes, and play an important role in the alternation of generations. For this study, we selected Adiantum flabellulatum as it has both ornamental and medicinal value and is also an indicator plant of acidic soil. The full-length transcriptome sequencing of its gametophytes was carried out using PacBio three-generation sequencing technology. A total of 354,228 transcripts were obtained, and 231,705 coding sequences (CDSs) were predicted, including 5,749 transcription factors (TFs), 2,214 transcription regulators (TRs) and 4,950 protein kinases (PKs). The transcripts annotated by non-redundant protein sequence database (NR), Kyoto encyclopedia of genes and genomes (KEGG), eukaryotic ortholog groups (KOG), Swissprot, protein family (Pfma), nucleotide sequence database (NT) and gene ontology (GO) were 251,501, 197,474, 193,630, 194,639, 195,956, 113,069 and 197,883, respectively. In addition, 138,995 simple sequence repeats (SSRs) and 111,793 long non-coding RNAs (lncRNAs) were obtained. We selected nine chlorophyll synthase genes for qRT-PCR, and the results showed that the full-length transcript sequences and the annotation information were reliable. This study can provide a reference gene set for subsequent gene expression quantification.

Four drug metabolism-related subgroups of pancreatic adenocarcinoma in prognosis, immune infiltration, and gene mutation

We aimed to screen the drug metabolism-related subgroups of pancreatic adenocarcinoma (PAAD) and to study the prognosis, clinical features, immune infiltration, and gene mutation differences of different subtypes in PAAD patients. All 181 cases of PAAD samples and clinical characteristics data were downloaded from The Cancer Genome Atlas (TCGA). After matching the drug metabolism-related genes downloaded from PMID 33202946 with the TCGA dataset, the drug metabolism-related genes were initially obtained. Besides, univariate Cox regression analysis was used to screen the drug metabolism genes related to the prognosis of PAAD. Moreover, the construction of the protein–protein interaction (PPI) network and gene ontology were performed. The four subgroups of PAAD obtained from unsupervised clustering analysis were systematically analyzed, including prognostic, GSVA, immune infiltration, and gene mutation analysis. A total of 83 drug metabolism genes related to the prognosis of PAAD were obtained and enriched in 16 pathways. The PPI network was composed of 248 relationship pairs. Four subgroups that can identify different subtypes of PPAD were obtained, and there were significant differences in survival and clinical characteristics, mutation types, and immune infiltration abundance between subgroups. A total of 17 different pathways among the four subgroups involved in cell cycle, response to stimulants such as drugs, and transmembrane transport. In this study, the four subgroups related to the drug metabolism of PAAD were comprehensively analyzed, and the important role of drug metabolism-related genes in the immune infiltration and prognosis of PAAD were emphasized.

Pan-Cancer Analysis Predicts FOXS1 as a Key Target in Prognosis and Tumor Immunotherapy

Purpose

Only a few studies have reported the role of FOXS1, a transcriptional factor, in the tumor development process. In this article, we investigate the function of FOXS1 in distinct neoplastic development and the tumor immune microenvironment (TIME).

Patients and Methods

The latent roles of FOXS1 in various tumors were prospected based on TCGA, GTEx, CCLE, GEPIA2, cBioPortal, TIMER, ImmuCellAI databases, GSVA datasets, GSEA datasets, and R packages. The expression difference, gene alteration, clinical characteristics, prognostic values, biological mechanism, potential pathways, tumor microenvironment, and immune cell infiltration related to FOXS1 were appraised.

Results

FOXS1 was strongly expressed in pan-cancer, and this gene was associated with low survival rates. FOXS1 was linked to many pathways that are cancer-promoting and immune-related. The expression of this transcriptional factor in cancers was positively related to immune cell infiltration, especially M2-like macrophages and Treg cells. In addition to that, FOXS1 demonstrated a positive relationship with many immune-suppression genes, such as TGFB1 and ARORA2A.

Conclusion

Our study identified an oncogenic effect of FOXS1, which may play a vital role as a prognosticative biological marker in pan-cancer. Exorbitant expression of FOXS1 is associated with high TAMs and Treg cells infiltration. These cells have an immunosuppressive function and promote the development of the immunosuppressive tumor microenvironment. The research of FOXS1 provided a potential drug target for tumor immunotherapy.

Estrogen receptor alpha and NFATc1 bind to a bone mineral density-associated SNP to repress WNT5B in osteoblasts

Genetic factors and estrogen deficiency contribute to the development of osteoporosis. The single-nucleotide polymorphism (SNP) rs2887571 is predicted from genome-wide association studies (GWASs) to associate with osteoporosis but has had an unknown mechanism. Analysis of osteoblasts from 110 different individuals who underwent joint replacement revealed that the genotype of rs2887571 correlates with WNT5B expression. Analysis of our ChIP-sequencing data revealed that SNP rs2887571 overlaps with an estrogen receptor alpha (ERα) binding site. Here we show that 17β-estradiol (E2) suppresses WNT5B expression and further demonstrate the mechanism of ERα binding at the enhancer containing rs2887571 to suppress WNT5B expression differentially in each genotype. ERα interacts with NFATc1, which is predicted to bind directly at rs2887571. CRISPR-Cas9 and ChIP-qPCR experiments confirm differential regulation of WNT5B between each allele. Homozygous GG has a higher binding affinity for ERα than homozygous AA and results in greater suppression of WNT5B expression. Functionally, WNT5B represses alkaline phosphatase expression and activity, decreasing osteoblast differentiation and mineralization. Furthermore, WNT5B increases interleukin-6 expression and suppresses E2-induced expression of alkaline phosphatase during osteoblast differentiation. We show that WNT5B suppresses the differentiation of osteoblasts via receptor tyrosine kinase-like orphan receptor 1/2 (ROR1/2), which activates DVL2/3/RAC1/CDC42/JNK/SIN3A signaling and inhibits β-catenin activity. Together, our data provide mechanistic insight into how ERα and NFATc1 regulate the non-coding SNP rs2887571, as well as the function of WNT5B on osteoblasts, which could provide alternative therapeutic targets for osteoporosis.

Bioinformatics Tools to Understand Notch

As a result of the culmination of data, and the fast-paced advancement of new research, all the biological information collected can make it difficult to sort data. This is oftentimes experienced when learning about the human genome. Fortunately, with the advancement of technology, the field of bioinformatics has emerged which has allowed for the creation of a variety of biological databases. These biological databases provide a condensed reservoir of organized information that is easy to use and topic-specific. Here, we provide a list of 39 biological databases that help break down the fundamental details of a gene. This chapter uses the NOTCH1 gene as an example to demonstrate how biological databases can be used to extract gene information. Five sections were created to highlight the major areas needed to build a comprehensive foundation of NOTCH1. The first section lists databases containing basic gene and protein product information. The next section consists of protein interactions and signaling pathway databases which are essential in understanding the biological processes a gene product is involved in. Gene expression and disease databases are the next two sections which are connected since disease results from the aberrant expression of a gene product. The last database section examines model organisms which serve a key role in the study of human genetic diseases. Using these databases, we can elucidate NOTCH1's gene/protein structure, expression, and vital physiological function through the Notch signaling pathway.

Coupling of serum CK20 and hyper-methylated CLIP4 as promising biomarker for colorectal cancer diagnosis: from bioinformatics screening to clinical validation

Colorectal cancer (CRC) is one of the most common and lethal malignancies. The identification of minimally invasive and precise biomarkers is an urgent need for the early diagnosis of CRC. Through bioinformatics analysis of 395 CRC tissues and 63 CRC cell lines, CK18, CK20, de-methylated HPDL and hyper-methylated CLIP4 were identified as candidate serum biomarkers. Then, a training cohort consisting of 60 CRC, 30 colorectal adenomas (CA) and 33 healthy controls and a validation cohort consisting of 60 CRC, 30 CA and 30 healthy controls were enrolled. In the training cohort, enzyme-linked immunosorbent assay (ELISA) showed that CK18 and CK20 were all significantly higher in CRC and CA. CK18 diagnosed CRC with 46.67% sensitivity and 87.3% specificity; CK20 diagnosed CRC with 28.33% sensitivity and 90.47% specificity. Methylation-specific PCR (MSP) indicated that de-methylated HPDL and hyper-methylated CLIP4 were significantly detected in CRC and CA. De-methylated HPDL diagnosed CRC with 36.67% sensitivity and 93.65% specificity and hyper-methylated CLIP4 with 73.33% sensitivity and 84.13% specificity. Random combined analysis suggested that CK20/hyper-methylated CLIP4 diagnosed CRC with 91.67% sensitivity and 82.54% specificity. In the validation cohort, CK20 diagnosed CRC with 36.7% sensitivity and 88.3% specificity and hyper-methylated CLIP4 with 80% sensitivity and 85% specificity. CK20/hyper-methylated CLIP4 diagnosed CRC with 95% sensitivity and 81.7% specificity. Compared with serum biomarkers reported before, CK20/hyper-methylated CLIP4 possessed the potential to be a new effective and precise diagnostic biomarker for CRC.

Virtual Gene Concept and a Corresponding Pragmatic Research Program in Genetical Data Science

Mendel proposed an experimentally verifiable paradigm of particle-based heredity that has been influential for over 150 years. The historical arguments have been reflected in the near past as Mendel's concept has been diversified by new types of omics data. As an effect of the accumulation of omics data, a virtual gene concept forms, giving rise to genetical data science. The concept integrates genetical, functional, and molecular features of the Mendelian paradigm. I argue that the virtual gene concept should be deployed pragmatically. Indeed, the concept has already inspired a practical research program related to systems genetics. The program includes questions about functionality of structural and categorical gene variants, about regulation of gene expression, and about roles of epigenetic modifications. The methodology of the program includes bioinformatics, machine learning, and deep learning. Education, funding, careers, standards, benchmarks, and tools to monitor research progress should be provided to support the research program.

cyPhyRNA-seq: a genome-scale RNA-seq method to detect active self-cleaving ribozymes by capturing RNAs with 2',3' cyclic phosphates and 5' hydroxyl ends

Self-cleaving ribozymes are catalytically active RNAs that cleave themselves into a 5′-fragment with a 2′,3′-cyclic phosphate and a 3′-fragment with a 5′-hydroxyl. They are widely applied for the construction of synthetic RNA devices and RNA-based therapeutics. However, the targeted discovery of self-cleaving ribozymes remains a major challenge. We developed a transcriptome-wide method, called cyPhyRNA-seq, to screen for ribozyme cleavage fragments in total RNA extract. This approach employs the specific ligation-based capture of ribozyme 5′-fragments using a variant of the Arabidopsis thaliana tRNA ligase we engineered. To capture ribozyme 3′-fragments, they are enriched from total RNA by enzymatic treatments. We optimized and enhanced the individual steps of cyPhyRNA-seq in vitro and in spike-in experiments. Then, we applied cyPhyRNA-seq to total RNA isolated from the bacterium Desulfovibrio vulgaris and detected self-cleavage of the three predicted type II hammerhead ribozymes, whose activity had not been examined to date. cyPhyRNA-seq can be used for the global analysis of active self-cleaving ribozymes with the advantage to capture both ribozyme cleavage fragments from total RNA. Especially in organisms harbouring many self-cleaving RNAs, cyPhyRNA-seq facilitates the investigation of cleavage activity. Moreover, this method has the potential to be used to discover novel self-cleaving ribozymes in different organisms.

Low-frequency and rare coding variants of NUS1 contribute to susceptibility and phenotype of Parkinson's disease

NUS1 has been recently identified as a candidate gene for Parkinson's disease (PD). Few studies have examined the association of NUS1 variants with PD susceptibility and phenotypes. In the first cohort, whole-exome sequencing was performed to identify variants in NUS1 exon-coding and exon-intron regions in 1542 cases and 1625 controls. 13 variants were totally detected, of which 10 rare variants and 3 low-frequency variants. Burden analysis showed that rare NUS1 variants significantly enriched in PD (p=0.016). We also performed a meta-analysis based on previous and our studies to correlate NUS1 mutations with PD susceptibility. Integrating our previous cohort (3210 cases and 2807 controls) and the first cohort identified the significant association of rs539668656 with PD risk (odds ratio (OR) = 2.82, p = 0.016). The genotype-phenotype association analysis showed that patients carrying rare variants, or rs539668656 were significantly associated with earlier onset age, depression, emotional impairment and severe disease condition. Our results support the role of NUS1 rare variants and rs539668656 towards PD susceptibility and phenotype.

Identification and characteristic analysis of enhancers across 13 major cancer types

Enhancers are often mutated and dysregulated in various diseases such as cancer. By integrating the function annotation of the mammalian genome (FANTOM) enhancers expression profiles and RNA-seq data from The Cancer Genome Atlas (TCGA) of 13 cancers and their corresponding para-cancerous tissues, we systematically identified a total of 4702 significantly differentially expressed (DE) enhancers. Furthermore, a total of 1036 DE genes regulated by DE enhancers were identified. It was found that in these 13 cancers, most (61.13%) enhancers were ubiquitously expressed, whereas DE enhancers were more likely to be tissue-specific expressed, and the DE genes regulated by DE enhancers were significantly enriched in cancer-related pathways. Finally, it was manifested that 74 single nucleotide polymorphisms (SNPs) were located in 37 DE enhancers, and these SNPs affected the gain and loss of functional transcription factor binding sites of 758 transcription factors, which were shown to be highly correlated with tumorigenesis and development.

Analysis of the Prognostic Value and Gene Expression Mechanism of SHOX2 in Lung Adenocarcinoma

Background: Detection of SHOX2 methylation has been used to assist in the early diagnosis of lung cancer in many hospitals as SHOX2 may be important in the tumorigenesis of lung cancer. However, there are few studies on the mRNA expression, methylation, and molecular mechanism of SHOX2 in lung cancer. We aimed to explore the role of SHOX2 in lung adenocarcinoma (LUAD).

Methods: First, we examined the differential expression of SHOX2 mRNA and methylation in cancerous and normal tissues using databases. Second, we analyzed the relationship between SHOX2 expression and common clinical parameters in LUAD patients. Third, we further explored the methylated level and its specific location of SHOX2 and the mainly factors of SHOX2 gene expression. Finally, we screened the correlatively expressed genes to analyze the pathways from the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes using DAVID.

Results: We found that the mRNA expression of SHOX2 was higher in multiple cancers, including LUAD and lung squamous cell carcinoma (LUSC), than in normal tissues. Among LUAD patients, SHOX2 expression was higher in patients of middle–young age, with smoking history, in advanced stages, and with nodal distant metastasis. In addition, our results showed that patients with high expression of SHOX2 are prone to recurrence, poor differentiation, and poor prognosis. Thus, we identified that SHOX2 might be an oncogene for LUAD progression. The main factor influencing the high expression of SHOX2 mRNA may be DNA methylation, followed by copy number variation (CNV), but not by gene mutations in LUAD. Unexpectedly, we found that SHOX2 undergoes hypomethylation in the gene body instead of hypermethylation in the promoter. Additionally, SHOX2 has cross talk in the PI3K–Akt signaling pathway and ECM–receptor interaction.

Conclusion:SHOX2 is highly expressed in most cancers. SHOX2 gene expression might be mainly regulated by methylation of its gene body in LUAD, and its high expression or hypomethylation indicates poor differentiation and poor prognosis. SHOX2 could be involved in PI3K–Akt and other important cancer-related signaling pathways to promote tumorigenesis.

Identification of a five-gene signature of the RGS gene family with prognostic value in ovarian cancer

The RGS (regulator of G protein signaling) gene family, which includes negative regulators of G protein-coupled receptors, comprises important drug targets for malignant tumors. It is thus of great significance to explore the value of RGS family genes for diagnostic and prognostic prediction in ovarian cancer. The RNA-seq, immunophenotype, and stem cell index data of pan-cancer, The Cancer Genome Atlas (TCGA) data, and GTEx data of ovarian cancer were downloaded from the UCSC Xena database. In the pan-cancer database, the expression level of RGS1, RGS18, RGS19, and RGS13 was positively correlated with stromal and immune cell scores. Cancer patients with high RGS18 expression were more sensitive to cyclophosphamide and nelarabine, whereas those with high RGS19 expression were more sensitive to cladribine and nelarabine. The relationship between RGS family gene expression and overall survival (OS) and progression-free survival (PFS) of ovarian cancer patients was analyzed using the KM-plotter database, RGS17, RGS16, RGS1, and RGS8 could be used as diagnostic biomarkers of the immune subtype of ovarian cancer, and RGS10 and RGS16 could be used as biomarkers to predict the clinical stage of this disease. Further, Lasso cox analysis identified a five-gene risk score (RGS11, RGS10, RGS13, RGS4, and RGS3). Multivariate COX analysis showed that the risk score was an independent prognostic factor for patients with ovarian cancer. Immunohistochemistry and the HPA protein database confirmed that the five-gene signature is overexpressed in ovarian cancer. GSEA showed that it is mainly involved in the ECM-receptor interaction, TGF-beta signaling pathway, Wnt signaling pathway, and chemokine signaling pathway, which promote the occurrence and development of ovarian cancer. The prediction model of ovarian cancer constructed using RGS family genes is of great significance for clinical decision making and the personalized treatment of patients with ovarian cancer.

Contribution of coding/non-coding variants in NUS1 to late-onset sporadic Parkinson's disease

INTRODUCTION

A recent study reported that rare variants in NUS1 were associated with Parkinson's disease (PD). We aimed to assess the relative contribution of rare and common coding/non-coding variants of NUS1 to late-onset PD patients (LOPD).

METHODS

Whole genome sequencing data were analyzed for target NUS1 regions, derived from a cohort of 1962 cases and 1279 controls. The genetic association analyses were performed using logistic regression analysis and Sequence Kernel association test. Expression quantitative trait loci (eQTL) analysis was conducted to further explore the association of variants with NUS1 expression based on the data from GTEx database.

RESULTS

We identified 18 rare coding variants. p.Y131C was first identified in LOPD. However, no significant burden of rare NUS1 coding variants in LOPD was found. The rare variant sets of two regulatory elements (GH06J117605 and GH06J117674) were significantly enriched in LOPD even after Bonferroni correction (adjusted P = 0.013; adjusted P = 0.010). Considering the joint effect of rare and common variants, all variant sets within GH06J117605 and GH06J117674 showed association with LOPD but were no longer significant after Bonferroni correction. None of the common variants within coding/non-coding regions were significant after Bonferroni correction. The eQTL results suggested these variants in GH06J117605 and GH06J117674 could potentially have eQTL effects on the brain tissues.

CONCLUSIONS

These findings provide novel insight into the role of NUS1 regulatory regions in the development of LOPD and indicate that the variants in regulatory elements of NUS1 may be associated with LOPD by influencing the gene expression level.

Systematic Analysis of the Expression and Prognostic Significance of P4HA1 in Pancreatic Cancer and Construction of a lncRNA-miRNA-P4HA1 Regulatory Axis

Objectives

Prolyl 4-hydroxylase subunit alpha 1 (P4HA1) plays a crucial role in modulating extracellular matrix component and promoting tumor progression by changing tumor adhesion, migration, and other biological behaviors in some cancers. However, its expression pattern, biological function, and underlying mechanism in pancreatic cancer remain largely unclear.

Materials and Methods

In this study, a set of bioinformatics tools were used to analyze the expression of P4HA1 and its prognostic value in pancreatic cancer. In addition, the mechanism through which P4HA1 promotes the progression of pancreatic cancer was explored by constructing a competing endogenous RNA (ceRNA) regulatory axis.

Results

It was found that the mRNA and protein expression of P4HA1 was significantly higher in pancreatic cancer tissues than in normal tissues. Its high P4HA1 expression correlated with poor clinicopathological features (T stage: P = 0.0078; N stage: P = 0.0124; TNM stage: P = 0.0013; pathological grade: P = 0.0108) and poor prognosis [OS: HR = 1, 95% CI (1-1.01), P = 0.00028; DSS: HR = 1, 95% CI (1-1.01), P = 0.00049; PFI: HR = 1.01, 95% CI (1.01-1.02), P = 0.0057; and DFI: HR = 1, 95% CI (1-1.01), P = 0.0034]. The LINC01503/miR-335-5p/P4HA1 axis might mediate the effects of P4HA1 in promoting the progression on pancreatic cancer.

Conclusions

Collectively, our findings suggest that high expression of P4HA1 may be used as a promising prognostic biomarker and could be considered for the development of a novel therapeutic strategy for pancreatic cancer in the future.

Analysis of N6-Methyladenosine Methyltransferase Reveals METTL14 and ZC3H13 as Tumor Suppressor Genes in Breast Cancer

Objectives

Recently, an increasing number of studies have revealed that N6-methyladenosine (m6A) functions as a significant post-transcriptional modification which plays a critical role in the occurrence and progression of enriched tumors by regulating coding and non-coding RNA biogenesis. However, the biological function of m6A in breast cancer remains largely unclear.

Materials and Methods

In this study, we used a series of bioinformatic databases and tools to jointly analyze the expression of m6A methylation transferases (METTL3, METTL14, WTAP, RBM15, RBM15B and ZC3H13) and investigate the prognostic value of METTL14 and ZC3H13 in breast cancer. Besides, we analyzed the downstream carcinogenic molecular mechanisms related to METTL14 and ZC3H13 and their relationship with immune infiltration in breast tumor tissues.

Results

The results showed that METTL14 and ZC3H13 were the down-regulated m6A methylation transferases in breast cancer. Survival outcome analysis suggested that abnormally low expression of METTL14 and ZC3H13 could predict unfavorable prognosis in four breast cancer subtypes. Moreover, their down-regulation was associated with ER-, PR- and triple-negative breast cancer patients, as well as tumor progression (increased Scarff, Bloom and Richardson grade status and Nottingham Prognostic Index classification). Co-expression analysis revealed that METTL14 and ZC3H13 had a strong positive correlation with APC, an antagonist of the Wnt signaling pathway, indicating they might cooperate in regulating proliferation, invasion, and metastasis of tumor cells. METTL14, ZC3H13, and APC expression levels had significant positive correlation with infiltrating levels of CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells, and negative correlation with Treg cells in breast cancer.

Conclusions

This study demonstrated that down-regulation of METTL14 and ZC3H13 which act as two tumor suppressor genes was found in breast cancer and predicted poor prognosis. Their abnormal expression promoted breast cancer invasion by affecting pathways related to tumor progression and mediating immunosuppression.

Hypoxia-Associated Prognostic Markers and Competing Endogenous RNA Co-Expression Networks in Breast Cancer

Objective

Many primary tumors have insufficient supply of molecular oxygen, called hypoxia. Hypoxia is one of the leading characteristics of solid tumors resulting in a higher risk of local failure and distant metastasis. It is quite necessary to investigate the hypoxia associated molecular hallmarks in breast cancer.

Materials and Methods

According to the published studies, we selected 13 hypoxia related gene expression signature to define the hypoxia status of breast cancer using ConsensusClusterPlus package based on the data from The Cancer Genome Atlas (TCGA). Subsequently, we characterized the infiltration of 24 immune cell types under different hypoxic conditions. Furthermore, the differentially expressed hypoxia associated microRNAs, mRNAs and related signaling pathways were analyzed and depicted. On this basis, a series of prognostic markers related to hypoxia were identified and ceRNA co-expression networks were constructed.

Results

Two subgroups (cluster1 and cluster2) were identified and the 13 hypoxia related gene signature were all up-regulated in cluster1. Thus, we defined the cluster1 as “hypoxic subgroup” compared with cluster2. The infiltration of CD8+ T cell and CD4+ T cell were lower in cluster1 while the nTreg cell and iTreg cell were higher, indicating that there was immunosuppressive status in cluster1. We observed widespread hypoxia-associated dysregulation of microRNAs and mRNAs. Next, a risk signature for predicting prognosis of breast cancer patients was established based on 12 dysregulated hypoxia associated prognostic genes. Two microRNAs, hsa-miR-210-3p and hsa-miR-190b, with the most significant absolute logFC value were related to unfavorable and better prognosis, respectively. Several long non-coding RNAs were predicted to be microRNA targets and positively correlated with two selected mRNAs, CPEB2 and BCL11A. Predictions based on the LINC00899/PSMG3-AS1/PAXIP1-AS1- hsa-miR-210-3p-CPEB2 and SNHG16- hsa-miR-190b-BCL11A ceRNA regulation networks indicated that the two genes might act as tumor suppressor and oncogene, respectively.

Conclusion

Hypoxia plays an important role in the initiation and progression of breast cancer. Our research provides potential mechanisms into molecular-level understanding of tumor hypoxia.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

An integrative proteogenomics approach reveals peptides encoded by annotated lincRNA in the mouse kidney inner medulla

Long noncoding RNAs (lncRNAs) are intracellular transcripts longer than 200 nucleotides and lack protein-coding information. A subclass of lncRNA known as long intergenic noncoding RNAs (lincRNAs) are transcribed from genomic regions that share no overlap with annotated protein-coding genes. Increasing evidence has shown that some annotated lincRNA transcripts do in fact contain open reading frames (ORFs) encoding functional short peptides in the cell. Few robust methods for lincRNA-encoded peptide identification have been reported, and the tissue-specific expression of these peptides has been largely unexplored. Here we propose an integrative workflow for lincRNA-encoded peptide discovery and test it on the mouse kidney inner medulla (IM). In brief, low molecular weight protein fractions were enriched from homogenate of IMs and trypsinized into shorter peptides, which were sequenced by high resolution liquid chromatography-tandem mass spectrometry (LC-MS/MS). To curate a hypothetical lincRNA-encoded peptide database for peptide-spectrum matching following LC-MS/MS, we performed RNA-Seq on IMs, computationally removed reads overlapping with annotated protein-coding genes, and remapped the remaining reads to a database of mouse noncoding transcripts to infer lincRNA expression. Expressed lincRNAs were searched for ORFs by an existing rule-based algorithm, and translated ORFs were used for peptide-spectrum matching. Peptides identified by LC-MS/MS were further evaluated by using several quality control criteria and bioinformatics methods. We discovered three novel lincRNA-encoded peptides, which are conserved in mouse, rat, and human. The workflow can be adapted for discovery of small protein-coding genes in any species or tissue where noncoding transcriptome information is available.

Systemic Expression Analysis Reveals Prognostic Significance of WIPI3 in Hepatocellular Carcinoma

Introduction

WD repeat domain phosphoinositide-interacting protein 3 (WIPI3) is a member of the WIPI protein family, autophagy marker, that is associated with the malignant progression of various human cancers, but its role in hepatocellular carcinoma (HCC) is still unclear.

Materials and Methods

Firstly, we collected the mRNA expression of WIPI3 in HCC through the platform of Oncomine, as well as the DNA copy number variations (CNVs), and verified it on human HCC cell line and the GEO database. Then, the subgroups and prognosis of HCC were performed by the UALCAN web tool. The mutation of WIPI3 was analyzed by cBioPortal. The coexpression of WIPI3 in HCC was identified from the LinkedOmics database, and function enrichment analysis was done using the LinkFinder module in LinkedOmics. Coexpression gene network was constructed through the STRING database, and the MCODE plug-in of which was used to build the gene modules; both of them were visualized by the Cytoscape software. Finally, the top modular genes in the same patient cohort were constructed through data mining in The Cancer Genome Atlas (TCGA) liver hepatocellular carcinoma (LIHC) by using the UCSC Xena browser.

Results

The results indicated that WIPI3 was frequently overexpressed in HCC, which could lead to a poor prognosis through the Kaplan-Meier (KM) analysis. Moreover, there existed mutations of WIPI3 in HCC, and the prognosis of WIPI3-altered group was significantly poor based on KM plotter data. Coexpression analysis showed that the coexpression gene of WIPI3 was associated with cell cycle and spliceosome. Further analysis suggested that WIPI3 and eukaryotic translation initiation factor 4A3 (EIF4A3) coordinately regulated the cancer cell cycle by spliceosome as a result of the strong positive correlation between them.

Conclusion

In summary, WIPI3 is constantly overexpressed in HCC tissues, resulting in a poor prognosis; therefore, we can identify it as an effective target for the treatment of HCC.

GCH1 variants contribute to the risk and earlier age-at-onset of Parkinson's disease: a two-cohort case-control study

Background

Common and rare variants of guanosine triphosphate cyclohydrolase 1 (GCH1) gene may play important roles in Parkinson’s disease (PD). However, there is a lack of comprehensive analysis of GCH1 genotypes, especially in non-coding regions. The aim of this study was to explore the genetic characteristics of GCH1, including rare and common variants in coding and non-coding regions, in a large population of PD patients in Chinese mainland, as well as the phenotypic characteristics of GCH1 variant carriers.

Methods

In the first cohort of this case-control study, we performed whole-exome sequencing in 1555 patients with early-onset or familial PD and 2234 healthy controls; then in the second cohort, whole-genome sequencing was performed in sporadic late-onset PD samples (1962 patients), as well as 1279 controls. Variants at target GCH1 regions were extracted, and then genetic and detailed phenotypic data were analyzed using regression models and the sequence kernel association test. We also performed a meta-analysis to correlate deleterious GCH1 variants with age at onset (AAO) in PD patients.

Results

For coding variants, we identified a significant burden of GCH1 deleterious variants in early-onset or familial PD cases compared to controls (1.2% vs 0.1%, P < 0.0001). In the analysis of possible regulatory variants in GCH1 non-coding regions, rs12323905 (P = 0.001, odds ratio = 1.19, 95%CI 1.07–1.32) was significantly associated with PD, and variant sets in untranslated regions and intron regions, GCH1 brain-specific expression quantitative trait loci, and two possible promoter/enhancer (GH14J054857 and GH14J054880) were suggestively associated with PD. Genotype-phenotype correlation analysis revealed that the carriers of GCH1 deleterious variants manifested younger AAO (P < 0.0001), and had milder motor symptoms, milder fatigue symptoms and more autonomic nervous dysfunctions. Meta-analysis of six studies demonstrated 6.4-year earlier onset in GCH1 deleterious variant carriers (P = 0.0009).

Conclusions

The results highlight the importance of deleterious variants and non-coding variants of GCH1 in PD in Chinese mainland and suggest that GCH1 mutation can influence the PD phenotype, which may help design experimental studies to elucidate the mechanisms of GCH1 in the pathogenesis of PD.

Abnormal expression and prognostic significance of EPB41L1 in kidney renal clear cell carcinoma based on data mining

Background

EPB41L1 gene (erythrocyte membrane protein band 4.1 like 1) encodes the protein 4.1N, a member of 4.1 family, playing a vital role in cell adhesion and migration, which is associated with the malignant progression of various human cancers. However, the expression and prognostic significance of EPB41L1 in kidney renal clear cell carcinoma (KIRC) remain to be investigated.

Methods

In this study, we collected the mRNA expression of EPB41L1 in KIRC through the Oncomine platform, and used the HPA database to perform the pathological tissue immunohistochemistry in patients. Then, the sub-groups and prognosis of KIRC were performed by UALCAN and GEPIA web-tool, respectively. Further, the mutation of EPB41L1 in KIRC was analyzed by c-Bioportal. The co-expression genes of EPB41L1 in KIRC were displayed from the LinkedOmics database, and function enrichment analysis was used by LinkFinder module in LinkedOmics. The function of EPB41L1 in cell adhesion and migration was confirmed by wound healing assay using 786-O cells in vitro. Co-expression gene network was constructed through the STRING database, and the MCODE plug-in of which was used to build the gene modules, both of them was visualized by Cytoscape software. Finally, the top modular genes in the same patient cohort were constructed through data mining in TCGA by using the UCSC Xena browser.

Results

The results indicated that EPB41L1 was down-expressed in KIRC, leading to a poor prognosis. Moreover, there is a mutation in the FERM domain of EPB41L1, but it has no significant effect on the prognosis of KIRC. The co-expressed genes of EPB41L1 were associated with cell adhesion and confirmed in vitro. Further analysis suggested that EPB41L1 and amyloid beta precursor protein (APP) were coordinated to regulated cancer cell adhesion, thereby increasing the incidence of cancer cell metastasis and tumor invasion.

Conclusions

In summary, EPB41L1 is constantly down-expressed in KIRC tissues, resulting a poor prognosis. Therefore, we suggest that it can be an effective biomarker for the diagnosis of KIRC.

A scalable SCENIC workflow for single-cell gene regulatory network analysis

This protocol explains how to perform a fast SCENIC analysis alongside standard best practices steps on single-cell RNA-sequencing data using software containers and Nextflow pipelines. SCENIC reconstructs regulons (i.e., transcription factors and their target genes) assesses the activity of these discovered regulons in individual cells and uses these cellular activity patterns to find meaningful clusters of cells. Here we present an improved version of SCENIC with several advances. SCENIC has been refactored and reimplemented in Python (pySCENIC), resulting in a tenfold increase in speed, and has been packaged into containers for ease of use. It is now also possible to use epigenomic track databases, as well as motifs, to refine regulons. In this protocol, we explain the different steps of SCENIC: the workflow starts from the count matrix depicting the gene abundances for all cells and consists of three stages. First, coexpression modules are inferred using a regression per-target approach (GRNBoost2). Next, the indirect targets are pruned from these modules using cis-regulatory motif discovery (cisTarget). Lastly, the activity of these regulons is quantified via an enrichment score for the regulon's target genes (AUCell). Nonlinear projection methods can be used to display visual groupings of cells based on the cellular activity patterns of these regulons. The results can be exported as a loom file and visualized in the SCope web application. This protocol is illustrated on two use cases: a peripheral blood mononuclear cell data set and a panel of single-cell RNA-sequencing cancer experiments. For a data set of 10,000 genes and 50,000 cells, the pipeline runs in <2 h.

Discovering candidate imprinted genes and imprinting control regions in the human genome

Background

Genomic imprinting is a process thereby a subset of genes is expressed in a parent-of-origin specific manner. This evolutionary novelty is restricted to mammals and controlled by genomic DNA segments known as Imprinting Control Regions (ICRs) and germline Differentially Methylated Regions (gDMRs). Previously, I showed that in the mouse genome, the fully characterized ICRs/gDMRs often includes clusters of 2 or more of a set of composite-DNA-elements known as ZFBS-morph overlaps.

Results

Because of the importance of the ICRs to regulating parent-of-origin specific gene expression, I developed a genome-wide strategy for predicting their positions in the human genome. My strategy consists of creating plots to display the density of ZFBS-morph overlaps along the entire chromosomal DNA sequences. In initial evaluations, I found that peaks in these plots pinpointed several of the known ICRs/gDMRs along the DNA in chromosomal bands. I deduced that in density-plots, robust peaks corresponded to actual or candidate ICRs in the DNA. By locating the genes in the vicinity of candidate ICRs, I could discover potential imprinting genes. Additionally, my assessments revealed a connection between several of the potential imprinted genes and human developmental anomalies. Examples include Leber congenital amaurosis 11, Coffin-Siris syndrome, progressive myoclonic epilepsy-10, microcephalic osteodysplastic primordial dwarfism type II, and microphthalmia, cleft lip and palate, and agenesis of the corpus callosum.

Conclusion

With plots displaying the density of ZFBS-morph overlaps, researchers could locate candidate ICRs and imprinted genes. Since the datafiles are available for download and display at the UCSC genome browser, it is possible to examine the plots in the context of Single nucleotide polymorphisms (SNPs) to design experiments to discover novel ICRs and imprinted genes in the human genome.

Molecular Mechanism of Aflatoxin-Induced Hepatocellular Carcinoma Derived from a Bioinformatics Analysis

Exposure to aflatoxin is considered to be one of the causes of hepatocellular carcinoma (HCC). With the development of bioinformation, we sought to reveal the occurrence and development of aflatoxin-induced HCC through data research. We identified differentially expressed genes (DEGs) of datasets GSE127791 (Aflatoxin-treated pluripotent stem cell derived human hepatocytes vs. controls) and GSE64041 (liver carcinoma with unknown cause vs. non-cancerous tissue) by GEO2R to find the common DEGs. Gene ontology (GO) and KEGG path enrichment analysis were used to annotate the function of DEGs. Hub genes were screened from identified DEGs by protein-protein interaction (PPI) network analysis. The prognostic value of hub genes in cancer databases were evaluated. We obtained 132 common DEGs and 11 hub genes. According to cluster analysis and protein co-expression networks, we screened out the key genes, histidine-rich glycoprotein (HRG) and phosphoenolpyruvate carboxykinase 2 (PCK2). Oncomine database and survival curve analysis showed that the decline in HRG and PCK2 expression in the development of HCC indicated poor prognosis. We speculated that the decreased expression of HRG and PCK2 after aflatoxin exposure to hepatocyte may be related to aflatoxin induced hepatocyte injury and carcinogenesis. In addition, the decreased expression of HRG and PCK2 in the occurrence and development of HCC suggests a poor prognosis of HCC.

ATPase copper transporter A, negatively regulated by miR-148a-3p, contributes to cisplatin resistance in breast cancer cells

BACKGROUND

Breast cancer is the leading cause of death among women. Cisplatin is an effective drug for breast cancer, but resistance often develops during long term chemotherapy. While the mechanism of chemotherapy resistance is still not fully understood.

METHODS

Survival analyses of ATP7A and ATP7B were used to evaluate their effects on the development of Breast invasive carcinoma (BRCA). Immunostaining, flow cytometry, and IC50 assay were utilized to examine the effects of ATP7A-siRNA combined with cisplatin on apoptosis in breast cancer cells. Q-PCR, western blotting, and dual-luciferase assay were employed to confirm ATP7A is a novel target gene of miR-148a-3p.

RESULTS

In this current study, we identified knocking-down ATP7A could enhance cytotoxicity treatment of cisplatin in breast cancer cells. We also demonstrated miR-148a-3p overexpression in BRCA cells increased the sensitivity to cisplatin, and subsequently enhanced DNA damage and apoptosis. Moreover, we found ATP7A is a novel target gene of miR-148a-3p. In brief, our results showed miR-148a could accelerate chemotherapy induced-apoptosis in breast cancer cells by inhibiting ATP7A expression.

CONCLUSIONS

Our results highlight that inhibition of ATP7A is a potential strategy for targeting breast cancer resistant to cisplatin, and we provided an interesting method to compare the involvement of various genes in the assessment of cisplatin resistance.

Systematic expression analysis of WEE family kinases reveals the importance of PKMYT1 in breast carcinogenesis

OBJECTIVES

Many cancer cells depend on G2 checkpoint mechanism regulated by WEE family kinases to maintain genomic integrity. The PKMYT1 gene, as a member of WEE family kinases, participates in G2 checkpoint surveillance and probably links with tumorigenesis, but its role in breast cancer remains largely unclear.

MATERIALS AND METHODS

In this study, we used a set of bioinformatic tools to jointly analyse the expression of WEE family kinases and investigate the prognostic value of PKMYT1 in breast cancer.

RESULTS

The results indicated that PKMYT1 is the only frequently overexpressed member of WEE family kinases in breast cancer. KM plotter data suggests that abnormally high expression of PKMYT1 predicts poor prognosis, especially for some subtypes, such as luminal A/B and triple-negative (TNBC) types. Moreover, the up-regulation of PKMYT1 was associated with HER2-positive (HER2+), basal-like (Basal-like), TNBC statuses and increased classifications of Scarff, Bloom and Richardson (SBR). Co-expression analysis showed PKMYT1 has a strong positive correlation with Polo-like kinase 1 (PLK1), implying they may cooperate in regulating cancer cell proliferation by synchronizing rapid cell cycle with high quality of genome maintenance.

CONCLUSIONS

Collectively, this study demonstrates that overexpression of PKMYT1 is always found in breast cancer and predicts unfavourable prognosis, implicating it as an appealing therapeutic target for breast carcinoma.

Computational Approaches to Predict the Non-canonical DNAs

Background:

Although most nucleotides in the genome form canonical double-stranded B-DNA, many repeated sequences transiently present as non-canonical conformations (non-B DNA) such as triplexes, quadruplexes, Z-DNA, cruciforms, and slipped/hairpins. Those noncanonical DNAs (ncDNAs) are not only associated with many genetic events such as replication, transcription, and recombination, but are also related to the genetic instability that results in the predisposition to disease. Due to the crucial roles of ncDNAs in cellular and genetic functions, various computational methods have been implemented to predict sequence motifs that generate ncDNA.

Objective:

Here, we review strategies for the identification of ncDNA motifs across the whole genome, which is necessary for further understanding and investigation of the structure and function of ncDNAs.

Conclusion:

There is a great demand for computational prediction of non-canonical DNAs that play key functional roles in gene expression and genome biology. In this study, we review the currently available computational methods for predicting the non-canonical DNAs in the genome. Current studies not only provide an insight into the computational methods for predicting the secondary structures of DNA but also increase our understanding of the roles of non-canonical DNA in the genome.

TNFAIP3 genetic polymorphisms reduce ankylosing spondylitis risk in Eastern Chinese Han population

This study was conducted to clarify the associations of tumor necrosis factor-α induced protein 3 (TNFAIP3) and TNFAIP3-interacting protein 1 (TNIP1) genetic polymorphisms with ankylosing spondylitis (AS) susceptibility. Five single nucleotide polymorphisms (SNPs) in TNFAIP3 gene and four in TNIP1 gene were genotyped in 667 AS patients and 667 matched healthy controls. Genotypes and haplotype analysis were conducted by using SPSS 23.0 and Haploview 4.2 software. The T allele and CT genotype in TNFAIP3 rs10499194 were significantly associated with a reduced AS risk (T allele vs. C allele, OR = 0.619, 95% CI = 0.430–0.889, P = 0.009; CT vs. CC, OR = 0.603, 95% CI = 0.416–0.875, P = 0.007). However, no association remained significant after Bonferroni correction. The rs13207033^A- rs10499194^T haplotype of TNFAIP3 conferred a protective effect on AS susceptibility. Stratification analyses suggested that rs10499194 polymorphism decreased the risk of AS in the male subgroup, subgroup aged ≥ 29, HLA-B27 positive subgroup as well as the subgroups of BASFI < 4 and BASDAI < 4 (all P < 0.05). Furthermore, the functional annotation suggested a potential function of rs10499194 mutation. Our results demonstrated that TNFAIP3 rs10499194 polymorphism may be associated with a reduced risk of AS.

Up-Regulation of PI 3-Kinases and the Activation of PI3K-Akt Signaling Pathway in Cancer Stem-Like Cells Through DNA Hypomethylation Mediated by the Cancer Microenvironment

Previously, we have succeeded in converting induced pluripotent stem cells (iPSCs) into cancer stem cells (CSCs) by treating the iPSCs with conditioned medium of Lewis lung carcinoma (LLC) cells. The converted CSCs, named miPS-LLCcm cells, exhibited the self-renewal, differentiation potential, and potential to form malignant tumors with metastasis. In this study, we further characterized miPS-LLCcm cells both in vivo and in vitro. The tumors formed by subcutaneous injection showed the structures with pathophysiological features consisting of undifferentiated and malignant phenotypes generally found in adenocarcinoma. Metastasis in the lung was also observed as nodule structures. Excising from the tumors, primary cultured cells from the tumor and the nodule showed self-renewal, differentiation potential as well as tumor forming ability, which are the essential characters of CSCs. We then characterized the epigenetic regulation occurring in the CSCs. By comparing the DNA methylation level of CG rich regions, the differentially methylated regions (DMRs) were evaluated in all stages of CSCs when compared with the parental iPSCs. In DMRs, hypomethylation was found superior to hypermethylation in the miPS-LLCcm cells and its derivatives. The hypo- and hypermethylated genes were used to nominate KEGG pathways related with CSC. As a result, several categories were defined in the KEGG pathways from which most related with cancers, significant and high expression of components was PI3K-AKT signaling pathway. Simultaneously, the AKT activation was also confirmed in the CSCs. The PI3K-Akt signaling pathway should be an important pathway for the CSCs established by the treatment with conditioned medium of LLC cells.

Datasets on the genomic positions of the MLL1 morphemes, the ZFP57 binding site, and ZFBS-Morph overlaps in the build mm9 of the mouse genome

While MLL1 activates gene expression in most tissues, ZFP57 represses transcription. MLL1 selectively interacts with a group of nonmethylated DNA sequences known as the MLL1 morphemes. ZFP57 associates with a methylated hexamer (ZFBS), dispersed in the genomic DNA segments known as Imprinted Control Regions (ICRs) and germline Differentially Methylated Regions (gDMRs), to maintain allele-specific gene repression. We have identified a set of composite DNA elements (ZFBS-Morph overlaps) that provides the sequence context of ZFBS in the canonical ICRs/gDMRs. This report provides tables listing the nucleotide sequences of the MLL1 morphemes and ZFBS-Morph overlaps. The report also offers links to the data repository at Purdue University, for downloading the positions of the MLL1 morphemes, the ZFP57 binding site, and the ZFBS-Morph overlaps in the mouse genome.

Imprinted control regions include composite DNA elements consisting of the ZFP57 binding site overlapping MLL1 morphemes

Mammalian genomes include DNA segments that are imprinted (CpG-methylated) only on one of the two parental chromosomes, leading to parent-of-origin-specific gene expression. The process is regulated by Imprinting Control Regions (ICRs) and germline Differentially Methylated Regions (gDMRs). Previously, ZFP57 was shown to recognize a methylated hexanucleotide in ICRs to maintain allele-specific gene repression. In Bioinformatics analyses, I found that the hexamer occurred frequently in mouse chromosomal DNA, suggesting that beside the ZFP57 binding site (ZFBS), ICRs contained sequence features with unknown characteristics. To identify such features, I examined chromosomal abundance of motifs in which the length of the hexamer was extended by one or several nucleotides. Results led to the discovery of a group of functionally significant composite DNA elements (ZFBS-Morph overlaps) that may play dual roles in the regulation of allele-specific gene expression. Importantly, results of genome-wide evaluations revealed that nearly 90% of the gDMRs included closely-spaced ZFBS-Morph overlaps.

Identifying low density lipoprotein cholesterol associated variants in the Annexin A2 (ANXA2) gene

BACKGROUND AND AIMS

Annexin-A2 (AnxA2) is an endogenous inhibitor of proprotein convertase subtilisin/kexin type-9 (PCSK9). The repeat-one (R1) domain of AnxA2 binds to PCSK9, blocking its ability to promote degradation of low-density lipoprotein cholesterol-receptors (LDL-R) and thereby regulate low-density lipoprotein cholesterol (LDL-C) levels. Here we identify variants in ANXA2 influencing LDL-C levels and we determine the molecular mechanisms of their effects.

RESULTS

The ANXA2 single nucleotide polymorphism (SNP) genotype-phenotype association was examined using the Second-Northwick-Park Heart Study (NPHSII) (n∼2700) and the UCL-LSHTM-Edinburgh-Bristol (UCLEB) consortium (n∼14,600). The ANXA2-R1 domain coding-SNP rs17845226 (V98L) associated with LDL-C, homozygotes for the minor allele having ≈18.8% higher levels of LDL-C (p = 0.004), and higher risk of coronary heart disease (CHD) (p = 0.04). The SNP is in modest linkage disequilibrium (r² > 0.5) with two intergenic SNPs, rs17191344 and rs11633032. Both SNPs showed allele-specific protein binding, and the minor alleles caused significant reduction in reporter gene expression (≈18%, p < 0.001). In the expression quantitative trait loci (eQTL) study, minor allele homozygotes have significantly lower levels of ANXA2-mRNA expression (p = 1.36 × 10^-05).

CONCLUSIONS

Both rs11633032 and rs17191344 SNPs are functional variants, where the minor alleles create repressor-binding protein sites for transcription factors that contribute to reduced ANXA2 gene expression. Lower AnxA2 levels could increase plasma levels of PCSK9 and thus increase LDL-C levels and risk of CHD. This supports, for the first time in humans, previous observations in mouse models that changes in the levels of AnxA2 directly influence plasma LDL-C levels, and thus implicate this protein as a potential therapeutic target for LDL-C lowering.

Identification of the Functional Variant(s) that Explain the Low-Density Lipoprotein Receptor (LDLR) GWAS SNP rs6511720 Association with Lower LDL-C and Risk of CHD

Background

The Low-Density Lipoprotein Receptor (LDLR) SNP rs6511720 (G>T), located in intron-1 of the gene, has been identified in genome-wide association studies (GWAS) as being associated with lower plasma levels of LDL-C and a lower risk of coronary heart disease (CHD). Whether or not rs6511720 is itself functional or a marker for a functional variant elsewhere in the gene is not known.

Methods

The association of LDLR SNP rs6511720 with incidence of CHD and levels of LDL-C was determined by reference to CARDIoGRAM, C4D and Global lipids genetics consortium (GLGC) data. SNP annotation databases were used to identify possible SNP function and prioritization. Luciferase reporter assays in the liver cell line Huh7 were used to measure the effect of variant genotype on gene expression. Electrophoretic Mobility Shift Assays (EMSAs) were used to identify the Transcription Factors (TFs) involved in gene expression regulation.

Results

The phenotype-genotype analysis showed that the rs6511720 minor allele is associated with lower level of LDL-C [beta = -0.2209, p = 3.85 x10^-262], and lower risk of CHD [log (OR) = 0.1155, p = 1.04 x10^-7]. Rs6511720 is in complete linkage. Rs6511720 is in complete linkage disequilibrium (LD) with three intron-1 SNPs (rs141787760, rs60173709, rs57217136). Luciferase reporter assays in Huh7 cells showed that the rare alleles of both rs6511720 and rs57217136 caused a significant increase in LDLR expression compared to the common alleles (+29% and +24%, respectively). Multiplex Competitor-EMSAs (MC-EMSA) identified that the transcription factor serum response element (SRE) binds to rs6511720, while retinoic acid receptor (RAR) and signal transducer and activator of transcription 1 (STAT1) bind to rs57217136.

Conclusion

Both LDLR rs6511720 and rs57217136 are functional variants. Both these minor alleles create enhancer-binding protein sites for TFs and may contribute to increased LDLR expression, which is consequently associated with reduced LDL-C levels and 12% lower CHD risk.

SCN8A encephalopathy: Research progress and prospects

On April 21, 2015, the first SCN8A Encephalopathy Research Group convened in Washington, DC, to assess current research into clinical and pathogenic features of the disorder and prepare an agenda for future research collaborations. The group comprised clinical and basic scientists and representatives of patient advocacy groups. SCN8A encephalopathy is a rare disorder caused by de novo missense mutations of the sodium channel gene SCN8A, which encodes the neuronal sodium channel Nav 1.6. Since the initial description in 2012, approximately 140 affected individuals have been reported in publications or by SCN8A family groups. As a result, an understanding of the severe impact of SCN8A mutations is beginning to emerge. Defining a genetic epilepsy syndrome goes beyond identification of molecular etiology. Topics discussed at this meeting included (1) comparison between mutations of SCN8A and the SCN1A mutations in Dravet syndrome, (2) biophysical properties of the Nav 1.6 channel, (3) electrophysiologic effects of patient mutations on channel properties, (4) cell and animal models of SCN8A encephalopathy, (5) drug screening strategies, (6) the phenotypic spectrum of SCN8A encephalopathy, and (7) efforts to develop a bioregistry. A panel discussion of gaps in bioregistry, biobanking, and clinical outcomes data was followed by a planning session for improved integration of clinical and basic science research. Although SCN8A encephalopathy was identified only recently, there has been rapid progress in functional analysis and phenotypic classification. The focus is now shifting from identification of the underlying molecular cause to the development of strategies for drug screening and prioritized patient care.

Systematic characterization and comparison of the CYP2C9 variability of the Orang Asli in Malaysia with 12 populations

We conducted a systematic characterization of CYP2C9 variants in 61 Orang Asli and 96 Singaporean Malays using the whole genome sequences data and compared the variants with the other 11 HapMap populations. The frequency of rs1057910 (CYP2C9*3) is the highest in the Orang Asli compared to other populations. Three alleles with clinical implication were detected in the Orang Asli while 2 were found in the Singaporean Malays. Large numbers of the Orang Asli are predicted to have reduced metabolic capacity and therefore they would require a lower dose of drugs which are metabolized by CYP2C9. They are also at increased risks of adverse effects and therapeutic failures. A large number of CYP2C9 variants in the Orang Asli were not in the Hardy Weinberg Equilibrium which could be due to small sample size or mutations that disrupt the equilibrium of allele frequencies. In conclusion, different polymorphism patterns, allele frequencies, genotype frequencies and LD blocks are observed between the Orang Asli, the Singaporean Malays and the other populations. The study provided new information on the genetic polymorphism of CYP2C9 which is important for the implementation of precision medicine for the Orang Asli.

sORFs.org: a repository of small ORFs identified by ribosome profiling

With the advent of ribosome profiling, a next generation sequencing technique providing a "snap-shot'' of translated mRNA in a cell, many short open reading frames (sORFs) with ribosomal activity were identified. Follow-up studies revealed the existence of functional peptides, so-called micropeptides, translated from these 'sORFs', indicating a new class of bio-active peptides. Over the last few years, several micropeptides exhibiting important cellular functions were discovered. However, ribosome occupancy does not necessarily imply an actual function of the translated peptide, leading to the development of various tools assessing the coding potential of sORFs. Here, we introduce sORFs.org (http://www.sorfs.org), a novel database for sORFs identified using ribosome profiling. Starting from ribosome profiling, sORFs.org identifies sORFs, incorporates state-of-the-art tools and metrics and stores results in a public database. Two query interfaces are provided, a default one enabling quick lookup of sORFs and a BioMart interface providing advanced query and export possibilities. At present, sORFs.org harbors 263 354 sORFs that demonstrate ribosome occupancy, originating from three different cell lines: HCT116 (human), E14_mESC (mouse) and S2 (fruit fly). sORFs.org aims to provide an extensive sORFs database accessible to researchers with limited bioinformatics knowledge, thus enabling easy integration into personal projects.

Long-term reduction of T-cell intracellular antigens reveals a transcriptome associated with extracellular matrix and cell adhesion components

Knockdown of T-cell intracellular antigens TIA1 and TIAR contributes to a cellular phenotype characterised by uncontrolled proliferation and tumorigenesis. Massive-scale poly(A+) RNA sequencing of TIA1 or TIAR-knocked down HeLa cells reveals transcriptome signatures comprising genes and functional categories potentially able to modulate several aspects of membrane dynamics associated with extracellular matrix and focal/cell adhesion events. The transcriptomic heterogeneity is the result of differentially expressed genes and RNA isoforms generated by alternative splicing and/or promoter usage. These results suggest a role for TIA proteins in the regulation and/or modulation of cellular homeostasis related to focal/cell adhesion, extracellular matrix and membrane and cytoskeleton dynamics.