Intronic regulatory elements determine the divergent expression patterns of AGAMOUS-LIKE6 subfamily members in Arabidopsis

The MADS-box genes SOC1 and AGL24 antagonize XAL2 functions in Arabidopsis thaliana root development

MADS-domain transcription factors play pivotal roles in numerous developmental processes in Arabidopsis thaliana. While their involvement in flowering transition and floral development has been extensively examined, their functions in root development remain relatively unexplored. Here, we explored the function and genetic interaction of three MADS-box genes (XAL2, SOC1 and AGL24) in primary root development. By analyzing loss-of-function and overexpression lines, we found that SOC1 and AGL24, both critical components in flowering transition, redundantly act as repressors of primary root growth as the loss of function of either SOC1 or AGL24 partially recovers the primary root growth, meristem cell number, cell production rate, and the length of fully elongated cells of the short-root mutant xal2-2. Furthermore, we observed that the simultaneous overexpression of AGL24 and SOC1 leads to short-root phenotypes, affecting meristem cell number and fully elongated cell size, whereas SOC1 overexpression is sufficient to affect columella stem cell differentiation. Additionally, qPCR analyses revealed that these genes exhibit distinct modes of transcriptional regulation in roots compared to what has been previously reported for aerial tissues. We identified 100 differentially expressed genes in xal2-2 roots by RNA-seq. Moreover, our findings revealed that the expression of certain genes involved in cell differentiation, as well as stress responses, which are either upregulated or downregulated in the xal2-2 mutant, reverted to WT levels in the absence of SOC1 or AGL24.

A transcriptional atlas identifies key regulators and networks for the development of spike tissues in barley

Grain number and size determine grain yield in crops and are closely associated with spikelet fertility and grain filling in barley (Hordeum vulgare). Abortion of spikelet primordia within individual barley spikes causes a 30%-50% loss in the potential number of grains during development from the awn primordium stage to the tipping stage, after that grain filling is the primary factor regulating grain size. To identify transcriptional signatures associated with spike development, we use a six-rowed barley cultivar (Morex) to develop a spatiotemporal transcriptome atlas containing 255 samples covering 17 stages and 5 positions along the spike. We identify several fundamental regulatory networks, in addition to key regulators of spike development and morphology. Specifically, we show HvGELP96, encoding a GDSL domain-containing protein, as a regulator of spikelet fertility and grain number. Our transcriptional atlas offers a powerful resource to answer fundamental questions in spikelet development and degeneration in barley.

Genome-wide identification and expression analysis of PYL family genes and functional characterization of GhPYL8D2 under drought stress in Gossypium hirsutum

Cotton is a crucial economic crop, serving as a natural fiber source for the textile industry. However, drought stress poses a significant threat to cotton fiber quality and productivity worldwide. Pyrabactin Resistance 1-Like (PYL) proteins, as abscisic acid (ABA) receptors, play a crucial role in adverse stress responses, but knowledge about the PYLs in cotton remains limited. In our study, we identified 40 GhPYL genes in Gossypium hirsutum through a genome-wide analysis of the cotton genome database. Our analysis revealed that the PYL family formed three distinct subfamilies with typical family characteristics in G. hirsutum. Additionally, through quantitative expression analysis, including transcriptome dataset and qRT-PCR, we found that all GhPYLs were expressed in all tissues of G. hirsutum, and all GhPYLs were differentially expressed under drought stress. Among them, GhPYL4A1, GhPY5D1, GhPY8D2, and a member of the type 2C protein phosphatases clade A family in Gossypium hirsutum (GhPP2CA), GhHAI2D, showed significant differences in expression levels within 12 h after stress treatment. Our protein interaction analysis and BiFC demonstrated the complex regulatory network between GhPYL family proteins and GhPP2CA proteins. We also found that there is an interaction between GhPYL8D2 and GhHAI2D, and through drought treatment of transgenic cotton, we found that GhPYL8D2 played a vital role in the response of G. hirsutum to drought through stomatal control via co-regulation with GhHAI2D. Our findings provide useful insights into the regulation of GhPYL family genes that occur in response to abiotic stresses in cotton.

Functional Divergence Analysis of AGL6 Genes in Prunus mume

The AGAMOUS-LIKE6 (AGL6) lineage is an important clade of MADS-box transcription factors that play essential roles in floral organ development. The genome of Prunus mume contains two homoeologous AGL6 genes that are replicated as gene fragments. In this study, two AGL6 homologs, PmAGL6-1 and PmAGL6-2, were cloned from P. mume and then functionally characterized. Sequence alignment and phylogenetic analyses grouped both genes into the AGL6 lineage. The expression patterns and protein-protein interaction patterns showed significant differences between the two genes. However, the ectopic expression of the two genes in Arabidopsis thaliana resulted in similar phenotypes, including the promotion of flowering, alteration of floral organ structure, participation in the formation of the floral meristem and promotion of pod bending. Therefore, gene duplication has led to some functional divergence of PmAGL6-1 and PmAGL6-2 but their functions are similar. We thus speculated that AGL6 genes play a crucial role in flower development in P. mume.

Dynamic genome evolution in a model fern

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Regulatory network for FOREVER YOUNG FLOWER-like genes in regulating Arabidopsis flower senescence and abscission

FOREVER YOUNG FLOWER (FYF) has been reported to play an important role in regulating flower senescence/abscission. Here, we functionally analyzed five Arabidopsis FYF-like genes, two in the FYF subgroup (FYL1/AGL71 and FYL2/AGL72) and three in the SOC1 subgroup (SOC1/AGL20, AGL19, and AGL14/XAL2), and showed their involvement in the regulation of flower senescence and/or abscission. We demonstrated that in FYF subgroup, FYF has both functions in suppressing flower senescence and abscission, FYL1 only suppresses flower abscission and FYL2 has been converted as an activator to promote flower senescence. In SOC1 subgroup, AGL19/AGL14/SOC1 have only one function in suppressing flower senescence. We also found that FYF-like proteins can form heterotetrameric complexes with different combinations of A/E functional proteins (such as AGL6 and SEP1) and AGL15/18-like proteins to perform their functions. These findings greatly expand the current knowledge behind the multifunctional evolution of FYF-like genes and uncover their regulatory network in plants.

Identification of key regulatory genes involved in the sporophyte and gametophyte development in Ginkgo biloba ovules revealed by in situ expression analyses

PREMISE

In Arabidopsis thaliana, the role of the most important key genes that regulate ovule development is widely known. In nonmodel species, and especially in gymnosperms, the ovule developmental processes are still quite obscure. In this study, we describe the putative roles of Ginkgo biloba orthologs of regulatory genes during ovule development. Specifically, we studied AGAMOUS (AG), AGAMOUS-like 6 (AGL6), AINTEGUMENTA (ANT), BELL1 (BEL1), Class III HD-Zip, and YABBY Ginkgo genes.

METHODS

We analyzed their expression domains through in situ hybridizations on two stages of ovule development: the very early stage that corresponds to the ovule primordium, still within wintering buds, and the late stage at pollination time.

RESULTS

GBM5 (Ginkgo ortholog of AG), GbMADS8 (ortholog of AGL6) and GbC3HDZ1-2-3 were expressed in both the stages of ovule development, while GbMADS1, GbAGL6-like genes (orthologs of AGL6), GbBEL1-2 and YABBY Ginkgo orthologs (GbiYAB1B and GbiYABC) seem mostly involved at pollination time. GbANTL1 was not expressed in the studied stages and was different from GbANTL2 and GbBEL1, which seem to be involved at both stages of ovule development. In Ginkgo, the investigated genes display patterns of expression only partially comparable to those of other studied seed plants.

CONCLUSIONS

The expression of most of these regulatory genes in the female gametophyte region at pollination time leads to suggest a communication between the sporophytic maternal tissue and the developing female gametophyte, as demonstrated for well-studied model angiosperms.

Genetic and Molecular Mechanisms Conferring Heat Stress Tolerance in Tomato Plants

