Intraspecific violation of genetic colinearity and its implications in maize

Construction of high-quality recombination maps with low-coverage genomic sequencing for joint linkage analysis in maize

BACKGROUND

A genome-wide association study (GWAS) is the foremost strategy used for finding genes that control human diseases and agriculturally important traits, but it often reports false positives. In contrast, its complementary method, linkage analysis, provides direct genetic confirmation, but with limited resolution. A joint approach, using multiple linkage populations, dramatically improves resolution and statistical power. For example, this approach has been used to confirm that many complex traits, such as flowering time controlling adaptation in maize, are controlled by multiple genes with small effects. In addition, genotyping by sequencing (GBS) at low coverage not only produces genotyping errors, but also results in large datasets, making the use of high-throughput sequencing technologies computationally inefficient or unfeasible.

RESULTS

In this study, we converted raw SNPs into effective recombination bins. The reduced bins not only retain the original information, but also correct sequencing errors from low-coverage genomic sequencing. To further increase the statistical power and resolution, we merged a new temperate maize nested association mapping (NAM) population derived in China (CN-NAM) with the existing maize NAM population developed in the US (US-NAM). Together, the two populations contain 36 families and 7,000 recombinant inbred lines (RILs). One million SNPs were generated for all the RILs with GBS at low coverage. We developed high-quality recombination maps for each NAM population to correct genotyping errors and improve the computational efficiency of the joint linkage analysis. The original one million SNPs were reduced to 4,932 and 5,296 recombination bins with average interval distances of 0.34 cM and 0.28 cM for CN-NAM and US-NAM, respectively. The quantitative trait locus (QTL) mapping for flowering time (days to tasseling) indicated that the high-density, recombination bin map improved resolution of QTL mapping by 50 % compared with that using a medium-density map. We also demonstrated that combining the CN-NAM and US-NAM populations improves the power to detect QTL by 50 % compared to single NAM population mapping. Among the QTLs mapped by joint usage of the US-NAM and CN-NAM maps, 25 % of the QTLs overlapped with known flowering-time genes in maize.

CONCLUSION

This study provides directions and resources for the research community, especially maize researchers, for future studies using the recombination bin strategy for joint linkage analysis. Available resources include efficient usage of low-coverage genomic sequencing, detailed positions for genes controlling maize flowering, and recombination bin maps and flowering- time data for both CN and US NAMs. Maize researchers even have the opportunity to grow both CN and US NAM populations to study the traits of their interest, as the seeds of both NAM populations are available from the seed repository in China and the US.

High-resolution genetic mapping of maize pan-genome sequence anchors

In addition to single-nucleotide polymorphisms, structural variation is abundant in many plant genomes. The structural variation across a species can be represented by a 'pan-genome', which is essential to fully understand the genetic control of phenotypes. However, the pan-genome's complexity hinders its accurate assembly via sequence alignment. Here we demonstrate an approach to facilitate pan-genome construction in maize. By performing 18 trillion association tests we map 26 million tags generated by reduced representation sequencing of 14,129 maize inbred lines. Using machine-learning models we select 4.4 million accurately mapped tags as sequence anchors, 1.1 million of which are presence/absence variations. Structural variations exhibit enriched association with phenotypic traits, indicating that it is a significant source of adaptive variation in maize. The ability to efficiently map ultrahigh-density pan-genome sequence anchors enables fine characterization of structural variation and will advance both genetic research and breeding in many crops.

New insight into a complex plant-fungal pathogen interaction

The coevolution of plants and microbes has shaped plant mechanisms that detect and repel pathogens. A newly identified plant gene confers partial resistance to a fungal pathogen not by preventing initial infection but by limiting its spread through the plant.

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.

Recombination in diverse maize is stable, predictable, and associated with genetic load

Among the fundamental evolutionary forces, recombination arguably has the largest impact on the practical work of plant breeders. Varying over 1,000-fold across the maize genome, the local meiotic recombination rate limits the resolving power of quantitative trait mapping and the precision of favorable allele introgression. The consequences of low recombination also theoretically extend to the species-wide scale by decreasing the power of selection relative to genetic drift, and thereby hindering the purging of deleterious mutations. In this study, we used genotyping-by-sequencing (GBS) to identify 136,000 recombination breakpoints at high resolution within US and Chinese maize nested association mapping populations. We find that the pattern of cross-overs is highly predictable on the broad scale, following the distribution of gene density and CpG methylation. Several large inversions also suppress recombination in distinct regions of several families. We also identify recombination hotspots ranging in size from 1 kb to 30 kb. We find these hotspots to be historically stable and, compared with similar regions with low recombination, to have strongly differentiated patterns of DNA methylation and GC content. We also provide evidence for the historical action of GC-biased gene conversion in recombination hotspots. Finally, using genomic evolutionary rate profiling (GERP) to identify putative deleterious polymorphisms, we find evidence for reduced genetic load in hotspot regions, a phenomenon that may have considerable practical importance for breeding programs worldwide.

Transposable elements contribute to activation of maize genes in response to abiotic stress

Transposable elements (TEs) account for a large portion of the genome in many eukaryotic species. Despite their reputation as “junk” DNA or genomic parasites deleterious for the host, TEs have complex interactions with host genes and the potential to contribute to regulatory variation in gene expression. It has been hypothesized that TEs and genes they insert near may be transcriptionally activated in response to stress conditions. The maize genome, with many different types of TEs interspersed with genes, provides an ideal system to study the genome-wide influence of TEs on gene regulation. To analyze the magnitude of the TE effect on gene expression response to environmental changes, we profiled gene and TE transcript levels in maize seedlings exposed to a number of abiotic stresses. Many genes exhibit up- or down-regulation in response to these stress conditions. The analysis of TE families inserted within upstream regions of up-regulated genes revealed that between four and nine different TE families are associated with up-regulated gene expression in each of these stress conditions, affecting up to 20% of the genes up-regulated in response to abiotic stress, and as many as 33% of genes that are only expressed in response to stress. Expression of many of these same TE families also responds to the same stress conditions. The analysis of the stress-induced transcripts and proximity of the transposon to the gene suggests that these TEs may provide local enhancer activities that stimulate stress-responsive gene expression. Our data on allelic variation for insertions of several of these TEs show strong correlation between the presence of TE insertions and stress-responsive up-regulation of gene expression. Our findings suggest that TEs provide an important source of allelic regulatory variation in gene response to abiotic stress in maize.

