#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Maize adaptation across temperate climates was obtained via expression of two florigen genes


Authors: Sara Castelletti aff001;  Aude Coupel-Ledru aff002;  Italo Granato aff002;  Carine Palaffre aff003;  Llorenç Cabrera-Bosquet aff002;  Chiara Tonelli aff001;  Stéphane D. Nicolas aff004;  François Tardieu aff002;  Claude Welcker aff002;  Lucio Conti aff001
Authors place of work: Department of Biosciences, University of Milan, Milan, Italy aff001;  LEPSE, INRAe, Univ. Montpellier, SupAgro, Montpellier, France aff002;  Unité Expérimentale du Maïs, INRAe, Univ. Bordeaux, Saint Martin de Hinx, France aff003;  GQE—Le Moulon, INRAe, Université Paris-Saclay, CNRS, AgroParisTech, Gif-sur-Yvette, France aff004
Published in the journal: Maize adaptation across temperate climates was obtained via expression of two florigen genes. PLoS Genet 16(7): e32767. doi:10.1371/journal.pgen.1008882
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008882

Summary

Expansion of the maize growing area was central for food security in temperate regions. In addition to the suppression of the short-day requirement for floral induction, it required breeding for a large range of flowering time that compensates the effect of South-North gradients of temperatures. Here we show the role of a novel florigen gene, ZCN12, in the latter adaptation in cooperation with ZCN8. Strong eQTLs of ZCN8 and ZCN12, measured in 327 maize lines, accounted for most of the genetic variance of flowering time in platform and field experiments. ZCN12 had a strong effect on flowering time of transgenic Arabidopsis thaliana plants; a path analysis showed that it directly affected maize flowering time together with ZCN8. The allelic composition at ZCN QTLs showed clear signs of selection by breeders. This suggests that florigens played a central role in ensuring a large range of flowering time, necessary for adaptation to temperate areas.

Keywords:

Plant genomics – Plant genetics – Heredity – Maize – Quantitative trait loci – Gene mapping – Genome-wide association studies – Leaves – Single nucleotide polymorphisms

Introduction

Maize was domesticated 9000 years ago in central Mexico from teosinte (Zea mays subsp. parviglumis), a tropical plant that only flowers under short days [1]. In contrast to its wild ancestor, modern maize shows a large geographic distribution from 0 to 50° N. Adaptation to temperate latitudes, with long days during summer, required breeders to first suppress the short-day requirement for floral induction, but a second adaptation was also required. Because the time to flowering strongly depends on temperature, it is shorter for a given non-photoperiodic genotype at lower than at higher latitudes of the temperate area. A long intrinsic cycle duration is necessary to counteract this effect in warmest areas, thereby avoiding decreases in cumulated photosynthesis and biomass. Conversely, a short intrinsic duration is needed in the coolest areas to avoid grain filling to occur in Autumn with low light and temperature [2]. Accordingly, the range of time to flowering in current germplasm typically ranges from 35 to 90 days after sowing in a given site [3].

About ninety quantitative trait loci (QTLs) in the maize genome account for the natural diversity in flowering time [48] and most of them explain a small proportion of the phenotypic variance, with prevalently additive effects [4]. The molecular characterization of a few large-effect alleles helped uncover key events underpinning the flowering adaptation of maize from tropical to temperate climates [9]. These alleles were found to regulate the expression of florigens, i.e. floral genes expressed in the leaf vasculature. Florigen proteins are members of the Phosphatidyl Ethanolamine Binding Proteins (PEBPs) family that include FLOWERING LOCUS T in Arabidopsis and Heading date 3a in rice that move systemically to promote flowering at the shoot apex [10]. The maize genome encodes several florigen genes, including ZEA CENTRORADIALIS 8 (ZCN8) that acts as major florigen [9,1113]. Its expression only occurs under short-day conditions in tropically adapted maize, but not in long days [13]. The uncoupling of ZCN8 expression from short day-dependent activation occurred in a stepwise manner through the selection of mutations in the ZCN8 expression regulatory processes. For example, independent transposon insertions at two related ZmCCT genes (9 and 10) reduce maize flowering photoperiod sensitivity [14,15]. The active (expressed) alleles of ZmCCT9 and 10 of tropical maize delay flowering under long days by repressing ZCN8 [13,16]. The gradual northbound expansion of maize was also driven by variations at cis-regulatory sequences of the ZCN8 promoter [9], which modify the binding landscape for different transcription factors [17]. Another ZCN8 regulatory node involves ZmMADS69, which in turn represses the ZCN8 negative regulator Related to apetala2.7 (Rap2.7) [18]. Both ZmMADS69 and Rap2.7 map within well-characterized QTLs of flowering time suggesting that allelic variations at these genes control part of the diversity of flowering time of temperate maize [1922]. Classic mutagenesis screens also identified major genes that affect flowering time, thereby helping to refine the ZCN8 regulatory pathway with additional components. These include Indeterminate1 (Id1) and delayed flowering 1 (dlf1) that encode, respectively, an upstream regulator required for ZCN8 expression [23,24] and a bZIP transcription factor presumably required for ZCN8 protein signalling at the shoot apex [13,25].

Although the above studies reveal a prominent role for ZCN8 activation in the control of flowering time in short vs long days, its contribution to the genetic variability of flowering time in non-photoperiodic temperate maize varieties is less clear. Indeed, it has been argued that the variability in flowering time in temperate maize is controlled by other pathways that act in parallel to or downstream of ZCN8 transcriptional activation [26]. Two elements challenge this view: (i) the downregulation of ZCN8 in temperate maize leads to a significant delay in flowering time (albeit less prominent compared to id1 or dlf1 mutants) [13,16,18], (ii) QTL mapping shows that flowering time is linked to allelic variations at a genomic region that harbours ZCN8 (also known as the Vgt2 QTL) [6,19,27,28].

A genome-wide study of temperate maize germplasm was necessary to reconcile these views on the potential role of florigen expression in the control of flowering time of temperate maize. Here we address this question by jointly measuring flowering time and the expression of ZCN8, and of two related genes with unknown variability, ZCN7 and ZCN12 [12], in a panel of 327 maize temperate lines. Our study revealed large variations in ZCN8 and ZCN12 transcript levels, which closely correlate to flowering time. Regulatory regions for ZCN8 and ZCN12 coincided with flowering time QTLs and we detected signatures of selection for some of these regions. Collectively, these results suggest that the adaptation of maize within temperate regions derives from the control of florigen pathways.

Results

ZCN8 levels highly correlate with flowering time in temperate maize

We first measured the time from germination to pollen shedding on florets (anthesis), to the emergence of silks from the ear (silking) and also the final leaf number which reflects the duration from germination to floral transition, in four experiments in a phenotyping platform and in the field. Large variability was observed in two platform experiments for anthesis (BLUES ranging from 56 to 99 d20°C), silking (from 57 to 101 d20°C) and leaf number (from 13 to 25) (S1A and S1B Fig), resulting in narrow-sense heritability up to 0.80 (S1 Table). Field experiments provided similar results to those obtained in the platform (Fig 1A, S1C Fig, S1 Table).

Fig. 1. ZCN8 contributes to flowering time variability in temperate maize lines.
<i>ZCN8</i> contributes to flowering time variability in temperate maize lines.
(A) Relationship between genotypic values for anthesis in field (2015–2016) and platform (2014) experiments. (B) Anthesis time plotted against ZCN8 expression after the floral transition (platform experiment in 2014). (C) Relationship between anthesis measured in the field (2015–2016) and ZCN8 expression in the platform experiment (2014). In (A), (B) and (C) n = 313 lines. Genotypic values are BLUEs. Anthesis time expressed in equivalent days at 20°C after germination. The red triangle represents genotypic values for the reference line B73. Correlation coefficients are displayed together with their significance (*** p < 10−4).

