A MITE insertion abolishes the AP3-3 self-maintenance regulatory loop in apetalous flowers of Nigella damascena.

MADS-box transcription factors are important regulators of floral organ identity through their binding to specific motifs, termed CArG, in the promoter of their target genes. Petal initiation and development depend on class A and B genes, but MADS-box genes of the APETALA3 (AP3) clade are key regulators of this process. In the early diverging eudicot Nigella damascena, an apetalous [T] morph is characterized by the lack of expression of the NdAP3-3 gene, with its expression being petal-specific in the wild-type [P] morph. All [T] morph plants are homozygous for an NdAP3-3 allele with a Miniature Inverted-repeat Transposable Element (MITE) insertion in the second intron of the gene. Here, we investigated to which extent the MITE insertion impairs regulation of the NdAP3-3 gene. We found that expression of NdAP3-3 is initiated in the [T] morph, but the MITE insertion prevents its positive self-maintenance by affecting the correct splicing of the mRNA. We also found specific CArG features in the promoter of the NdAP3-3 genes with petal-specific expression. However, they are not sufficient to drive expression only in petals of transgenic Arabidopsis, highlighting the existence of Nigella-specific cis/trans-acting factors in regulating AP3 paralogs.

Transcriptome Analysis Reveals Putative Target Genes of APETALA3-3 During Early Floral Development in Nigella damascena L.

Even though petals are homoplastic structures, their identity consistently involves genes of the APETALA3 (AP3) lineage. However, the extent to which the networks downstream of AP3 are conserved in species with petals of different evolutionary origins is unknown. In Ranunculaceae, the specificity of the AP3-III lineage offers a great opportunity to identify the petal gene regulatory network in a comparative framework. Using a transcriptomic approach, we investigated putative target genes of the AP3-III ortholog NdAP3-3 in Nigella damascena at early developmental stages when petal identity is determined, and we compared our data with that from selected eudicot species. We generated a de novo reference transcriptome to carry out a differential gene expression analysis between the wild-type and mutant NdAP3-3 genotypes differing by the presence vs. absence of petals at early stages of floral development. Among the 1,620 genes that were significantly differentially expressed between the two genotypes, functional annotation suggested a large involvement of nuclear activities, including regulation of transcription, and enrichment in processes linked to cell proliferation. Comparing with Arabidopsis data, we found that highly conserved genes between the two species are enriched in homologs of direct targets of the AtAP3 protein. Integrating AP3-3 binding site data from another Ranunculaceae species, Aquilegia coerulea, allowed us to identify a set of 18 putative target genes that were conserved between the three species. Our results suggest that, despite the independent evolutionary origin of petals in core eudicots and Ranunculaceae, a small conserved set of genes determines petal identity and early development in these taxa.

Common gardens in teosintes reveal the establishment of a syndrome of adaptation to altitude.

In plants, local adaptation across species range is frequent. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Genome scans are popular to uncover outlier loci potentially involved in the genetic architecture of local adaptation, however links between outliers and phenotypic variation are rarely addressed. Here we focused on adaptation of teosinte populations along two elevation gradients in Mexico that display continuous environmental changes at a short geographical scale. We used two common gardens, and phenotyped 18 traits in 1664 plants from 11 populations of annual teosintes. In parallel, we genotyped these plants for 38 microsatellite markers as well as for 171 outlier single nucleotide polymorphisms (SNPs) that displayed excess of allele differentiation between pairs of lowland and highland populations and/or correlation with environmental variables. Our results revealed that phenotypic differentiation at 10 out of the 18 traits was driven by local selection. Trait covariation along the elevation gradient indicated that adaptation to altitude results from the assembly of multiple co-adapted traits into a complex syndrome: as elevation increases, plants flower earlier, produce less tillers, display lower stomata density and carry larger, longer and heavier grains. The proportion of outlier SNPs associating with phenotypic variation, however, largely depended on whether we considered a neutral structure with 5 genetic groups (73.7%) or 11 populations (13.5%), indicating that population stratification greatly affected our results. Finally, chromosomal inversions were enriched for both SNPs whose allele frequencies shifted along elevation as well as phenotypically-associated SNPs. Altogether, our results are consistent with the establishment of an altitudinal syndrome promoted by local selective forces in teosinte populations in spite of detectable gene flow. Because elevation mimics climate change through space, SNPs that we found underlying phenotypic variation at adaptive traits may be relevant for future maize breeding.

Testing the link between genome size and growth rate in maize.

