Pleiotropic effects of the duplicate maize FLORICAULA/LEAFY genes zfl1 and zfl2 on traits under selection during maize domestication

Determination of single or paired-kernel-rows is controlled by two quantitative loci during maize domestication

KEY MESSAGE

qPEDS1, a major quantitative trait locus that determines kernel row number during domestication, harbors the proposed causal gene Zm00001d033675, which may affect jasmonic acid biosynthesis and determine the fate of spikelets. Maize domestication has achieved the production of maize with enlarged ears, enhancing grain productivity dramatically. Kernel row number (KRN), an important yield-related trait, has increased from two rows in teosinte to at least eight rows in modern maize. However, the genetic mechanisms underlying this process remain unclear. To understand KRN domestication, we developed a teosinte-maize BC2F7 population by introgressing teosinte into a maize background. We identified one line, Teosinte ear rank1 (Ter1), with only 5-7 kernel rows which is fewer than those in almost all maize inbred lines. We detected two quantitative trait loci underlying Ter1 and fine-mapped the major one to a 300-kb physical interval. Two candidate genes, Zm674 and Zm675, were identified from 26 maize reference genomes and teosinte bacterial artificial chromosome sequences. Finally, we proposed that Ter1 affects jasmonic acid biosynthesis in the developing ear to determine KRN by the fate of spikelets. This study provides novel insights into the genetic and molecular mechanisms underlying KRN domestication and candidates for de novo wild teosinte domestication.

Haplotype-based breeding: A new insight in crop improvement

Haplotype-based breeding (HBB) is one of the cutting-edge technologies in the realm of crop improvement due to the increasing availability of Single Nucleotide Polymorphisms identified by Next Generation Sequencing technologies. The complexity of the data can be decreased with fewer statistical tests and a lower probability of spurious associations by combining thousands of SNPs into a few hundred haplotype blocks. The presence of strong genomic regions in breeding lines of most crop species facilitates the use of haplotypes to improve the efficiency of genomic and marker-assisted selection. Haplotype-based breeding as a Genomic Assisted Breeding (GAB) approach harnesses the genome sequence data to pinpoint the allelic variation used to hasten the breeding cycle and circumvent the challenges associated with linkage drag. This review article demonstrates ways to identify candidate genes, superior haplotype identification, haplo-pheno analysis, and haplotype-based marker-assisted selection. The crop improvement strategies that utilize superior haplotypes will hasten the breeding progress to safeguard global food security.

The Function of Florigen in the Vegetative-to-Reproductive Phase Transition in and around the Shoot Apical Meristem

Plants undergo a series of developmental phases throughout their life-cycle, each characterized by specific processes. Three critical features distinguish these phases: the arrangement of primordia (phyllotaxis), the timing of their differentiation (plastochron) and the characteristics of the lateral organs and axillary meristems. Identifying the unique molecular features of each phase, determining the molecular triggers that cause transitions and understanding the molecular mechanisms underlying these transitions are keys to gleaning a complete understanding of plant development. During the vegetative phase, the shoot apical meristem (SAM) facilitates continuous leaf and stem formation, with leaf development as the hallmark. The transition to the reproductive phase induces significant changes in these processes, driven mainly by the protein FT (FLOWERING LOCUS T) in Arabidopsis and proteins encoded by FT orthologs, which are specified as 'florigen'. These proteins are synthesized in leaves and transported to the SAM, and act as the primary flowering signal, although its impact varies among species. Within the SAM, florigen integrates with other signals, culminating in developmental changes. This review explores the central question of how florigen induces developmental phase transition in the SAM. Future research may combine phase transition studies, potentially revealing the florigen-induced developmental phase transition in the SAM.

Studies on Candidate Genes Related to Flowering Time in a Multiparent Population of Maize Derived from Tropical and Temperate Germplasm

A comprehensive study on maize flowering traits, focusing on the regulation of flowering time and the elucidation of molecular mechanisms underlying the genes controlling flowering, holds the potential to significantly enhance our understanding of the associated regulatory gene network. In this study, three tropical maize inbreds, CML384, CML171, and CML444, were used, along with a temperate maize variety, Shen137, as parental lines to cross with Ye107. The resulting F1s underwent seven consecutive generations of self-pollination through the single-seed descent (SSD) method to develop a multiparent population. To investigate the regulation of maize flowering time-related traits and to identify loci and candidate genes, a genome-wide association study (GWAS) was conducted. GWAS analysis identified 556 SNPs and 12 candidate genes that were significantly associated with flowering time-related traits. Additionally, an analysis of the effect of the estimated breeding values of the subpopulations on flowering time was conducted to further validate the findings of the present study. Collectively, this study offers valuable insights into novel candidate genes, contributing to an improved understanding of maize flowering time-related traits. This information holds practical significance for future maize breeding programs aimed at developing high-yielding hybrids.

Green revolution to genome revolution: driving better resilient crops against environmental instability

Crop improvement programmes began with traditional breeding practices since the inception of agriculture. Farmers and plant breeders continue to use these strategies for crop improvement due to their broad application in modifying crop genetic compositions. Nonetheless, conventional breeding has significant downsides in regard to effort and time. Crop productivity seems to be hitting a plateau as a consequence of environmental issues and the scarcity of agricultural land. Therefore, continuous pursuit of advancement in crop improvement is essential. Recent technical innovations have resulted in a revolutionary shift in the pattern of breeding methods, leaning further towards molecular approaches. Among the promising approaches, marker-assisted selection, QTL mapping, omics-assisted breeding, genome-wide association studies and genome editing have lately gained prominence. Several governments have progressively relaxed their restrictions relating to genome editing. The present review highlights the evolutionary and revolutionary approaches that have been utilized for crop improvement in a bid to produce climate-resilient crops observing the consequence of climate change. Additionally, it will contribute to the comprehension of plant breeding succession so far. Investing in advanced sequencing technologies and bioinformatics will deepen our understanding of genetic variations and their functional implications, contributing to breakthroughs in crop improvement and biodiversity conservation.

Genomic-regions associated with cold stress tolerance in Asia-adapted tropical maize germplasm

Maize is gaining impetus in non-traditional and non-conventional seasons such as off-season, primarily due to higher demand and economic returns. Maize varieties directed for growing in the winter season of South Asia must have cold resilience as an important trait due to the low prevailing temperatures and frequent cold snaps observed during this season in most parts of the lowland tropics of Asia. The current study involved screening of a panel of advanced tropically adapted maize lines to cold stress during vegetative and flowering stage under field conditions. A suite of significant genomic loci (28) associated with grain yield along and agronomic traits such as flowering (15) and plant height (6) under cold stress environments. The haplotype regression revealed 6 significant haplotype blocks for grain yield under cold stress across the test environments. Haplotype blocks particularly on chromosomes 5 (bin5.07), 6 (bin6.02), and 9 (9.03) co-located to regions/bins that have been identified to contain candidate genes involved in membrane transport system that would provide essential tolerance to the plant. The regions on chromosome 1 (bin1.04), 2 (bin 2.07), 3 (bin 3.05-3.06), 5 (bin5.03), 8 (bin8.05-8.06) also harboured significant SNPs for the other agronomic traits. In addition, the study also looked at the plausibility of identifying tropically adapted maize lines from the working germplasm with cold resilience across growth stages and identified four lines that could be used as breeding starts in the tropical maize breeding pipelines.