Climate change is a major threat to global food security. Changes in climate can directly impact food systems by reducing the production and genetic diversity of crops and their wild relatives, thereby restricting future options for breeding improved varieties and reducing the ability to adapt crops to future challenges. The global surface temperature is predicted to rise by an average of 0.3°C during the next decade, and the Paris Agreement (Paris Climate Accords) aims to limit global warming to below an average of 2°C, preferably to 1.5°C compared to pre-industrial levels. Even if the goal of the Paris Agreement can be met, the predicted rise in temperatures will increase the likelihood of extreme weather events, including heatwaves, making heat stress (HS) a major global abiotic stress factor for many crops. HS can have adverse effects on plant morphology, physiology, and biochemistry during all stages of vegetative and reproductive development. In fruiting vegetables, even moderate HS reduces fruit set and yields, and high temperatures may result in poor fruit quality. In this review, we emphasize the effects of abiotic stress, especially at high temperatures, on crop plants, such as tomatoes, touching upon key processes determining plant growth and yield. Specifically, we investigated the molecular mechanisms involved in HS tolerance and the challenges of developing heat-tolerant tomato varieties. Finally, we discuss a strategy for effectively improving the heat tolerance of vegetable crops.

Transcriptional regulation in plants: Using omics data to crack the cis-regulatory code

Innovative omics technologies, advanced bioinformatics, and machine learning methods are rapidly becoming integral tools for plant functional genomics, with tremendous recent advances made in this field. In transcriptional regulation, an initial lag in the accumulation of plant omics data relative to that of animals stimulated the development of computational methods capable of extracting maximum information from the available data sets. Recent comprehensive studies of transcription factor-binding profiles in Arabidopsis and maize and the accumulation of uniformly processed omics data in public databases have brought plant biologists into the big leagues, with many cutting-edge methods available. Here, we summarize the state-of-the-art bioinformatics approaches used to predict or infer the cis-regulatory code behind transcriptional gene regulation, focusing on their plant research applications.

Challenges and Approaches to Crop Improvement Through C3-to-C4 Engineering

With a rapidly growing world population and dwindling natural resources, we are now facing the enormous challenge of increasing crop yields while simultaneously improving the efficiency of resource utilization. Introduction of C4 photosynthesis into C3 crops is widely accepted as a key strategy to meet this challenge because C4 plants are more efficient than C3 plants in photosynthesis and resource usage, particularly in hot climates, where the potential for productivity is high. Lending support to the feasibility of this C3-to-C4 engineering, evidence indicates that C4 photosynthesis has evolved from C3 photosynthesis in multiple lineages. Nevertheless, C3-to-C4 engineering is not an easy task, as several features essential to C4 photosynthesis must be introduced into C3 plants. One such feature is the spatial separation of the two phases of photosynthesis (CO₂ fixation and carbohydrate synthesis) into the mesophyll and bundle sheath cells, respectively. Another feature is the Kranz anatomy, characterized by a close association between the mesophyll and bundle sheath (BS) cells (1:1 ratio). These anatomical features, along with a C4-specific carbon fixation enzyme (PEPC), form a CO₂-concentration mechanism that ensures a high photosynthetic efficiency. Much effort has been taken in the past to introduce the C4 mechanism into C3 plants, but none of these attempts has met with success, which is in my opinion due to a lack of system-level understanding and manipulation of the C3 and C4 pathways. As a prerequisite for the C3-to-C4 engineering, I propose that not only the mechanisms that control the Kranz anatomy and cell-type-specific expression in C3 and C4 plants must be elucidated, but also a good understanding of the gene regulatory network underlying C3 and C4 photosynthesis must be achieved. In this review, I first describe the past and current efforts to increase photosynthetic efficiency in C3 plants and their limitations; I then discuss a systems approach to tackling down this challenge, some practical issues, and recent technical innovations that would help us to solve these problems.

Genome-Wide Identification and Analysis of the MADS-Box Gene Family in American Beautyberry (Callicarpa americana)

The MADS-box gene family encodes a number of transcription factors that play key roles in various plant growth and development processes from response to environmental cues to cell differentiation and organ identity, especially the floral organogenesis, as in the prominent ABCDE model of flower development. Recently, the genome of American beautyberry (Callicarpa americana) has been sequenced. It is a shrub native to the southern region of United States with edible purple-colored berries; it is a member of the Lamiaceae family, a family of medical and agricultural importance. Seventy-eight MADS-box genes were identified from 17 chromosomes of the C. americana assembled genome. Peptide sequences blast and analysis of phylogenetic relationships with MADS-box genes of Sesame indicum, Solanum lycopersicum, Arabidopsis thaliana, and Amborella trichopoda were performed. Genes were separated into 32 type I and 46 type II MADS-box genes. C. americana MADS-box genes were clustered into four groups: MIKC^C, MIKC*, Mα-type, and Mγ-type, while the Mβ-type group was absent. Analysis of the gene structure revealed that from 1 to 15 exons exist in C. americana MADS-box genes. The number of exons in type II MADS-box genes (5–15) greatly exceeded the number in type I genes (1–9). The motif distribution analysis of the two types of MADS-box genes showed that type II MADS-box genes contained more motifs than type I genes. These results suggested that C. americana MADS-box genes type II had more complex structures and might have more diverse functions. The role of MIKC-type MADS-box genes in flower and fruit development was highlighted when the expression profile was analyzed in different organs transcriptomes. This study is the first genome-wide analysis of the C. americana MADS-box gene family, and the results will further support any functional and evolutionary studies of C. americana MADS-box genes and serve as a reference for related studies of other plants in the medically important Lamiaceae family.

Tomato agamous-like6 parthenocarpy is facilitated by ovule integument reprogramming involving the growth regulator KLUH

Fruit set is established during and soon after fertilization of the ovules inside the quiescent ovary, but the signaling pathways involved remain obscure. The tomato (Solanum lycopersicum) CRISPR loss-of-function mutant of the transcription factor gene AGAMOUS-like6 (SlAGL6; slagl6CR-sg1) is capable of fertilization-independent setting of normal, yet seedless (parthenocarpic), fruit. To gain insight into the mechanism of fleshy fruit set, in this study, we investigated how slagl6CR-sg1 uncouples fruit set from fertilization. We found that mutant ovules were enlarged due to integument over-proliferation and failed to differentiate an endothelium, the integument's innermost layer, upon maturation. A causal relationship between slagl6 loss-of-function and these abnormal phenotypes is inferred from the observation that SlAGL6 is predominantly expressed in the immature ovule integument, and upon ovule maturation, its expression shifts to the endothelium. The transcriptome of unfertilized mutant ovules profoundly differs from that of wild-type and exhibits substantial overlap with the transcriptomes of fertilized ovules sporophytic tissues. One prominent upregulated gene was the fertilization-induced cytochrome P450 cell proliferation regulator SlKLUH. Indeed, ectopic overexpression of SlKLUH stimulated both integument growth in unfertilized ovules and parthenocarpy, suggesting that its suppression by SlAGL6 is paramount for preventing fertilization-independent fruit set. Taken together, our study informs on the transcriptional programs that are regulated by SlAGL6 and demonstrates that it acts from within the ovule integument to inhibit ovary growth beyond anthesis. That by suppressing components of the fertilization-induced ovule reprogramming underlying fruit set.

Potentials and Challenges of Genomics for Breeding Cannabis Cultivars

Cannabis (Cannabis sativa L.) is an influential yet controversial agricultural plant with a very long and prominent history of recreational, medicinal, and industrial usages. Given the importance of this species, we deepened some of the main challenges-along with potential solutions-behind the breeding of new cannabis cultivars. One of the main issues that should be fixed before starting new breeding programs is the uncertain taxonomic classification of the two main taxa (e.g., indica and sativa) of the Cannabis genus. We tried therefore to examine this topic from a molecular perspective through the use of DNA barcoding. Our findings seem to support a unique species system (C. sativa) based on two subspecies: C. sativa subsp. sativa and C. sativa subsp. indica. The second key issue in a breeding program is related to the dioecy behavior of this species and to the comprehension of those molecular mechanisms underlying flower development, the main cannabis product. Given the role of MADS box genes in flower identity, we analyzed and reorganized all the genomic and transcriptomic data available for homeotic genes, trying to decipher the applicability of the ABCDE model in Cannabis. Finally, reviewing the limits of the conventional breeding methods traditionally applied for developing new varieties, we proposed a new breeding scheme for the constitution of F₁ hybrids, without ignoring the indisputable contribution offered by genomics. In this sense, in parallel, we resumed the main advances in the genomic field of this species and, ascertained the lack of a robust set of SNP markers, provided a discriminant and polymorphic panel of SSR markers as a valuable tool for future marker assisted breeding programs.

Global Quantitative Mapping of Enhancers in Rice by STARR-seq

