Identification and fine mapping of quantitative trait loci for the number of vascular bundle in maize stem

QTL mapping and BSR-seq revealed loci and candidate genes associated with the sporadic multifoliolate phenotype in soybean (Glycine max)

KEY MESSAGE

The QTLs and candidate genes governing the multifoliolate phenotype were identified by combining linkage mapping with BSR-seq, revealing a possible interplay between genetics and the environment in soybean leaf development. Soybean, as a legume, is typified by trifoliolate leaves. Although multifoliolate leaves (compound leaves with more than three leaflets each) have been reported in soybean, including sporadic appearances in the first compound leaves in a recombinant inbred line (RIL) population from a cross between cultivated soybean C08 and wild soybean W05 from this study, the genetic basis of this phenomenon is still unclear. Here, we integrated quantitative trait locus (QTL) mapping with bulked segregant RNA sequencing (BSR-seq) to identify the genetic loci associated with the multifoliolate phenotype in soybean. Using linkage mapping, ten QTLs related to the multifoliolate trait were identified. Among these, a significant and major QTL, qMF-2-1 on chromosome 2 and consistently detected across biological replicates, explained more than 10% of the phenotypic variation. Together with BSR-seq analyses, which analyzed the RILs with the highest multifoliolate frequencies and those with the lowest frequencies as two distinct bulks, two candidate genes were identified: Glyma.06G204300 encoding the transcription factor TCP5, and Glyma.06G204400 encoding LONGIFOLIA 2 (LNG2). Transcriptome analyses revealed that stress-responsive genes were significantly differentially expressed between high-multifoliolate occurrence lines and low occurrence ones, indicating environmental factors probably influence the appearance of multifoliolate leaves in soybean through stress-responsive genes. Hence, this study offers new insights into the genetic mechanism behind the multifoliolate phenotype in soybean.

ZmCCT10-relayed photoperiod sensitivity regulates natural variation in the arithmetical formation of male germinal cells in maize

Extensive mutational screening studies have documented genes regulating anther and pollen development. Knowledge concerning how formation of male germinal cell is arithmetically controlled in natural populations, under different environmental conditions, is lacking. We counted pollen number within a single anther and a maize-teosinte BC₂ S₃ recombinant inbred line population to identify ZmCCT10 as a major determinant of pollen number variation. ZmCCT10 was originally identified as a photoperiod-sensitive negative regulator of flowering. ZmCCT10 inactivation, after transposon insertion within its promoter, is proposed to have accelerated maize spread toward higher latitudes, thus allowing temperate maize to flower under long-day conditions. We showed that the active ZmCCT10 allele decreased pollen formation. As different active and inactive ZmCCT10 alleles have been found in natural maize populations, this represents the first report of a gene controlling pollen number in a crop natural population. These findings suggest that higher pollen number, which provides a competitive advantage in open-pollinated populations, may have been one of the major driving forces for the selection of an inactive ZmCCT10 allele during tropical maize domestication. We provide evidence that ZmCCT10 has opposite effects on cell proliferation of archesporial and tapetum cells and it modulates expression of key regulators during early anther development.

QTL analyses of soybean root system architecture revealed genetic relationships with shoot-related traits

The genetic basis of soybean root system architecture (RSA) and the genetic relationship between shoot and RSA were revealed by integrating data from recombinant inbred population grafting and QTL mapping. Variations in root system architecture (RSA) affect the functions of roots and thus play vital roles in plant adaptations and agricultural productivity. The aim of this study was to unravel the genetic relationship between RSA traits and shoot-related traits in soybean. This study characterized RSA variability at seedling stage in a recombinant inbred population, derived from a cross between cultivated soybean C08 and wild soybean W05, and performed high-resolution quantitative trait locus (QTL) mapping. In total, 34 and 41 QTLs were detected for RSA-related and shoot-related traits, respectively, constituting eight QTL clusters. Significant QTL correspondence was found between shoot biomass and RSA-related traits, consistent with significant correlations between these phenotypes. RSA-related QTLs also overlapped with selection regions in the genome, suggesting the cultivar RSA could be a partial consequence of domestication. Using reciprocal grafting, we confirmed that shoot-derived signals affected root development and the effects were controlled by multiple loci. Meanwhile, RSA-related QTLs were found to co-localize with four soybean flowering-time loci. Consistent with the phenotypes of the parental lines of our RI population, diminishing the function of flowering controlling E1 family through RNA interference (RNAi) led to reduced root growth. This implies that the flowering time-related genes within the RSA-related QTLs are actually contributing to RSA. To conclude, this study identified the QTLs that determine RSA through controlling root growth indirectly via regulating shoot functions, and discovered superior alleles from wild soybean that could be used to improve the root structure in existing soybean cultivars.

