Zhiming Cai

A metagenome-wide association study of gut microbiota in type 2 diabetes

Assessment and characterization of gut microbiota has become a major research area in human disease, including type 2 diabetes, the most prevalent endocrine disease worldwide. To carry out analysis on gut microbial content in patients with type 2 diabetes, we developed a protocol for a metagenome-wide association study (MGWAS) and undertook a two-stage MGWAS based on deep shotgun sequencing of the gut microbial DNA from 345 Chinese individuals. We identified and validated approximately 60,000 type-2-diabetes-associated markers and established the concept of a metagenomic linkage group, enabling taxonomic species-level analyses. MGWAS analysis showed that patients with type 2 diabetes were characterized by a moderate degree of gut microbial dysbiosis, a decrease in the abundance of some universal butyrate-producing bacteria and an increase in various opportunistic pathogens, as well as an enrichment of other microbial functions conferring sulphate reduction and oxidative stress resistance. An analysis of 23 additional individuals demonstrated that these gut microbial markers might be useful for classifying type 2 diabetes.

Genome-wide association study in a Chinese Han population identifies nine new susceptibility loci for systemic lupus erythematosus

We performed a genome-wide association study (GWAS) of systemic lupus erythematosus (SLE) in a Chinese Han population by genotyping 1,047 cases and 1,205 controls using Illumina Human610-Quad BeadChips and replicating 78 SNPs in two additional cohorts (3,152 cases and 7,050 controls). We identified nine new susceptibility loci (ETS1, IKZF1, RASGRP3, SLC15A4, TNIP1, 7q11.23, 10q11.22, 11q23.3 and 16p11.2; 1.77 × 10−25 ≤ Pcombined ≤ 2.77 × 10−8) and confirmed seven previously reported loci (BLK, IRF5, STAT4, TNFAIP3, TNFSF4, 6q21 and 22q11.21; 5.17 × 10−42 ≤ Pcombined ≤ 5.18 × 10−12). Comparison with previous GWAS findings highlighted the genetic heterogeneity of SLE susceptibility between Chinese Han and European populations. This study not only advances our understanding of the genetic basis of SLE but also highlights the value of performing GWAS in diverse ancestral populations.

The genomic sequence of the Chinese hamster ovary (CHO)-K1 cell line

Chinese hamster ovary (CHO)–derived cell lines are the preferred host cells for the production of therapeutic proteins. Here we present a draft genomic sequence of the CHO-K1 ancestral cell line. The assembly comprises 2.45 Gb of genomic sequence, with 24,383 predicted genes. We associate most of the assembled scaffolds with 21 chromosomes isolated by microfluidics to identify chromosomal locations of genes. Furthermore, we investigate genes involved in glycosylation, which affect therapeutic protein quality, and viral susceptibility genes, which are relevant to cell engineering and regulatory concerns. Homologs of most human glycosylation-associated genes are present in the CHO-K1 genome, although 141 of these homologs are not expressed under exponential growth conditions. Many important viral entry genes are also present in the genome but not expressed, which may explain the unusual viral resistance property of CHO cell lines. We discuss how the availability of this genome sequence may facilitate genome-scale science for the optimization of biopharmaceutical protein production.

Frequent mutations of chromatin remodeling genes in transitional cell carcinoma of the bladder.

Transitional cell carcinoma (TCC) is the most common type of bladder cancer. Here we sequenced the exomes of nine individuals with TCC and screened all the somatically mutated genes in a prevalence set of 88 additional individuals with TCC with different tumor stages and grades. In our study, we discovered a variety of genes previously unknown to be mutated in TCC. Notably, we identified genetic aberrations of the chromatin remodeling genes (UTX, MLL-MLL3, CREBBP-EP300, NCOR1, ARID1A and CHD6) in 59% of our 97 subjects with TCC. Of these genes, we showed UTX to be altered substantially more frequently in tumors of low stages and grades, highlighting its potential role in the classification and diagnosis of bladder cancer. Our results provide an overview of the genetic basis of TCC and suggest that aberration of chromatin regulation might be a hallmark of bladder cancer.

Single-Cell Exome Sequencing Reveals Single-Nucleotide Mutation Characteristics of a Kidney Tumor