Enhancers activate transcription in a distance-, orientation-, and position-independent manner, which makes them difficult to be identified. Self-transcribing active regulatory region sequencing (STARR-seq) measures the enhancer activity of millions of DNA fragments in parallel. Here we used STARR-seq to generate a quantitative global map of rice enhancers. Most enhancers were mapped within genes, especially at the 5′ untranslated regions (5′UTR) and in coding sequences. Enhancers were also frequently mapped proximal to silent and lowly-expressed genes in transposable element (TE)-rich regions. Analysis of the epigenetic features of enhancers at their endogenous loci revealed that most enhancers do not co-localize with DNase I hypersensitive sites (DHSs) and lack the enhancer mark of histone modification H3K4me1. Clustering analysis of enhancers according to their epigenetic marks revealed that about 40% of identified enhancers carried one or more epigenetic marks. Repressive H3K27me3 was frequently enriched with positive marks, H3K4me3 and/or H3K27ac, which together label enhancers. Intergenic enhancers were also predicted based on the location of DHS regions relative to genes, which overlap poorly with STARR-seq enhancers. In summary, we quantitatively identified enhancers by functional analysis in the genome of rice, an important model plant. This work provides a valuable resource for further mechanistic studies in different biological contexts.

Origination and selection of ABCDE and AGL6 subfamily MADS-box genes in gymnosperms and angiosperms

BACKGROUND

The morphological diversity of flower organs is closely related to functional divergence within the MADS-box gene family. Bryophytes and seedless vascular plants have MADS-box genes but do not have ABCDE or AGAMOUS-LIKE6 (AGL6) genes. ABCDE and AGL6 genes belong to the subgroup of MADS-box genes. Previous works suggest that the B gene was the first ABCDE and AGL6 genes to emerge in plant but there are no mentions about the probable origin time of ACDE and AGL6 genes. Here, we collected ABCDE and AGL6 gene 381 protein sequences and 361 coding sequences from gymnosperms and angiosperms and reconstructed a complete Bayesian phylogeny of these genes. In this study, we want to clarify the probable origin time of ABCDE and AGL6 genes is a great help for understanding the role of the formation of the flower, which can decipher the forming order of MADS-box genes in the future.

RESULTS

These genes appeared to have been under purifying selection and their evolutionary rates are not significantly different from each other. Using the Bayesian evolutionary analysis by sampling trees (BEAST) tool, we estimated that: the mutation rate of the ABCDE and AGL6 genes was 2.617 × 10-3 substitutions/site/million years, and that B genes originated 339 million years ago (MYA), CD genes originated 322 MYA, and A genes shared the most recent common ancestor with E/AGL6 296 MYA, respectively.

CONCLUSIONS

The phylogeny of ABCDE and AGL6 genes subfamilies differed. The APETALA1 (AP1 or A gene) subfamily clustered into one group. The APETALA3/PISTILLATA (AP3/PI or B genes) subfamily clustered into two groups: the AP3 and PI clades. The AGAMOUS/SHATTERPROOF/SEEDSTICK (AG/SHP/STK or CD genes) subfamily clustered into a single group. The SEPALLATA (SEP or E gene) subfamily in angiosperms clustered into two groups: the SEP1/2/4 and SEP3 clades. The AGL6 subfamily clustered into a single group. Moreover, ABCDE and AGL6 genes appeared in the following order: AP3/PI → AG/SHP/STK → AGL6/SEP/AP1. In this study, we collected candidate sequences from gymnosperms and angiosperms. This study highlights important events in the evolutionary history of the ABCDE and AGL6 gene families and clarifies their evolutionary path.

Genomics of Flower Identity in Grapevine (Vitis vinifera L.)

The identity of the four characteristic whorls of typical eudicots, namely, sepals, petals, stamens, and carpels, is specified by the overlapping action of homeotic genes, whose single and combined contributions have been described in detail in the so-called ABCDE model. Continuous species-specific refinements and translations resulted in this model providing the basis for understanding the genetic and molecular mechanisms of flower development in model organisms, such as Arabidopsis thaliana and other main plant species. Although grapevine (Vitis vinifera L.) represents an extremely important cultivated fruit crop globally, studies related to the genetic determinism of flower development are still rare, probably because of the limited interest in sexual reproduction in a plant that is predominantly propagated asexually. Nonetheless, several studies have identified and functionally characterized some ABCDE orthologs in grapevine. The present study is intended to provide a comprehensive screenshot of the transcriptional behavior of 18 representative grapevine ABCDE genes encoding MADS-box transcription factors in a developmental kinetic process, from preanthesis to the postfertilization stage and in different flower organs, namely, the calyx, calyptra, anthers, filaments, ovary, and embryos. The transcript levels found were compared with the proposed model for Arabidopsis to evaluate their biological consistency. With a few exceptions, the results confirmed the expression pattern expected based on the Arabidopsis data.

Floral MADS-box protein interactions in the early diverging angiosperm Aristolochia fimbriata Cham. (Aristolochiaceae: Piperales)

Floral identity MADS-box A, B, C, D, E, and AGL6 class genes are predominantly single copy in Magnoliids, and predate the whole genome duplication (WGD) events in monocots and eudicots. By comparison with the model species Arabidopsis thaliana, the expression patterns of B-, C-, and D-class genes in stamen, carpel, and ovules are conserved in Aristolochia fimbriata, whereas A-, E-class, and AGL6 genes have different expression patterns. Nevertheless, the interactions of these proteins that act through multimeric complexes remain poorly known in early divergent angiosperms. This study evaluates protein interactions among all floral MADS-box A. fimbriata proteins using the Yeast Two Hybrid System (Y2H). We found no homodimers and less heterodimers formed by AfimFUL when compared to AfimAGL6, which allowed us to suggest AGL6 homodimers in combination with AfimSEP2 as the most likely tetramer in sepal identity. We found AfimAP3-AfimPI obligate heterodimers and AfimAG-AfimSEP2 protein interactions intact suggesting conserved stamen and carpel tetrameric complexes in A. fimbriata. We observed a broader interaction partner set for AfimSEP2 than for its paralog AfimSEP1. We show conserved and exclusive MADS-box protein interactions in A. fimbriata in comparison with other eudicot and monocot model species in order to establish plesiomorphic MADS-box protein floral networks in angiosperms.

Evolution and Diversification of FRUITFULL Genes in Solanaceae

Ecologically and economically important fleshy edible fruits have evolved from dry fruit numerous times during angiosperm diversification. However, the molecular mechanisms that underlie these shifts are unknown. In the Solanaceae there has been a major shift to fleshy fruits in the subfamily Solanoideae. Evidence suggests that an ortholog of FRUITFULL (FUL), a transcription factor that regulates cell proliferation and limits the dehiscence zone in the silique of Arabidopsis, plays a similar role in dry-fruited Solanaceae. However, studies have shown that FUL orthologs have taken on new functions in fleshy fruit development, including regulating elements of tomato ripening such as pigment accumulation. FUL belongs to the core eudicot euFUL clade of the angiosperm AP1/FUL gene lineage. The euFUL genes fall into two paralogous clades, euFULI and euFULII. While most core eudicots have one gene in each clade, Solanaceae have two: FUL1 and FUL2 in the former, and MBP10 and MBP20 in the latter. We characterized the evolution of the euFUL genes to identify changes that might be correlated with the origin of fleshy fruit in Solanaceae. Our analyses revealed that the Solanaceae FUL1 and FUL2 clades probably originated through an early whole genome multiplication event. By contrast, the data suggest that the MBP10 and MBP20 clades are the result of a later tandem duplication event. MBP10 is expressed at weak to moderate levels, and its atypical short first intron lacks putative transcription factor binding sites, indicating possible pseudogenization. Consistent with this, our analyses show that MBP10 is evolving at a faster rate compared to MBP20. Our analyses found that Solanaceae euFUL gene duplications, evolutionary rates, and changes in protein residues and expression patterns are not correlated with the shift in fruit type. This suggests deeper analyses are needed to identify the mechanism underlying the change in FUL ortholog function.

HISTONE DEACETYLASE 19 and the flowering time gene FD maintain reproductive meristem identity in an age-dependent manner

The shoot apical meristem (SAM) undergoes developmental transitions that include a shift from vegetative to reproductive growth. This transition is triggered by flowering time genes, which up-regulate floral meristem (FM) identity genes that, in turn, control flower development by activating floral organ identity genes. This cascade of transcriptional activation is refined by repression mechanisms that temporally and spatially restrict gene expression to ensure proper development. Here, we demonstrate that HISTONE DEACETYLASE 19 (HDA19) maintains the identity of the reproductive SAM, or inflorescence meristem (IM), late in Arabidopsis thaliana development. At late stages of growth, hda19 IMs display a striking patterning defect characterized by ectopic expression of floral organ identity genes and the replacement of flowers with individual stamenoid organs. We further show that the flowering time gene FD has a specific function in this regulatory process, as fd hastens the emergence of these patterning defects in hda19 growth. Our work therefore identifies a new role for FD in reproductive patterning, as FD regulates IM function together with HDA19 in an age-dependent fashion. To effect these abnormalities, hda19 and fd may accentuate the weakening of transcriptional repression that occurs naturally with reproductive meristem proliferation.