Identification of QTL underlying the main stem related traits in a doubled haploid barley population

Lodging reduces grain yield in cereal crops. The height, diameter and strength of stem are crucial for lodging resistance, grain yield, and photosynthate transport in barley. Understanding the genetic basis of stem benefits barley breeding. Here, we evaluated 13 stem related traits after 28 days of heading in a barley DH population in two consecutive years. Significant phenotypic correlations between lodging index (LI) and other stem traits were observed. Three mapping methods using the experimental data and the BLUP data, detected 27 stable and major QTLs, and 22 QTL clustered regions. Many QTLs were consistent with previously reported traits for grain filling rate, internodes, panicle and lodging resistance. Further, candidate genes were predicted for stable and major QTLs and were associated with plant development and adverse stress in the transition from vegetative stage to reproductive stage. This study provided potential genetic basis and new information for exploring barley stem morphology, and laid a foundation for map-based cloning and further fine mapping of these QTLs.

Comparative transcriptomic reveals the molecular mechanism of maize hybrid Zhengdan538 in response to water deficit

Heterosis, known as one of the most successful strategies for increasing grain yield and abiotic/biotic stress tolerance, has been widely exploited in maize breeding. However, the underlying molecular processes are still to be elucidated. The maize hybrid "Zhengdan538" shows high tolerance to drought stress. The transcriptomes of the seedling leaves of its parents, "ZhengA88" and "ZhengT22" and their reciprocal F₁ hybrid under well-watered and water deficit conditions, were analyzed by RNA sequencing (RNA-Seq). Transcriptome profiling of the reciprocal hybrid revealed 2994-4692 differentially expressed genes (DEGs) under well-watered and water-deficit conditions, which were identified by comparing with their parents. The reciprocal hybrid was more closely related to the parental line "ZhengT22" than to the parental line "ZhengA88" in terms of gene expression patterns under water-deficit condition. Furthermore, genes showed expression level dominance (ELD), especially the high-parental ELD (Class 3 and 5), accounted for the largest proportion of DEGs between the reciprocal F₁ hybrid and their parental lines under water deficit. These ELD genes mainly participated in photosynthesis, energy biosynthesis, and metabolism processes. The results indicated that ELD genes played important roles in hybrid tolerance to water deficit. Moreover, a set of important drought-responsive transcription factors were found to be encoded by the identified ELD genes and are thought to function in improving drought tolerance in maize hybrid plants. Our results provide a better understanding of the molecular mechanism of drought tolerance in hybrid maize.

Dissection of the genetic architecture of peduncle vascular bundle-related traits in maize by a genome-wide association study

The peduncle vascular system of maize is critical for the transport of photosynthetic products, nutrients, and water from the roots and leaves to the ear. Accordingly, it positively affects the grain yield. However, the genetic basis of peduncle vascular bundle (PVB)-related traits in maize remains unknown. Thus, 15 PVB-related traits of 386 maize inbred lines were investigated at three locations (Yongcheng, 17YC; Kaifeng, 20KF; and Yuanyang, 20YY). The repeatability for the 15 traits ranged from 35.53% to 92.13%. A genome-wide association study was performed and 69 non-redundant quantitative trait loci (QTL) were detected, including 9, 41, and 27 QTL identified at 17YC, 20KF, and 20YY, respectively. These QTL jointly explained 4.72% (SLL) to 37.30% (NSVB) of the phenotypic variation. Eight QTL were associated with the same trait at two locations. Furthermore, four pleiotropic QTL were identified. Moreover, one QTL (qPVB44), associated with NSVB_20KF, was co-localized with a previously reported locus related to kernel width, implying qPVB44 may affect the kernel width by modulating the number of small vascular bundles. Examinations of the 69 QTL identified 348 candidate genes that were classified in five groups. Additionally, 26 known VB-related homologous genes (e.g. VLN2, KNOX1, and UGT72B3) were detected in 20 of the 69 QTL. A comparison of the NSVB between a Zmvln2 EMS mutant and its wild type elucidated the function of the candidate gene ZmVLN2. These results are important for clarifying the genetic basis of PVB-related traits and may be useful for breeding new high-yielding maize cultivars.

