Kelvin X. Zhang

DNA demethylation dynamics in the human prenatal germline

Global DNA demethylation in humans is a fundamental process that occurs in pre-implantation embryos and reversion to naive ground state pluripotent stem cells (PSCs). However, the extent of DNA methylation reprogramming in human germline cells is unknown. Here, we performed whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq) of human prenatal germline cells from 53 to 137 days of development. We discovered that the transcriptome and methylome of human germline is distinct from both human PSCs and the inner cell mass (ICM) of human blastocysts. Using this resource to monitor the outcome of global DNA demethylation with reversion of primed PSCs to the naive ground state, we uncovered hotspots of ultralow methylation at transposons that are protected from demethylation in the germline and ICM. Taken together, the human germline serves as a valuable in vivo tool for monitoring the epigenome of cells that have emerged from a global DNA demethylation event.

The ontogeny of cKIT + human primordial germ cells proves to be a resource for human germ line reprogramming, imprint erasure and in vitro differentiation

The generation of research-quality, clinically relevant cell types in vitro from human pluripotent stem cells requires a detailed understanding of the equivalent human cell types. Here we analysed 134 human embryonic and fetal samples from 6 to 20 developmental weeks and identified the stages at which cKIT+ primordial germ cells (PGCs), the precursors of gametes, undergo whole-genome epigenetic reprogramming with global depletion of 5mC, H3K27me3 and H2A.Z, and the time at which imprint erasure is initiated and 5hmC is present. Using five alternative in vitro differentiation strategies combined with single-cell microfluidic analysis and a bona fide human cKIT+ PGC signature, we show the stage of cKIT+ PGC formation in the first 16 days of differentiation. Taken together, our study creates a resource of human germ line ontogeny that is essential for future studies aimed at in vitro differentiation and unveiling the mechanisms necessary to pass human DNA from one generation to the next.

Novel stem/progenitor cell population from murine tracheal submucosal gland ducts with multipotent regenerative potential.

The airway epithelium is in direct contact with the environment and therefore constantly at risk for injury. Basal cells (BCs) have been found to repair the surface epithelium (SE), but the contribution of other stem cell populations to airway epithelial repair has not been identified. We demonstrated that airway submucosal gland (SMG) duct cells, in addition to BCs, survived severe hypoxic‐ischemic injury. We developed a method to isolate duct cells from the airway. In vitro and in vivo models were used to compare the self‐renewal and differentiation potential of duct cells and BCs. We found that only duct cells were capable of regenerating SMG tubules and ducts, as well as the SE overlying the SMGs. SMG duct cells are therefore a multipotent stem cell for airway epithelial repair This is of importance to the field of lung regeneration as determining the repairing cell populations could lead to the identification of novel therapeutic targets and cell‐based therapies for patients with airway diseases. STEM CELLS 2011;29:1283–1293

Ig Superfamily Ligand and Receptor Pairs Expressed in Synaptic Partners in Drosophila.

Information processing relies on precise patterns of synapses between neurons. The cellular recognition mechanisms regulating this specificity are poorly understood. In the medulla of the Drosophila visual system, different neurons form synaptic connections in different layers. Here, we sought to identify candidate cell recognition molecules underlying this specificity. Using RNA sequencing (RNA-seq), we show that neurons with different synaptic specificities express unique combinations of mRNAs encoding hundreds of cell surface and secreted proteins. Using RNA-seq and protein tagging, we demonstrate that 21 paralogs of the Dpr family, a subclass of immunoglobulin (Ig)-domain containing proteins, are expressed in unique combinations in homologous neurons with different layer-specific synaptic connections. Dpr interacting proteins (DIPs), comprising nine paralogs of another subclass of Ig-containing proteins, are expressed in a complementary layer-specific fashion in a subset of synaptic partners. We propose that pairs of Dpr/DIP paralogs contribute to layer-specific patterns of synaptic connectivity.

Molecular Profiling of Premalignant Lesions in Lung Squamous Cell Carcinomas Identifies Mechanisms Involved in Stepwise Carcinogenesis

