IgA nephropathy(IgAN) is the most common primary glomerular disease in China. Primary infections always occur before IgAN. However, the pathology of IgAN is still unclear. Previously we found that LL37, a protein secreted by senescent cells, was specific for the progression of IgAN, and also played a role in the neutrophil function. So we hypothesized that the infiltration of neutrophils, inflammation factors, and aging markers , which were modulated by functional networks, induced the immune response and renal injury. RNA-Sequencing (RNA-seq) can be used to study the whole transcriptome and detect splicing variants that are expressed in a specific cell type or tissue. We separate glomerulus from the renal biopsy tissues. After RNA extraction, the sequences were analyzed with Illumina HiSeq 2000/2500. 381 genes with differential expression between the IgAN patients and the healthy controls were identified. Only PLAU, JUN, and FOS were related to DNA damage, telomere dysfunction-induced aging markers, neutrophil function and IgA nephropathy. The networks showed the possibility of these genes being connected. We conclude that DNA damage and telomere dysfunction could play important roles in IgA nephropathy. In addition, neutrophils are also important factors in this disease. The networks of these markers showed the mechanism pathways that are involved in the duration of the occurrence and progression of IgA nephropathy and might be a new therapeutic opportunity for disease treatment.

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r(2) = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer.

Intrauterine growth restriction (IUGR) affects the foetus and has a number of pathological consequences throughout life. Recent work has indicated that variations in DNA methylation might cause placental dysfunction, which may be associated with adverse pregnancy complications. Here, we investigated the promoter methylomes of placental shares from seven monochorionic (MC) twins with selective intrauterine growth restriction (sIUGR) using the healthy twin as an ideal control. Our work demonstrated that the IUGR placental shares harboured a distinct DNA hypomethylation pattern and that the methylation variations preferentially occurred in CpG island shores or non-CpG island promoters. The differentially methylated promoters could significantly separate the IUGR placental shares from the healthy ones. Ultra-performance liquid chromatography/tandem mass spectrometry (UPLC-MS/MS) further confirmed the genome-wide DNA hypomethylation and the lower level of hydroxymethylation statuses in the IUGR placental shares. The methylation variations of the LRAT and SLC19A1 promoters, which are involved in vitamin A metabolism and folate transportation, respectively, and the EFS promoter were further validated in an additional 12 pairs of MC twins with sIUGR. Although the expressions of LRAT, SLC19A1 and EFS were not affected, we still speculated that DNA methylation and hydroxymethylation might serve a functional role during in utero foetal development.

Runs of homozygosity (ROH) represents extended length of homozygotes on a long genomic distance. In oncology, it is known as loss of heterozygosity (LOH) if identified exclusively in cancer cell rather than in matched control cell. Studies have identified several genomic regions which show consistent ROH in different kinds of carcinoma. To query whether this consistency can be observed on broader spectrum, both in more cancer types and in wider genomic regions, we investigated ROH patterns in the National Cancer Institute 60 cancer cell line panel (NCI-60) and HapMap Caucasian healthy trio families. Using results from Affymetrix 500 K SNP arrays, we report a genome wide significant association of ROH regions between the NCI-60 and HapMap samples, with much a higher level of ROH (11 fold) in the cancer cell lines. Analysis shows that more severe ROH found in cancer cells appears to be the extension of existing ROH in healthy state. In the HapMap trios, the adult subgroup had a slightly but significantly higher level (1.02 fold) of ROH than did the young subgroup. For several ROH regions we observed the co-occurrence of fragile sites (FRAs). However, FRA on the genome wide level does not show a clear relationship with ROH regions.

