Differentiation of the maize subgenomes by genome dominance and both ancient and ongoing gene loss

Characterization of Rubisco activase genes in maize: an α-isoform gene functions alongside a β-isoform gene

Rubisco activase (RCA) catalyzes the activation of Rubisco in vivo and plays a crucial role in regulating plant growth. In maize (Zea mays), only β-form RCA genes have been cloned and characterized. In this study, a genome-wide survey revealed the presence of an α-form RCA gene and a β-form RCA gene in the maize genome, herein referred to as ZmRCAα and ZmRCAβ, respectively. An analysis of genomic DNA and complementary DNA sequences suggested that alternative splicing of the ZmRCAβ precursor mRNA (premRNA) at its 3' untranslated region could produce two distinctive ZmRCAβ transcripts. Analyses by electrophoresis and matrix-assisted laser desorption/ionization-tandem time-of-flight mass spectrometry showed that ZmRCAα and ZmRCAβ encode larger and smaller polypeptides of approximately 46 and 43 kD, respectively. Transcriptional analyses demonstrated that the expression levels of both ZmRCAα and ZmRCAβ were higher in leaves and during grain filling and that expression followed a specific cyclic day/night pattern. In 123 maize inbred lines with extensive genetic diversity, the transcript abundance and protein expression levels of these two RCA genes were positively correlated with grain yield. Additionally, both genes demonstrated a similar correlation with grain yield compared with three C₄ photosynthesis genes. Our data suggest that, in addition to the β-form RCA-encoding gene, the α-form RCA-encoding gene also contributes to the synthesis of RCA in maize and support the hypothesis that RCA genes may play an important role in determining maize productivity.

A near complete snapshot of the Zea mays seedling transcriptome revealed from ultra-deep sequencing

RNA-sequencing (RNA-seq) enables in-depth exploration of transcriptomes, but typical sequencing depth often limits its comprehensiveness. In this study, we generated nearly 3 billion RNA-Seq reads, totaling 341 Gb of sequence, from a Zea mays seedling sample. At this depth, a near complete snapshot of the transcriptome was observed consisting of over 90% of the annotated transcripts, including lowly expressed transcription factors. A novel hybrid strategy combining de novo and reference-based assemblies yielded a transcriptome consisting of 126,708 transcripts with 88% of expressed known genes assembled to full-length. We improved current annotations by adding 4,842 previously unannotated transcript variants and many new features, including 212 maize transcripts, 201 genes, 10 genes with undocumented potential roles in seedlings as well as maize lineage specific gene fusion events. We demonstrated the power of deep sequencing for large transcriptome studies by generating a high quality transcriptome, which provides a rich resource for the research community.

Origin, inheritance, and gene regulatory consequences of genome dominance in polyploids

Whole-genome duplications happen repeatedly in a typical flowering plant lineage. Following most ancient tetraploidies, the two subgenomes are distinguishable because one subgenome, the dominant subgenome, tends to have more genes than the other subgenome. Additionally, among retained pairs, the gene on the dominant subgenome tends to be expressed more than its recessive homeolog. Using comparative genomics, we show that genome dominance is heritable. The dominant subgenome of one postpolyploidy event remains dominant through a subsequent polyploidy event. We show that transposon-derived 24-nt RNAs target and cover the upstream region of retained genes preferentially when located on the recessive subgenome, and with little regard for a gene's level of expression. We hypothesize that small RNA (smRNA)-mediated silencing of transposons near genes causes position-effect down-regulation. Unlike 24-nt smRNA coverage, transposon coverage tracks gene expression, so not all transposons behave identically. We propose that successful ancient tetraploids begin as wide crosses between two lines, each evolved for different tradeoffs between transposon silencing and negative position effects on gene expression. We hypothesize that following a chaotic wide-cross/new tetraploid period, genes acquire their new expression balances based on differences in transposon coverage in the parents. We envision patches of silenceable transposon as quantitative cis-regulators of baseline transcription rate. Attractive solutions to heterosis and the C-value paradox are mentioned.

An interspecific fungal hybrid reveals cross-kingdom rules for allopolyploid gene expression patterns

Polyploidy, a state in which the chromosome complement has undergone an increase, is a major force in evolution. Understanding the consequences of polyploidy has received much attention, and allopolyploids, which result from the union of two different parental genomes, are of particular interest because they must overcome a suite of biological responses to this merger, known as “genome shock.” A key question is what happens to gene expression of the two gene copies following allopolyploidization, but until recently the tools to answer this question on a genome-wide basis were lacking. Here we utilize high throughput transcriptome sequencing to produce the first genome-wide picture of gene expression response to allopolyploidy in fungi. A novel pipeline for assigning sequence reads to the gene copies was used to quantify their expression in a fungal allopolyploid. We find that the transcriptional response to allopolyploidy is predominantly conservative: both copies of most genes are retained; over half the genes inherit parental gene expression patterns; and parental differential expression is often lost in the allopolyploid. Strikingly, the patterns of gene expression change are highly concordant with the genome-wide expression results of a cotton allopolyploid. The very different nature of these two allopolyploids implies a conserved, eukaryote-wide transcriptional response to genome merger. We provide evidence that the transcriptional responses we observe are mostly driven by intrinsic differences between the regulatory systems in the parent species, and from this propose a mechanistic model in which the cross-kingdom conservation in transcriptional response reflects conservation of the mutational processes underlying eukaryotic gene regulatory evolution. This work provides a platform to develop a universal understanding of gene expression response to allopolyploidy and suggests that allopolyploids are an exceptional system to investigate gene regulatory changes that have evolved in the parental species prior to allopolyploidization.

Organisms are complex biological systems that must continue to function even as their genomes evolve. While evolution is usually gradual, the formation of new species by the hybridization of different parents—allopolyploidization—occurs nearly instantaneously. A key question is what happens to expression of the two parental gene copies following genome merger. To determine this, we focused on a fungal allopolyploid from a group that dominates many of the world's pastoral economies. To investigate the fate of gene expression in this system, we developed a novel pipeline to assign high throughput RNA sequence reads to the two parental gene copies, thus allowing quantification of expression. We found transcriptional responses to be predominantly conservative: most gene copies either inherit parental expression patterns, or if differentially expressed in the parents, that difference is lost in the hybrid. Moreover, we identified an extraordinary level of concordance in the fate of genome-wide allopolyploid gene expression with that seen in cotton. The very different nature of these two allopolyploids suggests that there is a set of universal rules for the transcriptional response to genome merger. We propose a mechanistic model whereby this conserved response reflects similarities in mutational processes that underlie gene regulatory evolution.

Polyploids as a "model system" for the study of heterosis

Heterosis research over the past century has focused primarily on diploid plants and animals. This is despite the fact that most heterotic organisms contain polyploid events in their recent and/or ancient past and various important crop species are heterotic polyploids. We present an argument for the study of heterosis within polyploid systems and give examples of how their study can improve current hypotheses and generate new ones. Polyploid systems allow experiments not possible in diploids but the insights gained must be incorporated into models to explain heterosis at all levels.

Genome-wide identification, splicing, and expression analysis of the myosin gene family in maize (Zea mays)

