Patterns of chromosomal duplication in maize and their implications for comparative maps of the grasses

The biosynthesis of storage reserves and auxin is coordinated by a hierarchical regulatory network in maize endosperm

Grain filling in maize (Zea mays) is intricately linked to cell development, involving the regulation of genes responsible for the biosynthesis of storage reserves (starch, proteins, and lipids) and phytohormones. However, the regulatory network coordinating these biological functions remains unclear. In this study, we identified 1744 high-confidence target genes co-regulated by the transcription factors (TFs) ZmNAC128 and ZmNAC130 (ZmNAC128/130) through chromatin immunoprecipitation sequencing coupled with RNA-seq analysis in the zmnac128/130 loss-of-function mutants. We further constructed a hierarchical regulatory network using DNA affinity purification sequencing analysis of downstream TFs regulated by ZmNAC128/130. In addition to target genes involved in the biosynthesis of starch and zeins, we discovered novel target genes of ZmNAC128/130 involved in the biosynthesis of lipids and indole-3-acetic acid (IAA). Consistently, the number of oil bodies, as well as the contents of triacylglycerol, and IAA were significantly reduced in zmnac128/130. The hierarchical regulatory network centered by ZmNAC128/130 revealed a significant overlap between the direct target genes of ZmNAC128/130 and their downstream TFs, particularly in regulating the biosynthesis of storage reserves and IAA. Our results indicated that the biosynthesis of storage reserves and IAA is coordinated by a multi-TFs hierarchical regulatory network in maize endosperm.

The WIP6 transcription factor TOO MANY LATERALS specifies vein type in C4 and C3 grass leaves

Grass leaves are invariantly strap shaped with an elongated distal blade and a proximal sheath that wraps around the stem. Underpinning this shape is a scaffold of leaf veins, most of which extend in parallel along the proximo-distal leaf axis. Differences between species are apparent both in the vein types that develop and in the distance between veins across the medio-lateral leaf axis. A prominent engineering goal is to increase vein density in leaves of C3 photosynthesizing species to facilitate the introduction of the more efficient C4 pathway. Here, we discover that the WIP6 transcription factor TOO MANY LATERALS (TML) specifies vein rank in both maize (C4) and rice (C3). Loss-of-function tml mutations cause large lateral veins to develop in positions normally occupied by smaller intermediate veins, and TML transcript localization in wild-type leaves is consistent with a role in suppressing lateral vein development in procambial cells that form intermediate veins. Attempts to manipulate TML function in rice were unsuccessful because transgene expression was silenced, suggesting that precise TML expression is essential for shoot viability. This finding may reflect the need to prevent the inappropriate activation of downstream targets or, given that transcriptome analysis revealed altered cytokinin and auxin signaling profiles in maize tml mutants, the need to prevent local or general hormonal imbalances. Importantly, rice tml mutants display an increased occupancy of veins in the leaf, providing a step toward an anatomical chassis for C4 engineering. Collectively, a conserved mechanism of vein rank specification in grass leaves has been revealed.

Efficient Algorithms for Genomic Duplication Models

An important issue in evolutionary molecular biology is to discover genomic duplication episodes and their correspondence to the species tree. Existing approaches vary in the two fundamental aspects: the choice of evolutionary scenarios that model allowed locations of duplications in the species tree, and the rules of clustering gene duplications from gene trees into a single multiple duplication event. Here we study the method of clustering called minimum episodes for several models of allowed evolutionary scenarios with a focus on interval models in which every gene duplication has an interval consisting of allowed locations in the species tree. We present mathematical foundations for general genomic duplication problems. Next, we propose the first linear time and space algorithm for minimum episodes clustering jointly for any interval model and the algorithm for the most general model in which every evolutionary scenario is allowed. We also present a comparative study of different models of genomic duplication based on simulated and empirical datasets. We provided algorithms and tools that could be applied to solve efficiently minimum episodes clustering problems. Our comparative study helps to identify which model is the most reasonable choice in inferring genomic duplication events.

Dynamic gene amplification and function diversification of grass-specific O-methyltransferase gene family

The plant O-methyltransferases are dependent on S-Adenosyl-l-methionine, which can catalyze a variety of secondary metabolites. Here we identified different number of OMT genes from the respective grass genomes. Phylogenetic analysis showed that this OMT gene family is a grass-specific gene family that is different from COMT. Most of genes were expanded by tandem and segment duplication after the species split from their progenitor. Furthermore, genes from Group I and two clusters from group II are only present in Panicoideae, which included Bx10 and Bx7 involved in the benzoxazinoids pathway, suggesting these genes could participate in insect resistance in Panicoideae. Gene expression profiles showed that OMT genes were preferentially expressed in vegetative stages, especially in roots. These results revealed that this grass-specific OMT gene family could affect the development of vegetative stages, and be involved in the benzoxazinoids pathway or suberin biosynthesis that was different from COMT.

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.

Genome Halving with an Outgroup

Some genomes are known to have incurred a genome doubling (tetraploidization) event in their evolutionary history, and this is reflected today in patterns of duplicated segments scattered throughout their chromosomes. These duplications may be used as data to “halve” the genome, i.e. to reconstruct the ancestral genome at the moment of tetraploidization, but the solution is often highly non-unique. To resolve this problem, we adapt the genome halving algorithm of El-Mabrouk and Sankoff to take account of an external reference genome. We apply this to reconstruct the tetraploid ancestor of maize, using either rice or sorghum as the reference.

Evolutionary divergence of phytochrome protein function in Zea mays PIF3 signaling

Two maize phytochrome-interacting factor (PIF) basic helix-loop-helix (bHLH) family members, ZmPIF3.1 and ZmPIF3.2, were identified, cloned and expressed in vitro to investigate light-signaling interactions. A phylogenetic analysis of sequences of the maize bHLH transcription factor gene family revealed the extent of the PIF family, and a total of seven predicted PIF-encoding genes were identified from genes encoding bHLH family VIIa/b proteins in the maize genome. To investigate the role of maize PIFs in phytochrome signaling, full-length cDNAs for phytochromes PhyA2, PhyB1, PhyB2 and PhyC1 from maize were cloned and expressed in vitro as chromophorylated holophytochromes. We showed that ZmPIF3.1 and ZmPIF3.2 interact specifically with the Pfr form of maize holophytochrome B1 (ZmphyB1), showing no detectable affinity for the Pr form. Maize holophytochrome B2 (ZmphyB2) showed no detectable binding affinity for PIFs in either Pr or Pfr forms, but phyB Pfr from Arabidopsis interacted with ZmPIF3.1 similarly to ZmphyB1 Pfr. We conclude that subfunctionalization at the protein-protein interaction level has altered the role of phyB2 relative to that of phyB1 in maize. Since the phyB2 mutant shows photomorphogenic defects, we conclude that maize phyB2 is an active photoreceptor, without the binding of PIF3 seen in other phyB family proteins.