Objective. Bringing together genomics, transcriptomics, proteomics, and other -omics technologies is an important step towards developing highly personalized medicine. However, instrumentation has advances far beyond expectations and now we are able to generate data faster than it can be interpreted. Materials and Methods. We have developed PANDA (Pathway AND Annotation) Explorer, a visualization tool that integrates gene-level annotation in the context of biological pathways to help interpret complex data from disparate sources. PANDA is a web-based application that displays data in the context of well-studied pathways like KEGG, BioCarta, and PharmGKB. PANDA represents data/annotations as icons in the graph while maintaining the other data elements (i.e., other columns for the table of annotations). Custom pathways from underrepresented diseases can be imported when existing data sources are inadequate. PANDA also allows sharing annotations among collaborators. Results. In our first use case, we show how easy it is to view supplemental data from a manuscript in the context of a user's own data. Another use-case is provided describing how PANDA was leveraged to design a treatment strategy from the somatic variants found in the tumor of a patient with metastatic sarcomatoid renal cell carcinoma. Conclusion. PANDA facilitates the interpretation of gene-centric annotations by visually integrating this information with context of biological pathways. The application can be downloaded or used directly from our website: http://bioinformaticstools.mayo.edu/research/panda-viewer/.

Aims/Hypothesis: Transforming growth factor-beta (TGF-β1) plays an important regulatory role in the progression of chronic kidney failure. Further, damage to kidney glomerular mesangial cells is central to the progression of diabetic nephropathy. The aim of this study was to explore the genetic associations between mRNA, microRNA, and epigenetics in mesangial cells in response to TGF-β1. Methods: The regulatory effects of TGF-β1 on mesangial cells were investigated at different molecular levels by treating mesangial cells with TGF-β1 for 3 days followed by genome-wide miRNA, RNA, DNA methylation, and H3K27me3 expression profiling using next generation sequencing (NGS). Results: Our results provide the first comprehensive, computationally integrated report of RNA-Seq, miRNA-Seq, and epigenomic analyses across all genetic variations, confirming the occurrence of DNA methylation and H3K27me3 in response to TGF-β1. Our findings show that the expression of KLF7 and Gja4 are involved in TGF-β1 regulated DNA methylation. Our data also provide evidence of the association between epigenetic changes and the expression of genes closely related to TGF-β1 regulation. Conclusion: This study has advanced our current knowledge of mechanisms that contribute to the expression of TGF-β1-regulated genes involved in the pathogenesis of kidney disease. The molecular underpinnings of TGF-β1 stimulation of kidney cells was determined, thereby providing a robust platform for further target exploration.

Human mesenchymal stem cell (MSC)-based tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) gene delivery is regarded as an effective treatment for glioblastoma (GBM). However, adverse-free target site homing of the delivery vehicles to the tumor microsatellite nests is challenging, leading to erroneously sustained released of this suicide protein into the normal brain parenchyma; therefore, limiting off-target cytotoxicity and controlled expression of the suicide inductor is a prerequisite for the safe use of therapeutic stem cells.