The actin-based myosin system is essential for the organization and dynamics of the endomembrane system and transport network in plant cells. Plants harbour two unique myosin groups, class VIII and class XI, and the latter is structurally and functionally analogous to the animal and fungal class V myosin. Little is known about myosins in grass, even though grass includes several agronomically important cereal crops. Here, we identified 14 myosin genes from the genome of maize (Zea mays). The relatively larger sizes of maize myosin genes are due to their much longer introns, which are abundant in transposable elements. Phylogenetic analysis indicated that maize myosin genes could be classified into class VIII and class XI, with three and 11 members, respectively. Apart from subgroup XI-F, the remaining subgroups were duplicated at least in one analysed lineage, and the duplication events occurred more extensively in Arabidopsis than in maize. Only two pairs of maize myosins were generated from segmental duplication. Expression analysis revealed that most maize myosin genes were expressed universally, whereas a few members (XI-1, -6, and -11) showed an anther-specific pattern, and many underwent extensive alternative splicing. We also found a short transcript at the O1 locus, which conceptually encoded a headless myosin that most likely functions at the transcriptional level rather than via a dominant-negative mechanism at the translational level. Together, these data provide significant insights into the evolutionary and functional characterization of maize myosin genes that could transfer to the identification and application of homologous myosins of other grasses.

Practical halving; the Nelumbo nucifera evidence on early eudicot evolution

We present a stepwise optimal genome halving algorithm designed for large eukaryote genomes with largely single-copy genes, taking advantage of a signature pattern of paralog distribution in ancient polyploids. This is applied to the genome of Nelumbo nucifera, the sacred lotus, which is the descendant of a duplicated basal eudicot genome. In concert with the reconstructed ancestor of the grape, we investigate early events in eudicot evolution and show that the chromosome number of the common ancestor of lotus and grape was likely between 5 and 7. We show that the duplication of the ancestor of lotus and the triplication of the ancestor of grape were not closely preceded by any additional such event before the divergence of their two lineages.

Extensive translational regulation of gene expression in an allopolyploid (Glycine dolichocarpa)

All flowering plants have experienced repeated rounds of polyploidy (whole-genome duplication), which has in turn driven the evolution of novel phenotypes and ecological tolerances and been a major driver of speciation. The effects of polyploidy on gene expression have been studied extensively at the level of transcription and, to a much lesser extent, at the level of the steady state proteome, but not at the level of translation. We used polysome profiling by RNA-Seq to quantify translational regulation of gene expression in a recently formed (∼100,000 years ago) allotetraploid (Glycine dolichocarpa) closely related to the cultivated soybean (Glycine max). We show that there is a high level of concordance between the allopolyploid transcriptome and translatome overall but that at least one-quarter of the transcriptome is translationally regulated. We further show that translational regulation preferentially targets genes involved in transcription, translation, and photosynthesis, causes regional and possibly whole-chromosome shifts in expression bias between duplicated genes (homoeologs), and reduces transcriptional differences between the polyploid and its diploid progenitors, possibly attenuating misregulation resulting from genome merger and/or doubling. Finally, translational regulation correlates positively with long-term retention of homoeologs from a paleopolyploidy event, suggesting that it plays a significant role in polyploid evolution.

Insights into the maize pan-genome and pan-transcriptome

Genomes at the species level are dynamic, with genes present in every individual (core) and genes in a subset of individuals (dispensable) that collectively constitute the pan-genome. Using transcriptome sequencing of seedling RNA from 503 maize (Zea mays) inbred lines to characterize the maize pan-genome, we identified 8681 representative transcript assemblies (RTAs) with 16.4% expressed in all lines and 82.7% expressed in subsets of the lines. Interestingly, with linkage disequilibrium mapping, 76.7% of the RTAs with at least one single nucleotide polymorphism (SNP) could be mapped to a single genetic position, distributed primarily throughout the nonpericentromeric portion of the genome. Stepwise iterative clustering of RTAs suggests, within the context of the genotypes used in this study, that the maize genome is restricted and further sampling of seedling RNA within this germplasm base will result in minimal discovery. Genome-wide association studies based on SNPs and transcript abundance in the pan-genome revealed loci associated with the timing of the juvenile-to-adult vegetative and vegetative-to-reproductive developmental transitions, two traits important for fitness and adaptation. This study revealed the dynamic nature of the maize pan-genome and demonstrated that a substantial portion of variation may lie outside the single reference genome for a species.

Genome triplication drove the diversification of Brassica plants

The genus Brassica belongs to the plant family Brassicaceae, which includes many important crop species that are used as oilseed, condiments, or vegetables throughout the world. Brassica plants comprise many diverse species, and each species contains rich morphotypes showing extreme traits. Brassica species experienced an extra whole genome triplication (WGT) event compared with the model plant Arabidopsis thaliana. Whole genome sequencing of the Brassica species Brassica rapa, Brassica oleracea and others demonstrated that WGT plays an important role in the speciation and morphotype diversification of Brassica plants. Comparative genomic analysis based on the genome sequences of B. rapa and A. thaliana clearly identified the WGT event and further demonstrated that the translocated Proto-Calepine Karyotype (tPCK, n=7) was the diploid ancestor of the three subgenomes in B. rapa. Following WGT, subsequent extensive genome fractionation, block reshuffling and chromosome reduction accompanied by paleocentromere descent from the three tPCK subgenomes during the rediploidization process produced stable diploid species. Genomic rearrangement of the diploid species and their hybridization then contributed to Brassica speciation. The subgenome dominance effect and biased gene retention, such as the over-retention of auxin-related genes after WGT, promoted functional gene evolution and thus propelled the expansion of rich morphotypes in the Brassica species. In conclusion, the WGT event initiated subsequent genomic and gene-level evolution, which further drove Brassica speciation and created rich morphotypes in each species.

Genome-wide characterization of microRNA in foxtail millet (Setaria italica)

BACKGROUND

MicroRNAs (miRNAs) are a class of short non-coding, endogenous RNAs that play key roles in many biological processes in both animals and plants. Although many miRNAs have been identified in a large number of organisms, the miRNAs in foxtail millet (Setaria italica) have, until now, been poorly understood.

RESULTS

In this study, two replicate small RNA libraries from foxtail millet shoots were sequenced, and 40 million reads representing over 10 million unique sequences were generated. We identified 43 known miRNAs, 172 novel miRNAs and 2 mirtron precursor candidates in foxtail millet. Some miRNA*s of the known and novel miRNAs were detected as well. Further, eight novel miRNAs were validated by stem-loop RT-PCR. Potential targets of the foxtail millet miRNAs were predicted based on our strict criteria. Of the predicted target genes, 79% (351) had functional annotations in InterPro and GO analyses, indicating the targets of the miRNAs were involved in a wide range of regulatory functions and some specific biological processes. A total of 69 pairs of syntenic miRNA precursors that were conserved between foxtail millet and sorghum were found. Additionally, stem-loop RT-PCR was conducted to confirm the tissue-specific expression of some miRNAs in the four tissues identified by deep-sequencing.

CONCLUSIONS

We predicted, for the first time, 215 miRNAs and 447 miRNA targets in foxtail millet at a genome-wide level. The precursors, expression levels, miRNA* sequences, target functions, conservation, and evolution of miRNAs we identified were investigated. Some of the novel foxtail millet miRNAs and miRNA targets were validated experimentally.

Two evolutionarily distinct classes of paleopolyploidy

