Adaptation of the Invasive Plant Sphagneticola trilobata to Flooding Stress by Hybridization with Native Relatives.

Hybridization is common between invasive and native species and may produce more adaptive hybrids. The hybrid (Sphagneticola × guangdongensis) of Sphagneticola trilobata (an invasive species) and S. calendulacea (a native species) was found in South China. In this study, S. trilobata, S. calendulacea, and Sphagneticola × guangdongensis were used as research materials to explore their adaptability to flooding stress. Under flooding stress, the ethylene content and the expression of key enzyme genes related to ethylene synthesis in Sphagneticola × guangdongensis and S. calendulacea were significantly higher than those in S. trilobata. A large number of adventitious roots and aerenchyma were generated in Sphagneticola × guangdongensis and S. calendulacea. The contents of reactive oxygen species and malondialdehyde in Sphagneticola × guangdongensis and S. calendulacea were lower than those in S. trilobata, and the leaves of S. trilobata were the most severely damaged under flooding stress. The results indicate that hybridization catalyzed the tolerance of Sphagneticola × guangdongensis to flooding stress, and the responses of Sphagneticola × guangdongensis to flooding stress were more similar to that of its native parent. This suggests that hybridization with native relatives is an important way for invasive species to overcome environmental pressure and achieve invasion.

Genome-wide identification and evolution of the tubulin gene family in Camelina sativa.

BACKGROUND: Tubulins play crucial roles in numerous fundamental processes of plant development. In flowering plants, tubulins are grouped into α-, ß- and γ-subfamilies, while α- and ß-tubulins possess a large isotype diversity and gene number variations among different species. This circumstance leads to insufficient recognition of orthologous isotypes and significantly complicates extrapolation of obtained experimental results, and brings difficulties for the identification of particular tubulin isotype function. The aim of this research is to identify and characterize tubulins of an emerging biofuel crop Camelina sativa. RESULTS: We report comprehensive identification and characterization of tubulin gene family in C. sativa, including analyses of exon-intron organization, duplicated genes comparison, proper isotype designation, phylogenetic analysis, and expression patterns in different tissues. 17 α-, 34 ß- and 6 γ-tubulin genes were identified and assigned to a particular isotype. Recognition of orthologous tubulin isotypes was cross-referred, involving data of phylogeny, synteny analyses and genes allocation on reconstructed genomic blocks of Ancestral Crucifer Karyotype. An investigation of expression patterns of tubulin homeologs revealed the predominant role of N6 (A) and N7 (B) subgenomes in tubulin expression at various developmental stages, contrarily to general the dominance of transcripts of H7 (C) subgenome. CONCLUSIONS: For the first time a complete set of tubulin gene family members was identified and characterized for allohexaploid C. sativa species. The study demonstrates the comprehensive approach of precise inferring gene orthology. The applied technique allowed not only identifying C. sativa tubulin orthologs in model Arabidopsis species and tracking tubulin gene evolution, but also uncovered that A. thaliana is missing orthologs for several particular isotypes of α- and ß-tubulins.

Structure and dynamics of microbial communities associated with the resurrection plant Boea hygrometrica in response to drought stress.

