Major regulatory genes in maize contribute to standing variation in teosinte (Zea mays ssp. parviglumis)

In plants, many major regulatory genes that control plant growth and development have been identified and characterized. Despite a detailed knowledge of the function of these genes little is known about how they contribute to the natural variation for complex traits. To determine whether major regulatory genes of maize contribute to standing variation in Balsas teosinte we conducted association mapping in 584 Balsas teosinte individuals. We tested 48 markers from nine candidate regulatory genes against 13 traits for plant and inflorescence architecture. We identified significant associations using a mixed linear model that controls for multiple levels of relatedness. Ten associations involving five candidate genes were significant after correction for multiple testing, and two survive the conservative Bonferroni correction. zfl2, the maize homolog of FLORICAULA of Antirrhinum, was associated with plant height. zap1, the maize homolog of APETALA1 of Arabidopsis, was associated with inflorescence branching. Five SNPs in the maize domestication gene, teosinte branched1, were significantly associated with either plant or inflorescence architecture. Our data suggest that major regulatory genes in maize do play a role in the natural variation for complex traits in teosinte and that some of the minor variants we identified may have been targets of selection during domestication.     

Evidence for a natural allelic series at the Maize domestication Locus teosinte branched1

Despite numerous quantitative trait loci and association mapping studies, our understanding of the extent to which natural allelic series contribute to the variation for complex traits is limited. In this study, we investigate the occurrence of a natural allelic series for complex traits at the teosinte branched1 (tb1) gene in natural populations of teosinte (Zea mays ssp. parviglumis, Z. mays ssp. mexicana, and Z. diploperennis). Previously, tb1 was shown to confer large effects on both plant architecture and ear morphology between domesticated maize and teosinte; however, the effect of tb1 on trait variation in natural populations of teosinte has not been investigated. We compare the effects of nine teosinte alleles of tb1 that were introgressed into an isogenic maize inbred background. Our results provide evidence for a natural allelic series at tb1 for several complex morphological traits. The teosinte introgressions separate into three distinct phenotypic classes, which correspond to the taxonomic origin of the alleles. The effects of the three allelic classes also correspond to known morphological differences between the teosinte taxa. Our results suggest that tb1 contributed to the morphological diversification of teosinte taxa as well as to the domestication of maize.     

Using Association Mapping in Teosinte to Investigate the Function of Maize Selection-Candidate Genes

Background

Large-scale screens of the maize genome identified 48 genes that show the putative signature of artificial selection during maize domestication or improvement. These selection-candidate genes may act as quantitative trait loci (QTL) that control the phenotypic differences between maize and its progenitor, teosinte. The selection-candidate genes appear to be located closer in the genome to domestication QTL than expected by chance.

Methods and Findings

As a step toward defining the traits controlled by these genes, we performed phenotype-genotype association mapping in teosinte for 32 of the 48 plus three other selection-candidate genes. Our analyses assayed 32 phenotypic traits, many of which were altered during maize domestication or improvement. We observed several significant associations between SNPs in the selection-candidate genes and trait variation in teosinte. These included two associations that surpassed the Bonferroni correction and five instances where a gene significantly associated with the same trait in both of our association mapping panels. Despite these significant associations, when compared as a group the selection-candidate genes performed no better than randomly chosen genes.

Conclusions

Our results suggest association analyses can be helpful for identifying traits under the control of selection-candidate genes. Indeed, we present evidence for new functions for several selection-candidate genes. However, with the current set of selection-candidate genes and our association mapping strategy, we found very few significant associations overall and no more than we would have found with randomly chosen genes. We discuss possible reasons that a large number of significant genotype-phenotype associations were not discovered.     

The Role of cis Regulatory Evolution in Maize Domestication

Gene expression differences between divergent lineages caused by modification of cis regulatory elements are thought to be important in evolution. We assayed genome-wide cis and trans regulatory differences between maize and its wild progenitor, teosinte, using deep RNA sequencing in F₁ hybrid and parent inbred lines for three tissue types (ear, leaf and stem). Pervasive regulatory variation was observed with approximately 70% of ∼17,000 genes showing evidence of regulatory divergence between maize and teosinte. However, many fewer genes (1,079 genes) show consistent cis differences with all sampled maize and teosinte lines. For ∼70% of these 1,079 genes, the cis differences are specific to a single tissue. The number of genes with cis regulatory differences is greatest for ear tissue, which underwent a drastic transformation in form during domestication. As expected from the domestication bottleneck, maize possesses less cis regulatory variation than teosinte with this deficit greatest for genes showing maize-teosinte cis regulatory divergence, suggesting selection on cis regulatory differences during domestication. Consistent with selection on cis regulatory elements, genes with cis effects correlated strongly with genes under positive selection during maize domestication and improvement, while genes with trans regulatory effects did not. We observed a directional bias such that genes with cis differences showed higher expression of the maize allele more often than the teosinte allele, suggesting domestication favored up-regulation of gene expression. Finally, this work documents the cis and trans regulatory changes between maize and teosinte in over 17,000 genes for three tissues.     