We then measured florigen expression of the same lines, in samples collected from fully expanded leaves of plants grown in the phenotyping platform. Samples were collected in the early morning at two dates, when leaves 4 and 7 reached ligulation, representing bona fide pre and post floral transition [13]. Care was taken to precisely record the harvest time of each individual sample and to record the developmental stage of every sampled plant, which slightly differed between plants because of the genetic variability of the phyllochron [29]. Mixed models including a fixed effect of sampling date and time of day (both of them small, S2A–S2D Fig), together with spatial effects in the greenhouse, were used to extract genotypic values for each line (S3C and S3D Fig).

The ZCN8 expression displayed large genetic variability. The four biological replicates, analysed at two sampling dates, provided highly correlated expression levels (up to r = 0.88) when values were normalised by those of the housekeeping genes UBIQUITIN-CONJUGATING ENZYME E2 (UCe) and ELONGATION FACTOR 1 ALPHA (EF1a) (S3A and S3B Fig). ZCN8 accumulation significantly increased between the two sampling dates (S4B Fig), consistent with previous studies [13], with higher heritability in the second than in the first date (h2 = 0.64 and 0.58, respectively). Still, we could detect a significant correlation in ZCN8 accumulation between the two sampling dates (r = 0.53, S4A Fig). We thus focused on ZCN8 expression detected at the second sampling date (leaf 7 stage).

Anthesis time was significantly correlated with the BLUEs of ZCN8 (r = -0.6, Fig 1B) and so was the total leaf number (r = -0.65, S5A Fig). The ranking of genotypes for ZCN8 was confirmed in an independent platform experiment (year 2015) on a subset of genotypes (n = 4 replicates for 276 lines, r = 0.63, S5C Fig). ZCN8 transcript levels detected in the phenotyping platform were also highly correlated with flowering times measured in field experiments (Fig 1C and S5A and S5B Fig). For any given expression level at the second sampling date, lines characterised by a strong increase in expression flowered earlier than the other lines (r = -0.36, S4C Fig). This suggests that constitutively high levels of ZCN8 and an increase in ZCN8 expression during development both contribute to early flowering time in temperate maize.

ZCN12 is a novel regulator of maize flowering and it is co-expressed with ZCN8

Although ZCN12 levels have not been documented so far for their role in flowering time, they showed high heritability (narrow-sense h2 = 0.67, S1 Table) and displayed a strong correlation with anthesis time (r = -0.64, Fig 2A) and leaf number in all platform and field experiments (S6A Fig). ZCN12 expression also tightly correlated with that of ZCN8 (r = 0.76, Fig 2B and S6B Fig), suggesting common regulatory mechanisms. This was in contrast to the ZCN7 expression that displayed low heritability (narrow sense h2 = 0.11, S1 Table) and a weaker correlation with flowering time (r = -0.31, S6C Fig) than ZCN8, despite synteny [30] and shared biochemical functions [31]. ZCN7 and ZCN8 expressions were also weakly correlated (S6C Fig). The accumulation of ZmCONZ1, a putative upstream regulator of the florigen genes in maize [32], was not linked to ZCN8 or ZCN7 expressions, nor with flowering time (S6D Fig).

Fig. 2. ZCN12 is a novel maize florigen co-regulated with ZCN8.
<i>ZCN12</i> is a novel maize florigen co-regulated with <i>ZCN8</i>.
ZCN12 expression correlates with anthesis time (platform experiment 2014), (A), and with ZCN8 expression (B). Genotypic values are BLUEs. Pearson’s correlation coefficients are indicated with their significance (***, p < 10−4). The red triangle represents genotypic values for the reference line B73. In (A) and (B), n = 173. (C) Arabidopsis T1 plants overexpressing ZCN12 but not ZCN8 show an extreme early flowering phenotype compared to empty vector transformation (Vect.); plants were photographed 3 weeks after sowing, horizontal bar, 1 cm. (D) ZCN8 and ZCN12 T1 over-expressing plants are earlier flowering than wild type (p < 0.01, one-way ANOVA with post-hoc Tukey HSD test). Horizontal bars, mean values of rosette leaves number. n = 37, 38 and 14 for, respectively, 35S:HA:ZCN8 35S:HA:ZCN12 and wild type (empty vector). (E) RT-PCR analysis on ZCN8/12 transcripts derived from independent T2 transgenic plants of the indicated genotypes. Numbers on the left are kb based on DNA ladder (Mk) migration. Amplification of vectors (Vect) harbouring genomic versions of ZCN8 and ZCN12 afforded a positive control whereas untransformed wild type and no DNA reactions (-) were used as negative controls. Actin expression was used for normalization (lower panel). Amplifications were conducted for 25 cycles.

In order to demonstrate the role of ZCN12 on flowering time, we cloned the ZCN12 gene for over-expression studies in Arabidopsis thaliana. This demonstrated a strong florigenic activity for the ZCN12 protein. Transgenic lines overexpressing ZCN12 flowered significantly earlier compared with control plants (Fig 2C–2E). Sixteen of the 40 transgenic plants showed an extremely early flowering phenotype characterised by the conversion of the shoot apical meristem in a terminal flower, with no observable vegetative leaves. The same transformation for the ZCN8 gene also caused significantly earlier flowering, but with only two cases of extreme early flowering.

Direct roles of ZCN12 and ZCN8 in the control of flowering time

We then showed that ZCN12 had a direct effect on maize flowering time. Indeed, structural equation modelling rejected models in which ZCN12 had an upstream effect by controlling ZCN8 and/or ZCN7 expressions, which would in turn control anthesis (Fig 3A and 3B, S7A and S7B Fig). Conversely, our analysis supported models in which ZCN12 accumulation had a direct effect on anthesis time, with either ZCN12 being controlled by ZCN8, or ZCN8 being controlled by ZCN12 or co-variation of both genes (Fig 3C, S7C Fig). The latter models displayed good indices of fit (RMSEA < 0.05, comparative fit index > 0.95). They were considered as equally plausible, thereby supporting the existence of a two-tier florigen system in maize in which ZCN12 acts in parallel with ZCN8.

Fig. 3. ZCN12 directly affects anthesis.
<i>ZCN12</i> directly affects anthesis.
Three families of path diagrams were tested against the whole dataset (BLUEs, platform experiment of 2014). Models in which ZCN12 has an indirect effect on anthesis (A and B) were rejected, whereas models in which both ZCN8 and ZCN12 directly affect anthesis were accepted, with an indirect effect of ZCN7 (C). Arrows, linear functional relationships between anthesis and florigen expression. Simple-headed arrows, relationships are considered as causal by the model. Double-headed arrows, free correlations. Standardized path coefficients are indicated on each arrow with the level of significance (***: P < 0.001, ns: not significant). Other models tested are displayed in S7 Fig.

ZCN8 and ZCN12 regulators map to region controlling maize flowering time

