The genome of maize (Zea mays ssp. mays) consists mostly of transposable elements (TEs) and varies in size among lines. This variation extends to other species in the genus Zea: although maize and Zea luxurians diverged only ∼140,000 years ago, their genomes differ in size by ∼50%. We used paired-end Illumina sequencing to evaluate the potential contribution of TEs to the genome size difference between these two species. We aligned the reads both to a filtered gene set and to an exemplar database of unique repeats representing 1,514 TE families; ∼85% of reads mapped against TE repeats in both species. The relative contribution of TE families to the B73 genome was highly correlated with previous estimates, suggesting that reliable estimates of TE content can be obtained from short high-throughput sequencing reads, even at low coverage. Because we used paired-end reads, we could assess whether a TE was near a gene by determining if one paired read mapped to a TE and the second read mapped to a gene. Using this method, Class 2 DNA elements were found significantly more often in genic regions than Class 1 RNA elements, but Class 1 elements were found more often near other TEs. Overall, we found that both Class 1 and 2 TE families account for ∼70% of the genome size difference between B73 and luxurians. Interestingly, the relative abundance of TE families was conserved between species (r = 0.97), suggesting genome-wide control of TE content rather than family-specific effects.

Background

The species Zea mays includes both domesticated maize (ssp. mays) and its closest wild relatives known as the teosintes. While genetic and archaeological studies have provided a well-established history of Z. mays evolution, there is currently minimal description of its current and past distribution. Here, we implemented species distribution modeling using paleoclimatic models of the last interglacial (LI; ∼135,000 BP) and the last glacial maximum (LGM; ∼21,000 BP) to hindcast the distribution of Zea mays subspecies over time and to revisit current knowledge of its phylogeography and evolutionary history.

Methodology/principal findings

Using a large occurrence data set and the distribution modeling MaxEnt algorithm, we obtained robust present and past species distributions of the two widely distributed teosinte subspecies (ssps. parviglumis and mexicana) revealing almost perfect complementarity, stable through time, of their occupied distributions. We also investigated the present distributions of primitive maize landraces, which overlapped but were broader than those of the teosintes. Our data reinforced the idea that little historical gene flow has occurred between teosinte subspecies, but maize has served as a genetic bridge between them. We observed an expansion of teosinte habitat from the LI, consistent with population genetic data. Finally, we identified locations potentially serving as refugia for the teosintes throughout epochs of climate change and sites that should be targeted in future collections.

Conclusion/significance

The restricted and highly contrasting ecological niches of the wild teosintes differ substantially from domesticated maize. Variables determining the distributions of these taxa can inform future considerations of local adaptation and the impacts of climate change. Our assessment of the changing distributions of Zea mays taxa over time offers a unique glimpse into the history of maize, highlighting a strategy for the study of domestication that may prove useful for other species.