MAIN CONCLUSION: The resurrection plant Boea hygrometrica selectively recruits and assembles drought-specific microbial communities across the plant-soil compartments, which may benefit plant growth and fitness under extreme drought conditions. Plant-associated microbes are essential for facilitating plant growth and fitness under drought stress. The resurrection plant Boea hygrometrica in natural habitats with seasonal rainfall can survive rapid desiccation, yet their interaction with microbiomes under drought conditions remains unexplored. This study examined the bacterial and fungal microbiome structure and drought response across plant-soil compartments of B. hygrometrica by high-throughput amplicon sequencing of 16S rRNA gene and internal transcribed spacer. Our results demonstrated that the diversity, composition, and functional profile of the microbial community varied considerably across the plant-soil compartments and were strongly affected by drought stress. Bacterial and fungal diversity was significantly reduced from soil to endosphere and belowground to aboveground compartments. The compartment-specific enrichment of the dominant bacteria phylum Cyanobacteriota and genus Methylorubrum in leaf endosphere, genera Pseudonocardia in rhizosphere soil and Actinoplanes in root endosphere, and fungal phylum Ascomycota in the aboveground compartments and genera Knufia in root endosphere and Cladosporium in leaf endosphere composed part of the core microbiota with corresponding enrichment of beneficial functions for plant growth and fitness. Moreover, the recruitment of dominant microbial genera Sphingosinicella and Plectosphaerella, Ceratobasidiaceae mycorrhizal fungi, and numerous plant growth-promoting bacteria involving nutrient supply and auxin regulation was observed in desiccated B. hygrometrica plants. Our results suggest that the stable assembled drought-specific microbial community of B. hygrometrica may contribute to plant survival under extreme environments and provide valuable microbial resources for the microbe-mediated drought tolerance enhancement in crops.

Parental conflict driven regulation of endosperm cellularization by a family of Auxin Response Factors.

The endosperm is a reproductive tissue supporting embryo development. In most flowering plants, the initial divisions of endosperm nuclei are not succeeded by cellularization; this process occurs only after a specific number of mitotic cycles have taken place. The timing of cellularization significantly influences seed viability and size. Previous research implicated auxin as a key factor in initiating nuclear divisions and determining the timing of cellularization. Here we uncover the involvement of a family of clustered auxin response factors (cARFs) as dosage-sensitive regulators of endosperm cellularization. cARFs, maternally expressed and paternally silenced, are shown to induce cellularization, thereby restricting seed growth. Our findings align with the predictions of the parental conflict theory, suggesting that cARFs represent major molecular targets in this conflict. We further demonstrate a recurring amplification of cARFs in the Brassicaceae, suggesting an evolutionary response to parental conflict by reinforcing maternal control over endosperm cellularization. Our study highlights that antagonistic parental control on endosperm cellularization converges on auxin biosynthesis and signalling.

Reciprocal conversion between annual and polycarpic perennial flowering behavior in the Brassicaceae.

The development of perennial crops holds great promise for sustainable agriculture and food security. However, the evolution of the transition between perenniality and annuality is poorly understood. Here, using two Brassicaceae species, Crucihimalaya himalaica and Erysimum nevadense, as polycarpic perennial models, we reveal that the transition from polycarpic perennial to biennial and annual flowering behavior is a continuum determined by the dosage of three closely related MADS-box genes. Diversification of the expression patterns, functional strengths, and combinations of these genes endows species with the potential to adopt various life-history strategies. Remarkably, we find that a single gene among these three is sufficient to convert winter-annual or annual Brassicaceae plants into polycarpic perennial flowering plants. Our work delineates a genetic basis for the evolution of diverse life-history strategies in plants and lays the groundwork for the generation of diverse perennial Brassicaceae crops in the future.

Myrosinase isogenes in wasabi (Wasabia japonica Matsum) and their putative roles in glucosinolate metabolism.

