Kerry L. Bubb

Transcriptional outputs of the Caenorhabditis elegans forkhead protein DAF-16

In Caenorhabditis elegans, the forkhead protein DAF‐16 transduces insulin‐like signals that regulate larval development and adult lifespan. To identify DAF‐16‐dependent transcriptional alterations that occur in a long‐lived C. elegans strain, we used cDNA microarrays and genomic analysis to identify putative direct and indirect DAF‐16 transcriptional target genes. Our analysis suggests that DAF‐16 action regulates a wide range of physiological responses by altering the expression of genes involved in metabolism, energy generation and cellular stress responses. Furthermore, we observed a large overlap between DAF‐16‐dependent transcription and genes normally expressed in the long‐lived dauer larval stage. Finally, we examined the in vivo role of 35 of these target genes by RNA‐mediated interference and identified one gene encoding a putative protease that is necessary for the daf‐2 Age phenotype.

Analysis of xbx genes in C. elegans

Cilia and flagella are widespread eukaryotic subcellular components that are conserved from green algae to mammals. In different organisms they function in cell motility, movement of extracellular fluids and sensory reception. While the function and structural description of cilia and flagella are well established, there are many questions that remain unanswered. In particular, very little is known about the developmental mechanisms by which cilia are generated and shaped and how their components are assembled into functional machineries. To find genes involved in cilia development we used as a search tool a promoter motif, the X-box, which participates in the regulation of certain ciliary genes in the nematode Caenorhabditis elegans. By using a genome search approach for X-box promoter motif-containing genes (xbx genes) we identified a list of about 750 xbx genes (candidates). This list comprises some already known ciliary genes as well as new genes, many of which we hypothesize to be important for cilium structure and function. We derived a C. elegans X-box consensus sequence by in vivo expression analysis. We found that xbx gene expression patterns were dependent on particular X-box nucleotide compositions and the distance from the respective gene start. We propose a model where DAF-19, the RFX-type transcription factor binding to the X-box, is responsible for the development of a ciliary module in C. elegans, which includes genes for cilium structure, transport machinery, receptors and other factors.

Mapping and dynamics of regulatory DNA and transcription factor networks in A. thaliana.

Our understanding of gene regulation in plants is constrained by our limited knowledge of plant cis-regulatory DNA and its dynamics. We mapped DNase I hypersensitive sites (DHSs) in A. thaliana seedlings and used genomic footprinting to delineate ∼ 700,000 sites of in vivo transcription factor (TF) occupancy at nucleotide resolution. We show that variation associated with 72 diverse quantitative phenotypes localizes within DHSs. TF footprints encode an extensive cis-regulatory lexicon subject to recent evolutionary pressures, and widespread TF binding within exons may have shaped codon usage patterns. The architecture of A. thaliana TF regulatory networks is strikingly similar to that of animals in spite of diverged regulatory repertoires. We analyzed regulatory landscape dynamics during heat shock and photomorphogenesis, disclosing thousands of environmentally sensitive elements and enabling mapping of key TF regulatory circuits underlying these fundamental responses. Our results provide an extensive resource for the study of A. thaliana gene regulation and functional biology.

Scan of Human Genome Reveals No New Loci Under Ancient Balancing Selection

There has been much speculation as to what role balancing selection has played in evolution. In an attempt to identify regions, such as HLA, at which polymorphism has been maintained in the human population for millions of years, we scanned the human genome for regions of high SNP density. We found 16 regions that, outside of HLA and ABO, are the most highly polymorphic regions yet described; however, evidence for balancing selection at these sites is notably lacking—indeed, whole-genome simulations indicate that our findings are expected under neutrality. We propose that (i) because it is rarely stable, long-term balancing selection is an evolutionary oddity, and (ii) when a balanced polymorphism is ancient in origin, the requirements for detection by means of SNP data alone will rarely be met.

GC-Rich DNA Elements Enable Replication Origin Activity in the Methylotrophic Yeast Pichia pastoris

The well-studied DNA replication origins of the model budding and fission yeasts are A/T-rich elements. However, unlike their yeast counterparts, both plant and metazoan origins are G/C-rich and are associated with transcription start sites. Here we show that an industrially important methylotrophic budding yeast, Pichia pastoris, simultaneously employs at least two types of replication origins—a G/C-rich type associated with transcription start sites and an A/T-rich type more reminiscent of typical budding and fission yeast origins. We used a suite of massively parallel sequencing tools to map and dissect P. pastoris origins comprehensively, to measure their replication dynamics, and to assay the global positioning of nucleosomes across the genome. Our results suggest that some functional overlap exists between promoter sequences and G/C-rich replication origins in P. pastoris and imply an evolutionary bifurcation of the modes of replication initiation.

Profiling of accessible chromatin regions across multiple plant species and cell types reveals common gene regulatory principles and new control modules

