Genetic signals of origin, spread, and introgression in a large sample of maize landraces

Archaeological findings show the extent of primitive characteristics of maize in South America

Scientific research has suggested that maize spread from Mexico and arrived in lowland South America in a state of partial domestication. However, archaeological samples with primitive morphological characteristics that corroborate this finding have not been recorded in the region thus far. Unexpectedly, many samples were identified in the Peruaçu Valley with characteristics never previously observed in South America. These archaeological samples with primitive characteristics, which are the focus of this work, represent the furthest records from the center of origin of the species and the longest duration of the maintenance of such characteristics (between 1010 and 570 years before present). The findings of this study, including archaeological samples, native races, and samples of teosinte, attest to a long history of maize diversification in lowland South America.

Genomic basis determining root system architecture in maize

KEY MESSAGE

A total of 389 and 344 QTLs were identified by GWAS and QTL mapping explaining accumulatively 32.2-65.0% and 23.7-63.4% of phenotypic variation for 14 shoot-borne root traits using more than 1300 individuals across multiple field trails. Efficient nutrient and water acquisition from soils depends on the root system architecture (RSA). However, the genetic determinants underlying RSA in maize remain largely unexplored. In this study, we conducted a comprehensive genetic analysis for 14 shoot-borne root traits using 513 inbred lines and 800 individuals from four recombinant inbred line (RIL) populations at the mature stage across multiple field trails. Our analysis revealed substantial phenotypic variation for these 14 root traits, with a total of 389 and 344 QTLs identified through genome-wide association analysis (GWAS) and linkage analysis, respectively. These QTLs collectively explained 32.2-65.0% and 23.7-63.4% of the trait variation within each population. Several a priori candidate genes involved in auxin and cytokinin signaling pathways, such as IAA26, ARF2, LBD37 and CKX3, were found to co-localize with these loci. In addition, a total of 69 transcription factors (TFs) from 27 TF families (MYB, NAC, bZIP, bHLH and WRKY) were found for shoot-borne root traits. A total of 19 genes including PIN3, LBD15, IAA32, IAA38 and ARR12 and 19 GWAS signals were overlapped with selective sweeps. Further, significant additive effects were found for root traits, and pyramiding the favorable alleles could enhance maize root development. These findings could contribute to understand the genetic basis of root development and evolution, and provided an important genetic resource for the genetic improvement of root traits in maize.

Cytoplasmic genome contributions to domestication and improvement of modern maize

BACKGROUND

Studies on maize evolution and domestication are largely limited to the nuclear genomes, and the contribution of cytoplasmic genomes to selection and domestication of modern maize remains elusive. Maize cytoplasmic genomes have been classified into fertile (NA and NB) and cytoplasmic-nuclear male-sterility (CMS-S, CMS-C, and CMS-T) groups, but their contributions to modern maize breeding have not been systematically investigated.

RESULTS

Here we report co-selection and convergent evolution between nuclear and cytoplasmic genomes by analyzing whole genome sequencing data of 630 maize accessions modern maize and its relatives, including 24 fully assembled mitochondrial and chloroplast genomes. We show that the NB cytotype is associated with the expansion of modern maize to North America, gradually replaces the fertile NA cytotype probably through unequal division, and predominates in over 90% of modern elite inbred lines. The mode of cytoplasmic evolution is increased nucleotypic diversity among the genes involved in photosynthesis and energy metabolism, which are driven by selection and domestication. Furthermore, genome-wide association study reveals correlation of cytoplasmic nucleotypic variation with key agronomic and reproductive traits accompanied with the diversification of the nuclear genomes.

CONCLUSIONS

Our results indicate convergent evolution between cytoplasmic and nuclear genomes during maize domestication and breeding. These new insights into the important roles of mitochondrial and chloroplast genomes in maize domestication and improvement should help select elite inbred lines to improve yield stability and crop resilience of maize hybrids.

Maize pollen carry bacteria that suppress a fungal pathogen that enters through the male gamete fertilization route

In flowering plants, after being released from pollen grains, the male gametes use the style channel to migrate towards the ovary where they fertilize awaiting eggs. Environmental pathogens exploit the style passage, resulting in diseased progeny seed. The belief is that pollen also transmits pathogens into the style. By contrast, we hypothesized that pollen carries beneficial microbes that suppress environmental pathogens on the style passage. No prior studies have reported pollen-associated bacterial functions in any plant species. Here, bacteria were cultured from maize (corn) pollen encompassing wild ancestors and farmer-selected landraces from across the Americas, grown in a common field in Canada for one season. In total, 298 bacterial isolates were cultured, spanning 45 genera, 103 species, and 88 OTUs, dominated by Pantoea, Bacillus, Pseudomonas, Erwinia, and Microbacterium. Full-length 16S DNA-based taxonomic profiling showed that 78% of bacterial taxa from the major wild ancestor of maize (Parviglumis teosinte) were present in at least one cultivated landrace. The species names of the bacterial isolates were used to search the pathogen literature systematically; this preliminary evidence predicted that the vast majority of the pollen-associated bacteria analyzed are not maize pathogens. The pollen-associated bacteria were tested in vitro against a style-invading Fusarium pathogen shown to cause Gibberella ear rot (GER): 14 isolates inhibited this pathogen. Genome mining showed that all the anti-Fusarium bacterial species encode phzF, associated with biosynthesis of the natural fungicide, phenazine. To mimic the male gamete migration route, three pollen-associated bacterial strains were sprayed onto styles (silks), followed by Fusarium inoculation; these bacteria reduced GER symptoms and mycotoxin accumulation in progeny seed. Confocal microscopy was used to search for direct evidence that pollen-associated bacteria can defend living silks against Fusarium graminearum (Fg); bacterial strain AS541 (Kluyvera intermedia), isolated from pollen of ancestral Parviglumis, was observed to colonize the susceptible style/silk entry points of Fg (silk epidermis, trichomes, wounds). Furthermore, on style/silk tissue, AS541 colonized/aggregated on Fg hyphae, and was associated with Fg hyphal breaks. These results suggest that pollen has the potential to carry bacteria that can defend the style/silk passage against an environmental pathogen - a novel observation.

Conservation and diversity of the pollen microbiome of Pan-American maize using PacBio and MiSeq

Pollen is a vector for diversification, fitness-selection, and transmission of plant genetic material. The extent to which the pollen microbiome may contribute to host diversification is largely unknown, because pollen microbiome diversity within a plant species has not been reported, and studies have been limited to conventional short-read 16S rRNA gene sequencing (e.g., V4-MiSeq) which suffers from poor taxonomic resolution. Here we report the pollen microbiomes of 16 primitive and traditional accessions of maize (corn) selected by indigenous peoples across the Americas, along with the modern U.S. inbred B73. The maize pollen microbiome has not previously been reported. The pollen microbiomes were identified using full-length (FL) 16S rRNA gene PacBio SMRT sequencing compared to V4-MiSeq. The Pan-American maize pollen microbiome encompasses 765 taxa spanning 39 genera and 46 species, including known plant growth promoters, insect-obligates, plant pathogens, nitrogen-fixers and biocontrol agents. Eleven genera and 13 species composed the core microbiome. Of 765 taxa, 63% belonged to only four genera: 28% were Pantoea, 15% were Lactococcus, 11% were Pseudomonas, and 10% were Erwinia. Interestingly, of the 215 Pantoea taxa, 180 belonged to a single species, P. ananatis. Surprisingly, the diversity within P. ananatis ranged nearly 10-fold amongst the maize accessions analyzed (those with ≥3 replicates), despite being grown in a common field. The highest diversity within P. ananatis occurred in accessions that originated near the center of diversity of domesticated maize, with reduced diversity associated with the north-south migration of maize. This sub-species diversity was revealed by FL-PacBio but missed by V4-MiSeq. V4-MiSeq also mis-identified some dominant genera captured by FL-PacBio. The study, though limited to a single season and common field, provides initial evidence that pollen microbiomes reflect evolutionary and migratory relationships of their host plants.