BACKGROUND: Wasabi, a Brassicaceae member, is well-known for its unique pungent and hot flavor which is produced from glucosinolate (GSL) degradation. Myrosinase (MYR) is a principle enzyme catalyzing the primary conversion of GSLs to GSL hydrolysis products (GHPs) which is responsible for plant defense system and food quality. Due to the limited information in relation to MYRs present in wasabi (Wasabia japonica M.), this study aimed to identify the MYR isogenes in W. japonica and analyze their roles in relation to GSL metabolism. RESULTS: In results, WjMYRI-1 was abundantly expressed in all organs, whereas WjMYRI-2 showed only trace expression levels. WjMYRII was highly expressed in the aboveground tissues. Interestingly, WjMYRII expression was significantly upregulated by certain abiotic factors, such as methyl jasmonate (more than 40-fold in petioles and 15-fold in leaves) and salt (tenfold in leaves). Young leaves and roots contained 97.89 and 91.17 µmol‧g-1 of GSL, whereas less GSL was produced in mature leaves and petioles (38.36 and 44.79 µmol‧g-1, respectively). Similar pattern was observed in the accumulation of GHPs in various plant organs. Notably, despite the non-significant changes in GSL production, abiotic factors treated samples enhanced significantly GHP content. Pearson's correlation analysis revealed that WjMYRI-1 expression significantly correlated with GSL accumulation and GHP formation, suggesting the primary role of WjMYRI-1-encoding putative protein in GSL degradation. In contrast, WjMYRII expression level showed no correlation with GSL or GHP content, suggesting another physiological role of WjMYRII in stress-induced response. CONCLUSIONS: In conclusions, three potential isogenes (WjMYRI-1, WjMYRI-2, and WjMYRII) encoding for different MYR isoforms in W. japonica were identified. Our results provided new insights related to MYR and GSL metabolism which are important for the implications of wasabi in agriculture, food and pharmaceutical industry. Particularly, WjMYRI-1 may be primarily responsible for GSL degradation, whereas WjMYRII (clade II) may be involved in other regulatory pathways induced by abiotic factors.

Systematic characterization of CsbZIP transcription factors in Camelina sativa and functional analysis of CsbZIP-A12 mediating regulation of unsaturated fatty acid-enriched oil biosynthesis.

The basic leucine zipper (bZIP) transcription factors (TFs) function importantly in numerous life processes in plants. However, bZIP members and their biological roles remain unknown in Camelina sativa, a worldwide promising oil crop. Here, 220 CsbZIP proteins were identified in camelina and classified into thirteen groups. Two and 347 pairs of tandem and segmental duplication genes were detected to be underwent purification selection, with segmental duplication as the main driven-force of CsbZIP gene family expansion. Most CsbZIP genes displayed a tissue-specific expression pattern. Particularly, CsbZIP-A12 significantly positively correlated with many FA/oil biosynthesis-related genes, indicating CsbZIP-A12 may regulate lipid biosynthesis. Notably, yeast one-hybrid (Y1H), ß-Glucuronidase (GUS), dual-luciferase (LUC) and EMSA assays evidenced that CsbZIP-A12 located in nucleus interacted with the promoters of CsSAD2-3 and CsFAD3-3 genes responsible for unsaturated fatty acid (UFA) synthesis, thus activating their transcriptions. Overexpression of CsbZIP-A12 led to an increase of total lipid by 3.275 % compared to the control, followed with oleic and α-linolenic acid levels enhanced by 3.4 % and 5.195 %, and up-regulated the expressions of CsSAD2-3, CsFAD3-3 and CsPDAT2-3 in camelina seeds. Furthermore, heterogeneous expression of CsbZIP-A12 significantly up-regulated the expressions of NtSAD2, NtFAD3 and NtPDAT genes in tobacco plants, thereby improving the levels of total lipids and UFAs in both leaves and seeds without negative effects on other agronomic traits. Together, our findings suggest that CsbZIP-A12 upregulates FA/oil biosynthesis by activating CsSAD2-3 and CsFAD3-3 as well as possible other related genes. These data lay a foundation for further functional analyses of CsbZIPs, providing new insights into the TF-based lipid metabolic engineering to increase vegetable oil yield and health-beneficial quality in oilseeds.

Molecular mechanisms of self-incompatibility in Brassicaceae and Solanaceae.

Self-incompatibility (SI) is a mechanism for preventing self-fertilization in flowering plants. SI is controlled by a single S-locus with multiple haplotypes (S-haplotypes). When the pistil and pollen share the same S-haplotype, the pollen is recognized as self and rejected by the pistil. This review introduces our research on Brassicaceae and Solanaceae SI systems to identify the S-determinants encoded at the S-locus and uncover the mechanisms of self/nonself-discrimination and pollen rejection. The recognition mechanisms of SI systems differ between these families. A self-recognition system is adopted by Brassicaceae, whereas a collaborative nonself-recognition system is used by Solanaceae. Work by our group and subsequent studies indicate that plants have evolved diverse SI systems.