pATsi: Paralogs and Singleton Genes from Arabidopsis thaliana

Arabidopsis thaliana is widely accepted as a model species in plant biology. Its genome, due to its small size and diploidy, was the first to be sequenced among plants, making this species also a reference for plant comparative genomics. Nevertheless, the evolutionary mechanisms that shaped the Arabidopsis genome are still controversial. Indeed, duplications, translocations, inversions, and gene loss events that contributed to the current organization are difficult to be traced. A reliable identification of paralogs and single-copy genes is essential to understand these mechanisms. Therefore, we implemented a dedicated pipeline to identify paralog genes and classify single-copy genes into opportune categories. PATsi, a web-accessible database, was organized to allow the straightforward access to the paralogs organized into networks and to the classification of single-copy genes. This permits to efficiently explore the gene collection of Arabidopsis for evolutionary investigations and comparative genomics.

High level of microsynteny and purifying selection affect the evolution of WRKY family in Gramineae

The WRKY gene family, which encodes proteins in the regulation processes of diverse developmental stages, is one of the largest families of transcription factors in higher plants. In this study, by searching for interspecies gene colinearity (microsynteny) and dating the age distributions of duplicated genes, we found 35 chromosomal segments of subgroup I genes of WRKY family (WRKY I) in four Gramineae species (Brachypodium, rice, sorghum, and maize) formed eight orthologous groups. After a stepwise gene-by-gene reciprocal comparison of all the protein sequences in the WRKY I gene flanking areas, highly conserved regions of microsynteny were found in the four Gramineae species. Most gene pairs showed conserved orientation within syntenic genome regions. Furthermore, tandem duplication events played the leading role in gene expansion. Eventually, environmental selection pressure analysis indicated strong purifying selection for the WRKY I genes in Gramineae, which may have been followed by gene loss and rearrangement. The results presented in this study provide basic information of Gramineae WRKY I genes and form the foundation for future functional studies of these genes. High level of microsynteny in the four grass species provides further evidence that a large-scale genome duplication event predated speciation.

Development of two major resources for pea genomics: the GenoPea 13.2K SNP Array and a high-density, high-resolution consensus genetic map

Single nucleotide polymorphism (SNP) arrays represent important genotyping tools for innovative strategies in both basic research and applied breeding. Pea is an important food, feed and sustainable crop with a large (about 4.45 Gbp) but not yet available genome sequence. In the present study, 12 pea recombinant inbred line populations were genotyped using the newly developed GenoPea 13.2K SNP Array. Individual and consensus genetic maps were built providing insights into the structure and organization of the pea genome. Largely collinear genetic maps of 3918-8503 SNPs were obtained from all mapping populations, and only two of these exhibited putative chromosomal rearrangement signatures. Similar distortion patterns in different populations were noted. A total of 12 802 transcript-derived SNP markers placed on a 15 079-marker high-density, high-resolution consensus map allowed the identification of ohnologue-rich regions within the pea genome and the localization of local duplicates. Dense syntenic networks with sequenced legume genomes were further established, paving the way for the identification of the molecular bases of important agronomic traits segregating in the mapping populations. The information gained on the structure and organization of the genome from this research will undoubtedly contribute to the understanding of the evolution of the pea genome and to its assembly. The GenoPea 13.2K SNP Array and individual and consensus genetic maps are valuable genomic tools for plant scientists to strengthen pea as a model for genetics and physiology and enhance breeding.

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.

Biochemical characterization of the maize cytokinin dehydrogenase family and cytokinin profiling in developing maize plantlets in relation to the expression of cytokinin dehydrogenase genes

The cytokinin dehydrogenases (CKX; EC 1.5.99.12) are a protein family that maintains the endogenous levels of cytokinins in plants by catalyzing their oxidative degradation. The CKX family in maize (Zea mays L.) has thirteen members, only two of which--ZmCKX1 and ZmCKX10--have previously been characterized in detail. In this study, nine further maize CKX isoforms were heterologously expressed in Escherichia coli, purified by affinity and ion-exchange chromatography and biochemically characterized. ZmCKX6 and ZmCKX9 could only be expressed successfully after the removal of putative sequence-specific vacuolar sorting signals (LLPT and LPTS, respectively), suggesting that these proteins are localized to the vacuole. Substrate specificity analyses revealed that the CKX isoforms can be grouped into two subfamilies: members of the first strongly prefer cytokinin free bases while members of the second degrade a broad range of substrates. The most active isoform was found to be ZmCKX1. One of the studied isoforms, ZmCKX6, seemed to encode a nonfunctional enzyme due to a mutation in a conserved HFG protein domain at the C-terminus. Site-directed mutagenesis experiments revealed that this domain is essential for CKX activity. The roles of the maize CKX enzymes in the development of maize seedlings during the two weeks immediately after radicle emergence were also investigated. It appears that ZmCKX1 is a key regulator of active cytokinin levels in developing maize roots. However, the expression of individual CKX isoforms in the shoots varied and none of them seemed to have strong effects on the cytokinin pool.

Rapid genomic and transcriptomic alterations induced by wide hybridization: Chrysanthemum nankingense × Tanacetum vulgare and C. crassum × Crossostephium chinense (Asteraceae)

Background

Hybridization is a major driver of evolution in plants. In a number of plant species, the process of hybridization has been revealed to be accompanied by wide-ranging genetic and epigenetic alterations, some of which have consequences on gene transcripts. The Asteraceae family includes a number of polyploid species, and wide crossing is seen as a viable means of genetically improving ornamental species such as Chrysanthemum spp. However, the consequences of hybridization in this taxon have yet to be characterized.

Results