The parthenocarpic gene Pat-k is generated by a natural mutation of SlAGL6 affecting fruit development in tomato (Solanum lycopersicum L.)

BACKGROUND

Parthenocarpy is a desired trait in tomato because it can overcome problems with fruit setting under unfavorable environmental conditions. A parthenocarpic tomato cultivar, 'MPK-1', with a parthenocarpic gene, Pat-k, exhibits stable parthenocarpy that produces few seeds. Because 'MPK-1' produces few seeds, seedlings are propagated inefficiently via cuttings. It was reported that Pat-k is located on chromosome 1. However, the gene had not been isolated and the relationship between the parthenocarpy and low seed set in 'MPK-1' remained unclear. In this study, we isolated Pat-k to clarify the relationship between parthenocarpy and low seed set in 'MPK-1'.

RESULTS

Using quantitative trait locus (QTL) analysis for parthenocarpy and seed production, we detected a major QTL for each trait on nearly the same region of the Pat-k locus on chromosome 1. To isolate Pat-k, we performed fine mapping using an F₄ population following the cross between a non-parthenocarpic cultivar, 'Micro-Tom' and 'MPK-1'. The results showed that Pat-k was located in the 529 kb interval between two markers, where 60 genes exist. By using data from a whole genome re-sequencing and genome sequence analysis of 'MPK-1', we could identify that the SlAGAMOUS-LIKE 6 (SlAGL6) gene of 'MPK-1' was mutated by a retrotransposon insertion. The transcript level of SlAGL6 was significantly lower in ovaries of 'MPK-1' than a non-parthenocarpic cultivar. From these results, we could conclude that Pat-k is SlAGL6, and its down-regulation in 'MPK-1' causes parthenocarpy and low seed set. In addition, we observed abnormal micropyles only in plants homozygous for the 'MPK-1' allele at the Pat-k/SlAGL6 locus. This result suggests that Pat-k/SlAGL6 is also related to ovule formation and that the low seed set in 'MPK-1' is likely caused by abnormal ovule formation through down-regulation of Pat-k/SlAGL6.

CONCLUSIONS

Pat-k is identical to SlAGL6, and its down-regulation causes parthenocarpy and low seed set in 'MPK-1'. Moreover, down-regulation of Pat-k/SlAGL6 could cause abnormal ovule formation, leading to a reduction in the number of seeds.

Ovule identity mediated by pre-mRNA processing in Arabidopsis

Ovules are fundamental for plant reproduction and crop yield as they are the precursors of seeds. Therefore, ovule specification is a critical developmental program. In Arabidopsis thaliana, ovule identity is redundantly conferred by the homeotic D-class genes SHATTERPROOF1 (SHP1), SHP2 and SEEDSTICK (STK), phylogenetically related to the MADS-domain regulatory gene AGAMOUS (AG), essential in floral organ specification. Previous studies have shown that the HUA-PEP activity, comprised of a suite of RNA-binding protein (RBP) encoding genes, regulates AG pre-mRNA processing and thus flower patterning and organ identity. Here, we report that the HUA-PEP activity additionally governs ovule morphogenesis. Accordingly, in severe hua-pep backgrounds ovules transform into flower organ-like structures. These homeotic transformations are most likely due to the dramatic reduction in SHP1, SHP2 and STK activity. Our molecular and genome-wide profiling strategies revealed the accumulation of prematurely terminated transcripts of D-class genes in hua-pep mutants and reduced amounts of their respective functional messengers, which points to pre-mRNA processing misregulation as the origin of the ovule developmental defects in such backgrounds. RNA processing and transcription are coordinated by the RNA polymerase II (RNAPII) carboxyl-terminal domain (CTD). Our results show that HUA-PEP activity members can interact with the CTD regulator C-TERMINAL DOMAIN PHOSPHATASE-LIKE1 (CPL1), supporting a co-transcriptional mode of action for the HUA-PEP activity. Our findings expand the portfolio of reproductive developmental programs in which HUA-PEP activity participates, and further substantiates the importance of RNA regulatory mechanisms (pre-mRNA co-transcriptional regulation) for correct gene expression during plant morphogenesis.

I Want to (Bud) Break Free: The Potential Role of DAM and SVP-Like Genes in Regulating Dormancy Cycle in Temperate Fruit Trees

Bud dormancy is an adaptive process that allows trees to survive the hard environmental conditions that they experience during the winter of temperate climates. Dormancy is characterized by the reduction in meristematic activity and the absence of visible growth. A prolonged exposure to cold temperatures is required to allow the bud resuming growth in response to warm temperatures. In fruit tree species, the dormancy cycle is believed to be regulated by a group of genes encoding MADS-box transcription factors. These genes are called DORMANCY-ASSOCIATED MADS-BOX (DAM) and are phylogenetically related to the Arabidopsis thaliana floral regulators SHORT VEGETATIVE PHASE (SVP) and AGAMOUS-LIKE 24. The interest in DAM and other orthologs of SVP (SVP-like) genes has notably increased due to the publication of several reports suggesting their role in the control of bud dormancy in numerous fruit species, including apple, pear, peach, Japanese apricot, and kiwifruit among others. In this review, we briefly describe the physiological bases of the dormancy cycle and how it is genetically regulated, with a particular emphasis on DAM and SVP-like genes. We also provide a detailed report of the most recent advances about the transcriptional regulation of these genes by seasonal cues, epigenetics and plant hormones. From this information, we propose a tentative classification of DAM and SVP-like genes based on their seasonal pattern of expression. Furthermore, we discuss the potential biological role of DAM and SVP-like genes in bud dormancy in antagonizing the function of FLOWERING LOCUS T-like genes. Finally, we draw a global picture of the possible role of DAM and SVP-like genes in the bud dormancy cycle and propose a model that integrates these genes in a molecular network of dormancy cycle regulation in temperate fruit trees.

Arabidopsis thaliana MLO genes are expressed in discrete domains during reproductive development

MLOs in Plant Reproduction. The MILDEW RESISTANCE LOCUS-O (MLO) protein family, comprised of 15 members, plays roles in diverse cell-cell communication processes such as powdery mildew susceptibility, root thigmomorphogenesis, and pollen tube reception. The NORTIA (NTA, AtMLO7) gene is expressed in the synergid cells of the female gametophyte where it functions in intercellular communication with the pollen tube. Discrepancies between previously published promoter::GUS and promoter::gene-GUS constructs expression patterns led us to explore the regulation of NTA expression. Here we found via NTA_pro::gNTA-GUS truncations that sequences within the NTA gene negatively regulate its expression in the stomata and carpel walls. This led to the hypothesis that other MLO family members may also have additional regulatory sequences within the gene. MLO_pro::gMLO-GUS constructs were examined for each family member focusing specifically on flowers in order to determine whether other MLOs could play a role in reproductive cell-cell communication. Notably, several MLOs were expressed in the pollen, in the stigma, in the pollinated style, and in the synergids and central cell. These findings indicate that other MLOs in addition to NTA could play a role in reproduction. Previous studies on the MLO family showed that phylogenetically related MLOs had redundant functions in powdery mildew infection and root thigmomorphogenesis; however, MLO expression in reproductive tissues did not strictly follow phylogenetic relationships, indicating that MLOs from different evolutionary origins may have been recruited for function in sexual reproduction.

Tomato facultative parthenocarpy results from SlAGAMOUS-LIKE 6 loss of function

The extreme sensitivity of the microsporogenesis process to moderately high or low temperatures is a major hindrance for tomato (Solanum lycopersicum) sexual reproduction and hence year-round cropping. Consequently, breeding for parthenocarpy, namely, fertilization-independent fruit set, is considered a valuable goal especially for maintaining sustainable agriculture in the face of global warming. A mutant capable of setting high-quality seedless (parthenocarpic) fruit was found following a screen of EMS-mutagenized tomato population for yielding under heat stress. Next-generation sequencing followed by marker-assisted mapping and CRISPR/Cas9 gene knockout confirmed that a mutation in SlAGAMOUS-LIKE 6 (SlAGL6) was responsible for the parthenocarpic phenotype. The mutant is capable of fruit production under heat stress conditions that severely hamper fertilization-dependent fruit set. Different from other tomato recessive monogenic mutants for parthenocarpy, Slagl6 mutations impose no homeotic changes, the seedless fruits are of normal weight and shape, pollen viability is unaffected, and sexual reproduction capacity is maintained, thus making Slagl6 an attractive gene for facultative parthenocarpy. The characteristics of the analysed mutant combined with the gene's mode of expression imply SlAGL6 as a key regulator of the transition between the state of 'ovary arrest' imposed towards anthesis and the fertilization-triggered fruit set.