Lung squamous cell carcinoma (SCC) is thought to arise from premalignant lesions in the airway epithelium; therefore, studying these lesions is critical for understanding lung carcinogenesis. Previous microarray and sequencing studies designed to discover early biomarkers and therapeutic targets for lung SCC had limited success identifying key driver events in lung carcinogenesis, mostly due to the cellular heterogeneity of patient samples examined and the interindividual variability associated with difficult to obtain airway premalignant lesions and appropriate normal control samples within the same patient. We performed RNA sequencing on laser-microdissected representative cell populations along the SCC pathologic continuum of patient-matched normal basal cells, premalignant lesions, and tumor cells. We discovered transcriptomic changes and identified genomic pathways altered with initiation and progression of SCC within individual patients. We used immunofluorescent staining to confirm gene expression changes in premalignant lesions and tumor cells, including increased expression of SLC2A1, CEACAM5, and PTBP3 at the protein level and increased activation of MYC via nuclear translocation. Cytoband enrichment analysis revealed coordinated loss and gain of expression in chromosome 3p and 3q regions, respectively, during carcinogenesis. This is the first gene expression profiling study of airway premalignant lesions with patient-matched SCC tumor samples. Our results provide much needed information about the biology of premalignant lesions and the molecular changes that occur during stepwise carcinogenesis of SCC, and it highlights a novel approach for identifying some of the earliest molecular changes associated with initiation and progression of lung carcinogenesis within individual patients. Cancer Prev Res; 7(5); 487–95. ©2014 AACR.

Genome-wide analysis of dorsal and ventral transcriptomes of the Xenopus laevis gastrula

RNA sequencing has allowed high-throughput screening of differential gene expression in many tissues and organisms. Xenopus laevis is a classical embryological and cell-free extract model system, but its genomic sequence had been lacking due to difficulties arising from allotetraploidy. There is currently much excitement surrounding the release of the completed X. laevis genome (version 9.1) by the Joint Genome Institute (JGI), which provides a platform for genome-wide studies. Here we present a deep RNA-seq dataset of transcripts expressed in dorsal and ventral lips of the early Xenopus gastrula embryo using the new genomic information, which was further annotated by blast searches against the human proteome. Overall, our findings confirm previous results from differential screenings using other methods that uncovered classical dorsal genes such as Chordin, Noggin and Cerberus, as well as ventral genes such as Sizzled, Ventx, Wnt8 and Bambi. Complete transcriptome-wide tables of mRNAs suitable for data mining are presented, which include many novel dorsal- and ventral-specific genes. RNA-seq was very quantitative and reproducible, and allowed us to define dorsal and ventral signatures useful for gene set expression analyses (GSEA). As an example of a new gene, we present here data on an organizer-specific secreted protein tyrosine kinase known as Pkdcc (protein kinase domain containing, cytoplasmic) or Vlk (vertebrate lonesome kinase). Overexpression experiments indicate that Pkdcc can act as a negative regulator of Wnt/ β-catenin signaling independently of its kinase activity. We conclude that RNA-Seq in combination with the X. laevis complete genome now available provides a powerful tool for unraveling cell-cell signaling pathways during embryonic induction.

Activation of Notch1 synergizes with multiple pathways in promoting castration-resistant prostate cancer

Significance A high nuclear Notch homolog 1, translocation-associated (Notch1) intracellular domain level distinguishes high-risk prostate cancer and castration-resistant prostate cancer from benign and low/intermediate-risk prostate cancer. Chronic activation of Notch1 cooperates with multiple oncogenic pathways altered in early prostate cancer, including AKT, Myc, and Ras/Raf/MAPK, to promote progression to androgen ablation-resistant prostate adenocarcinoma. Metastatic castration-resistant prostate cancer (CRPC) is the primary cause of prostate cancer-specific mortality. Defining new mechanisms that can predict recurrence and drive lethal CRPC is critical. Here, we demonstrate that localized high-risk prostate cancer and metastatic CRPC, but not benign prostate tissues or low/intermediate-risk prostate cancer, express high levels of nuclear Notch homolog 1, translocation-associated (Notch1) receptor intracellular domain. Chronic activation of Notch1 synergizes with multiple oncogenic pathways altered in early disease to promote the development of prostate adenocarcinoma. These tumors display features of epithelial-to-mesenchymal transition, a cellular state associated with increased tumor aggressiveness. Consistent with its activation in clinical CRPC, tumors driven by Notch1 intracellular domain in combination with multiple pathways altered in prostate cancer are metastatic and resistant to androgen deprivation. Our study provides functional evidence that the Notch1 signaling axis synergizes with alternative pathways in promoting metastatic CRPC and may represent a new therapeutic target for advanced prostate cancer.

Aldehyde Dehydrogenase Activity Enriches for Proximal Airway Basal Stem Cells and Promotes Their Proliferation