Do large effect QTL fractionate? A case study at the maize domestication QTL teosinte branched1

Quantitative trait loci (QTL) mapping is a valuable tool for studying the genetic architecture of trait variation. Despite the large number of QTL studies reported in the literature, the identified QTL are rarely mapped to the underlying genes and it is usually unclear whether a QTL corresponds to one or multiple linked genes. Similarly, when QTL for several traits colocalize, it is usually unclear whether this is due to the pleiotropic action of a single gene or multiple linked genes, each affecting one trait. The domestication gene teosinte branched1 (tb1) was previously identified as a major domestication QTL with large effects on the differences in plant and ear architecture between maize and teosinte. Here we present the results of two experiments that were performed to determine whether the single gene tb1 explains all trait variation for its genomic region or whether the domestication QTL at tb1 fractionates into multiple linked QTL. For traits measuring plant architecture, we detected only one QTL per trait and these QTL all mapped to tb1. These results indicate that tb1 is the sole gene for plant architecture traits that segregates in our QTL mapping populations. For most traits related to ear morphology, we detected multiple QTL per trait in the tb1 genomic region, including a large effect QTL at tb1 itself plus one or two additional linked QTL. tb1 is epistatic to two of these additional QTL for ear traits. Overall, these results provide examples for both a major QTL that maps to a single gene, as well as a case in which a QTL fractionates into multiple linked QTL.     

The role of regulatory genes during maize domestication: evidence from nucleotide polymorphism and gene expression

We investigated DNA sequence variation in 72 candidate genes in maize landraces and the wild ancestor of maize, teosinte. The candidate genes were chosen because they exhibit very low sequence diversity among maize inbreds and have sequence homology to known regulatory genes. We observed signatures of selection in 17 candidate genes, indicating that they were potential targets of artificial selection during domestication. In addition, 21 candidate genes were identified as potential targets of natural selection in teosinte. A comparison of the proportion of selected genes between our regulatory genes and genes unfiltered for their potential function (but also with very low sequence diversity among maize inbreds) provided some weak evidence that regulatory genes are overrepresented among selected genes. We detected no significant association between the positions of genes identified as potential targets of selection during domestication and quantitative trait loci (QTL) responsible for maize domestication traits. However, a subset of these genes, those identified by sequence homology as kinase/phosphatase genes, significantly cluster with the domestication QTL. We also analyzed expression profiles of genes in distinct maize tissues and observed that domestication genes are expressed on average at a significantly higher level than neutral genes in reproductive organs, including kernels.     

Genetic Analysis of Central Carbon Metabolism Unveils an Amino Acid Substitution That Alters Maize NAD-Dependent Isocitrate Dehydrogenase Activity

Background

Central carbon metabolism (CCM) is a fundamental component of life. The participating genes and enzymes are thought to be structurally and functionally conserved across and within species. Association mapping utilizes a rich history of mutation and recombination to achieve high resolution mapping. Therefore, applying association mapping in maize (Zea mays ssp. mays), the most diverse model crop species, to study the genetics of CCM is a particularly attractive system.

Methodology/Principal Findings

We used a maize diversity panel to test the CCM functional conservation. We found heritable variation in enzyme activity for every enzyme tested. One of these enzymes was the NAD-dependent isocitrate dehydrogenase (IDH, E.C. 1.1.1.41), in which we identified a novel amino-acid substitution in a phylogenetically conserved site. Using candidate gene association mapping, we identified that this non-synonymous polymorphism was associated with IDH activity variation. The proposed mechanism for the IDH activity variation includes additional components regulating protein level. With the comparison of sequences from maize and teosinte (Zea mays ssp. Parviglumis), the maize wild ancestor, we found that some CCM genes had also been targeted for selection during maize domestication.

Conclusions/Significance

Our results demonstrate the efficacy of association mapping for dissecting natural variation in primary metabolic pathways. The considerable genetic diversity observed in maize CCM genes underlies heritable phenotypic variation in enzyme activities and can be useful to identify putative functional sites.     

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