Whole genome duplications (WGDs) occurred in the distant evolutionary history of many lineages and are particularly frequent in the flowering plant lineages. Following paleopolyploidization in plants, most duplicated genes are deleted by intrachromosomal recombination, a process referred to as fractionation. In the examples studied so far, genes are disproportionately lost from one of the parental subgenomes (biased fractionation) and the subgenome having lost the lowest number of genes is more expressed (genome dominance). In the present study, we analyzed the pattern of gene deletion and gene expression following the most recent WGD in banana (alpha event) and extended our analyses to seven other sequenced plant genomes: poplar, soybean, medicago, arabidopsis, sorghum, brassica, and maize. We propose a new class of ancient WGD, with Musa (alpha), poplar, and soybean as members, where genes are both deleted and expressed to an equal extent (unbiased fractionation and genome equivalence). We suggest that WGDs with genome dominance and biased fractionation (Class I) may result from ancient allotetraploidies, while WGDs without genome dominance or biased fractionation (Class II) may result from ancient autotetraploidies.

The maize (Zea mays L.) AUXIN/INDOLE-3-ACETIC ACID gene family: phylogeny, synteny, and unique root-type and tissue-specific expression patterns during development

The plant hormone auxin plays a key role in the coordination of many aspects of growth and development. AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) genes encode instable primary auxin responsive regulators of plant development that display a protein structure with four characteristic domains. In the present study, a comprehensive analysis of the 34 members of the maize Aux/IAA gene family was performed. Phylogenetic reconstructions revealed two classes of Aux/IAA proteins that can be distinguished by alterations in their domain III. Seven pairs of paralogous maize Aux/IAA proteins were discovered. Comprehensive root-type and tissue-specific expression profiling revealed unique expression patterns of the diverse members of the gene family. Remarkably, five of seven pairs of paralogous genes displayed highly correlated expression patterns in roots. All but one (ZmIAA23) tested maize Aux/IAA genes were auxin inducible, displaying two types of auxin induction within three hours of treatment. Moreover, 51 of 55 (93%) differential Aux/IAA expression patterns between different root-types followed the expression tendency: crown roots > seminal roots > primary roots > lateral roots. This pattern might imply root-type-specific regulation of Aux/IAA transcript abundance. In summary, the detailed analysis of the maize Aux/IAA gene family provides novel insights in the evolution and developmental regulation and thus the function of these genes in different root-types and tissues.

Practical aliquoting of flowering plant genomes

We pose the problem of dissecting an ancient polyploid genome into its constituent subgenomes despite fragmentation and noise caused by genome rearrangements and fractionation of multi-copy genes. We formulate this in terms of decomposition into "defective" k-partite graphs, distinguished by location within the genome. We devise and implement a clustering heuristic for solving realistic instances of the problem. An unusual focus of our method is the focus on prioritizing gene density or lack of gaps in the assembly of fragments into larger regions, rather than maximizing the number of genes. We validate the method against the grape genome in which the ancient core eudicot triplication is readily detectible and is already well known. We then analyze the tomato genome, whose proposed status as a descendant of a more recent Solanum hexaploid is controversial, and confirm this proposal. The solution reveals unexpected information about the evolution of the tomato.

The circadian clock-associated gene gigantea1 affects maize developmental transitions

The circadian clock is an internal timing mechanism that allows plants to make developmental decisions in accordance with environmental conditions. In model plants, circadian clock-associated gigantea (gi) genes are directly involved in control of growth and developmental transitions. The maize gigantea1 (gi1) gene is the more highly expressed of the two gi homeologs, and its function is uncharacterized. To understand the role of gi1 in the regulatory networks of the maize circadian clock system, gi1 mutants were evaluated for changes in flowering time, phase change and growth control. When grown in long-day (LD) photoperiods, gi1 mutants flowered earlier than non-mutant plants, but this difference was not apparent in short-day (SD) photoperiods. Therefore, gi1 participates in a pathway that suppresses flowering in LD photoperiods, but not in SD. Part of the underlying cause of early flowering was up-regulated expression of the FT-like floral activator gene zea mays centroradialis8 (zcn8) and the CONSTANS-like flowering regulatory gene constans of zea mays1 (conz1). gi1 mutants also underwent vegetative phase change earlier and grew taller than non-mutant plants. These findings indicate gi1 has a repressive function in multiple regulatory pathways that govern maize growth and development.

Gene dosage effects: nonlinearities, genetic interactions, and dosage compensation

High-throughput genomic analyses have shown that many mutations, including loss-of-function (LOF) mutations, are present in diseased as well as in healthy individuals. Gene dosage effects due to deletions, duplications, and LOF mutations provide avenues to explore oligo- and multigenic inheritance. Here, we focus on several mechanisms that mediate gene dosage effects and analyze biochemical interactions among multiple gene products that are sources of nonlinear relations connecting genotypes and phenotypes. We also explore potential mechanisms that compensate for gene dosage effects. Understanding these issues is critical to understanding why an individual bearing a few damaging mutations can be severely diseased, whereas others harboring tens of potentially deleterious mutations can appear quite healthy.

Fractionation, rearrangement and subgenome dominance

Motivation: Fractionation is arguably the greatest cause of gene order disruption following whole genome duplication, causing severe biases in chromosome rearrangement-based estimates of evolutionary divergence.

Results: We show how to correct for this bias almost entirely by means of a ‘consolidation’ algorithm for detecting and suitably transforming identifiable regions of fractionation. We characterize the process of fractionation and the performance of the algorithm through realistic simulations. We apply our method to a number of core eudicot genomes, we and by studying the fractionation regions detected, are able to address topical issues in polyploid evolution.

Availability and implementation: Code for the consolidation algorithm, and sample data, is available at: http://137.122.149.195/Software/Fractionation/fractionation.html

Contact:sankoff@uottawa.ca

Maize gene atlas developed by RNA sequencing and comparative evaluation of transcriptomes based on RNA sequencing and microarrays

Transcriptome analysis is a valuable tool for identification and characterization of genes and pathways underlying plant growth and development. We previously published a microarray-based maize gene atlas from the analysis of 60 unique spatially and temporally separated tissues from 11 maize organs [1]. To enhance the coverage and resolution of the maize gene atlas, we have analyzed 18 selected tissues representing five organs using RNA sequencing (RNA-Seq). For a direct comparison of the two methodologies, the same RNA samples originally used for our microarray-based atlas were evaluated using RNA-Seq. Both technologies produced similar transcriptome profiles as evident from high Pearson's correlation statistics ranging from 0.70 to 0.83, and from nearly identical clustering of the tissues. RNA-Seq provided enhanced coverage of the transcriptome, with 82.1% of the filtered maize genes detected as expressed in at least one tissue by RNA-Seq compared to only 56.5% detected by microarrays. Further, from the set of 465 maize genes that have been historically well characterized by mutant analysis, 427 show significant expression in at least one tissue by RNA-Seq compared to 390 by microarray analysis. RNA-Seq provided higher resolution for identifying tissue-specific expression as well as for distinguishing the expression profiles of closely related paralogs as compared to microarray-derived profiles. Co-expression analysis derived from the microarray and RNA-Seq data revealed that broadly similar networks result from both platforms, and that co-expression estimates are stable even when constructed from mixed data including both RNA-Seq and microarray expression data. The RNA-Seq information provides a useful complement to the microarray-based maize gene atlas and helps to further understand the dynamics of transcription during maize development.

The pattern of Phosphate transporter 1 genes evolutionary divergence in Glycine max L

BACKGROUND

