The completion of the human genome sequence has made possible genome-wide studies of retroviral DNA integration. Here we report an analysis of 3,127 integration site sequences from human cells. We compared retroviral vectors derived from human immunodeficiency virus (HIV), avian sarcoma-leukosis virus (ASLV), and murine leukemia virus (MLV). Effects of gene activity on integration targeting were assessed by transcriptional profiling of infected cells. Integration by HIV vectors, analyzed in two primary cell types and several cell lines, strongly favored active genes. An analysis of the effects of tissue-specific transcription showed that it resulted in tissue-specific integration targeting by HIV, though the effect was quantitatively modest. Chromosomal regions rich in expressed genes were favored for HIV integration, but these regions were found to be interleaved with unfavorable regions at CpG islands. MLV vectors showed a strong bias in favor of integration near transcription start sites, as reported previously. ASLV vectors showed only a weak preference for active genes and no preference for transcription start regions. Thus, each of the three retroviruses studied showed unique integration site preferences, suggesting that virus-specific binding of integration complexes to chromatin features likely guides site selection.

How growth regulators provoke context-specific signals is a fundamental question in developmental biology. In plants, both auxin and brassinosteroids (BRs) promote cell expansion, and it was thought that they activated this process through independent mechanisms. In this work, we describe a shared auxin:BR pathway required for seedling growth. Genetic, physiological, and genomic analyses demonstrate that response from one pathway requires the function of the other, and that this interdependence does not act at the level of hormone biosynthetic control. Increased auxin levels saturate the BR-stimulated growth response and greatly reduce BR effects on gene expression. Integration of these two pathways is downstream from BES1 and Aux/IAA proteins, the last known regulatory factors acting downstream of each hormone, and is likely to occur directly on the promoters of auxin:BR target genes. We have developed a new approach to identify potential regulatory elements acting in each hormone pathway, as well as in the shared auxin:BR pathway. We show that one element highly overrepresented in the promoters of auxin- and BR-induced genes is responsive to both hormones and requires BR biosynthesis for normal expression. This work fundamentally alters our view of BR and auxin signaling and describes a powerful new approach to identify regulatory elements required for response to specific stimuli.

Organisms are able to anticipate changes in the daily environment with an internal oscillator know as the circadian clock. Transcription is an important mechanism in maintaining these oscillations. Here we explore, using whole genome tiling arrays, the extent of rhythmic expression patterns genome-wide, with an unbiased analysis of coding and noncoding regions of the Arabidopsis genome.

Although pioneered by human geneticists as a potential solution to the challenging problem of finding the genetic basis of common human diseases, genome-wide association (GWA) studies have, owing to advances in genotyping and sequencing technology, become an obvious general approach for studying the genetics of natural variation and traits of agricultural importance. They are particularly useful when inbred lines are available, because once these lines have been genotyped they can be phenotyped multiple times, making it possible (as well as extremely cost effective) to study many different traits in many different environments, while replicating the phenotypic measurements to reduce environmental noise. Here we demonstrate the power of this approach by carrying out a GWA study of 107 phenotypes in Arabidopsis thaliana, a widely distributed, predominantly self-fertilizing model plant known to harbour considerable genetic variation for many adaptively important traits. Our results are dramatically different from those of human GWA studies, in that we identify many common alleles of major effect, but they are also, in many cases, harder to interpret because confounding by complex genetics and population structure make it difficult to distinguish true associations from false. However, a-priori candidates are significantly over-represented among these associations as well, making many of them excellent candidates for follow-up experiments. Our study demonstrates the feasibility of GWA studies in A. thaliana and suggests that the approach will be appropriate for many other organisms.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

