Organization and variability of the maize genome

Biochemical and Epigenetic Regulation of Glutamate Metabolism in Maize (Zea mays L.) Leaves under Salt Stress

The effect of salt stress (150 mM NaCl) on the expression of genes, methylation of their promoters, and enzymatic activity of glutamate dehydrogenase (GDH), glutamate decarboxylase (GAD), and the 2-oxoglutarate (2-OG)-dehydrogenase (2-OGDH) complex was studied in maize (Zea mays L.). GDH activity increased continuously under salt stress, being 3-fold higher after 24 h. This was accompanied by the appearance of a second isoform with lower electrophoretic mobility. The expression of the Gdh1 gene strongly increased after 6-12 h of incubation, which corresponded to the demethylation of its promoter, while Gdh2 gene expression slightly increased after 2-6 h and then decreased. GAD activity gradually increased in the first 12 h, and then returned to the control level. This corresponded to the increase of Gad expression and its demethylation. Salt stress led to a 2-fold increase in the activity of 2-OGDH during the first 6 h of NaCl treatment, then the activity returned to the control level. Expression of the genes Ogdh1 and Ogdh3 peaked after 1-2 h of incubation. After 6-8 h with NaCl, the expression of these genes declined below the control levels, which correlated with the higher methylation of their promoters. We conclude that salt stress causes a redirection of the 2-OG flux to the γ-aminobutyric acid shunt via its amination to glutamate, by altering the expression of the Gdh1 and Gdh2 genes, which likely promotes the assembly of the native GDH molecule having a different subunit composition and greater affinity for 2-OG.

Genetic Diversity within a Collection of Italian Maize Inbred Lines: A Resource for Maize Genomics and Breeding

Genetic diversity is fundamental for studying the complex architecture of the traits of agronomic importance, controlled by major and minor loci. Moreover, well-characterized germplasm collections are essential tools for dissecting and analyzing genetic and phenotypic diversity in crops. A panel of 360 entries, a subset of a larger collection maintained within the GenBank at CREA Bergamo, which includes the inbreds derived from traditional Italian maize open-pollinated (OP) varieties and advanced breeding ones (Elite Inbreds), was analyzed to identify SNP markers using the tGBS® genotyping-by-sequencing technology. A total of 797,368 SNPs were found during the initial analysis. Imputation and filtering processes were carried out based on the percentage of missing data, redundant markers, and rarest allele frequencies, resulting in a final dataset of 15,872 SNP markers for which a physical map position was identified. Using this dataset, the inbred panel was characterized for linkage disequilibrium (LD), genetic diversity, population structure, and genetic relationships. LD decay at a genome-wide level indicates that the collection is a suitable resource for association mapping. Population structure analyses, which were carried out with different clustering methods, showed stable grouping statistics for four groups, broadly corresponding to 'Insubria', 'Microsperma', and 'Scagliolino' genotypes, with a fourth group composed prevalently of elite accessions derived from Italian and US breeding programs. Based on these results, the CREA Italian maize collection, genetically characterized in this study, can be considered an important tool for the mapping and characterization of useful traits and associated loci/alleles, to be used in maize breeding programs.

On the usefulness of mock genomes to define heterotic pools, testers, and hybrid predictions in orphan crops

The advances in genomics in recent years have increased the accuracy and efficiency of breeding programs for many crops. Nevertheless, the adoption of genomic enhancement for several other crops essential in developing countries is still limited, especially for those that do not have a reference genome. These crops are more often called orphans. This is the first report to show how the results provided by different platforms, including the use of a simulated genome, called the mock genome, can generate in population structure and genetic diversity studies, especially when the intention is to use this information to support the formation of heterotic groups, choice of testers, and genomic prediction of single crosses. For that, we used a method to assemble a reference genome to perform the single-nucleotide polymorphism (SNP) calling without needing an external genome. Thus, we compared the analysis results using the mock genome with the standard approaches (array and genotyping-by-sequencing (GBS)). The results showed that the GBS-Mock presented similar results to the standard methods of genetic diversity studies, division of heterotic groups, the definition of testers, and genomic prediction. These results showed that a mock genome constructed from the population's intrinsic polymorphisms to perform the SNP calling is an effective alternative for conducting genomic studies of this nature in orphan crops, especially those that do not have a reference genome.

Megabase-scale presence-absence variation with Tripsacum origin was under selection during maize domestication and adaptation

BACKGROUND

Structural variants (SVs) significantly drive genome diversity and environmental adaptation for diverse species. Unlike the prevalent small SVs (< kilobase-scale) in higher eukaryotes, large-size SVs rarely exist in the genome, but they function as one of the key evolutionary forces for speciation and adaptation.

RESULTS

In this study, we discover and characterize several megabase-scale presence-absence variations (PAVs) in the maize genome. Surprisingly, we identify a 3.2 Mb PAV fragment that shows high integrity and is present as complete presence or absence in the natural diversity panel. This PAV is embedded within the nucleolus organizer region (NOR), where the suppressed recombination is found to maintain the PAV against the evolutionary variation. Interestingly, by analyzing the sequence of this PAV, we not only reveal the domestication trace from teosinte to modern maize, but also the footprints of its origin from Tripsacum, shedding light on a previously unknown contribution from Tripsacum to the speciation of Zea species. The functional consequence of the Tripsacum segment migration is also investigated, and environmental fitness conferred by the PAV may explain the whole segment as a selection target during maize domestication and improvement.

CONCLUSIONS

These findings provide a novel perspective that Tripsacum contributes to Zea speciation, and also instantiate a strategy for evolutionary and functional analysis of the "fossil" structure variations during genome evolution and speciation.

Genotyping-by-sequencing and SNP-arrays are complementary for detecting quantitative trait loci by tagging different haplotypes in association studies

BACKGROUND

Single Nucleotide Polymorphism (SNP) array and re-sequencing technologies have different properties (e.g. calling rate, minor allele frequency profile) and drawbacks (e.g. ascertainment bias). This lead us to study their complementarity and the consequences of using them separately or combined in diversity analyses and Genome-Wide Association Studies (GWAS). We performed GWAS on three traits (grain yield, plant height and male flowering time) measured in 22 environments on a panel of 247 F1 hybrids obtained by crossing 247 diverse dent maize inbred lines with a same flint line. The 247 lines were genotyped using three genotyping technologies (Genotyping-By-Sequencing, Illumina Infinium 50 K and Affymetrix Axiom 600 K arrays).

RESULTS

The effects of ascertainment bias of the 50 K and 600 K arrays were negligible for deciphering global genetic trends of diversity and for estimating relatedness in this panel. We developed an original approach based on linkage disequilibrium (LD) extent in order to determine whether SNPs significantly associated with a trait and that are physically linked should be considered as a single Quantitative Trait Locus (QTL) or several independent QTLs. Using this approach, we showed that the combination of the three technologies, which have different SNP distributions and densities, allowed us to detect more QTLs (gain in power) and potentially refine the localization of the causal polymorphisms (gain in resolution).