A pollen selection system links self and interspecific incompatibility in the Brassicaceae.

Self-incompatibility and recurrent transitions to self-compatibility have shaped the extant mating systems underlying the nonrandom mating critical for speciation in angiosperms. Linkage between self-incompatibility and speciation is illustrated by the shared pollen rejection pathway between self-incompatibility and interspecific unilateral incompatibility (UI) in the Brassicaceae. However, the pollen discrimination system that activates this shared pathway for heterospecific pollen rejection remains unknown. Here we show that Stigma UI3.1, the genetically identified stigma determinant of UI in Arabidopsis lyrata × Arabidopsis arenosa crosses, encodes the S-locus-related glycoprotein 1 (SLR1). Heterologous expression of A. lyrata or Capsella grandiflora SLR1 confers on some Arabidopsis thaliana accessions the ability to discriminate against heterospecific pollen. Acquisition of this ability also requires a functional S-locus receptor kinase (SRK), whose ligand-induced dimerization activates the self-pollen rejection pathway in the stigma. SLR1 interacts with SRK and interferes with SRK homomer formation. We propose a pollen discrimination system based on competition between basal or ligand-induced SLR1-SRK and SRK-SRK complex formation. The resulting SRK homomer levels would be sensed by the common pollen rejection pathway, allowing discrimination among conspecific self- and cross-pollen as well as heterospecific pollen. Our results establish a mechanistic link at the pollen recognition phase between self-incompatibility and interspecific incompatibility.

Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae.

Proteins containing domain of unknown function (DUF) are prevalent in eukaryotic genome. The DUF1216 proteins possess a conserved DUF1216 domain resembling to the mediator protein of Arabidopsis RNA polymerase II transcriptional subunit-like protein. The DUF1216 family are specifically existed in Brassicaceae, however, no comprehensive evolutionary analysis of DUF1216 genes have been performed. We performed a first comprehensive genome-wide analysis of DUF1216 proteins in Brassicaceae. Totally 284 DUF1216 genes were identified in 27 Brassicaceae species and classified into four subfamilies on the basis of phylogenetic analysis. The analysis of gene structure and conserved motifs revealed that DUF1216 genes within the same subfamily exhibited similar intron/exon patterns and motif composition. The majority members of DUF1216 genes contain a signal peptide in the N-terminal, and the ninth position of the signal peptide in most DUF1216 is cysteine. Synteny analysis revealed that segmental duplication is a major mechanism for expanding of DUF1216 genes in Brassica oleracea, Brassica juncea, Brassica napus, Lepidium meyneii, and Brassica carinata, while in Arabidopsis thaliana and Capsella rubella, tandem duplication plays a major role in the expansion of the DUF1216 gene family. The analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios for DUF1216 paralogous indicated that most of gene pairs underwent purifying selection. DUF1216 genes displayed a specifically high expression in reproductive tissues in most Brassicaceae species, while its expression in Brassica juncea was specifically high in root. Our studies offered new insights into the phylogenetic relationships, gene structures and expressional patterns of DUF1216 members in Brassicaceae, which provides a foundation for future functional analysis.

A High-Throughput Approach for Photosynthesis Studies in a Brassicaceae Panel.

The study of natural variations in photosynthesis in the Brassicaceae family offers the possibility of identifying mechanisms to enhance photosynthetic efficiency in crop plants. Indeed, this family, and particularly its tribe Brassiceae, has been shown to harbor species that have a higher-than-expected photosynthetic efficiency, possibly as a result of a complex evolutionary history. Over the past two decades, methods have been developed to measure photosynthetic efficiency based on chlorophyll fluorescence. Chlorophyll fluorescence measurements are performed with special cameras, such as the FluorCams, which can be included in robotic systems to create high-throughput phenotyping platforms. While these platforms have so far demonstrated high efficiency in measuring small model species like Arabidopsis thaliana, they have the drawback of limited adaptability to accommodate different plant sizes. As a result, the range of species that can be analyzed is restricted. This chapter presents our approach to analyze the photosynthetic parameters: ϕPSII and Fv/Fm for a panel of Brassicaceae species, including a high-photosynthesis species, Hirschfeldia incana, and the adaptations to the phenotyping platform that are required to accommodate this varied group of plants.