Plant architecture is optimized for the local light environment. In response to foliar shade or neighbor proximity (low red to far-red light), some plant species exhibit shade-avoiding phenotypes, including increased stem and hypocotyl growth, which increases the likelihood of outgrowing competitor plants. If shade persists, early flowering and the reallocation of growth resources to stem elongation ultimately affect the yield of harvestable tissues in crop species. Previous studies have shown that hypocotyl growth in low red to far-red shade is largely dependent on the photoreceptor phytochrome B and the phytohormone auxin. However, where shade is perceived in the plant and how auxin regulates growth spatially are less well understood. Using the oilseed and vegetable crop species Brassica rapa, we show that the perception of low red to far-red shade by the cotyledons triggers hypocotyl cell elongation and auxin target gene expression. Furthermore, we find that following shade perception, elevated auxin levels occur in a basipetal gradient away from the cotyledons and that this is coincident with a gradient of auxin target gene induction. These results show that cotyledon-generated auxin regulates hypocotyl elongation. In addition, we find in mature B. rapa plants that simulated shade does not affect seed oil composition but may affect seed yield. This suggests that in field settings where mutual shading between plants may occur, a balance between plant density and seed yield per plant needs to be achieved for maximum oil yield, while oil composition might remain constant.

Ethylene gas is essential for many developmental processes and stress responses in plants. EIN2 plays a key role in ethylene signalling but its function remains enigmatic. Here, we show that ethylene specifically elevates acetylation of histone H3K14 and the non-canonical acetylation of H3K23 in etiolated seedlings. The up-regulation of these two histone marks positively correlates with ethylene-regulated transcription activation, and the elevation requires EIN2. Both EIN2 and EIN3 interact with a SANT domain protein named EIN2 nuclear associated protein 1 (ENAP1), overexpression of which results in elevation of histone acetylation and enhanced ethylene-inducible gene expression in an EIN2-dependent manner. On the basis of these findings we propose a model where, in the presence of ethylene, the EIN2 C terminus contributes to downstream signalling via the elevation of acetylation at H3K14 and H3K23. ENAP1 may potentially mediate ethylene-induced histone acetylation via its interactions with EIN2 C terminus.

Phytochromes are red/far-red light receptors that function in photomorphogenesis of plants. Photoisomerization of phytochrome by red light leads to its translocation to the nucleus, where it regulates gene expression. We examined whether phytochrome is phosphorylated in response to light, and we report that phytochrome B (phyB)'s N terminus contains a region with a number of phosphoserines, threonines, and tyrosines. The light-dependent phosphorylation of tyrosine 104 (Y104) appears to play a negative role in phyB's activity, because a phosphomimic mutant, phyBY104E, is unable to complement any phyB-related phenotype, is defective in binding to its signaling partner PIF3, and fails to form stable nuclear bodies even though it retains normal photochemistry in vitro. In contrast, plants stably expressing a nonphosphorylatable mutant, phyBY104F, are hypersensitive to light. The proper response to changes in the light environment is crucial for plant survival, and our study brings tyrosine phosphorylation to the forefront of light-signaling mechanisms.

Numerous methods are available to profile several epigenetic marks, providing data with different genome coverage and resolution. Large epigenomic datasets are then generated, and often combined with other high-throughput data, including RNA-seq, ChIP-seq for transcription factors (TFs) binding and DNase-seq experiments. Despite the numerous computational tools covering specific steps in the analysis of large-scale epigenomics data, comprehensive software solutions for their integrative analysis are still missing. Multiple tools must be identified and combined to jointly analyze histone marks, TFs binding and other -omics data together with DNA methylation data, complicating the analysis of these data and their integration with publicly available datasets.

DNA microarrays are a well-established technology for measuring gene expression levels. Microarrays designed for this purpose use relatively few probes for each gene and are biased toward known and predicted gene structures. Recently, high-density oligonucleotide-based whole-genome microarrays have emerged as a preferred platform for genomic analysis beyond simple gene expression profiling. Potential uses for such whole-genome arrays include empirical annotation of the transcriptome, chromatin-immunoprecipitation-chip studies, analysis of alternative splicing, characterization of the methylome (the methylation state of the genome), polymorphism discovery and genotyping, comparative genome hybridization, and genome resequencing. Here we review different whole-genome microarray designs and applications of this technology to obtain a wide variety of genomic scale information.