CONCLUSIONS

Conceptually different technologies are complementary for detecting QTLs by tagging different haplotypes in association studies. Considering LD, marker density and the combination of different technologies (SNP-arrays and re-sequencing), the genotypic data available were most likely enough to well represent polymorphisms in the centromeric regions, whereas using more markers would be beneficial for telomeric regions.

The maize W22 genome provides a foundation for functional genomics and transposon biology

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

Sequence analysis of European maize inbred line F2 provides new insights into molecular and chromosomal characteristics of presence/absence variants

BACKGROUND

Maize is well known for its exceptional structural diversity, including copy number variants (CNVs) and presence/absence variants (PAVs), and there is growing evidence for the role of structural variation in maize adaptation. While PAVs have been described in this important crop species, they have been only scarcely characterized at the sequence level and the extent of presence/absence variation and relative chromosomal landscape of inbred-specific regions remain to be elucidated.

RESULTS

De novo genome sequencing of the French F2 maize inbred line revealed 10,044 novel genomic regions larger than 1 kb, making up 88 Mb of DNA, that are present in F2 but not in B73 (PAV). This set of maize PAV sequences allowed us to annotate PAV content and to analyze sequence breakpoints. Using PAV genotyping on a collection of 25 temperate lines, we also analyzed Linkage Disequilibrium in PAVs and flanking regions, and PAV frequencies within maize genetic groups.

CONCLUSIONS

We highlight the possible role of MMEJ-type double strand break repair in maize PAV formation and discover 395 new genes with transcriptional support. Pattern of linkage disequilibrium within PAVs strikingly differs from this of flanking regions and is in accordance with the intuition that PAVs may recombine less than other genomic regions. We show that most PAVs are ancient, while some are found only in European Flint material, thus pinpointing structural features that may be at the origin of adaptive traits involved in the success of this material. Characterization of such PAVs will provide useful material for further association genetic studies in European and temperate maize.

Integrative analysis of DNA methylation, mRNAs, and small RNAs during maize embryo dedifferentiation

BACKGROUND

Maize (Zea mays) is an important model crop for transgenic studies. However, genetic transformation of maize requires embryonic calli derived from immature embryo, and the impact of utilizing tissue culture methods on the maize epigenome is poorly understood. Here, we generated whole-genome MeDIP-seq data examining DNA methylation in dedifferentiated and normal immature maize embryos.

RESULTS

We observed that most of the dedifferentiated embryos exhibited a methylation increase compared to normal embryos. Increased methylation at promoters was associated with down-regulated protein-coding gene expression; however, the correlation was not strong. Analysis of the callus and immature embryos indicated that the methylation increase was induced during induction of embryonic callus, suggesting phenotypic consequences may be caused by perturbations in genomic DNA methylation levels. The correlation between the 21-24nt small RNAs and DNA methylation regions were investigated but only a statistically significant correlation for 24nt small RNAs was observed.

CONCLUSIONS

These data extend the significance of epigenetic changes during maize embryo callus formation, and the methylation changes might explain some of the observed embryonic callus variation in callus formation.

Variable Level of Dominance of Candidate Genes Controlling Drought Functional Traits in Maize Hybrids

Breeding maize for drought tolerance necessitates the knowledge on tolerant genotypes, molecular basis of drought tolerance mechanism, action, and expression pattern of genes. Studying the expression pattern and gene action of candidate genes during drought stress in the hybrids will help in choosing target genes for drought tolerance breeding. In the present investigation, a set of five hybrids and their seven parents with a variable level of tolerance to drought stress was selected to study the magnitude and the direction of 52 drought-responsive candidate genes distributed across various biological functions, viz., stomatal regulation, root development, detoxification, hormone signaling, photosynthesis, and sugar metabolism. The tolerant parents, HKI1105 and CML425, and their hybrid, ADWLH2, were physiologically active under drought stress, since vital parameters viz., chlorophyll, root length and relative water content, were on par with the respective well-watered control. All the genes were up-regulated in ADWLH2, many were down-regulated in HM8 and HM9, and most were down-regulated in PMH1 and PMH3 in the shoots and roots. The nature of the gene action was controlled by the parental combination rather than the parent per se. The differentially expressed genes in all five hybrids explained a mostly non-additive gene action over additivity, which was skewed toward any of the parental lines. Tissue-specific gene action was also noticed in many of the genes. The non-additive gene action is driven by genetic diversity, allele polymorphism, events during gene regulation, and small RNAs under the stress condition. Differential regulation and cross-talk of genes controlling various biological functions explained the basis of drought tolerance in subtropical maize hybrids. The nature of the gene action and the direction of the expression play crucial roles in designing introgression and hybrid breeding programmes to breed drought tolerant maize hybrids.

European flint landraces grown in situ reveal adaptive introgression from modern maize

We have investigated the role of selection in the determination of the detected levels of introgression from modern maize hybrid varieties into maize landraces still cultivated in situ in Italy. We exploited the availability of a historical collection of landraces undertaken before the introduction and widespread use of modern maize, to analyse genomic changes that have occurred in these maize landraces over 50 years of co-existence with hybrid varieties. We have combined a previously published SSR dataset (n=21) with an AFLP loci dataset (n=168) to provide higher resolution power and to obtain a more detailed picture. We show that selection pressures for adaptation have favoured new alleles introduced by migration from hybrids. This shows the potential for analysis of historical introgression even over this short period of 50 years, for an understanding of the evolution of the genome and for the identification of its functionally important regions. Moreover, this demonstrates that landraces grown in situ represent almost unique populations for use for such studies when the focus is on the domesticated plant. This is due to their adaptation, which has arisen from their dynamic evolution under a continuously changing agro-ecological environment, and their capture of new alleles from hybridisation. We have also identified loci for which selection has inhibited introgression from modern germplasm and has enhanced the distinction between landraces and modern maize. These loci indicate that selection acted in the past, during the formation of the flint and dent gene pools. In particular, the locus showing the strongest signals of selection is a Misfit transposable element. Finally, molecular characterisation of the same samples with two different molecular markers has allowed us to compare their performances. Although the genetic-diversity and population-structure analyses provide the same global qualitative pattern, which thus provides the same inferences, there are differences related to their natures and characteristics.

Genomic distribution of H3K9me2 and DNA methylation in a maize genome

