Identification of the genomic locations of duplicate nucleotide sequences in maize by analysis of restriction fragment length polymorphisms

Locus-specific paramutation in Zea mays is maintained by a PICKLE-like chromodomain helicase DNA-binding 3 protein controlling development and male gametophyte function

Paramutations represent directed and meiotically-heritable changes in gene regulation leading to apparent violations of Mendelian inheritance. Although the mechanism and evolutionary importance of paramutation behaviors remain largely unknown, genetic screens in maize (Zea mays) identify five components affecting 24 nucleotide RNA biogenesis as required to maintain repression of a paramutant purple plant1 (pl1) allele. Currently, the RNA polymerase IV largest subunit represents the only component also specifying proper development. Here we identify a chromodomain helicase DNA-binding 3 (CHD3) protein orthologous to Arabidopsis (Arabidopsis thaliana) PICKLE as another component maintaining both pl1 paramutation and normal somatic development but without affecting overall small RNA biogenesis. In addition, genetic tests show this protein contributes to proper male gametophyte function. The similar mutant phenotypes documented in Arabidopsis and maize implicate some evolutionarily-conserved gene regulation while developmental defects associated with the two paramutation mutants are largely distinct. Our results show that a CHD3 protein responsible for normal plant ontogeny and sperm transmission also helps maintain meiotically-heritable epigenetic regulatory variation for specific alleles. This finding implicates an intersection of RNA polymerase IV function and nucleosome positioning in the paramutation process.

Genes are switched “on” and “off” during normal development by regulating DNA accessibility within the chromosomes. How certain gene variants permanently maintain “off” states from one generation to the next remains unclear, but studies in multiple eukaryotes implicate roles for specific types of small RNAs, some of which define cytosine methylation patterns. In corn, these RNAs come from at least two RNA polymerase II-derived complexes sharing a common catalytic subunit (RPD1). Although RPD1 both controls the normal developmental switching of many genes and permanently maintains some of these “off” states across generations, how RPD1 function defines heritable DNA accessibility is unknown. We discovered that a protein (CHD3a) belonging to a group known to alter nucleosome positioning is also required to help maintain a heritable “off” state for one particular corn gene variant controlling both plant and flower color. We also found CHD3a necessary for normal plant development and sperm transmission consistent with the idea that proper nucleosome positioning defines evolutionarily-important gene expression patterns. Because both CHD3a and RPD1 maintain the heritable “off” state of a specific gene variant, their functions appear to be mechanistically linked.

A Malvaceae mystery: A mallow maelstrom of genome multiplications and maybe misleading methods?

Previous research suggests that Gossypium has undergone a 5- to 6-fold multiplication following its divergence from Theobroma. However, the number of events, or where they occurred in the Malvaceae phylogeny remains unknown. We analyzed transcriptomic and genomic data from representatives of eight of the nine Malvaceae subfamilies. Phylogenetic analysis of nuclear data placed Dombeya (Dombeyoideae) as sister to the rest of Malvadendrina clade, but the plastid DNA tree strongly supported Durio (Helicteroideae) in this position. Intraspecific Ks plots indicated that all sampled taxa, except Theobroma (Byttnerioideae), Corchorus (Grewioideae), and Dombeya (Dombeyoideae), have experienced whole genome multiplications (WGMs). Quartet analysis suggested WGMs were shared by Malvoideae-Bombacoideae and Sterculioideae-Tilioideae, but did not resolve whether these are shared with each other or Helicteroideae (Durio). Gene tree reconciliation and Bayesian concordance analysis suggested a complex history. Alternative hypotheses are suggested, each involving two independent autotetraploid and one allopolyploid event. They differ in that one entails an allopolyploid origin for the Durio lineage, whereas the other invokes an allopolyploid origin for Malvoideae-Bombacoideae. We highlight the need for more genomic information in the Malvaceae and improved methods to resolve complex evolutionary histories that may include allopolyploidy, incomplete lineage sorting, and variable rates of gene and genome evolution.

Simple Sequence Repeat (SSR) Genetic Linkage Map of D Genome Diploid Cotton Derived from an Interspecific Cross between Gossypium davidsonii and Gossypium klotzschianum

The challenge in tetraploid cotton cultivars is the narrow genetic base and therefore, the bottleneck is how to obtain interspecific hybrids and introduce the germplasm directly from wild cotton to elite cultivars. Construction of genetic maps has provided insight into understanding the genome structure, interrelationships between organisms in relation to evolution, and discovery of genes that carry important agronomic traits in plants. In this study, we generated an interspecific hybrid between two wild diploid cottons, Gossypium davidsonii and Gossypium klotzschianum, and genotyped 188 F2:3 populations in order to develop a genetic map. We screened 12,560 SWU Simple Sequence Repeat (SSR) primers and obtained 1000 polymorphic markers which accounted for only 8%. A total of 928 polymorphic primers were successfully scored and only 728 were effectively linked across the 13 chromosomes, but with an asymmetrical distribution. The map length was 1480.23 cM, with an average length of 2.182 cM between adjacent markers. A high percentage of the markers on the map developed, and for the physical map of G. raimondii, exhibited highly significant collinearity, with two types of duplication. High level of segregation distortion was observed. A total of 27 key genes were identified with diverse roles in plant hormone signaling, development, and defense reactions. The achievement of developing the F2:3 population and its genetic map constructions may be a landmark in establishing a new tool for the genetic improvement of cultivars from wild plants in cotton. Our map had an increased recombination length compared to other maps developed from other D genome cotton species.

Unraveling the molecules hidden in the gray shadows of quantitative disease resistance to pathogens

Una de las preguntas más desafiantes del fitomejoramiento y de la fitopatología molecular es ¿cuáles son las bases genéticas y moleculares de la resistencia cuantitativa a enfermedades?. El escaso conocimiento de cómo este tipo de resistencia funciona ha obstaculizado que los fitomejoradores la aprovecharlo plenamente. Para superar estos obstáculos se han desarrollado nuevas metodologías para el estudio de rasgos cuantitativos. Los enfoques como el mapeo genético, la identificación de loci de rasgos cuantitativos (QTL) y el mapeo por asociaciones, incluyendo el enfoque de genes candidatos y los estudios de asociación amplia del genoma, se han llevado a cabo históricamente para describir rasgos cuantitativos y por lo tanto para estudiar QDR. Además, se han proporcionado grandes avances en la obtención de datos fenotípicos cuantitativos para mejorar estos análisis. Recientemente, algunos genes asociados a QDR han sido clonados, lo que conduce a nuevas hipótesis sobre las bases moleculares de este tipo de resistencia. En esta revisión presentamos los avances más recientes sobre QDR y la correspondiente aplicación, que han permitido postular nuevas ideas que pueden ayudar a construir nuevos modelos. Algunas de las hipótesis presentadas aquí como posibles explicaciones para QDR están relacionadas con el nivel de expresión y el splicing alternativo de algunos genes relacionados con la defensa, la acción de "alelos débiles" de genes R, la presencia de variantes alélicas en los genes implicados en la respuesta de defensa y un papel central de quinasas o pseudoqinasas. Con la información recapitulada en esta revisión es posible concluir que la distinción conceptual entre resistencia cualitativa y cuantitativa puede ser cuestionada ya que ambos comparten importantes componentes.

Competitive Ability of Maize Pollen Grains Requires Paralogous Serine Threonine Protein Kinases STK1 and STK2