Abscisic acid (ABA) signal transduction has been proposed to utilize cytosolic Ca(2+) in guard cell ion channel regulation. However, genetic mutants in Ca(2+) sensors that impair guard cell or plant ion channel signaling responses have not been identified, and whether Ca(2+)-independent ABA signaling mechanisms suffice for a full response remains unclear. Calcium-dependent protein kinases (CDPKs) have been proposed to contribute to central signal transduction responses in plants. However, no Arabidopsis CDPK gene disruption mutant phenotype has been reported to date, likely due to overlapping redundancies in CDPKs. Two Arabidopsis guard cell-expressed CDPK genes, CPK3 and CPK6, showed gene disruption phenotypes. ABA and Ca(2+) activation of slow-type anion channels and, interestingly, ABA activation of plasma membrane Ca(2+)-permeable channels were impaired in independent alleles of single and double cpk3cpk6 mutant guard cells. Furthermore, ABA- and Ca(2+)-induced stomatal closing were partially impaired in these cpk3cpk6 mutant alleles. However, rapid-type anion channel current activity was not affected, consistent with the partial stomatal closing response in double mutants via a proposed branched signaling network. Imposed Ca(2+) oscillation experiments revealed that Ca(2+)-reactive stomatal closure was reduced in CDPK double mutant plants. However, long-lasting Ca(2+)-programmed stomatal closure was not impaired, providing genetic evidence for a functional separation of these two modes of Ca(2+)-induced stomatal closing. Our findings show important functions of the CPK6 and CPK3 CDPKs in guard cell ion channel regulation and provide genetic evidence for calcium sensors that transduce stomatal ABA signaling.

Direct reprogramming of fibroblasts to neurons induces widespread cellular and transcriptional reconfiguration. Here, we characterized global epigenomic changes during the direct reprogramming of mouse fibroblasts to neurons using whole-genome base-resolution DNA methylation (mC) sequencing. We found that the pioneer transcription factor Ascl1 alone is sufficient for inducing the uniquely neuronal feature of non-CG methylation (mCH), but co-expression of Brn2 and Mytl1 was required to establish a global mCH pattern reminiscent of mature cortical neurons. Ascl1 alone induced promoter CG methylation (mCG) of fibroblast specific genes, while BAM overexpression additionally targets a competing myogenic program and directs a more faithful conversion to neuronal cells. Ascl1 induces local demethylation at its binding sites. Surprisingly, co-expression with Brn2 and Mytl1 inhibited the ability of Ascl1 to induce demethylation, suggesting a contextual regulation of transcription factor - epigenome interaction. Finally, we found that de novo methylation by DNMT3A is required for efficient neuronal reprogramming.

Unlike animals, plants can pause their life cycle as dormant seeds. In both plants and animals, DNA methylation is involved in the regulation of gene expression and genome integrity. In animals, reprogramming erases and re-establishes DNA methylation during development. However, knowledge of reprogramming or reconfiguration in plants has been limited to pollen and the central cell. To better understand epigenetic reconfiguration in the embryo, which forms the plant body, we compared time-series methylomes of dry and germinating seeds to publicly available seed development methylomes.

A detailed understanding of abiotic stress tolerance in plants is essential to provide food security in the face of increasingly harsh climatic conditions. Glucosinolates (GLSs) are secondary metabolites found in the Brassicaceae that protect plants from herbivory and pathogen attack. Here we report that in Arabidopsis, aliphatic GLS levels are regulated by the auxin-sensitive Aux/IAA repressors IAA5, IAA6, and IAA19. These proteins act in a transcriptional cascade that maintains expression of GLS levels when plants are exposed to drought conditions. Loss of IAA5/6/19 results in reduced GLS levels and decreased drought tolerance. Further, we show that this phenotype is associated with a defect in stomatal regulation. Application of GLS to the iaa5,6,19 mutants restores stomatal regulation and normal drought tolerance. GLS action is dependent on the receptor kinase GHR1, suggesting that GLS may signal via reactive oxygen species. These results provide a novel connection between auxin signaling, GLS levels and drought response.