Little is known about the factors driving within species Genome Size (GS) variation. GS may be shaped indirectly by natural selection on development and adaptative traits. Because GS variation is particularly pronounced in maize, we have sampled 83 maize inbred lines from three well described genetic groups adapted to contrasted climate conditions: inbreds of tropical origin, Flint inbreds grown in temperate climates, and Dent inbreds distributed in the Corn Belt. As a proxy for growth rate, we measured the Leaf Elongation Rate maximum during nighttime (LERmax) as well as GS in all inbred lines. In addition we combined available and new nucleotide polymorphism data at 29,090 sites to characterize the genetic structure of our panel. We found significant variation for both LERmax and GS among groups defined by our genetic structuring. Tropicals displayed larger GS than Flints while Dents exhibited intermediate values. LERmax followed the opposite trend with greater growth rate in Flints than in Tropicals. In other words, LERmax and GS exhibited a significantly negative correlation (r = - 0.27). However, this correlation was driven by among-group variation rather than within-group variation-it was no longer significant after controlling for structure and kinship among inbreds. Our results indicate that selection on GS may have accompanied ancient maize diffusion from its center of origin, with large DNA content excluded from temperate areas. Whether GS has been targeted by more intense selection during modern breeding within groups remains an open question.

Molecular evolution accompanying functional divergence of duplicated genes along the plant starch biosynthesis pathway.

BACKGROUND: Starch is the main source of carbon storage in the Archaeplastida. The starch biosynthesis pathway (sbp) emerged from cytosolic glycogen metabolism shortly after plastid endosymbiosis and was redirected to the plastid stroma during the green lineage divergence. The SBP is a complex network of genes, most of which are members of large multigene families. While some gene duplications occurred in the Archaeplastida ancestor, most were generated during the sbp redirection process, and the remaining few paralogs were generated through compartmentalization or tissue specialization during the evolution of the land plants. In the present study, we tested models of duplicated gene evolution in order to understand the evolutionary forces that have led to the development of SBP in angiosperms. We combined phylogenetic analyses and tests on the rates of evolution along branches emerging from major duplication events in six gene families encoding sbp enzymes. RESULTS: We found evidence of positive selection along branches following cytosolic or plastidial specialization in two starch phosphorylases and identified numerous residues that exhibited changes in volume, polarity or charge. Starch synthases, branching and debranching enzymes functional specializations were also accompanied by accelerated evolution. However, none of the sites targeted by selection corresponded to known functional domains, catalytic or regulatory. Interestingly, among the 13 duplications tested, 7 exhibited evidence of positive selection in both branches emerging from the duplication, 2 in only one branch, and 4 in none of the branches. CONCLUSIONS: The majority of duplications were followed by accelerated evolution targeting specific residues along both branches. This pattern was consistent with the optimization of the two sub-functions originally fulfilled by the ancestral gene before duplication. Our results thereby provide strong support to the so-called "Escape from Adaptive Conflict" (EAC) model. Because none of the residues targeted by selection occurred in characterized functional domains, we propose that enzyme specialization has occurred through subtle changes in affinity, activity or interaction with other enzymes in complex formation, while the basic function defined by the catalytic domain has been maintained.

An APETALA3 homolog controls both petal identity and floral meristem patterning in Nigella damascena L. (Ranunculaceae).

Flower architecture mutants provide a unique opportunity to address the genetic origin of flower diversity. Here we study a naturally occurring floral dimorphism in Nigella damascena (Ranunculaceae), involving replacement of the petals by numerous sepal-like and chimeric sepal/stamen organs. We performed a comparative study of floral morphology and floral development, and characterized the expression of APETALA3 and PISTILLATA homologs in both morphs. Segregation analyses and gene silencing were used to determine the involvement of an APETALA3 paralog (NdAP3-3) in the floral dimorphism. We demonstrate that the complex floral dimorphism is controlled by a single locus, which perfectly co-segregates with the NdAP3-3 gene. This gene is not expressed in the apetalous morph and exhibits a particular expression dynamic during early floral development in the petalous morph. NdAP3-3 silencing in petalous plants perfectly phenocopies the apetalous morph. Our results show that NdAP3-3 is fully responsible for the complex N. damascena floral dimorphism, suggesting that it plays a role not only in petal identity but also in meristem patterning, possibly through regulation of perianth organ number and the perianth/stamen boundary.

Accelerated evolution and coevolution drove the evolutionary history of AGPase sub-units during angiosperm radiation.

