DNA sequence evidence for the segmental allotetraploid origin of maize

Use of mRNA-seq to discriminate contributions to the transcriptome from the constituent genomes of the polyploid crop species Brassica napus

Background

Polyploidy often results in considerable changes in gene expression, both immediately and over evolutionary time. New phenotypes often arise with polyploid formation and may contribute to the fitness of polyploids in nature or their selection for use in agriculture. Oilseed rape (Brassica napus) is widely used to study the process of polyploidy both in artificially resynthesised and natural forms. mRNA-Seq, a recently developed approach to transcriptome profiling using deep-sequencing technologies is an alternative to microarrays for the study of gene expression in a polyploid.

Results

Illumina mRNA-Seq is comparable to microarray analysis for transcript quantification but has increased sensitivity and, very importantly, the potential to distinguish between homoeologous genes in polyploids. Using a novel curing process, we adapted a reference sequence that was a consensus derived from ESTs from both Brassica A and C genomes to one containing separate A and C genome versions for each of the 94,558 original unigenes. We aligned reads from B. napus to this cured reference, finding 38% more reads mapping from resynthesised lines and 28% more reads mapping from natural lines. Where the A and C versions differed at single nucleotide positions, termed inter-homoeologue polymorphisms (IHPs), we were able to apportion expression in the polyploid between the A and C genome homoeologues. 43,761 unigenes contained at least one IHP, with a mean frequency of 10.5 per kb unigene sequence. 6,350 of the unigenes with IHPs were differentially expressed between homoeologous gene pairs in resynthesised B. napus. 3,212 unigenes showed a similar pattern of differential expression across a range of natural B. napus crop varieties and, of these, 995 were in common with resynthesised B. napus. Functional classification showed over-representation in gene ontology categories not associated with dosage-sensitivity.

Conclusion

mRNA-Seq is the method of choice for measuring transcript abundance in polyploids due to its ability to measure the contributions of homoeologues to gene expression. The identification of large numbers of differentially expressed genes in both a newly resynthesised polyploid and natural B. napus confirms that there are both immediate and long-term alterations in the expression of homoeologous gene pairs following polyploidy.

Genetic Fingerprinting of Wheat and Its Progenitors by Mitochondrial Gene orf256

orf256 is a wheat mitochondrial gene associated with cytoplasmic male sterility (CMS) that has different organization in various species. This study exploited the orf256 gene as a mitochondrial DNA marker to study the genetic fingerprint of Triticum and Aegilops species. PCR followed by sequencing of common parts of the orf256 gene were employed to determine the fingerprint and molecular evolution of Triticum and Aegilops species. Although many primer pairs were used, two pairs of orf256 specific primers (5:-94/C: 482, 5:253/C: 482), amplified DNA fragments of 576 bp and 230 bp respectively in all species were tested. A common 500 bp of nine species of Triticum and Aegilops were aligned and showed consistent results with that obtained from other similar chloroplast or nuclear genes. Base alignment showed that there were various numbers of base substitutions in all species compared to S. cereal (Sc) (the outgroup species). Phylogenetic relationship revealed similar locations and proximity on phylogenetic trees established using plastid and nuclear genes. The results of this study open a good route to use unknown function genes of mitochondria in studying the molecular relationships and evolution of wheat and complex plant genomes.

Microsatellite allele dose and configuration establishment (MADCE): an integrated approach for genetic studies in allopolyploids

BACKGROUND

Genetic studies in allopolyploid plants are challenging because of the presence of similar sub-genomes, which leads to multiple alleles and complex segregation ratios. In this study, we describe a novel method for establishing the exact dose and configuration of microsatellite alleles for any accession of an allopolyploid plant species. The method, named Microsatellite Allele Dose and Configuration Establishment (MADCE), can be applied to mapping populations and pedigreed (breeding) germplasm in allopolyploids.

RESULTS

Two case studies are presented to demonstrate the power and robustness of the MADCE method. In the mapping case, five microsatellites were analysed. These microsatellites amplified 35 different alleles based on size. Using MADCE, we uncovered 30 highly informative segregating alleles. A conventional approach would have yielded only 19 fully informative and six partially informative alleles. Of the ten alleles that were present in all progeny (and thereby ignored or considered homozygous when using conventional approaches), six were found to segregate by dosage when analysed with MADCE. Moreover, the full allelic configuration of the mapping parents could be established, including null alleles, homozygous loci, and alleles that were present on multiple homoeologues. In the second case, 21 pedigreed cultivars were analysed using MADCE, resulting in the establishment of the full allelic configuration for all 21 cultivars and a tracing of allele flow over multiple generations.

CONCLUSIONS

The procedure described in this study (MADCE) enhances the efficiency and information content of mapping studies in allopolyploids. More importantly, it is the first technique to allow the determination of the full allelic configuration in pedigreed breeding germplasm from allopolyploid plants. This enables pedigree-based marker-trait association studies the use of algorithms developed for diploid crops, and it may increase the effectiveness of LD-based association studies. The MADCE method therefore enables researchers to tackle many of the genotyping problems that arise when performing mapping, pedigree, and association studies in allopolyploids. We discuss the merits of MADCE in comparison to other marker systems in polyploids, including SNPs, and how MADCE could aid in the development of SNP markers in allopolyploids.

Dynamic evolution of herbivore-induced sesquiterpene biosynthesis in sorghum and related grass crops

Sorghum (Sorghum bicolor) plants damaged by insects emit a blend of volatiles, predominantly sesquiterpenes, that are implicated in attracting natural enemies of the attacking insects. To characterize sesquiterpene biosynthesis in sorghum, seven terpene synthase (TPS) genes, SbTPS1 through SbTPS7, were identified based on their evolutionary relatedness to known sesquiterpene synthase genes from maize and rice. While SbTPS6 and SbTPS7 encode truncated proteins, all other TPS genes were determined to encode functional sesquiterpene synthases. Both SbTPS1 and SbTPS2 produced the major products zingiberene, β-bisabolene and β-sesquiphellandrene, but with opposite ratios of zingiberene to β-sesquiphellandrene. SbTPS3 produced (E)-α-bergamotene and (E)-β-farnesene. SbTPS4 formed (E)-β-caryophyllene as the major product. SbTPS5 produced mostly (E)-α-bergamotene and (Z)-γ-bisabolene. Based on the genome sequences of sorghum, maize and rice and the sesquiterpene synthase genes they contain, collinearity analysis identified the orthologs of sorghum sesquiterpene synthase genes, except for SbTPS4, in maize and rice. Phylogenetic analysis implied that SbTPS1, SbTPS2 and SbTPS3, which exist as tandem repeats, evolved as a consequence of local gene duplication in a lineage-specific manner. Structural modeling and site-directed mutagenesis experiments revealed that three amino acids in the active site play critical roles in defining product specificity of SbTPS1, SbTPS2, SbTPS3 and their orthologs in maize and rice. The naturally occurring functional variations of sesquiterpene synthases within and between species suggest that multiple mechanisms, including lineage-specific gene duplication, subfunctionalization, neofunctionalization and pseudogenization of duplicated genes, have all played a role in the dynamic evolution of insect-induced sesquiterpene biosynthesis in grasses.

The Maize PIN Gene Family of Auxin Transporters

Auxin is a key regulator of plant development and its differential distribution in plant tissues, established by a polar cell to cell transport, can trigger a wide range of developmental processes. A few members of the two families of auxin efflux transport proteins, PIN-formed (PIN) and P-glycoprotein (ABCB/PGP), have so far been characterized in maize. Nine new Zea mays auxin efflux carriers PIN family members and two maize PIN-like genes have now been identified. Four members of PIN1 (named ZmPIN1a-d) cluster, one gene homologous to AtPIN2 (ZmPIN2), three orthologs of PIN5 (ZmPIN5a-c), one gene paired with AtPIN8 (ZmPIN8), and three monocot-specific PINs (ZmPIN9, ZmPIN10a, and ZmPIN10b) were cloned and the phylogenetic relationships between early-land plants, monocots, and eudicots PIN proteins investigated, including the new maize PIN proteins. Tissue-specific expression patterns of the 12 maize PIN genes, 2 PIN-like genes and ZmABCB1, an ABCB auxin efflux carrier, were analyzed together with protein localization and auxin accumulation patterns in normal conditions and in response to drug applications. ZmPIN gene transcripts have overlapping expression domains in the root apex, during male and female inflorescence differentiation and kernel development. However, some PIN family members have specific tissue localization: ZmPIN1d transcript marks the L1 layer of the shoot apical meristem and inflorescence meristem during the flowering transition and the monocot-specific ZmPIN9 is expressed in the root endodermis and pericycle. The phylogenetic and gene structure analyses together with the expression pattern of the ZmPIN gene family indicate that subfunctionalization of some maize PINs can be associated to the differentiation and development of monocot-specific organs and tissues and might have occurred after the divergence between dicots and monocots.