Cytokinin fulfills its diverse roles in planta through a series of transcriptional responses. We identify the in vivo DNA binding site profiles for three genetically redundant type-B ARABIDOPSIS RESPONSE REGULATORS (B-ARRs): ARR1, ARR10, and ARR12. The expression and genome-wide DNA binding locations of the three B-ARRs extensively overlap. Constructing a primary cytokinin response transcriptional network reveals a recurring theme of widespread cross-regulation between the components of the cytokinin pathway and other plant hormone pathways. The B-ARRs are found to have similar DNA binding motifs, though sequences flanking the core motif were degenerate. Cytokinin treatments amalgamate the three different B-ARRs motifs to identical DNA binding signatures (AGATHY, H(a/t/c), Y(t/c)) which suggests cytokinin may regulate binding activity of B-ARR family members. Furthermore, we find that WUSCHEL, a key gene required for apical meristem maintenance, is a cytokinin-dependent B-ARR target gene, demonstrating the importance of the cytokinin transcription factor network in shoot development.

A comprehensive characterization of neuronal cell types, their distributions, and patterns of connectivity is critical for understanding the properties of neural circuits and how they generate behaviors. Here we review the experiences of the BRAIN Initiative Cell Census Consortium, ten pilot projects funded by the U.S. BRAIN Initiative, in developing, validating, and scaling up emerging genomic and anatomical mapping technologies for creating a complete inventory of neuronal cell types and their connections in multiple species and during development. These projects lay the foundation for a larger and longer-term effort to generate whole-brain cell atlases in species including mice and humans.

The Encyclopedia of DNA Elements (ENCODE) project has established a genomic resource for mammalian development, profiling a diverse panel of mouse tissues at 8 developmental stages from 10.5 days after conception until birth, including transcriptomes, methylomes and chromatin states. Here we systematically examined the state and accessibility of chromatin in the developing mouse fetus. In total we performed 1,128 chromatin immunoprecipitation with sequencing (ChIP-seq) assays for histone modifications and 132 assay for transposase-accessible chromatin using sequencing (ATAC-seq) assays for chromatin accessibility across 72 distinct tissue-stages. We used integrative analysis to develop a unified set of chromatin state annotations, infer the identities of dynamic enhancers and key transcriptional regulators, and characterize the relationship between chromatin state and accessibility during developmental gene regulation. We also leveraged these data to link enhancers to putative target genes and demonstrate tissue-specific enrichments of sequence variants associated with disease in humans. The mouse ENCODE data sets provide a compendium of resources for biomedical researchers and achieve, to our knowledge, the most comprehensive view of chromatin dynamics during mammalian fetal development to date.

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

Strigolactones (SLs) are a diverse class of butenolide-bearing plant hormones associated with several processes of major agricultural concern. SLs initiate symbiosis between plants and arbuscular mycorrhizal fungi, cause germination of crop-devastating parasitic plants, and inhibit shoot branching in vascular plants. SLs are perceived by dual receptor-hydrolase proteins, and capturing the intact ligand inside the receptor remains a key challenge for structural biologists. In addition, many discovered SLs are hard to obtain and too unstable to work with. In a computer-based approach, we investigated the interaction of 20 different SL molecules with nine crystal structures of SL receptors. Our results suggest an important role of the active site for ligand binding and orientation, and that the parasitic plant Striga hermonthica has developed both promiscuous and type-specific SL receptors as part of its host recognition strategy.

Integrating large single-cell gene expression, chromatin accessibility and DNA methylation datasets requires general and scalable computational approaches. Here we describe online integrative non-negative matrix factorization (iNMF), an algorithm for integrating large, diverse and continually arriving single-cell datasets. Our approach scales to arbitrarily large numbers of cells using fixed memory, iteratively incorporates new datasets as they are generated and allows many users to simultaneously analyze a single copy of a large dataset by streaming it over the internet. Iterative data addition can also be used to map new data to a reference dataset. Comparisons with previous methods indicate that the improvements in efficiency do not sacrifice dataset alignment and cluster preservation performance. We demonstrate the effectiveness of online iNMF by integrating more than 1 million cells on a standard laptop, integrating large single-cell RNA sequencing and spatial transcriptomic datasets, and iteratively constructing a single-cell multi-omic atlas of the mouse motor cortex.