BACKGROUND AND AIMS: ADP-glucose pyrophosphorylase (AGPase) is a key enzyme of starch biosynthesis. In the green plant lineage, it is composed of two large (LSU) and two small (SSU) sub-units encoded by paralogous genes, as a consequence of several rounds of duplication. First, our aim was to detect specific patterns of molecular evolution following duplication events and the divergence between monocotyledons and dicotyledons. Secondly, we investigated coevolution between amino acids both within and between sub-units. METHODS: A phylogeny of each AGPase sub-unit was built using all gymnosperm and angiosperm sequences available in databases. Accelerated evolution along specific branches was tested using the ratio of the non-synonymous to the synonymous substitution rate. Coevolution between amino acids was investigated taking into account compensatory changes between co-substitutions. KEY RESULTS: We showed that SSU paralogues evolved under high functional constraints during angiosperm radiation, with a significant level of coevolution between amino acids that participate in SSU major functions. In contrast, in the LSU paralogues, we identified residues under positive selection (1) following the first LSU duplication that gave rise to two paralogues mainly expressed in angiosperm source and sink tissues, respectively; and (2) following the emergence of grass-specific paralogues expressed in the endosperm. Finally, we found coevolution between residues that belong to the interaction domains of both sub-units. CONCLUSIONS: Our results support the view that coevolution among amino acid residues, especially those lying in the interaction domain of each sub-unit, played an important role in AGPase evolution. First, within SSU, coevolution allowed compensating mutations in a highly constrained context. Secondly, the LSU paralogues probably acquired tissue-specific expression and regulatory properties via the coevolution between sub-unit interacting domains. Finally, the pattern we observed during LSU evolution is consistent with repeated sub-functionalization under 'Escape from Adaptive Conflict', a model rarely illustrated in the literature.

Epistatic interactions between Opaque2 transcriptional activator and its target gene CyPPDK1 control kernel trait variation in maize.

Association genetics is a powerful method to track gene polymorphisms responsible for phenotypic variation, since it takes advantage of existing collections and historical recombination to study the correlation between large genetic diversity and phenotypic variation. We used a collection of 375 maize (Zea mays ssp. mays) inbred lines representative of tropical, American, and European diversity, previously characterized for genome-wide neutral markers and population structure, to investigate the roles of two functionally related candidate genes, Opaque2 and CyPPDK1, on kernel quality traits. Opaque2 encodes a basic leucine zipper transcriptional activator specifically expressed during endosperm development that controls the transcription of many target genes, including CyPPDK1, which encodes a cytosolic pyruvate orthophosphate dikinase. Using statistical models that correct for population structure and individual kinship, Opaque2 polymorphism was found to be strongly associated with variation of the essential amino acid lysine. This effect could be due to the direct role of Opaque2 on either zein transcription, zeins being major storage proteins devoid of lysine, or lysine degradation through the activation of lysine ketoglutarate reductase. Moreover, we found that a polymorphism in the Opaque2 coding sequence and several polymorphisms in the CyPPDK1 promoter nonadditively interact to modify both lysine content and the protein-versus-starch balance, thus revealing the role in quantitative variation in plants of epistatic interactions between a transcriptional activator and one of its target genes.

Whole genome surveys of rice, maize and sorghum reveal multiple horizontal transfers of the LTR-retrotransposon Route66 in Poaceae.

BACKGROUND: Horizontal transfers (HTs) refer to the transmission of genetic material between phylogenetically distant species. Although most of the cases of HTs described so far concern genes, there is increasing evidence that some involve transposable elements (TEs) in Eukaryotes. The availability of the full genome sequence of two cereal species, (i.e. rice and Sorghum), as well as the partial genome sequence of maize, provides the opportunity to carry out genome-wide searches for TE-HTs in Poaceae. RESULTS: We have identified an LTR-retrotransposon, that we named Route66, with more than 95% sequence identity between rice and Sorghum. Using a combination of in silico and molecular approaches, we are able to present a substantial phylogenetic evidence that Route66 has been transferred horizontally between Panicoideae and several species of the genus Oryza. In addition, we show that it has remained active after these transfers. CONCLUSION: This study constitutes a new case of HTs for an LTR-retrotransposon and we strongly believe that this mechanism could play a major role in the life cycle of transposable elements. We therefore propose to integrate classe I elements into the previous model of transposable element evolution through horizontal transfers.

Patterns of molecular evolution associated with two selective sweeps in the Tb1-Dwarf8 region in maize.