Two teosintes made modern maize

The origins of maize were the topic of vigorous debate for nearly a century, but neither the current genetic model nor earlier archaeological models account for the totality of available data, and recent work has highlighted the potential contribution of a wild relative, Zea mays ssp. mexicana. Our population genetic analysis reveals that the origin of modern maize can be traced to an admixture between ancient maize and Zea mays ssp. mexicana in the highlands of Mexico some 4000 years after domestication began. We show that variation in admixture is a key component of maize diversity, both at individual loci and for additive genetic variation underlying agronomic traits. Our results clarify the origin of modern maize and raise new questions about the anthropogenic mechanisms underlying dispersal throughout the Americas.

Teosinte Pollen Drive guides maize diversification and domestication by RNAi

Meiotic drivers subvert Mendelian expectations by manipulating reproductive development to bias their own transmission. Chromosomal drive typically functions in asymmetric female meiosis, while gene drive is normally postmeiotic and typically found in males. Using single molecule and single-pollen genome sequencing, we describe Teosinte Pollen Drive, an instance of gene drive in hybrids between maize (Zea mays ssp. mays) and teosinte mexicana (Zea mays ssp. mexicana), that depends on RNA interference (RNAi). 22nt small RNAs from a non-coding RNA hairpin in mexicana depend on Dicer-Like 2 (Dcl2) and target Teosinte Drive Responder 1 (Tdr1), which encodes a lipase required for pollen viability. Dcl2, Tdr1, and the hairpin are in tight pseudolinkage on chromosome 5, but only when transmitted through the male. Introgression of mexicana into early cultivated maize is thought to have been critical to its geographical dispersal throughout the Americas, and a tightly linked inversion in mexicana spans a major domestication sweep in modern maize. A survey of maize landraces and sympatric populations of teosinte mexicana reveals correlated patterns of admixture among unlinked genes required for RNAi on at least 4 chromosomes that are also subject to gene drive in pollen from synthetic hybrids. Teosinte Pollen Drive likely played a major role in maize domestication and diversification, and offers an explanation for the widespread abundance of "self" small RNAs in the germlines of plants and animals.

Genomic diversity and population structure of teosinte (Zea spp.) and its conservation implications

The wild species of the genus Zea commonly named teosintes, comprise nine different taxa, distributed from northern Mexico to Costa Rica. Although this genus of plants has been extensively studied from a morphological, ecogeographical and genetic point of view, most contributions have been limited to the study of a few populations and taxa. To understand the great variability that exists between and within teosinte species, it is necessary to include the vast majority of known populations. In this context, the objective of this work was to evaluate the diversity and genomic structure of 276 teosinte populations. Molecular analyzes were performed with 3,604 plants and with data from 33,929 SNPs. The levels of genetic diversity by taxonomic group show a marked difference between species, races and sections, where the highest values of genomic diversity was found in ssp. parviglumis and ssp. mexicana. The lower values were obtained for the Luxuriantes section as well as ssp. huehuetenagensis of the section Zea. The results of structure show that there is a great genetic differentiation in all the taxonomic groups considered. For ssp. parviglumis and mexicana, which are the taxa with the largest number of populations, a marked genomic differentiation was found that is consistent with their geographic distribution patterns. These results showed a loss of diversity in several teosinte populations, making a strong case for further collection, and ex situ and in situ conservation. Also, this study highlights the importance of integrating genomic diversity and structure for the applications of conservation and management.

Genetic Incompatibilities and Evolutionary Rescue by Wild Relatives Shaped Grain Amaranth Domestication

Crop domestication and the subsequent expansion of crops have long been thought of as a linear process from a wild ancestor to a domesticate. However, evidence of gene flow from locally adapted wild relatives that provided adaptive alleles into crops has been identified in multiple species. Yet, little is known about the evolutionary consequences of gene flow during domestication and the interaction of gene flow and genetic load in crop populations. We study the pseudo-cereal grain amaranth that has been domesticated three times in different geographic regions of the Americas. We quantify the amount and distribution of gene flow and genetic load along the genome of the three grain amaranth species and their two wild relatives. Our results show ample gene flow between crop species and between crops and their wild relatives. Gene flow from wild relatives decreased genetic load in the three crop species. This suggests that wild relatives could provide evolutionary rescue by replacing deleterious alleles in crops. We assess experimental hybrids between the three crop species and found genetic incompatibilities between one Central American grain amaranth and the other two crop species. These incompatibilities might have created recent reproductive barriers and maintained species integrity today. Together, our results show that gene flow played an important role in the domestication and expansion of grain amaranth, despite genetic species barriers. The domestication of plants was likely not linear and created a genomic mosaic by multiple contributors with varying fitness effects for today's crops.

Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Singh A, Kumar A. Unlocking the potential of ancient hexaploid Indian dwarf wheat, Tritium sphaerococcum for grain quality improvement

Wild and ancient wheat are considered to be a rich source of nutrients and better stress tolerant, hence being re-considered for mainstreaming its cultivation by the farmers and bringing it back to the food basket. In the present study, thirty-four diverse accessions of Indian dwarf wheat, Triticum sphaerococcum conserved in the Indian National Genebank were evaluated for thirteen-grain quality parameters namely thousand-grain weight (TGW), hectolitre weight (HW), sedimentation value (Sed), grain hardness index (HI), protein (Pro), albumin (Alb), globulin (Glo), gliadin (Gli), glutenin (Glu), gluten, lysine (Lys), Fe²⁺ and Zn²⁺ content, and four antioxidant enzymes activities. Substantial variations were recorded for studied traits. TGW, HW, Sed, HI, Pro, Alb, Glo, Gli, Glu, Gluten, Lys, Fe²⁺, and Zn²⁺ varied from 26.50-45.55 g, 70.50-86.00 kg/hl, 24.00-38.00 ml, 40.49-104.90, 15.34-19.35%, 17.60-40.31 mg/g, 10.75-16.56 mg/g, 26.35-44.94 mg/g, 24.47-39.56 mg/g, 55.33-75.06 mg/g, 0.04-0.29%, 42.72-90.72 ppm, and 11.45-25.70 ppm, respectively. Among antioxidants, peroxidase (POX), catalase (CAT), glutathione reductase (GR), and superoxide dismutase (SOD) activity ranged from 0.06-0.60 unit/ml, 0.02-0.61 unit/ml, 0.11-2.26 unit/ml, and 0.14-0.97 unit/ml, respectively. Hardness Index was positively associated with Pro and Zn²⁺ content whereas Lys was negatively associated with gluten content. Likewise, gluten and Fe²⁺ content had a positive association with the major protein fraction i.e., Gli and Glu. Hierarchical cluster analysis grouped 34 accessions into four clusters and the major group had nine indigenous and eight exotic accessions. We also validated high GPC accessions and EC182958 (17.16%), EC187176 and EC182945 (16.16%), EC613057 (15.79%), IC634028 (15.72%) and IC533826 (15.01%) were confirmed with more than 15% GPC. Also, superior trait-specific accessions namely, EC187167, IC534021, EC613055, EC180066, and EC182959 for low gluten content and IC384530, EC313761, EC180063, IC397363, EC10494 for high iron content (>76.51) were identified that may be used in wheat quality improvement for nutritional security of mankind.