Transposable elements are mobile DNA elements that are a prevalent component of many eukaryotic genomes. While transposable elements can often have deleterious effects through insertions into protein-coding genes they may also contribute to regulatory variation of gene expression. There are a handful of examples in which specific transposon insertions contribute to regulatory variation of nearby genes, particularly in response to environmental stress. We sought to understand the genome-wide influence of transposable elements on gene expression responses to abiotic stress in maize, a plant with many families of transposable elements located in between genes. Our analysis suggests that a small number of maize transposable element families may contribute to the response of nearby genes to abiotic stress by providing stress-responsive enhancer-like functions. The specific insertions of transposable elements are often polymorphic within a species. Our data demonstrate that allelic variation for insertions of the transposable elements associated with stress-responsive expression can contribute to variation in the regulation of nearby genes. Thus novel insertions of transposable elements provide a potential mechanism for genes to acquire cis-regulatory influences that could contribute to heritable variation for stress response.

Comparative proteomic analysis of a membrane-enriched fraction from flag leaves reveals responses to chemical hybridization agent SQ-1 in wheat

The induction of wheat male fertile lines by using the chemical hybridizing agent SQ-1 (CHA-SQ-1) is an effective approach in the utilization of heterosis; however, the molecular basis of male fertility remains unknown. Wheat flag leaves are the initial receptors of CHA-SQ-1 and their membrane structure plays a vital role in response to CHA-SQ-1 stress. To investigate the response of wheat flag leaves to CHA-SQ-1 stress, we compared their quantitative proteomic profiles in the absence and presence of CHA-SQ-1. Our results indicated that wheat flag leaves suffered oxidative stress during CHA-SQ-1 treatments. Leaf O2 (-), H2O2, and malonaldehyde levels were significantly increased within 10 h after CHA-SQ-1 treatment, while the activities of major antioxidant enzymes such as superoxide dismutase, catalase, and guaiacol peroxidase were significantly reduced. Proteome profiles of membrane-enriched fraction showed a change in the abundance of a battery of membrane proteins involved in multiple biological processes. These variable proteins mainly impaired photosynthesis, ATP synthesis protein mechanisms and were involved in the response to stress. These results provide an explanation of the relationships between membrane proteomes and anther abortion and the practical application of CHA for hybrid breeding.

Genome evolution in maize: from genomes back to genes

Maize occupies dual roles as both (a) one of the big-three grain species (along with rice and wheat) responsible for providing more than half of the calories consumed around the world, and (b) a model system for plant genetics and cytogenetics dating back to the origin of the field of genetics in the early twentieth century. The long history of genetic investigation in this species combined with modern genomic and quantitative genetic data has provided particular insight into the characteristics of genes linked to phenotypes and how these genes differ from many other sequences in plant genomes that are not easily distinguishable based on molecular data alone. These recent results suggest that the number of genes in plants that make significant contributions to phenotype may be lower than the number of genes defined by current molecular criteria, and also indicate that syntenic conservation has been underemphasized as a marker for gene function.

Investigating the molecular genetic basis of heterosis for internode expansion in maize by microRNA transcriptomic deep sequencing

Heterosis has been used widely in the breeding of maize and other crops and plays an important role in increasing yield, improving quality, and enhancing stress resistance, but its molecular mechanism is far from clear. To determine whether microRNA (miRNA)-dependent gene regulation is responsible for heterosis of elongating internodes below the ear and ear height in maize, a deep-sequencing strategy was applied to the elite hybrid Xundan20, which is currently cultivated widely in China, and its two parents. RNA was extracted from the eighth internode because it shows clear internode length heterosis. A total of 99 conserved maize miRNAs were detected in both the hybrid and parental lines. Most of these miRNAs were expressed nonadditively in the hybrid compared with its parental lines. These results indicated that miRNAs might participate in heterosis during internode expansion in maize and exert an influence on ear and plant height via the repression of their target genes. In total, eight novel miRNAs belonging to four miRNA families were predicted in the expanding internode. Global repression of miRNAs in the hybrid, which might result in enhanced gene expression, might be one reason why the hybrid shows longer internodes and taller seedlings compared with its parental lines.

Polarized gene conversion at the bz locus of maize

Nucleotide diversity is greater in maize than in most organisms studied to date, so allelic pairs in a hybrid tend to be highly polymorphic. Most recombination events between such pairs of maize polymorphic alleles are crossovers. However, intragenic recombination events not associated with flanking marker exchange, corresponding to noncrossover gene conversions, predominate between alleles derived from the same progenitor. In these dimorphic heterozygotes, the two alleles differ only at the two mutant sites between which recombination is being measured. To investigate whether gene conversion at the bz locus is polarized, two large diallel crossing matrices involving mutant sites spread across the bz gene were performed and more than 2,500 intragenic recombinants were scored. In both diallels, around 90% of recombinants could be accounted for by gene conversion. Furthermore, conversion exhibited a striking polarity, with sites located within 150 bp of the start and stop codons converting more frequently than sites located in the middle of the gene. The implications of these findings are discussed with reference to recent data from genome-wide studies in other plants.

Heterosis in early maize ear inflorescence development: a genome-wide transcription analysis for two maize inbred lines and their hybrid

Heterosis, or hybrid vigor, contributes to superior agronomic performance of hybrids compared to their inbred parents. Despite its importance, little is known about the genetic and molecular basis of heterosis. Early maize ear inflorescences formation affects grain yield, and are thus an excellent model for molecular mechanisms involved in heterosis. To determine the parental contributions and their regulation during maize ear-development-genesis, we analyzed genome-wide digital gene expression profiles in two maize elite inbred lines (B73 and Mo17) and their F₁ hybrid using deep sequencing technology. Our analysis revealed 17,128 genes expressed in these three genotypes and 22,789 genes expressed collectively in the present study. Approximately 38% of the genes were differentially expressed in early maize ear inflorescences from heterotic cross, including many transcription factor genes and some presence/absence variations (PAVs) genes, and exhibited multiple modes of gene action. These different genes showing differential expression patterns were mainly enriched in five cellular component categories (organelle, cell, cell part, organelle part and macromolecular complex), five molecular function categories (structural molecule activity, binding, transporter activity, nucleic acid binding transcription factor activity and catalytic activity), and eight biological process categories (cellular process, metabolic process, biological regulation, regulation of biological process, establishment of localization, cellular component organization or biogenesis, response to stimulus and localization). Additionally, a significant number of genes were expressed in only one inbred line or absent in both inbred lines. Comparison of the differences of modes of gene action between previous studies and the present study revealed only a small number of different genes had the same modes of gene action in both maize seedlings and ear inflorescences. This might be an indication that in different tissues or developmental stages, different global expression patterns prevail, which might nevertheless be related to heterosis. Our results support the hypotheses that multiple molecular mechanisms (dominance and overdominance modes) contribute to heterosis.