serine threonine kinase1 (stk1) and serine threonine kinase2 (stk2) are closely related maize paralogous genes predicted to encode serine/threonine protein kinases. Pollen mutated in stk1 or stk2 competes poorly with normal pollen, pointing to a defect in pollen tube germination or growth. Both genes are expressed in pollen, but not in most other tissues. In germination media, STK1 and STK2 fluorescent fusion proteins localize to the plasma membrane of the vegetative cell. RNA-seq experiments identified 534 differentially expressed genes in stk1 mutant pollen relative to wild type. Gene ontology (GO) molecular functional analysis uncovered several differentially expressed genes with putative ribosome initiation and elongation functions, suggesting that stk1 might affect ribosome function. Of the two paralogs, stk1 may play a more important role in pollen development than stk2, as stk2 mutations have a smaller pollen transmission effect. However, stk2 does act as an enhancer of stk1 because the double mutant combination is only infrequently pollen-transmitted in double heterozygotes. We conclude that the stk paralogs play an essential role in pollen development.

Has the connection between polyploidy and diversification actually been tested?

Many major clades of angiosperms have several whole genome duplications (polyploidization events) in their distant past, suggesting that polyploidy drives or at least permits diversification. However, data on recently diverged groups are more equivocal, finding little evidence of elevated diversification following polyploidy. The discrepancy may be attributable at least in part to methodology. Many studies use indirect methods, such as chromosome numbers, genome size, and Ks plots, to test polyploidy, although these approaches can be misleading, and often lack sufficient resolution. A direct test of diversification following polyploidy requires a sequence-based approach that traces the history of nuclear genomes rather than species. These methods identify the point of coalescence of ancestral genomes, but may be misleading about the time and thus the extent of diversification. Limitations of existing methods mean that the connection between polyploidy and diversification has not been rigorously tested and remains unknown.

A small, differentially regulated family of farnesyl diphosphate synthases in maize (Zea mays) provides farnesyl diphosphate for the biosynthesis of herbivore-induced sesquiterpenes

Of the three functional FPPS identified in maize, fpps3 is induced by herbivory to produce FDP important for the formation of the volatile sesquiterpenes of plant defense. Sesquiterpenes are not only crucial for the growth and development of a plant but also for its interaction with the environment. The biosynthesis of sesquiterpenes proceeds over farnesyl diphosphate (FDP), which is either used as a substrate for protein prenylation, converted to squalene, or to volatile sesquiterpenes. To elucidate the regulation of sesquiterpene biosynthesis in maize, we identified and characterized the farnesyl diphosphate synthase (FPPS) gene family which consists of three genes. Synteny analysis indicates that fpps2 and fpps3 originate from a genome duplication in an ancient tetraploid ancestor. The three FPPSs encode active enzymes that produce predominantly FDP from the isopentenyl diphosphate and dimethylallyl diphosphate substrates. Only fpps1 and fpps3 are induced by elicitor treatment, but induced fpps1 levels are much lower and only increased to the amounts of fpps3 levels in intact leaves. Elicitor-induced fpps3 levels in leaves increase to more than 15-fold of background levels. In undamaged roots, transcript levels of fpps1 are higher than those of fpps3, but only fpps3 transcripts are induced in response to herbivory by Diabrotica virgifera virgifera. A kinetic of transcript abundance in response to herbivory in leaves provided further evidence that the regulation of fpps3 corresponds to that of tps23, a terpene synthase, that converts FDP to the volatile (E)-ß-caryophyllene. Our study indicates that the differential expression of fpps1 and fpps3 provides maize with FDP for both primary metabolism and terpene-based defenses. The expression of fpps3 seems to coincide with the herbivore-induced emission of volatile sesquiterpenes that were demonstrated to be important defense signals.

The role of reticulate evolution in creating innovation and complexity

Reticulate evolution encompasses processes that conflict with traditional Tree of Life efforts. These processes, horizontal gene transfer (HGT), gene and whole-genome duplications through allopolyploidization, are some of the main driving forces for generating innovation and complexity. HGT has a profound impact on prokaryotic and eukaryotic evolution. HGTs can lead to the invention of new metabolic pathways and the expansion and enhancement of previously existing pathways. It allows for organismal adaptation into new ecological niches and new host ranges. Although many HGTs appear to be selected for because they provide some benefit to their recipient lineage, other HGTs may be maintained by chance through random genetic drift. Moreover, some HGTs that may initially seem parasitic in nature can cause complexity to arise through pathways of neutral evolution. Another mechanism for generating innovation and complexity, occurring more frequently in eukaryotes than in prokaryotes, is gene and genome duplications, which often occur through allopolyploidizations. We discuss how these different evolutionary processes contribute to generating innovation and complexity.

The role of PIN auxin efflux carriers in polar auxin transport and accumulation and their effect on shaping maize development

In plants, proper seed development and the continuing post-embryonic organogenesis both require that different cell types are correctly differentiated in response to internal and external stimuli. Among internal stimuli, plant hormones and particularly auxin and its polar transport (PAT) have been shown to regulate a multitude of plant physiological processes during vegetative and reproductive development. Although our current auxin knowledge is almost based on the results from researches on the eudicot Arabidopsis thaliana, during the last few years, many studies tried to transfer this knowledge from model to crop species, maize in particular. Applications of auxin transport inhibitors, mutant characterization, and molecular and cell biology approaches, facilitated by the sequencing of the maize genome, allowed the identification of genes involved in auxin metabolism, signaling, and particularly in polar auxin transport. PIN auxin efflux carriers have been shown to play an essential role in regulating PAT during both seed and post-embryonic development in maize. In this review, we provide a summary of the recent findings on PIN-mediated polar auxin transport during maize development. Similarities and differences between maize and Arabidopsis are analyzed and discussed, also considering that their different plant architecture depends on the differentiation of structures whose development is controlled by auxins.

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.

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.

RNA-seq in grain unveils fate of neo- and paleopolyploidization events in bread wheat (Triticum aestivum L.)

Background

Whole genome duplication is a common evolutionary event in plants. Bread wheat (Triticum aestivum L.) is a good model to investigate the impact of paleo- and neoduplications on the organization and function of modern plant genomes.

Results

We performed an RNA sequencing-based inference of the grain filling gene network in bread wheat and identified a set of 37,695 non-redundant sequence clusters, which is an unprecedented resolution corresponding to an estimated half of the wheat genome unigene repertoire. Using the Brachypodium distachyon genome as a reference for the Triticeae, we classified gene clusters into orthologous, paralogous, and homoeologous relationships. Based on this wheat gene evolutionary classification, older duplicated copies (dating back 50 to 70 million years) exhibit more than 80% gene loss and expression divergence while recent duplicates (dating back 1.5 to 3 million years) show only 54% gene loss and 36 to 49% expression divergence.

Conclusions

We suggest that structural shuffling due to duplicated gene loss is a rapid process, whereas functional shuffling due to neo- and/or subfunctionalization of duplicates is a longer process, and that both shuffling mechanisms drive functional redundancy erosion. We conclude that, as a result of these mechanisms, half the gene duplicates in plants are structurally and functionally altered within 10 million years of evolution, and the diploidization process is completed after 45 to 50 million years following polyploidization.

B73-Mo17 near-isogenic lines demonstrate dispersed structural variation in maize

