Joshua G. Steffen

Multiple reference genomes and transcriptomes for Arabidopsis thaliana

Genetic differences between Arabidopsis thaliana accessions underlie the plant’s extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.

The Capsella rubella genome and the genomic consequences of rapid mating system evolution

The shift from outcrossing to selfing is common in flowering plants, but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible <200,000 years ago. We report a C. rubella reference genome sequence and compare RNA expression and polymorphism patterns between C. rubella and its outcrossing progenitor Capsella grandiflora. We found a clear shift in the expression of genes associated with flowering phenotypes, similar to that seen in Arabidopsis, in which self fertilization evolved about 1 million years ago. Comparisons of the two Capsella species showed evidence of rapid genome-wide relaxation of purifying selection in C. rubella without a concomitant change in transposable element abundance. Overall we document that the transition to selfing may be typified by parallel shifts in gene expression, along with a measurable reduction of purifying selection.

An atlas of over 90,000 conserved noncoding sequences provides insight into crucifer regulatory regions

Despite the central importance of noncoding DNA to gene regulation and evolution, understanding of the extent of selection on plant noncoding DNA remains limited compared to that of other organisms. Here we report sequencing of genomes from three Brassicaceae species (Leavenworthia alabamica, Sisymbrium irio and Aethionema arabicum) and their joint analysis with six previously sequenced crucifer genomes. Conservation across orthologous bases suggests that at least 17% of the Arabidopsis thaliana genome is under selection, with nearly one-quarter of the sequence under selection lying outside of coding regions. Much of this sequence can be localized to approximately 90,000 conserved noncoding sequences (CNSs) that show evidence of transcriptional and post-transcriptional regulation. Population genomics analyses of two crucifer species, A. thaliana and Capsella grandiflora, confirm that most of the identified CNSs are evolving under medium to strong purifying selection. Overall, these CNSs highlight both similarities and several key differences between the regulatory DNA of plants and other species.

AGL80 Is Required for Central Cell and Endosperm Development in Arabidopsis

During plant reproduction, the central cell of the female gametophyte becomes fertilized to produce the endosperm, a storage tissue that nourishes the developing embryo within the seed. The molecular mechanisms controlling the specification and differentiation of the central cell are poorly understood. We identified a female gametophyte mutant in Arabidopsis thaliana, fem111, that is affected in central cell development. In fem111 female gametophytes, the central cells nucleolus and vacuole fail to mature properly. In addition, endosperm development is not initiated after fertilization of fem111 female gametophytes. fem111 contains a T-DNA insertion in AGAMOUS-LIKE80 (AGL80). FEM111/AGL80 is a member of the MADS box family of genes that likely encode transcription factors. An AGL80–green fluorescent protein fusion protein is localized to the nucleus. Within the ovule and seed, FEM111/AGL80 is expressed exclusively in the central cell and uncellularized endosperm. FEM111/AGL80 expression is also detected in roots, leaves, floral stems, anthers, and young flowers by real-time RT-PCR. FEM111/AGL80 is required for the expression of two central cell–expressed genes, DEMETER and DD46, but not for a third central cell–expressed gene, FERTILIZATION-INDEPENDENT SEED2. Together, these data suggest that FEM111/AGL80 functions as a transcription factor within the central cell gene regulatory network and controls the expression of downstream genes required for central cell development and function.

The AGL62 MADS Domain Protein Regulates Cellularization during Endosperm Development in Arabidopsis

Endosperm, a storage tissue in the angiosperm seed, provides nutrients to the embryo during seed development and/or to the developing seedling during germination. A major event in endosperm development is the transition between the syncytial phase, during which the endosperm nuclei undergo many rounds of mitosis without cytokinesis, and the cellularized phase, during which cell walls form around the endosperm nuclei. The molecular processes controlling this phase transition are not understood. In agl62 seeds, the endosperm cellularizes prematurely, indicating that AGL62 is required for suppression of cellularization during the syncytial phase. AGL62 encodes a Type I MADS domain protein that likely functions as a transcription factor. During seed development, AGL62 is expressed exclusively in the endosperm. During wild-type endosperm development, AGL62 expression is strong during the syncytial phase and then declines abruptly just before cellularization. By contrast, in mutant seeds containing defects in some FERTILIZATION-INDEPENDENT SEED (FIS) class Polycomb group genes, the endosperm fails to cellularize and AGL62 expression fails to decline. Together, these data suggest that AGL62 suppresses cellularization during the syncytial phase of endosperm development and that endosperm cellularization is triggered via direct or indirect AGL62 inactivation by the FIS polycomb complex.

AGL61 Interacts with AGL80 and Is Required for Central Cell Development in Arabidopsis