Elucidating the patterns of pleiotropy and its biological relevance in maize

Pleiotropy-when a single gene controls two or more seemingly unrelated traits-has been shown to impact genes with effects on flowering time, leaf architecture, and inflorescence morphology in maize. However, the genome-wide impact of biological pleiotropy across all maize phenotypes is largely unknown. Here, we investigate the extent to which biological pleiotropy impacts phenotypes within maize using GWAS summary statistics reanalyzed from previously published metabolite, field, and expression phenotypes across the Nested Association Mapping population and Goodman Association Panel. Through phenotypic saturation of 120,597 traits, we obtain over 480 million significant quantitative trait nucleotides. We estimate that only 1.56-32.3% of intervals show some degree of pleiotropy. We then assess the relationship between pleiotropy and various biological features such as gene expression, chromatin accessibility, sequence conservation, and enrichment for gene ontology terms. We find very little relationship between pleiotropy and these variables when compared to permuted pleiotropy. We hypothesize that biological pleiotropy of common alleles is not widespread in maize and is highly impacted by nuisance terms such as population structure and linkage disequilibrium. Natural selection on large standing natural variation in maize populations may target wide and large effect variants, leaving the prevalence of detectable pleiotropy relatively low.

Interspecies transfer of RAMOSA1 orthologs and promoter cis sequences impacts maize inflorescence architecture

Grass inflorescences support floral structures that each bear a single grain, where variation in branch architecture directly impacts yield. The maize (Zea mays) RAMOSA1 (ZmRA1) transcription factor acts as a key regulator of inflorescence development by imposing branch meristem determinacy. Here, we show RA1 transcripts accumulate in boundary domains adjacent to spikelet meristems in sorghum (Sorghum bicolor, Sb) and green millet (Setaria viridis, Sv) inflorescences similar as in the developing maize tassel and ear. To evaluate the functional conservation of syntenic RA1 orthologs and promoter cis sequences in maize, sorghum, and setaria, we utilized interspecies gene transfer and assayed genetic complementation in a common inbred background by quantifying recovery of normal branching in highly ramified ra1-R mutants. A ZmRA1 transgene that includes endogenous upstream and downstream flanking sequences recovered normal tassel and ear branching in ra1-R. Interspecies expression of two transgene variants of the SbRA1 locus, modeled as the entire endogenous tandem duplication or just the nonframeshifted downstream copy, complemented ra1-R branching defects and induced unusual fasciation and branch patterns. The SvRA1 locus lacks conserved, upstream noncoding cis sequences found in maize and sorghum; interspecies expression of a SvRA1 transgene did not or only partially recovered normal inflorescence forms. Driving expression of the SvRA1 coding region by the ZmRA1 upstream region, however, recovered normal inflorescence morphology in ra1-R. These data leveraging interspecies gene transfer suggest that cis-encoded temporal regulation of RA1 expression is a key factor in modulating branch meristem determinacy that ultimately impacts grass inflorescence architecture.

Molecular genetic variation of animals and plants under domestication

Domesticated plants and animals played crucial roles as models for evolutionary change by means of natural selection and for establishing the rules of inheritance, originally proposed by Charles Darwin and Gregor Mendel, respectively. Here, we review progress that has been made during the last 35 y in unraveling the molecular genetic variation underlying the stunning phenotypic diversity in crops and domesticated animals that inspired Mendel and Darwin. We notice that numerous domestication genes, crucial for the domestication process, have been identified in plants, whereas animal domestication appears to have a polygenic background with no obvious "domestication genes" involved. Although model organisms, such as Drosophila and Arabidopsis, have replaced domesticated species as models for basic research, the latter are still outstanding models for evolutionary research because phenotypic change in these species represents an evolutionary process over thousands of years. A consequence of this is that some alleles contributing to phenotypic diversity have evolved by accumulating multiple changes in the same gene. The continued molecular characterization of crops and farm animals with ever sharper tools is essential for future food security.

Deciphering Pleiotropic Signatures of Regulatory SNPs in Zea mays L. Using Multi-Omics Data and Machine Learning Algorithms

Maize is one of the most widely grown cereals in the world. However, to address the challenges in maize breeding arising from climatic anomalies, there is a need for developing novel strategies to harness the power of multi-omics technologies. In this regard, pleiotropy is an important genetic phenomenon that can be utilized to simultaneously enhance multiple agronomic phenotypes in maize. In addition to pleiotropy, another aspect is the consideration of the regulatory SNPs (rSNPs) that are likely to have causal effects in phenotypic development. By incorporating both aspects in our study, we performed a systematic analysis based on multi-omics data to reveal the novel pleiotropic signatures of rSNPs in a global maize population. For this purpose, we first applied Random Forests and then Markov clustering algorithms to decipher the pleiotropic signatures of rSNPs, based on which hierarchical network models are constructed to elucidate the complex interplay among transcription factors, rSNPs, and phenotypes. The results obtained in our study could help to understand the genetic programs orchestrating multiple phenotypes and thus could provide novel breeding targets for the simultaneous improvement of several agronomic traits.

Quantitative and population genomics suggest a broad role of stay-green loci in the drought adaptation of sorghum

Drought is a major constraint on plant productivity globally. Sorghum [Sorghum bicolor (L.) Moench] landraces have evolved in drought-prone regions, but the genetics of their adaptation is poorly understood. Here we sought to identify novel drought-tolerance loci and test hypotheses on the role of known loci including those underlying stay-green (Stg) postflowering drought tolerance. We phenotyped 590 diverse sorghum accessions from West Africa in 10 environments, under field-based managed drought stress [preflowering water stress (WS1), postflowering water stress (WS2), and well-watered (WW)] and rainfed (RF) conditions over 4 yr. Days to 50% flowering (DFLo), aboveground dry biomass (DBM), plant height (PH), and plant grain yield components (including grain weight [GrW], panicle weight [PW] and grain number [GrN] per plant, and 1000-grain weight [TGrW]) were measured, and genome-wide association studies (GWAS) was conducted. Broad-sense heritability for biomass and plant grain yield was high (33-92%) across environments. There was a significant correlation between stress tolerance index (STI) for GrW per plant across WS1 and WS2. Genome-wide association studies revealed that SbZfl1 and SbCN12, orthologs of maize (Zea mays L.) flowering genes, likely underlie flowering time variation under these conditions. Genome-wide association studies further identified associations (n = 134; common between two GWAS models) for STI and drought effects on plant yield components including 16 putative pleiotropic associations. Thirty of the associations colocalized with Stg1, Stg2, Stg3, and Stg4 loci and had large effects. Seven lead associations, including some within Stg1, overlapped with positive selection outliers. Our findings reveal previously undescribed natural genetic variation for drought-tolerance-related traits and suggest a broad role of Stg loci in drought adaptation of sorghum.

Identification of a major QTL and genome-wide epistatic interactions for single vs. paired spikelets in a maize-teosinte F2 population