The Phosphate transporter 1 (PHT1) gene family has crucial roles in phosphate uptake, translocation, remobilization, and optimization of metabolic processes using of Pi. Gene duplications expand the size of gene families, and subfunctionalization of paralog gene pairs is a predominant tendency after gene duplications. To date, experimental evidence for the evolutionary relationships among different paralog gene pairs of a given gene family in soybean is limited.

RESULTS

All potential Phosphate transporter 1 genes in Glycine max L. (GmPHT1) were systematically analyzed using both bioinformatics and experimentation. The soybean PHT1 genes originated from four distinct ancestors prior to the Gamma WGT and formed 7 paralog gene pairs and a singleton gene. Six of the paralog gene pairs underwent subfunctionalization, and while GmPHT1;4 paralog gene experienced pseudogenization. Examination of long-term evolutionary changes, six GmPHT1 paralog gene pairs diverged at multiple levels, in aspects of spatio-temporal expression patterns and/or quanta, phosphates affinity properties, subcellular localization, and responses to phosphorus stress.

CONCLUSIONS

These characterized divergences occurred in tissue- and/or development-specific modes, or conditional modes. Moreover, they have synergistically shaped the evolutionary rate of GmPHT1 family, as well as maintained phosphorus homeostasis at cells and in the whole plant.

Genomic distribution of maize facultative heterochromatin marked by trimethylation of H3K27

Trimethylation of histone H3 Lys-27 (H3K27me3) plays a critical role in regulating gene expression during plant and animal development. We characterized the genome-wide distribution of H3K27me3 in five developmentally distinct tissues in maize (Zea mays) plants of two genetic backgrounds, B73 and Mo17. There were more substantial differences in the genome-wide profile of H3K27me3 between different tissues than between the two genotypes. The tissue-specific patterns of H3K27me3 were often associated with differences in gene expression among the tissues and most of the imprinted genes that are expressed solely from the paternal allele in endosperm are targets of H3K27me3. A comparison of the H3K27me3 targets in rice (Oryza sativa), maize, and Arabidopsis thaliana provided evidence for conservation of the H3K27me3 targets among plant species. However, there was limited evidence for conserved targeting of H3K27me3 in the two maize subgenomes derived from whole-genome duplication, suggesting the potential for subfunctionalization of chromatin regulation of paralogs. Genomic profiling of H3K27me3 in loss-of-function mutant lines for Maize Enhancer of zeste-like2 (Mez2) and Mez3, two of the three putative H3K27me3 methyltransferases present in the maize genome, suggested partial redundancy of this gene family for maintaining H3K27me3 patterns. Only a portion of the targets of H3K27me3 required Mez2 and/or Mez3, and there was limited evidence for functional consequences of H3K27me3 at these targets.

Watching the grin fade: tracing the effects of polyploidy on different evolutionary time scales

Polyploidy, or whole-genome duplication (WGD), is a recurrent mutation both in cell lineages and over evolutionary time. By globally changing the relationship between gene copy number and other cellular entities, it can induce dramatic changes at the cellular and phenotypic level. Perhaps surprisingly, then, the insights that these events can bring to understanding other cellular features are not as well appreciated as they could be. In this review, we draw on examples of polyploidy from animals, plants and yeast to explore how investigations of polyploid cells have improved our understanding of the cell cycle, biological network complexity, metabolic phenotypes and tumor biology. We argue that the study of polyploidy across organisms, cell types, and time scales serves not only as a window into basic cell biology, but also as a basis for a predictive biology with applications ranging from crop improvement to treating cancer.

Patterns of nucleotide asymmetries in plant and animal genomes

Symmetry in biology provides many intriguing puzzles to the scientist's mind. Chargaff's second parity rule states a symmetric distribution of oligonucleotides within a single strand of double-stranded DNA. While this rule has been verified in a wide range of microbial genomes, it still awaits explanation. In our study, we inquired into patterns of mono- and trinucleotide intra-strand parity in complex plant genomic sequences that became available during the last few years, and compared these to equally complex animal genomes. The degree and patterns of deviation from Chargaff's second rule were different between plant and animal species. We observed a universal inter-chromosomal homogeneity of mononucleotide skews in coding sequences of plant chromosomes, while the base composition of animal coding sequences differed between chromosomes even within a single species. We also found differences in the base composition of dicot introns in comparison to those of monocots. These genome-wide patterns were limited to genic regions and were not encountered in inter-genic sequences. We discuss the implications of our findings in relation to hypotheses about functional correlations of intra-strand parity which have hitherto been put forward. Furthermore, we propose more recent polyploidization and subsequent homogenization of homoeologues as a possible reason for more homogeneous skew patterns in plants.

The dynamics of LTR retrotransposon accumulation across 25 million years of panicoid grass evolution

Sample sequence analysis was employed to investigate the repetitive DNAs that were most responsible for the evolved variation in genome content across seven panicoid grasses with >5-fold variation in genome size and different histories of polyploidy. In all cases, the most abundant repeats were LTR retrotransposons, but the particular families that had become dominant were found to be different in the Pennisetum, Saccharum, Sorghum and Zea lineages. One element family, Huck, has been very active in all of the studied species over the last few million years. This suggests the transmittal of an active or quiescent autonomous set of Huck elements to this lineage at the founding of the panicoids. Similarly, independent recent activity of Ji and Opie elements in Zea and of Leviathan elements in Sorghum and Saccharum species suggests that members of these families with exceptional activation potential were present in the genome(s) of the founders of these lineages. In a detailed analysis of the Zea lineage, the combined action of several families of LTR retrotransposons were observed to have approximately doubled the genome size of Zea luxurians relative to Zea mays and Zea diploperennis in just the last few million years. One of the LTR retrotransposon amplification bursts in Zea may have been initiated by polyploidy, but the great majority of transposable element activations are not. Instead, the results suggest random activation of a few or many LTR retrotransposons families in particular lineages over evolutionary time, with some families especially prone to future activation and hyper-amplification.

Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol

Switchgrass (Panicum virgatum L.) is a perennial grass that has been designated as an herbaceous model biofuel crop for the United States of America. To facilitate accelerated breeding programs of switchgrass, we developed both an association panel and linkage populations for genome-wide association study (GWAS) and genomic selection (GS). All of the 840 individuals were then genotyped using genotyping by sequencing (GBS), generating 350 GB of sequence in total. As a highly heterozygous polyploid (tetraploid and octoploid) species lacking a reference genome, switchgrass is highly intractable with earlier methodologies of single nucleotide polymorphism (SNP) discovery. To access the genetic diversity of species like switchgrass, we developed a SNP discovery pipeline based on a network approach called the Universal Network-Enabled Analysis Kit (UNEAK). Complexities that hinder single nucleotide polymorphism discovery, such as repeats, paralogs, and sequencing errors, are easily resolved with UNEAK. Here, 1.2 million putative SNPs were discovered in a diverse collection of primarily upland, northern-adapted switchgrass populations. Further analysis of this data set revealed the fundamentally diploid nature of tetraploid switchgrass. Taking advantage of the high conservation of genome structure between switchgrass and foxtail millet (Setaria italica (L.) P. Beauv.), two parent-specific, synteny-based, ultra high-density linkage maps containing a total of 88,217 SNPs were constructed. Also, our results showed clear patterns of isolation-by-distance and isolation-by-ploidy in natural populations of switchgrass. Phylogenetic analysis supported a general south-to-north migration path of switchgrass. In addition, this analysis suggested that upland tetraploid arose from upland octoploid. All together, this study provides unparalleled insights into the diversity, genomic complexity, population structure, phylogeny, phylogeography, ploidy, and evolutionary dynamics of switchgrass.