Phenotypic variation on which selection can act during evolution may be caused by variation in activity level of developmental regulatory genes. In many cases, however, such genes affect multiple traits. This situation can lead to co-evolution of traits, or evolutionary constraint if some pleiotropic effects are detrimental. Here, we present an analysis of quantitative traits associated with gene copy number of two important maize regulatory genes, the duplicate FLORICAULA/LEAFY orthologs zfl1 and zfl2. We found statistically significant associations between several quantitative traits and copy number of both zfl genes in several maize genetic backgrounds. Despite overlap in traits associated with these duplicate genes, zfl1 showed stronger associations with flowering time, while zfl2 associated more strongly with branching and inflorescence structure traits, suggesting some divergence of function. Since zfl2 associates with quantitative variation for ear rank and also maps near a quantitative trait locus (QTL) on chromosome 2 controlling ear rank differences between maize and teosinte, we tested whether zfl2 might have been involved in the evolution of this trait using a QTL complementation test. The results suggest that zfl2 activity is important for the QTL effect, supporting zfl2 as a candidate gene for a role in morphological evolution of maize.     

Fine mapping of a major quantitative trait locus that regulates pod shattering in soybean

Legumes represent the second most important family of crop plants, accounting for ~27 % of the world’s crop production. While some legumes are grown as forages or vegetables, most crop legumes are grown for harvesting their nutritious seeds. The legume seeds are contained in the pod, which is composed of a single seed-bearing carpel that, when matures, splits open along two seams, a process called pod dehiscence or pod shattering. Pod shattering before or during harvest causes yield losses of grain legumes. Moreover, the dominant shattering trait of the wild progenitors is a limiting factor for efficient introgression of value-added traits into elite breeding lines. Knowledge of the genetic mechanisms underlying pod shattering will facilitate breeding of shattering-resistant varieties, expedite introgression of agronomically favorable traits from wild species to elite breeding lines, and enrich our understanding of the evolution of seed dispersal and crop domestication in diverse crop species. Here we report fine mapping of a major quantitative trait locus (designated as qPDH1) that regulates pod shattering in soybean (Glycine max). A combination of linkage and association mapping allowed us to delimit the qPDH1 locus within a 47-kb region on chromosome 16. The data reported here will facilitate positional cloning of the underlying gene and the development of breeder-friendly genetic markers for marker-assisted selection in soybean.     

A transposable element insertion disturbed starch synthase gene SSIIb in maize

Transposable elements, as the most active genetic factors, have driven genome evolution in maize and reshaped certain key loci responsible for maize domestication, exemplified by an inserted transposon in teosinte branched1 (tb1), which controls plant architecture. In this study, we detected an insertion of a transposable element in the second exon of the coding sequence of the maize starch synthase gene SSIIb, leading to a splicing modification and gene frameshift. This insertion provided a means of determining the function of SSIIb for starch synthesis during maize domestication. Association and quantitative trait locus (QTL) mappings showed that SSIIb was not associated with starch eating quality and total starch content of kernel, and two maize near-isogenic-line-like lines with and without the insertion of the transposable element further exhibited the same starch content of kernel and leaf; in addition, nucleotide diversity analysis revealed that maize SSIIb was not under selection during domestication. All these results demonstrated that maize SSIIb might serve as a very minor genetic factor or a functional redundancy gene in starch synthesis. Global BLAST showed that the maize genome harbored 1,387 copies of this transposable element, of which 135 copies were located in genic regions. At least three genes beside maize SSIIb were disturbed by this transposable element. Five patterns of transposition, according to the insertion sites close to or within genes such as maize SSIIb in this study, are under discussion and a large quantity of present/absent variations due to the insertion of varieties of transposable elements, discovered by revolutionary next-generation sequencing, would rapidly accelerate QTL and association mappings for maize domestication through candidate gene tactics in the near future.     

Fine Mapping of a QTL Associated with Kernel Row Number on Chromosome 1 of Maize

The genetic factors underlying changes in ear morphology, and particularly the inheritance of kernel row number (KRN), have been broadly investigated in diverse mapping populations in maize (Zea mays L.). In this study, we mapped a region on the long arm of chromosome 1 containing a QTL for KRN. This work was performed using a set of recombinant chromosome nearly isogenic lines (RCNILs) derived from a BC₂S₃ population produced using the inbred maize line W22 and teosinte (Zea mays ssp. parviglumis) as the parents. A set of 48 RCNILs was evaluated in the field during the summer of 2013 in order to perform the mapping. A QTL for KRN was found that explained approximately 51% of the phenotypic variance and had a 1.5-LOD confidence interval of 203 kb. Seven genes are described in this interval. One of these candidate genes may have been the target of domestication processes in maize and contributed to the shift from two kernel row ears in teosinte to a highly polystichous ear in maize.     

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds

Background

The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants.

Results

We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves.

Conclusions

We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Electronic supplementary material

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

Genetic basis of kernel nutritional traits during maize domestication and improvement

The nutritional traits of maize kernels are important for human and animal nutrition, and these traits have undergone selection to meet the diverse nutritional needs of humans. However, our knowledge of the genetic basis of selecting for kernel nutritional traits is limited. Here, we identified both single and epistatic quantitative trait loci (QTLs) that contributed to the differences of oil and carotenoid traits between maize and teosinte. Over half of teosinte alleles of single QTLs increased the values of the detected oil and carotenoid traits. Based on the pleiotropism or linkage information of the identified single QTLs, we constructed a trait–locus network to help clarify the genetic basis of correlations among oil and carotenoid traits. Furthermore, the selection features and evolutionary trajectories of the genes or loci underlying variations in oil and carotenoid traits revealed that these nutritional traits produced diverse selection events during maize domestication and improvement. To illustrate more, a mutator distance–relative transposable element (TE) in intron 1 of DXS2, which encoded a rate‐limiting enzyme in the methylerythritol phosphate pathway, was identified to increase carotenoid biosynthesis by enhancing DXS2 expression. This TE occurs in the grass teosinte, and has been found to have undergone selection during maize domestication and improvement, and is almost fixed in yellow maize. Our findings not only provide important insights into evolutionary changes in nutritional traits, but also highlight the feasibility of reintroducing back into commercial agricultural germplasm those nutritionally important genes hidden in wild relatives.     

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

Studies that investigated the genetic basis of source and sink related traits have been widely conducted. However, the vascular system that links source and sink received much less attention. When maize was domesticated from its wild ancestor, teosinte, the external morphology has changed dramatically; however, less is known for the internal anatomy changes. In this study, using a large maize‐teosinte experimental population, we performed a high‐resolution quantitative trait locus (QTL) mapping for the number of vascular bundle in the uppermost internode of maize stem. The results showed that vascular bundle number is dominated by a large number of small‐effect QTLs, in which a total of 16 QTLs that jointly accounts for 52.2% of phenotypic variation were detected, with no single QTL explaining more than 6% of variation. Different from QTLs for typical domestication traits, QTLs for vascular bundle number might not be under directional selection following domestication. Using Near Isogenic Lines (NILs) developed from heterogeneous inbred family (HIF), we further validated the effect of one QTL qVb9‐2 on chromosome 9 and fine mapped the QTL to a 1.8‐Mb physical region. This study provides important insights for the genetic architecture of vascular bundle number in maize stem and sets basis for cloning of qVb9‐2.     

Genetic Architecture of Ear Fasciation in Maize (Zea mays) under QTL Scrutiny

Maize ear fasciation

Knowledge of the genes affecting maize ear inflorescence may lead to better grain yield modeling. Maize ear fasciation, defined as abnormal flattened ears with high kernel row number, is a quantitative trait widely present in Portuguese maize landraces.

Material and Methods

Using a segregating population derived from an ear fasciation contrasting cross (consisting of 149 F_2:3 families) we established a two location field trial using a complete randomized block design. Correlations and heritabilities for several ear fasciation-related traits and yield were determined. Quantitative Trait Loci (QTL) involved in the inheritance of those traits were identified and candidate genes for these QTL proposed.

Results and Discussion

Ear fasciation broad-sense heritability was 0.73. Highly significant correlations were found between ear fasciation and some ear and cob diameters and row number traits. For the 23 yield and ear fasciation-related traits, 65 QTL were identified, out of which 11 were detected in both environments, while for the three principal components, five to six QTL were detected per environment. Detected QTL were distributed across 17 genomic regions and explained individually, 8.7% to 22.4% of the individual traits or principal components phenotypic variance. Several candidate genes for these QTL regions were proposed, such as bearded-ear1, branched silkless1, compact plant1, ramosa2, ramosa3, tasselseed4 and terminal ear1. However, many QTL mapped to regions without known candidate genes, indicating potential chromosomal regions not yet targeted for maize ear traits selection.

Conclusions

Portuguese maize germplasm represents a valuable source of genes or allelic variants for yield improvement and elucidation of the genetic basis of ear fasciation traits. Future studies should focus on fine mapping of the identified genomic regions with the aim of map-based cloning.     

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.