Maize ear carries paired spikelets, whereas the ear of its wild ancestor, teosinte, bears single spikelets. However, little is known about the genetic basis of the processes of transformation of single spikelets in teosinte ear to paired spikelets in maize ear. In this study, a two-ranked, paired-spikelets primitive maize and a two-ranked, single-spikelet teosinte were utilized to develop an F2 population, and quantitative trait locus (loci) (QTL) mapping for single vs. paired spikelets (PEDS) was performed. One major QTL (qPEDS3.1) for PEDS located on chromosome 3S was identified in the 162 F2 plants using the inclusive composite interval mapping of additive (ICIM-ADD) module, explaining 23.79% of the phenotypic variance. Out of the 409 F2 plants, 43 plants with PEDS = 0% and 43 plants with PEDS > 20% were selected for selective genotyping, and the QTL (qPEDS3.1) was detected again. Moreover, the QTL (qPEDS3.1) was validated in three environments, which explained 31.05%, 38.94%, and 23.16% of the phenotypic variance, respectively. In addition, 50 epistatic QTLs were detected in the 162 F2 plants using the two-locus epistatic QTL (ICIM-EPI) module; they were distributed on all 10 chromosomes and explained 94.40% of the total phenotypic variance. The results contribute to a better understanding of the genetic basis of domestication of paired spikelets and provide a genetic resource for future map-based cloning; in addition, the systematic dissection of epistatic interactions underlies a theoretical framework for overcoming epistatic effects on QTL fine mapping.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01276-x.

The relative impact of evolving pleiotropy and mutational correlation on trait divergence

Both pleiotropic connectivity and mutational correlations can restrict the decoupling of traits under divergent selection, but it is unknown which is more important in trait evolution. To address this question, we create a model that permits within-population variation in both pleiotropic connectivity and mutational correlation, and compare their relative importance to trait evolution. Specifically, we developed an individual-based stochastic model where mutations can affect whether a locus affects a trait and the extent of mutational correlations in a population. We find that traits can decouple whether there is evolution in pleiotropic connectivity or mutational correlation, but when both can evolve, then evolution in pleiotropic connectivity is more likely to allow for decoupling to occur. The most common genotype found in this case is characterized by having one locus that maintains connectivity to all traits and another that loses connectivity to the traits under stabilizing selection (subfunctionalization). This genotype is favored because it allows the subfunctionalized locus to accumulate greater effect size alleles, contributing to increasingly divergent trait values in the traits under divergent selection without changing the trait values of the other traits (genetic modularization). These results provide evidence that partial subfunctionalization of pleiotropic loci may be a common mechanism of trait decoupling under regimes of corridor selection.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

Dissecting the Genetic Basis of Flowering Time and Height Related-Traits Using Two Doubled Haploid Populations in Maize

In the field, maize flowering time and height traits are closely linked with yield, planting density, lodging resistance, and grain fill. To explore the genetic basis of flowering time and height traits in maize, we investigated six related traits, namely, days to anthesis (AD), days to silking (SD), the anthesis-silking interval (ASI), plant height (PH), ear height (EH), and the EH/PH ratio (ER) in two locations for two years based on two doubled haploid (DH) populations. Based on the two high-density genetic linkage maps, 12 and 22 quantitative trait loci (QTL) were identified, respectively, for flowering time and height-related traits. Of these, ten QTLs had overlapping confidence intervals between the two populations and were integrated into three consensus QTLs (qFT_YZ1a, qHT_YZ5a, and qHT_YZ7a). Of these, qFT_YZ1a, conferring flowering time, is located at 221.1-277.0 Mb on chromosome 1 and explained 10.0-12.5% of the AD and SD variation, and qHT_YZ5a, conferring height traits, is located at 147.4-217.3 Mb on chromosome 5 and explained 11.6-15.3% of the PH and EH variation. These consensus QTLs, in addition to the other repeatedly detected QTLs, provide useful information for further genetic studies and variety improvements in flowering time and height-related traits.

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.

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.

Harnessing Knowledge from Maize and Rice Domestication for New Crop Breeding

Crop domestication has fundamentally altered the course of human history, causing a shift from hunter-gatherer to agricultural societies and stimulating the rise of modern civilization. A greater understanding of crop domestication would provide a theoretical basis for how we could improve current crops and develop new crops to deal with environmental challenges in a sustainable manner. Here, we provide a comprehensive summary of the similarities and differences in the domestication processes of maize and rice, two major staple food crops that feed the world. We propose that maize and rice might have evolved distinct genetic solutions toward domestication. Maize and rice domestication appears to be associated with distinct regulatory and evolutionary mechanisms. Rice domestication tended to select de novo, loss-of-function, coding variation, while maize domestication more frequently favored standing, gain-of-function, regulatory variation. At the gene network level, distinct genetic paths were used to acquire convergent phenotypes in maize and rice domestication, during which different central genes were utilized, orthologous genes played different evolutionary roles, and unique genes or regulatory modules were acquired for establishing new traits. Finally, we discuss how the knowledge gained from past domestication processes, together with emerging technologies, could be exploited to improve modern crop breeding and domesticate new crops to meet increasing human demands.

Teosinte (Zea mays ssp parviglumis) growth and transcriptomic response to weed stress identifies similarities and differences between varieties and with modern maize varieties

Transcriptomic responses of plants to weed presence gives insight on the physiological and molecular mechanisms involved in the stress response. This study evaluated transcriptomic and morphological responses of two teosinte (Zea mays ssp parviglumis) (an ancestor of domesticated maize) lines (Ames 21812 and Ames 21789) to weed presence and absence during two growing seasons. Responses were compared after 6 weeks of growth in Aurora, South Dakota, USA. Plant heights between treatments were similar in Ames 21812, whereas branch number decreased when weeds were present. Ames 21789 was 45% shorter in weedy vs weed-free plots, but branch numbers were similar between treatments. Season-long biomass was reduced in response to weed stress in both lines. Common down-regulated subnetworks in weed-stressed plants were related to light, photosynthesis, and carbon cycles. Several unique response networks (e.g. aging, response to chitin) and gene sets were present in each line. Comparing transcriptomic responses of maize (determined in an adjacent study) and teosinte lines indicated three common gene ontologies up-regulated when weed-stressed: jasmonic acid response/signaling, UDP-glucosyl and glucuronyltransferases, and quercetin glucosyltransferase (3-O and 7-O). Overall, morphologic and transcriptomic differences suggest a greater varietal (rather than a conserved) response to weed stress, and implies multiple responses are possible. These findings offer insights into opportunities to define and manipulate gene expression of several different pathways of modern maize varieties to improve performance under weedy conditions.

Evolution of multifunctionality through a pleiotropic substitution in the innate immune protein S100A9

Multifunctional proteins are evolutionary puzzles: how do proteins evolve to satisfy multiple functional constraints? S100A9 is one such multifunctional protein. It potently amplifies inflammation via Toll-like receptor four and is antimicrobial as part of a heterocomplex with S100A8. These two functions are seemingly regulated by proteolysis: S100A9 is readily degraded, while S100A8/S100A9 is resistant. We take an evolutionary biochemical approach to show that S100A9 evolved both functions and lost proteolytic resistance from a weakly proinflammatory, proteolytically resistant amniote ancestor. We identify a historical substitution that has pleiotropic effects on S100A9 proinflammatory activity and proteolytic resistance but has little effect on S100A8/S100A9 antimicrobial activity. We thus propose that mammals evolved S100A8/S100A9 antimicrobial and S100A9 proinflammatory activities concomitantly with a proteolytic 'timer' to selectively regulate S100A9. This highlights how the same mutation can have pleiotropic effects on one functional state of a protein but not another, thus facilitating the evolution of multifunctionality.