Genome-wide association studies have generated over thousands of susceptibility loci for many human complex traits, and yet for most of these associations the true causal variants remain unknown. Tissue/cell type-specific prediction and prioritization of non-coding regulatory variants will facilitate the identification of causal variants and underlying pathogenic mechanisms for particular complex diseases and traits. By leveraging recent large-scale functional genomics/epigenomics data, we develop an intuitive web server, GWAS4D (http://mulinlab.tmu.edu.cn/gwas4d or http://mulinlab.org/gwas4d), that systematically evaluates GWAS signals and identifies context-specific regulatory variants. The updated web server includes six major features: (i) updates the regulatory variant prioritization method with our new algorithm; (ii) incorporates 127 tissue/cell type-specific epigenomes data; (iii) integrates motifs of 1480 transcriptional regulators from 13 public resources; (iv) uniformly processes Hi-C data and generates significant interactions at 5 kb resolution across 60 tissues/cell types; (v) adds comprehensive non-coding variant functional annotations; (vi) equips a highly interactive visualization function for SNP-target interaction. Using a GWAS fine-mapped set for 161 coronary artery disease risk loci, we demonstrate that GWAS4D is able to efficiently prioritize disease-causal regulatory variants.

Characterizing the binding preference of RNA-binding proteins (RBP) is essential for us to understand the interaction between an RBP and its RNA targets, and to decipher the mechanism of post-transcriptional regulation. Experimental methods have been used to generate protein-RNA binding data for a number of RBPs in vivo and in vitro. Utilizing the binding data, a couple of computational methods have been developed to detect the RNA sequence or structure preferences of the RBPs. However, the majority of RBPs have not yet been experimentally characterized and lack RNA binding data. For these poorly studied RBPs, the identification of their binding preferences cannot be performed by most existing computational methods because the experimental binding data are prerequisite to these methods.

As the cost of whole genome sequencing (WGS) decreases, clinical laboratories will be looking at broadly adopting this technology to screen for variants of clinical significance. To fully leverage this technology in a clinical setting, results need to be reported quickly, as the turnaround rate could potentially impact patient care. The latest sequencers can sequence a whole human genome in about 24 hours. However, depending on the computing infrastructure available, the processing of data can take several days, with the majority of computing time devoted to aligning reads to genomics regions that are to date not clinically interpretable. In an attempt to accelerate the reporting of clinically actionable variants, we have investigated the utility of a multi-step alignment algorithm focused on aligning reads and calling variants in genomic regions of clinical relevance prior to processing the remaining reads on the whole genome. This iterative workflow significantly accelerates the reporting of clinically actionable variants with no loss of accuracy when compared to genotypes obtained with the OMNI SNP platform or to variants detected with a standard workflow that combines Novoalign and GATK.

Structural variation (SV) represents a significant, yet poorly understood contribution to an individual's genetic makeup. Advanced next-generation sequencing technologies are widely used to discover such variations, but there is no single detection tool that is considered a community standard. In an attempt to fulfil this need, we developed an algorithm, SoftSearch, for discovering structural variant breakpoints in Illumina paired-end next-generation sequencing data. SoftSearch combines multiple strategies for detecting SV including split-read, discordant read-pair, and unmated pairs. Co-localized split-reads and discordant read pairs are used to refine the breakpoints.

Conventionally, overall gene expressions from microarrays are used to infer gene networks, but it is challenging to account splicing isoforms. High-throughput RNA Sequencing has made splice variant profiling practical. However, its true merit in quantifying splicing isoforms and isoform-specific exon expressions is not well explored in inferring gene networks. This study demonstrates SpliceNet, a method to infer isoform-specific co-expression networks from exon-level RNA-Seq data, using large dimensional trace. It goes beyond differentially expressed genes and infers splicing isoform network changes between normal and diseased samples. It eases the sample size bottleneck; evaluations on simulated data and lung cancer-specific ERBB2 and MAPK signaling pathways, with varying number of samples, evince the merit in handling high exon to sample size ratio datasets. Inferred network rewiring of well established Bcl-x and EGFR centered networks from lung adenocarcinoma expression data is in good agreement with literature. Gene level evaluations demonstrate a substantial performance of SpliceNet over canonical correlation analysis, a method that is currently applied to exon level RNA-Seq data. SpliceNet can also be applied to exon array data. SpliceNet is distributed as an R package available at http://www.jjwanglab.org/SpliceNet.

Chromatin immunoprecipitation (ChIP) followed by next-generation sequencing (ChIP-Seq) has been widely used to identify genomic loci of transcription factor (TF) binding and histone modifications. ChIP-Seq data analysis involves multiple steps from read mapping and peak calling to data integration and interpretation. It remains challenging and time-consuming to process large amounts of ChIP-Seq data derived from different antibodies or experimental designs using the same approach. To address this challenge, there is a need for a comprehensive analysis pipeline with flexible settings to accelerate the utilization of this powerful technology in epigenetics research.

Large body of evidences accumulated in clinical and epidemiological studies indicate that hearts of diabetic subjects are more sensitive to ischemia reperfusion injury (IRI), which results in a higher rate of mortality at post-operation than that of non-diabetes. However, experimental results are equivocal and point to either increased or decreased susceptibility of the diabetic hearts to IRI, especially at the early stage of the disease. The present study was designed to test the hypothesis that the duration/severity of the indexed ischemia is a major determinant of the vulnerability to myocardial IRI at early stage of diabetes.

Mammalian target of rapamycin (mTOR) represents a key downstream intermediate for a myriad of oncogenic receptor tyrosine kinases. In the case of the insulin-like growth factor (IGF) pathway, the mTOR complex (mTORC1) mediates IGF-1 receptor (IGF-1R)-induced estrogen receptor alpha (ERα) phosphorylation/activation and leads to increased proliferation and growth in breast cancer cells. As a result, the prevalence of mTOR inhibitors combined with hormonal therapy has increased in recent years. Conversely, activated mTORC1 provides negative feedback regulation of IGF signaling via insulin receptor substrate (IRS)-1/2 serine phosphorylation and subsequent proteasomal degradation. Thus, the IGF pathway may provide escape (e.g. de novo or acquired resistance) from mTORC1 inhibitors. It is therefore plausible that combined inhibition of mTORC1 and IGF-1R for select subsets of ER-positive breast cancer patients presents as a viable therapeutic option.

Epigenetic alterations, such as aberrant DNA methylation of promoter and enhancer regions, which lead to atypical gene expression, have been associated with carcinogenesis. In hepatocellular carcinoma (HCC), genome-wide analysis of methylation has only recently been used. For a better understanding of hepatocarcinogenesis, we applied an even higher resolution analysis of the promoter methylome to identify previously unknown regions and genes differentially methylated in HCC.

It remains challenging to predict regulatory variants in particular tissues or cell types due to highly context-specific gene regulation. By connecting large-scale epigenomic profiles to expression quantitative trait loci (eQTLs) in a wide range of human tissues/cell types, we identify critical chromatin features that predict variant regulatory potential. We present cepip, a joint likelihood framework, for estimating a variant's regulatory probability in a context-dependent manner. Our method exhibits significant GWAS signal enrichment and is superior to existing cell type-specific methods. Furthermore, using phenotypically relevant epigenomes to weight the GWAS single-nucleotide polymorphisms, we improve the statistical power of the gene-based association test.

Interest in preclinical drug development for ovarian cancer has stimulated development of patient-derived xenograft (PDX) or tumorgraft models. However, the unintended formation of human lymphoma in severe combined immunodeficiency (SCID) mice from Epstein-Barr virus (EBV)-infected human lymphocytes can be problematic. In this study, we have characterized ovarian cancer PDXs which developed human lymphomas and explore methods to suppress lymphoproliferative growth. Fresh human ovarian tumors from 568 patients were transplanted intraperitoneally in SCID mice. A subset of PDX models demonstrated atypical patterns of dissemination with mediastinal masses, hepatosplenomegaly, and CD45-positive lymphoblastic atypia without ovarian tumor engraftment. Expression of human CD20 but not CD3 supported a B-cell lineage, and EBV genomes were detected in all lymphoproliferative tumors. Immunophenotyping confirmed monoclonal gene rearrangements consistent with B-cell lymphoma, and global gene expression patterns correlated well with other human lymphomas. The ability of rituximab, an anti-CD20 antibody, to suppress human lymphoproliferation from a patient's ovarian tumor in SCID mice and prevent growth of an established lymphoma led to a practice change with a goal to reduce the incidence of lymphomas. A single dose of rituximab during the primary tumor heterotransplantation process reduced the incidence of CD45-positive cells in subsequent PDX lines from 86.3% (n = 117 without rituximab) to 5.6% (n = 160 with rituximab), and the lymphoma rate declined from 11.1% to 1.88%. Taken together, investigators utilizing PDX models for research should routinely monitor for lymphoproliferative tumors and consider implementing methods to suppress their growth.

Genome-wide methylation profiling has led to more comprehensive insights into gene regulation mechanisms and potential therapeutic targets. Illumina Human Methylation BeadChip is one of the most commonly used genome-wide methylation platforms. Similar to other microarray experiments, methylation data is susceptible to various technical artifacts, particularly batch effects. To date, little attention has been given to issues related to normalization and batch effect correction for this kind of data.

The rapid development of next generation sequencing (NGS) technology provides a new chance to extend the scale and resolution of genomic research. How to efficiently map millions of short reads to the reference genome and how to make accurate SNP calls are two major challenges in taking full advantage of NGS. In this article, we reviewed the current software tools for mapping and SNP calling, and evaluated their performance on samples from The Cancer Genome Atlas (TCGA) project. We found that BWA and Bowtie are better than the other alignment tools in comprehensive performance for Illumina platform, while NovoalignCS showed the best overall performance for SOLiD. Furthermore, we showed that next-generation sequencing platform has significantly lower coverage and poorer SNP-calling performance in the CpG islands, promoter and 5'-UTR regions of the genome. NGS experiments targeting for these regions should have higher sequencing depth than the normal genomic region.