We focused on a region encompassing a major maize domestication locus, Tb1, and a locus involved in the flowering time variation, Dwarf8 (D8), to investigate the consequences of two closely linked selective sweeps on nucleotide variation and gain some insights into maize geographical diffusion, through climate adaptation. First, we physically mapped D8 at approximately 300 kb 3' of Tb1. Second, we analyzed patterns of nucleotide variation at Tb1, D8, and seven short regions (400-700 bp) located in the Tb1-D8 region sequenced on a 40 maize inbred lines panel encompassing early-flowering temperate and late-flowering tropical lines. The pattern of polymorphism along the region is characterized by two valleys of depleted polymorphism while the region in between exhibits an appreciable amount of diversity. Our results reveal that a region approximately 100 kb upstream of the D8 gene exhibits hallmarks of divergent selection between temperate and tropical lines and is likely closer than the D8 gene to the target of selection for climate adaptation. Selection in the tropical lines appears more recent than in the temperate lines, suggesting an initial domestication of early-flowering maize. Simulation results indicate that the polymorphism pattern is consistent with two interfering selective sweeps at Tb1 and D8.

Key impact of Vgt1 on flowering time adaptation in maize: evidence from association mapping and ecogeographical information.

An association study conducted on 375 maize inbred lines indicates a strong relationship between Vgt1 polymorphisms and flowering time, extending former quantitative trait loci (QTL) mapping results. Analysis of allele frequencies in a landrace collection supports a key role of Vgt1 in maize altilatitudinal adaptation.

Maize adaptation to temperate climate: relationship between population structure and polymorphism in the Dwarf8 gene.

To investigate the genetic basis of maize adaptation to temperate climate, collections of 375 inbred lines and 275 landraces, representative of American and European diversity, were evaluated for flowering time under short- and long-day conditions. The inbred line collection was genotyped for 55 genomewide simple sequence repeat (SSR) markers. Comparison of inbred line population structure with that of landraces, as determined with 24 SSR loci, underlined strong effects of both historical and modern selection on population structure and a clear relationship with geographical origins. The late tropical groups and the early "Northern Flint" group from the northern United States and northern Europe exhibited different flowering times. Both collections were genotyped for a 6-bp insertion/deletion in the Dwarf8 (D8idp) gene, previously reported to be potentially involved in flowering time variation in a 102 American inbred panel. Among-group D8idp differentiation was much higher than that for any SSR marker, suggesting diversifying selection. Correcting for population structure, D8idp was associated with flowering time under long-day conditions, the deletion allele showing an average earlier flowering of 29 degree days for inbreds and 145 degree days for landraces. Additionally, the deletion allele occurred at a high frequency (>80%) in Northern Flint while being almost absent (<5%) in tropical materials. Altogether, these results indicate that Dwarf8 could be involved in maize climatic adaptation through diversifying selection for flowering time.

Molecular evolution of the Opaque-2 gene in Zea mays L.

The Opaque-2 gene (O2) in maize encodes a transcriptional activator that controls the expression of various genes during kernel development, particularly some of the most abundant endosperm storage protein genes. Compared to its wild relative teosinte, maize has bigger and heavier kernels, with an increased proportion of starch and an altered distribution of the various storage protein categories. The molecular evolution of the O2 gene was investigated in connection with its possible involvement in the domestication process. Most of the coding sequence and parts of introns, 5'UTR, and 3' noncoding regions were sequenced in a set of cultivated and teosinte accessions. One hundred six polymorphic sites (5.4%) and 72 insertions/deletions, located mostly in noncoding regions, were found. Molecular diversity was quite high (pi = 0.0138, theta = 0.0167) compared to that of other transcription factors in maize. The synonymous and nonsynonymous diversity patterns along the coding sequence suggested that different regions are submitted to different functional constraints. Such an evolution would probably be favored by the observed rapid decay of linkage disequilibrium with distance. Cultivated accessions retained about 70% of the diversity observed in teosintes. Purifying selection was detected in both maize and teosintes. No conclusive evidence was obtained for a role of the O2 gene in the domestication process.

QTLs for enzyme activities and soluble carbohydrates involved in starch accumulation during grain filling in maize.