Variation in allelic expression associated with a recombination hotspot in Zea mays

Gene expression is a complex process, requiring precise spatial and temporal regulation of transcription factor activity; however, modifications of individual cis- and trans-acting modules can be molded by natural selection to create a sizeable number of novel phenotypes. Results from decades of research indicate that developmental and phenotypic divergence among eukaryotic organisms is driven primarily by variation in levels of gene expression that are dictated by mutations, either in structural or regulatory regions, of genes. The relative contributions and interplay of cis- and trans-acting regulatory factors to this evolutionary process, however, remain poorly understood. Analysis of eight genes in the Bz1-Sh1 interval of Zea mays (maize) indicates significant allele-specific expression biases in at least one tissue for all genes, ranging from 1.3-fold to 36-fold. All detected effects were cis-regulatory in nature, although genetic background may also influence the level of expression bias and tissue specificity for some allelic combinations. Most allelic pairs exhibited the same direction and approximate intensity of bias across all four tissues; however, a subset of allelic pairs show alternating dominance across different tissue types or variation in the degree of bias in different tissues. In addition, the genes showing the most striking levels of allelic bias co-localize with a previously described recombination hotspot in this region, suggesting a naturally occurring genetic mechanism for creating regulatory variability for a subset of plant genes that may ultimately lead to evolutionary diversification.

HelitronScanner uncovers a large overlooked cache of Helitron transposons in many plant genomes

Transposons make up the bulk of eukaryotic genomes, but are difficult to annotate because they evolve rapidly. Most of the unannotated portion of sequenced genomes is probably made up of various divergent transposons that have yet to be categorized. Helitrons are unusual rolling circle eukaryotic transposons that often capture gene sequences, making them of considerable evolutionary importance. Unlike other DNA transposons, Helitrons do not end in inverted repeats or create target site duplications, so they are particularly challenging to identify. Here we present HelitronScanner, a two-layered local combinational variable (LCV) tool for generalized Helitron identification that represents a major improvement over previous identification programs based on DNA sequence or structure. HelitronScanner identified 64,654 Helitrons from a wide range of plant genomes in a highly automated way. We tested HelitronScanner's predictive ability in maize, a species with highly heterogeneous Helitron elements. LCV scores for the 5' and 3' termini of the predicted Helitrons provide a primary confidence level and element copy number provides a secondary one. Newly identified Helitrons were validated by PCR assays or by in silico comparative analysis of insertion site polymorphism among multiple accessions. Many new Helitrons were identified in model species, such as maize, rice, and Arabidopsis, and in a variety of organisms where Helitrons had not been reported previously to our knowledge, leading to a major upward reassessment of their abundance in plant genomes. HelitronScanner promises to be a valuable tool in future comparative and evolutionary studies of this major transposon superfamily.

The control of recombination in wheat by Ph1 and its use in breeding

Two meiotic processes have a major influence on the plant breeding, namely, the independent assortment of chromosomes, and recombination. The major chromosome pairing locus in hexaploid and tetraploid wheat, Ph1, has a significant effect on both these processes. This chapter reviews our current understanding of this locus and how mutants of it can be exploited for breeding purposes.

Inheritance patterns and stability of DNA methylation variation in maize near-isogenic lines

DNA methylation is a chromatin modification that contributes to epigenetic regulation of gene expression. The inheritance patterns and trans-generational stability of 962 differentially methylated regions (DMRs) were assessed in a panel of 71 near-isogenic lines (NILs) derived from maize (Zea mays) inbred lines B73 and Mo17. The majority of DMRs exhibit inheritance patterns that would be expected for local (cis) inheritance of DNA methylation variation such that DNA methylation level was coupled to local genotype. There are few examples of DNA methylation that exhibit trans-acting control or paramutation-like patterns. The cis-inherited DMRs provide an opportunity to study the stability of inheritance for DNA methylation variation. There was very little evidence for alterations of DNA methylation levels at these DMRs during the generations of the NIL population development. DNA methylation level was associated with local genotypes in nearly all of the >30,000 potential cases of inheritance. The majority of the DMRs were not associated with small RNAs. Together, our results suggest that a significant portion of DNA methylation variation in maize exhibits locally (cis) inherited patterns, is highly stable, and does not require active programming by small RNAs for maintenance. DNA methylation may contribute to heritable epigenetic information in many eukaryotic genomes. In this study, we have documented the inheritance patterns and trans-generational stability for nearly 1000 DNA methylation variants in a segregating maize population. At most loci studied, the DNA methylation differences are locally inherited and are not influenced by the other allele or other genomic regions. The inheritance of DNA methylation levels across generations is quite robust with almost no examples of unstable inheritance, suggesting that DNA methylation differences can be quite stably inherited, even in segregating populations.

Genomic prediction in maize breeding populations with genotyping-by-sequencing