Plant Inflorescence Architecture: The Formation, Activity, and Fate of Axillary Meristems

The above-ground plant body in different plant species can have very distinct forms or architectures that arise by recurrent redeployment of a finite set of building blocks-leaves with axillary meristems, stems or branches, and flowers. The unique architectures of plant inflorescences in different plant families and species, on which this review focuses, determine the reproductive success and yield of wild and cultivated plants. Major contributors to the inflorescence architecture are the activity and developmental trajectories adopted by axillary meristems, which determine the degree of branching and the number of flowers formed. Recent advances in genetic and molecular analyses in diverse flowering plants have uncovered both common regulatory principles and unique players and/or regulatory interactions that underlie inflorescence architecture. Modulating activity of these regulators has already led to yield increases in the field. Additional insight into the underlying regulatory interactions and principles will not only uncover how their rewiring resulted in altered plant form, but will also enhance efforts at optimizing plant architecture in desirable ways in crop species.

Large-scale metabolite quantitative trait locus analysis provides new insights for high-quality maize improvement

It is generally recognized that many favorable genes which were lost during domestication, including those related to both nutritional value and stress resistance, remain hidden in wild relatives. To uncover such genes in teosinte, an ancestor of maize, we conducted metabolite profiling in a BC₂ F₇ population generated from a cross between the maize wild relative (Zea mays ssp. mexicana) and maize inbred line Mo17. In total, 65 primary metabolites were quantified in four tissues (seedling-stage leaf, grouting-stage leaf, young kernel and mature kernel) with clear tissue-specific patterns emerging. Three hundred and fifty quantitative trait loci (QTLs) for these metabolites were obtained, which were distributed unevenly across the genome and included two QTL hotspots. Metabolite concentrations frequently increased in the presence of alleles from the teosinte genome while the opposite was observed for grain yield and shape trait QTLs. Combination of the multi-tissue transcriptome and metabolome data provided considerable insight into the metabolic variations between maize and its wild relatives. This study thus identifies favorable genes hidden in the wild relative which should allow us to balance high yield and quality in future modern crop breeding programs.

Over-expression of the photoperiod response regulator ZmCCT10 modifies plant architecture, flowering time and inflorescence morphology in maize

Maize originated as a tropical plant that required short days to transition from vegetative to reproductive development. ZmCCT10 [CO, CONSTANS, CO-LIKE and TIMING OF CAB1 (CCT) transcription factor family] is a regulator of photoperiod response and was identified as a major QTL controlling photoperiod sensitivity in maize. We modulated expression of ZmCCT10 in transgenic maize using two constitutive promoters with different expression levels. Transgenic plants over expressing ZmCCT10 with either promoter were delayed in their transition from vegetative to reproductive development but were not affected in their switch from juvenile-to-adult vegetative growth. Strikingly, transgenic plants containing the stronger expressing construct had a prolonged period of vegetative growth accompanied with dramatic modifications to plant architecture that impacted both vegetative and reproductive traits. These plants did not produce ears, but tassels were heavily branched. In more than half of the transgenic plants, tassels were converted into a branched leafy structure resembling phyllody, often composed of vegetative plantlets. Analysis of expression modules controlling the floral transition and meristem identity linked these networks to photoperiod dependent regulation, whereas phase change modules appeared to be photoperiod independent. Results from this study clarified the influence of the photoperiod pathway on vegetative and reproductive development and allowed for the fine-tuning of the maize flowering time model.

Genome-wide association analysis of stalk biomass and anatomical traits in maize

BACKGROUND

Maize stover is an important source of crop residues and a promising sustainable energy source in the United States. Stalk is the main component of stover, representing about half of stover dry weight. Characterization of genetic determinants of stalk traits provide a foundation to optimize maize stover as a biofuel feedstock. We investigated maize natural genetic variation in genome-wide association studies (GWAS) to detect candidate genes associated with traits related to stalk biomass (stalk diameter and plant height) and stalk anatomy (rind thickness, vascular bundle density and area).

RESULTS

Using a panel of 942 diverse inbred lines, 899,784 RNA-Seq derived single nucleotide polymorphism (SNP) markers were identified. Stalk traits were measured on 800 members of the panel in replicated field trials across years. GWAS revealed 16 candidate genes associated with four stalk traits. Most of the detected candidate genes were involved in fundamental cellular functions, such as regulation of gene expression and cell cycle progression. Two of the regulatory genes (Zmm22 and an ortholog of Fpa) that were associated with plant height were previously shown to be involved in regulating the vegetative to floral transition. The association of Zmm22 with plant height was confirmed using a transgenic approach. Transgenic lines with increased expression of Zmm22 showed a significant decrease in plant height as well as tassel branch number, indicating a pleiotropic effect of Zmm22.

CONCLUSION

Substantial heritable variation was observed in the association panel for stalk traits, indicating a large potential for improving useful stalk traits in breeding programs. Genome-wide association analyses detected several candidate genes associated with multiple traits, suggesting common regulatory elements underlie various stalk traits. Results of this study provide insights into the genetic control of maize stalk anatomy and biomass.

Selection Signatures Underlying Dramatic Male Inflorescence Transformation During Modern Hybrid Maize Breeding

Inflorescence capacity plays a crucial role in reproductive fitness in plants, and in production of hybrid crops. Maize is a monoecious species bearing separate male and female flowers (tassel and ear, respectively). The switch from open-pollinated populations of maize to hybrid-based breeding schemes in the early 20th century was accompanied by a dramatic reduction in tassel size, and the trend has continued with modern breeding over the recent decades. The goal of this study was to identify selection signatures in genes that may underlie this dramatic transformation. Using a population of 942 diverse inbred maize accessions and a nested association mapping population comprising three 200-line biparental populations, we measured 15 tassel morphological characteristics by manual and image-based methods. Genome-wide association studies identified 242 single nucleotide polymorphisms significantly associated with measured traits. We compared 41 unselected lines from the Iowa Stiff Stalk Synthetic (BSSS) population to 21 highly selected lines developed by modern commercial breeding programs, and found that tassel size and weight were reduced significantly. We assayed genetic differences between the two groups using three selection statistics: cross population extended haplotype homozogysity, cross-population composite likelihood ratio, and fixation index. All three statistics show evidence of selection at genomic regions associated with tassel morphology relative to genome-wide null distributions. These results support the tremendous effect, both phenotypic and genotypic, that selection has had on maize male inflorescence morphology.

Maize domestication and gene interaction