ADPglucose, the essential substrate for starch synthesis, is synthesized in maize by a pathway involving at least invertases, sucrose synthase, and ADPglucose pyrophosphorylase, as shown by the starch-deficient mutants, mn1, sh1, and bt2 or sh2, respectively. To improve understanding of the relationship between early grain-filling traits and carbohydrate composition in mature grain, QTLs linked to soluble invertase, sucrose synthase, and ADPglucose pyrophosphorylase activities and to starch, sucrose, fructose, and glucose concentrations were investigated. In order to take into account the specific time-course of each enzyme activity during grain filling, sampling was carried out at three periods (15, 25, and 35 d after pollination) on 100 lines from a recombinant inbred family, grown in the field. The MQTL method associated with QTL interaction analysis revealed numerous QTLs for all traits, but only one QTL was consistently observed at the three sampling periods. Some chromosome zones were heavily labelled, forming clusters of QTLs. Numerous possible candidate genes of the starch synthetic pathway co-located with QTLs. Four QTLs were found close to the locus Sh1 (bin 9.01) coding for the sucrose synthase. In order to confirm the importance of this locus, the CAPS polymorphism of the Sh1 gene was analysed in 45 genetically unrelated maize lines from various geographical origins. The DNA polymorphism was significantly associated with phenotypic traits related to grain filling (starch and amylose content, grain matter, and ADPglucose pyrophosphorylase activity at 35 DAP). Thus, the Sh1 locus could provide a physiologically pertinent marker for maize selection.

Nucleotide diversity of the ZmPox3 maize peroxidase gene: relationships between a MITE insertion in exon 2 and variation in forage maize digestibility.

BACKGROUND: Polymorphisms were investigated within the ZmPox3 maize peroxidase gene, possibly involved in lignin biosynthesis because of its colocalization with a cluster of QTL related to lignin content and cell wall digestibility. The purpose of this study was to identify, on the basis of 37 maize lines chosen for their varying degrees of cell wall digestibility and representative of temperate regions germplasm, ZmPox3 haplotypes or individual polymorphisms possibly associated with digestibility. RESULTS: Numerous haplotypes with high diversity were identified. Frequency of nucleotide changes was high with on average one SNP every 57 bp. Nucleotide diversity was not equally distributed among site categories: the estimated pi was on average eight times higher for silent sites than for non-synonymous sites. Numerous sites were in linkage disequilibrium that decayed with increasing physical distance. A zmPox3 mutant allele, carrying an insertion of a transposable element in the second exon, was found in lines derived from the early flint inbred line, F7. This element possesses many structural features of miniature inverted-repeat transposable elements (MITE). The mutant allele encodes a truncated protein lacking important functional sites. An ANOVA performed with a subset of 31 maize lines indicated that the transposable element was significantly associated with cell wall digestibility. This association was confirmed using an additional set of 25 flint lines related to F7. Moreover, RT-PCR experiments revealed a decreased amount of corresponding mRNA in plants with the MITE insertion. CONCLUSION: These results showed that ZmPox3 could possibly be involved in monolignol polymerisation, and that a deficiency in ZmPox3 peroxidase activity seemingly has a negative effect on cell wall digestibility. Also, genetic diversity analyses of ZmPox3 indicated that this peroxidase could be a relevant target for grass digestibility improvement using specific allele introgressions.

GENETIC DRIFT AND FOUNDER EFFECT IN NATIVE VERSUS INTRODUCED POPULATIONS OF AN INVADING PLANT, LYTHRUM SALICARIA (LYTHRACEAE).

There are few convincing examples of genetic drift at loci under selection in natural populations. The plant sexual polymorphism tristyly provides an opportunity to investigate genetic drift because stochastic processes interacting with frequency-dependent selection give rise to a diagnostic pattern of morph-frequency variation. A previous study of 102 Ontario populations of the introduced tristylous wetland herb Lythrum salicaria provided evidence for the role of stochastic processes during colonization. However, whether stochastic effects are greater in these recently introduced populations compared to native Eurasian populations remains unclear. The propensity of this species to invade disturbed habitats suggests that episodes of colonization and periods of small population size must also occur in the native range. A survey of 102 populations in southwestern France indicated reduced stochastic effects in native populations. Populations exhibited significantly lower morph loss than in Ontario (5% vs. 23%) and significantly higher values of morph evenness. The greater incidence of trimorphism in French populations was not associated with larger population sizes; populations were significantly smaller than those in Ontario (means: 266 vs. 487). Morph evenness was positively correlated with population size among French but not Ontario populations, providing further evidence of nonequilibrium conditions in introduced compared to native populations. The incidence of trimorphism was unexpectedly high in small native populations (N ≤ 25; 22 of 27 populations trimorphic). Computer simulations indicated that levels of gene flow on the order of m ≥ 0.05 can account for the maintenance of tristyly in small populations. The high connectivity of populations within the agricultural landscape typical of southwestern France may facilitate levels of gene flow sufficient to maintain trimorphism in small populations.