Joshua T. Burdick

Mapping determinants of human gene expression by regional and genome-wide association

To study the genetic basis of natural variation in gene expression, we previously carried out genome-wide linkage analysis and mapped the determinants of ∼1,000 expression phenotypes. In the present study, we carried out association analysis with dense sets of single-nucleotide polymorphism (SNP) markers from the International HapMap Project. For 374 phenotypes, the association study was performed with markers only from regions with strong linkage evidence; these regions all mapped close to the expressed gene. For a subset of 27 phenotypes, analysis of genome-wide association was performed with >770,000 markers. The association analysis with markers under the linkage peaks confirmed the linkage results and narrowed the candidate regulatory regions for many phenotypes with strong linkage evidence. The genome-wide association analysis yielded highly significant results that point to the same locations as the genome scans for about 50% of the phenotypes. For one candidate determinant, we carried out functional analyses and confirmed the variation in cis-acting regulatory activity. Our findings suggest that association studies with dense SNP maps will identify susceptibility loci or other determinants for some complex traits or diseases.

Common genetic variants account for differences in gene expression among ethnic groups

Variation in DNA sequence contributes to individual differences in quantitative traits, but in humans the specific sequence variants are known for very few traits. We characterized variation in gene expression in cells from individuals belonging to three major population groups. This quantitative phenotype differs significantly between European-derived and Asian-derived populations for 1,097 of 4,197 genes tested. For the phenotypes with the strongest evidence of cis determinants, most of the variation is due to allele frequency differences at cis-linked regulators. The results show that specific genetic variation among populations contributes appreciably to differences in gene expression phenotypes. Populations differ in prevalence of many complex genetic diseases, such as diabetes and cardiovascular disease. As some of these are probably influenced by the level of gene expression, our results suggest that allele frequency differences at regulatory polymorphisms also account for some population differences in prevalence of complex diseases.

Monozygotic Twins Reveal Germline Contribution to Allelic Expression Differences

Variation in the level of gene expression is a major determinant of a cells function and characteristics. Common allelic variants of genes can be expressed at different levels and thus contribute to phenotypic diversity. We have measured allelic expression differences at heterozygous loci in monozygotic twins and in unrelated individuals. We show that the extent of differential allelic expression is highly similar within monozygotic twin pairs for many loci, implying that allelic differences in gene expression are under genetic control. We also show that even subtle departures from equal allelic expression are often genetically determined.

A quantitative model of normal Caenorhabditis elegans embryogenesis and its disruption after stress.

The invariant lineage of Caenorhabditis elegans has powerful potential for quantifying developmental variability in normal and stressed embryos. Previous studies of division timing by automated lineage tracing suggested that variability in cell cycle timing is low in younger embryos, but manual lineage tracing of specific lineages suggested that variability may increase for later divisions. We developed improved automated lineage tracing methods that allowroutine lineage tracing through the last round of embryonic cell divisions and we applied these methods to trace the lineage of 18 wild-type embryos. Cell cycle lengths, division axes and cell positions are remarkably consistent among these embryos at all stages, with only slight increase in variability later in development. The resulting quantitative 4-dimensional model of embryogenesis provides a powerful reference dataset to identify defects in mutants or in embryos that have experienced environmental perturbations. We also traced the lineages of embryos imaged at higher temperatures to quantify the decay in developmental robustness under temperature stress. Developmental variability increases modestly at 25°C compared with 22°C and dramatically at 26°C, and we identify homeotic transformations in a subset of embryos grown at 26°C. The deep lineage tracing methods provide a powerful tool for analysis of normal development, gene expression and mutants and we provide a graphical user interface to allow other researchers to explore the average behavior of arbitrary cells in a reference embryo.

GenMapDB: a database of mapped human BAC clones