Natural haplotypes of FLM non-coding sequences fine-tune flowering time in ambient spring temperatures in Arabidopsis

Cool ambient temperatures are major cues determining flowering time in spring. The mechanisms promoting or delaying flowering in response to ambient temperature changes are only beginning to be understood. In Arabidopsis thaliana, FLOWERING LOCUS M (FLM) regulates flowering in the ambient temperature range and FLM is transcribed and alternatively spliced in a temperature-dependent manner. We identify polymorphic promoter and intronic sequences required for FLM expression and splicing. In transgenic experiments covering 69% of the available sequence variation in two distinct sites, we show that variation in the abundance of the FLM-ß splice form strictly correlate (R² = 0.94) with flowering time over an extended vegetative period. The FLM polymorphisms lead to changes in FLM expression (PRO2+) but may also affect FLM intron 1 splicing (INT6+). This information could serve to buffer the anticipated negative effects on agricultural systems and flowering that may occur during climate change.

DOI:http://dx.doi.org/10.7554/eLife.22114.001

Transcriptome Analysis of a Female-sterile Mutant (fsm) in Chinese Cabbage (Brassica campestris ssp. pekinensis)

Female-sterile mutants are ideal materials for studying pistil development in plants. Here, we identified a female-sterile mutant fsm in Chinese cabbage. This mutant, which exhibited stable inheritance, was derived from Chinese cabbage DH line 'FT' using a combination of isolated microspore culture and ethyl methanesulfonate mutagenesis. Compared with the wild-type line 'FT,' the fsm plants exhibited pistil abortion, and floral organs were also relatively smaller. Genetic analysis indicated that the phenotype of fsm is controlled by a single recessive nuclear gene. Morphological observations revealed that the presence of abnormal ovules in fsm likely influenced normal fertilization process, ultimately leading to female sterility. Comparative transcriptome analysis on the flower buds of 'FT' and fsm using RNA-Seq revealed a total of 1,872 differentially expressed genes (DEGs). Of these, a number of genes involved in pistil development were identified, such as PRETTY FEW SEEDS 2 (PFS2), temperature-induced lipocalin (TIL), AGAMOUS-LIKE (AGL), and HECATE (HEC). Furthermore, GO and KEGG pathway enrichment analyses of the DEGs suggested that a variety of biological processes and metabolic pathways are significantly enriched during pistil development. In addition, the expression patterns of 16 DEGs, including four pistil development-related genes and 12 floral organ development-related genes, were analyzed using qRT-PCR. A total of 31,272 single nucleotide polymorphisms were specifically detected in fsm. These results contribute to shed light on the regulatory mechanisms underlying pistil development in Chinese cabbage.

Bioinformatic prediction of transcription factor binding sites at promoter regions of genes for photoperiod and vernalization responses in model and temperate cereal plants

Background

Many genes involved in responses to photoperiod and vernalization have been characterized or predicted in Arabidopsis (Arabidopsis thaliana), Brachypodium (Brachypodium distachyon), wheat (Triticum aestivum) and barley (Hordeum vulgare). However, little is known about the transcription regulation of these genes, especially in the large, complex genomes of wheat and barley.

Results

We identified 68, 60, 195 and 61 genes that are known or postulated to control pathways of photoperiod (PH), vernalization (VE) and pathway integration (PI) in Arabidopsis, Brachypodium, wheat and barley for predicting transcription factor binding sites (TFBSs) in the promoters of these genes using the FIMO motif search tool of the MEME Suite. The initial predicted TFBSs were filtered to confirm the final numbers of predicted TFBSs to be 1066, 1379, 1528, and 789 in Arabidopsis, Brachypodium, wheat and barley, respectively. These TFBSs were mapped onto the PH, VE and PI pathways to infer about the regulation of gene expression in Arabidopsis and cereal species. The GC contents in promoters, untranslated regions (UTRs), coding sequences and introns were higher in the three cereal species than those in Arabidopsis. The predicted TFBSs were most abundant for two transcription factor (TF) families: MADS-box and CSD (cold shock domain). The analysis of publicly available gene expression data showed that genes with similar numbers of MADS-box and CSD TFBSs exhibited similar expression patterns across several different tissues and developmental stages. The intra-specific Tajima D-statistics of TFBS motif diversity showed different binding specificity among different TF families. The inter-specific Tajima D-statistics suggested faster TFBS divergence in TFBSs than in coding sequences and introns. Mapping TFBSs onto the PH, VE and PI pathways showed the predominance of MADS-box and CSD TFBSs in most genes of the four species, and the difference in the pathway regulations between Arabidopsis and the three cereal species.

Conclusion

Our approach to associating the key flowering genes with their potential TFs through prediction of putative TFBSs provides a framework to explore regulatory mechanisms of photoperiod and vernalization responses in flowering plants. The predicted TFBSs in the promoters of the flowering genes provide a basis for molecular characterization of transcription regulation in the large, complex genomes of important crop species, wheat and barley.

Electronic supplementary material

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

Towards genome-wide prediction and characterization of enhancers in plants

Enhancers are important cis-regulatory DNA elements that regulate transcription programs by recruiting transcription factors and directing them to the promoters of target genes in a cell-type/tissue-specific manner. The expression of a gene can be regulated by one or multiple enhancers at different developmental stages and/or in different tissues. Enhancers are difficult to identify because of their unpredictable positions relative to their cognate promoters. Remarkably, only a handful of enhancers have been identified in plant species largely due to the lack of general approaches for enhancer identification. Extensive genomic and epigenomic research in mammalian species has revealed that the genomic locations of enhancers can be predicted based on the binding sites of transcriptional co-factors and several distinct features associated with open chromatin. Here we review the methodologies used in enhancer prediction in mammalian species. We also review the recent applications of these methodologies in Arabidopsis thaliana and discuss the future directions of enhancer identification in plants. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Flower development: the evolutionary history and functions of the AGL6 subfamily MADS-box genes

AGL6 is an ancient subfamily of MADS-box genes found in both gymnosperms and angiosperms. Its functions remained elusive despite the fact that the MADS-box genes and the ABC model have been studied for >20 years. Nevertheless, recent discoveries in petunia, rice, and maize support its involvement in the 'E' function of floral development, very similar to the closely related AGL2 (SEPALLATA) subfamily which has been well characterized. The known functions of AGL6 span from ancient conserved roles to new functions acquired in specific plant families. The AGL6 genes are involved in floral meristem regulation, in floral organs, and ovule (integument) and seed development, and have possible roles in both male and female germline and gametophyte development. In grasses, they are also important for the development of the first whorl of the flower, whereas in Arabidopsis they may play additional roles before floral meristem formation. This review covers these recent insights and some other aspects that are not yet fully elucidated, which deserve more studies in the future.

Modulation of Ambient Temperature-Dependent Flowering in Arabidopsis thaliana by Natural Variation of FLOWERING LOCUS M

Plants integrate seasonal cues such as temperature and day length to optimally adjust their flowering time to the environment. Compared to the control of flowering before and after winter by the vernalization and day length pathways, mechanisms that delay or promote flowering during a transient cool or warm period, especially during spring, are less well understood. Due to global warming, understanding this ambient temperature pathway has gained increasing importance. In Arabidopsis thaliana, FLOWERING LOCUS M (FLM) is a critical flowering regulator of the ambient temperature pathway. FLM is alternatively spliced in a temperature-dependent manner and the two predominant splice variants, FLM-ß and FLM-δ, can repress and activate flowering in the genetic background of the A. thaliana reference accession Columbia-0. The relevance of this regulatory mechanism for the environmental adaptation across the entire range of the species is, however, unknown. Here, we identify insertion polymorphisms in the first intron of FLM as causative for accelerated flowering in many natural A. thaliana accessions, especially in cool (15°C) temperatures. We present evidence for a potential adaptive role of this structural variation and link it specifically to changes in the abundance of FLM-ß. Our results may allow predicting flowering in response to ambient temperatures in the Brassicaceae.