We then analysed to which extent genomic regions controlling florigens expression co-localised with regions controlling flowering time. A genome-wide association study (GWAS) identified 570 QTLs which largely overlapped between field and greenhouse experiments (S2 and S3 Tables), resulting in 214 meta-QTLs that harboured at least one QTL (Fig 4A) of flowering time, in either field or greenhouse experiments. An eQTL GWAS detected 52 meta-QTLs for ZCN8, among which half co-localised with flowering time QTLs. We also detected 85 meta-QTLs for ZCN12 and 38 for ZCN7 (Fig 4A). 12 meta-QTLs were shared between ZCN8 and ZCN12, of which 7 mapped to flowering time QTLs (Fig 4A). They also overlapped with flowering time QTLs detected in published studies [7] (S4 Table). ZCN eQTLs had a high contribution to the genetic variance of flowering time as shown by mixed-models applied to individual experiments (S5 Table). Indeed, in the platform experiment of 2014, the 19 ZCN eQTLs that co-localised with anthesis QTLs accounted for 57% of the genetic variance of anthesis in the platform experiment (vs 75% for all 43 QTLs of anthesis). Similar results were found when considering anthesis in the field, or total leaf number in each experiment (S5 Table).

Fig. 4. ZCN eQTLs are major contributors to the genetic variance of flowering time.
<i>ZCN</i> eQTLs are major contributors to the genetic variance of flowering time.
(A) ZCN-related QTLs overlap with flowering time QTLs. Venn diagram highlights ‘meta-QTLs’ shared among flowering-related traits and ZCN transcripts accumulation. “Flowering” meta-QTLs are defined as regions containing at least one QTL among anthesis, silking, total number of leaves in at least one experiment (platform, field, years). “ZCN-” meta-QTLs are defined as regions containing at least one ZCN- eQTL in at least one experiment (platform 2014–2015) and leaf stage. (B) Regional association plot for one meta-QTL on chromosome 3. Distribution of the -log10(pval) for all variants in the region. The dotted grey line corresponds to–log10(pval) = 4.5 for anthesis-related traits evaluated in the field (purple) and platform (orange) and 4 for the expression traits evaluated in the platform (green). (C), linkage disequilibrium (LD) heat map of all SNPs in the QTL showing the local LD (r2) between all the variants; Black lines represent the distribution of the SNPs and the blue boxes represents the genes mapped for the region. Red asterisks represent the position of the most significant SNPs. (See also S4 Table for candidate genes). Grey areas represent sub-regions harbouring the most significant SNPs and genes and having low LD between them on average.

We inspected the ZCN eQTLs for the presence of candidate genes or QTLs already known to control flowering time. In the Bin 3.05 hotspot region known for flowering time regulation [21,33], we identified five ZCN-related eQTLs which co-localised with at least one QTL of flowering time (S3 Table). An eQTL of ZCN8/12 (158.73–159.93Mb) was located at a genomic position compatible with the newly described transcriptional regulator ZmMADS69 (GRMZM2G171650) [18,21,34] (Fig 4B and 4C). Another eQTL region (Bin 8.05) comprised the gene GRMZM2G700665 (also known as Rap2.7) (S8B Fig and S4 Table), encoding a floral repressor [19] and negative regulator of ZCN8 [18].

Analysis of published whole-genome expression datasets [35] was used to check if transcripts whose expression correlated with ZCN8 were over-represented at ZCN eQTLs positions. We confirmed significant correlations between anthesis and ZCN8/12/7 expression (n = 181 lines, r = up to - 0.54) (S9A–S9C Fig). Furthermore, these datasets revealed strong patterns of co-regulation not only between ZCN8 and ZCN12 but also between ZCN8 and ZCN7 (n = 204 lines, r = up to 0.76) (S9D–S9F Fig). Additional 549 transcripts showed a significant correlation with ZCN8 expression levels (S6 Table). These were parsed in two clusters composed of genes with opposite trends to ZCN8 gene expression and correspondingly opposite association with anthesis (S9G Fig). Neither ZmMADS69 nor Rap2.7 nor any other gene in the eQTLs regions within Bin 3.05 were co-expressed with ZCN8, and co-expressed transcripts were not particularly enriched in eQTLs regions (S4 Table).

Signs of selection at ZCN8 and ZCN12 eQTLs to control maize flowering time

The genetic groups in the studied panel showed markedly different values of ZCNs expression levels and anthesis dates (S10A,S10B, S10D and S10E Fig), and clear co-expression of ZCN8 and ZCN12 within each group (S10C and S10F Fig). The observed differences in flowering time between admixture groups may be due to selection by past breeding of temperate maize lines. Consistent with this hypothesis, we found significant differences, between the seven genetic groups, in allelic distribution at QTLs showing co-location between ZCN8/12 expression and flowering time (Fig 5A). Alleles for increased expression of ZCN8 and ZCN12 were more frequent in genetic groups bred for early maturity (Fig 5A). Hence, a pool of early maturity alleles was represented within all genetic groups, but in different proportions, suggesting that some combinations were preferentially selected in different groups.

Fig. 5. Flowering time precocity of temperate maize derives from selection at ZCN8/12 regulatory regions.
Flowering time precocity of temperate maize derives from selection at <i>ZCN8/12</i> regulatory regions.
(A) Allelic distributions within genetic groups for 8 major QTLs. The SNPs displayed belong to the core set of 19 QTLs common to anthesis in the platform experiment in 2014, and ZCN transcript accumulations. At each locus the orange colour corresponds to the 'late allele' (i.e. conferring delayed anthesis), which is for all 8 SNPs the B73 allele. Genetic groups are listed according to the flowering time phenology (late to early, top to bottom). (B) Genome-wide distribution of FST Bayescan (FSTBS) along chromosome 3 and 8. SNPs under selection are highlighted in different colors to reflect the different level of confidence for selection (Black: low confidence, Orange: substantial, Red: strong). Regions shaded in blue correspond to meta-QTLs intervals.

We further tested for genome-wide deviation from neutrality in our panel for the chromosomes 3 and 8 harbouring major meta-QTLs. Fst outliers approach using BayeScan highlighted SNPs under selection in the ZCN-related QTLs and other regions that might contribute to flowering time (Fig 5B and S7 Table). When comparing all six populations (excluding the admixture group ‘diverse dent’), the BayeScan revealed SNPs under selection ranging from a log10(PO)—posterior odds—of 0.5 to 1 (substantially significant) to higher than 1 (highly significant). Most importantly, this analysis revealed SNPs under selection comprised in the ZCN-related QTLs regions, further indicating the contribution of ZCN8/12 regulation to flowering variability.

Discussion

ZCN12 and ZCN8 are co-regulated by trans-acting factors to control temperate maize flowering

Our study provides an example for how shifts in reproductive development that enable local adaptation can emerge through different florigen levels. Besides confirming the role of ZCN8 in flowering time regulation we also uncover the role of ZCN12 in this process. Interestingly ZCN12 is more distantly related to ZCN8 than ZCN7 [12] and located in a region which is not syntenic to ZCN8 [30], yet it underwent a similar regulatory trajectory to ZCN8 as shown by the strong co-variation across different lines and experiments.

GWAS approaches support the major contribution of trans-acting factors in driving the expression of ZCN8 and ZCN12 in temperate maize. We also detected a SNP associated with anthesis and silking located approximately 1.7 kbp upstream of ZCN8 (i.e. AX-91100620, S2 Table). However, this and eight other SNPs located in this putative regulatory region were not associated with ZCN8 transcript levels. It is unlikely that the absence of cis-type eQTLs for ZCN8 could depend on insufficient coverage of the genetic diversity, since they were discovered on a large set of diverse inbred lines: 32,758 accessions for 6 GBS markers, [36] and 30 accessions for 3 Affymetrix Axiom markers [34]. Multiple trans-acting genes and complex long-range regulatory interactions are thus more likely responsible for the accumulation of ZCN8/ZCN12. Supporting this hypothesis, we detected trans-eQTLs for ZCN8/ZCN12 that colocalize with two genomic regions—Vgt1 and Vgt3—carrying genes involved in flowering time variation (Fig 4B, S8B Fig). These regions include transcriptional regulators known to be under selection in maize latitudinal adaptation, namely ZmMADS69/Vgt3 [18,21,33] and Rap2.7/Vgt1 [19,22]. We did not detect any obvious candidate genes associated with photoperiodic signalling (e.g. ZmCCT10), suggesting its small contribution in ZCN8/12 regulation in our panel.