Recombinant inbred lines developed from the maize (Zea mays ssp. mays) inbreds B73 and Mo17 have been widely used to discover quantitative trait loci controlling a wide variety of phenotypic traits and as a resource to produce high-resolution genetic maps. These two parents were used to produce a set of near-isogenic lines (NILs) with small regions of introgression into both backgrounds. A novel array-based genotyping platform was used to score genotypes of over 7,000 loci in 100 NILs with B73 as the recurrent parent and 50 NILs with Mo17 as the recurrent parent. This population contains introgressions that cover the majority of the maize genome. The set of NILs displayed an excess of residual heterozygosity relative to the amount expected based on their pedigrees, and this excess residual heterozygosity is enriched in the low-recombination regions near the centromeres. The genotyping platform provided the ability to survey copy number variants that exist in more copies in Mo17 than in B73. The majority of these Mo17-specific duplications are located in unlinked positions throughout the genome. The utility of this population for the discovery and validation of quantitative trait loci was assessed through analysis of plant height variation.

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.

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.

The polychromatic Helitron landscape of the maize genome

Maize Helitron transposons are intriguing because of their notable ability to capture gene fragments and move them around the genome. To document more extensively their variability and their contribution to the remarkable genome structure variation of present-day maize, we have analyzed their composition, copy number, timing of insertion, and chromosomal distribution. First, we searched 2.4 Gb of sequences generated by the Maize Genome Sequencing Project with our HelitronFinder program. We identified 2,791 putative nonautonomous Helitrons and manually curated a subset of 272. The predicted Helitrons measure 11.9 kb on average and carry from zero to nine gene fragments, captured from 376 different genes. Although the diversity of Helitron gene fragments in maize is greater than in other species, more than one-third of annotated Helitrons carry fragments derived from just one of two genes. Most members in these two subfamilies inserted in the genome less than one million years ago. Second, we conducted a BLASTN search of the maize sequence database with queries from two previously described agenic Helitrons not detected by HelitronFinder. Two large subfamilies of Helitrons or Helitron-related transposons were identified. One subfamily, termed Cornucopious, consists of thousands of copies of an approximately 1.0-kb agenic Helitron that may be the most abundant transposon in maize. The second subfamily consists of >150 copies of a transposon-like sequence, termed Heltir, that has terminal inverted repeats resembling Helitron 3' termini. Nonautonomous Helitrons make up at least 2% of the maize genome and most of those tested show +/- polymorphisms among modern inbred lines.

Phenotypic, genetic and genomic consequences of natural and synthetic polyploidization of Nicotiana attenuata and Nicotiana obtusifolia

BACKGROUND AND METHODS

Polyploidy results in genetic turmoil, much of which is associated with new phenotypes that result in speciation. Five independent lines of synthetic allotetraploid N. x obtusiata (N x o) were created from crosses between the diploid N. attenuata (Na) (male) and N. obtusifolia (No) (female) and the autotetraploids of Na (NaT) and No (NoT) were synthesized. Their genetic, genomic and phenotypic changes were then compared with those of the parental diploid species (Na and No) as well as to the natural allotetraploids, N. quadrivalvis (Nq) and N. clevelandii (Nc), which formed 1 million years ago from crosses between ancient Na and No.

KEY RESULTS

DNA fingerprinting profiles (by UP-PCR) revealed that the five N x o lines shared similar but not identical profiles. Both synthetic and natural polyploidy showed a dosage effect on genome size (as measured in seeds); however, only Nq was associated with a genome upsizing. Phenotypic analysis revealed that at the cellular level, N x o lines had phenotypes intermediate of the parental phenotypes. Both allo- and autotetraploidization had a dosage effect on seed and dry biomass (except for NaT), but not on stalk height at first flower. Nc showed paternal (Na) cellular phenotypes but inherited maternal (No) biomass and seed mass, whereas Nq showed maternal (No) cellular phenotypes but inherited paternal (Na) biomass and seed mass patterns. Principal component analysis grouped Nq with N x o lines, due to similar seed mass, stalk height and genome size. These traits separated Nc, No and Na from Nq and N x o lines, whereas biomass distinguished Na from N x o and Nq lines, and NaT clustered closer to Nq and N x o lines than to Na.

CONCLUSIONS

Both allo- and autotetraploidy induce considerable morphological, genetic and genomic changes, many of which are retained by at least one of the natural polyploids. It is proposed that both natural and synthetic polyploids are well suited for studying the evolution of adaptive responses.

Floret-specific differences in gene expression and support for the hypothesis that tapetal degeneration of Zea mays L. occurs via programmed cell death

The maize (Zea mays) spikelet consists of two florets, each of which contains three developmentally synchronized anthers. Morphologically, the anthers in the upper and lower florets proceed through apparently similar developmental programs. To test for global differences in gene expression and to identify genes that are coordinately regulated during maize anther development, RNA samples isolated from upper and lower floret anthers at six developmental stages were hybridized to cDNA microarrays. Approximately 9% of the tested genes exhibited statistically significant differences in expression between anthers in the upper and lower florets. This finding indicates that several basic biological processes are differentially regulated between upper and lower floret anthers, including metabolism, protein synthesis and signal transduction. Genes that are coordinately regulated across anther development were identified via cluster analysis. Analysis of these results identified stage-specific, early in development, late in development and bi-phasic expression profiles. Quantitative RT-PCR analysis revealed that four genes whose homologs in other plant species are involved in programmed cell death are up-regulated just prior to the time the tapetum begins to visibly degenerate (i.e., the mid-microspore stage). This finding supports the hypothesis that developmentally normal tapetal degeneration occurs via programmed cell death.

The 'inner circle' of the cereal genomes

Early marker-based macrocolinearity studies between the grass genomes led to arranging their chromosomes into concentric 'crop circles' of synteny blocks that initially consisted of 30 rice-independent linkage groups representing the ancestral cereal genome structure. Recently, increased marker density and genome sequencing of several cereal genomes allowed the characterization of intragenomic duplications and their integration with intergenomic colinearity data to identify paleo-duplications and propose a model for the evolution of the grass genomes from a common ancestor. On the basis of these data an 'inner circle' comprising five ancestral chromosomes was defined providing a new reference for the grass chromosomes and new insights into their ancestral relationships and origin, as well as an efficient tool to design cross-genome markers for genetic studies.

Development of genic-microsatellite markers for sorghum staygreen QTL using a comparative genomic approach with rice

The already available comprehensive genome sequence information of model crops along with the transcriptomic resource from other crops provides an excellent opportunity for comparative genome analysis. We studied the synteny between each of the four major sorghum staygreen quantitative trait loci (QTL) regions with that in the rice genome and attempted to increase marker density around the QTL with genic-microsatellites from the sorghum transcriptomic resource using the rice genome as template. For each of the sorghum QTL regions, the reported RFLP markers were compiled, used for sequence similarity searches against the rice genome which identified syntenous regions on rice chromosome 1 for Stg1 and Stg2 QTL, on chromosome 9 for Stg3 QTL, and on chromosome 11 for Stg4 QTL. Using the Gramene genome browsing tool, 869 non-redundant sorghum expressed sequence tags (ESTs) were selected and 50 genic-microsatellites (18, 12, 15, and 5, for Stg1, Stg2, Stg3, and Stg4 QTL, respectively) could be developed. We could experimentally establish synteny of the Stg1, Stg2, Stg3, and Stg4 QTL regions with that of the rice genome by mapping ten polymorphic genic-microsatellite markers (20%) to the positions of the staygreen QTL. The simple strategy demonstrated in the present study could readily be extrapolated to other cereals of the Poaceae family. The markers developed in this study provide a basis for the isolation of genes underling these QTL using an association study or map-based gene isolation approach, and create an additional option for MAS of the staygreen trait in sorghum.