DNA methylation and dimethylation of lysine 9 of histone H3 (H3K9me2) are two chromatin modifications that can be associated with gene expression or recombination rate. The maize genome provides a complex landscape of interspersed genes and transposons. The genome-wide distribution of DNA methylation and H3K9me2 were investigated in seedling tissue for the maize inbred B73 and compared to patterns of these modifications observed in Arabidopsis thaliana. Most maize transposons are highly enriched for DNA methylation in CG and CHG contexts and for H3K9me2. In contrast to findings in Arabidopsis, maize CHH levels in transposons are generally low but some sub-families of transposons are enriched for CHH methylation and these families exhibit low levels of H3K9me2. The profile of modifications over genes reveals that DNA methylation and H3K9me2 is quite low near the beginning and end of genes. Although elevated CG and CHG methylation are found within gene bodies, CHH and H3K9me2 remain low. Maize has much higher levels of CHG methylation within gene bodies than observed in Arabidopsis and this is partially attributable to the presence of transposons within introns for some maize genes. These transposons are associated with high levels of CHG methylation and H3K9me2 but do not appear to prevent transcriptional elongation. Although the general trend is for a strong depletion of H3K9me2 and CHG near the transcription start site there are some putative genes that have high levels of these chromatin modifications. This study provides a clear view of the relationship between DNA methylation and H3K9me2 in the maize genome and how the distribution of these modifications is shaped by the interplay of genes and transposons.

Genomic localization of AtRE1 and AtRE2, copia-type retrotransposons, in natural variants of Arabidopsis thaliana

Retrotransposons are ubiquitous components of plant genomes. They affect genome organization, and can also affect the expression patterns of neighboring genes. Retrotransposons are therefore important elements for changing genomic information. To understand the evolution of the Arabidopsis genome, we examined the distribution of certain retrotransposons, AtRE1s and AtRE2s, in the genomes of 12 natural variants (accessions) of Arabidopsis thaliana. AtRE1 and AtRE2 are copia-type retrotransposons that are potentially active. Their copy numbers are low, and they are absent from the genomes of some accessions. We detected four loci with AtRE1s inserted in six accessions, and one locus with an insertion of a solo-LTR-like sequence derived from AtRE1 in two accessions. Seven loci with AtRE2s inserted were detected on eight accessions. These loci were distributed in euchromatic regions of chromosomes 1, 2, 3, and 4. The AtRE1 and AtRE2 sequences at some loci identified in this study have not been recorded in the database of the 1001 Genome project. The sequences of AtRE1s and those of AtRE2s in different accessions and at different loci were highly conserved. There was a complete or almost complete conservation of sequences of both long terminal repeats in each AtRE1 and in each AtRE2. These results suggest that AtRE1 and AtRE2 appeared quite recently in the Arabidopsis genome. Furthermore, sequence comparisons of AtRE1 and AtRE2 loci among accessions revealed the possibility that large deletions containing entire sequences of AtRE1 and AtRE2 have occurred in some accessions.

Comprehensive analysis of imprinted genes in maize reveals allelic variation for imprinting and limited conservation with other species

In plants, a subset of genes exhibit imprinting in endosperm tissue such that expression is primarily from the maternal or paternal allele. Imprinting may arise as a consequence of mechanisms for silencing of transposons during reproduction, and in some cases imprinted expression of particular genes may provide a selective advantage such that it is conserved across species. Separate mechanisms for the origin of imprinted expression patterns and maintenance of these patterns may result in substantial variation in the targets of imprinting in different species. Here we present deep sequencing of RNAs isolated from reciprocal crosses of four diverse maize genotypes, providing a comprehensive analysis that allows evaluation of imprinting at more than 95% of endosperm-expressed genes. We find that over 500 genes exhibit statistically significant parent-of-origin effects in maize endosperm tissue, but focused our analyses on a subset of these genes that had >90% expression from the maternal allele (69 genes) or from the paternal allele (108 genes) in at least one reciprocal cross. Over 10% of imprinted genes show evidence of allelic variation for imprinting. A comparison of imprinting in maize and rice reveals that 13% of genes with syntenic orthologs in both species exhibit conserved imprinting. Genes that exhibit conserved imprinting between maize and rice have elevated nonsynonymous to synonymous substitution ratios compared with other imprinted genes, suggesting a history of more rapid evolution. Together, these data suggest that imprinting only has functional relevance at a subset of loci that currently exhibit imprinting in maize.

Randomly detected genetically modified (GM) maize (Zea mays L.) near a transport route revealed a fragile 45S rDNA phenotype

Monitoring of genetically modified (GM) crops has been emphasized to prevent their potential effects on the environment and human health. Monitoring of the inadvertent dispersal of transgenic maize in several fields and transport routes in Korea was carried out by qualitative multiplex PCR, and molecular analyses were conducted to identify the events of the collected GM maize. Cytogenetic investigations through fluorescence in situ hybridization (FISH) of the GM maize were performed to check for possible changes in the 45S rDNA cluster because this cluster was reported to be sensitive to replication and transcription stress. Three GM maize kernels were collected from a transport route near Incheon port, Korea, and each was found to contain NK603, stacked MON863 x NK603, and stacked NK603 x MON810 inserts, respectively. Cytogenetic analysis of the GM maize containing the stacked NK603 x MON810 insert revealed two normal compact 5S rDNA signals, but the 45S rDNA showed a fragile phenotype, demonstrating a “beads-on-a-string” fragmentation pattern, which seems to be a consequence of genetic modification. Implications of the 45S rDNA cluster fragility in GM maize are also discussed.

Progress toward understanding heterosis in crop plants

Although heterosis, or hybrid vigor, is widely exploited in agriculture, a complete description of its molecular underpinnings has remained elusive despite extensive investigation. It appears that there is not a single, simple explanation for heterosis. Instead, it is likely that heterosis arises in crosses between genetically distinct individuals as a result of a diversity of mechanisms. Heterosis generally results from the action of multiple loci, and different loci affect heterosis for different traits and in different hybrids. Hence, multigene models are likely to prove most informative for understanding heterosis. Complementation of allelic variation, as well as complementation of variation in gene content and gene expression patterns, is likely to be an important contributor to heterosis. Epigenetic variation has the potential to interact in hybrid genotypes via novel mechanisms. Several other intriguing hypotheses are also under investigation. In crops, heterosis must be considered within the context of the genomic impacts of prior selection for agronomic traits.

Rogowsky P. PPR8522 encodes a chloroplast-targeted pentatricopeptide repeat protein necessary for maize embryogenesis and vegetative development

The pentatricopeptide repeat (PPR) domain is an RNA binding domain allowing members of the PPR superfamily to participate in post-transcriptional processing of organellar RNA. Loss of PPR8522 from maize (Zea mays) confers an embryo-specific (emb) phenotype. The emb8522 mutation was isolated in an active Mutator (Mu) population and co-segregation analysis revealed that it was tightly linked to a MuDR insertion in the first exon of PPR8522. Independent evidence that disruption of PPR8522 caused the emb phenotype was provided by fine mapping to a region of 116kb containing no other gene than PPR8522 and complementation of the emb8522 mutant by a PPR8522 cDNA. The deduced PPR8522 amino acid sequence of 832 amino acids contains 10 PPR repeats and a chloroplast target peptide, the function of which was experimentally demonstrated by transient expression in Nicotiana benthamiana. Whereas mutant endosperm is apparently normal, mutant embryos deviate from normal development as early as 3 days after pollination, are reduced in size, exhibit more or less severe morphological aberrations depending on the genetic background, and generally do not germinate. The emb8522 mutation is the first to associate the loss of a PPR gene with an embryo-lethal phenotype in maize. Analyses of mutant plantlets generated by embryo-rescue experiments indicate that emb8522 also affects vegetative plant growth and chloroplast development. The loss of chloroplast transcription dependent on plastid-encoded RNA polymerase is the likely cause for the lack of an organized thylakoid network and an albino, seedling-lethal phenotype.