Contents Summary 395 I. Introduction 395 II. The genetic basis of maize domestication 396 III. The tempo of maize domestication 401 IV. Genetic interactions and selection during maize domestication 401 V. Gene networks of maize domestication alleles 404 VI. Implications of gene interactions on evolution and selection404 VII. Conclusions 405 Acknowledgements 405 References 405 SUMMARY: Domestication is a tractable system for following evolutionary change. Under domestication, wild populations respond to shifting selective pressures, resulting in adaptation to the new ecological niche of cultivation. Owing to the important role of domesticated crops in human nutrition and agriculture, the ancestry and selection pressures transforming a wild plant into a domesticate have been extensively studied. In Zea mays, morphological, genetic and genomic studies have elucidated how a wild plant, the teosinte Z. mays subsp. parviglumis, was transformed into the domesticate Z. mays subsp. mays. Five major morphological differences distinguish these two subspecies, and careful genetic dissection has pinpointed the molecular changes responsible for several of these traits. But maize domestication was a consequence of more than just five genes, and regions throughout the genome contribute. The impacts of these additional regions are contingent on genetic background, both the interactions between alleles of a single gene and among alleles of the multiple genes that modulate phenotypes. Key genetic interactions include dominance relationships, epistatic interactions and pleiotropic constraint, including how these variants are connected in gene networks. Here, we review the role of gene interactions in generating the dramatic phenotypic evolution seen in the transition from teosinte to maize.

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.

Exploiting SPL genes to improve maize plant architecture tailored for high-density planting

Maize (Zea mays ssp. mays) is an agronomically important crop and also a classical genetic model for studying the regulation of plant architecture formation, which is a critical determinant of grain yield. Since the 1930s, increasing planting density has been a major contributing factor to the >7-fold increase in maize grain yield per unit land area in the USA, which is accompanied by breeding and utilization of cultivars characterized by high-density-tolerant plant architecture, including decreased ear height, lodging resistance, more upright leaves, reduced tassel branch number, and reduced anthesis-silking interval (ASI). Recent studies demonstrated that phytochrome-mediated red/far-red light signaling pathway and the miR156/SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) regulatory module co-ordinately regulate the shade avoidance response and diverse aspects of plant architecture in responding to shading in Arabidopsis. The maize genome contains 30 ZmSPL genes, and 18 of them are predicted as direct targets of zma-miR156s. Accumulating evidence indicates that ZmSPL genes play important roles in regulating maize flowering time, plant/ear height, tilling, leaf angle, tassel and ear architecture, and grain size and shape. Finally, we discuss ways to exploit maize SPL genes and downstream targets for improving maize plant architecture tailored for high-density planting.

Asymmetric expression patterns reveal a strong maternal effect and dosage compensation in polyploid hybrid fish

Background

Hybridization and polyploidization are regarded as the major driving forces in plant speciation, diversification, and ecological adaptation. Our knowledge regarding the mechanisms of duplicated-gene regulation following genomic merging or doubling is primarily derived from plants and is sparse for vertebrates.

Results

We successfully obtained an F1 generation (including allodiploid hybrids and triploid hybrids) from female Megalobrama amblycephala Yih (BSB, 2n = 48) × male Xenocypri davidi Bleeker (YB, 2n = 48). The duplicated-gene expression patterns of the two types of hybrids were explored using RNA-Seq data. In total, 5.44 × 10⁸ (69.32 GB) clean reads and 499,631 assembled unigenes were obtained from the testis transcriptomes. The sequence similarity analysis of 4265 orthologs revealed that the merged genomes were dominantly expressed in different ploidy hybrids. The differentially expressed genes in the two types of hybrids were asymmetric compared with those in both parents. Furthermore, the genome-wide expression level dominance (ELD) was biased toward the maternal BSB genome in both the allodiploid and triploid hybrids. In addition, the dosage-compensation mechanisms that reduced the triploid expression levels to the diploid state were determined in the triploid hybrids.

Conclusions

Our results indicate that divergent genomes undergo strong interactions and domination in allopolyploid offspring. Genomic merger has a greater effect on the gene-expression patterns than genomic doubling. The various expression mechanisms (including maternal effect and dosage compensation) in different ploidy hybrids suggest that the initial genomic merger and doubling play important roles in polyploidy adaptation and evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4883-7) contains supplementary material, which is available to authorized users.

Identification of minor effect QTLs for plant architecture related traits using super high density genotyping and large recombinant inbred population in maize (Zea mays)

BACKGROUND

Plant Architecture Related Traits (PATs) are of great importance for maize breeding, and mainly controlled by minor effect quantitative trait loci (QTLs). However, cloning or even fine-mapping of minor effect QTLs is very difficult in maize. Theoretically, large population and high density genetic map can be helpful for increasing QTL mapping resolution and accuracy, but such a possibility have not been actually tested.

RESULTS

Here, we employed a genotyping-by-sequencing (GBS) strategy to construct a linkage map with 16,769 marker bins for 1021 recombinant inbred lines (RILs). Accurately mapping of well studied genes P1, pl1 and r1 underlying silk color demonstrated the map quality. After QTL analysis, a total of 51 loci were mapped for six PATs. Although all of them belong to minor effect alleles, the lengths of the QTL intervals, with a minimum and median of 1.03 and 3.40 Mb respectively, were remarkably reduced as compared with previous reports using smaller size of population or small number of markers. Several genes with known function in maize were shown to be overlapping with or close neighboring to these QTL peaks, including na1, td1, d3 for plant height, ra1 for tassel branch number, and zfl2 for tassel length. To further confirm our mapping results, a plant height QTL, qPH1a, was verified by an introgression lines (ILs).

CONCLUSIONS

We demonstrated a method for high resolution mapping of minor effect QTLs in maize, and the resulted comprehensive QTLs for PATs are valuable for maize molecular breeding in the future.

The potential role of genetic assimilation during maize domestication

Domestication research has largely focused on identification of morphological and genetic differences between extant populations of crops and their wild relatives. Little attention has been paid to the potential effects of environment despite substantial known changes in climate from the time of domestication to modern day. In recent research, the exposure of teosinte (i.e., wild maize) to environments similar to the time of domestication, resulted in a plastic induction of domesticated phenotypes in teosinte. These results suggest that early agriculturalists may have selected for genetic mechanisms that cemented domestication phenotypes initially induced by a plastic response of teosinte to environment, a process known as genetic assimilation. To better understand this phenomenon and the potential role of environment in maize domestication, we examined differential gene expression in maize (Zea mays ssp. mays) and teosinte (Zea mays ssp. parviglumis) between past and present conditions. We identified a gene set of over 2000 loci showing a change in expression across environmental conditions in teosinte and invariance in maize. In fact, overall we observed both greater plasticity in gene expression and more substantial changes in co-expressionnal networks in teosinte across environments when compared to maize. While these results suggest genetic assimilation played at least some role in domestication, genes showing expression patterns consistent with assimilation are not significantly enriched for previously identified domestication candidates, indicating assimilation did not have a genome-wide effect.

Exploring and Harnessing Haplotype Diversity to Improve Yield Stability in Crops