Characterization of the monoterpene synthase gene tps26, the ortholog of a gene induced by insect herbivory in maize

Plants damaged by insects can synthesize and release volatile chemicals that attract natural enemies of the herbivore. The maize (Zea mays subsp. mays) terpene synthase gene stc1 is part of that indirect defense response, being induced in seedling blades in response to herbivory by beet army worm. Many genes in maize are duplicated because of a past whole-genome duplication event, and several of these orthologs display different expression patterns. We report here the isolation and characterization of tps26 and confirm by homology and synteny criteria that it is the ortholog of stc1. Prior genetic analysis revealed that the stc1 function is not duplicated, raising the interesting question of how the two orthologs have become differentiated in their expression. tps26 encodes a 633-amino acid protein that is highly conserved with STC1. Like stc1, tps26 is induced by wounding, but in the roots and leaf sheath, instead of the blade, and not in response to beet army worm feeding. tps26 maps near a quantitative trait locus for Southwestern corn borer resistance, making it a plausible candidate gene for that quantitative trait locus. However, while possessing highly polymorphic tps26 alleles, the resistant and susceptible parents of the mapping population do not differ in levels of tps26 expression. Moreover, tps26 is not induced specifically by Southwestern corn borer feeding. Therefore, although they share a wounding response, the stc1 and tps26 maize orthologs differ in their tissue specificity and their induction by insect herbivores. The N termini of STC1 and TPS26 are predicted to encode plastid transit peptides; fusion proteins of green fluorescent protein to either N terminus localized to the plastid, confirming that prediction. The mature proteins, but not the respective complete proteins, were active and synthesized a blend of monoterpenes, indicating that they are monoterpene synthases. A gene closely related to stc1/tps26 is found in the sorghum (Sorghum spp.) genome at a location that is not orthologous with stc1. The possible origin of stc1-like genes is discussed.

Transcript profiling by 3'-untranslated region sequencing resolves expression of gene families

Differences in gene expression underlie central questions in plant biology extending from gene function to evolutionary mechanisms and quantitative traits. However, resolving expression of closely related genes (e.g. alleles and gene family members) is challenging on a genome-wide scale due to extensive sequence similarity and frequently incomplete genome sequence data. We present a new expression-profiling strategy that utilizes long-read, high-throughput sequencing to capture the information-rich 3'-untranslated region (UTR) of messenger RNAs (mRNAs). Resulting sequences resolve gene-specific transcripts independent of a sequenced genome. Analysis of approximately 229,000 3'-anchored sequences from maize (Zea mays) ovaries identified 14,822 unique transcripts represented by at least two sequence reads. Total RNA from ovaries of drought-stressed wild-type and viviparous-1 mutant plants was used to construct a multiplex cDNA library. Each sample was labeled by incorporating one of 16 unique three-base key codes into the 3'-cDNA fragments, and combined samples were sequenced using a GS 20 454 instrument. Transcript abundance was quantified by frequency of sequences identifying each unique mRNA. At least 202 unique transcripts showed highly significant differences in abundance between wild-type and mutant samples. For a subset of mRNAs, quantitative differences were validated by real-time reverse transcription-polymerase chain reaction. The 3'-UTR profile resolved 12 unique cellulose synthase (CesA) transcripts in maize ovaries and identified previously uncharacterized members of a histone H1 gene family. In addition, this method resolved nearly identical paralogs, as illustrated by two auxin-repressed, dormancy-associated (Arda) transcripts, which showed reciprocal mRNA abundance in wild-type and mutant samples. Our results demonstrate the potential of 3'-UTR profiling for resolving gene- and allele-specific transcripts.

Distinct mechanisms govern the dosage-dependent and developmentally regulated resistance conferred by the maize Hm2 gene

The maize Hm2 gene provides protection against the leaf spot and ear mold disease caused by Cochliobolus carbonum race 1 (CCR1). In this regard, it is similar to Hm1, the better-known disease resistance gene of the maize-CCR1 pathosystem. However, in contrast to Hm1, which provides completely dominant resistance at all stages of plant development, Hm2-conferred resistance is only partially dominant and becomes fully effective only at maturity. To investigate why Hm2 behaves in this manner, we cloned it on the basis of its homology to Hm1. As expected, Hm2 is a duplicate of Hm1, although the protein it encodes is grossly truncated compared with HM1. The efficacy of Hm2 in conferring resistance improves gradually over time, changing from having little or no impact in seedling tissues to providing complete immunity at anthesis. The developmentally specified phenotype of Hm2 is not dictated transcriptionally, because the expression level of the gene, whether occurring constitutively or undergoing substantial and transient induction in response to infection, does not change with plant age. In contrast, however, the Hm2 transcript is much more abundant in plants homozygous for this gene compared with plants that contain only one copy of the gene, suggesting a transcriptional basis for the dosage-dependent nature of Hm2. Thus, different mechanisms seem to underlie the developmentally programmed versus the partially dominant resistance phenotype of Hm2.

Linkage mapping of domestication loci in a large maize teosinte backcross resource

An ultimate objective of QTL mapping is cloning genes responsible for quantitative traits. However, projects seldom go beyond segments <5 cM without subsequent breeding and genotyping lines to identify additional crossovers in a genomic region of interest. We report on a QTL analysis performed as a preliminary step in the development of a resource for map-based cloning of domestication and improvement genes in corn. A large backcross (BC)1 population derived from a cross between maize (Zea mays ssp. mays) and teosinte (ssp. parviglumis) was grown for the analysis. A total of 1749 progenies were genotyped for 304 markers and measured for 22 morphological traits. The results are in agreement with earlier studies showing a small number of genomic regions having greater impact on the morphological traits distinguishing maize and teosinte. Despite considerable power to detect epistasis, few QTL interactions were identified. To create a permanent resource, seed of BC1 plants was archived and 1000 BC2S6 BC1-derived lines are in development for fine mapping and cloning. The identification of four BC1 progeny with crossovers in a single gene, tb1, indicated that enough derived lines already exist to clone many QTL without the need to generate and identify additional crossovers.

Genetic diversity of African maize inbred lines revealed by SSR markers

Knowledge of genetic diversity (GD) and relationships among maize inbred lines is indispensable in a breeding program. Our objectives were to (1) investigate the level of genetic diversity among maize inbred lines and (2) assess their genetic structures by applying simple sequence repeat (SSR) markers. Fifty-six highland and mid-altitude maize inbred lines obtained from CIMMYT programs in Ethiopia and Zimbabwe were genotyped using 27 SSR loci. All of the genotypes studied could unequivocally be distinguished with the combination of the SSRs used. In total, 104 SSR alleles were identified, with a mean of 3.85 alleles per locus. The average polymorphism information content (PIC) was 0.58. GD expressed as Euclidean distance, varied from 0.28 to 0.73 with an average of 0.59. Cluster analysis using unweighted pair group method with arithmetic average (UPGMA) suggested five groups among the inbred lines. Most of the inbred lines adapted to the highlands and the mid-altitudes were positioned in different clusters with a few discrepancies. The pattern of groupings of the inbred lines was mostly consistent with available pedigree information. The variability detected using SSR markers could potentially contribute towards effective utilization of the inbred lines for the exploitation of heterosis and formation of genetically diverse source populations in Ethiopian maize improvement programs.