Recent advances in sequencing technologies have enabled large-scale surveys of genetic diversity in model species with a wholly or partly sequenced reference genome. However, thousands of key species, which are essential for food, health, energy, and ecology, do not have reference genomes. To accelerate their breeding cycle via marker assisted selection, high-throughput genotyping is required for these valuable species, in spite of the absence of reference genomes. Based on genotyping by sequencing (GBS) technology, we developed a new single nucleotide polymorphism (SNP) discovery protocol, the Universal Network-Enabled Analysis Kit (UNEAK), which can be widely used in any species, regardless of genome complexity or the availability of a reference genome. Here we test this protocol on switchgrass, currently the prime energy crop species in the United States of America. In addition to the discovery of over a million SNPs and construction of high-density linkage maps, we provide novel insights into the genome complexity, ploidy, phylogeny, and evolution of switchgrass. This is only the beginning: we believe UNEAK offers the key to the exploration and exploitation of the genetic diversity of thousands of non-model species.

Gene body methylation is conserved between plant orthologs and is of evolutionary consequence

DNA methylation is a common feature of eukaryotic genomes and is especially common in noncoding regions of plants. Protein coding regions of plants are often methylated also, but the extent, function, and evolutionary consequences of gene body methylation remain unclear. Here we investigate gene body methylation using an explicit comparative evolutionary approach. We generated bisulfite sequencing data from two tissues of Brachypodium distachyon and compared genic methylation patterns to those of rice (Oryza sativa ssp. japonica). Gene body methylation was strongly conserved between orthologs of the two species and affected a biased subset of long, slowly evolving genes. Because gene body methylation is conserved over evolutionary time, it shapes important features of plant genome evolution, such as the bimodality of G+C content among grass genes. Our results superficially contradict previous observations of high cytosine methylation polymorphism within Arabidopsis thaliana genes, but reanalyses of these data are consistent with conservation of methylation within gene regions. Overall, our results indicate that the methylation level is a long-term property of individual genes and therefore of evolutionary consequence.

Comparative analysis of syntenic genes in grass genomes reveals accelerated rates of gene structure and coding sequence evolution in polyploid wheat

Cycles of whole-genome duplication (WGD) and diploidization are hallmarks of eukaryotic genome evolution and speciation. Polyploid wheat (Triticum aestivum) has had a massive increase in genome size largely due to recent WGDs. How these processes may impact the dynamics of gene evolution was studied by comparing the patterns of gene structure changes, alternative splicing (AS), and codon substitution rates among wheat and model grass genomes. In orthologous gene sets, significantly more acquired and lost exonic sequences were detected in wheat than in model grasses. In wheat, 35% of these gene structure rearrangements resulted in frame-shift mutations and premature termination codons. An increased codon mutation rate in the wheat lineage compared with Brachypodium distachyon was found for 17% of orthologs. The discovery of premature termination codons in 38% of expressed genes was consistent with ongoing pseudogenization of the wheat genome. The rates of AS within the individual wheat subgenomes (21%-25%) were similar to diploid plants. However, we uncovered a high level of AS pattern divergence between the duplicated homeologous copies of genes. Our results are consistent with the accelerated accumulation of AS isoforms, nonsynonymous mutations, and gene structure rearrangements in the wheat lineage, likely due to genetic redundancy created by WGDs. Whereas these processes mostly contribute to the degeneration of a duplicated genome and its diploidization, they have the potential to facilitate the origin of new functional variations, which, upon selection in the evolutionary lineage, may play an important role in the origin of novel traits.

Accessing complex crop genomes with next-generation sequencing

Many important crop species have genomes originating from ancestral or recent polyploidisation events. Multiple homoeologous gene copies, chromosomal rearrangements and amplification of repetitive DNA within large and complex crop genomes can considerably complicate genome analysis and gene discovery by conventional, forward genetics approaches. On the other hand, ongoing technological advances in molecular genetics and genomics today offer unprecedented opportunities to analyse and access even more recalcitrant genomes. In this review, we describe next-generation sequencing and data analysis techniques that vastly improve our ability to dissect and mine genomes for causal genes underlying key traits and allelic variation of interest to breeders. We focus primarily on wheat and oilseed rape, two leading examples of major polyploid crop genomes whose size or complexity present different, significant challenges. In both cases, the latest DNA sequencing technologies, applied using quite different approaches, have enabled considerable progress towards unravelling the respective genomes. Our ability to discover the extent and distribution of genetic diversity in crop gene pools, and its relationship to yield and quality-related traits, is swiftly gathering momentum as DNA sequencing and the bioinformatic tools to deal with growing quantities of genomic data continue to develop. In the coming decade, genomic and transcriptomic sequencing, discovery and high-throughput screening of single nucleotide polymorphisms, presence-absence variations and other structural chromosomal variants in diverse germplasm collections will give detailed insight into the origins, domestication and available trait-relevant variation of polyploid crops, in the process facilitating novel approaches and possibilities for genomics-assisted breeding.

Automated conserved non-coding sequence (CNS) discovery reveals differences in gene content and promoter evolution among grasses

Conserved non-coding sequences (CNS) are islands of non-coding sequence that, like protein coding exons, show less divergence in sequence between related species than functionless DNA. Several CNSs have been demonstrated experimentally to function as cis-regulatory regions. However, the specific functions of most CNSs remain unknown. Previous searches for CNS in plants have either anchored on exons and only identified nearby sequences or required years of painstaking manual annotation. Here we present an open source tool that can accurately identify CNSs between any two related species with sequenced genomes, including both those immediately adjacent to exons and distal sequences separated by >12 kb of non-coding sequence. We have used this tool to characterize new motifs, associate CNSs with additional functions, and identify previously undetected genes encoding RNA and protein in the genomes of five grass species. We provide a list of 15,363 orthologous CNSs conserved across all grasses tested. We were also able to identify regulatory sequences present in the common ancestor of grasses that have been lost in one or more extant grass lineages. Lists of orthologous gene pairs and associated CNSs are provided for reference inbred lines of arabidopsis, Japonica rice, foxtail millet, sorghum, brachypodium, and maize.

Switchgrass genomic diversity, ploidy, and evolution: novel insights from a network-based SNP discovery protocol

Switchgrass (Panicum virgatum L.) is a perennial grass that has been designated as an herbaceous model biofuel crop for the United States of America. To facilitate accelerated breeding programs of switchgrass, we developed both an association panel and linkage populations for genome-wide association study (GWAS) and genomic selection (GS). All of the 840 individuals were then genotyped using genotyping by sequencing (GBS), generating 350 GB of sequence in total. As a highly heterozygous polyploid (tetraploid and octoploid) species lacking a reference genome, switchgrass is highly intractable with earlier methodologies of single nucleotide polymorphism (SNP) discovery. To access the genetic diversity of species like switchgrass, we developed a SNP discovery pipeline based on a network approach called the Universal Network-Enabled Analysis Kit (UNEAK). Complexities that hinder single nucleotide polymorphism discovery, such as repeats, paralogs, and sequencing errors, are easily resolved with UNEAK. Here, 1.2 million putative SNPs were discovered in a diverse collection of primarily upland, northern-adapted switchgrass populations. Further analysis of this data set revealed the fundamentally diploid nature of tetraploid switchgrass. Taking advantage of the high conservation of genome structure between switchgrass and foxtail millet (Setaria italica (L.) P. Beauv.), two parent-specific, synteny-based, ultra high-density linkage maps containing a total of 88,217 SNPs were constructed. Also, our results showed clear patterns of isolation-by-distance and isolation-by-ploidy in natural populations of switchgrass. Phylogenetic analysis supported a general south-to-north migration path of switchgrass. In addition, this analysis suggested that upland tetraploid arose from upland octoploid. All together, this study provides unparalleled insights into the diversity, genomic complexity, population structure, phylogeny, phylogeography, ploidy, and evolutionary dynamics of switchgrass.