ZmMADS69 and Rap2.7 did not show significant trends of co-expression with their trans-regulated targets ZCN8/12. It is thus possible that that there are many more genes of smaller effect contributing to ZCN variations which were not captured in our GWAS. These loci could affect the activity or accumulation of some of the abovementioned transcription factors. For example, microRNA172s are known to interfere with the translation of APETALA2 genes [37] (such as Rap2.7), indicating that transcript accumulation of Rap2.7 may not reflect its gene product levels. Another not mutually exclusive explanation may be that regulatory networks that impact on ZCNs levels are transient or highly tissue-specific, which would make it difficult to compare independent expression datasets.

We found limited numbers of genes located within eQTLs regions and co-expressed with ZCN8 and 12 (S4 Table). These included a putative RNA binding protein GRMZM2G107491 (underlying a ZCN12-specific eQTL), and two chloroplast-associated proteins (underlying different ZCN8 eQTLs), implying carbon metabolism in ZCN8 regulation. Carbon status is linked to the transcriptional regulation of the Arabidopsis florigen FLOWERING LOCUS T [38] independent of photoperiod conditions. Moreover, recent studies indicate that temperate maize flowering is associated with expression of genes involved in carbon partitioning (which provides a different molecular signature from photoperiodic flowering induction) [24]. We note that patterns of co-variations may not necessarily define causal, but also indirect relationships. Also, since ZCN8 levels are linked to flowering, our list of ZCN8 co-expressed genes would lead to overemphasizing the association between gene expression and flowering. Thus, co-expressed genes might also be associated with a different, more biologically distal phenotype. Overall, our data support the contribution of carbon metabolism upstream of ZCN8 and provide an initial framework for the identification of photoperiod independent mechanisms that shape ZCN8 levels across temperate maize varieties.

ZCN8/12 levels underpin independent events of selection for flowering adaptation

A large body of evidence highlights the importance of the regulation of gene expression in maize evolution and adaptation [39,40] and the contribution of the expressed genome to the heritable phenotypic variance [41]. It is hypothesised that temperate maize is the product of a northbound expansion through North America which was driven by the accumulation of an increasing number of early flowering alleles [42,43]. Rap2.7 and many other loci were under selection during the adaptation to high latitudes [22,27] and general patterns of differentiation for flowering time precocity closely followed the proportion of early alleles found in northern flints (earliest), corn belt dents (intermediate) and tropical maize genomes (late maturity). Besides revealing that at least some of these loci exert their effects on flowering time by affecting ZCNs levels, our study indicates that diversification in flowering can be imposed by combinatorial allelic effects in a relatively small number of ZCN-related eQTLs regions. While flowering time may represent the main trait under selection, these alleles (through their effect on florigen levels) may also affect other phenotypic outputs, for example florigen levels have been related to inflorescence development in other grasses [44]. We finally provide evidence for different signatures of selection in the main regions previously mapped in the genetic control of flowering time, thus confirming their potential contribution in the adaptation process.

Conclusion

We show that florigen genes play an essential role in the genetic variability of flowering time in those lines in which breeding suppressed photoperiodism. We reveal the role of a novel florigen, ZCN12, a non-syntenic paralog of ZCN8, whereas other related candidates had no clear role. Our study points to combinations of trans-type regulatory mechanisms underpinning selection for flowering time in different genetic groups through florigen levels.

Materials and methods

Plant material

A maize collection of 354 lines was assembled in the frame of the DROPS (European project FP7-244374) and Amaizing projects (French National Agency ANR-10-BTBR-01). It consists of historical and breeding American and European dent lines from the main founder group of modern hybrids [33]. Lines were genotyped using 50K Infinium HD Illumina array [45], a 600K Axiom Affymetrix array [34] and a set of 500K markers obtained by Genotyping by Sequencing (GBS) [33]. Genotyping data were filtered and missing data were imputed using Beagle V3.3.2 [46] after assembling three sub-matrices and removing duplicate loci based on their physical position. After quality control, 758,863 polymorphic single‐nucleotide polymorphisms (SNPs) were retained for GWAS analyses. All physical positions referred to hereafter are based on the B73 reference genome [47] RefGen_V2. Based on the 50K [33], seven genetic groups were identified using the ADMIXTURE software. These groups were constituted by: (i) 53 lines in the Non Stiff Stalk (Iodent) family traced by PH207, (ii) 18 European dent in the one traced by D06, (iii) 38 lines in the Lancaster family traced by Mo17, (iv) 17 in the one traced by Oh43, (v) 39 lines in the stiff stalk family traced by B73, (vi) 43 lines in the one traced by B14a, and (vii) 146 lines that did not fit into the six primary heterotic groups, such as W117, NC358, and F252 (referred to as “diverse dent”). B73 was considered as the reference line [47].

Assessment of flowering related traits in the field

The panel was assessed in two consecutive field experiments for flowering related traits. Experiments were conducted at Saint Martin de Hinx (INRA experimental station, 43°N, 1.3°W, France) in spring-summer 2015 and 2016, using alpha-lattice designs with two replicates under well-watered conditions piloted by soil sensors. Male (anthesis) and female (silking) flowering dates, as well as the total number of leaves (including the ones lost) produced were measured for all the lines. Air temperature was measured every hour in each experiment at 2m height over a reference grass canopy, and the progression of the crop cycle site was characterized via thermal time after emergence, expressed in equivalent days at 20°C (d20°C) [48].

Flowering-related traits measurements in the platform

Two subpanels of 327 lines and 276 lines were grown in the Phenoarch platform [49] (https://www6.montpellier.inrae.fr/lepse_eng/M3P) respectively in autumn 2014 and spring 2015 under natural fluctuating daylength conditions (with a minimum daylength of 14 hours). The two subpanels were optimized using the algorithm described in Rincent et al. [50] that allows to maximise genetic information by discarding individuals that might be more easily predictable. Five plants per line were grown under well-watered conditions and daily assessed for developmental traits. The number of visible and ligulated leaves of every plant was scored twice a week. Anthesis, silking (in d20°C) and the total number of leaves produced (including the ones lost) were assessed for all the lines.

Tissue sampling for transcripts expression analysis

Maize leaves were sampled at two developmental stages, ligulated leaf 4 and 7 (L4 and L7, respectively) from plants growing in the platform. The sampling day was decided based on B73 stage, which was used as a morphological reference. Four plants were sampled for each line. Samples deprived of the midvein were collected in the distal portion of fully expanded leaves (with clearly visible ligule) of each plant at approximately the same time of the day (within 90 minutes from sunrise). Sampling was conducted over three consecutive days for each sampling stage, and the fourth sampling was done soon after the third one. The collected material was kept in liquid nitrogen and the tubes were stored at -80°C.

RNA extraction and cDNA synthesis

The frozen samples were ground at high speed with glass beads in a TissueLyser II apparatus (Qiagen). Qiazol reagent (Qiagen) was added to individual tubes and total RNA was extracted using the Direct-zol™-96 kit (Zymo Research); on-column DNase I digestion was performed following the kit protocol. Total RNA was then quantified using a NanoDrop One instrument (ThermoFisher) and 1000 ng of each RNA sample were used for retrotranscription with the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). The cDNA samples were diluted 1:5 with water before proceeding to qPCR experiments.