The dimorphic diaspore model Aethionema arabicum (Brassicaceae): Distinct molecular and morphological control of responses to parental and germination temperatures.

Plants in habitats with unpredictable conditions often have diversified bet-hedging strategies that ensure fitness over a wider range of variable environmental factors. A striking example is the diaspore (seed and fruit) heteromorphism that evolved to maximize species survival in Aethionema arabicum (Brassicaceae) in which external and endogenous triggers allow the production of two distinct diaspores on the same plant. Using this dimorphic diaspore model, we identified contrasting molecular, biophysical, and ecophysiological mechanisms in the germination responses to different temperatures of the mucilaginous seeds (M+ seed morphs), the dispersed indehiscent fruits (IND fruit morphs), and the bare non-mucilaginous M- seeds obtained by pericarp (fruit coat) removal from IND fruits. Large-scale comparative transcriptome and hormone analyses of M+ seeds, IND fruits, and M- seeds provided comprehensive datasets for their distinct thermal responses. Morph-specific differences in co-expressed gene modules in seeds, as well as in seed and pericarp hormone contents, identified a role of the IND pericarp in imposing coat dormancy by generating hypoxia affecting abscisic acid (ABA) sensitivity. This involved expression of morph-specific transcription factors, hypoxia response, and cell wall remodeling genes, as well as altered ABA metabolism, transport, and signaling. Parental temperature affected ABA contents and ABA-related gene expression and altered IND pericarp biomechanical properties. Elucidating the molecular framework underlying the diaspore heteromorphism can provide insight into developmental responses to globally changing temperatures.

Exo84c-regulated degradation is involved in the normal self-incompatible response in Brassicaceae.

The self-incompatibility system evolves in angiosperms to promote cross-pollination by rejecting self-pollination. Here, we show the involvement of Exo84c in the SI response of both Brassica napus and Arabidopsis. The expression of Exo84c is specifically elevated in stigma during the SI response. Knocking out Exo84c in B. napus and SI Arabidopsis partially breaks down the SI response. The SI response inhibits both the protein secretion in papillae and the recruitment of the exocyst complex to the pollen-pistil contact sites. Interestingly, these processes can be partially restored in exo84c SI Arabidopsis. After incompatible pollination, the turnover of the exocyst-labeled compartment is enhanced in papillae. However, this process is perturbed in exo84c SI Arabidopsis. Taken together, our results suggest that Exo84c regulates the exocyst complex vacuolar degradation during the SI response. This process is likely independent of the known SI pathway in Brassicaceae to secure the SI response.

Overexpression of the Phosphatidylcholine:DiacylglycerolCholinephosphotransferase (PDCT) gene increases carbon flux toward triacylglycerol (TAG) synthesis in Camelinasativa seeds.