High-Throughput Phenotyping Accelerates the Dissection of the Phenotypic Variation and Genetic Architecture of Shank Vascular Bundles in Maize (Zea mays L.)

The vascular bundle of the shank is an important ‘flow’ organ for transforming maize biological yield to grain yield, and its microscopic phenotypic characteristics and genetic analysis are of great significance for promoting the breeding of new varieties with high yield and good quality. In this study, shank CT images were obtained using the standard process for stem micro-CT data acquisition at resolutions up to 13.5 μm. Moreover, five categories and 36 phenotypic traits of the shank including related to the cross-section, epidermis zone, periphery zone, inner zone and vascular bundle were analyzed through an automatic CT image process pipeline based on the functional zones. Next, we analyzed the phenotypic variations in vascular bundles at the base of the shank among a group of 202 inbred lines based on comprehensive phenotypic information for two environments. It was found that the number of vascular bundles in the inner zone (IZ_VB_N) and the area of the inner zone (IZ_A) varied the most among the different subgroups. Combined with genome-wide association studies (GWAS), 806 significant single nucleotide polymorphisms (SNPs) were identified, and 1245 unique candidate genes for 30 key traits were detected, including the total area of vascular bundles (VB_A), the total number of vascular bundles (VB_N), the density of the vascular bundles (VB_D), etc. These candidate genes encode proteins involved in lignin, cellulose synthesis, transcription factors, material transportation and plant development. The results presented here will improve the understanding of the phenotypic traits of maize shank and provide an important phenotypic basis for high-throughput identification of vascular bundle functional genes of maize shank and promoting the breeding of new varieties with high yield and good quality.

Genome-Wide Association Study of Vascular Bundle-Related Traits in Maize Stalk

The vascular bundle plays an important role in nutrient transportation in plants and exerts great influence on crop yield. Maize is widely used for food, feed, and fuel, producing the largest yield in the world. However, genes and molecular mechanism controlling vascular bundle-related traits in maize have largely remained undiscovered. In this study, a natural population containing 248 diverse maize inbred lines genotyped with high-throughput SNP markers was used for genome-wide association study. The results showed that broad variations existed for the vascular bundle-related traits which are subject to genetic structure and it was suitable for association analysis. In this study, we identified 15, 13, 2, 1, and 5 SNPs significantly associated with number of small vascular bundle, number of large vascular bundle, average area of single small vascular bundle, average area of single large vascular bundle, and cross-sectional area, respectively. The 210 candidate genes in the confidence interval can be classified into ten biological processes, three cellular components, and eight molecular functions. As for the Kyoto Encyclopedia of Genes and Genomes analysis of the candidate genes, a total of six pathways were identified. Finally, we found five genes related to vascular development, three genes related to cell wall, and two genes related to the mechanical strength of the stalk. Our results provide the further understanding of the genetic foundation of vascular bundle-related traits in maize stalk.

Natural Variation in Crops: Realized Understanding, Continuing Promise

Crops feed the world's population and shape human civilization. The improvement of crop productivity has been ongoing for almost 10,000 years and has evolved from an experience-based to a knowledge-driven practice over the past three decades. Natural alleles and their reshuffling are long-standing genetic changes that affect how crops respond to various environmental conditions and agricultural practices. Decoding the genetic basis of natural variation is central to understanding crop evolution and, in turn, improving crop breeding. Here, we review current advances in the approaches used to map the causal alleles of natural variation, provide refined insights into the genetics and evolution of natural variation, and outline how this knowledge promises to drive the development of sustainable agriculture under the dome of emerging technologies.

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.