TAQMAN and SYBR green assays design, implementation and analysis

The genomic and coding sequences of the targets GRMZM2G179264 (ZCN8), GRMZM2G141756 (ZCN7), GRMZM2G103666 (ZCN12) and GRMZM2G405368 (ZmCONZ1) and that of reference genes GRMZM2G102471 (UBIQUITIN CARRIER PROTEIN) [51], and GRMZM2G153541 (ELONGATION FACTOR 1Α) [52], were retrieved from the MaizeGDB website. Probe design was done by Bio-Rad specialists and PrimePCR Probe Custom Assay (Bio-Rad) were synthesized accordingly; the main features of probes and primers are listed in S8 Table. qPCR experiments were carried out in a CFX96 Real-time PCR System (Bio-Rad) using SsoAdvanced Universal Probes Supermix (Bio-Rad) and following manufacturer’s specifications. Amplification curves were evaluated using the instrument's software (Vs 3.1) and relative gene expression was calculated by applying the ΔCt method [5355]. Gene expression was estimated by averaging the ΔCt values computed against either GRMZM2G102471, GRMZM2G153541 or the mean of both as references. Expression in the 2015 experiment was evaluated using solely GRMZM2G153541 as a reference. ZCN12 expression levels derived from 173 of the 327 lines studied, because the detection assay was impaired by a three-nucleotide deletion in the 5’UTR of ZCN12 coinciding with the forward primer annealing site (S11 Fig). Nevertheless, these 173 lines distributed across the genetic groups of our panel in similar proportions as the initial 327 lines. This small deletion had no effects on ZCN12 transcript accumulation as further assayed by qPCR using the QuantiNova SYBR Green PCR Kit (Qiagen) on a subpopulation of lines characterised by the presence/absence of the deletion.

Mixed-models, BLUEs and heritability calculation

Analyses were performed using the R software version 3.4.3 [56]. Preliminary data exploration (ANOVA) was conducted to assess the effects of sampling day and time, as well as the actual developmental stage at sampling, on transcripts accumulation (S2 Fig). For all traits measured in the field and platform experiments (flowering-related traits and transcript accumulation), models were selected among several mixed-models to calculate genotypic means. Models were fitted with ASReml-R [57] (version 3), including the line and replicate effects as fixed, random spatial effects, spatially correlated errors, and other fixed effects (year for the field experiments; sampling date or sampling period for transcript accumulation in the platform). For each trait measured in the field and platform experiments (flowering-related traits and transcript accumulation), the best mixed-model was selected to estimate the best linear unbiased estimations (BLUEs) of the genotypic means, which were then used in the rest of the analyses. The same model, but considering the line effects as random, was used to estimate variance components, which were used to calculate broad-sense heritability. Narrow-sense heritability was estimated with a model assuming additive SNP effects using the R-package Heritability [58] (version 1.2) and a relatedness matrix as in Millet et al., [59].

Structural equation modelling

Structural equation modelling is a generalized method for the analysis of covariance relationships and is used to evaluate the fit of data to a priori causal hypotheses about the functioning of a system [6062]. These multivariate hypotheses are represented as graphical path models. Structural equation modelling then allows the assessment of the degree of fit between the observed and expected covariances structures, which is expressed as a goodness-of-fit χ2. Here, the aim was to impose a theoretical structure relating the direct and indirect relationships between ZCN8, ZCN12 and ZCN7 transcripts accumulation and flowering taking into account the results of bivariate correlations. BLUEs from the platform experiment (2014) were used and the assumption of univariate as well as multivariate normality were verified. Network structures tested were constructed as follow: (i) anthesis was placed downstream of genes expression based on previous knowledge (flowering is a consequence of florigens expression), (ii) then testing for all possible network structures between the 3 florigens acting upstream of anthesis. Models tested are detailed in the results section. Structural equation models were tested in R using the lavaan structural equation modelling package [63] (version 0.6–5), which uses the standard maximum likelihood estimator. A significant goodness-of-fit χ2 statistic indicates that the model does not fit the data. Once a model has not been rejected and considered biologically plausible, parameter estimates can be used to study direct, as well as indirect, effects of the variables. Standardized path coefficients quantify the strength of a relationship, whereas the effects of the other variables are held constant. Parameter estimates are tested for significance using z statistics. Root mean square error approximation (RMSEA) and comparative fit index (CFI) indices are used to assess the closeness of fit. Good models have a RMSEA < 0.05 and CFI > 0.95.

GWAS analysis

GWAS was performed on individual traits for each experiment as previously described [59]. Briefly, we used the single locus mixed model:

where Y is the vector of phenotypic values, μ the overall mean, X is the vector of SNP scores, β is the additive effect, and G and E represent random polygenic and residual effects. As in Rincent et al. [64], the variance-covariance matrix of G was determined by a genetic relatedness (or kinship) matrix, derived from all SNPs except those on the chromosome containing the SNP being tested. The SNP effects β were estimated by generalized least squares, and their significance (H0: β = 0) tested with an F-statistic. This model accounting for relatedness was found the best to control for confounding factors as compared to models accounting for structure or relatedness plus structure. Analyses were performed with FaST-LMM v2.07 [65]. Physical positions of significant SNPs were projected on the consensus genetic map for Dent genetic material [66]. Candidate SNPs distant less than 0.1 cM were considered as belonging to a common QTL, described via the most significant SNP in the QTL and the interval between all SNPs belonging to the QTL. Co-localizations between QTLs (across traits and experiments) were identified using the same 0.1cM window.

Regional plot definition

Annotated genes located within QTLs were identified according to maize annotation version 2 (MaizeGDB). For the entire QTL, the pairwise LD (R2) between all variants was estimated using the r-package snpStat [67] (version 1.33.0) and plotted with LDheatmap [68] (version 0.99–7) here shown as a triangle plot. Based on the LD pattern and the levels of significance of SNPs, the QTL were split into sub-regions.

Exploration of selective pressure for flowering-related traits and transcripts accumulation

Allelic distributions within the population were explored at SNPs representative of the core-set of QTLs and revealed strong contrasts across genetic groups. An FST study was then undertaken for identifying outlier loci that would differentiate the genetic groups more than the average genomic background. Bayescan version 2.1 [69] decomposes the Locus-population FST into a component specific for the population shared by all loci (beta) and a component specific for the locus shared by all the populations (alpha). Departure from neutrality at any given locus is assumed when alpha is significantly different from zero. Decisions about the chance of each locus being under selection are made based on posterior odds (PO) comparing the models with and without alpha component. We used as evidence of SNPs under selection markers with log10(PO) > 0.5 which corresponds to markers considered substantially significant [27].

Re-analysis of RNAseq data

Expression counts relative to RNA-seq experiments and flowering time [35] were downloaded from the Cyverse Discovery Environment under the directory: https://datacommons.cyverse.org/browse/iplant/home/shared/commons_repo/curated/Kremling_Nature3RNASeq282_March2018. We used data for post-sexual maturity leaves during the day (LMAD) because expression counts for ZCN7/8/12 at earlier stages (e.g. leaf 3) were extremely low and not available for all lines. We included transcripts accumulation dataset for 204 dent lines (with no biological replicates), of which 182 had also data available for anthesis (BLUPs). We discarded extreme outliers (respectively ZCN8 > 400, ZCN12 > 350 and ZCN7 > 150 expression counts). We thus obtained correlations between the expression counts of each annotated gene in the maize genome and ZCN8 was calculated (correlation coefficient (r) and associated p-value). The same method was applied to establish the relationship between the expression of each gene and anthesis data. Co-expression data were further reduced by applying a Bonferroni correction and analysed with a k-means clustering approach [70].