Camelinasativa has considerable promise as a dedicated industrial oilseed crop. Its oil-based blends have been tested and approved as liquid transportation fuels. Previously, we utilized metabolomic and transcriptomic profiling approaches and identified metabolic bottlenecks that control oil production and accumulation in seeds. Accordingly, we selected candidate genes for the metabolic engineering of Camelina. Here we targeted the overexpression of Camelina PDCT gene, which encodes the phosphatidylcholine: diacylglycerol cholinephosphotransferase enzyme. PDCT is proposed as a gatekeeper responsible for the interconversions of diacylglycerol (DAG) and phosphatidylcholine (PC) pools and has the potential to increase the levels of TAG in seeds. To confirm whether increased CsPDCT activity in developing Camelina seeds would enhance carbon flux toward increased levels of TAG and alter oil composition, we overexpressed the CsPDCT gene under the control of the seed-specific phaseolin promoter. Camelina transgenics exhibited significant increases in seed yield (19-56%), seed oil content (9-13%), oil yields per plant (32-76%), and altered polyunsaturated fatty acid (PUFA) content compared to their parental wild-type (WT) plants. Results from [14C] acetate labeling of Camelina developing embryos expressing CsPDCT in culture indicated increased rates of radiolabeled fatty acid incorporation into glycerolipids (up to 64%, 59%, and 43% higher in TAG, DAG, and PC, respectively), relative to WT embryos. We conclude that overexpression of PDCT appears to be a positive strategy to achieve a synergistic effect on the flux through the TAG synthesis pathway, thereby further increasing oil yields in Camelina.

The evolution of ephemeral flora in Xinjiang, China: insights from plastid phylogenomic analyses of Brassicaceae.

BACKGROUND: The ephemeral flora of northern Xinjiang, China, plays an important role in the desert ecosystems. However, the evolutionary history of this flora remains unclear. To gain new insights into its origin and evolutionary dynamics, we comprehensively sampled ephemeral plants of Brassicaceae, one of the essential plant groups of the ephemeral flora. RESULTS: We reconstructed a phylogenetic tree using plastid genomes and estimated their divergence times. Our results indicate that ephemeral species began to colonize the arid areas in north Xinjiang during the Early Miocene and there was a greater dispersal of ephemeral species from the surrounding areas into the ephemeral community of north Xinjiang during the Middle and Late Miocene, in contrast to the Early Miocene or Pliocene periods. CONCLUSIONS: Our findings, together with previous studies, suggest that the ephemeral flora originated in the Early Miocene, and species assembly became rapid from the Middle Miocene onwards, possibly attributable to global climate changes and regional geological events.

Evolutionary Analysis of Six Gene Families Part of the Reactive Oxygen Species (ROS) Gene Network in Three Brassicaceae Species.

Climate change is expected to intensify the occurrence of abiotic stress in plants, such as hypoxia and salt stresses, leading to the production of reactive oxygen species (ROS), which need to be effectively managed by various oxido-reductases encoded by the so-called ROS gene network. Here, we studied six oxido-reductases families in three Brassicaceae species, Arabidopsis thaliana as well as Nasturtium officinale and Eutrema salsugineum, which are adapted to hypoxia and salt stress, respectively. Using available and new genomic data, we performed a phylogenomic analysis and compared RNA-seq data to study genomic and transcriptomic adaptations. This comprehensive approach allowed for the gaining of insights into the impact of the adaptation to saline or hypoxia conditions on genome organization (gene gains and losses) and transcriptional regulation. Notably, the comparison of the N. officinale and E. salsugineum genomes to that of A. thaliana highlighted changes in the distribution of ohnologs and homologs, particularly affecting class III peroxidase genes (CIII Prxs). These changes were specific to each gene, to gene families subjected to duplication events and to each species, suggesting distinct evolutionary responses. The analysis of transcriptomic data has allowed for the identification of genes related to stress responses in A. thaliana, and, conversely, to adaptation in N. officinale and E. salsugineum.

Nudix hydrolase WvNUDX24 is involved in borneol biosynthesis in Wurfbainia villosa.