Photorespiration and the evolution of C4 photosynthesis

C(4) photosynthesis is one of the most convergent evolutionary phenomena in the biological world, with at least 66 independent origins. Evidence from these lineages consistently indicates that the C(4) pathway is the end result of a series of evolutionary modifications to recover photorespired CO(2) in environments where RuBisCO oxygenation is high. Phylogenetically informed research indicates that the repositioning of mitochondria in the bundle sheath is one of the earliest steps in C(4) evolution, as it may establish a single-celled mechanism to scavenge photorespired CO(2) produced in the bundle sheath cells. Elaboration of this mechanism leads to the two-celled photorespiratory concentration mechanism known as C(2) photosynthesis (commonly observed in C(3)-C(4) intermediate species) and then to C(4) photosynthesis following the upregulation of a C(4) metabolic cycle.

Escape from preferential retention following repeated whole genome duplications in plants

The well supported gene dosage hypothesis predicts that genes encoding proteins engaged in dose-sensitive interactions cannot be reduced back to single copies once all interacting partners are simultaneously duplicated in a whole genome duplication. The genomes of extant flowering plants are the result of many sequential rounds of whole genome duplication, yet the fraction of genomes devoted to encoding complex molecular machines does not increase as fast as expected through multiple rounds of whole genome duplications. Using parallel interspecies genomic comparisons in the grasses and crucifers, we demonstrate that genes retained as duplicates following a whole genome duplication have only a 50% chance of being retained as duplicates in a second whole genome duplication. Genes which fractionated to a single copy following a second whole genome duplication tend to be the member of a gene pair with less complex promoters, lower levels of expression, and to be under lower levels of purifying selection. We suggest the copy with lower levels of expression and less purifying selection contributes less to effective gene-product dosage and therefore is under less dosage constraint in future whole genome duplications, providing an explanation for why flowering plant genomes are not overrun with subunits of large dose-sensitive protein complexes.

ZCN8 encodes a potential orthologue of Arabidopsis FT florigen that integrates both endogenous and photoperiod flowering signals in maize

Higher plants use multiple perceptive measures to coordinate flowering time with environmental and endogenous cues. Physiological studies show that florigen is a mobile factor that transmits floral inductive signals from the leaf to the shoot apex. Arabidopsis FT protein is widely regarded as the archetype florigen found in diverse plant species, particularly in plants that use inductive photoperiods to flower. Recently, a large family of FT homologues in maize, the Zea CENTRORADIALIS (ZCN) genes, was described, suggesting that maize also contains FT-related proteins that act as a florigen. The product of one member of this large family, ZCN8, has several attributes that make it a good candidate as a maize florigen. Mechanisms underlying the floral transition in maize are less well understood than those of other species, partly because flowering in temperate maize is dependent largely on endogenous signals. The maize indeterminate1 (id1) gene is an important regulator of maize autonomous flowering that acts in leaves to mediate the transmission or production of florigenic signals. This study finds that id1 acts upstream of ZCN8 to control its expression, suggesting a possible new link to flowering in day-neutral maize. Moreover, in teosinte, a tropical progenitor of maize that requires short-day photoperiods to induce flowering, ZCN8 is highly up-regulated in leaves under inductive photoperiods. Finally, vascular-specific expression of ZCN8 in Arabidopsis complements the ft-1 mutation, demonstrating that leaf-specific expression of ZCN8 can induce flowering. These results suggest that ZCN8 may encode a florigen that integrates both endogenous and environmental signals in maize.

Sequencing crop genomes: approaches and applications

Many challenges face plant scientists, in particular those working on crop production, such as a projected increase in population, decrease in water and arable land, changes in weather patterns and predictability. Advances in genome sequencing and resequencing can and should play a role in our response to meeting these challenges. However, several barriers prevent rapid and effective deployment of these tools to a wide variety of crops. Because of the complexity of crop genomes, de novo sequencing with next-generation sequencing technologies is a process fraught with difficulties that then create roadblocks to the utilization of these genome sequences for crop improvement. Collecting rapid and accurate phenotypes in crop plants is a hindrance to integrating genomics with crop improvement, and advances in informatics are needed to put these tools in the hands of the scientists on the ground.

Molecular aspects of flower development in grasses

The grass family (Poaceae) of the monocotyledons includes about 10,000 species and represents one of the most important taxa among angiosperms. Their flower morphology is remarkably different from those of other monocotyledons and higher eudicots. The peculiar floral structure of grasses is the floret, which contains carpels and stamens, like eudicots, but lacks petals and sepals. The reproductive organs are surrounded by two lodicules, which correspond to eudicot petals, and by a palea and lemma, whose correspondence to eudicot organs remains controversial. The molecular and genetic analysis of floral morphogenesis and organ specification, primarily performed in eudicot model species, led to the ABCDE model of flower development. Several genes required for floral development in grasses correspond to class A, B, C, D, and E genes of eudicots, but others appear to have unique and diversified functions. In this paper, we outline the present knowledge on the evolution and diversification of grass genes encoding MIKC-type MADS-box transcription factors, based on information derived from studies in rice, maize, and wheat. Moreover, we review recent advances in studying the genes involved in the control of flower development and the extent of structural and functional conservation of these genes between grasses and eudicots.

Evolutionary divergence and limits of conserved non-coding sequence detection in plant genomes

The discovery of regulatory motifs embedded in upstream regions of plants is a particularly challenging bioinformatics task. Previous studies have shown that motifs in plants are short compared with those found in vertebrates. Furthermore, plant genomes have undergone several diversification mechanisms such as genome duplication events which impact the evolution of regulatory motifs. In this article, a systematic phylogenomic comparison of upstream regions is conducted to further identify features of the plant regulatory genomes, the component of genomes regulating gene expression, to enable future de novo discoveries. The findings highlight differences in upstream region properties between major plant groups and the effects of divergence times and duplication events. First, clear differences in upstream region evolution can be detected between monocots and dicots, thus suggesting that a separation of these groups should be made when searching for novel regulatory motifs, particularly since universal motifs such as the TATA box are rare. Second, investigating the decay rate of significantly aligned regions suggests that a divergence time of ~100 mya sets a limit for reliable conserved non-coding sequence (CNS) detection. Insights presented here will set a framework to help identify embedded motifs of functional relevance by understanding the limits of bioinformatics detection for CNSs.

Loss of centromeric histone H3 (CENH3) from centromeres precedes uniparental chromosome elimination in interspecific barley hybrids

Uniparental chromosome elimination occurs in several interspecific hybrids of plants. We studied the mechanism underlying selective elimination of the paternal chromosomes during the early development of Hordeum vulgare × Hordeum bulbosum embryos. The following conclusions regarding the role of the centromere-specific histone H3 variant (CENH3) in the process of chromosome elimination were drawn: (i) centromere inactivity of H. bulbosum chromosomes triggers the mitosis-dependent process of uniparental chromosome elimination in unstable H. vulgare × H. bulbosum hybrids; (ii) centromeric loss of CENH3 protein rather than uniparental silencing of CENH3 genes causes centromere inactivity; (iii) in stable species combinations, cross-species incorporation of CENH3 occurs despite centromere-sequence differences, and not all CENH3 variants get incorporated into centromeres if multiple CENH3s are present in species combinations; and (iv) diploid barley species encode two CENH3 variants, the proteins of which are intermingled within centromeres throughout mitosis and meiosis.

The use of 2D-electrophoresis and de novo sequencing to characterize inter- and intra-cultivar protein polymorphisms in an allopolyploid crop

