Karen M. McGinnis

The Phaeodactylum genome reveals the evolutionary history of diatom genomes.

Diatoms are photosynthetic secondary endosymbionts found throughout marine and freshwater environments, and are believed to be responsible for around one-fifth of the primary productivity on Earth. The genome sequence of the marine centric diatom Thalassiosira pseudonana was recently reported, revealing a wealth of information about diatom biology. Here we report the complete genome sequence of the pennate diatom Phaeodactylum tricornutum and compare it with that of T. pseudonana to clarify evolutionary origins, functional significance and ubiquity of these features throughout diatoms. In spite of the fact that the pennate and centric lineages have only been diverging for 90 million years, their genome structures are dramatically different and a substantial fraction of genes (∼40%) are not shared by these representatives of the two lineages. Analysis of molecular divergence compared with yeasts and metazoans reveals rapid rates of gene diversification in diatoms. Contributing factors include selective gene family expansions, differential losses and gains of genes and introns, and differential mobilization of transposable elements. Most significantly, we document the presence of hundreds of genes from bacteria. More than 300 of these gene transfers are found in both diatoms, attesting to their ancient origins, and many are likely to provide novel possibilities for metabolite management and for perception of environmental signals. These findings go a long way towards explaining the incredible diversity and success of the diatoms in contemporary oceans.

The Arabidopsis SLEEPY1 Gene Encodes a Putative F-Box Subunit of an SCF E3 Ubiquitin Ligase

The Arabidopsis SLY1 (SLEEPY1) gene positively regulates gibberellin (GA) signaling. Positional cloning of SLY1 revealed that it encodes a putative F-box protein. This result suggests that SLY1 is the F-box subunit of an SCF E3 ubiquitin ligase that regulates GA responses. The DELLA domain protein RGA (repressor of ga1-3) is a repressor of GA response that appears to undergo GA-stimulated protein degradation. RGA is a potential substrate of SLY1, because sly1 mutations cause a significant increase in RGA protein accumulation even after GA treatment. This result suggests SCFSLY1-targeted degradation of RGA through the 26S proteasome pathway. Further support for this model is provided by the observation that an rga null allele partially suppresses the sly1-10 mutant phenotype. The predicted SLY1 amino acid sequence is highly conserved among plants, indicating a key role in GA response.

An RNA-dependent RNA polymerase is required for paramutation in maize

Paramutation is an allele-dependent transfer of epigenetic information, which results in the heritable silencing of one allele by another. Paramutation at the b1 locus in maize is mediated by unique tandem repeats that communicate in trans to establish and maintain meiotically heritable transcriptional silencing. The mop1 (mediator of paramutation1) gene is required for paramutation, and mop1 mutations reactivate silenced Mutator elements. Plants carrying mutations in the mop1 gene also stochastically exhibit pleiotropic developmental phenotypes. Here we report the map-based cloning of mop1, an RNA-dependent RNA polymerase gene (RDRP), most similar to the RDRP in plants that is associated with the production of short interfering RNA (siRNA) targeting chromatin. Nuclear run-on assays reveal that the tandem repeats required for b1 paramutation are transcribed from both strands, but siRNAs were not detected. We propose that the mop1 RDRP is required to maintain a threshold level of repeat RNA, which functions in trans to establish and maintain the heritable chromatin states associated with paramutation.

Genome-wide analyses of alternative splicing in plants: Opportunities and challenges

Alternative splicing (AS) creates multiple mRNA transcripts from a single gene. While AS is known to contribute to gene regulation and proteome diversity in animals, the study of its importance in plants is in its early stages. However, recently available plant genome and transcript sequence data sets are enabling a global analysis of AS in many plant species. Results of genome analysis have revealed differences between animals and plants in the frequency of alternative splicing. The proportion of plant genes that have one or more alternative transcript isoforms is approximately 20%, indicating that AS in plants is not rare, although this rate is approximately one-third of that observed in human. The majority of plant AS events have not been functionally characterized, but evidence suggests that AS participates in important plant functions, including stress response, and may impact domestication and trait selection. The increasing availability of plant genome sequence data will enable larger comparative analyses that will identify functionally important plant AS events based on their evolutionary conservation, determine the influence of genome duplication on the evolution of AS, and discover plant-specific cis-elements that regulate AS. This review summarizes recent analyses of AS in plants, discusses the importance of further analysis, and suggests directions for future efforts.

Computational identification and functional predictions of long noncoding RNA in Zea mays.

Background Computational analysis of cDNA sequences from multiple organisms suggests that a large portion of transcribed DNA does not code for a functional protein. In mammals, noncoding transcription is abundant, and often results in functional RNA molecules that do not appear to encode proteins. Many long noncoding RNAs (lncRNAs) appear to have epigenetic regulatory function in humans, including HOTAIR and XIST. While epigenetic gene regulation is clearly an essential mechanism in plants, relatively little is known about the presence or function of lncRNAs in plants. Methodology/Principal Findings To explore the connection between lncRNA and epigenetic regulation of gene expression in plants, a computational pipeline using the programming language Python has been developed and applied to maize full length cDNA sequences to identify, classify, and localize potential lncRNAs. The pipeline was used in parallel with an SVM tool for identifying ncRNAs to identify the maximal number of ncRNAs in the dataset. Although the available library of sequences was small and potentially biased toward protein coding transcripts, 15% of the sequences were predicted to be noncoding. Approximately 60% of these sequences appear to act as precursors for small RNA molecules and may function to regulate gene expression via a small RNA dependent mechanism. ncRNAs were predicted to originate from both genic and intergenic loci. Of the lncRNAs that originated from genic loci, ∼20% were antisense to the host gene loci. Conclusions/Significance Consistent with similar studies in other organisms, noncoding transcription appears to be widespread in the maize genome. Computational predictions indicate that maize lncRNAs may function to regulate expression of other genes through multiple RNA mediated mechanisms.