Responses of Maize Internode to Water Deficit Are Different at the Biochemical and Histological Levels

Maize feeding value is strongly linked to plant digestibility. Cell wall composition and structure can partly explain cell wall digestibility variations, and we recently showed that tissue lignification and lignin spatial distribution also contribute to cell wall digestibility variations. Although the genetic determinism of digestibility and cell wall composition has been studied for more than 20 years, little is available concerning that of tissue lignification. Moreover, maize yield is negatively impacted by water deficit, and we newly highlighted the impact of water deficit on cell wall digestibility and composition together with tissue lignification. Consequently, the aim of this study was to explore the genetic mechanisms of lignin distribution in link with cell wall composition and digestibility under contrasted water regimes. Maize internodes from a recombinant inbred line (RIL) population grown in field trials with contrasting irrigation scenarios were biochemically and histologically quantified. Results obtained showed that biochemical and histological traits have different response thresholds to water deficit. Histological profiles were therefore only modified under pronounced water deficit, while most of the biochemical traits responded whatever the strength of the water deficit. Three main clusters of quantitative trait locus (QTL) for histological traits were detected. Interestingly, overlap between the biochemical and histological clusters is rare, and one noted especially colocalizations between histological QTL/clusters and QTL for p-coumaric acid content. These findings reinforce the suspected role of tissue p-coumaroylation for both the agronomic properties of plants as well as their digestibility.

Dissecting the phenotypic components and genetic architecture of maize stem vascular bundles using high-throughput phenotypic analysis

High-throughput phenotyping is increasingly becoming an important tool for rapid advancement of genetic gain in breeding programmes. Manual phenotyping of vascular bundles is tedious and time-consuming, which lags behind the rapid development of functional genomics in maize. More robust and automated techniques of phenotyping vascular bundles traits at high-throughput are urgently needed for large crop populations. In this study, we developed a standard process for stem micro-CT data acquisition and an automatic CT image process pipeline to obtain vascular bundle traits of stems including geometry-related, morphology-related and distribution-related traits. Next, we analysed the phenotypic variation of stem vascular bundles between natural population subgroup (480 inbred lines) based on 48 comprehensively phenotypic information. Also, the first database for stem micro-phenotypes, MaizeSPD, was established, storing 554 pieces of basic information of maize inbred lines, 523 pieces of experimental information, 1008 pieces of CT scanning images and processed images, and 24 192 pieces of phenotypic data. Combined with genome-wide association studies (GWASs), a total of 1562 significant single nucleotide polymorphism (SNPs) were identified for 30 stem micro-phenotypic traits, and 84 unique genes of 20 traits such as VBNum, VBAvArea and PZVBDensity were detected. Candidate genes identified by GWAS mainly encode enzymes involved in cell wall metabolism, transcription factors, protein kinase and protein related to plant signal transduction and stress response. The results presented here will advance our knowledge about phenotypic trait components of stem vascular bundles and provide useful information for understanding the genetic controls of vascular bundle formation and development.

QTG-Finder2: A Generalized Machine-Learning Algorithm for Prioritizing QTL Causal Genes in Plants

Linkage mapping has been widely used to identify quantitative trait loci (QTL) in many plants and usually requires a time-consuming and labor-intensive fine mapping process to find the causal gene underlying the QTL. Previously, we described QTG-Finder, a machine-learning algorithm to rationally prioritize candidate causal genes in QTLs. While it showed good performance, QTG-Finder could only be used in Arabidopsis and rice because of the limited number of known causal genes in other species. Here we tested the feasibility of enabling QTG-Finder to work on species that have few or no known causal genes by using orthologs of known causal genes as the training set. The model trained with orthologs could recall about 64% of Arabidopsis and 83% of rice causal genes when the top 20% ranked genes were considered, which is similar to the performance of models trained with known causal genes. The average precision was 0.027 for Arabidopsis and 0.029 for rice. We further extended the algorithm to include polymorphisms in conserved non-coding sequences and gene presence/absence variation as additional features. Using this algorithm, QTG-Finder2, we trained and cross-validated Sorghum bicolor and Setaria viridis models. The S. bicolor model was validated by causal genes curated from the literature and could recall 70% of causal genes when the top 20% ranked genes were considered. In addition, we applied the S. viridis model and public transcriptome data to prioritize a plant height QTL and identified 13 candidate genes. QTL-Finder2 can accelerate the discovery of causal genes in any plant species and facilitate agricultural trait improvement.