Polyploidy and allopolyploidy have played an important role in the evolution of many plants and crops. Several techniques exist to characterize allopolyploid varieties. Analyzing the consequences of genomic reorganization at the gDNA level is a prerequisite but a better insight into the consequences for the phenotype is also primordial. As such, protein polymorphism analysis is important in understanding plant and crop biodiversity and is a driving force behind crop improvement. Our strategy to analyze protein isoforms and to detect possible gene silencing or deletion in bananas was based on protein analysis. Bananas are a good representative of a complex allopolyploid and important crop. We combined two-dimensional electrophoresis (2DE) and 2D DIGE with de novo MS/MS sequence determination to characterize a range of triploid varieties. Via Principal Component Analysis (PCA) and hierarchical clustering we were able to blindly classify the different varieties according to their presumed genome constitution. We report for the first time the application of an automated approach for the derivatization of peptides for facilitated MS/MS de novo sequence determination. We conclude that the proteome does not always correspond to the presumed genome formulae and that proteomics is a powerful tool to characterize varieties. The observations at the protein level provide good indications for a more complex genome structure and genomic rearrangement in some banana varieties.

High time for a roll call: gene duplication and phylogenetic relationships of TCP-like genes in monocots

BACKGROUND AND AIMS

The TCP family is an ancient group of plant developmental transcription factors that regulate cell division in vegetative and reproductive structures and are essential in the establishment of flower zygomorphy. In-depth research on eudicot TCPs has documented their evolutionary and developmental role. This has not happened to the same extent in monocots, although zygomorphy has been critical for the diversification of Orchidaceae and Poaceae, the largest families of this group. Investigating the evolution and function of TCP-like genes in a wider group of monocots requires a detailed phylogenetic analysis of all available sequence information and a system that facilitates comparing genetic and functional information.

METHODS

The phylogenetic relationships of TCP-like genes in monocots were investigated by analysing sequences from the genomes of Zea mays, Brachypodium distachyon, Oryza sativa and Sorghum bicolor, as well as EST data from several other monocot species.

KEY RESULTS

All available monocot TCP-like sequences are associated in 20 major groups with an average identity ≥64 % and most correspond to well-supported clades of the phylogeny. Their sequence motifs and relationships of orthology were documented and it was found that 67 % of the TCP-like genes of Sorghum, Oryza, Zea and Brachypodium are in microsyntenic regions. This analysis suggests that two rounds of whole genome duplication drove the expansion of TCP-like genes in these species.

CONCLUSIONS

A system of classification is proposed where putative or recognized monocot TCP-like genes are assigned to a specific clade of PCF-, CIN- or CYC/tb1-like genes. Specific biases in sequence data of this family that must be tackled when studying its molecular evolution and phylogeny are documented. Finally, the significant retention of duplicated TCP genes from Zea mays is considered in the context of balanced gene drive.

BARREN STALK FASTIGIATE1 is an AT-hook protein required for the formation of maize ears

Ears are the seed-bearing inflorescences of maize (Zea mays) plants and represent a crucial component of maize yield. The first step in the formation of ears is the initiation of axillary meristems in the axils of developing leaves. In the classic maize mutant barren stalk fastigiate1 (baf1), first discovered in the 1950s, ears either do not form or, if they do, are partially fused to the main stalk. We positionally cloned Baf1 and found that it encodes a transcriptional regulator containing an AT-hook DNA binding motif. Single coorthologs of Baf1 are found in syntenic regions of brachypodium (Brachypodium distachyon), rice (Oryza sativa), and sorghum (Sorghum bicolor), suggesting that the gene is likely present in all cereal species. Protein-protein interaction assays suggest that BAF1 is capable of forming homodimers and heterodimers with other members of the AT-hook family. Another transcriptional regulator required for ear initiation is the basic helix-loop-helix protein BARREN STALK1 (BA1). Genetic and expression analyses suggest that Baf1 is required to reach a threshold level of Ba1 expression for the initiation of maize ears. We propose that Baf1 functions in the demarcation of a boundary region essential for the specification of a stem cell niche.

Gene space dynamics during the evolution of Aegilops tauschii, Brachypodium distachyon, Oryza sativa, and Sorghum bicolor genomes

Nine different regions totaling 9.7 Mb of the 4.02 Gb Aegilops tauschii genome were sequenced using the Sanger sequencing technology and compared with orthologous Brachypodium distachyon, Oryza sativa (rice), and Sorghum bicolor (sorghum) genomic sequences. The ancestral gene content in these regions was inferred and used to estimate gene deletion and gene duplication rates along each branch of the phylogenetic tree relating the four species. The total gene number in the extant Ae. tauschii genome was estimated to be 36,371. The gene deletion and gene duplication rates and total gene numbers in the four genomes were used to estimate the total gene number in each node of the phylogenetic tree. The common ancestor of the Brachypodieae and Triticeae lineages was estimated to have had 28,558 genes, and the common ancestor of the Panicoideae, Ehrhartoideae, and Pooideae subfamilies was estimated to have had 27,152 or 28,350 genes, depending on the ancestral gene scenario. Relative to the Brachypodieae and Triticeae common ancestor, the gene number was reduced in B. distachyon by 3,026 genes and increased in Ae. tauschii by 7,813 genes. The sum of gene deletion and gene duplication rates, which reflects the rate of gene synteny loss, was correlated with the rate of structural chromosome rearrangements and was highest in the Ae. tauschii lineage and lowest in the rice lineage. The high rate of gene space evolution in the Ae. tauschii lineage accounts for the fact that, contrary to the expectations, the level of synteny between the phylogenetically more related Ae. tauschii and B. distachyon genomes is similar to the level of synteny between the Ae. tauschii genome and the genomes of the less related rice and sorghum. The ratio of gene duplication to gene deletion rates in these four grass species closely parallels both the total number of genes in a species and the overall genome size. Because the overall genome size is to a large extent a function of the repeated sequence content in a genome, we suggest that the amount and activity of repeated sequences are important factors determining the number of genes in a genome.

Expression of CENH3 alleles in synthesized allopolyploid Oryza species

Synthesized allopolyploids are valuable materials for comparative analyses of two or more distinct genomes, such as the expression changes (activation, inactivation or differential expression) of orthologous genes following allopolyploidization. CENH3 is a centromere- specific histone H3 variant and has been regarded as a central component in kinetochore formation and centromere function. In this study, interspecific hybrids of Oryza genus (AA × CC, AA × CCDD) and their backcross progenies were produced, and the genome constitutions were identified as AC, ACC, ACD, AACD, or AA(CD) by Genomic in situ hybridization (GISH). We further cloned and sequenced the CENH3 genes from O. sativa (AA), O. officinalis (CC) and O. latifolia (CCDD). Sequencing of RT-PCR products revealed that CENH3_C2 and CENH3_D, the two CENH3 alleles from O. latifolia, showed polymorphism in several sites, while CENH3_C2 and CENH3_C1 from O. officinalis were different at only two amino acids positions. Moreover, we found that the CENH3 genes from both parents are expressed in interspecific hybrids and their progenies. Specifically, based on our cDNA sequencing data, the ratio of expression level between CENH3_A and CENH3_C1 was approximately 1 in AC and 0.5 in ACC genomes, respectively. As a result, the CENH3 expression patterns shed more light on the inter-coordination between varied centromeric DNA sequences and highly conserved kinetochore protein in synthesized allopolyploids of Oryza genus.

Changes in expression pattern of the teosinte branched1-like genes in the Zingiberales provide a mechanism for evolutionary shifts in symmetry across the order

PREMISE OF THE STUDY

Floral symmetry is a trait of key importance when considering floral diversification because it is thought to play a significant role in plant-pollinator interactions. The CYCLOIDEA/TEOSINTE BRANCHED1 (CYC/TB1)-like genes have been implicated in the development and evolution of floral symmetry in numerous lineages. We thus chose to investigate a possible role for these genes in the evolution of floral symmetry within petaloid monocots, using the order Zingiberales as a model system. In the Zingiberales, evolutionary shifts in symmetry have occurred in all floral whorls, making the order ideal for studying the evolution of this ecologically significant trait.

METHODS

We analyzed TB1-like (TBL) genes from taxa spanning the order in a phylogenetic context. Using RNA in situ hybridization, we examined the expression of two TBL genes in Costus spicatus (Costaceae) and Heliconia stricta (Heliconiaceae), taxa with divergent floral symmetry patterns.

KEY RESULTS