Repeat associated small RNAs vary among parents and following hybridization in maize

Small RNAs (sRNAs) are hypothesized to contribute to hybrid vigor because they maintain genome integrity, contribute to genetic diversity, and control gene expression. We used Illumina sequencing to assess how sRNA populations vary between two maize inbred lines (B73 and Mo17) and their hybrid. We sampled sRNAs from the seedling shoot apex and the developing ear, two rapidly growing tissues that program the greater growth of maize hybrids. We found that parental differences in siRNAs primarily originate from repeat regions. Although the maize genome contains greater number and complexity of repeats compared with Arabidopsis or rice, we confirmed that, like these simpler plant genomes, 24-nt siRNAs whose abundance differs between maize parents also show a trend of down-regulation following hybridization. Surprisingly, hybrid vigor is fully maintained when 24-nt siRNAs are globally reduced by mutation of the RNA-dependent RNA polymerase 2 encoded by modifier of paramutation1 (mop1). We also discovered that 21-22-nt siRNAs derived from a number of distinct retrotransposon families differentially accumulate between B73 and Mo17 as well as their hybrid. Thus, maize possesses a unique source of genetic variation for regulating transposons and genes at a genomic scale, which may contribute to its high degree of observed heterosis.

A large maize (Zea mays L.) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome

SNP genotyping arrays have been useful for many applications that require a large number of molecular markers such as high-density genetic mapping, genome-wide association studies (GWAS), and genomic selection. We report the establishment of a large maize SNP array and its use for diversity analysis and high density linkage mapping. The markers, taken from more than 800,000 SNPs, were selected to be preferentially located in genes and evenly distributed across the genome. The array was tested with a set of maize germplasm including North American and European inbred lines, parent/F1 combinations, and distantly related teosinte material. A total of 49,585 markers, including 33,417 within 17,520 different genes and 16,168 outside genes, were of good quality for genotyping, with an average failure rate of 4% and rates up to 8% in specific germplasm. To demonstrate this array's use in genetic mapping and for the independent validation of the B73 sequence assembly, two intermated maize recombinant inbred line populations - IBM (B73×Mo17) and LHRF (F2×F252) - were genotyped to establish two high density linkage maps with 20,913 and 14,524 markers respectively. 172 mapped markers were absent in the current B73 assembly and their placement can be used for future improvements of the B73 reference sequence. Colinearity of the genetic and physical maps was mostly conserved with some exceptions that suggest errors in the B73 assembly. Five major regions containing non-colinearities were identified on chromosomes 2, 3, 6, 7 and 9, and are supported by both independent genetic maps. Four additional non-colinear regions were found on the LHRF map only; they may be due to a lower density of IBM markers in those regions or to true structural rearrangements between lines. Given the array's high quality, it will be a valuable resource for maize genetics and many aspects of maize breeding.

Heritable epigenetic variation among maize inbreds

Epigenetic variation describes heritable differences that are not attributable to changes in DNA sequence. There is the potential for pure epigenetic variation that occurs in the absence of any genetic change or for more complex situations that involve both genetic and epigenetic differences. Methylation of cytosine residues provides one mechanism for the inheritance of epigenetic information. A genome-wide profiling of DNA methylation in two different genotypes of Zea mays (ssp. mays), an organism with a complex genome of interspersed genes and repetitive elements, allowed the identification and characterization of examples of natural epigenetic variation. The distribution of DNA methylation was profiled using immunoprecipitation of methylated DNA followed by hybridization to a high-density tiling microarray. The comparison of the DNA methylation levels in the two genotypes, B73 and Mo17, allowed for the identification of approximately 700 differentially methylated regions (DMRs). Several of these DMRs occur in genomic regions that are apparently identical by descent in B73 and Mo17 suggesting that they may be examples of pure epigenetic variation. The methylation levels of the DMRs were further studied in a panel of near-isogenic lines to evaluate the stable inheritance of the methylation levels and to assess the contribution of cis- and trans- acting information to natural epigenetic variation. The majority of DMRs that occur in genomic regions without genetic variation are controlled by cis-acting differences and exhibit relatively stable inheritance. This study provides evidence for naturally occurring epigenetic variation in maize, including examples of pure epigenetic variation that is not conditioned by genetic differences. The epigenetic differences are variable within maize populations and exhibit relatively stable trans-generational inheritance. The detected examples of epigenetic variation, including some without tightly linked genetic variation, may contribute to complex trait variation.

Heritable variation within a species provides the basis for natural and artificial selection. A substantial portion of heritable variation is based on alterations in DNA sequence among individuals and is termed genetic variation. There is also evidence for epigenetic variation, which refers to heritable differences that are not caused by DNA sequence changes. Methylation of cytosine residues provides one molecular mechanism for epigenetic variation in many eukaryotic species. The genome-wide distribution of DNA methylation was assessed in two different inbred genotypes of maize to identify differentially methylated regions that may contribute to epigenetic variation. There are hundreds of genomic regions that have differences in DNA methylation levels in these two different genotypes, including methylation differences in regions without genetic variation. By studying the inheritance of the differential methylation in near-isogenic progeny of the two inbred lines, it is possible to demonstrate relatively stable inheritance of epigenetic variation, even in the absence of DNA sequence changes. The epigenetic variation among individuals of the same species may provide important contributions to phenotypic variation within a species even in the absence of genetic differences.

A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species

Advances in next generation technologies have driven the costs of DNA sequencing down to the point that genotyping-by-sequencing (GBS) is now feasible for high diversity, large genome species. Here, we report a procedure for constructing GBS libraries based on reducing genome complexity with restriction enzymes (REs). This approach is simple, quick, extremely specific, highly reproducible, and may reach important regions of the genome that are inaccessible to sequence capture approaches. By using methylation-sensitive REs, repetitive regions of genomes can be avoided and lower copy regions targeted with two to three fold higher efficiency. This tremendously simplifies computationally challenging alignment problems in species with high levels of genetic diversity. The GBS procedure is demonstrated with maize (IBM) and barley (Oregon Wolfe Barley) recombinant inbred populations where roughly 200,000 and 25,000 sequence tags were mapped, respectively. An advantage in species like barley that lack a complete genome sequence is that a reference map need only be developed around the restriction sites, and this can be done in the process of sample genotyping. In such cases, the consensus of the read clusters across the sequence tagged sites becomes the reference. Alternatively, for kinship analyses in the absence of a reference genome, the sequence tags can simply be treated as dominant markers. Future application of GBS to breeding, conservation, and global species and population surveys may allow plant breeders to conduct genomic selection on a novel germplasm or species without first having to develop any prior molecular tools, or conservation biologists to determine population structure without prior knowledge of the genome or diversity in the species.