Genome-Wide Association Study Discovers Novel Germplasm Resources and Genetic Loci with Resistance to Gibberella Ear Rot Caused by Fusarium graminearum

Gibberella ear rot (GER) in maize caused by Fusarium graminearum is one of the most devastating maize diseases reducing grain yield and quality worldwide. The genetic bases of maize GER resistance remain largely unknown. Using artificial inoculation across multiple environments, the GER severity of an association panel consisting of 316 diverse inbred lines was observed with wide phenotypic variation. In the association panel, a genome-wide association study using a general linear model identified 69 single-nucleotide polymorphisms (SNPs) significantly associated with GER resistance at the threshold of 2.04 × 10-5, and the average phenotypic variation explained (PVE) of these SNPs was 5.09%. We also conducted a genome-wide association study analysis using a mixed linear model at a threshold of 1.0 × 10-4, and 16 significantly associated SNPs with an average PVE of 4.73% were detected. A combined general linear model and mixed linear model method obtained 10 co-localized significantly associated SNPs linked to GER resistance, including the most significant SNP (PZE-105079915) with the greatest PVE value, 9.07%, at bin 5.05 following 10 candidate genes. These findings are significant for the exploration of the complicated genetic variations in maize GER resistance. The regions and genes identified herein provide a list of candidate targets for further investigation, in addition to the elite germplasm resources that can be used for breeding GER resistance in maize.

Genotyping of DNA pools identifies untapped landraces and genomic regions to develop next-generation varieties

Landraces, that is, traditional varieties, have a large diversity that is underexploited in modern breeding. A novel DNA pooling strategy was implemented to identify promising landraces and genomic regions to enlarge the genetic diversity of modern varieties. As proof of concept, DNA pools from 156 American and European maize landraces representing 2340 individuals were genotyped with an SNP array to assess their genome-wide diversity. They were compared to elite cultivars produced across the 20th century, represented by 327 inbred lines. Detection of selective footprints between landraces of different geographic origin identified genes involved in environmental adaptation (flowering times, growth) and tolerance to abiotic and biotic stress (drought, cold, salinity). Promising landraces were identified by developing two novel indicators that estimate their contribution to the genome of inbred lines: (i) a modified Roger's distance standardized by gene diversity and (ii) the assignation of lines to landraces using supervised analysis. It showed that most landraces do not have closely related lines and that only 10 landraces, including famous landraces as Reid's Yellow Dent, Lancaster Surecrop and Lacaune, cumulated half of the total contribution to inbred lines. Comparison of ancestral lines directly derived from landraces with lines from more advanced breeding cycles showed a decrease in the number of landraces with a large contribution. New inbred lines derived from landraces with limited contributions enriched more the haplotype diversity of reference inbred lines than those with a high contribution. Our approach opens an avenue for the identification of promising landraces for pre-breeding.

Haplotype-phased and chromosome-level genome assembly of Puccinia polysora, a giga-scale fungal pathogen causing southern corn rust

Rust fungi are characterized by large genomes with high repeat content and have two haploid nuclei in most life stages, which makes achieving high-quality genome assemblies challenging. Here, we described a pipeline using HiFi reads and Hi-C data to assemble a gigabase-sized fungal pathogen, Puccinia polysora f.sp. zeae, to haplotype-phased and chromosome-scale. The final assembled genome is 1.71 Gbp, with ~850 Mbp and 18 chromosomes in each haplotype, being currently one of the two giga-scale fungi assembled to chromosome level. Transcript-based annotation identified 47,512 genes for the dikaryotic genome with a similar number for each haplotype. A high level of interhaplotype variation was found with 10% haplotype-specific BUSCO genes, 5.8 SNPs/kbp, and structural variation accounting for 3% of the genome size. The P. polysora genome displayed over 85% repeat contents, with genome-size expansion and copy number increasing of species-specific orthogroups. Interestingly, these features did not affect overall synteny with other Puccinia species having smaller genomes. Fine-time-point transcriptomics revealed seven clusters of coexpressed secreted proteins that are conserved between two haplotypes. The fact that candidate effectors interspersed with all genes indicated the absence of a "two-speed genome" evolution in P. polysora. Genome resequencing of 79 additional isolates revealed a clonal population structure of P. polysora in China with low geographic differentiation. Nevertheless, a minor population differentiated from the major population by having mutations on secreted proteins including AvrRppC, indicating the ongoing virulence to evade recognition by RppC, a major resistance gene in Chinese corn cultivars. The high-quality assembly provides valuable genomic resources for future studies on disease management and the evolution of P. polysora.

Large haplotypes highlight a complex age structure within the maize pan-genome

The genomes of maize and other eukaryotes contain stable haplotypes in regions of low recombination. These regions, including centromeres, long heterochromatic blocks, and rDNA arrays, have been difficult to analyze with respect to their diversity and origin. Greatly improved genome assemblies are now available that enable comparative genomics over these and other nongenic spaces. Using 26 complete maize genomes, we developed methods to align intergenic sequences while excluding genes and regulatory regions. The centromere haplotypes (cenhaps) extend for megabases on either side of the functional centromere regions and appear as evolutionary strata, with haplotype divergence/coalescence times dating as far back as 450 thousand years ago (kya). Application of the same methods to other low recombination regions (heterochromatic knobs and rDNA) and all intergenic spaces revealed that deep coalescence times are ubiquitous across the maize pan-genome. Divergence estimates vary over a broad timescale with peaks at ∼16 and 300 kya, reflecting a complex history of gene flow among diverging populations and changes in population size associated with domestication. Cenhaps and other long haplotypes provide vivid displays of this ancient diversity.

A common resequencing-based genetic marker data set for global maize diversity

Maize (Zea mays ssp. mays) populations exhibit vast ranges of genetic and phenotypic diversity. As sequencing costs have declined, an increasing number of projects have sought to measure genetic differences between and within maize populations using whole-genome resequencing strategies, identifying millions of segregating single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels). Unlike older genotyping strategies like microarrays and genotyping by sequencing, resequencing should, in principle, frequently identify and score common genetic variants. However, in practice, different projects frequently employ different analytical pipelines, often employ different reference genome assemblies and consistently filter for minor allele frequency within the study population. This constrains the potential to reuse and remix data on genetic diversity generated from different projects to address new biological questions in new ways. Here, we employ resequencing data from 1276 previously published maize samples and 239 newly resequenced maize samples to generate a single unified marker set of approximately 366 million segregating variants and approximately 46 million high-confidence variants scored across crop wild relatives, landraces as well as tropical and temperate lines from different breeding eras. We demonstrate that the new variant set provides increased power to identify known causal flowering-time genes using previously published trait data sets, as well as the potential to track changes in the frequency of functionally distinct alleles across the global distribution of modern maize.

The role of transposon inverted repeats in balancing drought tolerance and yield-related traits in maize

The genomic basis underlying the selection for environmental adaptation and yield-related traits in maize remains poorly understood. Here we carried out genome-wide profiling of the small RNA (sRNA) transcriptome (sRNAome) and transcriptome landscapes of a global maize diversity panel under dry and wet conditions and uncover dozens of environment-specific regulatory hotspots. Transgenic and molecular studies of Drought-Related Environment-specific Super eQTL Hotspot on chromosome 8 (DRESH8) and ZmMYBR38, a target of DRESH8-derived small interfering RNAs, revealed a transposable element-mediated inverted repeats (TE-IR)-derived sRNA- and gene-regulatory network that balances plant drought tolerance with yield-related traits. A genome-wide scan revealed that TE-IRs associate with drought response and yield-related traits that were positively selected and expanded during maize domestication. These results indicate that TE-IR-mediated posttranscriptional regulation is a key molecular mechanism underlying the tradeoff between crop environmental adaptation and yield-related traits, providing potential genomic targets for the breeding of crops with greater stress tolerance but uncompromised yield.