A comparison of open chromatin landscapes reveals commonalities in transcriptional regulation across species and identifies a transcription factor cascade in the Arabidopsis root hair. The transcriptional regulatory structure of plant genomes remains poorly defined relative to animals. It is unclear how many cis-regulatory elements exist, where these elements lie relative to promoters, and how these features are conserved across plant species. We employed the assay for transposase-accessible chromatin (ATAC-seq) in four plant species (Arabidopsis thaliana, Medicago truncatula, Solanum lycopersicum, and Oryza sativa) to delineate open chromatin regions and transcription factor (TF) binding sites across each genome. Despite 10-fold variation in intergenic space among species, the majority of open chromatin regions lie within 3 kb upstream of a transcription start site in all species. We find a common set of four TFs that appear to regulate conserved gene sets in the root tips of all four species, suggesting that TF-gene networks are generally conserved. Comparative ATAC-seq profiling of Arabidopsis root hair and non-hair cell types revealed extensive similarity as well as many cell-type-specific differences. Analyzing TF binding sites in differentially accessible regions identified a MYB-driven regulatory module unique to the hair cell, which appears to control both cell fate regulators and abiotic stress responses. Our analyses revealed common regulatory principles among species and shed light on the mechanisms producing cell-type-specific transcriptomes during development.

Detecting Disease-Causing Mutations in the Human Genome by Haplotype Matching

Comparisons between haplotypes from affected patients and the human reference genome are frequently used to identify candidates for disease-causing mutations, even though these alignments are expected to reveal a high level of background neutral polymorphism. This limits the scope of genetic studies to relatively small genomic intervals, because current methods for distinguishing potential causal mutations from neutral variation are inefficient. Here we describe a new strategy for detecting mutations that is based on comparing affected haplotypes with closely matched control sequences from healthy individuals, rather than with the human reference genome. We use theory, simulation, and a real data set to show that this approach is expected to reduce the number of sequence variants that must be subjected to follow-up analysis by at least a factor of 20 when closely matched control sequences are selected from a reference panel with as few as 100 control genomes. We also define a reference data resource that would allow efficient application of this strategy to large critical intervals across the genome.

Regulatory DNA in A. thaliana can tolerate high levels of sequence divergence

Variation in regulatory DNA is thought to drive evolution. Cross-species comparisons of regulatory DNA have provided evidence for both weak purifying selection and substantial turnover in regulatory regions. However, disruption of transcription factor binding sites can affect the expression of neighboring genes. Thus, the base-pair level functional annotation of regulatory DNA has proven challenging. Here, we explore regulatory DNA variation and its functional consequences in genetically diverse strains of the plant Arabidopsis thaliana, which largely maintain the positional homology of regulatory DNA. Using chromatin accessibility to delineate regulatory DNA genome-wide, we find that 15% of approximately 50,000 regulatory sites varied in accessibility among strains. Some of these accessibility differences are associated with extensive underlying sequence variation, encompassing many deletions and dramatically hypervariable sequence. For the majority of such regulatory sites, nearby gene expression was similar, despite this large genetic variation. However, among all regulatory sites, those with both high levels of sequence variation and differential chromatin accessibility are the most likely to reside near genes with differential expression among strains. Unexpectedly, the vast majority of regulatory sites that differed in chromatin accessibility among strains show little variation in the underlying DNA sequence, implicating variation in upstream regulators.

Mapping and dynamics of regulatory DNA during seed development

The genome is reprogrammed during development to produce diverse cell types, largely through altered expression and activity of key transcription factors. The accessibility and critical function of epidermal cells have made them a model for connecting transcriptional events to development in a range of model systems. In Arabidopsis and many other plants, fertilization triggers differentiation of specialized epidermal cells called the seed coat that have a unique morphology caused by large extracellular deposits of pectin. Here, we used DNase I-seq to generate regulatory landscapes of Arabidopsis thaliana seeds at two points in development. To enrich for seed coat signals, we used the INTACT method for capturing nuclei from GL2-expressing cells from whole siliques. Over half of the regulatory (i.e., accessible) bases identified in seeds at these two developmental time points had not been previously identified in seven-day-old seedling. We identified over 3000 regions that were developmentally dynamic (i.e., differentially accessible) at the two timepoints and found that genes neighboring these regions were enriched for gene ontology terms such as development, regulation and pigment. Genes neighboring differentially accessible regions were also significantly more likely to have differing expression at these two timepoints. The differentially accessible regions themselves were enriched for motifs belonging to transcription factors involved in development and environmental response, notably the TCP and MYB families at the earlier and later time points, respectively. These findings are consistent with previous studies, highlighting the dramatic changes that take place in the chromatin landscape during development, and indicating the importance of extending the approach described here to other cell types and developmental processes.

Developmental and conditional dynamics of gene expression in single root cells of A. thaliana

Single-cell RNA-seq can yield high-resolution cell-type-specific expression signatures that reveal new cell types and the developmental trajectories of cell lineages. Here, we apply this approach to A. thaliana root cells to capture gene expression in 3,121 root cells. We analyze these data with Monocle 3, which orders single cell transcriptomes in an unsupervised manner and uses machine learning to reconstruct single-cell developmental trajectories along pseudotime. We identify hundreds of genes with cell-type-specific expression, with pseudotime analysis of several cell lineages revealing both known and novel genes that are expressed along a developmental trajectory. We identify transcription factor motifs that are enriched in early and late cells, together with the corresponding candidate transcription factors that likely drive the observed expression patterns. Finally, by applying heat stress to whole seedlings, we address the longstanding question of possible heterogeneity among cell types in the response to an abiotic stress. Although the response of canonical heat shock genes dominates expression across cell types, subtle but significant differences in other genes can be detected among cell types. Taken together, our results demonstrate that single-cell transcriptomics holds promise for studying plant development and plant physiology with unprecedented resolution.