Cytological diploidization and rapid genome changes of the newly synthesized allotetraploids Cucumis x hytivus

We used a newly synthesized allotetraploid between C. sativus (2n = 2x = 14, n gametic chromosome number, x haploid chromosome number) and C. hystrix (2n = 2x = 24) to study the genomic events in its early generations. Results from cytological characterization of the F(1) and the allotetraploid progenies showed that the rate of bivalents in meiotic metaphase I of the F(1) was greatly improved by chromosome doubling, and further improved during the selfing process of allopolyploid resulting into relatively diploid-like meiosis. Extensive genomic changes were detected by amplified fragment length polymorphism analysis. The changes mainly involved loss of parental restriction fragments and gaining of novel fragments. The total detectable changes were from 11.1 to 32.1%, and the frequency of losing parental fragments was much higher than that of gaining novel fragments. Some of the changes were initiated as early as in the F(1) hybrid, whereas others occurred after chromosome doubling (polyploid formation). No significant differences were detected in the reciprocal F(1) hybrids and S(0) generations. But the data showed that the frequency of sequence losing in C. sativus was about two times higher than in the C. hystrix. Our results demonstrated that the sequence elimination was the major event of genomic changes, and it might provide the physical basis for the diploid-like meiotic behavior in the diploidization of the newly formed allopolyploids. Moreover, the results suggest that the sequence elimination was not caused by cytoplasmic factors, and might relate to genomic recombination and to the numbers of parental chromosome.

Annotation and expression profile analysis of 2073 full-length cDNAs from stress-induced maize (Zea mays L.) seedlings

Full-length cDNAs are very important for genome annotation and functional analysis of genes. The number of full-length cDNAs from maize (Zea mays L.) remains limited. Here we report the construction of a full-length enriched cDNA library from osmotically stressed maize seedlings by using the modified CAP trapper method. From this library, 2073 full-length cDNAs were collected and further analyzed by sequencing from both the 5'- and 3'-ends. A total of 1728 (83.4%) sequences did not match known maize mRNA and full-length cDNA sequences in the GenBank database and represent new full-length genes. After alignment of the 2073 full-length cDNAs with 448 maize BAC sequences, it was found that 84 full-length cDNAs could be mapped to the BACs. Of these, 43 genes (51.2%) have been correctly annotated from the BAC clones, 37 genes (44.0%) have been annotated with a different exon-intron structure from our cDNA, and four genes (4.76%) had no annotations in the TIGR database. Expression analysis of 2073 full-length maize cDNAs using a cDNA macroarray led to the identification of 79 genes upregulated by stress treatments and 329 downregulated genes. Of the 79 stress-inducible genes, 30 genes contain ABRE, DRE, MYB, MYC core sequences or other abiotic-responsive cis-acting elements in their promoters. These results suggest that these cis-acting elements and the corresponding transcription factors take part in plant responses to osmotic stress either cooperatively or independently. Additionally, the data suggest that an ethylene signaling pathway may be involved in the maize response to drought stress.

Evolution of proteinase inhibitor defenses in North American allopolyploid species of Nicotiana

We studied the jasmonate (JA)-elicited trypsin-proteinase inhibitor (TPI) anti-herbivore defense system in North American Nicotiana to understand how complex polygenetic traits evolve after allopolyploidy speciation. N. quadrivalvis (Nq) and N. clevelandii (Nc) are allotetraploid descendant species of the ancestors of the diploid species N. attenuata (Na) and N. obtusifolia (No). From cDNA, intron and promoter sequence analyses, and Southern blotting, we deduced that only the maternally derived No TPI genes were retained in the tetraploid genomes (Nq, Nc), whereas the sequences of the paternal Na ancestor were deleted. The number of TPI repeats in different Nicotiana taxa was independent of phylogenetic associations. In Na, TPI activity and mRNA transcript accumulation as well as JA levels increased dramatically above wound-induced levels when the oral secretions (OS) from Manduca sexta larvae were introduced into wounds. This OS-mediated amplification of defense signaling and downstream response was also found in the tetraploid genomes but was absent from No; in No, OS treatment suppresses TPI mRNA accumulation and activity and does not increase JA accumulation. Hence, the tetraploids retained components of Na's signaling system, but lost Na's TPI genes and used No's TPI genes to retain a functional TPI defense system, underscoring the genomic flexibility that enables complex polygenic traits to be retained in allopolyploid species.

Genetic map-based analysis of genome structure in the homosporous fern Ceratopteris richardii

Homosporous ferns have extremely high chromosome numbers relative to flowering plants, but the species with the lowest chromosome numbers show gene expression patterns typical of diploid organisms, suggesting that they may be diploidized ancient polyploids. To investigate the role of polyploidy in fern genome evolution, and to provide permanent genetic resources for this neglected group, we constructed a high-resolution genetic linkage map of the homosporous fern model species, Ceratopteris richardii (n = 39). Linkage map construction employed 488 doubled haploid lines (DHLs) that were genotyped for 368 RFLP, 358 AFLP, and 3 isozyme markers. Forty-one linkage groups were recovered, with average spacing between markers of 3.18 cM. Most loci (approximately 76%) are duplicated and most duplicates occur on different linkage groups, indicating that as in other eukaryotic genomes, gene duplication plays a prominent role in shaping the architecture of fern genomes. Although past polyploidization is a potential mechanism for the observed abundance of gene duplicates, a wide range in the number of gene duplicates as well as the absence of large syntenic regions consisting of duplicated gene copies implies that small-scale duplications may be the primary mode of gene duplication in C. richardii. Alternatively, evidence of past polyploidization(s) may be masked by extensive chromosomal rearrangements as well as smaller-scale duplications and deletions following polyploidization(s).

Synteny conservation between the Prunus genome and both the present and ancestral Arabidopsis genomes

Background

Due to the lack of availability of large genomic sequences for peach or other Prunus species, the degree of synteny conservation between the Prunus species and Arabidopsis has not been systematically assessed. Using the recently available peach EST sequences that are anchored to Prunus genetic maps and to peach physical map, we analyzed the extent of conserved synteny between the Prunus and the Arabidopsis genomes. The reconstructed pseudo-ancestral Arabidopsis genome, existed prior to the proposed recent polyploidy event, was also utilized in our analysis to further elucidate the evolutionary relationship.

Results

We analyzed the synteny conservation between the Prunus and the Arabidopsis genomes by comparing 475 peach ESTs that are anchored to Prunus genetic maps and their Arabidopsis homologs detected by sequence similarity. Microsyntenic regions were detected between all five Arabidopsis chromosomes and seven of the eight linkage groups of the Prunus reference map. An additional 1097 peach ESTs that are anchored to 431 BAC contigs of the peach physical map and their Arabidopsis homologs were also analyzed. Microsyntenic regions were detected in 77 BAC contigs. The syntenic regions from both data sets were short and contained only a couple of conserved gene pairs. The synteny between peach and Arabidopsis was fragmentary; all the Prunus linkage groups containing syntenic regions matched to more than two different Arabidopsis chromosomes, and most BAC contigs with multiple conserved syntenic regions corresponded to multiple Arabidopsis chromosomes. Using the same peach EST datasets and their Arabidopsis homologs, we also detected conserved syntenic regions in the pseudo-ancestral Arabidopsis genome. In many cases, the gene order and content of peach regions was more conserved in the ancestral genome than in the present Arabidopsis region. Statistical significance of each syntenic group was calculated using simulated Arabidopsis genome.