Genotyping-by-sequencing (GBS) technologies have proven capacity for delivering large numbers of marker genotypes with potentially less ascertainment bias than standard single nucleotide polymorphism (SNP) arrays. Therefore, GBS has become an attractive alternative technology for genomic selection. However, the use of GBS data poses important challenges, and the accuracy of genomic prediction using GBS is currently undergoing investigation in several crops, including maize, wheat, and cassava. The main objective of this study was to evaluate various methods for incorporating GBS information and compare them with pedigree models for predicting genetic values of lines from two maize populations evaluated for different traits measured in different environments (experiments 1 and 2). Given that GBS data come with a large percentage of uncalled genotypes, we evaluated methods using nonimputed, imputed, and GBS-inferred haplotypes of different lengths (short or long). GBS and pedigree data were incorporated into statistical models using either the genomic best linear unbiased predictors (GBLUP) or the reproducing kernel Hilbert spaces (RKHS) regressions, and prediction accuracy was quantified using cross-validation methods. The following results were found: relative to pedigree or marker-only models, there were consistent gains in prediction accuracy by combining pedigree and GBS data; there was increased predictive ability when using imputed or nonimputed GBS data over inferred haplotype in experiment 1, or nonimputed GBS and information-based imputed short and long haplotypes, as compared to the other methods in experiment 2; the level of prediction accuracy achieved using GBS data in experiment 2 is comparable to those reported by previous authors who analyzed this data set using SNP arrays; and GBLUP and RKHS models with pedigree with nonimputed and imputed GBS data provided the best prediction correlations for the three traits in experiment 1, whereas for experiment 2 RKHS provided slightly better prediction than GBLUP for drought-stressed environments, and both models provided similar predictions in well-watered environments.

InsertionMapper: a pipeline tool for the identification of targeted sequences from multidimensional high throughput sequencing data

Background

The advent of next-generation high-throughput technologies has revolutionized whole genome sequencing, yet some experiments require sequencing only of targeted regions of the genome from a very large number of samples. These regions can be amplified by PCR and sequenced by next-generation methods using a multidimensional pooling strategy. However, there is at present no available generalized tool for the computational analysis of target-enriched NGS data from multidimensional pools.

Results

Here we present InsertionMapper, a pipeline tool for the identification of targeted sequences from multidimensional high throughput sequencing data. InsertionMapper consists of four independently working modules: Data Preprocessing, Database Modeling, Dimension Deconvolution and Element Mapping. We illustrate InsertionMapper with an example from our project 'New reverse genetics resources for maize’, which aims to sequence-index a collection of 15,000 independent insertion sites of the transposon Ds in maize. Identified sequences are validated by PCR assays. This pipeline tool is applicable to similar scenarios requiring analysis of the tremendous output of short reads produced in NGS sequencing experiments of targeted genome sequences.

Conclusions

InsertionMapper is proven efficacious for the identification of target-enriched sequences from multidimensional high throughput sequencing data. With adjustable parameters and experiment configurations, this tool can save great computational effort to biologists interested in identifying their sequences of interest within the huge output of modern DNA sequencers. InsertionMapper is freely accessible at https://sourceforge.net/p/insertionmapper and http://bo.csam.montclair.edu/du/insertionmapper.

Genomic resources for gene discovery, functional genome annotation, and evolutionary studies of maize and its close relatives

Maize is one of the most important food crops and a key model for genetics and developmental biology. A genetically anchored and high-quality draft genome sequence of maize inbred B73 has been obtained to serve as a reference sequence. To facilitate evolutionary studies in maize and its close relatives, much like the Oryza Map Alignment Project (OMAP) (www.OMAP.org) bacterial artificial chromosome (BAC) resource did for the rice community, we constructed BAC libraries for maize inbred lines Zheng58, Chang7-2, and Mo17 and maize wild relatives Zea mays ssp. parviglumis and Tripsacum dactyloides. Furthermore, to extend functional genomic studies to maize and sorghum, we also constructed binary BAC (BIBAC) libraries for the maize inbred B73 and the sorghum landrace Nengsi-1. The BAC/BIBAC vectors facilitate transfer of large intact DNA inserts from BAC clones to the BIBAC vector and functional complementation of large DNA fragments. These seven Zea Map Alignment Project (ZMAP) BAC/BIBAC libraries have average insert sizes ranging from 92 to 148 kb, organellar DNA from 0.17 to 2.3%, empty vector rates between 0.35 and 5.56%, and genome equivalents of 4.7- to 8.4-fold. The usefulness of the Parviglumis and Tripsacum BAC libraries was demonstrated by mapping clones to the reference genome. Novel genes and alleles present in these ZMAP libraries can now be used for functional complementation studies and positional or homology-based cloning of genes for translational genomics.

Genetics and consequences of crop domestication

Phenotypic variation has been manipulated by humans during crop domestication, which occurred primarily between 3000 and 10000 years ago in the various centers of origin around the world. The process of domestication has profound consequences on crops, where the domesticate has moderately reduced genetic diversity relative to the wild ancestor across the genome, and severely reduced diversity for genes targeted by domestication. The question that remains is whether reduction in genetic diversity has affected crop production today. A case study in maize ( Zea mays ) demonstrates the application of understanding relationships between genetic diversity and phenotypic diversity in the wild ancestor and the domesticate. As an outcrossing species, maize has tremendous genetic variation. The complementary combination of genome-wide association mapping (GWAS) approaches, large HapMap data sets, and germplasm resources is leading to important discoveries of the relationship between genetic diversity and phenotypic variation and the impact of domestication on trait variation.

Comprehensive genotyping of the USA national maize inbred seed bank

BACKGROUND

Genotyping by sequencing, a new low-cost, high-throughput sequencing technology was used to genotype 2,815 maize inbred accessions, preserved mostly at the National Plant Germplasm System in the USA. The collection includes inbred lines from breeding programs all over the world.

RESULTS

The method produced 681,257 single-nucleotide polymorphism (SNP) markers distributed across the entire genome, with the ability to detect rare alleles at high confidence levels. More than half of the SNPs in the collection are rare. Although most rare alleles have been incorporated into public temperate breeding programs, only a modest amount of the available diversity is present in the commercial germplasm. Analysis of genetic distances shows population stratification, including a small number of large clusters centered on key lines. Nevertheless, an average fixation index of 0.06 indicates moderate differentiation between the three major maize subpopulations. Linkage disequilibrium (LD) decays very rapidly, but the extent of LD is highly dependent on the particular group of germplasm and region of the genome. The utility of these data for performing genome-wide association studies was tested with two simply inherited traits and one complex trait. We identified trait associations at SNPs very close to known candidate genes for kernel color, sweet corn, and flowering time; however, results suggest that more SNPs are needed to better explore the genetic architecture of complex traits.

CONCLUSIONS

The genotypic information described here allows this publicly available panel to be exploited by researchers facing the challenges of sustainable agriculture through better knowledge of the nature of genetic diversity.

Comparative proteomic analysis of embryos between a maize hybrid and its parental lines during early stages of seed germination