Gene cloning, transformation and growth of Arabidopsis

ZCN12 and ZCN8 were amplified from genomic DNA from B73 with oligonucleotides lcm116, 117, 118 and 119 as detailed in S7 Table. These primers incorporate sequences for Gateway-based cloning (Invitrogen) and include start and stop codons. PCR products were subject to a second round of PCR amplification with adapters art77 and art78. Gel-purified PCR products were cloned into Gateway pDONR221 vector through BP-mediated recombination and sequenced. ZCN12 and ZCN8 entry clones were finally recombined into the binary destination vector pEarleyGate201 [71] via LR clonase. We used the floral dip method to transform wild-type Columbia Arabidopsis and primary transformants were selected on soil by continuous applications of Basta. T1 plants that yielded enough seeds were analysed in T2 to monitor transgene-derived expression by RT-PCR. Plants were grown in a climatic chamber under long-day conditions (16 hours of light, 8 hours of dark), a mean daily temperature of 23 degrees (night-day) and 60% humidity, light was provided by fluorescent tubes and PAR (photosynthetically active radiation) was approximately 90 μmol m-2 s-1.

Supporting information

S1 Fig [a]
A wide flowering time variability conserved across different environments.

S2 Fig [a]
transcript accumulation is not consistently affected by day of sampling, time of sampling within a day, or actual leaf stage at sampling.

S3 Fig [a]
Three methods of normalization of transcript accumulation show highly consistent results.

S4 Fig [c]
accumulation increases between stage ligulated leaves 4 and 7 in a genotype-dependent fashion.

S5 Fig [a]
Strong correlation between accumulation and flowering time is confirmed across different environments.

S6 Fig [a]
Strong correlation between accumulation and maize flowering is confirmed across different environments whereas and have a weak or no effect.

S7 Fig [a]
Path diagrams indicate possible direct routes of to activate flowering.

S8 Fig [a]
eQTLs are major contributors to the genetic variance of flowering time in the platform experiment of 2014.

S9 Fig [b]
Strong patterns of co-regulation are verified in independent datasets.

S10 Fig [a]
and accumulation follows the flowering time precocity of different admixture groups.

S11 Fig [a]
A polymorphism in the region affects amplification efficiency.

S12 Fig [1a]
Graphical abstract of the experiments conducted in this study and of the associated datasets.

S1 Table [xlsx]
Statistical indicators of the genotypic variability for flowering related traits and florigens transcripts accumulation in individual experiments.

S2 Table [chr]
Complete set of QTLs for flowering related traits and florigens transcript accumulation.

S3 Table [xlsx]
Set of ‘meta-QTLs’ for flowering-related traits and florigen transcript accumulation.

S4 Table [ds]
Core set of QTLs for flowering related traits and florigens transcript accumulation.

S5 Table [g]

S6 Table [r]
List of coregulated genes and their relation to flowering time.

S7 Table [po]
Fst analysis for SNPs under selection.

S8 Table [xlsx]
Primers and assays information used in this study.


Zdroje

1. Matsuoka Y, Vigouroux Y, Goodman MM, Sanchez G. J, Buckler E, Doebley J. A single domestication for maize shown by multilocus microsatellite genotyping. Proc Natl Acad Sci. 2002;99: 6080–6084. doi: 10.1073/pnas.052125199 11983901

2. Parent B, Leclere M, Lacube S, Semenov MA, Welcker C, Martre P, et al. Maize yields over Europe may increase in spite of climate change, with an appropriate use of the genetic variability of flowering time. Proc Natl Acad Sci. 2018;115: 10642–10647. doi: 10.1073/pnas.1720716115 30275304

3. Colasanti J, Muszynski M. The Maize Floral Transition. Handbook of Maize: Its Biology. 2009. pp. 41–55. doi: 10.1007/978-0-387-79418-1_3

4. Buckler ES, Holland JB, Bradbury PJ, Acharya CB, Brown PJ, Browne C, et al. The genetic architecture of maize flowering time. Science. 2009;325: 714–718. doi: 10.1126/science.1174276 19661422

5. Li D, Wang X, Zhang X, Chen Q, Xu G, Xu D, et al. The genetic architecture of leaf number and its genetic relationship to flowering time in maize. New Phytol. 2016;210: 256–268. doi: 10.1111/nph.13765 26593156

6. Chardon F, Virlon B, Moreau L, Falque M, Joets J, Decousset L, et al. Genetic architecture of flowering time in maize as inferred from quantitative trait loci meta-analysis and synteny conservation with the rice genome. Genetics. 2004;168: 2169–2185. doi: 10.1534/genetics.104.032375 15611184

7. Li YX, Li C, Bradbury PJ, Liu X, Lu F, Romay CM, et al. Identification of genetic variants associated with maize flowering time using an extremely large multi-genetic background population. Plant J. 2016;86: 391–402. doi: 10.1111/tpj.13174 27012534

8. Romero Navarro JA, Willcox M, Burgueño J, Romay C, Swarts K, Trachsel S, et al. A study of allelic diversity underlying flowering-time adaptation in maize landraces. Nat Genet. 2017;49: 476–480. doi: 10.1038/ng.3784 28166212

9. Guo L, Wang X, Zhao M, Huang C, Li C, Li D, et al. Stepwise cis-Regulatory Changes in ZCN8 Contribute to Maize Flowering-Time Adaptation. Curr Biol. 2018;28: 3005–3015.e4. doi: 10.1016/j.cub.2018.07.029 30220503

10. Kobayashi Y, Weigel D. Move on up, it’s time for change—mobile signals controlling photoperiod-dependent flowering. Genes Dev. 2007;21: 2371–2384. doi: 10.1101/gad.1589007 17908925

11. Lazakis CM, Coneva V, Colasanti J. ZCN8 encodes a potential orthologue of Arabidopsis FT florigen that integrates both endogenous and photoperiod flowering signals in maize. J Exp Bot. 2011;62: 4833–4842. doi: 10.1093/jxb/err129 21730358

12. Danilevskaya ON, Meng X, Hou Z, Ananiev E V, Simmons CR. A genomic and expression compendium of the expanded PEBP gene family from maize. Plant Physiol. 2008;146: 250–264. doi: 10.1104/pp.107.109538 17993543

13. Meng X, Muszynski MG, Danilevskaya ON. The FT-like ZCN8 gene functions as a floral activator and is involved in photoperiod sensitivity in maize. Plant Cell. 2011;23: 942–960. doi: 10.1105/tpc.110.081406 21441432

14. Yang Q, Li Z, Li W, Ku L, Wang C, Ye J, et al. CACTA-like transposable element in ZmCCT attenuated photoperiod sensitivity and accelerated the postdomestication spread of maize. Proc Natl Acad Sci. 2013;110: 16969–16974. doi: 10.1073/pnas.1310949110 24089449

15. Hung H-Y, Shannon LM, Tian F, Bradbury PJ, Chen C, Flint-Garcia SA, et al. ZmCCT and the genetic basis of day-length adaptation underlying the postdomestication spread of maize. Proc Natl Acad Sci. 2012;109: E1913–21. doi: 10.1073/pnas.1203189109 22711828

16. Stephenson E, Estrada S, Meng X, Ourada J, Muszynski MG, Habben JE, et al. Over-expression of the photoperiod response regulator ZmCCT10 modifies plant architecture, flowering time and inflorescence morphology in maize. PLoS One. 2019;14. doi: 10.1371/journal.pone.0203728