Conclusion

We report here the result of the first extensive analysis of the conserved microsynteny using DNA sequences across the Prunus genome and their Arabidopsis homologs. Our study also illustrates that both the ancestral and present Arabidopsis genomes can provide a useful resource for marker saturation and candidate gene search, as well as elucidating evolutionary relationships between species.

Divergence of duplicated genes in maize: evolution of contrasting targeting information for enzymes in the porphyrin pathway

The divergence of sequence and expression pattern of duplicated genes provides a means for genetic innovation to occur without sacrificing an essential function. The cpx1 and cpx2 genes of maize are a singular example of duplicated genes that have diverged by deletion and creation of protein targeting information. The cpx genes encode coproporphyrinogen III oxidase ('coprogen oxidase'), which catalyzes a step in the synthesis of chlorophyll and heme. In plants, this enzyme has been found exclusively in the plastids. The cpx1 and cpx2 genes encode almost identical, catalytically active enzymes with distinctive N-terminal peptide sequences. The cpx1 gene encodes the expected plastid transit peptide, but this region is deleted from the cpx2 gene. While the 5' regions of both messenger RNAs are highly similar, the cpx2 gene has an open-reading frame that could encode a new targeting signal. GFP fused with CPX1 localized to the plastids. In contrast, the GFP fusion with CPX2 did not target plastids and appeared to localize to mitochondria. Both cpx genes are expressed ubiquitously but, based on mutant phenotype, they seem to have discrete biological roles. Seedlings homozygous for a null mutation in the cpx1 gene completely lack chlorophyll and develop necrotic lesions in the light. However, the mutant seedlings and callus cultures will grow in tissue culture in the dark, implying that they retain a capacity to produce heme. We discuss models for the evolution of the cpx genes and possible roles of mitochondrion-localized coprogen oxidase activity in maize.

Quantitative trait loci for cell wall components in recombinant inbred lines of maize (Zea mays L.) II: leaf sheath tissue

While maize silage is a significant feed component in animal production operations, little information is available on the genetic bases of fiber and lignin concentrations in maize, which are negatively correlated with digestibility. Fiber is composed largely of cellulose, hemicellulose and lignin, which are the primary components of plant cell walls. Variability for these traits in maize germplasm has been reported, but the sources of the variation and the relationships between these traits in different tissues are not well understood. In this study, 191 recombinant inbred lines of B73 (low-intermediate levels of cell wall components, CWCs) x De811 (high levels of CWCs) were analyzed for quantitative trait loci (QTL) associated with CWCs in the leaf sheath. Samples were harvested from plots at two locations in 1998 and one in 1999 and assayed for neutral detergent fiber (NDF), acid detergent fiber (ADF) and acid detergent lignin (ADL). QTL were detected on all ten chromosomes, most in tissue specific clusters in concordance with the high genotypic correlations for CWCs within the same tissue. Adjustment of NDF for its subfraction, ADF, revealed that most of the genetic variation in NDF was probably due to variation in ADF. The low to moderate genotypic correlations for the same CWC across leaf sheath and stalk tissues indicate that some genes for CWCs may only be expressed in certain tissues. Many of the QTL herein were detected in other populations, and some are linked to candidate genes for cell wall carbohydrate biosynthesis.

Cytogenetic mapping in maize

Cytogenetic maps depict the location and order of markers along chromosomes. Cytogenetic maps are important in genome research as they relate the genetic data and molecular sequences to the morphological features of chromosomes. In this paper, we discuss various methods used in cytogenetic mapping in maize, with special reference to fluorescence in situ hybridization (FISH) of single-copy sequences on meiotic pachytene chromosomes.

Allopolyploidy--a shaping force in the evolution of wheat genomes

Recent studies have shown that allopolyploidy accelerates genome evolution in wheat in two ways: (1) allopolyploidization triggers rapid genome changes (revolutionary changes) through the instantaneous generation of a variety of cardinal genetic and epigenetic alterations, and (2) the allopolyploid condition facilitates sporadic genomic changes during the life of the species (evolutionary changes) that are not attainable at the diploid level. The revolutionary changes comprise (1) non-random elimination of coding and non-coding DNA sequences, (2) epigenetic changes such as DNA methylation of coding and non-coding DNA leading, among others, to gene silencing, (3) activation of genes and retroelements which in turn alters the expression of adjacent genes. These highly reproducible changes occur in the F1 hybrids or in the first generation(s) of the nascent allopolyploids and were similar to those that occurred twice in nature: first in the formation of allotetraploid wheat (approximately 0.5 million years ago) and second in the formation of hexaploid wheat (approximately 10,000 years ago). Elimination of non-coding sequences from one of the two homoeologous pairs in tetraploids and from two homoeologous pairs in hexaploids, augments the differentiation of homoeologous chromosomes at the polyploid level, thus providing the physical basis for the diploid-like meiotic behavior of allopolyploid wheat. Regulation of gene expression may lead to improved inter-genomic interactions. Gene inactivation brings about rapid diploidization while activation of genes through demethylation or through transcriptional activation of retroelements altering the expression of adjacent genes, leads to novel expression patterns. The evolutionary changes comprise (1) horizontal inter-genomic transfer of chromosome segments between the constituent genomes, (2) production of recombinant genomes through hybridization and introgression between different allopolyploid species or, more seldom, between allopolyploids and diploids, and (3) mutations. These phenomena, emphasizing the plasticity of the genome with regards to both structure and function, might improve the adaptability of the newly formed allopolyploids and facilitate their rapid and successful establishment in nature.

Genome evolution of allopolyploids: a process of cytological and genetic diploidization

Allopolyploidy is a prominent mode of speciation in higher plants. Due to the coexistence of closely related genomes, a successful allopolyploid must have the ability to invoke and maintain diploid-like behavior, both cytologically and genetically. Recent studies on natural and synthetic allopolyploids have raised many discrepancies. Most species have displayed non-Mendelian behavior in the allopolyploids, but others have not. Some species have demonstrated rapid genome changes following allopolyploid formation, while others have conserved progenitor genomes. Some have displayed directed, non-random genome changes, whereas others have shown random changes. Some of the genomic changes have appeared in the F1 hybrids, which have been attributed to the union of gametes from different progenitors, while other changes have occurred during or after genome doubling. Although these observations provide significant novel insights into the evolution of allopolyploids, the overall mechanisms of the event are still elusive. It appears that both genetic and epigenetic operations are involved in the diploidization process of allopolyploids. Overall, genetic and epigenetic variations are often associated with the activities of repetitive sequences and transposon elements. Specifically, genomic sequence elimination and chromosome rearrangement are probably the major forces guiding cytological diploidization. Gene non-functionalization, sub-functionalization, neo-functionalization, as well as other kinds of epigenetic modifications, are likely the leading factors promoting genetic diploidization.