Amplified fragment length polymorphism (AFLP), methylation sensitive amplification polymorphism (MSAP) and cDNA-AFLP profiling of the two intergeneric hybrids C. nankingense × Tanacetum vulgare and C. crassum × Crossostephium chinense were employed to characterize, respectively, the genomic, epigenomic and transcriptomic changes induced by the hybridization event. The hybrids’ AFLP profiles included both the loss of specific parental fragments and the gain of fragments not present in either parent’s profile. About 10% of the paternal fragments were not inherited by the hybrid, while the corresponding rate for the maternal parent fragments was around 4–5%. The novel fragments detected may have arisen either due to heterozygosity in one or other parent, or as a result of a deletion event following the hybridization. Around one half of the cDNA-AFLP fragments were common to both parents, about 30% were specific to the female parent, and somewhat under 20% specific to the male parent; the remainder (2.9-4.7%) of the hybrids’ fragments were not present in either parent’s profile. The MSAP fingerprinting demonstrated that the hybridization event also reduced the amount of global cytosine methylation, since > 50% of the parental fragments were methylated, while the corresponding frequencies for the two hybrids were 48.5% and 50.4%.

Conclusions

Combining two different Asteraceae genomes via hybridization clearly induced a range of genomic and epigenomic alterations, some of which had an effect on the transcriptome. The rapid genomic and transcriptomic alterations induced by hybridization may accelerate the evolutionary process among progenies.

The genomics of incompatibility factors and sex determination in hybridizing species of Cottus (Pisces)

Cottus rhenanus and Cottus perifretum have formed hybrid lineages and narrow hybrid zones that can be best explained through the action of natural selection. However, the underlying selective forces as well as their genomic targets are not well understood. This study identifies genomic regions in the parental species that cause hybrid incompatibilities and tests whether these manifest in a sex-specific manner to learn about processes that affect natural hybridization in Cottus. Interspecific F2 crosses were analyzed for 255 markers for genetic mapping and to detect transmission distortion as a sign for genetic incompatibilities. The Cottus map consists of 24 linkage groups with a total length of 1575.4 cM. A male heterogametic (XY) sex determination region was found on different linkage groups in the two parental species. Genetic incompatibilities were incomplete, varied among individuals and populations and were not associated with the heterogametic sex. The variance between populations and individuals makes it unlikely that there are species-specific incompatibility loci that could affect the gene pool of natural hybrids in a simple and predictable way. Conserved synteny with sequenced fish genomes permits to genetically study the Cottus genome through the transfer of genomic information from the model fish species. Homology relationships of candidate genomic regions in Cottus indicate that sex determination is not based on the same genomic regions found in other fish species. This suggests a fast evolutionary turnover of the genetic basis of sex determination that, together with the small size of the heterogametic regions, may contribute to the absence of fitness effects related to the Haldane's rule.

Generation of tandem direct duplications by reversed-ends transposition of maize ac elements

Tandem direct duplications are a common feature of the genomes of eukaryotes ranging from yeast to human, where they comprise a significant fraction of copy number variations. The prevailing model for the formation of tandem direct duplications is non-allelic homologous recombination (NAHR). Here we report the isolation of a series of duplications and reciprocal deletions isolated de novo from a maize allele containing two Class II Ac/Ds transposons. The duplication/deletion structures suggest that they were generated by alternative transposition reactions involving the termini of two nearby transposable elements. The deletion/duplication breakpoint junctions contain 8 bp target site duplications characteristic of Ac/Ds transposition events, confirming their formation directly by an alternative transposition mechanism. Tandem direct duplications and reciprocal deletions were generated at a relatively high frequency (~0.5 to 1%) in the materials examined here in which transposons are positioned nearby each other in appropriate orientation; frequencies would likely be much lower in other genotypes. To test whether this mechanism may have contributed to maize genome evolution, we analyzed sequences flanking Ac/Ds and other hAT family transposons and identified three small tandem direct duplications with the structural features predicted by the alternative transposition mechanism. Together these results show that some class II transposons are capable of directly inducing tandem sequence duplications, and that this activity has contributed to the evolution of the maize genome.

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.

Identification, characterization and interpretation of single-nucleotide sequence variation in allopolyploid crop species

An understanding of nature and extent of nucleotide sequence variation is required for programmes of discovery and characterization of single nucleotide polymorphisms (SNPs), which provide the most versatile class of molecular genetic marker. A majority of higher plant species are polyploids, and allopolyploidy, because of hybrid formation between closely related taxa, is very common. Mutational variation may arise both between allelic (homologous) sequences within individual subgenomes and between homoeologous sequences among subgenomes, in addition to paralogous variation between duplicated gene copies. Successful SNP validation in allopolyploids depends on differentiation of the sequence variation classes. A number of biological factors influence the feasibility of discrimination, including degree of gene family complexity, inbreeding or outbreeding reproductive habit, and the level of knowledge concerning progenitor diploid species. In addition, developments in high-throughput DNA sequencing and associated computational analysis provide general solutions for the genetic analysis of allopolyploids. These issues are explored in the context of experience from a range of allopolyploid species, representing grain (wheat and canola), forage (pasture legumes and grasses), and horticultural (strawberry) crop. Following SNP discovery, detection in routine genotyping applications also presents challenges for allopolyploids. Strategies based on either design of subgenome-specific SNP assays through homoeolocus-targeted polymerase chain reaction (PCR) amplification, or detection of incremental changes in nucleotide variant dosage, are described.

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.

QTL involved in the partial restoration of male fertility of C-type cytoplasmic male sterility in maize

Partial restoration of male fertility limits the use of C-type cytoplasmic male sterility (C-CMS) for the production of hybrid seeds in maize. Nevertheless, the genetic basis of the trait is still unknown. Therefore, the aim to this study was to identify genomic regions that govern partial restoration by means of a QTL analysis carried out in an F(2) population (n = 180). This population was derived from the Corn Belt inbred lines B37C and K55. F(2)BC(1) progenies were phenotyped at three locations in Switzerland. Male fertility was rated according to the quality and number of anthers as well as the anthesis-silking interval. A weak effect of environment on the expression of partial restoration was reflected by high heritabilities of all fertility-related traits. Partial restoration was inherited like an oligogenic trait. Three major QTL regions were found consistently across environments in the chromosomal bins 2.09, 3.06 and 7.03. Therefore, a marker-assisted counter-selection of partial restoration is promising. Minor QTL regions were found on chromosomes 3, 4, 5, 6 and 8. A combination of partial restorer alleles at different QTL can lead to full restoration of fertility. The maternal parent was clearly involved in the partial restoration, because the restorer alleles at QTL in bins 2.09, 6.04 and 7.03 originated from B37. The three major QTL regions collocated with other restorer genes of maize, a phenomenon, which seems to be typical for restorer genes. Therefore, a study of the clusters of restorer genes in maize could lead to a better understanding of their evolution and function. In this respect, the long arm of chromosome 2 is particularly interesting, because it harbors restorer genes for the three major CMS systems (C, T and S) of maize.