Strategies for silencing and escape: the ancient struggle between transposable elements and their hosts

Over the past several years, there has been an explosion in our understanding of the mechanisms by which plant transposable elements (TEs) are epigenetically silenced and maintained in an inactive state over long periods of time. This highly efficient process results in vast numbers of inactive TEs; indeed, the majority of many plant genomes are composed of these quiescent elements. This observation has led to the rather static view that TEs represent an essentially inert portion of plant genomes. However, recent work has demonstrated that TE silencing is a highly dynamic process that often involves transcription of TEs at particular times and places during plant development. Plants appear to use transcripts from silenced TEs as an ongoing source of information concerning the mobile portion of the genome. In contrast to our understanding of silencing pathways, we know relatively little about the ways in which TEs evade silencing. However, vast differences in TE content between even closely related plant species suggest that they are often wildly successful at doing so. Here, we discuss TE activity in plants as the result of a constantly shifting balance between host strategies for TE silencing and TE strategies for escape and amplification.

Digital gene expression signatures for maize development

Genome-wide expression signatures detect specific perturbations in developmental programs and contribute to functional resolution of key regulatory networks. In maize (Zea mays) inflorescences, mutations in the RAMOSA (RA) genes affect the determinacy of axillary meristems and thus alter branching patterns, an important agronomic trait. In this work, we developed and tested a framework for analysis of tag-based, digital gene expression profiles using Illumina's high-throughput sequencing technology and the newly assembled B73 maize reference genome. We also used a mutation in the RA3 gene to identify putative expression signatures specific to stem cell fate in axillary meristem determinacy. The RA3 gene encodes a trehalose-6-phosphate phosphatase and may act at the interface between developmental and metabolic processes. Deep sequencing of digital gene expression libraries, representing three biological replicate ear samples from wild-type and ra3 plants, generated 27 million 20- to 21-nucleotide reads with frequencies spanning 4 orders of magnitude. Unique sequence tags were anchored to 3'-ends of individual transcripts by DpnII and NlaIII digests, which were multiplexed during sequencing. We mapped 86% of nonredundant signature tags to the maize genome, which associated with 37,117 gene models and unannotated regions of expression. In total, 66% of genes were detected by at least nine reads in immature maize ears. We used comparative genomics to leverage existing information from Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) in functional analyses of differentially expressed maize genes. Results from this study provide a basis for the analysis of short-read expression data in maize and resolved specific expression signatures that will help define mechanisms of action for the RA3 gene.

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

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

Overexpression of the epidermis-specific homeodomain-leucine zipper IV transcription factor Outer Cell Layer1 in maize identifies target genes involved in lipid metabolism and cuticle biosynthesis

Transcription factors of the homeodomain-leucine zipper IV (HD-ZIP IV) family play crucial roles in epidermis-related processes. To gain further insight into the molecular function of OUTER CELL LAYER1 (OCL1), 14 target genes up- or down-regulated in transgenic maize (Zea mays) plants overexpressing OCL1 were identified. The 14 genes all showed partial coexpression with OCL1 in maize organs, and several of them shared preferential expression in the epidermis with OCL1. They encoded proteins involved in lipid metabolism, defense, envelope-related functions, or cuticle biosynthesis and include ZmWBC11a (for white brown complex 11a), an ortholog of AtWBC11 involved in the transport of wax and cutin molecules. In support of the annotations, OCL1-overexpressing plants showed quantitative and qualitative changes of cuticular wax compounds in comparison with wild-type plants. An increase in C24 to C28 alcohols was correlated with the transcriptional up-regulation of ZmFAR1, coding for a fatty acyl-coenzyme A reductase. Transcriptional activation of ZmWBC11a by OCL1 was likely direct, since transactivation in transiently transformed maize kernels was abolished by a deletion of the activation domain in OCL1 or mutations in the L1 box, a cis-element bound by HD-ZIP IV transcription factors. Our data demonstrate that, in addition to AP2/EREBP and MYB-type transcription factors, members of the HD-ZIP IV family contribute to the transcriptional regulation of genes involved in cuticle biosynthesis.

Spatial configuration of transposable element Ac termini affects their ability to induce chromosomal breakage in maize

Composite or closely linked maize (Zea mays) Ac/Ds transposable elements can induce chromosome breakage, but the precise configurations of Ac/Ds elements that can lead to chromosome breakage are not completely defined. Here, we determined the structures and chromosome breakage properties of 15 maize p1 alleles: each allele contains a fixed fractured Ac (fAc) element and a closely linked full-length Ac at various flanking sites. Our results show that pairs of Ac/fAc elements in which the termini of different elements are in direct or reverse orientation can induce chromosome breakage. By contrast, no chromosome breakage is observed with alleles containing pairs of Ac/fAc elements in which the external termini of the paired elements can function as a macrotransposon. Among the structures that can lead to chromosome breaks, breakage frequency is inversely correlated with the distance between the interacting Ac/Ds termini. These results provide new insight into the mechanism of transposition-induced chromosome breakage, which is one outcome of the chromosome-restructuring ability of alternative transposition events.

Fine mapping and haplotype structure analysis of a major flowering time quantitative trait locus on maize chromosome 10

Flowering time is a major adaptive trait in plants and an important selection criterion for crop species. In maize, however, little is known about its molecular basis. In this study, we report the fine mapping and characterization of a major quantitative trait locus located on maize chromosome 10, which regulates flowering time through photoperiod sensitivity. This study was performed in near-isogenic material derived from a cross between the day-neutral European flint inbred line FV286 and the tropical short-day inbred line FV331. Recombinant individuals were identified among a large segregating population and their progenies were scored for flowering time. Combined genotypic characterization led to delimit the QTL to an interval of 170 kb and highlighted an unbalanced recombination pattern. Two bacterial artificial chromosomes (BACs) covering the region were analyzed to identify putative candidate genes, and synteny with rice, sorghum, and brachypodium was investigated. A gene encoding a CCT domain protein homologous to the rice Ghd7 heading date regulator was identified, but its causative role was not demonstrated and deserves further analyses. Finally, an association study showed a strong level of linkage disequilibrium over the region and highlighted haplotypes that could provide useful information for the exploitation of genetic resources and marker-assisted selection in maize.

Maize inbreds exhibit high levels of copy number variation (CNV) and presence/absence variation (PAV) in genome content