Gene enrichment in maize with hypomethylated partial restriction (HMPR) libraries

A new technology was developed to assist gene-enrichment sequencing of any complex plant genome, employing maize as the test organism. Hypomethylated partial restriction (HMPR) libraries were constructed by using independent partial digestions with methylation-sensitive restriction enzymes HpaII (5'-CCGG-3') and HpyCH4IV (5'-ACGT-3'). Fragments of 1-4 kb were purified and cloned, followed by sequence analysis of >2000 clones from 10 separate libraries. Organellar clones comprised approximately 10% of each library but were useful in showing that no chimeric clones were generated and that digestion efficiencies were 10%-25% in different libraries. Four separate HMPR libraries, analyzed in detail, exhibited very similar degrees of gene enrichment and repeat depletion. Known gene homologies were found in approximately 25% of the HMPR clones, compared with <4% in clones from a fully random set of unfiltered maize shotgun sequences. This six- to sevenfold enrichment for genes compares favorably with the best previous gene enrichment techniques in maize, High Cot analysis and methylation filtration. Compared with High Cot and methylation filtration, HMPR is exceptional in depleting retrotransposons' content to the lowest level yet observed (<5%, compared with >70% for unfiltered maize sequences) and in providing an unmatched enrichment for the "unknown" sequences that contain promoters, introns, and other gene-adjacent regions.

Cis-effects on meiotic recombination across distinct a1-sh2 intervals in a common Zea genetic background

Genetic distances across the a1-sh2 interval varied threefold in three near-isogenic stocks that carry structurally distinct teosinte A1 Sh2 haplotypes (from Z. mays spp. mexicana Chalco, Z. mays spp. parviglumis, and Z. luxurians) and a common maize a1::rdt sh2 haplotype. In each haplotype >85% of recombination events resolved in the proximal 10% of the approximately 130-kb a1-sh2 interval. Even so, significant differences in the distributions of recombination breakpoints were observed across subintervals among haplotypes. Each of the three previously detected recombination hot spots was detected in at least one of the three teosinte haplotypes and two of these hot spots were not detected in at least one teosinte haplotype. Moreover, novel hot spots were detected in two teosinte haplotypes. Due to the near-isogenic nature of the three stocks, the observed variation in the distribution of recombination events is the consequence of cis-modifications. Although generally negatively correlated with rates of recombination per megabase, levels of sequence polymorphisms do not fully account for the nonrandom distribution of recombination breakpoints. This study also suggests that estimates of linkage disequilibrium must be interpreted with caution when considering whether a gene has been under selection.

Quantitative trait loci for cell-wall components in recombinant inbred lines of maize (Zea mays L.) I: stalk tissue

Maize silage is a significant energy source for animal production operations, and the efficiency of the conversion of forage into animal mass is an important consideration when selecting cultivars for use as feed. Fiber and lignin are negatively correlated with digestibility of feed, so the development of forage with reduced levels of these cell-wall components (CWCs) is desirable. While variability for fiber and lignin is present in maize germplasm, traditional selection has focused on the yield of the ear rather than the forage quality of the whole plant, and little information is available concerning the genetics of fiber and lignin. The objectives of this study were to map quantitative trait loci (QTLs) for fiber and lignin in the maize stalk and compare them with QTLs from other populations. Stalk samples were harvested from 191 recombinant inbred lines (RILs) of B73 (an inbred line with low-to-intermediate levels of CWCs) x De811 (an inbred line with high levels of CWCs) at two locations in 1998 and one in 1999 and assayed for neutral detergent fiber (NDF), acid detergent fiber (ADF), and acid detergent lignin (ADL). The QTLs were detected on nine chromosomes, mostly clustered in concordance with the high genetic correlations between NDF and ADF. Adjustment of NDF for ADF and ADF for ADL revealed that most of the variability for CWCs in this population is in ADF. Many of the QTLs detected in this study have also been detected in other populations, and several are linked to candidate genes for cellulose or starch biosynthesis. The genetic information obtained in this study should be useful to breeding efforts aimed at improving the quality of maize silage.

Characterization of two GL8 paralogs reveals that the 3-ketoacyl reductase component of fatty acid elongase is essential for maize (Zea mays L.) development

Prior analyses established that the maize (Zea mays L.) gl8a gene encodes 3-ketoacyl reductase, a component of the fatty acid elongase required for the biosynthesis of very long chain fatty acids (VLCFAs). A paralogous gene, gl8b, has been identified that is 96% identical to gl8a. The gl8a and gl8b genes map to syntenic chromosomal regions, have similar, but not identical, expression patterns, and encode proteins that are 97% identical. Both of these genes are required for the normal accumulation of cuticular waxes on seedling leaves. The chemical composition of the cuticular waxes from gl8a and gl8b mutants indicates that these genes have at least overlapping, if not redundant, functions in cuticular wax biosynthesis. Although gl8a and gl8b double mutant kernels have endosperms that cannot be distinguished from wild-type siblings, these kernels are non-viable because their embryos fail to undergo normal development. Double mutant kernels accumulate substantially reduced levels of VLCFAs. VLCFAs are components of a variety of compounds, for example, cuticular waxes, suberin, and sphingolipids. Consistent with their essential nature in yeast, the accumulation of the ceramide moiety of sphingolipids is substantially reduced and their fatty acid composition altered in gl8a and gl8b double mutant kernels relative to wild-type kernels. Hence, we hypothesize that sphingolipids or other VLCFA-containing compounds are essential for normal embryo development.

Engrailed-ZmOCL1 fusions cause a transient reduction of kernel size in maize

ZmOCL1 is the founding member of the ZmOCL (Outer Cell Layer) family encoding putative transcription factors of the HD-ZIP IV class. It is expressed in the L1 cell layer of the embryo and several other tissues of maize. After determination of the intron/exon structure a mutator insertion was isolated in the upstream region. No notable phenotypes and wildtype levels of ZmOCL1 transcript were observed in homozygous mutant plants. In contrast transgenic plants carrying a fusion of the repressor domain of the Drosophila Engrailed gene with the DNA binding and dimerisation domains of ZmOCL1 showed a transient reduction of embryo, endosperm and kernel size that was most obvious around 15 DAP. An inverse relationship was observed between the degree of size reduction and the expression level of the transcript. In reciprocal crosses the size reduction was only observed when the transgenic plants were used as females and no expression of male transmitted transgenes was detected. Smaller kernels resembled younger kernels of wild-type siblings indicating that interference with ZmOCL1 function leads to an overall slow-down of early kernel development. Based on marker gene analysis ZmOCL1 may act via a modification of gibberellin levels. Phylogenetic analyses based on the intron/exon structure and sequence similarities of ZmOCL1 and other HD-ZIP IV proteins from maize, rice and Arabidopsis helped to identify orthologues and suggested an evolution in the function of individual genes after the divergence of monocots and dicots.

Identification and characterization of endoplasmic reticulum-associated degradation proteins differentially affected by endoplasmic reticulum stress

