Search

Isogenic individuals can display seemingly stochastic phenotypic differences, limiting the accuracy of genotype-to-phenotype predictions. The extent of this phenotypic variation depends in part on genetic background, raising questions about the genes involved in controlling stochastic phenotypic variation. Focusing on early seedling traits in Arabidopsisthaliana, we found that hypomorphs of the cuticle-related geneLTP2greatly increased variation in seedling phenotypes, including hypocotyl length, gravitropism and cuticle permeability. Manyltp2hypocotyls were significantly shorter than wild-type hypocotyls while others resembled the wild type. Differences in epidermal properties and gene expression betweenltp2seedlings with long and short hypocotyls suggest a loss of cuticle integrity as the primary determinant of the observed phenotypic variation. We identified environmental conditions that reveal or mask the increased variation inltp2hypomorphs, and found that increased expression of its closest paralogLTP1is necessary forltp2phenotypes. Our results illustrate how decreased expression of a single gene can generate starkly increased phenotypic variation in isogenic individuals in response to an environmental challenge.

Argonaute 1 (AGO1), the principal protein component of microRNA-mediated regulation, plays a key role in plant growth and development. AGO1 physically interacts with the chaperone HSP90, which buffers cryptic genetic variation in plants and animals. We sought to determine whether genetic perturbation of AGO1 in Arabidopsis thaliana would also reveal cryptic genetic variation, and if so, whether AGO1-dependent loci overlap with those dependent on HSP90. To address these questions, we introgressed a hypomorphic mutant allele of AGO1 into a set of mapping lines derived from the commonly used Arabidopsis strains Col-0 and Ler. Although we identified several cases in which AGO1 buffered genetic variation, none of the AGO1-dependent loci overlapped with those buffered by HSP90 for the traits assayed. We focused on one buffered locus where AGO1 perturbation uncoupled the traits days to flowering and rosette leaf number, which are otherwise closely correlated. Using a bulk segregant approach, we identified a non-functional Ler hua2 mutant allele as the causal AGO1-buffered polymorphism. Introduction of a non-functional hua2 allele into a Col-0 ago1 mutant background recapitulated the Ler-dependent ago1 phenotype, implying that coupling of these traits involves different molecular players in these closely related strains. Taken together, our findings demonstrate that even though AGO1 and HSP90 buffer genetic variation in the same traits, these robustness regulators interact epistatically with different genetic loci, suggesting that higher-order epistasis is uncommon.Article SummaryArgonaute 1 (AGO1), a key player in plant development, interacts with the chaperone HSP90 which buffers environmental and genetic variation. We found that AGO1 buffers environmental and genetic variation in the same traits; however, AGO1-dependent and HSP90-dependent loci do not overlap. Detailed analysis of a buffered locus found that a non-functional HUA2 allele decouples days to flowering and rosette leaf number in an AGO1-dependent manner, suggesting the AGO1-dependent buffering acts at the network level.

Transcriptional control is exerted primarily through the binding of transcription factor proteins to regulatory elements in DNA. By virtue of eukaryotic DNA being complexed with histones, transcription factor binding to DNA alters or eliminates histone-DNA contacts, leading to increased accessibility of the DNA region to nuclease enzymes. This hypersensitivity to nuclease digestion has been used to define DNA binding events and regulatory elements across genomes, and to compare these attributes between cell types or conditions. These approaches make it possible to define the regulatory elements in a genome as well as to predict the regulatory networks of transcription factors and their target genes in a given cell state. As these chromatin accessibility assays are increasingly used, it is important to consider how to analyze the resulting data to avoid artifactual results or misinterpretation. In this review, we focus on some of the key technical and computational caveats associated with plant chromatin accessibility data, including strategies for sample preparation, sequencing, read mapping, and downstream analyses.

Here, Zhou et al. took advantage of Saccharomyces cerevisiae and its well-characterized mating and invasion pathways to explain how transcription factor variants alter motif recognition and gene expression, and ultimately, organismal phenotypes..., Amino acid substitutions are commonly found in human transcription factors, yet the functional consequences of much of this variation remain unknown, even in well-characterized DNA-binding domains. Here, we examine how six single-amino acid variants in the DNA-binding domain of Ste12—a yeast transcription factor regulating mating and invasion—alter Ste12 genome binding, motif recognition, and gene expression to yield markedly different phenotypes. Using a combination of the “calling-card” method, RNA sequencing, and HT-SELEX (high throughput systematic evolution of ligands by exponential enrichment), we find that variants with dissimilar binding and expression profiles can converge onto similar cellular behaviors. Mating-defective variants led to decreased expression of distinct subsets of genes necessary for mating. Hyper-invasive variants also decreased expression of subsets of genes involved in mating, but increased the expression of other subsets of genes associated with the cellular response to osmotic stress. While single-amino acid changes in the coding region of this transcription factor result in complex regulatory reconfiguration, the major phenotypic consequences for the cell appear to depend on changes in the expression of a small number of genes with related functions.