17. Alter P, Bircheneder S, Zhou L-Z, Schlüter U, Gahrtz M, Sonnewald U, et al. Flowering Time-Regulated Genes in Maize Include the Transcription Factor ZmMADS1. Plant Physiol. 2016;172: 389–404. doi: 10.1104/pp.16.00285 27457125

18. Liang Y, Liu Q, Wang X, Huang C, Xu G, Hey S, et al. Zm MADS 69 functions as a flowering activator through the ZmRap2.7‐ ZCN 8 regulatory module and contributes to maize flowering time adaptation. New Phytol. 2019;221: 2335–2347. doi: 10.1111/nph.15512 30288760

19. Salvi S, Sponza G, Morgante M, Tomes D, Niu X, Fengler KA, et al. Conserved noncoding genomic sequences associated with a flowering-time quantitative trait locus in maize. Proc Natl Acad Sci. 2007;104: 11376–11381. doi: 10.1073/pnas.0704145104 17595297

20. Castelletti S, Tuberosa R, Pindo M, Salvi S. A MITE Transposon Insertion Is Associated with Differential Methylation at the Maize Flowering Time QTL Vgt1. G3; Genes|Genomes|Genetics. 2014;4: 805–812. doi: 10.1534/g3.114.010686 24607887

21. Hirsch CN, Foerster JM, Johnson JM, Sekhon RS, Muttoni G, Vaillancourt B, et al. Insights into the Maize Pan-Genome and Pan-Transcriptome. Plant Cell. 2014;26: 121–135. doi: 10.1105/tpc.113.119982 24488960

22. Ducrocq S, Madur D, Veyrieras J-BB, Camus-Kulandaivelu L, Kloiber-Maitz M, Presterl T, et al. Key Impact of Vgt1 on Flowering Time Adaptation in Maize: Evidence From Association Mapping and Ecogeographical Information. Genetics. 2008;178: 2433–2437. doi: 10.1534/genetics.107.084830 18430961

23. Colasanti J, Yuan Z, Sundaresan V. The indeterminate gene encodes a zinc finger protein and regulates a leaf-generated signal required for the transition to flowering in maize. Cell. 1998;93: 593–603. doi: 10.1016/s0092-8674(00)81188-5 9604934

24. Minow MAA, Ávila LM, Turner K, Ponzoni E, Mascheretti I, Dussault FM, et al. Distinct gene networks modulate floral induction of autonomous maize and photoperiod-dependent teosinte. J Exp Bot. 2018;69: 2937–2952. doi: 10.1093/jxb/ery110 29688423

25. Muszynski MG, Dam T, Li B, Shirbroun DM, Hou Z, Bruggemann E, et al. delayed flowering1 Encodes a Basic Leucine Zipper Protein That Mediates Floral Inductive Signals at the Shoot Apex in Maize. Plant Physiol. 2006;142: 1523–1536. doi: 10.1104/pp.106.088815 17071646

26. Dong Z, Danilevskaya O, Abadie T, Messina C, Coles N, Cooper M. A gene regulatory network model for Floral transition of the shoot apex in maize and its dynamic modeling. PLoS One. 2012;7: e43450. doi: 10.1371/journal.pone.0043450 22912876

27. Bouchet S, Servin B, Bertin P, Madur D, Combes V, Dumas F, et al. Adaptation of Maize to Temperate Climates: Mid-Density Genome-Wide Association Genetics and Diversity Patterns Reveal Key Genomic Regions, with a Major Contribution of the Vgt2 (ZCN8) Locus. PLoS One. 2013;8: e71377. doi: 10.1371/journal.pone.0071377 24023610

28. Romay MC, Millard MJ, Glaubitz JC, Peiffer JA, Swarts KL, Casstevens TM, et al. Comprehensive genotyping of the USA national maize inbred seed bank. Genome Biol. 2013;14: R55. doi: 10.1186/gb-2013-14-6-r55 23759205

29. Millet EJ, Kruijer W, Coupel-Ledru A, Alvarez Prado S, Cabrera-Bosquet L, Lacube S, et al. Genomic prediction of maize yield across European environmental conditions. Nature Genetics. 2019. doi: 10.1038/s41588-019-0414-y

30. Zhang Y, Ngu DW, Carvalho D, Liang Z, Qiu Y, Roston RL, et al. Differentially Regulated Orthologs in Sorghum and the Subgenomes of Maize. Plant Cell. 2017;29: 1938–1951. doi: 10.1105/tpc.17.00354 28733421

31. Mascheretti I, Turner K, Brivio RS, Hand A, Colasanti J, Rossi V. Florigen-Encoding Genes of Day-Neutral and Photoperiod-Sensitive Maize Are Regulated by Different Chromatin Modifications at the Floral Transition. Plant Physiol. 2015;168: 1351–1363. doi: 10.1104/pp.15.00535 26084920

32. Miller TA, Muslin EH, Dorweiler JE. A maize CONSTANS-like gene, conz1, exhibits distinct diurnal expression patterns in varied photoperiods. Planta. 2008;227: 1377–1388. doi: 10.1007/s00425-008-0709-1 18301915

33. Negro SS, Millet EJ, Madur D, Bauland C, Combes V, Welcker C, et al. Genotyping-by-sequencing and SNP-arrays are complementary for detecting quantitative trait loci by tagging different haplotypes in association studies. BMC Plant Biol. 2019;19: 318. doi: 10.1186/s12870-019-1926-4 31311506

34. Unterseer S, Bauer E, Haberer G, Seidel M, Knaak C, Ouzunova M, et al. A powerful tool for genome analysis in maize: Development and evaluation of the high density 600 k SNP genotyping array. BMC Genomics. 2014;15. doi: 10.1186/1471-2164-15-823

35. Kremling KAG, Chen S-Y, Su M-H, Lepak NK, Romay MC, Swarts KL, et al. Dysregulation of expression correlates with rare-allele burden and fitness loss in maize. Nature. 2018;555: 520–523. doi: 10.1038/nature25966 29539638

36. Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, et al. TASSEL-GBS: A high capacity genotyping by sequencing analysis pipeline. PLoS One. 2014. doi: 10.1371/journal.pone.0090346

37. Chen X. A MicroRNA as a Translational Repressor of APETALA2 in Arabidopsis Flower Development. Science. 2004. doi: 10.1126/science.1088060

38. Wahl V, Ponnu J, Schlereth A, Arrivault S, Langenecker T, Franke A, et al. Regulation of flowering by trehalose-6-phosphate signaling in Arabidopsis thaliana. Science. 2013;339: 704–707. doi: 10.1126/science.1230406 23393265

39. Lin H ying, Liu Q, Li X, Yang J, Liu S, Huang Y, et al. Substantial contribution of genetic variation in the expression of transcription factors to phenotypic variation revealed by eRD-GWAS. Genome Biol. 2017. doi: 10.1186/s13059-017-1328-6

40. Swanson-Wagner R, Briskine R, Schaefer R, Hufford MB, Ross-Ibarra J, Myers CL, et al. Reshaping of the maize transcriptome by domestication. Proc Natl Acad Sci. 2012;109: 11878–11883. doi: 10.1073/pnas.1201961109 22753482

41. Rodgers-Melnick E, Vera DL, Bass HW, Buckler ES. Open chromatin reveals the functional maize genome. Proc Natl Acad Sci. 2016;113: E3177–84. doi: 10.1073/pnas.1525244113 27185945

42. Tenaillon MI, Charcosset A. A European perspective on maize history. Comptes Rendus—Biologies. 2011. pp. 221–228. doi: 10.1016/j.crvi.2010.12.015 21377617