GenMapDB (http://genomics.med.upenn.edu/genmapdb) is a repository of human bacterial artificial chromosome (BAC) clones mapped by our laboratory to sequence-tagged site markers. Currently, GenMapDB contains over 3000 mapped clones that span 19 chromosomes, chromosomes 2, 4, 5, 9-22, X and Y. This database provides positional information about human BAC clones from the RPCI-11 human male BAC library. It also contains restriction fragment analysis data and end sequences of the clones. GenMapDB is freely available to the public. The main purpose of GenMapDB is to organize the mapping data and to allow the research community to search for mapped BAC clones that can be used in gene mapping studies and chromosomal mutation analysis projects.

Senataxin Mutation Reveals How R-Loops Promote Transcription by Blocking DNA Methylation at Gene Promoters

R-loops are three-stranded nucleic acid structures found abundantly and yet often viewed as by-products of transcription. Studying cells from patients with a motor neuron disease (amyotrophic lateral sclerosis 4 [ALS4]) caused by a mutation in senataxin, we uncovered how R-loops promote transcription. In ALS4 patients, the senataxin mutation depletes R-loops with a consequent effect on gene expression. With fewer R-loops in ALS4 cells, the expression of BAMBI, a negative regulator of transforming growth factor β (TGF-β), is reduced; that then leads to the activation of the TGF-β pathway. We uncovered that genome-wide R-loops influence promoter methylation of over 1,200 human genes. DNA methyl-transferase 1 favors binding to double-stranded DNA over R-loops. Thus, in forming R-loops, nascent RNA blocks DNA methylation and promotes further transcription. Hence, our results show that nucleic acid structures, in addition to sequences, influence the binding and activity of regulatory proteins.

The Lipocalin LPR-1 Cooperates with LIN-3/EGF Signaling To Maintain Narrow Tube Integrity in Caenorhabditis elegans

Lipocalins are secreted cup-shaped glycoproteins that bind sterols, fatty acids, and other lipophilic molecules. Lipocalins have been implicated in a wide array of processes related to lipophilic cargo transport, sequestration, and signaling, and several are used as biomarkers for human disease, but the functions of most lipocalins remain poorly understood. Here we show that the Caenorhabditis elegans lipocalin LPR-1 is required to maintain apical membrane integrity and a continuous lumen in two narrow unicellular tubes, the excretory duct and pore, during a period of rapid lumen elongation. LPR-1 fusion protein is expressed by the duct and pore and accumulates both intracellularly and in apical extracellular compartments, but it can also function cell nonautonomously when provided from outside of the excretory system. lpr-1 mutant defects can be rescued by increased signaling through the epidermal growth factor (EGF)-Ras-extracellular signal regulated kinase (ERK) pathway, which promotes the more elongated duct vs. less elongated pore tube fate. Spatial and temporal rescue experiments indicate that Ras signaling acts within the duct and pore tubes during or prior to cell fate determination to bypass the requirement for LPR-1. lpr-1 mutations did not disrupt LIN-3/EGF-dependent duct-fate specification, prevent functioning of any specific LIN-3/EGF isoform, or alter LET-23/EGFR localization, and reduced signaling did not phenocopy or enhance lpr-1 mutant defects. These data suggest that LPR-1 protects lumen integrity through a LIN-3/EGF-independent mechanism, but that increased signaling upregulates some target(s) that can compensate for lpr-1 absence.

Overlapping cell population expression profiling and regulatory inference in C. elegans

BackgroundUnderstanding gene expression across the diverse metazoan cell types during development is critical to understanding their function and regulation. However, most cell types have not been assayed for expression genome-wide.ResultsWe applied a novel approach we term “Profiling of Overlapping Populations of cells (POP-Seq)” to assay differential expression across all embryonic cells in the nematode Caenorhabditis elegans. In this approach, we use RNA-seq to define the transcriptome of diverse partially overlapping FACS-sorted cell populations. This identified thousands of transcripts differentially expressed across embryonic cells. Hierarchical clustering analysis identified over 100 sets of coexpressed genes corresponding to distinct patterns of cell type specific expression. We identified thousands of candidate regulators of these clusters based on enrichment of transcription factor motifs and experimentally determined binding sites.ConclusionsOur analysis provides new insight into embryonic gene regulation, and provides a resource for improving our knowledge of tissue-specific expression and its regulation throughout C. elegans development.

Deconvolution of gene expression from cell populations across the C. elegans lineage

BackgroundKnowledge of when and in which cells each gene is expressed across multicellular organisms is critical in understanding both gene function and regulation of cell type diversity. However, methods for measuring expression typically involve a trade-off between imaging-based methods, which give the precise location of a limited number of genes, and higher throughput methods such as RNA-seq, which include all genes, but are more limited in their resolution to apply to many tissues. We propose an intermediate method, which estimates expression in individual cells, based on high-throughput measurements of expression from multiple overlapping groups of cells. This approach has particular benefits in organisms such as C. elegans where invariant developmental patterns make it possible to define these overlapping populations of cells at single-cell resolution.ResultWe implement several methods to deconvolve the gene expression in individual cells from population-level data and determine the accuracy of these estimates on simulated data from the C. elegans embryo.ConclusionThese simulations suggest that a high-resolution map of expression in the C. elegans embryo may be possible with expression data from as few as 30 cell populations.

Human proteins that interact with RNA/DNA hybrids

RNA/DNA hybrids form when RNA hybridizes with its template DNA generating a three-stranded structure known as the R-loop. Knowledge of how they form and resolve, as well as their functional roles, is limited. Here, by pull-down assays followed by mass spectrometry, we identified 803 proteins that bind to RNA/DNA hybrids. Because these proteins were identified using in vitro assays, we confirmed that they bind to R-loops in vivo. They include proteins that are involved in a variety of functions, including most steps of RNA processing. The proteins are enriched for K homology (KH) and helicase domains. Among them, more than 300 proteins preferred binding to hybrids than double-stranded DNA. These proteins serve as starting points for mechanistic studies to elucidate what RNA/DNA hybrids regulate and how they are regulated.