Morphological and genetic diversity of maize landraces along an altitudinal gradient in the Southern Andes

Maize (Zea mays ssp. mays) is a major cereal crop worldwide and is traditionally or commercially cultivated almost all over the Americas. The North-Western Argentina (NWA) region constitutes one of the main diversity hotspots of the Southern Andes, with contrasting landscapes and a large number of landraces. Despite the extensive collections performed by the "Banco Activo de Germoplasma INTA Pergamino, Argentina" (BAP), most of them have not been characterized yet. Here we report the morphological and molecular evaluation of 30 accessions collected from NWA, along an altitudinal gradient between 1120 and 2950 meters above sea level (masl). Assessment of morphological variation in a common garden allowed the discrimination of two groups, which differed mainly in endosperm type and overall plant size. Although the groups retrieved by the molecular analyses were not consistent with morphological clusters, they showed a clear pattern of altitudinal structuring. Affinities among accessions were not in accordance with racial assignments. Overall, our results revealed that there are two maize gene pools co-existing in NWA, probably resulting from various waves of maize introduction in pre-Columbian times as well as from the adoption of modern varieties by local farmers. In conclusion, the NWA maize landraces preserved at the BAP possess high morphological and molecular variability. Our results highlight their potential as a source of diversity for increasing the genetic basis of breeding programs and provide useful information to guide future sampling and conservation efforts.

ABA-inducible DEEPER ROOTING 1 improves adaptation of maize to water deficiency

Root architecture remodelling is critical for forage moisture in water-limited soil. DEEPER ROOTING 1 (DRO1) in Oryza, Arabidopsis, and Prunus has been reported to improve drought avoidance by promoting roots to grow downward and acquire water from deeper soil. In the present study, we found that ZmDRO1 responded more strongly to abscisic acid (ABA)/drought induction in Zea mays ssp. mexicana, an ancestral species of cultivated maize, than in B73. It was proposed that this is one of the reasons why Zea mays ssp. mexicana has a more noticeable change in the downward direction angle of the root and fewer biomass penalties under water-deficient conditions. Thus, a robust, synthetic ABA/drought-inducible promoter was used to control the expression of ZmDRO1^B73 in Arabidopsis and cultivated maize for drought-resistant breeding. Interestingly, ABA-inducible ZmDRO1 promoted a larger downward root angle and improved grain yield by more than 40% under water-limited conditions. Collectively, these results revealed that different responses to ABA/drought induction of ZmDRO1 confer different drought avoidance abilities, and we demonstrated the application of ZmDRO1 via an ABA-inducible strategy to alter the root architecture of modern maize to improve drought adaptation in the field.

Maize stomatal responses against the climate change

Drought and heat, in the context of climate change, are expected to increase in many agricultural areas across the globe. Among current abiotic stresses, they are the most limiting factors that influence crop growth and productivity. Maize is one of most widely produced crops of the world, being the first in grain production with a yield that exceeded 1.1 billion tons in 2021. Despite its wide distribution in semi-arid regions, it is highly vulnerable to climate change, which triggers important losses in its productivity. This article explores how maize yield may persevere through climate change by focusing on the stomatal regulation of gas exchange. The emerging picture unravels that maize copes with drought stress by reducing stomatal size and stomatal pore area, and increasing stomatal density, which, in turn, reduces transpiration and photosynthetic rate. When drought and heat co-occur, heat enhances stomatal response to drought stress. To avoid plant heat damage, the decline in stomatal aperture could trigger the expansion of the distance of action, from the longitudinal leaf veins, of ZmSHR1, which might act to positively regulate ZmSPCHs/ZmICE1 heterodimers, increasing the stomatal density. Only when drought is not very severe, elevated CO2 levels reduce yield losses. The knowledge of the upcoming climate changes together with the prediction of the developmental and physiological stomatal responses will allow not only to anticipate maize yield in the next years, but also to contribute to the correct decision-making in the management of this important crop.

Quantitative evaluation of nonlinear methods for population structure visualization and inference

Population structure (also called genetic structure and population stratification) is the presence of a systematic difference in allele frequencies between subpopulations in a population as a result of nonrandom mating between individuals. It can be informative of genetic ancestry, and in the context of medical genetics, it is an important confounding variable in genome-wide association studies. Recently, many nonlinear dimensionality reduction techniques have been proposed for the population structure visualization task. However, an objective comparison of these techniques has so far been missing from the literature. In this article, we discuss the previously proposed nonlinear techniques and some of their potential weaknesses. We then propose a novel quantitative evaluation methodology for comparing these nonlinear techniques, based on populations for which pedigree is known a priori either through artificial selection or simulation. Based on this evaluation metric, we find graph-based algorithms such as t-SNE and UMAP to be superior to principal component analysis, while neural network-based methods fall behind.

The genomic signature of wild-to-crop introgression during the domestication of scarlet runner bean (Phaseolus coccineus L.)

The scarlet runner bean (Phaseolus coccineus) is one of the five domesticated Phaseolus species. It is cultivated in small-scale agriculture in the highlands of Mesoamerica for its dry seeds and immature pods, and unlike the other domesticated beans, P. coccineus is an open-pollinated legume. Contrasting with its close relative, the common bean, few studies focusing on its domestication history have been conducted. Demographic bottlenecks associated with domestication might reduce genetic diversity and facilitate the accumulation of deleterious mutations. Conversely, introgression from wild relatives could be a source of variation. Using Genotyping by Sequencing data (79,286 single-nucleotide variants) from 237 cultivated and wild samples, we evaluated the demographic history of traditional varieties from different regions of Mexico and looked for evidence of introgression between sympatric wild and cultivated populations. Traditional varieties have high levels of diversity, even though there is evidence of a severe initial genetic bottleneck followed by a population expansion. Introgression from wild to domesticated populations was detected, which might contribute to the recovery of the genetic variation. Introgression has occurred at different times: constantly in the center of Mexico; recently in the North West; and anciently in the South. Several factors are acting together to increase and maintain genetic diversity in P. coccineus cultivars, such as demographic expansion and introgression. Wild relatives represent a valuable genetic resource and have played a key role in scarlet runner bean evolution via introgression into traditional varieties.

Identification of Fusarium verticillioides Resistance Alleles in Three Maize Populations With Teosinte Gene Introgression

Fusarium ear rot (FER) is a common fungal disease in maize (Zea mays L.) caused by Fusarium verticillioides. Resistant germplasm resources for FER are rare in cultivated maize; however, teosintes (Z. mays ssp. parviglumis and Z. mays ssp. diploperennis), which are wild-type species of maize, have the potential to offer a novel source of resistance alleles to enhance pathogen resistance in modern maize. Therefore, the aim of this study was to identify favorable alleles that confer significant levels of resistance toward FER. Three populations of BC₂F₈ recombinant inbred lines (RILs) were developed by crossing two different teosintes, Z. diploperennis and Z. parviglumis, with maize inbred lines B73 and Zheng58, and were screened for FER resistance. We found that Z. diploperennis and Z. parviglumis had higher resistance toward F. verticillioides in the leaves than B73 and Zheng58. However, the resistance toward F. verticillioides in the leaf and ear was unrelated among RILs. FER resistance was positively correlated with grain yield in the B73 × diploperennis (BD) and Zheng58 × parviglumis (ZP) populations, partly because the quantitative trait loci (QTLs) of FER resistance and yield traits were located close together. Four coincident QTLs (qFERbd5.177, qFERbd10.140, qFERzp4.066, and qFERzp5.116) and two highly reliable resistance-yield synergistic QTLs (qFERbd10.140 and qFERzp4.066) were identified in the BD and ZP populations, opening up the possibility of breeding for FER resistance without reducing yield.