43. Swarts K, Gutaker RM, Benz B, Blake M, Bukowski R, Holland J, et al. Genomic estimation of complex traits reveals ancient maize adaptation to temperate North America. Science. 2017;357: 512–515. doi: 10.1126/science.aam9425 28774930

44. Boden SA, Cavanagh C, Cullis BR, Ramm K, Greenwood J, Jean Finnegan E, et al. Ppd-1 is a key regulator of inflorescence architecture and paired spikelet development in wheat. Nat Plants. 2015;1: 14016. doi: 10.1038/nplants.2014.16 27246757

45. Ganal MW, Durstewitz G, Polley A, Bérard A, Buckler ES, Charcosset A, et al. A large maize (zea mays L.) SNP genotyping array: Development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome. PLoS One. 2011;6. doi: 10.1371/journal.pone.0028334

46. Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007. doi: 10.1086/521987

47. Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, et al. The B73 maize genome: Complexity, diversity, and dynamics. Science. 2009;326: 1112–1115. doi: 10.1126/science.1178534 19965430

48. Parent B, Turc O, Gibon Y, Stitt M, Tardieu F. Modelling temperature-compensated physiological rates, based on the co-ordination of responses to temperature of developmental processes. J Exp Bot. 2010;61: 2057–2069. doi: 10.1093/jxb/erq003 20194927

49. Cabrera-Bosquet L, Fournier C, Brichet N, Welcker C, Suard B, Tardieu F. High-throughput estimation of incident light, light interception and radiation-use efficiency of thousands of plants in a phenotyping platform. New Phytol. 2016;212: 269–281. doi: 10.1111/nph.14027 27258481

50. Rincent R, Laloë D, Nicolas S, Altmann T, Brunel D, Revilla P, et al. Maximizing the reliability of genomic selection by optimizing the calibration set of reference individuals: Comparison of methods in two diverse groups of maize inbreds (Zea mays L.). Genetics. 2012. doi: 10.1534/genetics.112.141473

51. Manoli A, Sturaro A, Trevisan S, Quaggiotti S, Nonis A. Evaluation of candidate reference genes for qPCR in maize. J Plant Physiol. 2012;169: 807–815. doi: 10.1016/j.jplph.2012.01.019 22459324

52. Lin Y, Zhang C, Lan H, Gao S, Liu H, Liu J, et al. Validation of potential reference genes for qPCR in maize across abiotic stresses, hormone treatments, and tissue types. PLoS One. 2014;9. doi: 10.1371/journal.pone.0095445

53. Winterhalter WE, Fedorka KM. Sex-specific variation in the emphasis, inducibility and timing of the post-mating immune response in Drosophila melanogaster. Proc R Soc B Biol Sci. 2009;276: 1109–1117. doi: 10.1098/rspb.2008.1559

54. Steger K, Wilhelm J, Konrad L, Stalf T, Greb R, Diemer T, et al. Both protamine-1 to protamine-2 mRNA ratio and Bcl2 mRNA content in testicular spermatids and ejaculated spermatozoa discriminate between fertile and infertile men. Hum Reprod. 2008;23: 11–16. doi: 10.1093/humrep/dem363 18003625

55. Laumanns IP, Fink L, Wilhelm J, Wolff JC, Mitnacht-Kraus R, Graef-Hoechst S, et al. The noncanonical WNT pathway is operative in idiopathic pulmonary arterial hypertension. Am J Respir Cell Mol Biol. 2009;40: 683–691. doi: 10.1165/rcmb.2008-0153OC 19029018

56. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/. 2013. doi: 10.1348/000712608X366867

57. Butler DG, Cullis BR, Gilmour AR, Gogel BJ. ASReml-R reference manual mixed models for S language environments. VSN International Ltd, Hemel Hempstead, HP1 1ES, UK www.vsni.co.uk. 2009.

58. Kruijer W, Boer MP, Malosetti M, Flood PJ, Engel B, Kooke R, et al. Marker-based estimation of heritability in immortal populations. Genetics. 2014;199: 379–398. doi: 10.1534/genetics.114.167916 25527288

59. Millet E, Welcker C, Kruijer W, Negro S, Nicolas S, Praud S, et al. Genome-wide analysis of yield in Europe: allelic effects as functions of drought and heat scenarios. Plant Physiol. 2016; pp.00621.2016. doi: 10.1104/pp.16.00621

60. Shipley B. Cause and Correlation in Biology. Cause and Correlation in Biology. 2016. doi: 10.1017/cbo9781139979573

61. McArdle J. Structural Equation Modeling: Applications in Ecological and Evolutionary Biology. Edited by Pugesek Bruce H, Adrian Tomer, and, Alexander von Eye. Cambridge and New York: Cambridge University Press. $75.00. xiii + 409 p; ill.; index. ISBN: 0–521–. The Quarterly Review of Biology. 2004. pp. 330–330. doi: 10.1086/425822

62. Lamb E, Shirtliffe S, May W. Structural equation modeling in the plant sciences: An example using yield components in oat. Can J Plant Sci. 2011;91: 603–619. doi: 10.4141/cjps2010-035

63. Rosseel Y. lavaan: An R Package for Structural Equation Modeling. J Stat Softw. 2012. doi: 10.18637/jss.v048.i02

64. Rincent R, Moreau L, Monod H, Kuhn E, Melchinger AE, Malvar RA, et al. Recovering power in association mapping panels with variable levels of linkage disequilibrium. Genetics. 2014;197: 375–387. doi: 10.1534/genetics.113.159731 24532779

65. Lippert C, Listgarten J, Liu Y, Kadie CM, Davidson RI, Heckerman D. FaST linear mixed models for genome-wide association studies. Nat Methods. 2011;8: 833–835. doi: 10.1038/nmeth.1681 21892150

66. Giraud H, Lehermeier C, Bauer E, Falque M, Segura V, Bauland C, et al. Linkage disequilibrium with linkage analysis of multiline crosses reveals different multiallelic QTL for hybrid performance in the flint and dent heterotic groups of maize. Genetics. 2014;198: 1717–1734. doi: 10.1534/genetics.114.169367 25271305

67. Clayton D. Differences between snpStats and snpMatrix. Differences. 2011.

68. Shin J-H, Blay S, Graham J, McNeney B. LDheatmap: An R Function for Graphical Display of Pairwise Linkage Disequilibria Between Single Nucleotide Polymorphisms. J Stat Softw. 2015;16. doi: 10.18637/jss.v016.c03

69. Foll M, Gaggiotti O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: A Bayesian perspective. Genetics. 2008;180: 977–993. doi: 10.1534/genetics.108.092221 18780740

70. MacQueen J. Some Methods for classification and Analysis of Multivariate Observations. 5th Berkeley Symposium on Mathematical Statistics and Probability 1967. 1967. citeulike-article-id:6083430

71. Earley KW, Haag JR, Pontes O, Opper K, Juehne T, Song K, et al. Gateway-compatible vectors for plant functional genomics and proteomics. Plant Journal. 2006. pp. 616–629. doi: 10.1111/j.1365-313X.2005.02617.x 16441352

72. Pradal C, Dufour-Kowalski S, Boudon F, Fournier C, Godin C. OpenAlea: A visual programming and component-based software platform for plant modelling. Funct Plant Biol. 2008. doi: 10.1071/FP08084

73. Fournier C, Andrieu B. ADEL-maize: An L-system based model for the integration of growth processes from the organ to the canopy. Application to regulation of morphogenesis by light availability. Agronomie. 1999. doi: 10.1051/agro:19990311


Článek vyšel v časopise

PLOS Genetics


2020 Číslo 7
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#