Following the domestication of maize over the past approximately 10,000 years, breeders have exploited the extensive genetic diversity of this species to mold its phenotype to meet human needs. The extent of structural variation, including copy number variation (CNV) and presence/absence variation (PAV), which are thought to contribute to the extraordinary phenotypic diversity and plasticity of this important crop, have not been elucidated. Whole-genome, array-based, comparative genomic hybridization (CGH) revealed a level of structural diversity between the inbred lines B73 and Mo17 that is unprecedented among higher eukaryotes. A detailed analysis of altered segments of DNA conservatively estimates that there are several hundred CNV sequences among the two genotypes, as well as several thousand PAV sequences that are present in B73 but not Mo17. Haplotype-specific PAVs contain hundreds of single-copy, expressed genes that may contribute to heterosis and to the extraordinary phenotypic diversity of this important crop.

Deep sampling of the Palomero maize transcriptome by a high throughput strategy of pyrosequencing

BACKGROUND

In-depth sequencing analysis has not been able to determine the overall complexity of transcriptional activity of a plant organ or tissue sample. In some cases, deep parallel sequencing of Expressed Sequence Tags (ESTs), although not yet optimized for the sequencing of cDNAs, has represented an efficient procedure for validating gene prediction and estimating overall gene coverage. This approach could be very valuable for complex plant genomes. In addition, little emphasis has been given to efforts aiming at an estimation of the overall transcriptional universe found in a multicellular organism at a specific developmental stage.

RESULTS

To explore, in depth, the transcriptional diversity in an ancient maize landrace, we developed a protocol to optimize the sequencing of cDNAs and performed 4 consecutive GS20-454 pyrosequencing runs of a cDNA library obtained from 2 week-old Palomero Toluqueño maize plants. The protocol reported here allowed obtaining over 90% of informative sequences. These GS20-454 runs generated over 1.5 Million reads, representing the largest amount of sequences reported from a single plant cDNA library. A collection of 367,391 quality-filtered reads (30.09 Mb) from a single run was sufficient to identify transcripts corresponding to 34% of public maize ESTs databases; total sequences generated after 4 filtered runs increased this coverage to 50%. Comparisons of all 1.5 Million reads to the Maize Assembled Genomic Islands (MAGIs) provided evidence for the transcriptional activity of 11% of MAGIs. We estimate that 5.67% (86,069 sequences) do not align with public ESTs or annotated genes, potentially representing new maize transcripts. Following the assembly of 74.4% of the reads in 65,493 contigs, real-time PCR of selected genes confirmed a predicted correlation between the abundance of GS20-454 sequences and corresponding levels of gene expression.

CONCLUSION

A protocol was developed that significantly increases the number, length and quality of cDNA reads using massive 454 parallel sequencing. We show that recurrent 454 pyrosequencing of a single cDNA sample is necessary to attain a thorough representation of the transcriptional universe present in maize, that can also be used to estimate transcript abundance of specific genes. This data suggests that the molecular and functional diversity contained in the vast native landraces remains to be explored, and that large-scale transcriptional sequencing of a presumed ancestor of the modern maize varieties represents a valuable approach to characterize the functional diversity of maize for future agricultural and evolutionary studies.

Transposons that clean up after themselves

A transposon in the germline genome of the ciliate Oxytricha uses its transposase to remove itself, as well as other germline-limited DNA, from the differentiating somatic genome during development.

Too many ends: aberrant transposition

A recent study by Zhang and colleagues published in the March 15, 2009, issue of Genes & Development (pp. 755-765) demonstrates that maize Ac/Ds transposons mediate translocations and other rearrangements through aberrant execution of the normal transposition process. Ac transposase uses one end from each of two neighboring elements in these events, which may happen more commonly than previously thought. In genomes where there can be many transposon ends scattered across all the chromosomes, such mistakes can have important consequences.

Genome-wide and organ-specific landscapes of epigenetic modifications and their relationships to mRNA and small RNA transcriptomes in maize

Maize (Zea mays) has an exceptionally complex genome with a rich history in both epigenetics and evolution. We report genomic landscapes of representative epigenetic modifications and their relationships to mRNA and small RNA (smRNA) transcriptomes in maize shoots and roots. The epigenetic patterns differed dramatically between genes and transposable elements, and two repressive marks (H3K27me3 and DNA methylation) were usually mutually exclusive. We found an organ-specific distribution of canonical microRNAs (miRNAs) and endogenous small interfering RNAs (siRNAs), indicative of their tissue-specific biogenesis. Furthermore, we observed that a decreasing level of mop1 led to a concomitant decrease of 24-nucleotide siRNAs relative to 21-nucleotide miRNAs in a tissue-specific manner. A group of 22-nucleotide siRNAs may originate from long-hairpin double-stranded RNAs and preferentially target gene-coding regions. Additionally, a class of miRNA-like smRNAs, whose putative precursors can form short hairpins, potentially targets genes in trans. In summary, our data provide a critical analysis of the maize epigenome and its relationships to mRNA and smRNA transcriptomes.

Methylation-sensitive linking libraries enhance gene-enriched sequencing of complex genomes and map DNA methylation domains

Background

Many plant genomes are resistant to whole-genome assembly due to an abundance of repetitive sequence, leading to the development of gene-rich sequencing techniques. Two such techniques are hypomethylated partial restriction (HMPR) and methylation spanning linker libraries (MSLL). These libraries differ from other gene-rich datasets in having larger insert sizes, and the MSLL clones are designed to provide reads localized to "epigenetic boundaries" where methylation begins or ends.

Results

A large-scale study in maize generated 40,299 HMPR sequences and 80,723 MSLL sequences, including MSLL clones exceeding 100 kb. The paired end reads of MSLL and HMPR clones were shown to be effective in linking existing gene-rich sequences into scaffolds. In addition, it was shown that the MSLL clones can be used for anchoring these scaffolds to a BAC-based physical map. The MSLL end reads effectively identified epigenetic boundaries, as indicated by their preferential alignment to regions upstream and downstream from annotated genes. The ability to precisely map long stretches of fully methylated DNA sequence is a unique outcome of MSLL analysis, and was also shown to provide evidence for errors in gene identification. MSLL clones were observed to be significantly more repeat-rich in their interiors than in their end reads, confirming the correlation between methylation and retroelement content. Both MSLL and HMPR reads were found to be substantially gene-enriched, with the SalI MSLL libraries being the most highly enriched (31% align to an EST contig), while the HMPR clones exhibited exceptional depletion of repetitive DNA (to ~11%). These two techniques were compared with other gene-enrichment methods, and shown to be complementary.

Conclusion

MSLL technology provides an unparalleled approach for mapping the epigenetic status of repetitive blocks and for identifying sequences mis-identified as genes. Although the types and natures of epigenetic boundaries are barely understood at this time, MSLL technology flags both approximate boundaries and methylated genes that deserve additional investigation. MSLL and HMPR sequences provide a valuable resource for maize genome annotation, and are a uniquely valuable complement to any plant genome sequencing project. In order to make these results fully accessible to the community, a web display was developed that shows the alignment of MSLL, HMPR, and other gene-rich sequences to the BACs; this display is continually updated with the latest ESTs and BAC sequences.