Clear cell renal cell carcinoma (ccRCC) is the most common kidney cancer and has very few mutations that are shared between different patients. To better understand the intratumoral genetics underlying mutations of ccRCC, we carried out single-cell exome sequencing on a ccRCC tumor and its adjacent kidney tissue. Our data indicate that this tumor was unlikely to have resulted from mutations in VHL and PBRM1. Quantitative population genetic analysis indicates that the tumor did not contain any significant clonal subpopulations and also showed that mutations that had different allele frequencies within the population also had different mutation spectrums. Analyses of these data allowed us to delineate a detailed intratumoral genetic landscape at a single-cell level. Our pilot study demonstrates that ccRCC may be more genetically complex than previously thought and provides information that can lead to new ways to investigate individual tumors, with the aim of developing more effective cellular targeted therapies.

MicroRNA Expression Signatures of Bladder Cancer Revealed by Deep Sequencing

Background MicroRNAs (miRNAs) are a class of small noncoding RNAs that regulate gene expression. They are aberrantly expressed in many types of cancers. In this study, we determined the genome-wide miRNA profiles in bladder urothelial carcinoma by deep sequencing. Methodology/Principal Findings We detected 656 differentially expressed known human miRNAs and miRNA antisense sequences (miRNA*s) in nine bladder urothelial carcinoma patients by deep sequencing. Many miRNAs and miRNA*s were significantly upregulated or downregulated in bladder urothelial carcinoma compared to matched histologically normal urothelium. hsa-miR-96 was the most significantly upregulated miRNA and hsa-miR-490-5p was the most significantly downregulated one. Upregulated miRNAs were more common than downregulated ones. The hsa-miR-183, hsa-miR-200b∼429, hsa-miR-200c∼141 and hsa-miR-17∼92 clusters were significantly upregulated. The hsa-miR-143∼145 cluster was significantly downregulated. hsa-miR-182, hsa-miR-183, hsa-miR-200a, hsa-miR-143 and hsa-miR-195 were evaluated by Real-Time qPCR in a total of fifty-one bladder urothelial carcinoma patients. They were aberrantly expressed in bladder urothelial carcinoma compared to matched histologically normal urothelium (p<0.001 for each miRNA). Conclusions/Significance To date, this is the first study to determine genome-wide miRNA expression patterns in human bladder urothelial carcinoma by deep sequencing. We found that a collection of miRNAs were aberrantly expressed in bladder urothelial carcinoma compared to matched histologically normal urothelium, suggesting that they might play roles as oncogenes or tumor suppressors in the development and/or progression of this cancer. Our data provide novel insights into cancer biology.

Integrated Profiling of MicroRNAs and mRNAs: MicroRNAs Located on Xq27.3 Associate with Clear Cell Renal Cell Carcinoma

Background With the advent of second-generation sequencing, the expression of gene transcripts can be digitally measured with high accuracy. The purpose of this study was to systematically profile the expression of both mRNA and miRNA genes in clear cell renal cell carcinoma (ccRCC) using massively parallel sequencing technology. Methodology The expression of mRNAs and miRNAs were analyzed in tumor tissues and matched normal adjacent tissues obtained from 10 ccRCC patients without distant metastases. In a prevalence screen, some of the most interesting results were validated in a large cohort of ccRCC patients. Principal Findings A total of 404 miRNAs and 9,799 mRNAs were detected to be differentially expressed in the 10 ccRCC patients. We also identified 56 novel miRNA candidates in at least two samples. In addition to confirming that canonical cancer genes and miRNAs (including VEGFA, DUSP9 and ERBB4; miR-210, miR-184 and miR-206) play pivotal roles in ccRCC development, promising novel candidates (such as PNCK and miR-122) without previous annotation in ccRCC carcinogenesis were also discovered in this study. Pathways controlling cell fates (e.g., cell cycle and apoptosis pathways) and cell communication (e.g., focal adhesion and ECM-receptor interaction) were found to be significantly more likely to be disrupted in ccRCC. Additionally, the results of the prevalence screen revealed that the expression of a miRNA gene cluster located on Xq27.3 was consistently downregulated in at least 76.7% of ∼50 ccRCC patients. Conclusions Our study provided a two-dimensional map of the mRNA and miRNA expression profiles of ccRCC using deep sequencing technology. Our results indicate that the phenotypic status of ccRCC is characterized by a loss of normal renal function, downregulation of metabolic genes, and upregulation of many signal transduction genes in key pathways. Furthermore, it can be concluded that downregulation of miRNA genes clustered on Xq27.3 is associated with ccRCC.