In order to meet future food, feed, fiber, and bioenergy demands, global yields of all major crops need to be increased significantly. At the same time, the increasing frequency of extreme weather events such as heat and drought necessitates improvements in the environmental resilience of modern crop cultivars. Achieving sustainably increase yields implies rapid improvement of quantitative traits with a very complex genetic architecture and strong environmental interaction. Latest advances in genome analysis technologies today provide molecular information at an ultrahigh resolution, revolutionizing crop genomic research, and paving the way for advanced quantitative genetic approaches. These include highly detailed assessment of population structure and genotypic diversity, facilitating the identification of selective sweeps and signatures of directional selection, dissection of genetic variants that underlie important agronomic traits, and genomic selection (GS) strategies that not only consider major-effect genes. Single-nucleotide polymorphism (SNP) markers today represent the genotyping system of choice for crop genetic studies because they occur abundantly in plant genomes and are easy to detect. SNPs are typically biallelic, however, hence their information content compared to multiallelic markers is low, limiting the resolution at which SNP-trait relationships can be delineated. An efficient way to overcome this limitation is to construct haplotypes based on linkage disequilibrium, one of the most important features influencing genetic analyses of crop genomes. Here, we give an overview of the latest advances in genomics-based haplotype analyses in crops, highlighting their importance in the context of polyploidy and genome evolution, linkage drag, and co-selection. We provide examples of how haplotype analyses can complement well-established quantitative genetics frameworks, such as quantitative trait analysis and GS, ultimately providing an effective tool to equip modern crops with environment-tailored characteristics.

Analysis of the genetic architecture of maize ear and grain morphological traits by combined linkage and association mapping

Using combined linkage and association mapping, 26 stable QTL and six stable SNPs were detected across multiple environments for eight ear and grain morphological traits in maize. One QTL, PKS2, might play an important role in maize yield improvement. In the present study, one bi-parental population and an association panel were used to identify quantitative trait loci (QTL) for eight ear and grain morphological traits. A total of 108 QTL related to these traits were detected across four environments using an ultra-high density bin map constructed using recombinant inbred lines (RILs) derived from a cross between Ye478 and Qi319, and 26 QTL were identified in more than two environments. Furthermore, 64 single nucleotide polymorphisms (SNPs) were found to be significantly associated with the eight ear and grain morphological traits (-log₁₀(P) > 4) in an association panel of 240 maize inbred lines. Combining the two mapping populations, a total of 17 pleiotropic QTL/SNPs (pQTL/SNPs) were associated with various traits across multiple environments. PKS2, a stable locus influencing kernel shape identified on chromosome 2 in a genome-wide association study (GWAS), was within the QTL confidence interval defined by the RILs. The candidate region harbored a short 13-Kb LD block encompassing four SNPs (SYN11386, PHM14783.16, SYN11392, and SYN11378). In the association panel, 13 lines derived from the hybrid PI78599 possessed the same allele as Qi319 at the PHM14783.16 (GG) locus, with an average value of 0.21 for KS, significantly lower than that of the 34 lines derived from Ye478 that carried a different allele (0.25, P < 0.05). Therefore, further fine mapping of PKS2 will provide valuable information for understanding the genetic components of grain yield and improving molecular marker-assisted selection (MAS) in maize.

Surprisingly Low Limits of Selection in Plant Domestication

Current debate concerns the pace at which domesticated plants emerged from cultivated wild populations and how many genes were involved. Using an individual-based model, based on the assumptions of Haldane and Maynard Smith, respectively, we estimate that a surprisingly low number of 50–100 loci are the most that could be under selection in a cultivation regime at the selection strengths observed in the archaeological record. This finding is robust to attempts to rescue populations from extinction through selection from high standing genetic variation, gene flow, and the Maynard Smith-based model of threshold selection. Selective sweeps come at a cost, reducing the capacity of plants to adapt to new environments, which may contribute to the explanation of why selective sweeps have not been detected more frequently and why expansion of the agrarian package during the Neolithic was so frequently associated with collapse.

Expanding Maize Genetic Resources with Predomestication Alleles: Maize-Teosinte Introgression Populations

Teosinte ( subsp. H. H. Iltis & Doebley) has greater genetic diversity than maize inbreds and landraces ( subsp. ). There are, however, limited genetic resources to efficiently evaluate and tap this diversity. To broaden resources for genetic diversity studies in maize, we developed and evaluated 928 near-isogenic introgression lines (NILs) from 10 teosinte accessions in the B73 background. Joint linkage analysis of the 10 introgression populations identified several large-effect quantitative trait loci (QTL) for days to anthesis (DTA), kernel row number (KRN), and 50-kernel weight (Wt50k). Our results confirm prior reports of kernel domestication loci and identify previously uncharacterized QTL with a range of allelic effects enabling future research into the genetic basis of these traits. Additionally, we used a targeted set of NILs to validate the effects of a KRN QTL located on chromosome 2. These introgression populations offer novel tools for QTL discovery and validation as well as a platform for initiating fine mapping.

SUPPRESSOR OF APICAL DOMINANCE1 of Sporisorium reilianum Modulates Inflorescence Branching Architecture in Maize and Arabidopsis

The biotrophic fungus Sporisorium reilianum causes head smut of maize (Zea mays) after systemic plant colonization. Symptoms include the formation of multiple female inflorescences at subapical nodes of the stalk because of loss of apical dominance. By deletion analysis of cluster 19-1, the largest genomic divergence cluster in S. reilianum, we identified a secreted fungal effector responsible for S. reilianum-induced loss of apical dominance, which we named SUPPRESSOR OF APICAL DOMINANCE1 (SAD1). SAD1 transcript levels were highly up-regulated during biotrophic fungal growth in all infected plant tissues. SAD1-green fluorescent protein fusion proteins expressed by recombinant S. reilianum localized to the extracellular hyphal space. Transgenic Arabidopsis (Arabidopsis thaliana)-expressing green fluorescent protein-SAD1 displayed an increased number of secondary rosette-leaf branches. This suggests that SAD1 manipulates inflorescence branching architecture in maize and Arabidopsis through a conserved pathway. Using a yeast (Saccharomyces cerevisiae) two-hybrid library of S. reilianum-infected maize tissues, we identified potential plant interaction partners that had a predicted function in ubiquitination, signaling, and nuclear processes. Presence of SAD1 led to an increase of the transcript levels of the auxin transporter PIN-FORMED1 in the root and a reduction of the branching regulator TEOSINTE BRANCHED1 in the stalk. This indicates a role of SAD1 in regulation of apical dominance by modulation of branching through increasing transcript levels of the auxin transporter PIN1 and derepression of bud outgrowth.

KRN4 Controls Quantitative Variation in Maize Kernel Row Number

Kernel row number (KRN) is an important component of yield during the domestication and improvement of maize and controlled by quantitative trait loci (QTL). Here, we fine-mapped a major KRN QTL, KRN4, which can enhance grain productivity by increasing KRN per ear. We found that a ~3-Kb intergenic region about 60 Kb downstream from the SBP-box gene Unbranched3 (UB3) was responsible for quantitative variation in KRN by regulating the level of UB3 expression. Within the 3-Kb region, the 1.2-Kb Presence-Absence variant was found to be strongly associated with quantitative variation in KRN in diverse maize inbred lines, and our results suggest that this 1.2-Kb transposon-containing insertion is likely responsible for increased KRN. A previously identified A/G SNP (S35, also known as Ser220Asn) in UB3 was also found to be significantly associated with KRN in our association-mapping panel. Although no visible genetic effect of S35 alone could be detected in our linkage mapping population, it was found to genetically interact with the 1.2-Kb PAV to modulate KRN. The KRN4 was under strong selection during maize domestication and the favorable allele for the 1.2-Kb PAV and S35 has been significantly enriched in modern maize improvement process. The favorable haplotype (Hap1) of 1.2-Kb-PAV-S35 was selected during temperate maize improvement, but is still rare in tropical and subtropical maize germplasm. The dissection of the KRN4 locus improves our understanding of the genetic basis of quantitative variation in complex traits in maize.