Polycomb-group mediated epigenetic mechanisms through plant evolution

Polycomb Group (PcG) proteins form an epigenetic "memory system", conserved in both plants and animals, controlling global gene expression during development via histone modifications. The role of PcG proteins in plants was primarily explored in Arabidopsis thaliana, where PcG regulation of developmental processes was demonstrated throughout the plant life cycle. Our knowledge about the PcG machinery in terrestrial plants other than Arabidopsis began to accumulate only in recent years. In this review we summarize recent emerging data on the evolution and diversification of PcG mechanisms in various phyla, from early-diverging plants, including members of the Chlorophyte algae, through bryophytes and flowering plants. We describe the compositions of the PcG gene families, their so-far studied expression profiles, and finally summarize commonalities vs. differences among PcG functions across the various species. This article is part of a Special Issue entitled: Epigenetic control of cellular and developmental processes in plants.

Factors that contribute to variation in evolutionary rate among Arabidopsis genes

Surprisingly, few studies have described evolutionary rate variation among plant nuclear genes, with little investigation of the causes of rate variation. Here, we describe evolutionary rates for 11,492 ortholog pairs between Arabidopsis thaliana and A. lyrata and investigate possible contributors to rate variation among these genes. Rates of evolution at synonymous sites vary along chromosomes, suggesting that mutation rates vary on genomic scales, perhaps as a function of recombination rate. Rates of evolution at nonsynonymous sites correlate most strongly with expression patterns, but they also vary as to whether a gene is duplicated and retained after a whole-genome duplication (WGD) event. WGD genes evolve more slowly, on average, than nonduplicated genes and non-WGD duplicates. We hypothesize that levels and patterns of expression are not only the major determinants that explain nonsynonymous rate variation among genes but also a critical determinant of gene retention after duplication.

The selection and use of sorghum (Sorghum propinquum) bacterial artificial chromosomes as cytogenetic FISH probes for maize (Zea mays L.)

The integration of genetic and physical maps of maize is progressing rapidly, but the cytogenetic maps lag behind, with the exception of the pachytene fluorescence in situ hybridization (FISH) maps of maize chromosome 9. We sought to produce integrated FISH maps of other maize chromosomes using Core Bin Marker loci. Because these 1 Kb restriction fragment length polymorphism (RFLP) probes are below the FISH detection limit, we used BACs from sorghum, a small-genome relative of maize, as surrogate clones for FISH mapping. We sequenced 151 maize RFLP probes and compared in silico BAC selection methods to that of library filter hybridization and found the latter to be the best. BAC library screening, clone verification, and single-clone selection criteria are presented along with an example of transgenomic BAC FISH mapping. This strategy has been used to facilitate the integration of RFLP and FISH maps in other large-genome species.

Structure and expression of the maize (Zea mays L.) SUN-domain protein gene family: evidence for the existence of two divergent classes of SUN proteins in plants

BACKGROUND

The nuclear envelope that separates the contents of the nucleus from the cytoplasm provides a surface for chromatin attachment and organization of the cortical nucleoplasm. Proteins associated with it have been well characterized in many eukaryotes but not in plants. SUN (Sad1p/Unc-84) domain proteins reside in the inner nuclear membrane and function with other proteins to form a physical link between the nucleoskeleton and the cytoskeleton. These bridges transfer forces across the nuclear envelope and are increasingly recognized to play roles in nuclear positioning, nuclear migration, cell cycle-dependent breakdown and reformation of the nuclear envelope, telomere-led nuclear reorganization during meiosis, and karyogamy.

RESULTS

We found and characterized a family of maize SUN-domain proteins, starting with a screen of maize genomic sequence data. We characterized five different maize ZmSUN genes (ZmSUN1-5), which fell into two classes (probably of ancient origin, as they are also found in other monocots, eudicots, and even mosses). The first (ZmSUN1, 2), here designated canonical C-terminal SUN-domain (CCSD), includes structural homologs of the animal and fungal SUN-domain protein genes. The second (ZmSUN3, 4, 5), here designated plant-prevalent mid-SUN 3 transmembrane (PM3), includes a novel but conserved structural variant SUN-domain protein gene class. Mircroarray-based expression analyses revealed an intriguing pollen-preferred expression for ZmSUN5 mRNA but low-level expression (50-200 parts per ten million) in multiple tissues for all the others. Cloning and characterization of a full-length cDNA for a PM3-type maize gene, ZmSUN4, is described. Peptide antibodies to ZmSUN3, 4 were used in western-blot and cell-staining assays to show that they are expressed and show concentrated staining at the nuclear periphery.

CONCLUSIONS

The maize genome encodes and expresses at least five different SUN-domain proteins, of which the PM3 subfamily may represent a novel class of proteins with possible new and intriguing roles within the plant nuclear envelope. Expression levels for ZmSUN1-4 are consistent with basic cellular functions, whereas ZmSUN5 expression levels indicate a role in pollen. Models for possible topological arrangements of the CCSD-type and PM3-type SUN-domain proteins are presented.

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.

Comparative analysis of the Brassica oleracea genetic map and the Arabidopsis thaliana genome