In spite of commercial use of heterosis in agriculture, the molecular basis of heterosis is poorly understood. It was observed that maize hybrid Zong3/87-1 exhibited an earlier onset or heterosis in radicle emergence. To get insights into the underlying mechanism of heterosis in radicle emergence, differential proteomic analysis between hybrid and its parental lines was performed. In total, the number of differentially expressed protein spots between hybrid and its parental lines in dry and 24 h imbibed seed embryos were 134 and 191, respectively, among which 47.01% (63/134) and 34.55% (66/191) protein spots displayed nonadditively expressed pattern. Remarkably, 54.55% of nonadditively accumulated proteins in 24 h imbibed seed embryos displayed above or equal to the level of the higher parent patterns. Moreover, 155 differentially expressed protein spots were identified, which were grouped into eight functional classes, including transcription & translation, energy & metabolism, signal transduction, disease & defense, storage protein, transposable element, cell growth & division and unclassified proteins. In addition, one of the upregulated proteins in F₁ hybrids was ZmACT2, a homolog of Arabidopsis thaliana ACT7 (AtACT7). Expressing ZmACT2 driven by the AtACT7 promoter partially complemented the low germination phenotype in the Atact7 mutant. These results indicated that hybridization between two parental lines can cause changes in the expression of a variety of proteins, and it is concluded that the altered pattern of gene expression at translational level in the hybrid may be responsible for the observed heterosis.

Aluminum tolerance in maize is associated with higher MATE1 gene copy number

Genome structure variation, including copy number variation and presence/absence variation, comprises a large extent of maize genetic diversity; however, its effect on phenotypes remains largely unexplored. Here, we describe how copy number variation underlies a rare allele that contributes to maize aluminum (Al) tolerance. Al toxicity is the primary limitation for crop production on acid soils, which make up 50% of the world's potentially arable lands. In a recombinant inbred line mapping population, copy number variation of the Al tolerance gene multidrug and toxic compound extrusion 1 (MATE1) is the basis for the quantitative trait locus of largest effect on phenotypic variation. This expansion in MATE1 copy number is associated with higher MATE1 expression, which in turn results in superior Al tolerance. The three MATE1 copies are identical and are part of a tandem triplication. Only three maize inbred lines carrying the three-copy allele were identified from maize and teosinte diversity panels, indicating that copy number variation for MATE1 is a rare, and quite likely recent, event. These maize lines with higher MATE1 copy number are also Al-tolerant, have high MATE1 expression, and originate from regions of highly acidic soils. Our findings show a role for copy number variation in the adaptation of maize to acidic soils in the tropics and suggest that genome structural changes may be a rapid evolutionary response to new environments.

Characterization of genetic diversity and linkage disequilibrium of ZmLOX4 and ZmLOX5 loci in maize

Maize (Zea mays L.) lipoxygenases (ZmLOXs) are well recognized as important players in plant defense against pathogens, especially in cross kingdom lipid communication with pathogenic fungi. This study is among the first to investigate genetic diversity at important gene paralogs ZmLOX4 and ZmLOX5. Sequencing of these genes in 400 diverse maize lines showed little genetic diversity and low linkage disequilibrium in the two genes. Importantly, we identified one inbred line in which ZmLOX5 has a disrupted open reading frame, a line missing ZmLOX5, and five lines with a duplication of ZmLOX5. Tajima's D test suggests that both ZmLOX4 and ZmLOX5 have been under neutral selection. Further investigation of haplotype data revealed that within the ZmLOX family members only ZmLOX12, a monocot specific ZmLOX, showed strong linkage disequilibrium that extends further than expected in maize. Linkage disequilibrium patterns at these loci of interest are crucial for future candidate gene association mapping studies. ZmLOX4 and ZmLOX5 mutations and copy number variants are under further investigation for crop improvement.

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

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

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

Mendelian and non-Mendelian regulation of gene expression in maize

Transcriptome variation plays an important role in affecting the phenotype of an organism. However, an understanding of the underlying mechanisms regulating transcriptome variation in segregating populations is still largely unknown. We sought to assess and map variation in transcript abundance in maize shoot apices in the intermated B73×Mo17 recombinant inbred line population. RNA–based sequencing (RNA–seq) allowed for the detection and quantification of the transcript abundance derived from 28,603 genes. For a majority of these genes, the population mean, coefficient of variation, and segregation patterns could be predicted by the parental expression levels. Expression quantitative trait loci (eQTL) mapping identified 30,774 eQTL including 96 trans-eQTL “hotspots,” each of which regulates the expression of a large number of genes. Interestingly, genes regulated by a trans-eQTL hotspot tend to be enriched for a specific function or act in the same genetic pathway. Also, genomic structural variation appeared to contribute to cis-regulation of gene expression. Besides genes showing Mendelian inheritance in the RIL population, we also found genes whose expression level and variation in the progeny could not be predicted based on parental difference, indicating that non-Mendelian factors also contribute to expression variation. Specifically, we found 145 genes that show patterns of expression reminiscent of paramutation such that all the progeny had expression levels similar to one of the two parents. Furthermore, we identified another 210 genes that exhibited unexpected patterns of transcript presence/absence. Many of these genes are likely to be gene fragments resulting from transposition, and the presence/absence of their transcripts could influence expression levels of their ancestral syntenic genes. Overall, our results contribute to the identification of novel expression patterns and broaden the understanding of transcriptional variation in plants.

Phenotypes are determined by the expression of genes, the environment, and the interaction of gene expression and the environment. However, a complete understanding of the inheritance of and genome-wide regulation of gene expression is lacking. One approach, called expression quantitative trait locus (eQTL) mapping provides the opportunity to examine the genome-wide inheritance and regulation of gene expression. In this paper, we conducted high-throughput sequencing of gene transcripts to examine gene expression in the shoot apex of a maize biparental mapping population. We quantified expression levels from 28,603 genes in the population and showed that the vast majority of genes exhibited the expected pattern of Mendelian inheritance. We genetically mapped the expression patterns and identified genomic regions associated with gene expression. Notably, we detected gene expression patterns that exhibited non-Mendelian inheritance. These included 145 genes that exhibited expression patterns in the progeny that were similar to only one of the parents and 210 genes with unexpected presence/absence expression patterns. The findings of non-Mendelian inheritance underscore the complexity of gene expression and provide a framework for understanding these complexities.

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.

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

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