Maize (Zea mays L.) is one of the world's most important sources of calories for humans. With an expanding global population, the demands for maize-derived food, feed, and fuel are rapidly increasing. To meet these needs, geneticists and breeders are facing the challenge of enhancing grain yield through genetic improvement of maize germplasm. Understanding the genetic basis of grain yield is necessary to guide breeding efforts towards the development of high-yielding hybrids. Kernel row number (KRN) in maize is one of the most important yield components and a significant breeding target. Over the last few decades, many genes that determine inflorescence development and architecture have been identified and characterized. The formation of kernel rows is an integral part of the development of the female inflorescence in maize. Nevertheless, the genetic basis and molecular regulation of quantitative variation in KRN is poorly understood. This study provides experimental evidence for the hypothesis that variation in intergenic regions can regulate quantitative variation of important grain yield-related traits, and also provides tools for improving KRN in maize.

Comparative transcriptomics uncovers alternative splicing changes and signatures of selection from maize improvement

Background

Alternative splicing (AS) is an important regulatory mechanism that greatly contributes to eukaryotic transcriptome diversity. A substantial amount of evidence has demonstrated that AS complexity is relevant to eukaryotic evolution, development, adaptation, and complexity. In this study, six teosinte and ten maize transcriptomes were sequenced to analyze AS changes and signatures of selection in maize domestication and improvement.

Results

In maize and teosinte, 13,593 highly conserved genes, including 12,030 multiexonic genes, were detected. By identifying AS isoforms from mutliexonic genes, we found that AS types were not significantly different between maize and teosinte. In addition, the two main AS types (intron retention and alternative acceptor) contributed to more than 60% of the AS events in the two species, but the average unique AS events per each alternatively spliced gene in maize (4.12) was higher than that in teosinte (2.26). Moreover, 94 genes generating 98 retained introns with transposable element (TE) sequences were detected in maize, which is far more than 9 retained introns with TEs detected in teosinte. This indicates that TE insertion might be an important mechanism for intron retention in maize. Additionally, the AS levels of 3864 genes were significantly different between maize and teosinte. Of these, 151 AS level-altered genes that are involved in transcriptional regulation and in stress responses are located in regions that have been targets of selection during maize improvement. These genes were inferred to be putatively improved genes.

Conclusions

We suggest that both maize and teosinte share similar AS mechanisms, but more genes have increased AS complexity during domestication from teosinte to maize. Importantly, a subset of AS level-increased genes that encode transcription factors and stress-responsive proteins may have been selected during maize improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1582-5) contains supplementary material, which is available to authorized users.

Evolution. Comment on "A promiscuous intermediate underlies the evolution of LEAFY DNA binding specificity"

Sayou et al. (Reports, 7 February 2014, p. 645) proposed a new model for evolution of transcription factors without gene duplication, using LEAFY as an archetype. Their proposal contradicts the evolutionary history of plants and ignores evidence that LEAFY evolves through gene duplications. Within their data set, we identified a moss with multiple LEAFY orthologs, which contests their model and supports that LEAFY evolves through duplications.

New insights into the regulation of inflorescence architecture

The architecture of inflorescences displays the spatiotemporal arrangement of flowers and determines plant reproductive success through affecting fruit set and plant interaction with biotic or abiotic factors. Flowering plants have evolved a remarkable diversity of inflorescence branching patterns, which is largely governed by developmental decisions in inflorescence meristems and their derived meristems between maintenance of indeterminacy and commitment to the floral fate. Recent findings suggest that regulation of inflorescence architecture is mediated by flowering time genes, Arabidopsis LSH1 and Oryza G1 (ALOG) family genes, and the interaction between the auxin pathway and floral meristem regulators. In this review, we discuss how the relevant new players and mechanisms account for the development of appropriate inflorescence structures in flowering plants in response to environmental and developmental signals.

Analyses of sequence polymorphism and haplotype diversity of LEAFY genes revealed post-domestication selection in the Chinese elite maize inbred lines

Post-domestication selection refers to the artificial selection on the loci controlling important agronomic traits during the process of genetic improvement in a population. The maize genes Zfl1 and Zfl2, duplicate orthologs of Arabidopsis LEAFY, are key regulators in plant branching, inflorescence and flower development, and reproduction. In this study, the full gene sequences of Zfl1 and Zfl2 from 62 Chinese elite inbred lines were amplified to evaluate their nucleotide polymorphisms and haplotype diversities. A total of 254 and 192 variants that included SNPs and indels were identified from the full sequences of Zfl1 and Zfl2, respectively. Although most of the variants were found to be located in the non-coding regions, the polymorphisms of CDS sequences classified Zfl1 into 16 haplotypes encoding 16 different proteins and Zfl2 into 18 haplotypes encoding eight different proteins. The population of Huangzaosi and its derived lines showed statistically significant signals of post-domestication selection on the Zfl1 CDS sequences, as well as lower nucleotide polymorphism and haplotype diversity than the whole set. However, the Zfl2 locus was only selected for in the heterotic group Reid. Further evidence revealed that at least 17 recombination events contributed to the genetic and haplotype diversities at the Zfl1 locus and 16 recombination events at the Zfl2 locus.

CACTA-like transposable element in ZmCCT attenuated photoperiod sensitivity and accelerated the postdomestication spread of maize

The postdomestication adaptation of maize to longer days required reduced photoperiod sensitivity to optimize flowering time. We performed a genome-wide association study and confirmed that ZmCCT, encoding a CCT domain-containing protein, is associated with the photoperiod response. In early-flowering maize we detected a CACTA-like transposable element (TE) within the ZmCCT promoter that dramatically reduced flowering time. TE insertion likely occurred after domestication and was selected as maize adapted to temperate zones. This process resulted in a strong selective sweep within the TE-related block of linkage disequilibrium. Functional validations indicated that the TE represses ZmCCT expression to reduce photoperiod sensitivity, thus accelerating maize spread to long-day environments.

Designed transcriptional regulators for trait development

Development is largely controlled by proteins that regulate gene expression at the level of transcription. These regulatory proteins, the genes that control them, and the genes that they control, are organized in a hierarchical structure of complex interactions. Altering the expression of genes encoding regulatory proteins controlling critical nodes in this hierarchy has potential for dramatic phenotypic modification. Constitutive over-expression of genes encoding regulatory proteins in transgenic plants has resulted in agronomically interesting phenotypes along with developmental abnormalities. For trait development, the magnitude and timing of expression of genes encoding key regulatory proteins will need to be precisely controlled and targeted to specific cells and tissues at certain developmental timepoints. Such control is made possible by designed transcriptional regulators which are fusions of engineered DNA binding proteins and activator or repressor domains. Expression of genes encoding such designed transcriptional regulators enable the selective modulation of endogenous gene expression. Genes encoding proteins controlling regulatory networks are prime targets for up- or down-regulation via such designed transcriptional regulators.

Anatomical and transcriptional dynamics of maize embryonic leaves during seed germination