The Pic19 NBS-LRR gene family members are closely linked to Scmv1, but not involved in maize resistance to sugarcane mosaic virus

Sugarcane mosaic virus (SCMV) is the causal pathogen for a severe mosaic virus disease of maize worldwide. In our previous research, the maize resistance gene analog (RGA) Pic19 and its three cognate BAC contigs were mapped to the same region as the SCMV resistance gene Scmv1. Here we report the isolation and characterization of the Pic19R gene family members from the inbred line FAP1360A, which shows complete resistance to SCMV. Two primer pairs were designed based on the conserved regions among the known Pic19 paralogs and used for rapid amplification of cDNA ends of FAP1360A. Six full-length cDNAs, corresponding to the Pic19R-1 to -6 paralogs, were obtained. Three of them (Pic19R-1 to -3) had uninterrupted coding sequences and were, therefore, regarded as candidates for the Scmv1 gene. A total of 18 positive BAC clones harboring the Pic19R-2 to -5 paralogs were obtained from the FAP1360A BAC library and assembled into two BAC contigs. Two markers, tagging Pic19R-2 and -3 and Pic19R-4, were developed and used to genotype a high-resolution mapping population segregating solely for the Scmv1 locus. Although closely linked, none of these three Pic19R paralogs co-segregated with the Scmv1 locus. Analysis of the Pic19R family indicated that the Pic19R-1 paralog is identical to the known Rxo1 gene conferring resistance to rice bacterial streak disease and none of the other Pic19R paralogs seems to be involved in resistance to SCMV.

Display of a maize cDNA library on baculovirus infected insect cells

BACKGROUND

Maize is a good model system for cereal crop genetics and development because of its rich genetic heritage and well-characterized morphology. The sequencing of its genome is well advanced, and new technologies for efficient proteomic analysis are needed. Baculovirus expression systems have been used for the last twenty years to express in insect cells a wide variety of eukaryotic proteins that require complex folding or extensive posttranslational modification. More recently, baculovirus display technologies based on the expression of foreign sequences on the surface of Autographa californica (AcMNPV) have been developed. We investigated the potential of a display methodology for a cDNA library of maize young seedlings.

RESULTS

We constructed a full-length cDNA library of young maize etiolated seedlings in the transfer vector pAcTMVSVG. The library contained a total of 2.5 x 10(5) independent clones. Expression of two known maize proteins, calreticulin and auxin binding protein (ABP1), was shown by western blot analysis of protein extracts from insect cells infected with the cDNA library. Display of the two proteins in infected insect cells was shown by selective biopanning using magnetic cell sorting and demonstrated proof of concept that the baculovirus maize cDNA display library could be used to identify and isolate proteins.

CONCLUSION

The maize cDNA library constructed in this study relies on the novel technology of baculovirus display and is unique in currently published cDNA libraries. Produced to demonstrate proof of principle, it opens the way for the development of a eukaryotic in vivo display tool which would be ideally suited for rapid screening of the maize proteome for binding partners, such as proteins involved in hormone regulation or defence.

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.

Gene expression analyses in maize inbreds and hybrids with varying levels of heterosis

BACKGROUND

Heterosis is the superior performance of F1 hybrid progeny relative to the parental phenotypes. Maize exhibits heterosis for a wide range of traits, however the magnitude of heterosis is highly variable depending on the choice of parents and the trait(s) measured. We have used expression profiling to determine whether the level, or types, of non-additive gene expression vary in maize hybrids with different levels of genetic diversity or heterosis.

RESULTS

We observed that the distributions of better parent heterosis among a series of 25 maize hybrids generally do not exhibit significant correlations between different traits. Expression profiling analyses for six of these hybrids, chosen to represent diversity in genotypes and heterosis responses, revealed a correlation between genetic diversity and transcriptional variation. The majority of differentially expressed genes in each of the six different hybrids exhibited additive expression patterns, and approximately 25% exhibited statistically significant non-additive expression profiles. Among the non-additive profiles, approximately 80% exhibited hybrid expression levels between the parental levels, approximately 20% exhibited hybrid expression levels at the parental levels and ~1% exhibited hybrid levels outside the parental range.

CONCLUSION

We have found that maize inbred genetic diversity is correlated with transcriptional variation. However, sampling of seedling tissues indicated that the frequencies of additive and non-additive expression patterns are very similar across a range of hybrid lines. These findings suggest that heterosis is probably not a consequence of higher levels of additive or non-additive expression, but may be related to transcriptional variation between parents. The lack of correlation between better parent heterosis levels for different traits suggests that transcriptional diversity at specific sets of genes may influence heterosis for different traits.

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.

Characterization of paralogous protein families in rice

BACKGROUND

High gene numbers in plant genomes reflect polyploidy and major gene duplication events. Oryza sativa, cultivated rice, is a diploid monocotyledonous species with a ~390 Mb genome that has undergone segmental duplication of a substantial portion of its genome. This, coupled with other genetic events such as tandem duplications, has resulted in a substantial number of its genes, and resulting proteins, occurring in paralogous families.

RESULTS

Using a computational pipeline that utilizes Pfam and novel protein domains, we characterized paralogous families in rice and compared these with paralogous families in the model dicotyledonous diploid species, Arabidopsis thaliana. Arabidopsis, which has undergone genome duplication as well, has a substantially smaller genome (~120 Mb) and gene complement compared to rice. Overall, 53% and 68% of the non-transposable element-related rice and Arabidopsis proteins could be classified into paralogous protein families, respectively. Singleton and paralogous family genes differed substantially in their likelihood of encoding a protein of known or putative function; 26% and 66% of singleton genes compared to 73% and 96% of the paralogous family genes encode a known or putative protein in rice and Arabidopsis, respectively. Furthermore, a major skew in the distribution of specific gene function was observed; a total of 17 Gene Ontology categories in both rice and Arabidopsis were statistically significant in their differential distribution between paralogous family and singleton proteins. In contrast to mammalian organisms, we found that duplicated genes in rice and Arabidopsis tend to have more alternative splice forms. Using data from Massively Parallel Signature Sequencing, we show that a significant portion of the duplicated genes in rice show divergent expression although a correlation between sequence divergence and correlation of expression could be seen in very young genes.

CONCLUSION

Collectively, these data suggest that while co-regulation and conserved function are present in some paralogous protein family members, evolutionary pressures have resulted in functional divergence with differential expression patterns.

Genome size and GC content evolution of Festuca: ancestral expansion and subsequent reduction

BACKGROUND AND AIMS

Plant evolution is well known to be frequently associated with remarkable changes in genome size and composition; however, the knowledge of long-term evolutionary dynamics of these processes still remains very limited. Here a study is made of the fine dynamics of quantitative genome evolution in Festuca (fescue), the largest genus in Poaceae (grasses).

METHODS

Using flow cytometry (PI, DAPI), measurements were made of DNA content (2C-value), monoploid genome size (Cx-value), average chromosome size (C/n-value) and cytosine + guanine (GC) content of 101 Festuca taxa and 14 of their close relatives. The results were compared with the existing phylogeny based on ITS and trnL-F sequences.