Assessing the impacts of agricultural drought (SPI/SPEI) on maize and wheat yields across Hungary

This study examined the physical properties of agricultural drought (i.e., intensity, duration, and severity) in Hungary from 1961 to 2010 based on the Standardized Precipitation Index (SPI) and the Standardized Precipitation Evapotranspiration Index (SPEI). The study analyzed the interaction between drought and crop yield for maize and wheat using standardized yield residual series (SYRS), and the crop-drought resilient factor (CDRF). The results of both SPI and SPEI (-3, -6) showed that the western part of Hungary has significantly more prone to agricultural drought than the eastern part of the country. Drought frequency analysis reveals that the eastern, northern, and central parts of Hungary were the most affected regions. Drought analysis also showed that drought was particularly severe in Hungary during 1970–1973, 1990–1995, 2000–2003, and 2007. The yield of maize was more adversely affected than wheat especially in the western and southern regions of Hungary (1961–2010). In general, maize and wheat yields were severely non-resilient (CDRF < 0.8) in the central and western part of the country. The results suggest that drought events are a threat to the attainment of the second Sustainable Development Goals (SDG-2). Therefore, to ensure food security in Hungary and in other parts of the world, drought resistant crop varieties need to be developed to mitigate the adverse effects of climate change on agricultural production.

Maize dispersal patterns associated with different types of endosperm and migration of indigenous groups in lowland South America

BACKGROUND AND AIMS

The lowlands of South America appear to be remarkably important in the evolutionary history of maize, due to new evidence that suggests that maize dispersed from Mexico and arrived in this region in a state of partial domestication. This study aimed to identify dispersal patterns of maize genetic diversity in this part of the continent.

METHODS

A total of 170 maize accessions were characterized with 4398 single nucleotide polymorphisms (SNPs) and analysed to determine if maize dispersal was associated with types of endosperm and indigenous language families.

KEY RESULTS

Four genetic groups were identified in the discriminant analysis of principal components and five groups in the cluster analysis (neighbour-joining method). The groups were structured according to the predominance of endosperm types (popcorn, floury, flint/semi-flint). Spatial principal component analysis of genetic variation suggests different dispersal patterns for each endosperm type and can be associated with hypotheses of expansions of different indigenous groups.

CONCLUSIONS

From a possible origin in Southwestern Amazonia, different maize dispersal routes emerged: (1) towards Northern Amazonia, which continued towards the Caatinga and south-eastern Atlantic Forest (Floury); (2) towards Southern Brazil, passing through the Cerrado and Southern Atlantic Forest reaching the Pampa region (Floury); and (3) along the Atlantic Coast, following Tupi movements originating from two separate expansions: one (Tupinamba) from north to south, and the other (Guarani) in the opposite direction, from south to north (flint, floury and popcorn).

Demonstration of local adaptation in maize landraces by reciprocal transplantation

Populations are locally adapted when they exhibit higher fitness than foreign populations in their native habitat. Maize landrace adaptations to highland and lowland conditions are of interest to researchers and breeders. To determine the prevalence and strength of local adaptation in maize landraces, we performed a reciprocal transplant experiment across an elevational gradient in Mexico. We grew 120 landraces, grouped into four populations (Mexican Highland, Mexican Lowland, South American Highland, South American Lowland), in Mexican highland and lowland common gardens and collected phenotypes relevant to fitness and known highland-adaptive traits such as anthocyanin pigmentation and macrohair density. 67k DArTseq markers were generated from field specimens to allow comparisons between phenotypic patterns and population genetic structure. We found phenotypic patterns consistent with local adaptation, though these patterns differ between the Mexican and South American populations. Quantitative trait differentiation (Q ST) was greater than neutral allele frequency differentiation (F ST) for many traits, signaling directional selection between pairs of populations. All populations exhibited higher fitness metric values when grown at their native elevation, and Mexican landraces had higher fitness than South American landraces when grown in these Mexican sites. As environmental distance between landraces' native collection sites and common garden sites increased, fitness values dropped, suggesting landraces are adapted to environmental conditions at their natal sites. Correlations between fitness and anthocyanin pigmentation and macrohair traits were stronger in the highland site than the lowland site, supporting their status as highland-adaptive. These results give substance to the long-held presumption of local adaptation of New World maize landraces to elevation and other environmental variables across North and South America.

Genetic diversity, distribution and domestication history of the neglected GGAtAt genepool of wheat

We present a comprehensive survey of cytogenetic and genomic diversity of the GGA^tA^t genepool of wheat, thereby unlocking these plant genetic resources for wheat improvement. Wheat yields are stagnating around the world and new sources of genes for resistance or tolerances to abiotic traits are required. In this context, the tetraploid wheat wild relatives are among the key candidates for wheat improvement. Despite its potential huge value for wheat breeding, the tetraploid GGA^tA^t genepool is largely neglected. Understanding the population structure, native distribution range, intraspecific variation of the entire tetraploid GGA^tA^t genepool and its domestication history would further its use for wheat improvement. The paper provides the first comprehensive survey of genomic and cytogenetic diversity sampling the full breadth and depth of the tetraploid GGA^tA^t genepool. According to the results obtained, the extant GGA^tA^t genepool consists of three distinct lineages. We provide detailed insights into the cytogenetic composition of GGA^tA^t wheats, revealed group- and population-specific markers and show that chromosomal rearrangements play an important role in intraspecific diversity of T. araraticum. The origin and domestication history of the GGA^tA^t lineages is discussed in the context of state-of-the-art archaeobotanical finds. We shed new light on the complex evolutionary history of the GGA^tA^t wheat genepool and provide the basis for an increased use of the GGA^tA^t wheat genepool for wheat improvement. The findings have implications for our understanding of the origins of agriculture in southwest Asia.

Genetic dissection of maize plant architecture using a novel nested association mapping population

The leaf angle (LA), plant height (PH), and ear height (EH) are key plant architectural traits influencing maize (Zea mays L.) yield. However, their genetic determinants have not yet been well-characterized. Here, we developed a maize advanced backcross-nested association mapping population in Henan Agricultural University (HNAU-NAM1) comprised of 1,625 BC₁ F₄ /BC₂ F₄ lines. These were obtained by crossing a diverse set of 12 representative inbred lines with the common GEMS41 line, which were then genotyped using the MaizeSNP9.4K array. Genetic diversity and phenotypic distribution analyses showed considerable levels of genetic variation. We obtained 18-88 quantitative trait loci (QTLs) associated with LA, PH, and EH by using three complementary mapping methods, named as separate linkage mapping, joint linkage mapping, and genome-wide association studies. Our analyses enabled the identification of ten QTL hot-spot regions associated with the three traits, which were distributed on nine different chromosomes. We further selected 13 major QTLs that were simultaneously detected by three methods and deduced the candidate genes, of which eight were not reported before. The newly constructed HNAU-NAM1 population in this study will further broaden our insights into understanding of genetic regulation of plant architecture, thus will help to improve maize yield and provide an invaluable resource for maize functional genomics and breeding research.