Spreading of heterochromatin is limited to specific families of maize retrotransposons

Transposable elements (TEs) have the potential to act as controlling elements to influence the expression of genes and are often subject to heterochromatic silencing. The current paradigm suggests that heterochromatic silencing can spread beyond the borders of TEs and influence the chromatin state of neighboring low-copy sequences. This would allow TEs to condition obligatory or facilitated epialleles and act as controlling elements. The maize genome contains numerous families of class I TEs (retrotransposons) that are present in moderate to high copy numbers, and many are found in regions near genes, which provides an opportunity to test whether the spreading of heterochromatin from retrotransposons is prevalent. We have investigated the extent of heterochromatin spreading into DNA flanking each family of retrotransposons by profiling DNA methylation and di-methylation of lysine 9 of histone 3 (H3K9me2) in low-copy regions of the maize genome. The effects of different retrotransposon families on local chromatin are highly variable. Some retrotransposon families exhibit enrichment of heterochromatic marks within 800–1,200 base pairs of insertion sites, while other families exhibit very little evidence for the spreading of heterochromatic marks. The analysis of chromatin state in genotypes that lack specific insertions suggests that the heterochromatin in low-copy DNA flanking retrotransposons often results from the spreading of silencing marks rather than insertion-site preferences. Genes located near TEs that exhibit spreading of heterochromatin tend to be expressed at lower levels than other genes. Our findings suggest that a subset of retrotransposon families may act as controlling elements influencing neighboring sequences, while the majority of retrotransposons have little effect on flanking sequences.

Transposable elements comprise a substantial portion of many eukaryotic genomes. These mobile fragments of DNA can directly mutate genes through insertions into coding regions but may also affect the gene regulation through nearby insertions. There is evidence that the majority of transposable elements are epigenetically silenced, and in some cases this silencing may spread to neighboring sequences. This spreading of heterochromatin could create a significant fitness tradeoff between transposon silencing and gene expression. The maize genome has a complex organization with many genes flanked by retrotransposons, providing an opportunity to study the interaction of retrotransposons and genes. To survey the prevalence of heterochromatin spreading associated with different retrotransposon families, we profiled the spread of heterochromatin into nearby low copy sequences for 150 high copy retrotransposon families. While many retrotransposons exhibit little to no spreading of heterochromatin, there are some retrotransposon families that do exhibit spreading. Genes located near retrotransposons that spread heterochromatin have lower expression levels. The families of retrotransposons that spread heterochromatin marks to nearby low-copy sequences may have increased fitness costs for the host genome due to their suppression of genes located near insertions.

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

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

Paramutation: a trans-homolog interaction affecting heritable gene regulation

Paramutation describes both the process and results of trans-sensing between chromosomes that causes specific heritable changes in gene regulation. RNA molecules are implicated in mediating similar events in maize, mouse, and Drosophila. Changes in both small RNA profiles and cytosine methylation patterns in Arabidopsis hybrids represent a potential molecular equivalent to the interactions responsible for paramutations. Despite a seemingly unifying feature of RNA-directed changes, both recent and historical works show that paramutations in maize require plant-specific proteins and lack expected hallmarks of a trans-effect mediated solely by RNAs. Recent examples of nearby transposons affecting RNA polymerase II functions lead to an opinion that paramutations represent an emergent property of the transcriptional dynamics ongoing in plant genomes between repetitious features and nearby genes.

Heterosis: Many Genes, Many Mechanisms—End the Search for an Undiscovered Unifying Theory

Heterosis is the increase in vigor that is observed in progenies of matings of diverse individuals from different species, isolated populations, or selected strains within species or populations. Heterosis has been of immense economic value in agriculture and has important implications regarding the fitness and fecundity of individuals in natural populations. Genetic models based on complementation of deleterious alleles, especially in the context of linkage and epistasis, are consistent with many observed manifestations of heterosis. The search for the genes and alleles that underlie heterosis, as well as for broader allele-independent, genomewide mechanisms, has encompassed many species and systems. Common themes across these studies indicate that sequence diversity is necessary but not sufficient to produce heterotic phenotypes, and that the molecular pathways that produce heterosis involve chromatin modification, transcriptional control, translation and protein processing, and interactions between and within developmental and biochemical pathways. Taken together, there are many and diverse molecular mechanisms that translate DNA into phenotype, and it is the combination of all these mechanisms across many genes that produce heterosis in complex traits.

Genomic organization and sequence dynamics of the AvrPiz-t locus in Magnaporthe oryzae

Plants utilize multiple layers of defense mechanisms to fight against the invasion of diverse pathogens. The R gene mediates resistance, in most cases, dependent on the co-existence of its cognate pathogen-derived avirulence (Avr) gene. The rice blast R gene Piz-t corresponds in gene-for-gene fashion to the Magnaporthe oryzae Avr gene AvrPiz-t. In this study, we determined and compared the genomic sequences surrounding the AvrPiz-t gene in both avirulent and virulent isolates, designating as AvrPiz-t-ZB15 and avrPiz-t-70-15 regions, respectively. The sequence of the AvrPiz-t-ZB15 region is 120966 bp whereas avrPiz-t-70-15 is 146292 bp in length. The extreme sequence similarity and good synteny in gene order and content along with the absence of two predicted genes in the avrPiz-t-70-15 region were observed in the predicted protein-coding regions in the AvrPiz-t locus. Nevertheless, frequent presence/absence and highly dynamic organization of transposable elements (TEs) were identified, representing the major variation of the AvrPiz-t locus between different isolates. Moreover, TEs constitute 27.3% and 43.2% of the genomic contents of the AvrPiz-t-ZB15 and avrPiz-t-70-15 regions, respectively, indicating that TEs contribute largely to the organization and evolution of AvrPiz-t locus. The findings of this study suggest that M. oryzae could benefit in an evolutionary sense from the presence of active TEs in genes conferring avirulence and provide an ability to rapidly change and thus to overcome host R genes.

Dynamic evolution of bz orthologous regions in the Andropogoneae and other grasses