Analysis of the bread wheat genome using whole-genome shotgun sequencing

Bread wheat (Triticum aestivum) is a globally important crop, accounting for 20 per cent of the calories consumed by humans. Major efforts are underway worldwide to increase wheat production by extending genetic diversity and analysing key traits, and genomic resources can accelerate progress. But so far the very large size and polyploid complexity of the bread wheat genome have been substantial barriers to genome analysis. Here we report the sequencing of its large, 17-gigabase-pair, hexaploid genome using 454 pyrosequencing, and comparison of this with the sequences of diploid ancestral and progenitor genomes. We identified between 94,000 and 96,000 genes, and assigned two-thirds to the three component genomes (A, B and D) of hexaploid wheat. High-resolution synteny maps identified many small disruptions to conserved gene order. We show that the hexaploid genome is highly dynamic, with significant loss of gene family members on polyploidization and domestication, and an abundance of gene fragments. Several classes of genes involved in energy harvesting, metabolism and growth are among expanded gene families that could be associated with crop productivity. Our analyses, coupled with the identification of extensive genetic variation, provide a resource for accelerating gene discovery and improving this major crop.

Centromere retention and loss during the descent of maize from a tetraploid ancestor

Although centromere function is highly conserved in eukaryotes, centromere sequences are highly variable. Only a few centromeres have been sequenced in higher eukaryotes because of their repetitive nature, thus hindering study of their structure and evolution. Conserved single-copy sequences in pericentromeres (CSCPs) of sorghum and maize were found to be diagnostic characteristics of adjacent centromeres. By analyzing comparative map data and CSCP sequences of sorghum, maize, and rice, the major evolutionary events related to centromere dynamics were discovered for the maize lineage after its divergence from a common ancestor with sorghum. (i) Remnants of ancient CSCP regions were found for the 10 lost ancestral centromeres, indicating that two ancient homeologous chromosome pairs did not contribute any centromeres to the current maize genome, whereas two other pairs contributed both of their centromeres. (ii) Five cases of long-distance, intrachromosome movement of CSCPs were detected in the retained centromeres, with inversion the major process involved. (iii) The 12 major chromosomal rearrangements that led to maize chromosome number reduction from 20 to 10 were uncovered. (iv) In addition to whole chromosome insertion near (but not always into) other centromeres, translocation and fusion were found to be important mechanisms underlying grass chromosome number reduction. (v) Comparison of chromosome structures confirms the polyploid event that led to the tetraploid ancestor of modern maize.

Downy mildew resistance induced by Trichoderma harzianum T39 in susceptible grapevines partially mimics transcriptional changes of resistant genotypes

BACKGROUND

Downy mildew, caused by Plasmopara viticola, is one of the most severe diseases of grapevine and is commonly controlled by fungicide treatments. The beneficial microorganism Trichoderma harzianum T39 (T39) can induce resistance to downy mildew, although the molecular events associated with this process have not yet been elucidated in grapevine. A next generation RNA sequencing (RNA-Seq) approach was used to study global transcriptional changes associated with resistance induced by T39 in Vitis vinifera Pinot Noir leaves. The long-term aim was to develop strategies to optimize the use of this agent for downy mildew control.

RESULTS

More than 14.8 million paired-end reads were obtained for each biological replicate of T39-treated and control leaf samples collected before and 24 h after P. viticola inoculation. RNA-Seq analysis resulted in the identification of 7,024 differentially expressed genes, highlighting the complex transcriptional reprogramming of grapevine leaves during resistance induction and in response to pathogen inoculation. Our data show that T39 has a dual effect: it directly modulates genes related to the microbial recognition machinery, and it enhances the expression of defence-related processes after pathogen inoculation. Whereas several genes were commonly affected by P. viticola in control and T39-treated plants, opposing modulation of genes related to responses to stress and protein metabolism was found. T39-induced resistance partially inhibited some disease-related processes and specifically activated defence responses after P. viticola inoculation, causing a significant reduction of downy mildew symptoms.

CONCLUSIONS

The global transcriptional analysis revealed that defence processes known to be implicated in the reaction of resistant genotypes to downy mildew were partially activated by T39-induced resistance in susceptible grapevines. Genes identified in this work are an important source of markers for selecting novel resistance inducers and for the analysis of environmental conditions that might affect induced resistance mechanisms.

Insights from paleogenomic and population studies into the consequences of dosage sensitive gene expression in plants

Classical studies of plant phenotypes of individuals with whole or partial genome dosage changes led to the concept of genomic balance. Subsequent studies of gene expression in ploidy and aneuploidy series showed a greater number of modulations in aneuploid plants than with whole genome changes leading to the idea that gene expression processes were modulated by stoichiometric changes of interacting regulatory factors. Recent studies of genomic sequences and copy number variants in populations reveal different fates of duplicate genes depending on whole genome or segmental duplication. Following polyploidy formation, members of macromolecular complexes persist in the evolutionary lineage longer than random genes and a complementary pattern is found for segmental duplications in that there is an underrepresentation of members of macromolecular complexes. These and other studies described suggest there are negative fitness consequences when an imbalance occurs for members of macromolecular complexes including regulatory functions.

Changes in genome content generated via segregation of non-allelic homologs

A careful analysis of two maize recombinant inbred lines (RILs) relative to their inbred parents revealed the presence of several hundred apparently de novo copy number variants (CNVs). These changes in genome content were validated via both PCR and whole exome-array capture-and-sequencing experiments. One hundred and eighty-five genomic regions, which overlap with 38 high-confidence genes, exhibited apparently de novo copy number variation (CNV) in these two RILs and in many instances the same apparently de novo CNV events were observed in multiple RILs. Further analyses revealed that these recurrent apparently de novo CNVs were caused by segregation of single-copy homologous sequences that are located in non-allelic positions in the two parental inbred lines. F(1) individuals derived from these inbred lines will be hemizygous for each of these non-allelic homologs but RIL genotypes will contain these sequences at zero, one or two genomic loci. Hence, the segregation of non-allelic homologs may contribute to transgressive segregation. Indeed, statistical associations between phenotypic quantitative trait loci and genomic losses were observed for two of 14 tested pairs of non-allelic homologs.

The banana (Musa acuminata) genome and the evolution of monocotyledonous plants