Spreading of Heterochromatin Is Limited to Specific Families of Maize Retrotransposons

Transposable elements (TEs) have the potential to act as controlling elements to influence the expression of genes and are often subject to heterochromatic silencing. The current paradigm suggests that heterochromatic silencing can spread beyond the borders of TEs and influence the chromatin state of neighboring low-copy sequences. This would allow TEs to condition obligatory or facilitated epialleles and act as controlling elements. The maize genome contains numerous families of class I TEs (retrotransposons) that are present in moderate to high copy numbers, and many are found in regions near genes, which provides an opportunity to test whether the spreading of heterochromatin from retrotransposons is prevalent. We have investigated the extent of heterochromatin spreading into DNA flanking each family of retrotransposons by profiling DNA methylation and di-methylation of lysine 9 of histone 3 (H3K9me2) in low-copy regions of the maize genome. The effects of different retrotransposon families on local chromatin are highly variable. Some retrotransposon families exhibit enrichment of heterochromatic marks within 800–1,200 base pairs of insertion sites, while other families exhibit very little evidence for the spreading of heterochromatic marks. The analysis of chromatin state in genotypes that lack specific insertions suggests that the heterochromatin in low-copy DNA flanking retrotransposons often results from the spreading of silencing marks rather than insertion-site preferences. Genes located near TEs that exhibit spreading of heterochromatin tend to be expressed at lower levels than other genes. Our findings suggest that a subset of retrotransposon families may act as controlling elements influencing neighboring sequences, while the majority of retrotransposons have little effect on flanking sequences.

Transgene-induced RNA interference as a tool for plant functional genomics

RNA interference (RNAi) is a powerful tool for functional genomics in a number of species. The logistics and procedures for doing high-throughput RNAi to investigate the functions of large numbers of genes in Arabidopsis thaliana and in Zea mays are described. Publicly available plasmid vectors that facilitate the stable chromosomal integration of inverted repeat transgenes that trigger RNAi have been used to generate more than 50 independent transgenic lines each in Arabidopsis and maize. Analysis of mRNA abundance of the targeted genes in independent lines transformed with distinct constructs indicates that the success of RNAi-induced silencing is gene dependent. mRNA levels were not detectably reduced for some genes, but were dramatically reduced for a number of genes targeted. A common pattern was that multiple independent lines transgenic for the same construct showed the same extent of silencing. This chapter describes the procedures used to generate and test transgenic lines mediating RNAi in Arabidopsis and maize.

Genetic Perturbation of the Maize Methylome

Genetic analyses of maize genes in DNA methylation pathways reveal differences between maize and Arabidopsis, including evidence that DNA methylation is required for growth and development in maize. DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

Assessing the Efficiency of RNA Interference for Maize Functional Genomics

A large-scale functional genomics project was initiated to study the function of chromatin-related genes in maize (Zea mays). Transgenic lines containing short gene segments in inverted repeat orientation designed to reduce expression of target genes by RNA interference (RNAi) were isolated, propagated, and analyzed in a variety of assays. Analysis of the selectable marker expression over multiple generations revealed that most transgenes were transmitted faithfully, whereas some displayed reduced transmission or transgene silencing. A range of target-gene silencing efficiencies, from nondetectable silencing to nearly complete silencing, was revealed by semiquantitative reverse transcription-PCR analysis of transcript abundance for the target gene. In some cases, the RNAi construct was able to cause a reduction in the steady-state RNA levels of not only the target gene, but also another closely related gene. Correlation of silencing efficiency with expression level of the target gene and sequence features of the inverted repeat did not reveal any factors capable of predicting the silencing success of a particular RNAi-inducing construct. The frequencies of success of this large-scale project in maize, together with parameters for optimization at various steps, should serve as a useful framework for designing future RNAi-based functional genomics projects in crop plants.

Accessible DNA and Relative Depletion of H3K9me2 at Maize Loci Undergoing RNA-Directed DNA Methylation

Only a small fraction of the maize genome undergoes RNA-directed DNA methylation and, for several characteristics, the chromatin in these areas resembles euchromatin more than heterochromatin. RNA-directed DNA methylation (RdDM) in plants is a well-characterized example of RNA interference-related transcriptional gene silencing. To determine the relationships between RdDM and heterochromatin in the repeat-rich maize (Zea mays) genome, we performed whole-genome analyses of several heterochromatic features: dimethylation of lysine 9 and lysine 27 (H3K9me2 and H3K27me2), chromatin accessibility, DNA methylation, and small RNAs; we also analyzed two mutants that affect these processes, mediator of paramutation1 and zea methyltransferase2. The data revealed that the majority of the genome exists in a heterochromatic state defined by inaccessible chromatin that is marked by H3K9me2 and H3K27me2 but that lacks RdDM. The minority of the genome marked by RdDM was predominantly near genes, and its overall chromatin structure appeared more similar to euchromatin than to heterochromatin. These and other data indicate that the densely staining chromatin defined as heterochromatin differs fundamentally from RdDM-targeted chromatin. We propose that small interfering RNAs perform a specialized role in repressing transposons in accessible chromatin environments and that the bulk of heterochromatin is incompatible with small RNA production.