We further investigated genome macrosynteny for Brassica species and Arabidopsis thaliana. This work aimed at comparative map construction for B. oleracea and A. thaliana chromosomes based on 160 known A. thaliana probes: 147 expressed sequence tags (ESTs) and 13 full-length cDNA clones. Based on an in silico study of the A. thaliana genome, most of the selected ESTs (83%) represented unique or low-copy genes. We identified conserved segments by the visual inspection of comparative data with a priori assumptions, and established their significance with the LineUp algorithm. Evaluation of the number of B. oleracea gene copies per A. thaliana EST revealed a fixed upward trend. We established a segregation distortion pattern for all genetic loci, with particular consideration of the type of selection (gametic or zygotic), and discuss its possible impact on genetic map construction. Consistent with previous reports, we found evidence for numerous chromosome rearrangements and the genome fragment replication of B. oleracea that have taken place since the divergence of the two species. Also, we found that over 54% of the B. oleracea genome is covered by 24 segments conserved with the A. thaliana genome. The average conserved segment is composed of 5 loci covering 19.3 cM in the B. oleracea genetic map and 2.42 Mb in the A. thaliana physical map. We have also attempted to use a unified system of conserved blocks (previously described) to verify our results and perform a comprehensive comparison with other Brassica species.

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.

Construction of a high-resolution genetic linkage map and comparative genome analysis for the reef-building coral Acropora millepora

A high-resolution genetic linkage map for the coral Acropora millepora is constructed and compared with other metazoan genomes, revealing syntenic blocks.

Background

Worldwide, coral reefs are in decline due to a range of anthropogenic disturbances, and are now also under threat from global climate change. Virtually nothing is currently known about the genetic factors that might determine whether corals adapt to the changing climate or continue to decline. Quantitative genetics studies aiming to identify the adaptively important genomic loci will require a high-resolution genetic linkage map. The phylogenetic position of corals also suggests important applications for a coral genetic map in studies of ancestral metazoan genome architecture.

Results

We constructed a high-resolution genetic linkage map for the reef-building coral Acropora millepora, the first genetic map reported for any coral, or any non-Bilaterian animal. More than 500 single nucleotide polymorphism (SNP) markers were developed, most of which are transferable in populations from Orpheus Island and Great Keppel Island. The map contains 429 markers (393 gene-based SNPs and 36 microsatellites) distributed in 14 linkage groups, and spans 1,493 cM with an average marker interval of 3.4 cM. Sex differences in recombination were observed in a few linkage groups, which may be caused by haploid selection. Comparison of the coral map with other metazoan genomes (human, nematode, fly, anemone and placozoan) revealed synteny regions.

Conclusions

Our study develops a framework that will be essential for future studies of adaptation in coral and it also provides an important resource for future genome sequence assembly and for comparative genomics studies on the evolution of metazoan genome structure.

Detailed analysis of a contiguous 22-Mb region of the maize genome

Most of our understanding of plant genome structure and evolution has come from the careful annotation of small (e.g., 100 kb) sequenced genomic regions or from automated annotation of complete genome sequences. Here, we sequenced and carefully annotated a contiguous 22 Mb region of maize chromosome 4 using an improved pseudomolecule for annotation. The sequence segment was comprehensively ordered, oriented, and confirmed using the maize optical map. Nearly 84% of the sequence is composed of transposable elements (TEs) that are mostly nested within each other, of which most families are low-copy. We identified 544 gene models using multiple levels of evidence, as well as five miRNA genes. Gene fragments, many captured by TEs, are prevalent within this region. Elimination of gene redundancy from a tetraploid maize ancestor that originated a few million years ago is responsible in this region for most disruptions of synteny with sorghum and rice. Consistent with other sub-genomic analyses in maize, small RNA mapping showed that many small RNAs match TEs and that most TEs match small RNAs. These results, performed on approximately 1% of the maize genome, demonstrate the feasibility of refining the B73 RefGen_v1 genome assembly by incorporating optical map, high-resolution genetic map, and comparative genomic data sets. Such improvements, along with those of gene and repeat annotation, will serve to promote future functional genomic and phylogenomic research in maize and other grasses.

Characterization of new maize genes putatively involved in cytokinin metabolism and their expression during osmotic stress in relation to cytokinin levels

Plant hormones, cytokinins (CKs), have been for a long time considered to be involved in plant responses to stress. However, their exact roles in processes linked to stress signalization and acclimatization to adverse environmental conditions are unknown. In this study, expression profiles of the entire gene families of CK biosynthetic and degradation genes in maize (Zea mays) during development and stress responses are described. Transcript abundance of particular genes is discussed in relation to the levels of different CK metabolites. Salt and osmotic stresses induce expression of some CK biosynthetic genes in seedlings of maize, leading to a moderate increase of active forms of CKs lasting several days during acclimatization to stress. A direct effect of CKs to mediate activation of stress responses does not seem to be possible due to the slow changes in metabolite levels. However, expression of genes involved in cytokinin signal transduction is uniformly down-regulated within 0.5 h of stress induction by an unknown mechanism. cis-Zeatin and its derivatives were found to be the most abundant CKs in young maize seedlings. We demonstrate that levels of this zeatin isomer are significantly enhanced during early stress response and that it originates independently from de novo biosynthesis in stressed tissues, possibly by elevated specific RNA degradation. By enhancing their CK levels, plants could perhaps undergo a reduction of growth rates maintained by abscisic acid accumulation in stressed tissues. A second role for cytokinin receptors in sensing turgor response is hypothesized besides their documented function in CK signaling.

Understanding meiosis and the implications for crop improvement

Over the past 50 years, the understanding of meiosis has aged like a fine bottle of wine: the complexity is developing but the wine itself is still young. While emphasis in the plant kingdom has been placed on the model diploids Arabidopsis (Arabidopsis thaliana L.) and rice (Orzya sativa L.), our research has mainly focussed on the polyploid, bread wheat (Triticum aestivum L.). Bread wheat is an important food source for nearly two-thirds of the world's population. While creating new varieties can be achieved using existing or advanced breeding lines, we would also like to introduce beneficial traits from wild related species. However, expanding the use of non-adapted and wild germplasm in cereal breeding programs will depend on the ability to manipulate the cellular process of meiosis. Three important and tightly-regulated events that occur during early meiosis are chromosome pairing, synapsis and recombination. Which key genes control these events in meiosis (and how they do so) remains to be completely answered, particularly in crops such as wheat. Although the majority of published findings are from model organisms including yeast (Saccharomyces cerevisiae) and the nematode Caenorhabditis elegans, information from the plant kingdom has continued to grow in the past decade at a steady rate. It is with this new knowledge that we ask how meiosis will contribute to the future of cereal breeding. Indeed, how has it already shaped cereal breeding as we know it today?

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.