KEY RESULTS

The divergence of the fescue lineage from related Poeae was predated by about a 2-fold monoploid genome and chromosome size enlargement, and apparent GC content enrichment. The backward reduction of these parameters, running parallel in both main evolutionary lineages of fine-leaved and broad-leaved fescues, appears to diverge among the existing species groups. The most dramatic reductions are associated with the most recently and rapidly evolving groups which, in combination with recent intraspecific genome size variability, indicate that the reduction process is probably ongoing and evolutionarily young. This dynamics may be a consequence of GC-rich retrotransposon proliferation and removal. Polyploids derived from parents with a large genome size and high GC content (mostly allopolyploids) had smaller Cx- and C/n-values and only slightly deviated from parental GC content, whereas polyploids derived from parents with small genome and low GC content (mostly autopolyploids) generally had a markedly increased GC content and slightly higher Cx- and C/n-values.

CONCLUSIONS

The present study indicates the high potential of general quantitative characters of the genome for understanding the long-term processes of genome evolution, testing evolutionary hypotheses and their usefulness for large-scale genomic projects. Taken together, the results suggest that there is an evolutionary advantage for small genomes in Festuca.

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.

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.

Towards the molecular basis of heterosis

Heterosis describes the superior performance of heterozygous hybrid plants over their homozygous parental inbred lines. Despite the rediscovery of this phenomenon a century ago and its paramount agronomic importance, the genetic and molecular basis of heterosis remains enigmatic. Recently, various pioneer studies described differences in genome organization and gene expression of hybrids and their parental inbred lines. At the genomic level, a significant loss of colinearity at many loci between different inbred lines of maize was observed. At the level of gene expression, complex transcriptional networks specific for different developmental stages and tissues were monitored in maize (Zea mays), rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana). Integration of this complex expression data might contribute to improve our understanding of the molecular basis of heterosis.

The evolution of transposon repeat-induced point mutation in the genome of Colletotrichum cereale: reconciling sex, recombination and homoplasy in an ''asexual" pathogen

Mobile transposable elements are among the primary drivers of the evolution of eukaryotic genomes. For fungi, repeat-induced point mutation (RIP) silencing minimizes deleterious effects of transposons by mutating multicopy DNA during meiosis. In this study we identify five transposon species from the mitosporic fungus Colletotrichum cereale and report the signature pattern of RIP acting in a lineage-specific manner on 21 of 35 unique transposon copies, providing the first evidence for sexual recombination for this species. Sequence analysis of genomic populations of the retrotransposon Ccret2 showed repeated rounds of RIP mutation acting on different copies of the element. In the RIPped Ccret2 population, there were multiple inferences of incongruence primarily attributed to RIP-induced homoplasy. This study supports the view that the sequence variability of transposon populations in filamentous fungi reflects the activities of evolutionary processes that fall outside of typical phylogenetic or population genetic reconstructions.

Allele-specific expression patterns reveal biases and embryo-specific parent-of-origin effects in hybrid maize

We employed allele-specific expression (ASE) analyses to document biased allelic expression in maize (Zea mays). A set of 316 quantitative ASE assays were used to profile the relative allelic expression in seedling tissue derived from five maize hybrids. The different hybrids included in this study exhibit a range of heterosis levels; however, we did not observe differences in the frequencies of allelic bias. Allelic biases in gene expression were consistently observed for approximately 50% of the genes assayed in hybrid seedlings. The relative proportion of genes that exhibit cis- or trans-acting regulatory variation was very similar among the different genotypes. The cis-acting regulatory variation was more prevalent and resulted in greater expression differences than trans-acting regulatory variation for these genes. The ASE assays were further used to compare the relative expression of the B73 and Mo17 alleles in three tissue types (seedling, immature ear, and embryo) derived from reciprocal hybrids. These comparisons provided evidence for tissue-specific cis-acting variation and for a slight maternal expression bias in approximately 20% of genes in embryo tissue. Collectively, these data provide evidence for prevalent cis-acting regulatory variation that contributes to biased allelic expression between genotypes and between tissues.

Mapping of genome-wide resistance gene analogs (RGAs) in maize (Zea mays L.)

Isolation and mapping of genome-wide resistance (R) gene analogs (RGAs) is of importance in identifying candidate(s) for a particular resistance gene/QTL. Here we reported our result in mapping totally 228 genome-wide RGAs in maize. By developing RGA-tagged markers and subsequent genotyping a population consisting of 294 recombinant inbred lines (RILs), 67 RGAs were genetically mapped on maize genome. Meanwhile, in silico mapping was conducted to anchor 113 RGAs by comparing all 228 RGAs to those anchored EST and BAC/BAC-end sequences via tblastx search (E-value < 10(-20)). All RGAs from different mapping efforts were integrated into the existing SSR linkage map. After accounting for redundancy, the resultant RGA linkage map was composed of 153 RGAs that were mapped onto 172 loci on maize genome, and the mapped RGAs accounted for approximate three quarters of the genome-wide RGAs in maize. The extensive co-localizations were observed between mapped RGAs and resistance gene/QTL loci, implying the usefulness of this RGA linkage map in R gene cloning via candidate gene approach.

Highly asymmetric rice genomes

Background

Individuals in the same species are assumed to share the same genomic set. However, it is not unusual to find an orthologous gene only in small subset of the species, and recent genomic studies suggest that structural rearrangements are very frequent between genomes in the same species. Two recently sequenced rice genomes Oryza sativa L. var. Nipponbare and O. sativa L. var. 93-11 provide an opportunity to systematically investigate the extent of the gene repertoire polymorphism, even though the genomic data of 93-11 derived from whole-short-gun sequencing is not yet as complete as that of Nipponbare.

Results

We compared gene contents and the genomic locations between two rice genomes. Our conservative estimates suggest that at least 10% of the genes in the genomes were either under presence/absence polymorphism (5.2%) or asymmetrically located between genomes (4.7%). The proportion of these "asymmetric genes" varied largely among gene groups, in which disease resistance (R) genes and the RLK kinase gene group had 11.6 and 7.8 times higher proportion of asymmetric genes than housekeeping genes (Myb and MADS). The significant difference in the proportion of asymmetric genes among gene groups suggests that natural selection is responsible for maintaining genomic asymmetry. On the other hand, the nucleotide diversity in 17 R genes under presence/absence polymorphism was generally low (average nucleotide diversity = 0.0051).

Conclusion

The genomic symmetry was disrupted by 10% of asymmetric genes, which could cause genetic variation through more unequal crossing over, because these genes had no allelic counterparts to pair and then they were free to pair with homologues at non-allelic loci, during meiosis in heterozygotes. It might be a consequence of diversifying selection that increased the structural divergence among genomes, and of purifying selection that decreased nucleotide divergence in each R gene locus.