Field Screening of Rice Germplasm (Oryza sativa L. ssp. japonica) Based on Days to Flowering for Drought Escape

Terminal drought stress is one of the restrictive factors in rice production and is expected to upsurge under the current situation of climate change. The study evaluated the performance of 2030 rice genotypes under continuous drought stress conditions based on days to flowering (DF). The genotypes under augmented randomized complete block design were sown in May/June of 2017 and 2018 in the field with movable rainout that resulted in huge genetic diversity among the accessions. Descriptive statistics confirmed clear variation among accessions on growth duration, plant height to leaf, plant height to panicle, and germination percentage. Correlation, chemometric, and agglomerative hierarchical cluster analyses were performed that categorized the germplasm into 10 groups. Genotypes in clusters VIII and IX (drought-resistant) revealed better agronomic performance in terms of reduced days to flowering, but conversely taller plant height and higher maturity (%) under severe stress. Genotypes in clusters IV, V, and X were discovered to be drought-susceptible. The screened genotypes like Longjing 12, Longdun 102, Yanjing 22, Liaojing 27, Xiaohongbandao, Songjing 17, and Zaoshuqingsen can be utilized in rice breeding improvement programs for drought tolerance in terms of severe continuous drought, as well as terminal drought stress.

TeoNAM: A Nested Association Mapping Population for Domestication and Agronomic Trait Analysis in Maize

Recombinant inbred lines (RILs) are an important resource for mapping genes controlling complex traits in many species. While RIL populations have been developed for maize, a maize RIL population with multiple teosinte inbred lines as parents has been lacking. Here, we report a teosinte nested association mapping (TeoNAM) population, derived from crossing five teosinte inbreds to the maize inbred line W22. The resulting 1257 BC1S4 RILs were genotyped with 51,544 SNPs, providing a high-density genetic map with a length of 1540 cM. On average, each RIL is 15% homozygous teosinte and 8% heterozygous. We performed joint linkage mapping (JLM) and a genome-wide association study (GWAS) for 22 domestication and agronomic traits. A total of 255 QTL from JLM were identified, with many of these mapping near known genes or novel candidate genes. TeoNAM is a useful resource for QTL mapping for the discovery of novel allelic variation from teosinte. TeoNAM provides the first report that PROSTRATE GROWTH1, a rice domestication gene, is also a QTL associated with tillering in teosinte and maize. We detected multiple QTL for flowering time and other traits for which the teosinte allele contributes to a more maize-like phenotype. Such QTL could be valuable in maize improvement.

Evolutionary Metabolomics Identifies Substantial Metabolic Divergence between Maize and Its Wild Ancestor, Teosinte

Maize (Zea mays subsp mays) was domesticated from its wild ancestor, teosinte (Zea mays subsp parviglumis). Maize's distinct morphology and adaptation to diverse environments required coordinated changes in various metabolic pathways. However, how the metabolome was reshaped since domestication remains poorly understood. Here, we report a comprehensive assessment of divergence in the seedling metabolome between maize and teosinte. In total, 461 metabolites exhibited significant divergence due to selection. Interestingly, teosinte and tropical and temperate maize, representing major stages of maize evolution, targeted distinct sets of metabolites. Alkaloids, terpenoids, and lipids were specifically targeted in the divergence between teosinte and tropical maize, while benzoxazinoids were specifically targeted in the divergence between tropical and temperate maize. To identify genetic factors controlling metabolic divergence, we assayed the seedling metabolome of a large maize-by-teosinte cross population. We show that the recent metabolic divergence between tropical and temperate maize tended to have simpler genetic architecture than the divergence between teosinte and tropical maize. Through integrating transcriptome data, we identified candidate genes contributing to metabolic divergence, many of which were under selection at the nucleotide and transcript levels. Through overexpression or mutant analysis, we verified the roles of Flavanone 3-hydroxylase1, Purple aleurone1, and maize terpene synthase1 in the divergence of their related biosynthesis pathways. Our findings not only provide important insights into domestication-associated changes in the metabolism but also highlight the power of combining omics data for trait dissection.