tasselseed1 is a lipoxygenase affecting jasmonic acid signaling in sex determination of maize

Sex determination in maize is controlled by a developmental cascade leading to the formation of unisexual florets derived from an initially bisexual floral meristem. Abortion of pistil primordia in staminate florets is controlled by a tasselseed-mediated cell death process. We positionally cloned and characterized the function of the sex determination gene tasselseed1 (ts1). The TS1 protein encodes a plastid-targeted lipoxygenase with predicted 13-lipoxygenase specificity, which suggests that TS1 may be involved in the biosynthesis of the plant hormone jasmonic acid. In the absence of a functional ts1 gene, lipoxygenase activity was missing and endogenous jasmonic acid concentrations were reduced in developing inflorescences. Application of jasmonic acid to developing inflorescences rescued stamen development in mutant ts1 and ts2 inflorescences, revealing a role for jasmonic acid in male flower development in maize.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

BACKGROUND

Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers.

RESULTS

The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci (1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions.

CONCLUSION

The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

Involvement of the MADS-box gene ZMM4 in floral induction and inflorescence development in maize

The switch from vegetative to reproductive growth is marked by the termination of vegetative development and the adoption of floral identity by the shoot apical meristem (SAM). This process is called the floral transition. To elucidate the molecular determinants involved in this process, we performed genome-wide RNA expression profiling on maize (Zea mays) shoot apices at vegetative and early reproductive stages using massively parallel signature sequencing technology. Profiling revealed significant up-regulation of two maize MADS-box (ZMM) genes, ZMM4 and ZMM15, after the floral transition. ZMM4 and ZMM15 map to duplicated regions on chromosomes 1 and 5 and are linked to neighboring MADS-box genes ZMM24 and ZMM31, respectively. This gene order is syntenic with the vernalization1 locus responsible for floral induction in winter wheat (Triticum monococcum) and similar loci in other cereals. Analyses of temporal and spatial expression patterns indicated that the duplicated pairs ZMM4-ZMM24 and ZMM15-ZMM31 are coordinately activated after the floral transition in early developing inflorescences. More detailed analyses revealed ZMM4 expression initiates in leaf primordia of vegetative shoot apices and later increases within elongating meristems acquiring inflorescence identity. Expression analysis in late flowering mutants positioned all four genes downstream of the floral activators indeterminate1 (id1) and delayed flowering1 (dlf1). Overexpression of ZMM4 leads to early flowering in transgenic maize and suppresses the late flowering phenotype of both the id1 and dlf1 mutations. Our results suggest ZMM4 may play roles in both floral induction and inflorescence development.

Physical and genetic structure of the maize genome reflects its complex evolutionary history

Maize (Zea mays L.) is one of the most important cereal crops and a model for the study of genetics, evolution, and domestication. To better understand maize genome organization and to build a framework for genome sequencing, we constructed a sequence-ready fingerprinted contig-based physical map that covers 93.5% of the genome, of which 86.1% is aligned to the genetic map. The fingerprinted contig map contains 25,908 genic markers that enabled us to align nearly 73% of the anchored maize genome to the rice genome. The distribution pattern of expressed sequence tags correlates to that of recombination. In collinear regions, 1 kb in rice corresponds to an average of 3.2 kb in maize, yet maize has a 6-fold genome size expansion. This can be explained by the fact that most rice regions correspond to two regions in maize as a result of its recent polyploid origin. Inversions account for the majority of chromosome structural variations during subsequent maize diploidization. We also find clear evidence of ancient genome duplication predating the divergence of the progenitors of maize and rice. Reconstructing the paleoethnobotany of the maize genome indicates that the progenitors of modern maize contained ten chromosomes.

As a cash crop and a model biological system, maize is of great public interest. To facilitate maize molecular breeding and its basic biology research, we built a high-resolution physical map with two different fingerprinting methods on the same set of bacterial artificial chromosome clones. The physical map was integrated to a high-density genetic map and further serves as a framework for the maize genome-sequencing project. Comparative genomics showed that the euchromatic regions between rice and maize are very conserved. Physically we delimited these conserved regions and thus detected many genome rearrangements. We defined extensively the duplication blocks within the maize genome. These blocks allowed us to reconstruct the chromosomes of the maize progenitor. We detected that maize genome has experienced two rounds of genome duplications, an ancient one before maize–rice divergence and a recent one after tetraploidization.

Allopolyploidization-accommodated genomic sequence changes in triticale

BACKGROUND

Allopolyploidization is one of the major evolutionary modes of plant speciation. Recent interest in studying allopolyploids has provided significant novel insights into the mechanisms of allopolyploid formation. Compelling evidence indicates that genetic and/or epigenetic changes have played significant roles in shaping allopolyploids, but rates and modes of the changes can be very different among various species. Triticale (x Triticosecale) is an artificial species that has been used to study the evolutionary course of complex allopolyploids due to its recent origin and availability of a highly diversified germplasm pool. Scope This review summarizes recent genomics studies implemented in hexaploid and octoploid triticales and discusses the mechanisms of the changes and compares the major differences between genomic changes in triticale and other allopolyploid species.

CONCLUSIONS

Molecular studies have indicated extensive non-additive sequence changes or modifications in triticale, and the degree of variation appears to be higher than in other allopolyploid species. The data indicate that at least some sequence changes are non-random, and appear to be a function of genome relations, ploidy levels and sequence types. Specifically, the rye parental genome demonstrated a higher level of changes than the wheat genome. The frequency of lost parental bands was much higher than the frequency of gained novel bands, suggesting that sequence modification and/or elimination might be a major force causing genome variation in triticale. It was also shown that 68 % of the total changes occurred immediately following wide hybridization, but before chromosome doubling. Genome evolution following chromosome doubling occurred more slowly at a very low rate and the changes were mainly observed in the first five or so generations. The data suggest that cytoplasm and relationships between parental genomes are key factors in determining the direction, amount, timing and rate of genomic sequence variation that occurred during inter-generic allopolyploidization in this system.

Integrated cytogenetic map of mitotic metaphase chromosome 9 of maize: resolution, sensitivity, and banding paint development

To study the correlation of the sequence positions on the physical DNA finger print contig (FPC) map and cytogenetic maps of pachytene and somatic maize chromosomes, sequences located along the chromosome 9 FPC map approximately every 10 Mb were selected to place on maize chromosomes using fluorescent in situ hybridization (FISH). The probes were produced as pooled polymerase chain reaction products based on sequences of genetic markers or repeat-free portions of mapped bacterial artificial chromosome (BAC) clones. Fifteen probes were visualized on chromosome 9. The cytological positions of most sequences correspond on the pachytene, somatic, and FPC maps except some probes at the pericentromeric regions. Because of unequal condensation of mitotic metaphase chromosomes, being lower at pericentromeric regions and higher in the arms, probe positions are displaced to the distal ends of both arms. The axial resolution of FISH on somatic chromosome 9 varied from 3.3 to 8.2 Mb, which is 12-30 times lower than on pachytene chromosomes. The probe collection can be used as chromosomal landmarks or as a "banding paint" for the physical mapping of sequences including transgenes and BAC clones and for studying chromosomal rearrangements.

A genomic and expression compendium of the expanded PEBP gene family from maize

The phosphatidylethanolamine-binding proteins (PEBPs) represent an ancient protein family found across the biosphere. In animals they are known to act as kinase and serine protease inhibitors controlling cell growth and differentiation. In plants the most extensively studied PEBP genes, the Arabidopsis (Arabidopsis thaliana) FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) genes, function, respectively, as a promoter and a repressor of the floral transition. Twenty-five maize (Zea mays) genes that encode PEBP-like proteins, likely the entire gene family, were identified and named Zea mays CENTRORADIALIS (ZCN), after the first described plant PEBP gene from Antirrhinum. The maize family is expanded relative to eudicots (typically six to eight genes) and rice (Oryza sativa; 19 genes). Genomic structures, map locations, and syntenous relationships with rice were determined for 24 of the maize ZCN genes. Phylogenetic analysis assigned the maize ZCN proteins to three major subfamilies: TFL1-like (six members), MOTHER OF FT AND TFL1-like (three), and FT-like (15). Expression analysis demonstrated transcription for at least 21 ZCN genes, many with developmentally specific patterns and some having alternatively spliced transcripts. Expression patterns and protein structural analysis identified maize candidates likely having conserved gene function of TFL1. Expression patterns and interaction of the ZCN8 protein with the floral activator DLF1 in the yeast (Saccharomyces cerevisiae) two-hybrid assay strongly supports that ZCN8 plays an orthologous FT function in maize. The expression of other ZCN genes in roots, kernels, and flowers implies their involvement in diverse developmental processes.