We identified Zingiberales-specific gene duplications as well as a duplication in the TBL gene lineage that predates the diversification of commelinid monocots. Shifts in TBL gene expression were associated with evolutionary shifts in floral symmetry and stamen abortion. ZinTBL1a expression was found in the posterior (adaxial) staminode of H. stricta and in the abaxial staminodial labellum of C. spicatus. ZinTBL2 expression was strongest in the anterior (abaxial) sepals of H. stricta and in the adaxial fertile stamen of C. spicatus.

CONCLUSIONS

This study adds to the growing body of evidence that CYC/TB1-like genes have been repeatedly recruited throughout the course of evolution to generate bilateral floral symmetry (zygomorphy).

Phenotypic and gene expression analyses of a ploidy series of maize inbred Oh43

Polyploidization has repeatedly occurred during plant evolution. Although autopolyploidy is the best model to characterize the polyploidization effects in a highly controlled manner, there are limited studies on autopolyploids compared to allopolyploids. To improve our understanding of autopolyploidy effects in maize, we developed an inbred Oh43 ploidy series consisting of the diploid (2X), tetraploid (4X) and hexaploid (6X) lines and compared their phenotypes and gene expression in the mature adult leaf tissue. Our phenotypic study showed that plants of higher ploidy exhibit increased cell size but slower growth rate, later flowering, fewer tassel branches, reduced stature and fertility. Two-dimensional difference gel electrophoresis (2D DIGE) and gel electrophoresis followed by liquid chromatography and mass spectrometry (GeLC-MS) assays of the leaf proteomes revealed ~40 and 26% quantitative differentially expressed (DE) proteins, respectively, at the per genome level. A small number of qualitative DE proteins were also identified in the GeLC-MS assay. The majority of the quantitative DE proteins found in the 2D DIGE assay were present in either the 4X versus 6X or the 2X versus 6X comparison but not the 2X versus 4X comparison. Aneuploidy in some 6X plants might contribute to the more extensive changes of gene expression per genome in the 6X. Most changes of the protein expression per genome are less than twofold. Less than 5% of the DE genes exhibit a positive or negative continuous correlation through the ploidy series between their protein expression per genome, and the genome copy number. Hence, in the Oh43 ploidy series, expression for most proteins in a cell increases linearly with ploidy.

Back to the sea twice: identifying candidate plant genes for molecular evolution to marine life

Background

Seagrasses are a polyphyletic group of monocotyledonous angiosperms that have adapted to a completely submerged lifestyle in marine waters. Here, we exploit two collections of expressed sequence tags (ESTs) of two wide-spread and ecologically important seagrass species, the Mediterranean seagrass Posidonia oceanica (L.) Delile and the eelgrass Zostera marina L., which have independently evolved from aquatic ancestors. This replicated, yet independent evolutionary history facilitates the identification of traits that may have evolved in parallel and are possible instrumental candidates for adaptation to a marine habitat.

Results

In our study, we provide the first quantitative perspective on molecular adaptations in two seagrass species. By constructing orthologous gene clusters shared between two seagrasses (Z. marina and P. oceanica) and eight distantly related terrestrial angiosperm species, 51 genes could be identified with detection of positive selection along the seagrass branches of the phylogenetic tree. Characterization of these positively selected genes using KEGG pathways and the Gene Ontology uncovered that these genes are mostly involved in translation, metabolism, and photosynthesis.

Conclusions

These results provide first insights into which seagrass genes have diverged from their terrestrial counterparts via an initial aquatic stage characteristic of the order and to the derived fully-marine stage characteristic of seagrasses. We discuss how adaptive changes in these processes may have contributed to the evolution towards an aquatic and marine existence.

Karyotype of Zea luxurians and Z. mays subsp. mays using FISH/DAPI, and analysis of meiotic behavior of hybrids

The karyotypes of Zea luxurians and a race of maize from northwestern Argentina are described and compared using 4′,6-diamidino-2-phenylindole (DAPI) banding and fluorescent in situ hybridization (FISH) to localize the 180 bp knobs. The meiotic behavior of the F₁ artificial hybrids Z. luxurians × maize is also analyzed to determine the genomic relationships between both species. Neocentromere activity at knobs in the meiosis of the hybrids is particularly discussed. The meiotic behavior and the high pollen sterility of the hybrid revealed genetical and (or) chromosomal divergences, leading to postzygotic reproductive isolation among their parents. Here, we propose that maize shows lower genomic affinity to Z. luxurians than to other species of the genus with 2n = 20.

Dose-sensitivity, conserved non-coding sequences, and duplicate gene retention through multiple tetraploidies in the grasses

Whole genome duplications, or tetraploidies, are an important source of increased gene content. Following whole genome duplication, duplicate copies of many genes are lost from the genome. This loss of genes is biased both in the classes of genes deleted and the subgenome from which they are lost. Many or all classes are genes preferentially retained as duplicate copies are engaged in dose sensitive protein-protein interactions, such that deletion of any one duplicate upsets the status quo of subunit concentrations, and presumably lowers fitness as a result. Transcription factors are also preferentially retained following every whole genome duplications studied. This has been explained as a consequence of protein-protein interactions, just as for other highly retained classes of genes. We show that the quantity of conserved noncoding sequences (CNSs) associated with genes predicts the likelihood of their retention as duplicate pairs following whole genome duplication. As many CNSs likely represent binding sites for transcriptional regulators, we propose that the likelihood of gene retention following tetraploidy may also be influenced by dose-sensitive protein-DNA interactions between the regulatory regions of CNS-rich genes - nicknamed bigfoot genes - and the proteins that bind to them. Using grass genomes, we show that differential loss of CNSs from one member of a pair following the pre-grass tetraploidy reduces its chance of retention in the subsequent maize lineage tetraploidy.

DNA Methylation in Genomes of Several Annual Herbaceous and Woody Perennial Plants of Varying Ploidy as Detected by MSAP

Polyploidization is known to accompany altered DNA methylation in higher plants, which plays an important role in gene expression regulation and maintaining genome stability. While the characteristics of DNA methylation in different polyploid plants are still to be elucidated; here, status of genomic DNA methylation in a series of diploid, triploid, and tetraploid annual herbaceous plants (watermelon and Salvia) and woody perennials (pear, Poplar, and loquat) were explored by methylation-specific amplified polymorphism analysis. The results indicated that levels of DNA methylation in triploid watermelon and Salvia were lower than their diploid parents. In triploid Poplar and pear, higher levels of DNA methylation were detected, and no significant difference was observed between triploid and tetraploid in all tested materials. Further data analysis suggested that about half of the total detected sites underwent changes of DNA methylation patterns in triploid watermelons and Salvia, as well as an obvious trend towards demethylation. However, the changes of DNA methylation patterns in three triploid woody perennials were only 17.54-33.40%. This implied that the characteristics of DNA methylation are significantly different during the polyploidization of different plant species. Furthermore, the results suggested that the level of DNA methylation was nonlinearly related to the ploidy level, and triploid plants displayed more interesting DNA methylation status. The characteristics and possible functions of DNA methylation in different ploidy series are further discussed.

The amplification and evolution of orthologous 22-kDa α-prolamin tandemly arrayed genes in coix, sorghum and maize genomes

Tandemly arrayed genes (TAGs) account for about one-third of the duplicated genes in eukaryotic genomes. They provide raw genetic material for biological evolution, and play important roles in genome evolution. The 22-kDa prolamin genes in cereal genomes represent typical TAG organization, and provide the good material to investigate gene amplification of TAGs in closely related grass genomes. Here, we isolated and sequenced the Coix 22-kDa prolamin (coixin) gene cluster (283 kb), and carried out a comparative analysis with orthologous 22-kDa prolamin gene clusters from maize and sorghum. The 22-kDa prolamin gene clusters descended from orthologous ancestor genes, but underwent independent gene amplification paths after the separation of these species, therefore varied dramatically in sequence and organization. Our analysis indicated that the gene amplification model of 22-kDa prolamin gene clusters can be divided into three major stages. In the first stage, rare gene duplications occurred from the ancestor gene copy accidentally. In the second stage, rounds of gene amplification occurred by unequal crossing over to form tandem gene array(s). In the third stage, gene array was further diverged by other genomic activities, such as transposon insertions, segmental rearrangements, etc. Unlike their highly conserved sequences, the amplified 22-kDa prolamin genes diverged rapidly at their expression capacities and expression levels. Such processes had no apparent correlation to age or order of amplified genes within TAG cluster, suggesting a fast evolving nature of TAGs after gene amplification. These results provided insights into the amplification and evolution of TAG families in grasses.