A domestication-associated reduction in K+ -preferring HKT transporter activity underlies maize shoot K+ accumulation and salt tolerance

Maize was domesticated from Balsas teosinte c. 10 000 yr ago. Previous studies have suggested that increased tolerance to environmental stress occurred during maize domestication. However, the underlying genetic basis remains largely unknown. We used a maize (W22)-teosinte recombinant inbred line (RIL) to investigate the salt wild-type tolerance aspects of maize domestication. We revealed that ZmHKT2 is a major QTL that regulates K⁺ homeostasis in saline soils. ZmHKT2 encodes a K⁺ -preferring HKT family transporter and probably reduces shoot K⁺ content by removing K⁺ ions from root-to-shoot flowing xylem sap, ZmHKT2 deficiency increases xylem sap and shoot K⁺ concentrations, and increases salt tolerance. A coding sequence polymorphism in the ZmHKT2^W22 allele (SNP389-G) confers an amino acid variant ZmHKT2 that increases xylem sap K⁺ concentration, thereby increasing shoot K⁺ content and salt tolerance. Additional analyses showed that SNP389-G first existed in teosinte (allele frequency 56% in assayed accessions), then swept through the maize population (allele frequency 98%), and that SNP389-G probably underwent positive selection during maize domestication. We conclude that a domestication-associated reduction in K⁺ transport activity in ZmHKT2 underlies maize shoot K⁺ content and salt tolerance, and propose that CRISPR-based editing of ZmHKT2 might provide a feasible strategy for improving maize salt tolerance.

ZmMADS69 functions as a flowering activator through the ZmRap2.7-ZCN8 regulatory module and contributes to maize flowering time adaptation

Flowering time is a major determinant of the local adaptation of plants. Although numerous loci affecting flowering time have been mapped in maize, their underlying molecular mechanisms and roles in adaptation remain largely unknown. Here, we report the identification and characterization of MADS-box transcription factor ZmMADS69 that functions as a flowering activator through the ZmRap2.7-ZCN8 regulatory module and contributes to adaptation. We show that ZmMADS69 underlies a quantitative trait locus controlling the difference in flowering time between maize and its wild ancestor, teosinte. Maize ZmMADS69 allele is expressed at a higher level at floral transition and confers earlier flowering than the teosinte allele under long days and short days. Overexpression of ZmMADS69 causes early flowering, while a transposon insertion mutant of ZmMADS69 exhibits delayed flowering. ZmMADS69 shows pleiotropic effects for multiple traits of agronomic importance. ZmMADS69 functions upstream of the flowering repressor ZmRap2.7 to downregulate its expression, thereby relieving the repression of the florigen gene ZCN8 and causing early flowering. Population genetic analyses showed that ZmMADS69 was a target of selection and may have played an important role as maize spread from the tropics to temperate zones. Our findings provide important insights into the regulation and adaptation of flowering time.

High-throughput micro-phenotyping measurements applied to assess stalk lodging in maize (Zea mays L.)

BACKGROUND

The biomechanical properties of maize stalks largely determine their lodging resistance, which affects crop yield per unit area. However, the quantitative and qualitative relationship between micro-phenotypes and the biomechanics of maize stalks is still under examined. In particular, the roles of the number, geometry, and distribution of vascular bundles of stalks in maize lodging resistance remain unclear. Research on these biomechanical properties will benefit from high-resolution micro-phenotypic image acquisition capabilities, which have been improved by modern X-ray imaging devices such as micro-CT and the development of micro-phenotyping analysis software. Hence, high-throughput image analysis and accurate quantification of anatomical phenotypes of stalks are necessary.

RESULTS