Single cell transcriptomic profiling can map developmental trajectories and assess cell-type–specific changes during heat shock at single cell resolution., Single cell RNA sequencing 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 Arabidopsis (Arabidopsis 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. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. 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.

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 functions of epidermal cells have made them a model for connecting transcriptional events to development in a range of model systems. In Arabidopsis thaliana and many other plants, fertilization triggers differentiation of specialized epidermal seed coat cells that have a unique morphology caused by large extracellular deposits of polysaccharides. Here, we used DNase I-seq to generate regulatory landscapes of A. thaliana seeds at two critical time points in seed coat maturation (4 and 7 DPA), enriching for seed coat cells with the INTACT method. We found over 3,000 developmentally dynamic regulatory DNA elements and explored their relationship with nearby gene expression. The dynamic regulatory elements were enriched for motifs for several transcription factors families; most notably the TCP family at the earlier time point and the MYB family at the later one. To assess the extent to which the observed regulatory sites in seeds added to previously known regulatory sites in A. thaliana, we compared our data to 11 other data sets generated with 7-day-old seedlings for diverse tissues and conditions. Surprisingly, over a quarter of the regulatory, i.e. accessible, bases observed in seeds were novel. Notably, plant regulatory landscapes from different tissues, cell types, or developmental stages were more dynamic than those generated from bulk tissue in response to environmental perturbations, highlighting the importance of extending studies of regulatory DNA to single tissues and cell types during development.

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 toA. thalianaroot 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. We assess and interpret changes in total RNA expression along developmental trajectories and show that trajectory branch points mark developmental decisions. 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.

Robustness to both genetic and environmental change is an emergent feature of living systems. Loss of phenotypic robustness can be associated with increased penetrance of genetic variation. In model organisms and in humans, the phenotypic consequences of standing genetic variation can be buffered by the molecular chaperone HSP90. However, it has been argued that HSP90 has the opposite effect on newly introduced genetic variation. To test the buffering effect of HSP90 on new mutations, we introduced vast numbers of mutations into wild-type and HSP90-reduced plants and assessed embryonic lethality and early seedling phenotypes for thousands of offspring. Although the levels of newly introduced mutations were similar in the two backgrounds, the HSP90-reduced plants showed a significantly greater frequency of embryonic lethality and severe phenotypic abnormalities, consistent with higher penetrance and expressivity of newly introduced genetic variation. We further demonstrate that some mutant phenotypes were heritable in an HSP90-dependent manner, and we map candidate HSP90-dependent polymorphisms. Moreover, both sequence and phenotypic analyses of wild-type and HSP90-reduced plants suggest that the HSP90-dependent phenotypes are largely due the newly introduced mutations rather than to an increased mutation rate in HSP90-reduced plants. Taken together, our results support a model in which HSP90 buffers newly introduced mutations, and the phenotypic consequences of such mutations outweigh those of mutations arising de novo in response to HSP90 perturbation.

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 functions of epidermal cells have made them a model for connecting transcriptional events to development in a range of model systems. In Arabidopsis thaliana and many other plants, fertilization triggers differentiation of specialized epidermal seed coat cells that have a unique morphology caused by large extracellular deposits of pectin. Here, we used DNase I-seq to generate regulatory landscapes of A. thaliana seeds at two critical time points in seed coat maturation, enriching for seed coat cells with the INTACT method. We found over 3000 developmentally dynamic regulatory DNA elements and explored their relationship with nearby gene expression. The dynamic regulatory elements were enriched for motifs for several transcription factors families; most notably the TCP family at the earlier time point and the MYB family at the later one. To assess the extent to which the observed regulatory sites in seeds added to previously known regulatory sites in A. thaliana, we compared our data to 11 other data sets generated with seven-day-old seedlings for diverse tissues and conditions. Surprisingly, over a quarter of the regulatory, i.e. accessible, bases observed in seeds were novel. Notably, in this comparison, development exerted a stronger effect on the plant regulatory landscape than extreme environmental perturbations, highlighting the importance of extending studies of regulatory landscapes to other tissues and cell types during development.

Understanding gene regulatory networks in plants requires knowledge of cis- regulatory DNA, trans -acting factors, and their dynamics across development and in response to stimuli. Active cis -regulatory elements are hypersensitive to cleavage by the endonuclease DNase I. Motifs within DNase I hypersensitive sites indicate potential trans -acting factor occupancy and, when combined with DNase I cleavage data, can be used to construct provisional regulatory networks. Several recent studies have applied genome-wide DNase I hypersensitivity mapping to Arabidopsis thaliana and rice, generating chromatin accessibility landscapes for an array of tissues, cell types, and treatments. Here, we discuss these studies, with an emphasis on building regulatory networks and possible directions for improvements.

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 plantArabidopsis 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.

Dynamic Imprinting Gene imprinting—the silencing of either a maternally derived or paternally derived gene allele—is controlled in large part by DNA methylation. In plants, imprinting occurs in the endosperm, which nourishes the embryonic plant. Gehring et al. (p. 1447 ) and Hsieh et al. (p. 1451 ) analyzed the dynamics of DNA methylation in the endosperm and embryo of Arabidopsis and found extensive demethylation in the endosperm, suggesting that many imprinted genes are likely to exist. Gehring et al. characterized five imprinted genes in detail. Four of the 10 known imprinted genes are related homeodomain transcription factors. Furthermore, 5′ sequences demethylated in several of the genes were found to be derived from transposable elements, which supports the idea that imprinting arose as a by-product of silencing invading DNA.

Despite decade-long efforts, the genetic underpinnings of many complex traits and diseases remain largely elusive. It is increasingly recognized that a purely additive model, upon which most genome-wide association studies (GWAS) rely, is insufficient. Although thousands of significant trait-associated loci have been identified, purely additive models leave much of the inferred genetic variance unexplained. Several factors have been invoked to explain the ‘missing heritability’, including epistasis. Accounting for all possible epistatic interactions is computationally complex and requires very large samples. Here, we propose a simple two-state epistasis model, in which individuals show either high or low variant penetrance with respect to a certain trait. The use of this model increases the power to detect additive trait-associated loci. We show that this model is consistent with current GWAS results and improves fit with heritability observations based on twin studies. We suggest that accounting for variant penetrance will significantly increase our power to identify underlying additive loci.

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.

We define ‘genetic individuality’ as intraspecies variation that has substantial heritability and involves traits that are sufficiently common that they can be observed in any modest–sized sampling of individuals. We propose that genetic individuality is largely shaped by the combinatory shuffling of a modest number of genes, each of which exists as a family of functionally and structurally diverged alleles. Unequivocal examples of such allele families are found at the O–antigen–biosynthetic locus inPseudomonas aeruginosaand the human leucocyte antigen locus in humans. We examine characteristic features of these allele families and explore the possibility that genetic loci with similar characteristics can be recognized in a whole–genome scan of human genetic variation.

Human chromosome 7 has historically received prominent attention in the human genetics community, primarily related to the search for the cystic fibrosis gene and the frequent cytogenetic changes associated with various forms of cancer. Here we present more than 153 million base pairs representing 99.4% of the euchromatic sequence of chromosome 7, the first metacentric chromosome completed so far. The sequence has excellent concordance with previously established physical and genetic maps, and it exhibits an unusual amount of segmentally duplicated sequence (8.2%), with marked differences between the two arms. Our initial analyses have identified 1,150 protein-coding genes, 605 of which have been confirmed by complementary DNA sequences, and an additional 941 pseudogenes. Of genes confirmed by transcript sequences, some are polymorphic for mutations that disrupt the reading frame.

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.

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.

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., Author Summary Genome duplication in eukaryotes initiates at loci called replication origins. Origins in most budding and fission yeasts are A/T-rich DNA sequences, while metazoan origins are G/C-rich and are often associated with promoters. Here we have globally mapped replication origins and nucleosome positions in an industrially important methylotrophic yeast, Pichia pastoris. We show that P. pastoris has two general classes of origins—A/T-rich origins resembling those of most other yeasts, and a novel, G/C-rich class, that appear more robust and are associated with promoters. P. pastoris is the first known species using two kinds of origins and the first known budding yeast to use a G/C-rich origin motif. Additionally, the G/C-rich motif matches one of the motifs annotated as binding sites of the human Hsf1 transcriptional regulator suggesting that in this species there may be a link between transcriptional regulation and DNA replication initiation.

Summary 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.

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.