Pervasive gene content variation and copy number variation in maize and its undomesticated progenitor

Individuals of the same species are generally thought to have very similar genomes. However, there is growing evidence that structural variation in the form of copy number variation (CNV) and presence-absence variation (PAV) can lead to variation in the genome content of individuals within a species. Array comparative genomic hybridization (CGH) was used to compare gene content and copy number variation among 19 diverse maize inbreds and 14 genotypes of the wild ancestor of maize, teosinte. We identified 479 genes exhibiting higher copy number in some genotypes (UpCNV) and 3410 genes that have either fewer copies or are missing in the genome of at least one genotype relative to B73 (DownCNV/PAV). Many of these DownCNV/PAV are examples of genes present in B73, but missing from other genotypes. Over 70% of the CNV/PAV examples are identified in multiple genotypes, and the majority of events are observed in both maize and teosinte, suggesting that these variants predate domestication and that there is not strong selection acting against them. Many of the genes affected by CNV/PAV are either maize specific (thus possible annotation artifacts) or members of large gene families, suggesting that the gene loss can be tolerated through buffering by redundant functions encoded elsewhere in the genome. While this structural variation may not result in major qualitative variation due to genetic buffering, it may significantly contribute to quantitative variation.

Impact of natural genetic variation on the transcriptome of autotetraploid Arabidopsis thaliana

Polyploidy, the presence of more than two complete sets of chromosomes in an organism, has significantly shaped the genomes of angiosperms during evolution. Two forms of polyploidy are often considered: allopolyploidy, which originates from interspecies hybrids, and autopolyploidy, which originates from intraspecies genome duplication events. Besides affecting genome organization, polyploidy generates other genetic effects. Synthetic allopolyploid plants exhibit considerable transcriptome alterations, part of which are likely caused by the reunion of previously diverged regulatory hierarchies. In contrast, autopolyploids have relatively uniform genomes, suggesting lower alteration of gene expression. To evaluate the impact of intraspecies genome duplication on the transcriptome, we generated a series of unique Arabidopsis thaliana autotetraploids by using different ecotypes. A. thaliana autotetraploids show transcriptome alterations that strongly depend on their parental genome composition and include changed expression of both new genes and gene groups previously described from allopolyploid Arabidopsis. Alterations in gene expression are stable, nonstochastic, developmentally specific, and associated with changes in DNA methylation. We propose that Arabidopsis possesses an inherent and heritable ability to sense and respond to elevated, yet balanced chromosome numbers. The impact of natural variation on alteration of autotetraploid gene expression stresses its potential importance in the evolution and breeding of plants.

Ancient genome duplications during the evolution of kiwifruit (Actinidia) and related Ericales

BACKGROUND AND AIMS

To assess the number and phylogenetic distribution of large-scale genome duplications in the ancestry of Actinidia, publicly available expressed sequenced tags (ESTs) for members of the Actinidiaceae and related Ericales, including tea (Camellia sinensis), were analysed.

METHODS

Synonymous divergences (K(s)) were calculated for all duplications within gene families and examined for evidence of large-scale duplication events. Phylogenetic comparisons for a selection of orthologues among several related species in Ericales and two outgroups permitted placement of duplication events in relation to lineage divergences. Gene ontology (GO) categories were analysed for each whole-genome duplication (WGD) and the whole transcriptome.

KEY RESULTS

Evidence for three ancient WGDs in Actinidia was found. Analyses of paleologue GO categories indicated a different pattern of retained genes for each genome duplication, but a pattern consistent with the dosage-balance hypothesis among all retained paleologues.

CONCLUSIONS

This study provides evidence for one independent WGD in the ancestry of Actinidia (Ad-alpha), a WGD shared by Actinidia and Camellia (Ad-beta), and the well-established At-gamma WGD that occurred prior to the divergence of all taxa examined. More ESTs in other taxa are needed to elucidate which groups in Ericales share the Ad-beta or Ad-alpha duplications and their impact on diversification.

Development of highly polymorphic SNP markers from the complexity reduced portion of maize [Zea mays L.] genome for use in marker-assisted breeding

The duplicated and the highly repetitive nature of the maize genome has historically impeded the development of true single nucleotide polymorphism (SNP) markers in this crop. Recent advances in genome complexity reduction methods coupled with sequencing-by-synthesis technologies permit the implementation of efficient genome-wide SNP discovery in maize. In this study, we have applied Complexity Reduction of Polymorphic Sequences technology (Keygene N.V., Wageningen, The Netherlands) for the identification of informative SNPs between two genetically distinct maize inbred lines of North and South American origins. This approach resulted in the discovery of 1,123 putative SNPs representing low and single copy loci. In silico and experimental (Illumina GoldenGate (GG) assay) validation of putative SNPs resulted in mapping of 604 markers, out of which 188 SNPs represented 43 haplotype blocks distributed across all ten chromosomes. We have determined and clearly stated a specific combination of stringent criteria (>0.3 minor allele frequency, >0.8 GenTrainScore and >0.5 Chi_test100 score) necessary for the identification of highly polymorphic and genetically stable SNP markers. Due to these criteria, we identified a subset of 120 high-quality SNP markers to leverage in GG assay-based marker-assisted selection projects. A total of 32 high-quality SNPs represented 21 haplotypes out of 43 identified in this study. The information on the selection criteria of highly polymorphic SNPs in a complex genome such as maize and the public availability of these SNP assays will be of great value for the maize molecular genetics and breeding community.

Identification and characterization of the maize arogenate dehydrogenase gene family

In plants, the amino acids tyrosine and phenylalanine are synthesized from arogenate by arogenate dehydrogenase and arogenate dehydratase, respectively, with the relative flux to each being tightly controlled. Here the characterization of a maize opaque endosperm mutant (mto140), which also shows retarded vegetative growth, is described The opaque phenotype co-segregates with a Mutator transposon insertion in an arogenate dehydrogenase gene (zmAroDH-1) and this led to the characterization of the four-member family of maize arogenate dehydrogenase genes (zmAroDH-1-zmAroDH-4) which share highly similar sequences. A Mutator insertion at an equivalent position in AroDH-3, the most closely related family member to AroDH-1, is also associated with opaque endosperm and stunted vegetative growth phenotypes. Overlapping but differential expression patterns as well as subtle mutant effects on the accumulation of tyrosine and phenylalanine in endosperm, embryo, and leaf tissues suggest that the functional redundancy of this gene family provides metabolic plasticity for the synthesis of these important amino acids. mto140/arodh-1 seeds shows a general reduction in zein storage protein accumulation and an elevated lysine phenotype typical of other opaque endosperm mutants, but it is distinct because it does not result from quantitative or qualitative defects in the accumulation of specific zeins but rather from a disruption in amino acid biosynthesis.

Variation in antiherbivore defense responses in synthetic Nicotiana allopolyploids correlates with changes in uniparental patterns of gene expression

We examined the expression of Nicotiana attenuata (Na) and Nicotiana obtusifolia (No) herbivore-induced genes in synthetic autopolyploids (NaT and NoT) and five independent allopolyploid Nicotiana x obtusiata (Nxo) lines to understand how the expression of genes regulating complex polygenetic defense traits is altered in the early stages of allopolyploid hybridization. In Na, applying Manduca sexta oral secretions (OS) to wounds rapidly increased the transcript accumulation of wound-induced protein kinase (WIPK), lipoxygenase 3 (LOX3), nonexpressor of pathogenesis-related 1 (NPR1), and jasmonate-resistant 4 (JAR4) genes; these were correlated with increases in accumulation of jasmonic acid (JA), jasmonate-isoleucine, and trypsin protease inhibitors (TPIs). In No, OS elicitation reduced NPR1 transcripts and increased the level of salicylic acid (SA) that appeared to antagonize JA and JA-mediated defenses. OS elicited Nxo lines, accumulated high levels of the uniparental transcript of WIPK, LOX3, JAR4, and TPI, but low levels of both parental NPR1 transcripts that in turn were correlated with an increase in SA and a decrease in JA levels, suggesting SA/JA antagonism in the allopolyploid crosses. Methyl jasmonate treatment of Nxo lines elicited transcripts of both parental LOX3, JAR4, and TPIs, demonstrating that the uniparental pattern observed after OS elicitation was not due to gene inactivation. TPIs were induced at different levels among Nxo lines; some lines expressed high levels comparable to Na, others low levels similar to No, suggesting that synthetic neoallopolyploids rapidly readjust the expression of their parental defensive genes to generate diverse antiherbivore responses. Changes in the expression of key genes and posttranscriptional events likely facilitate adaptive radiations during allopolyploid speciation events.