Ancient Artworks and Crocus Genetics Both Support Saffron's Origin in Early Greece

Saffron crocus (Crocus sativus) is a male-sterile, triploid flower crop, and source of the spice and colorant saffron. For over three millennia, it was cultivated across the Mediterranean, including ancient Greece, Persia, and other cultures, later spreading all over the world. Despite saffron crocus' early omnipresence, its origin has been the matter of a century-old debate, in terms of area and time as well as parental species contribution. While remnants of the ancient arts, crafts, and texts still provide hints on its origin, modern genetics has the potential to efficiently follow these leads, thus shedding light on new possible lines of descent. In this review, we follow ancient arts and recent genetics to trace the evolutionary origin of saffron crocus. We focus on the place and time of saffron domestication and cultivation, and address its presumed autopolyploid origin involving cytotypes of wild Crocus cartwrightianus. Both ancient arts from Greece, Iran, and Mesopotamia as well as recent cytogenetic and comparative next-generation sequencing approaches point to saffron's emergence and domestication in ancient Greece, showing how both disciplines converge in tracing its origin.

Genetic Architecture of Grain Yield-Related Traits in Sorghum and Maize

Grain size, grain number per panicle, and grain weight are crucial determinants of yield-related traits in cereals. Understanding the genetic basis of grain yield-related traits has been the main research object and nodal in crop science. Sorghum and maize, as very close C4 crops with high photosynthetic rates, stress tolerance and large biomass characteristics, are extensively used to produce food, feed, and biofuels worldwide. In this review, we comprehensively summarize a large number of quantitative trait loci (QTLs) associated with grain yield in sorghum and maize. We placed great emphasis on discussing 22 fine-mapped QTLs and 30 functionally characterized genes, which greatly hinders our deep understanding at the molecular mechanism level. This review provides a general overview of the comprehensive findings on grain yield QTLs and discusses the emerging trend in molecular marker-assisted breeding with these QTLs.

Population genomics of Zea species identifies selection signatures during maize domestication and adaptation

BACKGROUND

Maize (Zea mays L. ssp. mays) was domesticated from teosinte (Zea mays ssp. parviglumis) about 9000 years ago in southwestern Mexico and adapted to a range of environments worldwide. Researchers have depicted the maize domestication and adaptation processes over the past two decades, but efforts have been limited either in sample size or genetic diversity. To better understand these processes, we conducted a genome-wide survey of 982 maize inbred lines and 190 teosinte accessions using over 40,000 single-nucleotide polymorphism markers.

RESULTS

Population structure, principal component analysis, and phylogenetic trees all confirmed the evolutionary relationship between maize and teosinte, and determined the evolutionary lineage of all species within teosinte. Shared haplotype analysis showed similar levels of ancestral alleles from Zea mays ssp. parviglumis and Zea mays ssp. mexicana in maize. Scans for selection signatures identified 394 domestication sweeps by comparing wild and cultivated maize and 360 adaptation sweeps by comparing tropical and temperate maize. Permutation tests revealed that the public association signals for flowering time were highly enriched in the domestication and adaptation sweeps. Genome-wide association study identified 125 loci significantly associated with flowering-time traits, ten of which identified candidate genes that have undergone selection during maize adaptation.

CONCLUSIONS

In this study, we characterized the history of maize domestication and adaptation at the population genomic level and identified hundreds of domestication and adaptation sweeps. This study extends the molecular mechanism of maize domestication and adaptation, and provides resources for basic research and genetic improvement in maize.

Proteomic analysis reveals the metabolic versatility of Amycolatopsis sp. BX17: A strain native from milpa agroecosystem soil

Amycolatopsis sp. BX17 is an actinobacterium isolated from milpa soils, which antagonizes the phytopathogenic fungus Fusarium graminearum. Metabolites secreted by the actinobacterium cultured in glucose-free medium inhibited 100% of the mycelial growth of F. graminearum RH1, while the inhibition rate was 65% in medium supplemented with 20 g/L glucose. With the aim of studying how the metabolism of strain BX17 is modulated by glucose as the main carbon source, media with 0 and 20 g/L glucose were selected to analyze the intracellular proteins by quantitative label-free proteomic analysis. Data are available via ProteomeXchange with identifier PXD028644. Proteins identified in bacteria cultured in medium without glucose were involved in glutamate metabolism, the Krebs cycle and the shikimate pathway, suggesting that amino acids are metabolized to synthesize antifungal compounds. In glucose-containing medium, carbon flux was directed mainly toward the synthesis of energy and cell growth. This study shows the metabolic versatility of Amycolatopsis BX17, and strengthens its potential use in designing biotechnological strategies for phytopathogen control. SIGNIFICANCE: Amycolatopsis BX17 is a bacterium isolated from milpa agroecosystems that antagonizes the phytopathogenic fungus Fusarium graminearum. Currently, there is scarce information about the metabolism involved in the biosynthesis of antifungal agents by this genus. We used a label-free proteomic approach to identify the differences in metabolic routes for antifungal biosynthesis in Amycolatopsis BX17 grown in media with 0 and 20 g/L glucose. Taken together the results suggest that the BX17 strain could be synthesizing the antifungal metabolite(s) from the Shikimate pathway through the synthesis and degradation of the amino acid tyrosine, which is a known precursor of glycopeptides with antibiotic and antifungal activity. While the lower antifungal activity of the metabolites secreted by Amycolatopsis BX17 when grown in a medium with glucose as the main carbon source, may be correlated with a lower synthesis of antifungal compounds, due to the directing of carbon flux toward metabolic pathways involved with energy synthesis and cell growth. Likewise, it is possible that the bacteria synthesize other compounds with biological activity, such as glycopeptides with antibiotic activity. These findings are relevant because they represent the first stage to understand the metabolic regulation involved in the biosynthesis of antifungal metabolites by the genus Amycolatopsis. Finally, improving our understanding of the metabolic regulation involved in the biosynthesis of antifungal metabolites is essential to design of strategies in agricultural biotechnology for phytopathogen control.

Molecular Parallelism Underlies Convergent Highland Adaptation of Maize Landraces

Convergent phenotypic evolution provides some of the strongest evidence for adaptation. However, the extent to which recurrent phenotypic adaptation has arisen via parallelism at the molecular level remains unresolved, as does the evolutionary origin of alleles underlying such adaptation. Here, we investigate genetic mechanisms of convergent highland adaptation in maize landrace populations and evaluate the genetic sources of recurrently selected alleles. Population branch excess statistics reveal substantial evidence of parallel adaptation at the level of individual single-nucleotide polymorphism (SNPs), genes, and pathways in four independent highland maize populations. The majority of convergently selected SNPs originated via migration from a single population, most likely in the Mesoamerican highlands, while standing variation introduced by ancient gene flow was also a contributor. Polygenic adaptation analyses of quantitative traits reveal that alleles affecting flowering time are significantly associated with elevation, indicating the flowering time pathway was targeted by highland adaptation. In addition, repeatedly selected genes were significantly enriched in the flowering time pathway, indicating their significance in adapting to highland conditions. Overall, our study system represents a promising model to study convergent evolution in plants with potential applications to crop adaptation across environmental gradients.

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.

Advances in Genomics Approaches Shed Light on Crop Domestication