Our anatomical analysis revealed that a dry maize seed contains four to five embryonic leaves at different developmental stages. Rudimentary kranz structure (KS) is apparent in the first leaf with a substantial density, but its density decreases toward younger leaves. Upon imbibition, leaf expansion occurs rapidly with new KSs initiated from the palisade-like ground meristem cells in the middle of the leaf. In parallel to the anatomical analysis, we obtained the time course transcriptomes for the embryonic leaves in dry and imbibed seeds every 6 h up to hour 72. Over this time course, the embryonic leaves exhibit transcripts of 30,255 genes at a level that can be regarded as "expressed." In dry seeds, ∼25,500 genes are expressed, showing functional enrichment in transcription, RNA processing, protein synthesis, primary metabolic pathways, and calcium transport. During the 72-h time course, ∼13,900 genes, including 590 transcription factor genes, are differentially expressed. Indeed, by 30 h postimbibition, ∼2,200 genes expressed in dry seeds are already down-regulated, and ∼2,000 are up-regulated. Moreover, the top 1% expressed genes at 54 h or later are very different from those before 30 h, reflecting important developmental and physiological transitions. Interestingly, clusters of genes involved in hormone metabolism, signaling, and responses are differentially expressed at various time points and TF gene expression is also modular and stage specific. Our dataset provides an opportunity for hypothesizing the timing of regulatory actions, particularly in the context of KS development.

Quantitative variation in maize kernel row number is controlled by the FASCIATED EAR2 locus

Domestication of cereal crops, such as maize, wheat and rice, had a profound influence on agriculture and the establishment of human civilizations. One major improvement was an increase in seed number per inflorescence, which enhanced yield and simplified harvesting and storage. The ancestor of maize, teosinte, makes 2 rows of kernels, and modern varieties make ∼8-20 rows. Kernel rows are initiated by the inflorescence shoot meristem, and shoot meristem size is controlled by a feedback loop involving the CLAVATA signaling proteins and the WUSCHEL transcription factor. We present a hypothesis that variation in inflorescence meristem size affects kernel row number (KRN), with the potential to increase yield. We also show that variation in the CLAVATA receptor-like protein FASCIATED EAR2 leads to increased inflorescence meristem size and KRN. These findings indicate that modulation of fundamental stem cell proliferation control pathways has the potential to enhance crop yields.

The plant ontology as a tool for comparative plant anatomy and genomic analyses

The Plant Ontology (PO; http://www.plantontology.org/) is a publicly available, collaborative effort to develop and maintain a controlled, structured vocabulary ('ontology') of terms to describe plant anatomy, morphology and the stages of plant development. The goals of the PO are to link (annotate) gene expression and phenotype data to plant structures and stages of plant development, using the data model adopted by the Gene Ontology. From its original design covering only rice, maize and Arabidopsis, the scope of the PO has been expanded to include all green plants. The PO was the first multispecies anatomy ontology developed for the annotation of genes and phenotypes. Also, to our knowledge, it was one of the first biological ontologies that provides translations (via synonyms) in non-English languages such as Japanese and Spanish. As of Release #18 (July 2012), there are about 2.2 million annotations linking PO terms to >110,000 unique data objects representing genes or gene models, proteins, RNAs, germplasm and quantitative trait loci (QTLs) from 22 plant species. In this paper, we focus on the plant anatomical entity branch of the PO, describing the organizing principles, resources available to users and examples of how the PO is integrated into other plant genomics databases and web portals. We also provide two examples of comparative analyses, demonstrating how the ontology structure and PO-annotated data can be used to discover the patterns of expression of the LEAFY (LFY) and terpene synthase (TPS) gene homologs.

A gene regulatory network model for floral transition of the shoot apex in maize and its dynamic modeling

The transition from the vegetative to reproductive development is a critical event in the plant life cycle. The accurate prediction of flowering time in elite germplasm is important for decisions in maize breeding programs and best agronomic practices. The understanding of the genetic control of flowering time in maize has significantly advanced in the past decade. Through comparative genomics, mutant analysis, genetic analysis and QTL cloning, and transgenic approaches, more than 30 flowering time candidate genes in maize have been revealed and the relationships among these genes have been partially uncovered. Based on the knowledge of the flowering time candidate genes, a conceptual gene regulatory network model for the genetic control of flowering time in maize is proposed. To demonstrate the potential of the proposed gene regulatory network model, a first attempt was made to develop a dynamic gene network model to predict flowering time of maize genotypes varying for specific genes. The dynamic gene network model is composed of four genes and was built on the basis of gene expression dynamics of the two late flowering id1 and dlf1 mutants, the early flowering landrace Gaspe Flint and the temperate inbred B73. The model was evaluated against the phenotypic data of the id1 dlf1 double mutant and the ZMM4 overexpressed transgenic lines. The model provides a working example that leverages knowledge from model organisms for the utilization of maize genomic information to predict a whole plant trait phenotype, flowering time, of maize genotypes.

CORNET 2.0: integrating plant coexpression, protein-protein interactions, regulatory interactions, gene associations and functional annotations

To enable easy access and interpretation of heterogeneous and scattered data, we have developed a user-friendly tool for data mining and integration in Arabidopsis, named CORNET. This tool allows the browsing of microarray data, the construction of coexpression and protein-protein interaction (PPI) networks and the exploration of diverse functional annotations. Here, we present the new functionalities of CORNET 2.0 for data integration in plants. First of all, CORNET allows the integration of regulatory interaction datasets accessible through the new transcription factor (TF) tool that can be used in combination with the coexpression tool or the PPI tool. In addition, we have extended the PPI tool to enable the analysis of gene-gene associations from AraNet as well as newly identified PPIs. Different search options are implemented to enable the construction of networks centered around multiple input genes or proteins. New functional annotation resources are included to retrieve relevant literature, phenotypes, plant ontology and biological pathways. We have also extended CORNET to attain the construction of coexpression and PPI networks in the crop species maize. Networks and associated evidence of the majority of currently available data types are visualized in Cytoscape. CORNET is available at https://bioinformatics.psb.ugent.be/cornet.

Reshaping of the maize transcriptome by domestication

Through domestication, humans have substantially altered the morphology of Zea mays ssp. parviglumis (teosinte) into the currently recognizable maize. This system serves as a model for studying adaptation, genome evolution, and the genetics and evolution of complex traits. To examine how domestication has reshaped the transcriptome of maize seedlings, we used expression profiling of 18,242 genes for 38 diverse maize genotypes and 24 teosinte genotypes. We detected evidence for more than 600 genes having significantly different expression levels in maize compared with teosinte. Moreover, more than 1,100 genes showed significantly altered coexpression profiles, reflective of substantial rewiring of the transcriptome since domestication. The genes with altered expression show a significant enrichment for genes previously identified through population genetic analyses as likely targets of selection during maize domestication and improvement; 46 genes previously identified as putative targets of selection also exhibit altered expression levels and coexpression relationships. We also identified 45 genes with altered, primarily higher, expression in inbred relative to outcrossed teosinte. These genes are enriched for functions related to biotic stress and may reflect responses to the effects of inbreeding. This study not only documents alterations in the maize transcriptome following domestication, identifying several genes that may have contributed to the evolution of maize, but highlights the complementary information that can be gained by combining gene expression with population genetic analyses.