A transgenomic cytogenetic sorghum (Sorghum propinquum) bacterial artificial chromosome fluorescence in situ hybridization map of maize (Zea mays L.) pachytene chromosome 9, evidence for regions of genome hyperexpansion

A cytogenetic FISH map of maize pachytene-stage chromosome 9 was produced with 32 maize marker-selected sorghum BACs as probes. The genetically mapped markers used are distributed along the linkage maps at an average spacing of 5 cM. Each locus was mapped by means of multicolor direct FISH with a fluorescently labeled probe mix containing a whole-chromosome paint, a single sorghum BAC clone, and the centromeric sequence, CentC. A maize-chromosome-addition line of oat was used for bright unambiguous identification of the maize 9 fiber within pachytene chromosome spreads. The locations of the sorghum BAC-FISH signals were determined, and each new cytogenetic locus was assigned a centiMcClintock position on the short (9S) or long (9L) arm. Nearly all of the markers appeared in the same order on linkage and cytogenetic maps but at different relative positions on the two. The CentC FISH signal was localized between cdo17 (at 9L.03) and tda66 (at 9S.03). Several regions of genome hyperexpansion on maize chromosome 9 were found by comparative analysis of relative marker spacing in maize and sorghum. This transgenomic cytogenetic FISH map creates anchors between various maps of maize and sorghum and creates additional tools and information for understanding the structure and evolution of the maize genome.

Genome halving with an outgroup

Some genomes are known to have incurred a genome doubling (tetraploidization) event in their evolutionary history, and this is reflected today in patterns of duplicated segments scattered throughout their chromosomes. These duplications may be used as data to "halve" the genome, i.e. to reconstruct the an cestral genome at the moment of tetraploidization, but the solution is often highly non-unique. To resolve this problem, we adapt the genome halving algorithm of El-Mabrouk and Sankoff to take account of an external reference genome. We apply this to reconstruct the tetraploid ancestor of maize, using either rice or sorghum as the reference.

A BAC pooling strategy combined with PCR-based screenings in a large, highly repetitive genome enables integration of the maize genetic and physical maps

Background

Molecular markers serve three important functions in physical map assembly. First, they provide anchor points to genetic maps facilitating functional genomic studies. Second, they reduce the overlap required for BAC contig assembly from 80 to 50 percent. Finally, they validate assemblies based solely on BAC fingerprints. We employed a six-dimensional BAC pooling strategy in combination with a high-throughput PCR-based screening method to anchor the maize genetic and physical maps.

Results

A total of 110,592 maize BAC clones (~ 6x haploid genome equivalents) were pooled into six different matrices, each containing 48 pools of BAC DNA. The quality of the BAC DNA pools and their utility for identifying BACs containing target genomic sequences was tested using 254 PCR-based STS markers. Five types of PCR-based STS markers were screened to assess potential uses for the BAC pools. An average of 4.68 BAC clones were identified per marker analyzed. These results were integrated with BAC fingerprint data generated by the Arizona Genomics Institute (AGI) and the Arizona Genomics Computational Laboratory (AGCoL) to assemble the BAC contigs using the FingerPrinted Contigs (FPC) software and contribute to the construction and anchoring of the physical map. A total of 234 markers (92.5%) anchored BAC contigs to their genetic map positions. The results can be viewed on the integrated map of maize [1,2].

Conclusion

This BAC pooling strategy is a rapid, cost effective method for genome assembly and anchoring. The requirement for six replicate positive amplifications makes this a robust method for use in large genomes with high amounts of repetitive DNA such as maize. This strategy can be used to physically map duplicate loci, provide order information for loci in a small genetic interval or with no genetic recombination, and loci with conflicting hybridization-based information.

Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids

Polyploidy, or whole-genome duplication (WGD), is an important genomic feature for all eukaryotes, especially many plants and some animals. The common occurrence of polyploidy suggests an evolutionary advantage of having multiple sets of genetic material for adaptive evolution. However, increased gene and genome dosages in autopolyploids (duplications of a single genome) and allopolyploids (combinations of two or more divergent genomes) often cause genome instabilities, chromosome imbalances, regulatory incompatibilities, and reproductive failures. Therefore, new allopolyploids must establish a compatible relationship between alien cytoplasm and nuclei and between two divergent genomes, leading to rapid changes in genome structure, gene expression, and developmental traits such as fertility, inbreeding, apomixis, flowering time, and hybrid vigor. Although the underlying mechanisms for these changes are poorly understood, some themes are emerging. There is compelling evidence that changes in DNA sequence, cis- and trans-acting effects, chromatin modifications, RNA-mediated pathways, and regulatory networks modulate differential expression of homoeologous genes and phenotypic variation that may facilitate adaptive evolution in polyploid plants and domestication in crops.

Subfunctionalization of PhyB1 and PhyB2 in the control of seedling and mature plant traits in maize

Phytochromes are the primary red/far-red photoreceptors of higher plants, mediating numerous developmental processes throughout the life cycle, from germination to flowering. In seed plants, phytochromes are encoded by a small gene family with each member performing both distinct and redundant roles in mediating physiological responses to light cues. Studies in both eudicot and monocot species have defined a central role for phytochrome B in mediating responses to light in the control of several agronomically important traits, including plant height, transitions to flowering and axillary branch meristem development. Here we characterize Mutator-induced alleles of PhyB1 and a naturally occurring deletion allele of PhyB2 in Zea mays (maize). Using single and double mutants, we show that the highly similar PhyB1 and PhyB2 genes encode proteins with both overlapping and non-redundant functions that control seedling and mature plant traits. PHYB1 and PHYB2 regulate elongation of sheath and stem tissues of mature plants and contribute to the light-mediated regulation of PhyA and Cab gene transcripts. However, PHYB1 and not PHYB2 contributes significantly to the inhibition of mesocotyl elongation under red light, whereas PHYB2 and to a lesser extent PHYB1 mediate the photoperiod-dependent floral transition. This sub functionalization of PHYB activities in maize has probably occurred since the tetraploidization of maize, and may contribute to flowering time variation in modern-day varieties.

Comparative genome mapping among Picea glauca, P. mariana x P. rubens and P. abies, and correspondence with other Pinaceae

A composite linkage map was constructed from four individual maps for the conifer Picea glauca (Moench) Voss, from anonymous and gene-specific markfers (714 AFLPs, 38 SSRs, and 53 ESTPs). A total of 12 linkage groups were delineated with an average marker density of 2.7 cM. Macro-synteny and macro-colinearity comparisons with two other composite linkage maps developed for the species complex P. mariana (Mill.) B.S.P. x P. rubens Sarg., and for P. abies (L.) Karst. revealed an identical number of linkage groups and a remarkable conservation of the gene content and gene order of linkage groups over the million years since the split between these taxa. Identical gene order among taxa was observed for 10 of the 12 assembled composite linkage groups. The discovery of one breakdown in synteny between P. glauca and the other two taxa indicated the occurrence of an inter-chromosomal rearrangement involving an insertional translocation. Analysis of marker colinearity also revealed a putative segmental duplication. The combined information from these three Picea genomes validated and improved large-scale genome comparisons at the inter-generic level in the family Pinaceae by allowing for the identification of 11 homoeologous linkage groups between Picea and Pinus, and nine such groups between Picea and Pseudotsuga menziesii. Notably, the analysis of synteny among the three genera revealed a putative case of chromosomal fission and an inter-chromosomal rearrangement in the genome of P. menziesii. Both of these changes are inter-connected, indicating much instability in this part of the P. menziesii genome. Overall, the macro-structure of the Pinaceae genome was well conserved, which is notable given the Cretaceous origin of its main lineages.

Statistical inference of chromosomal homology based on gene colinearity and applications to Arabidopsis and rice

Background

The identification of chromosomal homology will shed light on such mysteries of genome evolution as DNA duplication, rearrangement and loss. Several approaches have been developed to detect chromosomal homology based on gene synteny or colinearity. However, the previously reported implementations lack statistical inferences which are essential to reveal actual homologies.

Results

In this study, we present a statistical approach to detect homologous chromosomal segments based on gene colinearity. We implement this approach in a software package ColinearScan to detect putative colinear regions using a dynamic programming algorithm. Statistical models are proposed to estimate proper parameter values and evaluate the significance of putative homologous regions. Statistical inference, high computational efficiency and flexibility of input data type are three key features of our approach.

Conclusion

We apply ColinearScan to the Arabidopsis and rice genomes to detect duplicated regions within each species and homologous fragments between these two species. We find many more homologous chromosomal segments in the rice genome than previously reported. We also find many small colinear segments between rice and Arabidopsis genomes.

Genomic affinities between maize and Zea perennis using classical and molecular cytogenetic methods (GISH-FISH)

In this study we have analysed and compared the genomic composition, meiotic behaviour, and meiotic affinities of Zea perennis and Zea mays ssp. mays. To do so we studied the parental taxa and the interspecific hybrid Zea perennis x Zea mays ssp. mays, using classical cytogenetic methods, as well as GISH and FISH. GISH enabled us to recognize the genomic source of each chromosome involved in the meiotic configurations of this hybrid, and established the genomic affinities between their parental species. The results obtained here reinforce the hypothesis of the amphiploid origin of Zea perennis and, together with previous research, indicate that the chromosomes with divergent repetitive sequences in maize and Zea luxurians could be the remnants of a relict parental genome not shared with Zea perennis.

New evidence from Sinapis alba L. for ancestral triplication in a crucifer genome

We present clear evidence of ancestral genome triplication in Sinapis alba, a close relative of the cultivated Brassica species. Exceptionally high levels of heterozygosity in the parents of an F1 intercross permitted the mapping of an estimated 87% of all detected restriction fragment length polymorphism (RFLP) loci, with each RFLP probe typically detecting 2 or 3 loci. These duplicated loci were arranged in 8 triplicated homologous linkage blocks and 2 small, duplicated, homologous linkage blocks covering the majority of the S. alba genome. Several large-scale inversions and translocations appear to have rearranged the order of loci within homologous blocks. The role of successive polyploidization events on the evolution of crucifer species is discussed.

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).