Long intergenic non-coding RNA TUG1 is overexpressed in urothelial carcinoma of the bladder

Long intergenic non‐coding RNAs (lincRNAs) are a class of non‐coding RNAs that regulate gene expression via chromatin reprogramming. Taurine Up‐regulated Gene 1 (TUG1) is a lincRNA that is associated with chromatin‐modifying complexes and plays roles in gene regulation. In this study, we determined the expression patterns of TUG1 and the cell proliferation inhibition and apoptosis induced by silencing TUG1 in urothelial carcinoma of the bladder.

Inducing Cell Proliferation Inhibition, Apoptosis, and Motility Reduction by Silencing Long Noncoding Ribonucleic Acid Metastasis-associated Lung Adenocarcinoma Transcript 1 in Urothelial Carcinoma of the Bladder

OBJECTIVEnTo study the expression patterns of long noncoding ribonucleic acid (RNA) metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) and the cell proliferation inhibition, apoptosis, and motility changes induced by silencing MALAT1 in urothelial carcinoma of the bladder.nnnMATERIALS AND METHODSnThe expression levels of MALAT1 were determined using real-time quantitative polymerase chain reaction in cancerous tissues and paired normal tissues in a total of 36 patients with urothelial carcinoma of the bladder. Expression differences were analyzed according to the grade and stage. Bladder urothelial carcinoma T24 and 5637 cells were transfected with MALAT1 small interfering RNA or negative control small interfering RNA. The cell proliferation changes of the transfected bladder urothelial carcinoma cells were determined using the MTT assay. Apoptosis caused by silencing MALAT1 was evaluated using the flow cytometry assay and enzyme-linked immunosorbent assay. The motility changes induced by silencing MALAT1 were measured using the wound healing assay.nnnRESULTSnMALAT1 was upregulated in bladder urothelial carcinoma compared with matched normal urothelium (P=.008). The MALAT1 expression levels were greater in high-grade carcinomas than in low-grade carcinoma (P=.001). The MALAT1 expression levels were greater in invasive carcinoma than in noninvasive carcinoma (P=.018). Cell proliferation inhibition, increased apoptosis, and decreased motility were observed in MALAT1 small interfering RNA-transfected bladder urothelial carcinoma T24 and 5637 cells.nnnCONCLUSIONnMALAT1 plays an oncogenic role in urothelial carcinoma of the bladder. Silencing MALAT1 is a potential novel therapeutic approach for this cancer.

Synthetic miRNA-Mowers Targeting miR-183-96-182 Cluster or miR-210 Inhibit Growth and Migration and Induce Apoptosis in Bladder Cancer Cells

Background MicroRNAs (miRNAs) function as endogenous regulators of biological behaviors of human cancers. Several natural non-coding RNAs are reported to inhibit miRNAs by base-pairing interactions. These phenomena raise questions about the ability of artificial device to regulate miRNAs. The purpose of this study is to create synthetic devices that target a single miRNA or a miRNA cluster and to ascertain their therapeutic effects on the phenotypes of bladder cancer cells. Methodology/Principal Findings Tandem bulged miRNA binding sites were inserted into the 3′ untranslated region (UTR) of the SV-40 promoter-driven Renilla luciferase gene to construct two “miRNA-mowers” for suppression of miR-183-96-182 cluster or miR-210. A third device with tandem repeat sequences not complementary to any known miRNA was generated as an untargeted-control. In functional analyses, bladder cancer T24 and UM-UC-3 cells were transfected with each of the three devices, followed by assays for detection of their impacts. Luciferase assays indicated that the activities of the luciferase reporters in the miRNA-mowers were decreased to 30–50% of the untargeted-control. Using Real-Time qPCR, the expression levels of the target miRNAs were shown to be reduced 2-3-fold by the corresponding miRNA-mower. Cell growth, apoptosis, and migration were tested by MTT assay, flow cytometry assay, and in vitro scratch assay, respectively. Cell growth inhibition, increased apoptosis, and decreased motility were observed in miRNA-mowers-transfected bladder cancer cells. Conclusions/Significance Not only a single target miRNA but also the whole members of a target miRNA cluster can be blocked using this modular design strategy. Anti-cancer effects are induced by the synthetic miRNA-mowers in the bladder cancer cell lines. miR-183/96/182 cluster and miR-210 are shown to play oncogenic roles in bladder cancer. A potentially useful synthetic biology platform for miRNA loss-of-function study and cancer treatment has been established in this work.