We have updated VesselParser version 1.0 to version 2.0 and have improved its performance, accuracy, and computation strategies. Anatomical characteristics of the second and third stalk internodes of the cultivars 'Jingke968' and 'Jingdan38' were analyzed using VesselParser 2.0. The relationships between lodging resistance and anatomical phenotypes of stalks between the two different maize varieties were investigated. The total area of vascular bundles in the peripheral layer, auxiliary axis diameter, and total area of vascular bundles were revealed to have the highest correlation with mechanical properties, and anatomical phenotypes of maize stalk were better predictors of mechanical properties than macro features observed optically from direct measurement, such as diameter and perimeter.

CONCLUSIONS

This study demonstrates the utility of VesselParser 2.0 in assessing stalk mechanical properties. The combination of anatomical phenotypes and mechanical behavior research provides unique insights into the problem of stalk lodging, showing that micro phenotypes of vascular bundles are good predictors of maize stalk mechanical properties that may be important indices for the evaluation and identification of the biomechanical properties to improve lodging resistance of future maize varieties.

Stepwise cis-Regulatory Changes in ZCN8 Contribute to Maize Flowering-Time Adaptation

Maize (Zea mays ssp. mays) was domesticated in southwestern Mexico ∼9,000 years ago from its wild ancestor, teosinte (Zea mays ssp. parviglumis) [1]. From its center of origin, maize experienced a rapid range expansion and spread over 90° of latitude in the Americas [2-4], which required a novel flowering-time adaptation. ZEA CENTRORADIALIS 8 (ZCN8) is the maize florigen gene and has a central role in mediating flowering [5, 6]. Here, we show that ZCN8 underlies a major quantitative trait locus (QTL) (qDTA8) for flowering time that was consistently detected in multiple maize-teosinte experimental populations. Through association analysis in a large diverse panel of maize inbred lines, we identified a SNP (SNP-1245) in the ZCN8 promoter that showed the strongest association with flowering time. SNP-1245 co-segregated with qDTA8 in maize-teosinte mapping populations. We demonstrate that SNP-1245 is associated with differential binding by the flowering activator ZmMADS1. SNP-1245 was a target of selection during early domestication, which drove the pre-existing early flowering allele to near fixation in maize. Interestingly, we detected an independent association block upstream of SNP-1245, wherein the early flowering allele that most likely originated from Zea mays ssp. mexicana introgressed into the early flowering haplotype of SNP-1245 and contributed to maize adaptation to northern high latitudes. Our study demonstrates how independent cis-regulatory variants at a gene can be selected at different evolutionary times for local adaptation, highlighting how complex cis-regulatory control mechanisms evolve. Finally, we propose a polygenic map for the pre-Columbian spread of maize throughout the Americas.

Fine analysis of a genomic region involved in resistance to Mediterranean corn borer

BACKGROUND

Sesamia nonagrioides Lefebvere (Mediterranean corn borer, MCB) is the main pest of maize in the Mediterranean area. QTL for MCB stalk tunneling and grain yield under high MCB infestation had been located at bin 8.03-8.05 (4-21 cM and 10-30 cM respectively) in a previous analysis of the EP42 x EP39 RILs mapping population. The objective of the present work was to study with higher resolution those QTL, and validating and estimating with higher precision their locations and effects. To achieve this objective, we developed a set of 38 heterogeneous inbred families (HIFs) which were near-homozygous in the genome, except in the region under study. The HIFs were evaluated in multiple environments under artificial infestation with MCB and genotyped with SNPs.

RESULTS

The QTL for grain yield under high infestation was confirmed with higher precision and improved reliability at 112.6-116.9 Mb. On the contrary, the location of the QTL for stalk tunneling was not validated probably due to the fixation of some genomic regions during the development of the HIFs. Our study confirmed that the co-localization of the QTL for stalk tunneling and grain yield in the previous study was due to linked genes, not to pleiotropic effects. So, the QTL for grain yield can be used for improving grain yield without undesirable effect on stalk tunneling.

CONCLUSIONS

The HIF analysis is useful for validating QTL and for conducting deeper studies in traits related to corn borer resistance.

Complex genetic architecture underlies maize tassel domestication