Relationship, evolutionary fate and function of two maize co-orthologs of rice GW2 associated with kernel size and weight

BACKGROUND

In rice, the GW2 gene, found on chromosome 2, controls grain width and weight. Two homologs of this gene, ZmGW2-CHR4 and ZmGW2-CHR5, have been found in maize. In this study, we investigated the relationship, evolutionary fate and putative function of these two maize genes.

RESULTS

The two genes are located on duplicated maize chromosomal regions that show co-orthologous relationships with the rice region containing GW2. ZmGW2-CHR5 is more closely related to the sorghum counterpart than to ZmGW2-CHR4. Sequence comparisons between the two genes in eight diverse maize inbred lines revealed that the functional protein domain of both genes is completely conserved, with no non-synonymous polymorphisms identified. This suggests that both genes may have conserved functions, a hypothesis that was further confirmed through linkage, association, and expression analyses. Linkage analysis showed that ZmGW2-CHR4 is located within a consistent quantitative trait locus (QTL) for one-hundred kernel weight (HKW). Association analysis with a diverse panel of 121 maize inbred lines identified one single nucleotide polymorphism (SNP) in the promoter region of ZmGW2-CHR4 that was significantly associated with kernel width (KW) and HKW across all three field experiments examined in this study. SNPs or insertion/deletion polymorphisms (InDels) in other regions of ZmGW2-CHR4 and ZmGW2-CHR5 were also found to be significantly associated with at least one of the four yield-related traits (kernel length (KL), kernel thickness (KT), KW and HKW). None of the polymorphisms in either maize gene are similar to each other or to the 1 bp InDel causing phenotypic variation in rice. Expression levels of both maize genes vary over ear and kernel developmental stages, and the expression level of ZmGW2-CHR4 is significantly negatively correlated with KW.

CONCLUSIONS

The sequence, linkage, association and expression analyses collectively showed that the two maize genes represent chromosomal duplicates, both of which function to control some of the phenotypic variation for kernel size and weight in maize, as does their counterpart in rice. However, the different polymorphisms identified in the two maize genes and in the rice gene indicate that they may cause phenotypic variation through different mechanisms.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

Evidence for the involvement of Globosa-like gene duplications and expression divergence in the evolution of floral morphology in the Zingiberales

*The MADS box transcription factor family has long been identified as an important contributor to the control of floral development. It is often hypothesized that the evolution of floral development across angiosperms and within specific lineages may occur as a result of duplication, functional diversification, and changes in regulation of MADS box genes. Here we examine the role of Globosa (GLO)-like genes, members of the B-class MADS box gene lineage, in the evolution of floral development within the monocot order Zingiberales. *We assessed changes in perianth and stamen whorl morphology in a phylogenetic framework. We identified GLO homologs (ZinGLO1-4) from 50 Zingiberales species and investigated the evolution of this gene lineage. Expression of two GLO homologs was assessed in Costus spicatus and Musa basjoo. *Based on the phylogenetic data and expression results, we propose several family-specific losses and gains of GLO homologs that appear to be associated with key morphological changes. The GLO-like gene lineage has diversified concomitant with the evolution of the dimorphic perianth and the staminodial labellum. *Duplications and expression divergence within the GLO-like gene lineage may have played a role in floral diversification in the Zingiberales.

Following tetraploidy in maize, a short deletion mechanism removed genes preferentially from one of the two homologs

Following genome duplication and selfish DNA expansion, maize used a heretofore unknown mechanism to shed redundant genes and functionless DNA with bias toward one of the parental genomes.

Previous work in Arabidopsis showed that after an ancient tetraploidy event, genes were preferentially removed from one of the two homeologs, a process known as fractionation. The mechanism of fractionation is unknown. We sought to determine whether such preferential, or biased, fractionation exists in maize and, if so, whether a specific mechanism could be implicated in this process. We studied the process of fractionation using two recently sequenced grass species: sorghum and maize. The maize lineage has experienced a tetraploidy since its divergence from sorghum approximately 12 million years ago, and fragments of many knocked-out genes retain enough sequence similarity to be easily identifiable. Using sorghum exons as the query sequence, we studied the fate of both orthologous genes in maize following the maize tetraploidy. We show that genes are predominantly lost, not relocated, and that single-gene loss by deletion is the rule. Based on comparisons with orthologous sorghum and rice genes, we also infer that the sequences present before the deletion events were flanked by short direct repeats, a signature of intra-chromosomal recombination. Evidence of this deletion mechanism is found 2.3 times more frequently on one of the maize homeologs, consistent with earlier observations of biased fractionation. The over-fractionated homeolog is also a greater than 3-fold better target for transposon removal, but does not have an observably higher synonymous base substitution rate, nor could we find differentially placed methylation domains. We conclude that fractionation is indeed biased in maize and that intra-chromosomal or possibly a similar illegitimate recombination is the primary mechanism by which fractionation occurs. The mechanism of intra-chromosomal recombination explains the observed bias in both gene and transposon loss in the maize lineage. The existence of fractionation bias demonstrates that the frequency of deletion is modulated. Among the evolutionary benefits of this deletion/fractionation mechanism is bulk DNA removal and the generation of novel combinations of regulatory sequences and coding regions.

All genomes can accumulate dispensable DNA in the form of duplications of individual genes or even partial or whole genome duplications. Genomes also can accumulate selfish DNA elements. Duplication events specifically are often followed by extensive gene loss. The maize genome is particularly extreme, having become tetraploid 10 million years ago and played host to massive transposon amplifications. We compared the genome of sorghum (which is homologous to the pre-tetraploid maize genome) with the two identifiable parental genomes retained in maize. The two maize genomes differ greatly: one of the parental genomes has lost 2.3 times more genes than the other, and the selfish DNA regions between genes were even more frequently lost, suggesting maize can distinguish between the parental genomes present in the original tetraploid. We show that genes are actually lost, not simply relocated. Deletions were rarely longer than a single gene, and occurred between repeated DNA sequences, suggesting mis-recombination as a mechanism of gene removal. We hypothesize an epigenetic mechanism of genome distinction to account for the selective loss. To the extent that the rate of base substitutions tracks time, we neither support nor refute claims of maize allotetraploidy. Finally, we explain why it makes sense that purifying selection in mammals does not operate at all like the gene and genome deletion program we describe here.

On the tetraploid origin of the maize genome

Data from cytological and genetic mapping studies suggest that maize arose as a tetraploid. Two previous studies investigating the most likely mode of maize origin arrived at different conclusions. Gaut and Doebley proposed a segmental allotetraploid origin of the maize genome and estimated that the two maize progenitors diverged at 20.5 million years ago (mya). In a similar study, using larger data set, Brendel and colleagues suggested a single genome duplication at 16 mya. One of the key components of such analyses is to examine sequence divergence among strictly orthologous genes. In order to identify such genes, Lai and colleagues sequenced five duplicated chromosomal regions from the maize genome and the orthologous counterparts from the sorghum genome. They also identified the orthologous regions in rice. Using positional information of genetic components, they identified 11 orthologous genes across the two duplicated regions of maize, and the sorghum and rice regions. Swigonova et al. analyzed the 11 orthologues, and showed that all five maize chromosomal regions duplicated at the same time, supporting a tetraploid origin of maize, and that the two maize progenitors diverged from each other at about the same time as each of them diverged from sorghum, about 11.9 mya.

Genome-wide distribution of transposed Dissociation elements in maize