Crop domestication occurred ~10,000–12,000 years ago when humans shifted from a hunter–gatherer to an agrarian society. Crops were domesticated by selecting the traits in wild plant species that were suitable for human use. Research is crucial to elucidate the mechanisms and processes involved in modern crop improvement and breeding. Recent advances in genomics have revolutionized our understanding of crop domestication. In this review, we summarized cutting-edge crop domestication research by presenting its (1) methodologies, (2) current status, (3) applications, and (4) perspectives. Advanced genomics approaches have clarified crop domestication processes and mechanisms, and supported crop improvement.

The influence of maize genotype on the rhizosphere eukaryotic community

The microbiota colonizing the rhizosphere contributes to plant growth, productivity, carbon sequestration and phytoremediation. Several studies address plant-associated bacteria; however, few studies analyze the effect of plant genotype on the eukaryotic community. Here, we analyzed the eukaryotic composition of maize rhizosphere from three different plant landraces and one inbred line grown in the same soil (common garden approach). This experimental design, coupled with 18S rDNA gene amplicon sequencing, allowed us to test the influence of maize and its genotype on the rhizosphere's eukaryotic community. We found that plant growth modified the eukaryotic community in soil, as diversity comparisons between maize rhizosphere and unplanted soil revealed significantly different eukaryotic composition. Various genera of nematodes and fungi, predominantly bacterial feeding nematodes and mycorrhizal fungi among other taxa, were increased in the rhizosphere samples. We also observed that maize genotype differentially shaped the relative abundance of the following fungal families in the rhizosphere: Acaulosporaceae, Aspergillaceae, Chaetomiaceae, Claroideoglomeraceae, Corticiaceae, Mortierellaceae, Trichocomaceae and Trichomeriaceae. Thus, plant genotype has a selective influence on establishing fungal communities in the rhizosphere. This study emphasizes the importance of an integrated consideration of plant genetics for future agricultural applications of microbes to crops.

Single-molecule long-read sequencing reveals extensive genomic and transcriptomic variation between maize and its wild relative teosinte (Zea mays ssp. parviglumis)

Teosinte (Zea mays ssp. parviglumis), the wild progenitor of maize (Zea mays L.), is an important germplasm resource for improvement of modern maize lines. However, we have limited genetic and genomic information about teosinte and lack state-of-the-art tools to annotate transcriptomes assembled by single-molecule long-read sequencing without a reference genome. Here, we employed single-molecule long-read sequencing of cDNA libraries from five tissues of the teosinte inbred line TIL11 and identified 70,044 nonredundant transcript isoforms. We devised a state-of-the-art, machine learning-based bioinformatics pipeline DenovoAS_Finder to annotate the TIL11 transcriptome without a complete reference genome with an accuracy of up to 91%, providing a robust gene classifier of complex genomes. Additionally, we constructed a draft TIL11 genome with 16,633 high-quality contigs and a N50 of 112 kb by Nanopore sequencing. Genes from families that expanded from teosinte to maize were significantly enriched in the gene ontology (GO) term "RNA modification pathway" and had more transcript isoforms in TIL11 than in the maize inbred line B73. Genes showed collinearity between TIL11 and B73, and intergenic regions were extensively altered by transposable elements. Our study furthers the understanding of maize domestication and provides a resource for the utilization of wild germplasm in maize breeding.

Joint analysis of days to flowering reveals independent temperate adaptations in maize

Domesticates are an excellent model for understanding biological consequences of rapid climate change. Maize (Zea mays ssp. mays) was domesticated from a tropical grass yet is widespread across temperate regions today. We investigate the biological basis of temperate adaptation in diverse structured nested association mapping (NAM) populations from China, Europe (Dent and Flint) and the United States as well as in the Ames inbred diversity panel, using days to flowering as a proxy. Using cross-population prediction, where high prediction accuracy derives from overall genomic relatedness, shared genetic architecture, and sufficient diversity in the training population, we identify patterns in predictive ability across the five populations. To identify the source of temperate adapted alleles in these populations, we predict top associated genome-wide association study (GWAS) identified loci in a Random Forest Classifier using independent temperate-tropical North American populations based on lines selected from Hapmap3 as predictors. We find that North American populations are well predicted (AUC equals 0.89 and 0.85 for Ames and USNAM, respectively), European populations somewhat well predicted (AUC equals 0.59 and 0.67 for the Dent and Flint panels, respectively) and that the Chinese population is not predicted well at all (AUC is 0.47), suggesting an independent adaptation process for early flowering in China. Multiple adaptations for the complex trait days to flowering in maize provide hope for similar natural systems under climate change.

Genetic analysis of three maize husk traits by QTL mapping in a maize-teosinte population

BACKGROUND

The maize husk consists of numerous leafy layers and plays vital roles in protecting the ear from pathogen infection and dehydration. Teosinte, the wild ancestor of maize, has about three layers of small husk outer covering the ear. Although several quantitative trait loci (QTL) underlying husk morphology variation have been reported, the genetic basis of husk traits between teosinte and maize remains unclear.

RESULTS

A linkage population including 191 BC₂F₈ inbred lines generated from the maize line Mo17 and the teosinte line X26-4 was used to identify QTL associated with three husk traits: i.e., husk length (HL), husk width (HW) and the number of husk layers (HN). The best linear unbiased predictor (BLUP) depicted wide phenotypic variation and high heritability of all three traits. The HL exhibited greater correlation with HW than HN. A total of 4 QTLs were identified including 1, 1, 2, which are associated with HL, HW and HN, respectively. The proportion of phenotypic variation explained by these QTLs was 9.6, 8.9 and 8.1% for HL, HN and HW, respectively.

CONCLUSIONS

The QTLs identified in this study will pave a path to explore candidate genes regulating husk growth and development, and benefit the molecular breeding program based on molecular marker-assisted selection to cultivate maize varieties with an ideal husk morphology.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

Flowering time regulation model revisited by pooled sequencing of mass selection populations

Maize is one of the most broadly cultivated crops throughout the world, and flowering time is a major adaptive trait for its diffusion. The biggest challenge in understanding maize flowering genetic architecture is that the trait is confounded with population structure. To eliminate the effect, we revisited the flower time genetic network by using a tropical maize population Pop32, which was under mass selection for adaptation to early flowering time in China for six generations from tropical to temperate regions. The days to anthesis (DTA) of the initial (Pop32C0), intermedia (Pop32C3), and final population (Pop32C5) was 90.77, 84.63, and 79.72 days on average, respectively. To examine the genetic mechanism and identify the genetic loci underlying this rapid change in flowering time of Pop32, we bulked 30 individuals from C0, C3, and C5 to conduct the whole genome sequencing. And we finally identified 4,973,810 high-quality single nucleotide polymorphisms (SNPs) and 6,517 genes with allele frequency significantly changed during the artificial improvement process. We speculate that these genes might participate in the adaptive improvement process and control flowering time. To identify the candidate genes for flowering time from the gene set with allele frequency changed, we carried out weighted gene co-expression network analysis (WGCNA), and identified four co-expression modules that highly associated with the flowering time development, as well as constructed the co-expression network of key flowering time genes. Gene Ontology (GO) enrichment analysis revealed that the GO terms photosynthesis/light reaction, carbohydrate binding, auxin mediated signaling pathway, response to temperature stimulus that are closely connected with flowering time. Furthermore, targeted GWAS revealed the genes are significantly connected with the flowering time. qRT-PCR of four candidate genes GRMZM2G019879, GRMZM2G055905, GRMZM2G058158, and GRMZM2G171365 showed that their expression level is similar to the flowering time genes, which playing a key role in maize flowering time transition. This study revealed that the changes of flowering time in mass selection process may be strongly associated with the variations of allele frequency changes, and we identified some important candidate genes for flowering time, which will provide a new insight for the rapid improvement of maize important agronomic traits and promote the gene cloning of maize flowering time.