K-homology nuclear ribonucleoproteins regulate floral organ identity and determinacy in arabidopsis

Post-transcriptional control is nowadays considered a main checking point for correct gene regulation during development, and RNA binding proteins actively participate in this process. Arabidopsis thaliana FLOWERING LOCUS WITH KH DOMAINS (FLK) and PEPPER (PEP) genes encode RNA-binding proteins that contain three K-homology (KH)-domain, the typical configuration of Poly(C)-binding ribonucleoproteins (PCBPs). We previously demonstrated that FLK and PEP interact to regulate FLOWERING LOCUS C (FLC), a central repressor of flowering time. Now we show that FLK and PEP also play an important role in the maintenance of the C-function during floral organ identity by post-transcriptionally regulating the MADS-box floral homeotic gene AGAMOUS (AG). Previous studies have indicated that the KH-domain containing protein HEN4, in concert with the CCCH-type RNA binding protein HUA1 and the RPR-type protein HUA2, facilitates maturation of the AG pre-mRNA. In this report we show that FLK and PEP genetically interact with HEN4, HUA1, and HUA2, and that the FLK and PEP proteins physically associate with HUA1 and HEN4. Taken together, these data suggest that HUA1, HEN4, PEP and FLK are components of the same post-transcriptional regulatory module that ensures normal processing of the AG pre-mRNA. Our data better delineates the roles of PEP in plant development and, for the first time, links FLK to a morphogenetic process.

Unlike animals, angiosperms (flowering plants) lack a germline that is set-aside early in embryo development. Contrariwise, reproductive success relies on the formation of flowers during adult life, which provide the germ cells and the means for fertilization. Therefore, timing of flowering and flower organ morphogenesis are critical developmental operations that must be finely regulated and coordinated to complete reproduction. Arabidopsis thaliana FLOWERING LOCUS WITH KH DOMAINS (FLK) and PEPPER (PEP) encode two KH-domain RNA-binding proteins phylogenetically related to human proteins characterized by their high developmental versatility. FLK and PEP modulate the mRNA expression of the MADS-box gene FLOWERING LOCUS C, key in flowering control. In this work we have found that FLK and PEP also play a pivotal role in flower organogenesis by post-transcriptionally regulating the MADS-box floral organ identity gene AGAMOUS (AG). Interestingly, FLK and PEP physically interact with proteins involved in AG pre-mRNA processing to secure correct AG function in the floral meristem and flower. Taken together, our results reveal the existence of a post-transcriptional regulatory activity controlling key master genes for floral timing and flower morphogenesis, which might be instrumental for coordinating both developmental phases.

Flower Development and Perianth Identity Candidate Genes in the Basal Angiosperm Aristolochia fimbriata (Piperales: Aristolochiaceae)

Aristolochia fimbriata (Aristolochiaceae: Piperales) exhibits highly synorganized flowers with a single convoluted structure forming a petaloid perianth that surrounds the gynostemium, putatively formed by the congenital fusion between stamens and the upper portion of the carpels. Here we present the flower development and morphology of A. fimbriata, together with the expression of the key regulatory genes that participate in flower development, particularly those likely controlling perianth identity. A. fimbriata is a member of the magnoliids, and thus gene expression detected for all ABCE MADS-box genes in this taxon, can also help to elucidate patterns of gene expression prior the independent duplications of these genes in eudicots and monocots. Using both floral development and anatomy in combination with the isolation of MADS-box gene homologs, gene phylogenetic analyses and expression studies (both by reverse transcription PCR and in situ hybridization), we present hypotheses on floral organ identity genes involved in the formation of this bizarre flower. We found that most MADS-box genes were expressed in vegetative and reproductive tissues with the exception of AfimSEP2, AfimAGL6, and AfimSTK transcripts that are only found in flowers and capsules but are not detected in leaves. Two genes show ubiquitous expression; AfimFUL that is found in all floral organs at all developmental stages as well as in leaves and capsules, and AfimAG that has low expression in leaves and is found in all floral organs at all stages with a considerable reduction of expression in the limb of anthetic flowers. Our results indicate that expression of AfimFUL is indicative of pleiotropic roles and not of a perianth identity specific function. On the other hand, expression of B-class genes, AfimAP3 and AfimPI, suggests their conserved role in stamen identity and corroborates that the perianth is sepal and not petal-derived. Our data also postulates an AGL6 ortholog as a candidate gene for sepal identity in the Aristolochiaceae and provides testable hypothesis for a modified ABCE model in synorganized magnoliid flowers.

Analysis of CACTA transposases reveals intron loss as major factor influencing their exon/intron structure in monocotyledonous and eudicotyledonous hosts

Background

CACTA elements are DNA transposons and are found in numerous organisms. Despite their low activity, several thousand copies can be identified in many genomes. CACTA elements transpose using a ‘cut-and-paste’ mechanism, which is facilitated by a DDE transposase. DDE transposases from CACTA elements contain, despite their conserved function, different exon numbers among various CACTA families. While earlier studies analyzed the ancestral history of the DDE transposases, no studies have examined exon loss and gain with a view of mechanisms that could drive the changes.

Results

We analyzed 64 transposases from different CACTA families among monocotyledonous and eudicotyledonous host species. The annotation of the exon/intron boundaries showed a range from one to six exons. A robust multiple sequence alignment of the 64 transposases based on their protein sequences was created and used for phylogenetic analysis, which revealed eight different clades. We observed that the exon numbers in CACTA transposases are not specific for a host genome. We found that ancient CACTA lineages diverged before the divergence of monocotyledons and eudicotyledons. Most exon/intron boundaries were found in three distinct regions among all the transposases, grouping 63 conserved intron/exon boundaries.

Conclusions

We propose a model for the ancestral CACTA transposase gene, which consists of four exons, that predates the divergence of the monocotyledons and eudicotyledons. Based on this model, we propose pathways of intron loss or gain to explain the observed variation in exon numbers. While intron loss appears to have prevailed, a putative case of intron gain was nevertheless observed.

AGAMOUS-LIKE13, a putative ancestor for the E functional genes, specifies male and female gametophyte morphogenesis

Arabidopsis AGL13 is a member of the AGL6 clade of the MADS box gene family. GUS activity was specifically detected from the initiation to maturation of both pollen and ovules in AGL13:GUS Arabidopsis. The sterility of the flower with defective pollen and ovules was found in AGL13 RNAi knockdown and AGL13 + SRDX dominant-negative mutants. These results indicate that AGL13 acts as an activator in regulation of early initiation and further development of pollen and ovules. The production of similar floral organ defects in the severe AGL13 + SRDX and SEP2 + SRDX plants and the similar enhancement of AG nuclear localization efficiency by AGL13 and SEP3 proteins suggest a similar function for AGL13 and E functional SEP proteins. Additional fluorescence resonance energy transfer (FRET) analysis indicated that, similar to SEP proteins, AGL13 is able to interact with AG to form quartet-like complexes (AGL13-AG)2 and interact with AG-AP3-PI to form a higher-order heterotetrameric complex (AGL13-AG-AP3-PI). Through these complexes, AGL13 and AG could regulate the expression of similar downstream genes involved in pollen morphogenesis, anther cell layer formation and the ovule development. AGL13 also regulates AG/AP3/PI expression by positive regulatory feedback loops and suppresses its own expression through negative regulatory feedback loops by activating AGL6, which acts as a repressor of AGL13. Our data suggest that AGL13 is likely a putative ancestor for the E functional genes which specifies male and female gametophyte morphogenesis in plants during evolution.

Impact of an intragenic retrotransposon on the structural integrity and evolution of a major isoprenoid biosynthesis pathway gene in Hevea brasiliensis

Isoprenoids belong to a large family of structurally and functionally different natural compounds found universally from prokaryotes to higher animals and plants. In Hevea brasiliensis, the commercially important cis-polyisoprene (rubber) is synthesised as part of its defence mechanism in addition to other common isoprenoids like phytosterols, growth hormones etc. Farnesyl diphosphate synthase (FDPS) is a key enzyme in this process which catalyses the conversion of isoprene units into polyisoprene. Although prior sequence information is available, the structural variants of the FDPS gene presently existing in Hevea population are largely unknown. Since gene structure has a major role in gene regulation, extensive sequence analysis of this gene from different genotypes was carried out to identify the prevailing structural variants. We identified several SNPs and large indels which were associated with a partial transposable element (TE). Modification of key regulatory motifs and splice sites induced by the retroelement was also identified in the first intron. Screening of popular rubber clones, wild germplasm accessions and Hevea species revealed that the retroelement is responsible for the generation of new alleles with varying degrees of sequence homology. Segregation analysis of a progeny population confirmed that the alleles are not paralogs and are inherited in a Mendelian mode. Our findings suggest that the first intron of the FDPS gene has been subjected to various chromosomal rearrangements due to the interaction of a retrotransposon, resulting in novel alleles which may substantially contribute towards the evolution of this major gene in rubber. Moreover, the results indicate the possible existence of a retrotransposon-mediated epigenetic gene regulatory mechanism in Hevea.