Bananas (Musa spp.), including dessert and cooking types, are giant perennial monocotyledonous herbs of the order Zingiberales, a sister group to the well-studied Poales, which include cereals. Bananas are vital for food security in many tropical and subtropical countries and the most popular fruit in industrialized countries. The Musa domestication process started some 7,000 years ago in Southeast Asia. It involved hybridizations between diverse species and subspecies, fostered by human migrations, and selection of diploid and triploid seedless, parthenocarpic hybrids thereafter widely dispersed by vegetative propagation. Half of the current production relies on somaclones derived from a single triploid genotype (Cavendish). Pests and diseases have gradually become adapted, representing an imminent danger for global banana production. Here we describe the draft sequence of the 523-megabase genome of a Musa acuminata doubled-haploid genotype, providing a crucial stepping-stone for genetic improvement of banana. We detected three rounds of whole-genome duplications in the Musa lineage, independently of those previously described in the Poales lineage and the one we detected in the Arecales lineage. This first monocotyledon high-continuity whole-genome sequence reported outside Poales represents an essential bridge for comparative genome analysis in plants. As such, it clarifies commelinid-monocotyledon phylogenetic relationships, reveals Poaceae-specific features and has led to the discovery of conserved non-coding sequences predating monocotyledon-eudicotyledon divergence.

Homoploid hybrid origin of Yucca gloriosa: intersectional hybrid speciation in Yucca (Agavoideae, Asparagaceae)

There is a growing appreciation for the importance of hybrid speciation in angiosperm evolution. Here, we show that Yucca gloriosa (Asparagaceae: Agavoideae) is the product of intersectional hybridization between Y. aloifolia and Y. filamentosa. These species, all named by Carl Linnaeus, exist in sympatry along the southeastern Atlantic coast of the United States. Yucca gloriosa was found to share a chloroplast haplotype with Y. aloifolia in all populations sampled. In contrast, nuclear gene-based microsatellite markers in Y. gloriosa are shared with both parents. The hybrid origin of Y. gloriosa is supported by multilocus analyses of the nuclear microsatellite markers including principal coordinates analysis (PCO), maximum-likelihood hybrid index scoring (HINDEX), and Bayesian cluster analysis (STRUCTURE). The putative parental species share only one allele at a single locus, suggesting there is little to no introgressive gene flow occurring between these species and Y. gloriosa. At the same time, diagnostic markers are segregating in Y. gloriosa populations. Lack of variation in the chloroplast of Y. aloifolia, the putative maternal parent, makes it difficult to rule out multiple hybrid origins of Y. gloriosa, but allelic variation at nuclear loci can be explained by a single hybrid origin of Y. gloriosa. Overall, these data provide strong support for the homoploid hybrid origin of Y. gloriosa.

Gene balance hypothesis: connecting issues of dosage sensitivity across biological disciplines

We summarize, in this review, the evidence that genomic balance influences gene expression, quantitative traits, dosage compensation, aneuploid syndromes, population dynamics of copy number variants and differential evolutionary fate of genes after partial or whole-genome duplication. Gene balance effects are hypothesized to result from stoichiometric differences among members of macromolecular complexes, the interactome, and signaling pathways. The implications of gene balance are discussed.

Reliable transient transformation of intact maize leaf cells for functional genomics and experimental study

Maize (Zea mays) transformation routinely produces stable transgenic lines essential for functional genomics; however, transient expression of target proteins in maize cells is not yet routine. Such techniques are critical for rapid testing of transgene constructs and for experimental studies. Here, we report bombardment methods that depend on leaf developmental stage and result in successful expression with broad applications. Fluorescent marker genes were constructed and bombarded into five developmental regions in a growing maize leaf. Expression efficiency was highest in the basal-most 3 cm above the ligule of an approximately 50-cm growing adult leaf. Straightforward dissection procedures provide access to the receptive leaf regions, increasing efficiency from less than one transformant per cm(2) to over 21 transformants per cm(2). Successful expression was routine for proteins from full genomic sequences driven by native regulatory regions and from complementary DNA sequences driven by the constitutive maize polyubiquitin promoter and a heterologous terminator. Four tested fusion proteins, maize PROTEIN DISULFIDE ISOMERASE-Yellow Fluorescent Protein, GLOSSY8a-monomeric Red Fluorescent Protein and maize XYLOSYLTRANSFERASE, and maize Rho-of-Plants7-monomeric Teal Fluorescent Protein, localized as predicted in the endoplasmic reticulum, Golgi, and plasma membrane, respectively. Localization patterns were similar between transient and stable modes of expression, and cotransformation was equally successful. Coexpression was also demonstrated by transiently transforming cells in a stable line expressing a second marker protein, thus increasing the utility of a single stable transformant. Given the ease of dissection procedures, this method replaces heterologous expression assays with a more direct, native, and informative system, and the techniques will be useful for localization, colocalization, and functional studies.

Very low rate of gene conversion in the yeast genome

Gene duplication is a major driver of organismal adaptation and evolution and plays an important role in multiple human diseases. Whole-genome analyses have shown similar and high rates of gene duplication across a variety of eukaryotic species. Most of these studies, however, did not address the possible impact of interlocus gene conversion (IGC) on the evolution of gene duplicates. Because IGC homogenizes pairs of duplicates, widespread conversion would cause gene duplication events that happened long ago to appear more recent, resulting in artificially high estimates of duplication rates. Although the majority of genome-wide studies (including in the budding yeast Saccharomyces cerevisiae [Scer]) point to levels of IGC between paralogs ranging from 2% to 18%, Gao and Innan (Gao LZ, Innan H. 2004. Very low gene duplication rate in the yeast genome. Science 306:1367-1370.) found that gene conversion in yeast affected >80% of paralog pairs. If conversion rates really are this high, it would imply that the rate of gene duplication in eukaryotes is much lower than previously reported. In this work, we apply four different methodologies-including one approach that closely mirrors Gao and Innan's method-to estimate the level of IGC in Scer. Our analyses point to a maximum conversion level of 13% between paralogs in this species, in close agreement with most estimates of IGC in eukaryotes. We also show that the exceedingly high levels of conversion found previously derive from application of an accurate method to an inappropriate data set. In conclusion, our work provides the most striking evidence to date supporting the reduced incidence of IGC among Scer paralogs and sets up a framework for future analyses in other eukaryotes.

Decoding plant and animal genome plasticity from differential paleo-evolutionary patterns and processes

Continuing advances in genome sequencing technologies and computational methods for comparative genomics currently allow inferring the evolutionary history of entire plant and animal genomes. Based on the comparison of the plant and animal genome paleohistory, major differences are unveiled in 1) evolutionary mechanisms (i.e., polyploidization versus diploidization processes), 2) genome conservation (i.e., coding versus noncoding sequence maintenance), and 3) modern genome architecture (i.e., genome organization including repeats expansion versus contraction phenomena). This article discusses how extant animal and plant genomes are the result of inherently different rates and modes of genome evolution resulting in relatively stable animal and much more dynamic and plastic plant genomes.

Genic and nongenic contributions to natural variation of quantitative traits in maize

The complex genomes of many economically important crops present tremendous challenges to understand the genetic control of many quantitative traits with great importance in crop production, adaptation, and evolution. Advances in genomic technology need to be integrated with strategic genetic design and novel perspectives to break new ground. Complementary to individual-gene–targeted research, which remains challenging, a global assessment of the genomic distribution of trait-associated SNPs (TASs) discovered from genome scans of quantitative traits can provide insights into the genetic architecture and contribute to the design of future studies. Here we report the first systematic tabulation of the relative contribution of different genomic regions to quantitative trait variation in maize. We found that TASs were enriched in the nongenic regions, particularly within a 5-kb window upstream of genes, which highlights the importance of polymorphisms regulating gene expression in shaping the natural variation. Consistent with these findings, TASs collectively explained 44%–59% of the total phenotypic variation across maize quantitative traits, and on average, 79% of the explained variation could be attributed to TASs located in genes or within 5 kb upstream of genes, which together comprise only 13% of the genome. Our findings suggest that efficient, cost-effective genome-wide association studies (GWAS) in species with complex genomes can focus on genic and promoter regions.