The central cell of the female gametophyte plays a role in pollen tube guidance and in regulating the initiation of endosperm development. Following fertilization, the central cell gives rise to the seeds endosperm, which nourishes the developing embryo within the seed. The molecular mechanisms controlling specification and differentiation of the central cell are poorly understood. We identified AGL61 in a screen for transcription factor genes expressed in the female gametophyte. AGL61 encodes a Type I MADS domain protein, which likely functions as a transcription factor. Consistent with this, an AGL61-green fluorescent protein fusion protein is localized to the nucleus. In the context of the ovule and seed, AGL61 is expressed exclusively in the central cell and early endosperm. agl61 female gametophytes are affected in the central cell specifically. The morphological defects include an overall reduction in size of the central cell and a reduced or absent central cell vacuole. When fertilized with wild-type pollen, agl61 central cells fail to give rise to endosperm. In addition, synergid- and antipodal-expressed genes are ectopically expressed in agl61 central cells. The expression pattern and mutant phenotype of AGL61 are similar to those of AGL80, suggesting that AGL61 may function as a heterodimer with AGL80 within the central cell; consistent with this, AGL61 and AGL80 interact in yeast two-hybrid assays. Together, these data suggest that AGL61 functions as a transcription factor and controls the expression of downstream genes during central cell development.

Identification of genes expressed in the angiosperm female gametophyte

Until recently, identification of gene regulatory networks controlling the development of the angiosperm female gametophyte has presented a significant challenge to the plant biology community. The angiosperm female gametophyte is fairly inaccessible because it is a highly reduced structure relative to the sporophyte and is embedded within multiple layers of the sporophytic tissue of the ovule. Moreover, although mutations affecting the female gametophyte can be readily isolated, their analysis can be difficult because most affect genes involved in basic cellular processes that are also required in the diploid sporophyte. In recent years, expression-based approaches in multiple species have begun to uncover gene sets expressed in specific female gametophyte cells as a means of identifying regulatory networks controlling cell differentiation in the female gametophyte. Here, recent efforts to identify and analyse gene expression programmes in the Arabidopsis female gametophyte are reviewed.

RNA-Seq analysis of laser-capture microdissected cells of the developing central starchy endosperm of maize.

Endosperm is a product of double fertilization, and provides nutrients and signals to the embryo during seed development in flowering plants. Early stages of endosperm development are critical for the development of its storage capacity through synthesis and accumulation of starch and storage proteins. Here we report on the isolation and sequencing of mRNAs from the central portion of the starchy endosperm of Zea mays (maize) B73 at 6 days after pollination. We detected a high level of correlation among the four biological replicates of RNAs isolated using laser-capture microdissection of the cell type. Because the assayed developmental stage precedes the synthesis and accumulation of the major storage proteins and starch in the endosperm, our dataset likely include mRNAs for genes that are involved in control and establishment of these storage programs. The mRNA-Seq data has been deposited in Gene Expression Omnibus (accession number GSE58504).

SUPPRESSOR OF FRIGIDA (SUF4) supports gamete fusion via regulating arabidopsis EC1 gene expression

In Arabidopsis, gamete fusion requires the C2H2 transcription factor SUF4, which regulates the expression of the EGG CELL1 gene family. The EGG CELL1 (EC1) gene family of Arabidopsis (Arabidopsis thaliana) comprises five members that are specifically expressed in the egg cell and redundantly control gamete fusion during double fertilization. We investigated the activity of all five EC1 promoters in promoter-deletion studies and identified SUF4 (SUPPRESSOR OF FRIGIDA4), a C2H2 transcription factor, as a direct regulator of the EC1 gene expression. In particular, we demonstrated that SUF4 binds to all five Arabidopsis EC1 promoters, thus regulating their expression. The down-regulation of SUF4 in homozygous suf4-1 ovules results in reduced EC1 expression and delayed sperm fusion, which can be rescued by expressing SUF4-β-glucuronidase under the control of the SUF4 promoter. To identify more gene products able to regulate EC1 expression together with SUF4, we performed coexpression studies that led to the identification of MOM1 (MORPHEUS’ MOLECULE1), a component of a silencing mechanism that is independent of DNA methylation marks. In mom1-3 ovules, both SUF4 and EC1 genes are down-regulated, and EC1 genes show higher levels of histone 3 lysine-9 acetylation, suggesting that MOM1 contributes to the regulation of SUF4 and EC1 gene expression.

Genomic Rearrangements in Arabidopsis Considered as Quantitative Traits

Structural Rearrangements can have unexpected effects on quantitative phenotypes. Surprisingly, these rearrangements can also be considered as... To understand the population genetics of structural variants and their effects on phenotypes, we developed an approach to mapping structural variants that segregate in a population sequenced at low coverage. We avoid calling structural variants directly. Instead, the evidence for a potential structural variant at a locus is indicated by variation in the counts of short-reads that map anomalously to that locus. These structural variant traits are treated as quantitative traits and mapped genetically, analogously to a gene expression study. Association between a structural variant trait at one locus, and genotypes at a distant locus indicate the origin and target of a transposition. Using ultra-low-coverage (0.3×) population sequence data from 488 recombinant inbred Arabidopsis thaliana genomes, we identified 6502 segregating structural variants. Remarkably, 25% of these were transpositions. While many structural variants cannot be delineated precisely, we validated 83% of 44 predicted transposition breakpoints by polymerase chain reaction. We show that specific structural variants may be causative for quantitative trait loci for germination and resistance to infection by the fungus Albugo laibachii, isolate Nc14. Further we show that the phenotypic heritability attributable to read-mapping anomalies differs from, and, in the case of time to germination and bolting, exceeds that due to standard genetic variation. Genes within structural variants are also more likely to be silenced or dysregulated. This approach complements the prevalent strategy of structural variant discovery in fewer individuals sequenced at high coverage. It is generally applicable to large populations sequenced at low-coverage, and is particularly suited to mapping transpositions.