The maize (Zea mays) transposable element Dissociation (Ds) was mobilized for large-scale genome mutagenesis and to study its endogenous biology. Starting from a single donor locus on chromosome 10, over 1500 elements were distributed throughout the genome and positioned on the maize physical map. Genetic strategies to enrich for both local and unlinked insertions were used to distribute Ds insertions. Global, regional, and local insertion site trends were examined. We show that Ds transposed to both linked and unlinked sites and displayed a nonuniform distribution on the genetic map around the donor r1-sc:m3 locus. Comparison of Ds and Mutator insertions reveals distinct target preferences, which provide functional complementarity of the two elements for gene tagging in maize. In particular, Ds displays a stronger preference for insertions within exons and introns, whereas Mutator insertions are more enriched in promoters and 5'-untranslated regions. Ds has no strong target site consensus sequence, but we identified properties of the DNA molecule inherent to its local structure that may influence Ds target site selection. We discuss the utility of Ds for forward and reverse genetics in maize and provide evidence that genes within a 2- to 3-centimorgan region flanking Ds insertions will serve as optimal targets for regional mutagenesis.

The psychedelic genes of maize redundantly promote carbohydrate export from leaves

Whole-plant carbohydrate partitioning involves the assimilation of carbon in leaves and its translocation to nonphotosynthetic tissues. This process is fundamental to plant growth and development, but its regulation is poorly understood. To identify genes controlling carbohydrate partitioning, we isolated mutants that are defective in exporting fixed carbon from leaves. Here we describe psychedelic (psc), a new mutant of maize (Zea mays) that is perturbed in carbohydrate partitioning. psc mutants exhibit stable, discrete chlorotic and green regions within their leaves. psc chlorotic tissues hyperaccumulate starch and soluble sugars, while psc green tissues appear comparable to wild-type leaves. The psc chlorotic and green tissue boundaries are usually delineated by larger veins, suggesting that translocation of a mobile compound through the veins may influence the tissue phenotype. psc mutants display altered biomass partitioning, which is consistent with reduced carbohydrate export from leaves to developing tissues. We determined that the psc mutation is unlinked to previously characterized maize leaf carbohydrate hyperaccumulation mutants. Additionally, we found that the psc mutant phenotype is inherited as a recessive, duplicate-factor trait in some inbred lines. Genetic analyses with other maize mutants with variegated leaves and impaired carbohydrate partitioning suggest that Psc defines an independent pathway. Therefore, investigations into the psc mutation have uncovered two previously unknown genes that redundantly function to regulate carbohydrate partitioning in maize.

Comparative sequence analyses indicate that Coffea (Asterids) and Vitis (Rosids) derive from the same paleo-hexaploid ancestral genome

The complete sequence of Vitis vinifera revealed that the rosid clade derives from a hexaploid ancestor. At present, no analysis of complete genome sequence is available for an asterid, the other large eudicot clade, which includes the economically important species potato, tomato and coffee. To elucidate the genomic history of asterids, we compared the sequence of an 800 kb region of diploid Coffea genome to the orthologous regions of V. vinifera, Populus trichocarpa and Arabidopsis thaliana. We found a very high level of collinearity between around 80 genes of the three rosid species and Coffea. Collinearity comparisons between orthologous and paralogous regions indicates that (1) the Coffea (and consequently all asterids) and rosids share the same hexaploid ancestor; (2) the diploidization process (loss of duplicated and redundant copies from the whole genome duplication) was very advanced in the most recent common ancestor of rosids and asterids. Finally, no additional polyploidization events were detected in the Coffea lineage. Differences in gene loss rates were detected among the three rosid species and linked to the divergence in protein sequences.

Universal global imprints of genome growth and evolution--equivalent length and cumulative mutation density

Background

Segmental duplication is widely held to be an important mode of genome growth and evolution. Yet how this would affect the global structure of genomes has been little discussed.

Methods/Principal Findings

Here, we show that equivalent length, or , a quantity determined by the variance of fluctuating part of the distribution of the -mer frequencies in a genome, characterizes the latter's global structure. We computed the s of 865 complete chromosomes and found that they have nearly universal but (-dependent) values. The differences among the of a chromosome and those of its coding and non-coding parts were found to be slight.

Conclusions

We verified that these non-trivial results are natural consequences of a genome growth model characterized by random segmental duplication and random point mutation, but not of any model whose dominant growth mechanism is not segmental duplication. Our study also indicates that genomes have a nearly universal cumulative “point” mutation density of about 0.73 mutations per site that is compatible with the relatively low mutation rates of (15)10/site/Mya previously determined by sequence comparison for the human and E. coli genomes.

Dating the origins of polyploidy events

is a widespread speciation mechanism, particularly in plants. Estimating the time of origin of polyploid species is important for understanding issues such as gene loss and changes in regulation and expression among homoeologous copies that coexist in a single genome owing to polyploidy. Polyploid species can originate in various ways; the effects of mode of origin, genetic system, and sampling on estimates of the age of polyploid origin using distances between alleles of polyploids and their diploid progenitors, or between homoeologous loci in a polyploid genome, are explored. Even in the simplest cases, simulations confirm that different loci are expected to give very different estimates of the date of origin. The time of polyploid origin is at least as old as the time estimated from comparison of an allele sampled from the polyploid with the most closely related allele in the diploid progenitor. The polyploidy literature often does not make clear the longstanding observation that the divergence of homoeologous copies in an allopolyploid tracks the divergence of diploid species, not the origin of the polyploid. Estimating the date of origin of a polyploid is difficult, and in some circumstances impossible. Skepticism about dates of polyploid origins is clearly warranted.

Making a functional diploid: from polysomic to disomic inheritance

One little understood feature of polyploid speciation is the transition from polysomic to disomic inheritance, and much recent attention has focused on the role of pairing genes in this process. Using computer simulations we studied the effects of mutations, chromosomal inversions, chiasma, neofunctionalization, subfunctionalization and selection on the evolution of disomic inheritance in a polyploid over 10 000 generations. We show that: the evolution of pairing genes is not essential for the establishment of disomic inheritance, as genetic drift, coupled with a threshold for homologue pairing fidelity, is sufficient to explain the transition from polysomic to disomic inheritance; high rates of recombination increase the number of generations required for disomic inheritance to become established; both neofunctionalization and subfunctionalization speed up the transition to disomic inheritance. The data suggest that during polyploid species establishment, selection will favour reduced chiasma number and/or more focused distribution. The data also suggest a new role for subfunctionalization in that it can drive disomic inheritance. The evolution of subfunctionalization in genes across the genome will then act to maintain genes in syntenic blocks and may explain why such regions are so highly conserved.

Genomic and expression plasticity of polyploidy

Polyploidy or whole genome duplication (WGD) occurs throughout the evolutionary history of many plants and some animals, including crops such as wheat, cotton, and sugarcane. Recent studies have documented rapid and dynamic changes in genomic structure and gene expression in plant polyploids, which reflects genomic and functional plasticity of duplicate genes and genomes in plants. Common features of uniparental gene regulation and nonadditive gene expression in regulatory pathways responsive to growth, development, and stresses in many polyploids have led to the conclusion that epigenetic mechanisms including chromatin modifications and small RNAs play central roles in shaping molecular and phenotypic novelty that may be selected and domesticated in many polyploid plants and crops.

Insights from the comparison of plant genome sequences

The next decade will see essentially completed sequences for multiple branches of virtually all angiosperm clades that include major crops and/or botanical models. These sequences will provide a powerful framework for relating genome-level events to aspects of morphological and physiological variation that have contributed to the colonization of much of the planet by angiosperms. Clarification of the fundamental angiosperm gene set, its arrangement, lineage-specific variations in gene repertoire and arrangement, and the fates of duplicated gene pairs will advance knowledge of functional and regulatory diversity and perhaps shed light on adaptation by lineages to whole-genome duplication, which is a distinguishing feature of angiosperm evolution. Better understanding of the relationships among angiosperm genomes promises to provide a firm foundation upon which to base translational genomics: the leveraging of hard-won structural and functional genomic information from crown botanical models to dissect novel and, in some cases, economically important features in many additional organisms.

A genome-wide characterization of microRNA genes in maize