Genome structure exhibits remarkable plasticity within Zea mays. To examine how haplotype structure has evolved within the Andropogoneae tribe, we have analyzed the bz gene-rich region of maize (Zea mays), the Zea teosintes mays ssp. mexicana, luxurians and diploperennis, Tripsacum dactyloides, Coix lacryma-jobi and Sorghum propinquum. We sequenced and annotated BAC clones from these species and re-annotated the orthologous Sorghum bicolor region. Gene colinearity in the region is well conserved within the genus Zea. However, the orthologous regions of Coix and Sorghum exhibited several micro-rearrangements relative to Zea, including addition, truncation and deletion of genes. The stc1 gene, involved in the production of a terpenoid insect defense signal, is evolving particularly fast, and its progressive disappearance from some species is occurring by microhomology-mediated recombination. LTR retrotransposons are the main contributors to the dynamic evolution of the bz region. Common transposon insertion sites occur among haplotypes from different Zea mays sub-species, but not outside the species. As in Zea, different patterns of interspersion between genes and retrotransposons are observed in Sorghum. We estimate that the mean divergence times between maize and Tripsacum, Coix and Sorghum are 8.5, 12.1 and 12.4 million years ago, respectively, and that between Coix and Sorghum is 9.3 million years ago. A comparison of the bz orthologous regions of Zea, Sorghum and Coix with those of Brachypodium, Setaria and Oryza allows us to infer how the region has evolved by addition and deletion of genes in the approximately 50 million years since these genera diverged from a common progenitor.

Specific changes in total and mitochondrial proteomes are associated with higher levels of heterosis in maize hybrids

The phenomenon of hybrid vigor (heterosis) has long been harnessed by plant breeders to improve world food production. However, the changes that are essential for heterotic responses and the mechanisms responsible for heterosis remain undefined. Large increases in biomass and yield in high-heterosis hybrids suggest that alterations in bioenergetic processes may contribute to heterosis. Progeny from crosses between various inbred lines vary in the extent of vigor observed. Field-grown maize F₁ hybrids that consistently exhibited either low or high heterosis across a variety of environments were examined for changes in proteins that may be correlated with increased plant vigor and yield. Unpollinated ears at the time of flowering (ear shoots) were selected for the studies because they are metabolically active, rich in mitochondria, and the sizes of the ears are diagnostic of yield heterosis. Total protein and mitochondrial proteomes were compared among low- and higher-heterosis hybrids. Two-dimensional difference gel electrophoresis was used to identify allelic and/or isoform differences linked to heterosis. Identification of differentially regulated spots by mass spectrometry revealed proteins involved in stress responses as well as primary carbon and protein metabolism. Many of these proteins were identified in multiple spots, but analysis of their abundances by label-free mass spectrometry suggested that most of the expression differences were due to isoform variation rather than overall protein amount. Thus, our proteomics studies suggest that expression of specific alleles and/or post-translational modification of specific proteins correlate with higher levels of heterosis.

Transcriptional activation of Brassica napus β-ketoacyl-ACP synthase II with an engineered zinc finger protein transcription factor

Targeted gene regulation via designed transcription factors has great potential for precise phenotypic modification and acceleration of novel crop trait development. Canola seed oil composition is dictated largely by the expression of genes encoding enzymes in the fatty acid biosynthetic pathway. In the present study, zinc finger proteins (ZFPs) were designed to bind DNA sequences common to two canola β-ketoacyl-ACP Synthase II (KASII) genes downstream of their transcription start site. Transcriptional activators (ZFP-TFs) were constructed by fusing these ZFP DNA-binding domains to the VP16 transcriptional activation domain. Following transformation using Agrobacterium, transgenic events expressing ZFP-TFs were generated and shown to have elevated KASII transcript levels in the leaves of transgenic T(0) plants when compared to 'selectable marker only' controls as well as of T(1) progeny plants when compared to null segregants. In addition, leaves of ZFP-TF-expressing T(1) plants contained statistically significant decreases in palmitic acid (consistent with increased KASII activity) and increased total C18. Similarly, T(2) seed displayed statistically significant decreases in palmitic acid, increased total C18 and reduced total saturated fatty acid contents. These results demonstrate that designed ZFP-TFs can be used to regulate the expression of endogenous genes to elicit specific phenotypic modifications of agronomically relevant traits in a crop species.

Mutator System Derivatives Isolated from Sugarcane Genome Sequence

Mutator-like transposase is the most represented transposon transcript in the sugarcane transcriptome. Phylogenetic reconstructions derived from sequenced transcripts provided evidence that at least four distinct classes exist (I-IV) and that diversification among these classes occurred early in Angiosperms, prior to the divergence of Monocots/Eudicots. The four previously described classes served as probes to select and further sequence six BAC clones from a genomic library of cultivar R570. A total of 579,352 sugarcane base pairs were produced from these "Mutator system" BAC containing regions for further characterization. The analyzed genomic regions confirmed that the predicted structure and organization of the Mutator system in sugarcane is composed of two true transposon lineages, each containing a specific terminal inverted repeat and two transposase lineages considered to be domesticated. Each Mutator transposase class displayed a particular molecular structure supporting lineage specific evolution. MUSTANG, previously described domesticated genes, are located in syntenic regions across Sacharineae and, as expected for a host functional gene, posses the same gene structure as in other Poaceae. Two sequenced BACs correspond to hom(eo)logous locus with specific retrotransposon insertions that discriminate sugarcane haplotypes. The comparative studies presented, add information to the Mutator systems previously identified in the maize and rice genomes by describing lineage specific molecular structure and genomic distribution pattern in the sugarcane genome. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s12042-012-9104-y) contains supplementary material, which is available to authorized users.

Variation of presence/absence genes among Arabidopsis populations

BACKGROUND

Gene presence/absence (P/A) polymorphisms are commonly observed in plants and are important in individual adaptation and species differentiation. Detecting their abundance, distribution and variation among individuals would help to understand the role played by these polymorphisms in a given species. The recently sequenced 80 Arabidopsis genomes provide an opportunity to address these questions.

RESULTS

By systematically investigating these accessions, we identified 2,407 P/A genes (or 8.9%) absent in one or more genomes, averaging 444 absent genes per accession. 50.6% of P/A genes belonged to multi-copy gene families, or 31.0% to clustered genes. However, the highest proportion of P/A genes, outnumbered in singleton genes, was observed in the regions near centromeres. In addition, a significant correlation was observed between the P/A gene frequency among the 80 accessions and the diversity level at P/A loci. Furthermore, the proportion of P/A genes was different among functional gene categories. Finally, a P/A gene tree showed a diversified population structure in the worldwide Arabidopsis accessions.