Preferential retention of circadian clock genes during diploidization following whole genome triplication in Brassica rapa

Much has been learned about the architecture and function of the circadian clock of Arabidopsis thaliana, a model for plant circadian rhythms. Circadian rhythms contribute to evolutionary fitness, suggesting that circadian rhythmicity may also contribute to agricultural productivity. Therefore, we extend our study of the plant circadian clock to Brassica rapa, an agricultural crop. Since its separation from Arabidopsis, B. rapa has undergone whole genome triplication and subsequent diploidization that has involved considerable gene loss. We find that circadian clock genes are preferentially retained relative to comparison groups of their neighboring genes, a set of randomly chosen genes, and a set of housekeeping genes broadly conserved in eukaryotes. The preferential retention of clock genes is consistent with the gene dosage hypothesis, which predicts preferential retention of highly networked or dose-sensitive genes. Two gene families encoding transcription factors that play important roles in the plant core oscillator--the PSEUDO-RESPONSE REGULATORS, including TIMING OF CAB EXPRESSION1, and the REVEILLE family, including CIRCADIAN CLOCK ASSOCIATED1 and LATE ELONGATED HYPOCOTYL--exhibit preferential retention consistent with the gene dosage hypothesis, but a third gene family, including ZEITLUPE, that encodes F-Box proteins that regulate posttranslational protein stability offers an exception.

Emerging knowledge from genome sequencing of crop species

Extensive insights into the genome composition, organization, and evolution have been gained from the plant genome sequencing and annotation ongoing projects. The analysis of crop genomes provided surprising evidences with important implications in plant origin and evolution: genome duplication, ancestral re-arrangements and unexpected polyploidization events opened new doors to address fundamental questions related to species proliferation, adaptation, and functional modulations. Detailed paleogenomic analysis led to many speculation on how chromosomes have been shaped over time in terms of gene content and order. The completion of the genome sequences of several major crops, prompted to a detailed identification and annotation of transposable elements: new hypothesis related to their composition, chromosomal distribution, insertion models, amplification rate, and evolution patterns are coming up. Availability of full genome sequence of several crop species as well as from many accessions within species is providing new keys for biodiversity exploitation and interpretation. Re-sequencing is enabling high-throughput genotyping to identify a wealth of SNP and afterward to produce haplotype maps necessary to accurately associate molecular variation to phenotype. Conservation genomics is emerging as a powerful tool to explain adaptation, genetic drift, natural selection, hybridization and to estimate genetic variation, fitness and population's viability.

Development of pachytene FISH maps for six maize chromosomes and their integration with other maize maps for insights into genome structure variation

Integrated cytogenetic pachytene fluorescence in situ hybridization (FISH) maps were developed for chromosomes 1, 3, 4, 5, 6, and 8 of maize using restriction fragment length polymorphism marker-selected Sorghum propinquum bacterial artificial chromosomes (BACs) for 19 core bin markers and 4 additional genetic framework loci. Using transgenomic BAC FISH mapping on maize chromosome addition lines of oats, we found that the relative locus position along the pachytene chromosome did not change as a function of total arm length, indicative of uniform axial contraction along the fibers during mid-prophase for tested loci on chromosomes 4 and 5. Additionally, we cytogenetically FISH mapped six loci from chromosome 9 onto their duplicated syntenic regions on chromosomes 1 and 6, which have varying amounts of sequence divergence, using sorghum BACs homologous to the chromosome 9 loci. We found that successful FISH mapping was possible even when the chromosome 9 selective marker had no counterpart in the syntenic block. In total, these 29 FISH-mapped loci were used to create the most extensive pachytene FISH maps to date for these six maize chromosomes. The FISH-mapped loci were then merged into one composite karyotype for direct comparative analysis with the recombination nodule-predicted cytogenetic, genetic linkage, and genomic physical maps using the relative marker positions of the loci on all the maps. Marker colinearity was observed between all pair-wise map comparisons, although marker distribution patterns varied widely in some cases. As expected, we found that the recombination nodule-based predictions most closely resembled the cytogenetic map positions overall. Cytogenetic and linkage map comparisons agreed with previous studies showing a decrease in marker spacing in the peri-centromeric heterochromatin region on the genetic linkage maps. In fact, there was a general trend with most loci mapping closer towards the telomere on the linkage maps than on the cytogenetic maps, regardless of chromosome number or maize inbred line source, with just some of the telomeric loci exempted. Finally and somewhat surprisingly, we observed considerable variation between the relative arm positions of loci when comparing our cytogenetic FISH map to the B73 genomic physical maps, even where comparisons were to a B73-derived cytogenetic map. This variation is more evident between different chromosome arms, but less so within a given arm, ruling out any type of inbred-line dependent global features of linear deoxyribonucleic acid compared with the meiotic fiber organization. This study provides a means for analyzing the maize genome structure by producing new connections for integrating the cytogenetic, linkage, and physical maps of maize.

Biased gene fractionation and dominant gene expression among the subgenomes of Brassica rapa

Polyploidization, both ancient and recent, is frequent among plants. A “two-step theory" was proposed to explain the meso-triplication of the Brassica “A" genome: Brassica rapa. By accurately partitioning of this genome, we observed that genes in the less fractioned subgenome (LF) were dominantly expressed over the genes in more fractioned subgenomes (MFs: MF1 and MF2), while the genes in MF1 were slightly dominantly expressed over the genes in MF2. The results indicated that the dominantly expressed genes tended to be resistant against gene fractionation. By re-sequencing two B. rapa accessions: a vegetable turnip (VT117) and a Rapid Cycling line (L144), we found that genes in LF had less non-synonymous or frameshift mutations than genes in MFs; however mutation rates were not significantly different between MF1 and MF2. The differences in gene expression patterns and on-going gene death among the three subgenomes suggest that “two-step" genome triplication and differential subgenome methylation played important roles in the genome evolution of B. rapa.

Grass microRNA gene paleohistory unveils new insights into gene dosage balance in subgenome partitioning after whole-genome duplication

The recent availability of plant genome sequences, combined with a robust evolutionary scenario of the modern monocot and eudicot karyotypes from their diploid ancestors, offers an opportunity to gain insights into microRNA (miRNA) gene paleohistory in plants. Characterization and comparison of miRNAs and associated protein-coding targets in plants allowed us to unravel (1) contrasted genome conservation patterns of miRNAs in monocots and eudicots after whole-genome duplication (WGD), (2) an ancestral miRNA founder pool in the monocot genomes dating back to 100 million years ago, (3) miRNA subgenome dominance during the post-WGD diploidization process with selective miRNA deletion complemented with possible transposable element-mediated return flows, and (4) the miRNA/target interaction-directed differential loss/retention of miRNAs following the gene dosage balance rule. Together, our data suggest that overretained miRNAs in grass genomes may be implicated in connected gene regulations for stress responses, which is essential for plant adaptation and useful for crop variety innovation.

In silico archeogenomics unveils modern plant genome organisation, regulation and evolution

Increasing access to plant genome sequences as well as high resolution gene-based genetic maps have recently offered the opportunity to compare modern genomes and model their evolutionary history from their reconstructed founder ancestors on an unprecedented scale. In silico paleogenomic data have revealed the evolutionary forces that have shaped present-day genomes and allowed us to gain insight into how they are organised and regulated today.