MicroRNAs (miRNAs) are small, non-coding RNAs that play essential roles in plant growth, development, and stress response. We conducted a genome-wide survey of maize miRNA genes, characterizing their structure, expression, and evolution. Computational approaches based on homology and secondary structure modeling identified 150 high-confidence genes within 26 miRNA families. For 25 families, expression was verified by deep-sequencing of small RNA libraries that were prepared from an assortment of maize tissues. PCR-RACE amplification of 68 miRNA transcript precursors, representing 18 families conserved across several plant species, showed that splice variation and the use of alternative transcriptional start and stop sites is common within this class of genes. Comparison of sequence variation data from diverse maize inbred lines versus teosinte accessions suggest that the mature miRNAs are under strong purifying selection while the flanking sequences evolve equivalently to other genes. Since maize is derived from an ancient tetraploid, the effect of whole-genome duplication on miRNA evolution was examined. We found that, like protein-coding genes, duplicated miRNA genes underwent extensive gene-loss, with approximately 35% of ancestral sites retained as duplicate homoeologous miRNA genes. This number is higher than that observed with protein-coding genes. A search for putative miRNA targets indicated bias towards genes in regulatory and metabolic pathways. As maize is one of the principal models for plant growth and development, this study will serve as a foundation for future research into the functional roles of miRNA genes.

Maize centromere structure and evolution: sequence analysis of centromeres 2 and 5 reveals dynamic Loci shaped primarily by retrotransposons

We describe a comprehensive and general approach for mapping centromeres and present a detailed characterization of two maize centromeres. Centromeres are difficult to map and analyze because they consist primarily of repetitive DNA sequences, which in maize are the tandem satellite repeat CentC and interspersed centromeric retrotransposons of maize (CRM). Centromeres are defined epigenetically by the centromeric histone H3 variant, CENH3. Using novel markers derived from centromere repeats, we have mapped all ten centromeres onto the physical and genetic maps of maize. We were able to completely traverse centromeres 2 and 5, confirm physical maps by fluorescence in situ hybridization (FISH), and delineate their functional regions by chromatin immunoprecipitation (ChIP) with anti-CENH3 antibody followed by pyrosequencing. These two centromeres differ substantially in size, apparent CENH3 density, and arrangement of centromeric repeats; and they are larger than the rice centromeres characterized to date. Furthermore, centromere 5 consists of two distinct CENH3 domains that are separated by several megabases. Succession of centromere repeat classes is evidenced by the fact that elements belonging to the recently active recombinant subgroups of CRM1 colonize the present day centromeres, while elements of the ancestral subgroups are also found in the flanking regions. Using abundant CRM and non-CRM retrotransposons that inserted in and near these two centromeres to create a historical record of centromere location, we show that maize centromeres are fluid genomic regions whose borders are heavily influenced by the interplay of retrotransposons and epigenetic marks. Furthermore, we propose that CRMs may be involved in removal of centromeric DNA (specifically CentC), invasion of centromeres by non-CRM retrotransposons, and local repositioning of the CENH3.

Centromeres tend to be the last regions to be assembled in genome projects, as their mapping is hampered by their characteristically high repeat DNA content and lack of genetic recombination. Using unique markers derived from these repeat-rich regions, we were able to generate and annotate physical maps of two maize centromeres. Functional centromeres are defined not so much by their primary DNA sequence as by the presence of CENH3, a special histone that replaces canonical histone H3 in centromeric nucleosomes. Little is known about how deposition of CENH3 is regulated, or about the interplay between centromeric repeats and CENH3. By graphing the density of CENH3 nucleosomes onto the physical map, we delineated the functional centromeres in today's maize genome. We then used the large number of LTR retrotransposon insertions, for which the corn genome is well known, as “archeological evidence” to reconstruct the historic centromere boundaries. This was possible because i) some retrotransposon families of maize (CRM) appear to possess a unique ability to preferentially target centromeres during integration and ii) insertion times of individual retrotransposons can be calculated. Here we show that the centromere boundaries in maize have changed over time and are heavily influenced by centromeric and non-centromeric repeats.

A single molecule scaffold for the maize genome

About 85% of the maize genome consists of highly repetitive sequences that are interspersed by low-copy, gene-coding sequences. The maize community has dealt with this genomic complexity by the construction of an integrated genetic and physical map (iMap), but this resource alone was not sufficient for ensuring the quality of the current sequence build. For this purpose, we constructed a genome-wide, high-resolution optical map of the maize inbred line B73 genome containing >91,000 restriction sites (averaging 1 site/∼23 kb) accrued from mapping genomic DNA molecules. Our optical map comprises 66 contigs, averaging 31.88 Mb in size and spanning 91.5% (2,103.93 Mb/∼2,300 Mb) of the maize genome. A new algorithm was created that considered both optical map and unfinished BAC sequence data for placing 60/66 (2,032.42 Mb) optical map contigs onto the maize iMap. The alignment of optical maps against numerous data sources yielded comprehensive results that proved revealing and productive. For example, gaps were uncovered and characterized within the iMap, the FPC (fingerprinted contigs) map, and the chromosome-wide pseudomolecules. Such alignments also suggested amended placements of FPC contigs on the maize genetic map and proactively guided the assembly of chromosome-wide pseudomolecules, especially within complex genomic regions. Lastly, we think that the full integration of B73 optical maps with the maize iMap would greatly facilitate maize sequence finishing efforts that would make it a valuable reference for comparative studies among cereals, or other maize inbred lines and cultivars.

The maize genome contains abundant repeats interspersed by low-copy, gene-coding sequences that make it a challenge to sequence; consequently, current BAC sequence assemblies average 11 contigs per clone. The iMap deals with such complexity by the judicious integration of IBM genetic and B73 physical maps, but the B73 genome structure could differ from the IBM population because of genetic recombination and subsequent rearrangements. Accordingly, we report a genome-wide, high-resolution optical map of maize B73 genome that was constructed from the direct analysis of genomic DNA molecules without using genetic markers. The integration of optical and iMap resources with comparisons to FPC maps enabled a uniquely comprehensive and scalable assessment of a given BAC's sequence assembly, its placement within a FPC contig, and the location of this FPC contig within a chromosome-wide pseudomolecule. As such, the overall utility of the maize optical map for the validation of sequence assemblies has been significant and demonstrates the inherent advantages of single molecule platforms. Construction of the maize optical map represents the first physical map of a eukaryotic genome larger than 400 Mb that was created de novo from individual genomic DNA molecules.

The physical and genetic framework of the maize B73 genome

Maize is a major cereal crop and an important model system for basic biological research. Knowledge gained from maize research can also be used to genetically improve its grass relatives such as sorghum, wheat, and rice. The primary objective of the Maize Genome Sequencing Consortium (MGSC) was to generate a reference genome sequence that was integrated with both the physical and genetic maps. Using a previously published integrated genetic and physical map, combined with in-coming maize genomic sequence, new sequence-based genetic markers, and an optical map, we dynamically picked a minimum tiling path (MTP) of 16,910 bacterial artificial chromosome (BAC) and fosmid clones that were used by the MGSC to sequence the maize genome. The final MTP resulted in a significantly improved physical map that reduced the number of contigs from 721 to 435, incorporated a total of 8,315 mapped markers, and ordered and oriented the majority of FPC contigs. The new integrated physical and genetic map covered 2,120 Mb (93%) of the 2,300-Mb genome, of which 405 contigs were anchored to the genetic map, totaling 2,103.4 Mb (99.2% of the 2,120 Mb physical map). More importantly, 336 contigs, comprising 94.0% of the physical map (∼1,993 Mb), were ordered and oriented. Finally we used all available physical, sequence, genetic, and optical data to generate a golden path (AGP) of chromosome-based pseudomolecules, herein referred to as the B73 Reference Genome Sequence version 1 (B73 RefGen_v1).

Maize has been a cultural icon and staple food crop of Americans since the discovery of the new world in 1492. Contemporary society is now faced with growing demands for food and fuel in the face of global climate change and the potential for increased disease pressure. To provide a comprehensive foundation to systematically understand maize biology with the goal of breeding higher yielding, disease-resistant, and drought-tolerant cultivars, our consortium sequenced the B73 genome of maize. In this study, we used a comprehensive physical and genetic framework map to develop a minimum tiling path (MTP) of over 16,000 BAC clones across the genome. The MTP was generated dynamically and integrated numerous data types, such as in-coming genome sequence, over 8,000 sequence-based genetic markers, and the maize optical map. This allowed us to genetically anchor, order, and orient the majority of the maize physical map and genome sequence to the genetic map. Post-genome sequencing, we constructed a golden path (AGP) of sequence-based pseudomolecules representing the ten chromosomes of the maize B73 genome (B73 RefGen_v1). This unprecedented integration of genetic, physical, and genomic sequence into one framework will greatly facilitate all aspects of plant biological research.