Both basal and submucosal gland (SMG) duct stem cells of the airway epithelium are capable of sphere formation in the in vitro sphere assay, although the efficiency at which this occurs is very low. We sought to improve this efficiency of sphere formation by identifying subpopulations of airway basal stem cells (ABSC) and SMG duct cells based on their aldehyde dehydrogenase (ALDH) activity. ALDH(hi) ABSCs and SMG duct cells were highly enriched for the population of cells that could make spheres, while the co-culture of ALDH(hi) differentiated cells with the ALDH(hi) ABSCs increased their sphere-forming efficiency. Specific ALDH agonists and antagonists were used to show that airway specific ALDH isozymes are important for ABSC proliferation. Pathway analysis of gene expression profiling of ALDH(hi) and ALDH(lo) ABSCs revealed a significant upregulation of the arachidonic acid (AA) metabolism pathway in ALDH(hi) ABSCs. We confirmed the importance of this pathway in the metabolism of proliferating ALDH(hi) ABSCs using bioenergetics studies as well as agonists and antagonists of the AA pathway. These studies could lead to the development of novel strategies for altering ABSC proliferation in the airway epithelium.

Spemann organizer transcriptome induction by early beta-catenin, Wnt, Nodal, and Siamois signals in Xenopus laevis

Significance We present a genome-wide study of the signals responsible for the early induction of the body axis in the following experimental conditions: β-catenin morpholino; Wnt, Siamois, and Cerberus mRNAs; LiCl treatment; and dorsal-ventral regenerating half-embryos bisected at gastrula. Comparing 46 RNA-seq libraries, we uncovered the genetic networks that initiate dorsal-ventral patterning and Spemann’s organizer formation. We defined an early β-catenin signature that has only minor overlap with recently published late zygotic Wnt signatures. The relation of these early steps of development to endomesodermal germ layer induction was studied by overexpressing the growth factor antagonist Cerberus. This study offers a rich resource for understanding the earliest inductive events in the body plan of a model vertebrate embryo. The earliest event in Xenopus development is the dorsal accumulation of nuclear β-catenin under the influence of cytoplasmic determinants displaced by fertilization. In this study, a genome-wide approach was used to examine transcription of the 43,673 genes annotated in the Xenopus laevis genome under a variety of conditions that inhibit or promote formation of the Spemann organizer signaling center. Loss of function of β-catenin with antisense morpholinos reproducibly reduced the expression of 247 mRNAs at gastrula stage. Interestingly, only 123 β-catenin targets were enriched on the dorsal side and defined an early dorsal β-catenin gene signature. These genes included several previously unrecognized Spemann organizer components. Surprisingly, only 3 of these 123 genes overlapped with the late Wnt signature recently defined by two other groups using inhibition by Dkk1 mRNA or Wnt8 morpholinos, which indicates that the effects of β-catenin/Wnt signaling in early development are exquisitely regulated by stage-dependent mechanisms. We analyzed transcriptome responses to a number of treatments in a total of 46 RNA-seq libraries. These treatments included, in addition to β-catenin depletion, regenerating dorsal and ventral half-embryos, lithium chloride treatment, and the overexpression of Wnt8, Siamois, and Cerberus mRNAs. Only some of the early dorsal β-catenin signature genes were activated at blastula whereas others required the induction of endomesoderm, as indicated by their inhibition by Cerberus overexpression. These comprehensive data provide a rich resource for analyzing how the dorsal and ventral regions of the embryo communicate with each other in a self-organizing vertebrate model embryo.

FamAnn: an automated variant annotation pipeline to facilitate target discovery for family-based sequencing studies

FamAnn is an automated variant annotation pipeline designed for facilitating target discovery for family-based sequencing studies. It can apply a different inheritance pattern or a de novo mutations discovery model to each family and select single nucleotide variants and small insertions and deletions segregating in each family or shared by multiple families. It also provides a variety of variant annotations and retains and annotates all transcripts hit by a single variant. Excel-compatible outputs including all annotated variants segregating in each family or shared by multiple families will be provided for users to prioritize variants based on their customized thresholds. A list of genes that harbor the segregating variants will be provided as well for possible pathway/network analyses. FamAnn uses the de facto community standard Variant Call Format as the input format and can be applied to whole exome, genome or targeted resequencing data. AVAILABILITY https://sites.google.com/site/famannotation/home CONTACT: jianchaoyao@gmail.com, kelvinzhang@mednet.ucla.edu, mccombie@cshl.edu Supplementary information: Supplementary data are available at Bioinformatics online.