Homeotic Genes and the ABCDE Model for Floral Organ Formation in Wheat

Floral organ formation has been the subject of intensive study for over 20 years, particularly in the model dicot species Arabidopsis thaliana. These studies have led to the establishment of a general model for the development of floral organs in higher plants, the so-called ABCDE model, in which floral whorl-specific combinations of class A, B, C, D, or E genes specify floral organ identity. In Arabidopsis, class A, B, C, D, E genes encode MADS-box transcription factors except for the class A gene APETALA2. Mutation of these genes induces floral organ homeosis. In this review, I focus on the roles of these homeotic genes in bread wheat (Triticum aestivum), particularly with respect to the ABCDE model. Pistillody, the homeotic transformation of stamens into pistil-like structures, occurs in cytoplasmic substitution (alloplasmic) wheat lines that have the cytoplasm of the related wild species Aegilops crassa. This phenomenon is a valuable tool for analysis of the wheat ABCDE model. Using an alloplasmic line, the wheat ortholog of DROOPING LEAF (TaDL), a member of the YABBY gene family, has been shown to regulate pistil specification. Here, I describe the current understanding of the ABCDE model for floral organ formation in wheat.

Potential functional replacement of the plastidic acetyl-CoA carboxylase subunit (accD) gene by recent transfers to the nucleus in some angiosperm lineages

Eukaryotic cells originated when an ancestor of the nucleated cell engulfed bacterial endosymbionts that gradually evolved into the mitochondrion and the chloroplast. Soon after these endosymbiotic events, thousands of ancestral prokaryotic genes were functionally transferred from the endosymbionts to the nucleus. This process of functional gene relocation, now rare in eukaryotes, continues in angiosperms. In this article, we show that the chloroplastic acetyl-CoA carboxylase subunit (accD) gene that is present in the plastome of most angiosperms has been functionally relocated to the nucleus in the Campanulaceae. Surprisingly, the nucleus-encoded accD transcript is considerably smaller than the plastidic version, consisting of little more than the carboxylase domain of the plastidic accD gene fused to a coding region encoding a plastid targeting peptide. We verified experimentally the presence of a chloroplastic transit peptide by showing that the product of the nuclear accD fused to green fluorescent protein was imported in the chloroplasts. The nuclear gene regulatory elements that enabled the erstwhile plastidic gene to become functional in the nuclear genome were identified, and the evolution of the intronic and exonic sequences in the nucleus is described. Relocation and truncation of the accD gene is a remarkable example of the processes underpinning endosymbiotic evolution.

Comparative genome analysis of two Cryptosporidium parvum isolates with different host range

Parasites of the genus Cryptosporidium infect the intestinal and gastric epithelium of different vertebrate species. Some of the many Cryptosporidium species described to date differ with respect to host range; whereas some species' host range appears to be narrow, others have been isolated from taxonomically unrelated vertebrates. To begin to investigate the genetic basis of Cryptosporidium host specificity, the genome of a Cryptosporidium parvum isolate belonging to a sub-specific group found exclusively in humans was sequenced and compared to the reference C. parvum genome representative of the zoonotic group. Over 12,000 single-nucleotide polymorphisms (SNPs), or 1.4 SNP per kilobase, were identified. The genome distribution of SNPs was highly heterogeneous, but non-synonymous and silent SNPs were similarly distributed. On many chromosomes, the most highly divergent regions were located near the ends. Genes in the most diverged regions were almost twice as large as the genome-wide average. Transporters, and ABC transporters in particular, were over-represented among these genes, as were proteins with predicted signal peptide. Possibly reflecting the presence of regulatory sequences, the distribution of intergenic SNPs differed according to the function of the downstream open reading frame. A 3-way comparison of the newly sequenced anthroponotic C. parvum, the reference zoonotic C. parvum and the human parasite Cryptosporidium hominis identified genetic loci where the anthroponotic C. parvum sequence is more similar to C. hominis than to the zoonotic C. parvum reference. Because C. hominis and anthroponotic C. parvum share a similar host range, this unexpected observation suggests that proteins encoded by these genes may influence the host range.

DOF-binding sites additively contribute to guard cell-specificity of AtMYB60 promoter

BACKGROUND

We previously demonstrated that the Arabidopsis thaliana AtMYB60 protein is an R2R3MYB transcription factor required for stomatal opening. AtMYB60 is specifically expressed in guard cells and down-regulated at the transcriptional levels by the phytohormone ABA.

RESULTS

To investigate the molecular mechanisms governing AtMYB60 expression, its promoter was dissected through deletion and mutagenesis analyses. By studying different versions of AtMYB60 promoter::GUS reporter fusions in transgenic plants we were able to demonstrate a modular organization for the AtMYB60 promoter. Particularly we defined: a minimal promoter sufficient to confer guard cell-specific activity to the reporter gene; the distinct roles of different DOF-binding sites organised in a cluster in the minimal promoter in determining guard cell-specific expression; the promoter regions responsible for the enhancement of activity in guard cells; a promoter region responsible for the negative transcriptional regulation by ABA. Moreover from the analysis of single and multiple mutants we could rule out the involvement of a group of DOF proteins, known as CDFs, already characterised for their involvement in flowering time, in the regulation of AtMYB60 expression.

CONCLUSIONS

These findings shed light on the regulation of gene expression in guard cells and provide new promoter modules as useful tools for manipulating gene expression in guard cells, both for physiological studies and future biotechnological applications.

Specific expression of LATERAL SUPPRESSOR is controlled by an evolutionarily conserved 3' enhancer

Aerial plant architecture is largely based on the activity of axillary meristems (AMs), initiated in the axils of leaves. The Arabidopsis gene LATERAL SUPPRESSOR (LAS), which is expressed in well-defined domains at the adaxial boundary of leaf primordia, is a key regulator of AM formation. The precise definition of organ boundaries is an essential step for the formation of new organs in general and for meristem initiation; however, mechanisms leading to these specific patterns are not well understood. To increase understanding of how the highly specific transcript accumulation in organ boundary regions is established, we investigated the LAS promoter. Analysis of deletion constructs revealed that an essential enhancer necessary for complementation is situated about 3.2 kb downstream of the LAS open reading frame. This enhancer is sufficient to confer promoter specificity as upstream sequences in LAS could be replaced by non-specific promoters, such as the 35S minimal promoter. Further promoter swapping experiments using the PISTILLATA or the full 35S promoter demonstrated that the LAS 3' enhancer also has suppressor functions, largely overwriting the activity of different 5' promoters. Phylogenetic analyses suggest that LAS function and regulation are evolutionarily highly conserved. Homologous elements in downstream regulatory sequences were found in all LAS orthologs, including grasses. Transcomplementation experiments demonstrated the functional conservation of non-coding sequences between Solanum lycopersicum (tomato) and Arabidopsis. In summary, our results show that a highly conserved enhancer/suppressor element is the main regulatory module conferring the boundary-specific expression of LAS.

Transcriptome analysis of the Arabidopsis megaspore mother cell uncovers the importance of RNA helicases for plant germline development

Germ line specification is a crucial step in the life cycle of all organisms. For sexual plant reproduction, the megaspore mother cell (MMC) is of crucial importance: it marks the first cell of the plant "germline" lineage that gets committed to undergo meiosis. One of the meiotic products, the functional megaspore, subsequently gives rise to the haploid, multicellular female gametophyte that harbours the female gametes. The MMC is formed by selection and differentiation of a single somatic, sub-epidermal cell in the ovule. The transcriptional network underlying MMC specification and differentiation is largely unknown. We provide the first transcriptome analysis of an MMC using the model plant Arabidopsis thaliana with a combination of laser-assisted microdissection and microarray hybridizations. Statistical analyses identified an over-representation of translational regulation control pathways and a significant enrichment of DEAD/DEAH-box helicases in the MMC transcriptome, paralleling important features of the animal germline. Analysis of two independent T-DNA insertion lines suggests an important role of an enriched helicase, MNEME (MEM), in MMC differentiation and the restriction of the germline fate to only one cell per ovule primordium. In heterozygous mem mutants, additional enlarged MMC-like cells, which sometimes initiate female gametophyte development, were observed at higher frequencies than in the wild type. This closely resembles the phenotype of mutants affected in the small RNA and DNA-methylation pathways important for epigenetic regulation. Importantly, the mem phenotype shows features of apospory, as female gametophytes initiate from two non-sister cells in these mutants. Moreover, in mem gametophytic nuclei, both higher order chromatin structure and the distribution of LIKE HETEROCHROMATIN PROTEIN1 were affected, indicating epigenetic perturbations. In summary, the MMC transcriptome sets the stage for future functional characterization as illustrated by the identification of MEM, a novel gene involved in the restriction of germline fate.