Maize (Zea mays) tassels underwent profound morphological changes during maize domestication and improvement. Although a number of genes affecting maize inflorescence development have been identified, the genetic basis of the morphological changes in maize tassels since domestication is not well understood. Here, using a large population of 866 maize-teosinte BC₂ S₃ recombinant inbred lines genotyped using 19 838 single nucleotide polymorphism (SNP) markers, we performed high-resolution quantitative trait locus (QTL) mapping for five tassel morphological traits. We showed that the five tassel traits were associated with different genetic architecture features. Known genes for maize inflorescence development identified by mutagenesis were significantly enriched in the tassel trait QTLs, and many of these genes, including ramosa1 (ra1), barren inflorescence2 (bif2), unbranched2 (ub2), zea floricaula leafy2 (zfl2) and barren stalk fastigiate1 (baf1), showed evidence of selection. An in-depth nucleotide diversity analysis at the bif2 locus identified strong selection signatures in the 5'-regulatory region. We also found that several known flowering time genes co-localized with tassel trait QTLs. A further association analysis indicated that the maize photoperiod gene ZmCCT was significantly associated with tassel size variation. Using near-isogenic lines, we narrowed down a major-effect QTL for tassel length, qTL9-1, to a 513-kb physical region. These results provide important insights into the genetic architecture that controls maize tassel evolution.

The genetic architecture of leaf number and its genetic relationship to flowering time in maize

The number of leaves and their distributions on plants are critical factors determining plant architecture in maize (Zea mays), and leaf number is frequently used as a measure of flowering time, a trait that is key to local environmental adaptation. Here, using a large set of 866 maize-teosinte BC2 S3 recombinant inbred lines genotyped by using 19,838 single nucleotide polymorphism markers, we conducted a comprehensive genetic dissection to assess the genetic architecture of leaf number and its genetic relationship to flowering time. We demonstrated that the two components of total leaf number, the number of leaves above (LA) and below (LB) the primary ear, were under relatively independent genetic control and might be subject to differential directional selection during maize domestication and improvement. Furthermore, we revealed that flowering time and leaf number are commonly regulated at a moderate level. The pleiotropy of the genes ZCN8, dlf1 and ZmCCT on leaf number and flowering time were validated by near-isogenic line analysis. Through fine mapping, qLA1-1, a major-effect locus that specifically affects LA, was delimited to a region with severe recombination suppression derived from teosinte. This study provides important insights into the genetic basis of traits affecting plant architecture and adaptation. The genetic independence of LA from LB enables the optimization of leaf number for ideal plant architecture breeding in maize.

Modeling the Morphometric Evolution of the Maize Shoot Apical Meristem

The maize (Zea mays subsp. mays L.) shoot apical meristem (SAM) is a self-replenishing pool of stem cells that produces all above-ground plant tissues. Improvements in image acquisition and processing techniques have allowed high-throughput, quantitative genetic analyses of SAM morphology. As with other large-scale phenotyping efforts, meaningful descriptions of genetic architecture depend on the collection of relevant measures. In this study, we tested two quantitative image processing methods to describe SAM morphology within the genus Zea, represented by 33 wild relatives of maize and 841 lines from a domesticated maize by wild teosinte progenitor (MxT) backcross population, along with previously reported data from several hundred diverse maize inbred lines. Approximating the MxT SAM as a paraboloid derived eight parabolic estimators of SAM morphology that identified highly overlapping quantitative trait loci (QTL) on eight chromosomes, which implicated previously identified SAM morphology candidate genes along with new QTL for SAM morphological variation. Using a Fourier-transform related method of comprehensive shape analysis, we detected cryptic SAM shape variation that identified QTL on six chromosomes. We found that Fourier transform shape descriptors and parabolic estimation measures are highly correlated and identified similar QTL. Analysis of shoot apex contours from 73 anciently diverged plant taxa further suggested that parabolic shape may be a universal feature of plant SAMs, regardless of evolutionary clade. Future high-throughput examinations of SAM morphology may benefit from the ease of acquisition and phenotypic fidelity of modeling the SAM as a paraboloid.