Deciphering the Genetic Diversity of Landraces With High-Throughput SNP Genotyping of DNA Bulks: Methodology and Application to the Maize 50k Array

Genebanks harbor original landraces carrying many original favorable alleles for mitigating biotic and abiotic stresses. Their genetic diversity remains, however, poorly characterized due to their large within genetic diversity. We developed a high-throughput, cheap and labor saving DNA bulk approach based on single-nucleotide polymorphism (SNP) Illumina Infinium HD array to genotype landraces. Samples were gathered for each landrace by mixing equal weights from young leaves, from which DNA was extracted. We then estimated allelic frequencies in each DNA bulk based on fluorescent intensity ratio (FIR) between two alleles at each SNP using a two step-approach. We first tested either whether the DNA bulk was monomorphic or polymorphic according to the two FIR distributions of individuals homozygous for allele A or B, respectively. If the DNA bulk was polymorphic, we estimated its allelic frequency by using a predictive equation calibrated on FIR from DNA bulks with known allelic frequencies. Our approach: (i) gives accurate allelic frequency estimations that are highly reproducible across laboratories, (ii) protects against false detection of allele fixation within landraces. We estimated allelic frequencies of 23,412 SNPs in 156 landraces representing American and European maize diversity. Modified Roger's genetic Distance between 156 landraces estimated from 23,412 SNPs and 17 simple sequence repeats using the same DNA bulks were highly correlated, suggesting that the ascertainment bias is low. Our approach is affordable, easy to implement and does not require specific bioinformatics support and laboratory equipment, and therefore should be highly relevant for large-scale characterization of genebanks for a wide range of species.

Archaeological Central American maize genomes suggest ancient gene flow from South America

Maize is a global food staple with great economic and cultural importance. Archaeogenomic studies have revealed a process of protracted maize domestication and multiple waves of human-mediated dispersal in the Americas. Maize first arrived in South America as a partial domesticate, where the domestication syndrome became independently fixed and improved varieties developed away from the influence of wild gene flow. We demonstrate that hybrids of some of these improved varieties were likely reintroduced back to Central America. We hypothesize that this backflow of South American genetic material may have contributed to the development of a more productive staple, which was related to the growth and aggregation of human populations, and the formation of more complex social and political structures regionally.

Maize (Zea mays ssp. mays) domestication began in southwestern Mexico ∼9,000 calendar years before present (cal. BP) and humans dispersed this important grain to South America by at least 7,000 cal. BP as a partial domesticate. South America served as a secondary improvement center where the domestication syndrome became fixed and new lineages emerged in parallel with similar processes in Mesoamerica. Later, Indigenous cultivators carried a second major wave of maize southward from Mesoamerica, but it has been unclear until now whether the deeply divergent maize lineages underwent any subsequent gene flow between these regions. Here we report ancient maize genomes (2,300–1,900 cal. BP) from El Gigante rock shelter, Honduras, that are closely related to ancient and modern maize from South America. Our findings suggest that the second wave of maize brought into South America hybridized with long-established landraces from the first wave, and that some of the resulting newly admixed lineages were then reintroduced to Central America. Direct radiocarbon dates and cob morphological data from the rock shelter suggest that more productive maize varieties developed between 4,300 and 2,500 cal. BP. We hypothesize that the influx of maize from South America into Central America may have been an important source of genetic diversity as maize was becoming a staple grain in Central and Mesoamerica.

Gene regulatory effects of a large chromosomal inversion in highland maize

Chromosomal inversions play an important role in local adaptation. Inversions can capture multiple locally adaptive functional variants in a linked block by repressing recombination. However, this recombination suppression makes it difficult to identify the genetic mechanisms underlying an inversion's role in adaptation. In this study, we used large-scale transcriptomic data to dissect the functional importance of a 13 Mb inversion locus (Inv4m) found almost exclusively in highland populations of maize (Zea mays ssp. mays). Inv4m was introgressed into highland maize from the wild relative Zea mays ssp. mexicana, also present in the highlands of Mexico, and is thought to be important for the adaptation of these populations to cultivation in highland environments. However, the specific genetic variants and traits that underlie this adaptation are not known. We created two families segregating for the standard and inverted haplotypes of Inv4m in a common genetic background and measured gene expression effects associated with the inversion across 9 tissues in two experimental conditions. With these data, we quantified both the global transcriptomic effects of the highland Inv4m haplotype, and the local cis-regulatory variation present within the locus. We found diverse physiological effects of Inv4m across the 9 tissues, including a strong effect on the expression of genes involved in photosynthesis and chloroplast physiology. Although we could not confidently identify the causal alleles within Inv4m, this research accelerates progress towards understanding this inversion and will guide future research on these important genomic features.

A Tale of Maize, Palm, and Pine: Changing Socio-Ecological Interactions from Pre-Classic Maya to the Present Day in Belize

The environmental impact of the ancient Maya, and subsequent ecological recovery following the Terminal Classic decline, have been the key foci of research into socio-ecological interactions in the Yucatán peninsula. These foci, however, belie the complex pattern of resource exploitation and agriculture associated with post-Classic Maya societies and European colonisation. We present a high-resolution, 1200-year record of pollen and charcoal data from a 52-cm short core extracted from New River Lagoon, near to the European settlement of Indian Church, northern Belize. This study complements and extends a previous 3500-year reconstruction of past environmental change, located 1-km north of the new record and adjacent to the ancient Maya site of Lamanai. This current study shows a mixed crop production and palm agroforestry management strategy of the ancient Maya, which corroborates previous evidence at Lamanai. Comparison of the two records suggests that core agricultural and agroforestry activities shifted southwards, away from the centre of Lamanai, beginning at the post-Classic period. The new record also demonstrates that significant changes in land-use were not associated with drought at the Terminal Classic (ca. CE 1000) or the European Encounter (ca. CE 1500), but instead resulted from social and cultural change in the post-Classic period (CE 1200) and new economies associated with the British timber trade (CE 1680). The changes in land-use documented in two adjacent records from the New River Lagoon underline the need to reconstruct human–environment interactions using multiple, spatially, and temporally diverse records.

Symbiotic interactions between chickpea (Cicer arietinum L.) genotypes and Mesorhizobium strains

Legume genotype (GL) x rhizobium genotype (GR) interaction in chickpea was studied using a genetically diverse set of accessions and rhizobium strains in modified Leonard Jars. A subset of effective GL x GR combinations was subsequently evaluated in a pot experiment to identify combinations of chickpea genotypes and rhizobium strains with stable and superior symbiotic performance. A linear mixed model was employed to analyse the occurrence of GL x GR interaction and an additive main effects and multiplicative interaction (AMMI) model was used to study patterns in the performance of genotype-strain combinations. We found statistically significant interaction in jars in terms of symbiotic effectiveness that was entirely due to the inclusion of one of the genotypes, ICC6263. No interaction was found in a subsequent pot experiment. The presence of two genetic groups (Kabuli and Desi genepools) did not affect interaction with Mesorhizobium strains. With the exception of a negative interaction with genotype ICC6263 in the jar experiment, the type strain Mesorhizobium ciceri LMG 14989 outperformed or equalled other strains on all chickpea genotypes in both jar and pot experiments. Similar to earlier reports in common bean, our results suggest that efforts to find more effective strains may be more rewarding than aiming for identification of superior combinations of strains and genotypes.