The disposal of misfolded proteins from the lumen of the endoplasmic reticulum (ER) is one of the quality control mechanisms present in the protein secretory pathway. Through ER-associated degradation, misfolded substrates are targeted to the cytosol where they are degraded by the proteasome. We have identified four maize (Zea mays) Der1-like genes (Zm Derlins) that encode homologs of Der1p, a yeast (Saccharomyces cerevisiae) protein implicated in ER-associated degradation. Zm Derlins are capable of functionally complementing a yeast Der1 deletion mutant. Such complementation indicates that the Der1p function is conserved among species. Zm Derlin genes are expressed at low levels throughout the plant, but appear prevalent in tissues with high activity of secretory protein accumulation, including developing endosperm cells. Expression of three of the four Zm Derlin genes increases during ER stress, with Zm Derlin1-1 showing the strongest induction. Subcellular fractionation experiments localized Zm Derlin proteins to the membrane fraction of microsomes. In maize endosperm, Zm Derlin proteins were found primarily associated with ER-derived protein bodies regardless of the presence of an ER stress response.

DNA rearrangement in orthologous orp regions of the maize, rice and sorghum genomes

The homeologous Orp1 and Orp2 regions of maize and the orthologous regions in sorghum and rice were compared by generating sequence data for >486 kb of genomic DNA. At least three genic rearrangements differentiate the maize Orp1 and Orp2 segments, including an insertion of a single gene and two deletions that removed one gene each, while no genic rearrangements were detected in the maize Orp2 region relative to sorghum. Extended comparison of the orthologous Orp regions of sorghum and japonica rice uncovered numerous genic rearrangements and the presence of a transposon-rich region in rice. Only 11 of 27 genes (40%) are arranged in the same order and orientation between sorghum and rice. Of the 8 genes that are uniquely present in the sorghum region, 4 were found to have single-copy homologs in both rice and Arabidopsis, but none of these genes are located near each other, indicating frequent gene movement. Further comparison of the Orp segments from two rice subspecies, japonica and indica, revealed that the transposon-rich region is both an ancient and current hotspot for retrotransposon accumulation and genic rearrangement. We also identify unequal gene conversion as a mechanism for maize retrotransposon rearrangement.

The inheritance and evolution of leaf pigmentation and pubescence in teosinte

To investigate the genetic mechanisms that underlie morphological evolution in natural populations, we employed QTL mapping to dissect the inheritance of leaf sheath characters that distinguish Chalco from Balsas teosinte. Abundant macrohairs (trichomes) and intense anthocyanin accumulation are found in Chalco teosinte sheaths whereas Balsas teosinte leaf sheaths are green and glabrous. These character states may represent adaptations to the cooler highland (Chalco) vs. warmer middle-elevation (Balsas) climates. QTL mapping in multiple populations revealed a mix of major- and minor-effect QTL affecting both sheath color (anthocyanin) and macrohair abundance. The major QTL for macrohairs accounts for 52% of the parental difference. Epistatic interactions were detected between the major-effect QTL and multiple other QTL for both traits, accounting for substantial portions of phenotypic variance. Developmental analyses suggest that regulatory program changes underlie the phenotypic differences. Sheath anthocyanin QTL are clearly associated with b1 and a3, both of which are regulators of anthocyanin biosynthesis. Our findings suggest that changes in a small number of QTL can lead to morphological evolution by modulating existing developmental programs.

Duplication and DNA segmental loss in the rice genome: implications for diploidization

* Large-scale duplication events have been recently uncovered in the rice genome, but different interpretations were proposed regarding the extent of the duplications. * Through analysing the 370 Mb genome sequences assembled into 12 chromosomes of Oryza sativa subspecies indica, we detected 10 duplicated blocks on all 12 chromosomes that contained 47% of the total predicted genes. Based on the phylogenetic analysis, we inferred that this was a result of a genome duplication that occurred c. 70 million years ago, supporting the polyploidy origin of the rice genome. In addition, a segmental duplication was also identified involving chromosomes 11 and 12, which occurred c. 5 million years ago. * Following the duplications, there have been large-scale chromosomal rearrangements and deletions. About 30-65% of duplicated genes were lost shortly after the duplications, leading to a rapid diploidization. * Together with other lines of evidence, we propose that polyploidization is still an ongoing process in grasses of polyploidy origins.

Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives

The cultivated Brassica species are the group of crops most closely related to Arabidopsis thaliana (Arabidopsis). They represent models for the application in crops of genomic information gained in Arabidopsis and provide an opportunity for the investigation of polyploid genome formation and evolution. The scientific literature contains contradictory evidence for the dynamics of the evolution of polyploid genomes. We aimed at overcoming the inherent complexity of Brassica genomes and clarify the effects of polyploidy on the evolution of genome microstructure in specific segments of the genome. To do this, we have constructed bacterial artificial chromosome (BAC) libraries from genomic DNA of B. rapa subspecies trilocularis (JBr) and B. napus var Tapidor (JBnB) to supplement an existing BAC library from B. oleracea. These allowed us to analyse both recent polyploidization (under 10,000 years in B. napus) and more ancient polyploidization events (ca. 20 Myr for B. rapa and B. oleracea relative to Arabidopsis), with an analysis of the events occurring on an intermediate time scale (over the ca. 4 Myr since the divergence of the B. rapa and B. oleracea lineages). Using the Arabidopsis genome sequence and clones from the JBr library, we have analysed aspects of gene conservation and microsynteny between six regions of the genome of B. rapa with the homoeologous regions of the genomes of B. oleracea and Arabidopsis. Extensive divergence of gene content was observed between the B. rapa paralogous segments and their homoeologous segments within the genome of Arabidopsis. A pattern of interspersed gene loss was identified that is similar, but not identical, to that observed in B. oleracea. The conserved genes show highly conserved collinearity with their orthologues across genomes, but a small number of species-specific rearrangements were identified. Thus the evolution of genome microstructure is an ongoing process. Brassica napus is a recently formed polyploid resulting from the hybridization of B. rapa (containing the Brassica A genome) and B. oleracea (containing the Brassica C genome). Using clones from the JBnB library, we have analysed the microstructure of the corresponding segments of the B. napus genome. The results show that there has been little or no change to the microstructure of the analysed segments of the Brassica A and C genomes as a consequence of the hybridization event forming natural B. napus. The observations indicate that, upon polyploid formation, these segments of the genome did not undergo a burst of evolution discernible at the scale of microstructure.

Characterization of the maize endosperm transcriptome and its comparison to the rice genome

The cereal endosperm is a major organ of the seed and an important component of the world's food supply. To understand the development and physiology of the endosperm of cereal seeds, we focused on the identification of genes expressed at various times during maize endosperm development. We constructed several cDNA libraries to identify full-length clones and subjected them to a twofold enrichment. A total of 23,348 high-quality sequence-reads from 5'- and 3'-ends of cDNAs were generated and assembled into a unigene set representing 5326 genes with paired sequence-reads. Additional sequencing yielded a total of 3160 (59%) completely sequenced, full-length cDNAs. From 5326 unigenes, 4139 (78%) can be aligned with 5367 predicted rice genes and by taking only the "best hit" be mapped to 3108 positions on the rice genome. The 22% unigenes not present in rice indicate a rapid change of gene content between rice and maize in only 50 million years. Differences in rice and maize gene numbers also suggest that maize has lost a large number of duplicated genes following tetraploidization. The larger number of gene copies in rice suggests that as many as 30% of its genes arose from gene amplification, which would extrapolate to a significant proportion of the estimated 44,027 candidate genes of its entire genome. Functional classification of the maize endosperm unigene set indicated that more than a fourth of the novel functionally assignable genes found in this study are involved in carbohydrate metabolism, consistent with its role as a storage organ.