Borneol, camphor, and bornyl acetate are highly promising monoterpenoids widely used in medicine, flavor, food, and chemical applications. Bornyl diphosphate (BPP) serves as a common precursor for the biosynthesis of these monoterpenoids. Although bornyl diphosphate synthase (BPPS) that catalyzes the cyclization of geranyl diphosphate (GPP) to BPP has been identified in multiple plants, the enzyme responsible for the hydrolysis of BPP to produce borneol has not been reported. Here, we conducted in vitro and in vivo functional characterization to identify the Nudix hydrolase WvNUDX24 from W. villosa, which specifically catalyzes the hydrolysis of BPP to generate bornyl phosphate (BP), and then BP forms borneol under the action of phosphatase. Subcellular localization experiments indicated that the hydrolysis of BPP likely occurs in the cytoplasm. Furthermore, site-directed mutagenesis experiments revealed that four critical residues (R84, S96, P98, and G99) for the hydrolysis activity of WvNUDX24. Additionally, the functional identification of phosphatidic acid phosphatase (PAP) demonstrated that WvPAP5 and WvPAP10 were able to hydrolyze geranylgeranyl diphosphate (GGPP) and farnesyl diphosphate (FPP) to generate geranylgeranyl phosphate (GGP) and farnesyl phosphate (FP), respectively, but could not hydrolyze BPP, GPP, and neryl diphosphate (NPP) to produce corresponding monophosphate products. These findings highlight the essential role of WvNUDX24 in the first step of BPP hydrolysis to produce borneol and provide genetic elements for the production of BPP-related terpenoids through plant metabolic engineering and synthetic biology.

Transposable elements contribute to the establishment of the glycine shuttle in Brassicaceae species.

C3 -C4 intermediate photosynthesis has evolved at least five times convergently in the Brassicaceae, despite this family lacking bona fide C4 species. The establishment of this carbon concentrating mechanism is known to require a complex suite of ultrastructural modifications, as well as changes in spatial expression patterns, which are both thought to be underpinned by a reconfiguration of existing gene-regulatory networks. However, to date, the mechanisms which underpin the reconfiguration of these gene networks are largely unknown. In this study, we used a pan-genomic association approach to identify genomic features that could confer differential gene expression towards the C3 -C4 intermediate state by analysing eight C3 species and seven C3 -C4 species from five independent origins in the Brassicaceae. We found a strong correlation between transposable element (TE) insertions in cis-regulatory regions and C3 -C4 intermediacy. Specifically, our study revealed 113 gene models in which the presence of a TE within a gene correlates with C3 -C4 intermediate photosynthesis. In this set, genes involved in the photorespiratory glycine shuttle are enriched, including the glycine decarboxylase P-protein whose expression domain undergoes a spatial shift during the transition to C3 -C4 photosynthesis. When further interrogating this gene, we discovered independent TE insertions in its upstream region which we conclude to be responsible for causing the spatial shift in GLDP1 gene expression. Our findings hint at a pivotal role of TEs in the evolution of C3 -C4 intermediacy, especially in mediating differential spatial gene expression.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae).

KEY MESSAGE: Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plants feature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604 bp and 46.83%, 288,122 bp and 44.79%, and 287,054 bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editing sites. The Mli mitogenome features four conformations likely mediated by the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

Differential expression and evolutionary diversification of RNA helicases in Boechera sexual and apomictic reproduction.

In higher plants, sexual reproduction is characterized by meiosis of the first cells of the germlines, and double fertilization of the egg and central cell after gametogenesis. In contrast, in apomicts of the genus Boechera, meiosis is omitted or altered and only the central cell requires fertilization, while the embryo forms parthenogenetically from the egg cell. To deepen the understanding of the transcriptional basis underlying these differences, we applied RNA-seq to compare expression in reproductive tissues of different Boechera accessions. This confirmed previous evidence of an enrichment of RNA helicases in plant germlines. Furthermore, few RNA helicases were differentially expressed in female reproductive ovule tissues harboring mature gametophytes from apomictic and sexual accessions. For some of these genes, we further found evidence for a complex recent evolutionary history. This included a homolog of Arabidopsis thaliana FASCIATED STEM4 (FAS4). In contrast to AtFAS4, which is a single-copy gene, FAS4 is represented by three homologs in Boechera, suggesting a potential for subfunctionalization to modulate reproductive development. To gain first insights into functional roles of FAS4, we studied Arabidopsis lines carrying mutant alleles. This identified the crucial importance of AtFAS4 for reproduction, as we observed developmental defects and arrest during male and female gametogenesis.