The Arabidopsis SOC1-like genes AGL42, AGL71 and AGL72 promote flowering in the shoot apical and axillary meristems

The floral transition is the switch from vegetative development to flowering. Proper timing of the floral transition is regulated by different pathways and is critical for the reproductive success of plants. Some of the flowering pathways are controlled by environmental signals such as photoperiod and vernalization, others by autonomous signals such as the developmental state of the plant and hormones, particularly gibberellin. SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1) acts in Arabidopsis as an integrative centre of these pathways, promoting the floral transition. In this work, we show that AGAMOUS-LIKE 42 (AGL42), AGAMOUS-LIKE 71 (AGL71) and AGAMOUS-LIKE 72 (AGL72), which encode MADS-box transcription factors phylogenetically closely related to SOC1, are also involved in the floral transition. They promote flowering at the shoot apical and axillary meristems and seem to act through a gibberellin-dependent pathway. Furthermore SOC1 directly controls the expression of AGL42, AGL71 and AGL72 to balance the expression level of these SOC1-like genes. Our data reveal roles for AGL42, AGL71 and AGL72 in the floral transition, demonstrate their genetic interactions with SOC1 and suggest that their roles differ in the apical and axillary meristems.

BAC-recombineering for studying plant gene regulation: developmental control and cellular localization of SnRK1 kinase subunits

Recombineering, permitting precise modification of genes within bacterial artificial chromosomes (BACs) through homologous recombination mediated by lambda phage-encoded Red proteins, is a widely used powerful tool in mouse, Caenorhabditis and Drosophila genetics. As Agrobacterium-mediated transfer of large DNA inserts from binary BACs and TACs into plants occurs at low frequency, recombineering is so far seldom exploited in the analysis of plant gene functions. We have constructed binary plant transformation vectors, which are suitable for gap-repair cloning of genes from BACs using recombineering methods previously developed for other organisms. Here we show that recombineering facilitates PCR-based generation of precise translational fusions between coding sequences of fluorescent reporter and plant proteins using galK-based exchange recombination. The modified target genes alone or as part of a larger gene cluster can be transferred by high-frequency gap-repair into plant transformation vectors, stably maintained in Agrobacterium and transformed without alteration into plants. Versatile application of plant BAC-recombineering is illustrated by the analysis of developmental regulation and cellular localization of interacting AKIN10 catalytic and SNF4 activating subunits of Arabidopsis Snf1-related (SnRK1) protein kinase using in vivo imaging. To validate full functionality and in vivo interaction of tagged SnRK1 subunits, it is demonstrated that immunoprecipitated SNF4-YFP is bound to a kinase that phosphorylates SnRK1 candidate substrates, and that the GFP- and YFP-tagged kinase subunits co-immunoprecipitate with endogenous wild type AKIN10 and SNF4.

RD20, a stress-inducible caleosin, participates in stomatal control, transpiration and drought tolerance in Arabidopsis thaliana

Plants overcome water deficit conditions by combining molecular, biochemical and morphological changes. At the molecular level, many stress-responsive genes have been isolated, but knowledge of their physiological functions remains fragmentary. Here, we report data for RD20, a stress-inducible Arabidopsis gene that belongs to the caleosin family. As for other caleosins, we showed that RD20 localized to oil bodies. Although caleosins are thought to play a role in the degradation of lipids during seed germination, induction of RD20 by dehydration, salt stress and ABA suggests that RD20 might be involved in processes other than germination. Using plants carrying the promoter RD20::uidA construct, we show that RD20 is expressed in leaves, guard cells and flowers, but not in root or in mature seeds. Water deficit triggers a transient increase in RD20 expression in leaves that appeared predominantly dependent on ABA signaling. To assess the biological significance of these data, a functional analysis using rd20 knock-out and overexpressing complemented lines cultivated either in standard or in water deficit conditions was performed. The rd20 knock-out plants present a higher transpiration rate that correlates with enhanced stomatal opening and a reduced tolerance to drought as compared with the wild type. These results support a role for RD20 in drought tolerance through stomatal control under water deficit conditions.

Evolutionary trends in the floral transcriptome: insights from one of the basalmost angiosperms, the water lily Nuphar advena (Nymphaeaceae)

Current understanding of floral developmental genetics comes primarily from the core eudicot model Arabidopsis thaliana. Here, we explore the floral transcriptome of the basal angiosperm, Nuphar advena (water lily), for insights into the ancestral developmental program of flowers. We identify several thousand Nuphar genes with significantly upregulated floral expression, including homologs of the well-known ABCE floral regulators, deployed in broadly overlapping transcriptional programs across floral organ categories. Strong similarities in the expression profiles of different organ categories in Nuphar flowers are shared with the magnoliid Persea americana (avocado), in contrast to the largely organ-specific transcriptional cascades evident in Arabidopsis, supporting the inference that this is the ancestral condition in angiosperms. In contrast to most eudicots, floral organs are weakly differentiated in Nuphar and Persea, with staminodial intermediates between stamens and perianth in Nuphar, and between stamens and carpels in Persea. Consequently, the predominantly organ-specific transcriptional programs that characterize Arabidopsis flowers (and perhaps other eudicots) are derived, and correlate with a shift towards morphologically distinct floral organs, including differentiated sepals and petals, and a perianth distinct from stamens and carpels. Our findings suggest that the genetic regulation of more spatially discrete transcriptional programs underlies the evolution of floral morphology.

Evolutionary analysis of regulatory sequences (EARS) in plants

Identification of regulatory sequences within non-coding regions of DNA is an essential step towards elucidation of gene networks. This approach constitutes a major challenge, however, as only a very small fraction of non-coding DNA is thought to contribute to gene regulation. The mapping of regulatory regions traditionally involves the laborious construction of promoter deletion series which are then fused to reporter genes and assayed in transgenic organisms. Bioinformatic methods can be used to scan sequences for matches for known regulatory motifs, however these methods are currently hampered by the relatively small amount of such motifs and by a high false-discovery rate. Here, we demonstrate a robust and highly sensitive, in silico method to identify evolutionarily conserved regions within non-coding DNA. Sequence conservation within these regions is taken as evidence for evolutionary pressure against mutations, which is suggestive of functional importance. We test this method on a small set of well characterised promoters, and show that it successfully identifies known regulatory regions. We further show that these evolutionarily conserved sequences contain clusters of transcription binding sites, often described as regulatory modules. A version of the tool optimised for the analysis of plant promoters is available online at http://wsbc.warwick.ac.uk/ears/main.php.

Expression divergence of the AGL6 MADS domain transcription factor lineage after a core eudicot duplication suggests functional diversification

BACKGROUND

Because of their known role as transcriptional regulators of key plant developmental processes, the diversification of MADS-box gene function is thought to be a major driving force in the developmental evolution of plants. Yet the function of some MADS-box gene subfamilies has remained elusive thus far. One such lineage, AGL6, has now been functionally characterized in three angiosperm species, but a phylogenetic framework for comparison of AGL6 gene function is currently missing.

RESULTS

Based on phylogenetic analyses of newly isolated and EST-based sequences, we describe the duplication history of the AGL6 subfamily in angiosperms. Our analyses provide support for four ancient duplications in the evolution of the AGL6 lineage: one at the base of core eudicots resulting in euAGL6 and AGL6-like gene clades, one during basal angiosperm diversification and two in monocot evolution. To investigate whether the spatial domains in which AGL6 genes function have diverged after duplication, we use quantitative Real Time PCR. We show that the core eudicot AGL6-like clade acquired expression in vegetative tissues, while its paralog euAGL6 remains predominantly confined to reproductive tissues.

CONCLUSIONS

These and previous data lead us to propose that the AGL6 lineage in core eudicots, in addition to functions related to the expression in reproductive structures, may have acquired a function in developmental transitions of vegetative shoots.