CONCLUSIONS

An estimate of P/A genes and their frequency distribution in the worldwide Arabidopsis accessions was obtained. Our results suggest that there are diverse mechanisms to generate or maintain P/A genes, by which individuals and functionally different genes can selectively maintain P/A polymorphisms for a specific adaptation.

Maize HapMap2 identifies extant variation from a genome in flux

Whereas breeders have exploited diversity in maize for yield improvements, there has been limited progress in using beneficial alleles in undomesticated varieties. Characterizing standing variation in this complex genome has been challenging, with only a small fraction of it described to date. Using a population genetics scoring model, we identified 55 million SNPs in 103 lines across pre-domestication and domesticated Zea mays varieties, including a representative from the sister genus Tripsacum. We find that structural variations are pervasive in the Z. mays genome and are enriched at loci associated with important traits. By investigating the drivers of genome size variation, we find that the larger Tripsacum genome can be explained by transposable element abundance rather than an allopolyploid origin. In contrast, intraspecies genome size variation seems to be controlled by chromosomal knob content. There is tremendous overlap in key gene content in maize and Tripsacum, suggesting that adaptations from Tripsacum (for example, perennialism and frost and drought tolerance) can likely be integrated into maize.

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.

Maize (Zea mays L.) genome diversity as revealed by RNA-sequencing

Maize is rich in genetic and phenotypic diversity. Understanding the sequence, structural, and expression variation that contributes to phenotypic diversity would facilitate more efficient varietal improvement. RNA based sequencing (RNA-seq) is a powerful approach for transcriptional analysis, assessing sequence variation, and identifying novel transcript sequences, particularly in large, complex, repetitive genomes such as maize. In this study, we sequenced RNA from whole seedlings of 21 maize inbred lines representing diverse North American and exotic germplasm. Single nucleotide polymorphism (SNP) detection identified 351,710 polymorphic loci distributed throughout the genome covering 22,830 annotated genes. Tight clustering of two distinct heterotic groups and exotic lines was evident using these SNPs as genetic markers. Transcript abundance analysis revealed minimal variation in the total number of genes expressed across these 21 lines (57.1% to 66.0%). However, the transcribed gene set among the 21 lines varied, with 48.7% expressed in all of the lines, 27.9% expressed in one to 20 lines, and 23.4% expressed in none of the lines. De novo assembly of RNA-seq reads that did not map to the reference B73 genome sequence revealed 1,321 high confidence novel transcripts, of which, 564 loci were present in all 21 lines, including B73, and 757 loci were restricted to a subset of the lines. RT-PCR validation demonstrated 87.5% concordance with the computational prediction of these expressed novel transcripts. Intriguingly, 145 of the novel de novo assembled loci were present in lines from only one of the two heterotic groups consistent with the hypothesis that, in addition to sequence polymorphisms and transcript abundance, transcript presence/absence variation is present and, thereby, may be a mechanism contributing to the genetic basis of heterosis.

The turbulent life of Sirevirus retrotransposons and the evolution of the maize genome: more than ten thousand elements tell the story

Sireviruses are one of the three genera of Copia long terminal repeat (LTR) retrotransposons, exclusive to and highly abundant in plants, and with a unique, among retrotransposons, genome structure. Yet, perhaps due to the few references to the Sirevirus origin of some families, compounded by the difficulty in correctly assigning retrotransposon families into genera, Sireviruses have hardly featured in recent research. As a result, analysis at this key level of classification and details of their colonization and impact on plant genomes are currently lacking. Recently, however, it became possible to accurately assign elements from diverse families to this genus in one step, based on highly conserved sequence motifs. Hence, Sirevirus dynamics in the relatively obese maize genome can now be comprehensively studied. Overall, we identified >10 600 intact and approximately 28 000 degenerate Sirevirus elements from a plethora of families, some brought into the genus for the first time. Sireviruses make up approximately 90% of the Copia population and it is the only genus that has successfully infiltrated the genome, possibly by experiencing intense amplification during the last 600 000 years, while being constantly recycled by host mechanisms. They accumulate in chromosome-distal gene-rich areas, where they insert in between gene islands, mainly in preferred zones within their own genomes. Sirevirus LTRs are heavily methylated, while there is evidence for a palindromic consensus target sequence. This work brings Sireviruses in the spotlight, elucidating their lifestyle and history, and suggesting their crucial role in the current genomic make-up of maize, and possibly other plant hosts.

Genome evolution in diploid and tetraploid Coffea species as revealed by comparative analysis of orthologous genome segments

Sequence comparison of orthologous regions enables estimation of the divergence between genomes, analysis of their evolution and detection of particular features of the genomes, such as sequence rearrangements and transposable elements. Despite the economic importance of Coffea species, little genomic information is currently available. Coffea is a relatively young genus that includes more than one hundred diploid species and a single tetraploid species. Three Coffea orthologous regions of 470-900 kb were analyzed and compared: both subgenomes of allotetraploid Coffea arabica (contributed by the diploid species Coffea eugenioides and Coffea canephora) and the genome of diploid C. canephora. Sequence divergence was calculated on global alignments or on coding and non-coding sequences separately. A search for transposable elements detected 43 retrotransposons and 198 transposons in the sequences analyzed. Comparative insertion analysis made it possible to locate 165 TE insertions in the phylogenetic tree of the three genomes/subgenomes. In the tetraploid C. arabica, a homoeologous non-reciprocal transposition (HNRT) was detected and characterized: a 50 kb region of the C. eugenioides derived subgenome replaced the C. canephora derived counterpart. Comparative sequence analysis on three Coffea genomes/subgenomes revealed almost perfect gene synteny, low sequence divergence and a high number of shared transposable elements. Compared to the results of similar analysis in other genera (Aegilops/Triticum and Oryza), Coffea genomes/subgenomes appeared to be dramatically less diverged, which is consistent with the relatively recent radiation of the Coffea genus. Based on nucleotide substitution frequency, the HNRT was dated at 10,000-50,000 years BP, which is also the most recent estimation of the origin of C. arabica.