Functional identification and characterization of the Brassica napus transcription factor gene BnAP2, the ortholog of Arabidopsis thaliana APETALA2

Knockout of stigmatic ascorbate peroxidase 1 (APX1) delays pollen rehydration and germination by mediating ROS homeostasis in Brassica napus L

The successful interaction between pollen and stigma is a critical process for plant sexual reproduction, involving a series of intricate molecular and physiological events. After self-compatible pollination, a significant reduction in reactive oxygen species (ROS) production has been observed in stigmas, which is essential for pollen grain rehydration and subsequent pollen tube growth. Several scavenging enzymes tightly regulate ROS homeostasis. However, the potential role of these ROS-scavenging enzymes in the pollen-stigma interaction in Brassica napus remains unclear. Here, we showed that the activity of ascorbate peroxidase (APX), an enzyme that plays a crucial role in the detoxification of hydrogen peroxide (H2O2), was modulated depending on the compatibility of pollination in B. napus. We then identified stigma-expressed APX1s and generated pentuple mutants of APX1s using CRISPR/Cas9 technology. After compatible pollination, the BnaAPX1 pentuple mutants accumulated higher levels of H2O2 in the stigma, while the overexpression of BnaA09.APX1 resulted in lower levels of H2O2. Furthermore, the knockout of BnaAPX1 delayed the compatible response-mediated pollen rehydration and germination, which was consistent with the effects of a specific APX inhibitor, ρ-Aminophenol, on compatible pollination. In contrast, the overexpression of BnaA09.APX1 accelerated pollen rehydration and germination after both compatible and incompatible pollinations. However, delaying and promoting pollen rehydration and germination did not affect the seed set after compatible and incompatible pollination in APX1 pentuple mutants and overexpression lines, respectively. Our results demonstrate the fundamental role of BnaAPX1 in pollen rehydration and germination by regulating ROS homeostasis during the pollen-stigma interaction in B. napus.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Genome-wide analysis of the DREB family genes and functional identification of the involvement of BrDREB2B in abiotic stress in wucai (Brassica campestris L.)

Dehydration responsive element binding protein (DREB) is a significant transcription factor class known to be implicated in abiotic stresses. In this study, we systematically conducted a genome-wide identification and expression analysis of the DREB gene family, including gene structures, evolutionary relationships, chromosome distribution, conserved domains, and expression patterns. A total of 65 DREB family gene members were identified in Chinese cabbage (Brassica rapa L.) and were classified into five subgroups based on phylogenetic analysis. Through analysis of the conserved domains of BrDREB family genes, only one exon existed in the gene structure. Through the analysis of cis-acting elements, these genes were mainly involved in hormone regulation and adversity stress. In order to identify the function of BrDREB2B, overexpressed transgenic Arabidopsis was constructed. After different stress treatments, the germination rate, root growth, survival rate, and various plant physiological indicators were measured. The results showed that transgenic Arabidopsis thaliana plants overexpressing BrDREB2B exhibited enhanced tolerance to salt, heat and drought stresses. Taken together, our results are the first to report the BrDREB2B gene response to drought and heat stresses in Chinese cabbage and provide a basis for further studies to determine the function of BrDREBs in response to abiotic stresses.

SMRT and Illumina RNA-Seq Identifies Potential Candidate Genes Related to the Double Flower Phenotype and Unveils SsAP2 as a Key Regulator of the Double-Flower Trait in Sagittaria sagittifolia

Double flowers are one of the important objectives of ornamental plant breeding. Sagittaria sagittifolia is an aquatic herb in the Alismataceae family that is widely used as an ornamental plant in gardens. However, the reference genome has not been published, and the molecular regulatory mechanism of flower formation remains unclear. In this study, single molecule real-time (SMRT) sequencing technology combined with Illumina RNA-Seq was used to perform a more comprehensive analysis of S. sagittifolia for the first time. We obtained high-quality full-length transcripts, including 53,422 complete open reading frames, and identified 5980 transcription factors that belonged to 67 families, with many MADS-box genes involved in flower formation being obtained. The transcription factors regulated by plant hormone signals played an important role in the development of double flowers. We also identified an AP2 orthologous gene, SsAP2, with a deletion of the binding site for miR172, that overexpressed SsAP2 in S. sagittifolia and exhibited a delayed flowering time and an increased number of petals. This study is the first report of a full-length transcriptome of S. sagittifolia. These reference transcripts will be valuable resources for the analysis of gene structures and sequences, which provide a theoretical basis for the molecular regulatory mechanism governing the formation of double flowers.

Temporal Control of Seed Development in Dicots: Molecular Bases, Ecological Impact and Possible Evolutionary Ramifications

In flowering plants, seeds serve as organs of both propagation and dispersal. The developing seed passes through several consecutive stages, following a conserved general outline. The overall time needed for a seed to develop, however, may vary both within and between plant species, and these temporal developmental properties remain poorly understood. In the present paper, we summarize the existing data for seed development alterations in dicot plants. For genetic mutations, the reported cases were grouped in respect of the key processes distorted in the mutant specimens. Similar phenotypes arising from the environmental influence, either biotic or abiotic, were also considered. Based on these data, we suggest several general trends of timing alterations and how respective mechanisms might add to the ecological plasticity of the families considered. We also propose that the developmental timing alterations may be perceived as an evolutionary substrate for heterochronic events. Given the current lack of plausible models describing timing control in plant seeds, the presented suggestions might provide certain insights for future studies in this field.

Differential expression of SlKLUH controlling fruit and seed weight is associated with changes in lipid metabolism and photosynthesis-related genes

The sizes of plant organs such as fruit and seed are crucial yield components. Tomato KLUH underlies the locus fw3.2, an important regulator of fruit and seed weight. However, the mechanism by which the expression levels of KLUH affect organ size is poorly understood. We found that higher expression of SlKLUH increased cell proliferation in the pericarp within 5 d post-anthesis in tomato near-isogenic lines. Differential gene expression analyses showed that lower expression of SlKLUH was associated with increased expression of genes involved in lipid metabolism. Lipidomic analysis revealed that repression of SlKLUH mainly increased the contents of certain non-phosphorus glycerolipids and phospholipids and decreased the contents of four unknown lipids. Co-expression network analyses revealed that lipid metabolism was possibly associated with but not directly controlled by SlKLUH, and that this gene instead controls photosynthesis-related processes. In addition, many transcription factors putatively involved in the KLUH pathway were identified. Collectively, we show that SlKLUH regulates fruit and seed weight which is associated with altered lipid metabolism. The results expand our understanding of fruit and seed weight regulation and offer a valuable resource for functional studies of candidate genes putatively involved in regulation of organ size in tomato and other crops.

Evolution of the Brassicaceae-specific MS5-Like family and neofunctionalization of the novel MALE STERILITY 5 gene essential for male fertility in Brassica napus

New genes (or lineage-specific genes) can facilitate functional innovations. MALE STERILITY 5 (MS5) in Brassica napus is a fertility-related new gene, which has two wild-type alleles (BnMS5a and BnMS5c ) and two mutant alleles (BnMS5b and BnMS5d ) that could induce male sterility. Here, we studied the history and functional evolution of MS5 homologs in plants by phylogenetic analysis and molecular genetic experiments. We identified 727 MS5 homologs and found that they define a Brassicaceae-specific gene family that has expanded partly via multiple tandem gene duplications and also probably transpositions. The MS5 in B. napus is inherited from a basic diploid ancestor of B. rapa. Molecular genetic experiments indicate that BnMS5a and BnMS5c are functionally distinct in B. napus and that BnMS5d can inhibit BnMS5a in B. napus in a dosage-dependent manner. The BnMS5a protein can move in coordination with meiotic telomeres and interact with the nuclear envelope protein SUN1, with a possible crucial role in meiotic chromosome behavior. In summary, BnMS5 belongs to a Brassicaceae-specific new gene family, and has gained a novel function that is essential for male fertility in B. napus through neofunctionalization that has likely occurred since the origin of B. rapa.

Rolling Circle Amplification (RCA)-Mediated Genome-Wide ihpRNAi Mutant Library Construction in Brassica napus

With the successful completion of genomic sequencing for Brassica napus, identification of novel genes, determination of functions performed by genes, and exploring the molecular mechanisms underlying important agronomic traits were challenged. Mutagenesis-based functional genomics techniques including chemical, physical, and insertional mutagenesis have been used successfully in the functional characterization of genes. However, these techniques had their disadvantages and inherent limitations for allopolyploid Brassica napus, which contained a large number of homologous and redundant genes. Long intron-spliced hairpin RNA (ihpRNA) constructs which contained inverted repeats of the target gene separated by an intron, had been shown to be very effective in triggering RNAi in plants. In the present study, the genome-wide long ihpRNA library of B. napus was constructed with the rolling circle amplification (RCA)-mediated technology. Using the phytoene desaturase (PDS) gene as a target control, it was shown that the RCA-mediated long ihpRNA construct was significantly effective in triggering gene silence in B. napus. Subsequently, the resultant long ihpRNA library was transformed into B. napus to produce corresponding RNAi mutants. Among the obtained transgenic ihpRNA population of B. napus, five ihpRNA lines with observable mutant phenotypes were acquired including alterations in the floral model and the stamen development. The target genes could be quickly identified using specific primers. These results showed that the RCA-mediated ihpRNA construction method was effective for the genome-wide long ihpRNA library of B. napus, therefore providing a platform for study of functional genomics in allopolyploid B. napus.

Suppression of B function by chimeric repressor gene-silencing technology (CRES-T) reduces the petaloid tepal identity in transgenic Lilium sp

Some monocotyledonous plants, including liliaceous, amaryllidaceous and iridaceous ones, produce flowers with petaloid tepals in whorls 1 and 2 organs. For explaining the molecular mechanism of two-layered petaloid tepal development, the modified ABC model has been proposed, in which B class genes are expressed in whorl 1 organs as well as in whorls 2 and 3 organs. We have previously obtained results strongly support the modified ABC model by chimeric repressor gene-silencing technology (CRES-T)-mediated suppression of B function in the liliaceous plant Tricyrtis sp. In the present study, we introduced a CRES-T construct derived from the B class gene of Tricyrtis sp. (TrihDEFa-SRDX) into Lilium sp. in order to examine the effect of suppressing B function on the floral organ identity. Flowers of transgenic plants did not open fully and had pale pink-colored tepals with decreased numbers of papillae on the adaxial side in whorls 1 and 2 compared with those of non-transgenic plants. No apparent morphological alterations were observed in whorls 3 and 4 organs. Both the amount of total anthocyanins and the expression levels of endogenous flavonoid biosynthesis-related genes (LhMYB12, LhbHLH2, LhCHS, LhF3H, LhF3’H, LhDFR and LhANS) decreased in whorls 1 and 2 organs of transgenic plants compared with non-transgenic plants. In addition, the expression levels of endogenous B class genes (LFDEF, LFGLOA and LFGLOB) decreased in transgenic plants and the level was negatively correlated with the degree of morphological alteration. Thus suppression of B function may reduce the identity of petaloid tepals in whorls 1 and 2 of transgenic Lilium sp.

Genome-Wide Identification and Characterization of APETALA2/Ethylene-Responsive Element Binding Factor Superfamily Genes in Soybean Seed Development

Glycine max is one of the most important grain and oil crops, and improvement of seed yield is one of the major objectives in soybean breeding. The AP2/ERF superfamily members are involved in regulating flower and seed development in many species, and therefore play key roles in seed yield. However, it is still unknown that how many AP2/ERF members were presented in the G. max genome and whether these AP2/ERF family members function in flower and seed development in G. max. Here, we identified 380 AP2/ERF superfamily genes in the G. max genome. Phylogenetic analysis showed that 323 members were grouped into the ERF family, and 49 into the AP2 family. Among the AP2 family, 14 members of the euAP2 lineage showed high identity with their orthologs, and eight member of the ANT lineage were expressed highly in the seeds. Furthermore, seven of them (GmAP2-1 to GmAP2-7) were successfully cloned and over-expressed in Arabidopsis thaliana. The transgenic Arabidopsis plants over-expressing these GmAP2 genes flowered earlier relative to the wild type control. The seed length and width, and seed area of these over-expression lines were increased compared with the wild type, and seed weight of over-expression lines of GmAP2-1, GmAP2-4, GmAP2-5, and GmAP2-6 were greater than those of the wild type. Furthermore, the seed number per silique of the over-expression lines for GmAP2 genes were not affected except GmAP2-5. Collectively, GmAP2-1, GmAP2-4, and GmAP2-6 played important roles in regulating seed weight by affecting seed length, width and area, and further controlling seed yield.

Defective APETALA2 Genes Lead to Sepal Modification in Brassica Crops

Many vegetable and oilseed crops belong to Brassica species. The seed production of these crops is hampered often by abnormal floral organs, especially under the conditions of abiotic conditions. However, the molecular reasons for these abnormal floral organs remains poorly understood. Here, we report a novel pistil-like flower mutant of B. rapa. In the flower of this mutant, the four sepals are modified to one merged carpel that look like a ring in the sepal positions, enveloping some abnormal stamens and a pistil, and resulting in poor seed production. This novel mutant is named sepal-carpel modification (scm). DNA sequencing showed that the BrAP2a gene, the ortholog of Arabidopsis APETALA2 (AP2) that specifies sepal identity, losses the function of in scm mutant due to a 119-bp repeated sequence insertion that resulted in an early transcription termination. BrAP2b, the paralog of BrAP2a featured two single-nucleotide substitutions that cause a single amino acid substitution in the highly conserved acidic serine-rich transcriptional activation domain. Each of the two BrAP2 genes rescues the sepal defective phenotype of the ap2-5 mutant of Arabidopsis. Furthermore, the knockout mutation of the corresponding BnAP2 genes of oilseed rape (B. napus) by CRISPR/Cas9-mediated genome editing system resulted in scm-like phenotype. These results suggest that BrAP2 gene plays a key role in sepal modification. Our finding provides an insight into molecular mechanism underlying morphological modification of floral organs and is useful for genetic manipulation of flower modification and improvement of seed production of Brassica crops.

Induction of telomere-mediated chromosomal truncation and behavior of truncated chromosomes in Brassica napus

Engineered minichromosomes could be stably inherited and serve as a platform for simultaneously transferring and stably expressing multiple genes. Chromosomal truncation mediated by repeats of telomeric sequences is a promising approach for the generation of minichromosomes. In the present work, direct repetitive sequences of Arabidopsis telomere were used to study telomere-mediated truncation of chromosomes in Brassica napus. Transgenes containing alien Arabidopsis telomere were successfully obtained, and Southern blotting and fluorescence in situ hybridization (FISH) results show that the transgenes resulted in successful chromosomal truncation in B. napus. In addition, truncated chromosomes were inherited at rates lower than that predicted by Mendelian rules. To determine the potential manipulations and applications of the engineered chromosomes, such as the stacking of multiple transgenes and the Cre/lox and FRT/FLP recombination systems, both amenable to genetic manipulations through site-specific recombination in somatic cells, were tested for their ability to undergo recombination in B. napus. These results demonstrate that alien Arabidopsis telomere is able to mediate chromosomal truncation in B. napus. This technology would be feasible for chromosomal engineering and for studies on chromosome structure and function in B. napus.

High Density Linkage Map Construction and QTL Detection for Three Silique-Related Traits in Orychophragmus violaceus Derived Brassica napus Population

Seeds per silique (SS), seed weight (SW), and silique length (SL) are important determinant traits of seed yield potential in rapeseed (Brassica napus L.), and are controlled by naturally occurring quantitative trait loci (QTLs). Mapping QTLs to narrow chromosomal regions provides an effective means of characterizing the genetic basis of these complex traits. Orychophragmus violaceus is a crucifer with long siliques, many SS, and heavy seeds. A novel B. napus introgression line with many SS was previously selected from multiple crosses (B. rapa ssp. chinesis × O. violaceus) × B. napus. In present study, a doubled haploid (DH) population with 167 lines was established from a cross between the introgression line and a line with far fewer SS, in order to detect QTLs for silique-related traits. By screening with a Brassica 60K single nucleotide polymorphism (SNP) array, a high-density linkage map consisting of 1,153 bins and spanning a cumulative length of 2,209.1 cM was constructed, using 12,602 high-quality polymorphic SNPs in the DH population. The average recombination bin densities of the A and C subgenomes were 1.7 and 2.4 cM, respectively. 45 QTLs were identified for the three traits in all, which explained 4.0-34.4% of the total phenotypic variation; 20 of them were integrated into three unique QTLs by meta-analysis. These unique QTLs revealed a significant positive correlation between SS and SL and a significant negative correlation between SW and SS, and were mapped onto the linkage groups A05, C08, and C09. A trait-by-trait meta-analysis revealed eight, four, and seven consensus QTLs for SS, SW, and SL, respectively, and five major QTLs (cqSS.A09b, cqSS.C09, cqSW.A05, cqSW.C09, and cqSL.C09) were identified. Five, three, and four QTLs for SS, SW, and SL, respectively, might be novel QTLs because of the existence of alien genetic loci for these traits in the alien introgression. Thirty-eight candidate genes underlying nine QTLs for silique-related traits were identified.

AP2/ERF Transcription Factor, Ii049, Positively Regulates Lignan Biosynthesis in Isatis indigotica through Activating Salicylic Acid Signaling and Lignan/Lignin Pathway Genes

Lignans, such as lariciresinol and its derivatives, have been identified as effective antiviral ingredients in Isatis indigotica. Evidence suggests that the APETALA2/ethylene response factor (AP2/ERF) family might be related to the biosynthesis of lignans in I. indigotica. However, the special role played by the AP2/ERF family in the metabolism and its underlying putative mechanism still need to be elucidated. One novel AP2/ERF gene, named Ii049, was isolated and characterized from I. indigotica in this study. The quantitative real-time PCR analysis revealed that Ii049 was expressed highest in the root and responded to methyl jasmonate, salicylic acid (SA) and abscisic acid treatments to various degrees. Subcellular localization analysis indicated that Ii049 protein was localized in the nucleus. Knocking-down the expression of Ii049 caused a remarkable reduction of lignan/lignin contents and transcript levels of genes involved in the lignan/lignin biosynthetic pathway. Ii049 bound to the coupled element 1, RAV1AAT and CRTAREHVCBF2 motifs of genes IiPAL and IiCCR, the key structural genes in the lignan/lignin pathway. Furthermore, Ii049 was also essential for SA biosynthesis, and SA induced lignan accumulation in I. indigotica. Notably, the transgenic I. indigotica hairy roots overexpressing Ii049 showed high expression levels of lignan/lignin biosynthetic genes and SA content, resulting in significant accumulation of lignan/lignin. The best-engineered line (OVX049-10) produced 425.60 μg·g^-1 lariciresinol, an 8.3-fold increase compared with the wild type production. This study revealed the function of Ii049 in regulating lignan/lignin biosynthesis, which had the potential to increase the content of valuable lignan/lignin in economically significant medicinal plants.

Suppression of B function strongly supports the modified ABCE model in Tricyrtis sp. (Liliaceae)

B class MADS-box genes play important roles in petal and stamen development. Some monocotyledonous species, including liliaceous ones, produce flowers with petaloid tepals in whorls 1 and 2. A modified ABCE model has been proposed to explain the molecular mechanism of development of two-layered petaloid tepals. However, direct evidence for this modified ABCE model has not been reported to date. To clarify the molecular mechanism determining the organ identity of two-layered petaloid tepals, we used chimeric repressor gene-silencing technology (CRES-T) to examine the suppression of B function in the liliaceous ornamental Tricyrtis sp. Transgenic plants with suppressed B class genes produced sepaloid tepals in whorls 1 and 2 instead of the petaloid tepals as expected. In addition, the stamens of transgenic plants converted into pistil-like organs with ovule- and stigma-like structures. This report is the first to describe the successful suppression of B function in monocotyledonous species with two-layered petaloid tepals, and the results strongly support the modified ABCE model.

AfAP2-1, An Age-Dependent Gene of Aechmea fasciata, Responds to Exogenous Ethylene Treatment

The Bromeliaceae family is one of the most morphologically diverse families with a pantropical distribution. To schedule an appropriate flowering time for bromeliads, ethylene is commonly used to initiate flower development in adult plants. However, the mechanism by which ethylene induces flowering in adult bromeliads remains unknown. Here, we identified an APETALA2 (AP2)-like gene, AfAP2-1, in Aechmea fasciata. AfAP2-1 contains two AP2 domains and is a nuclear-localized protein. It functions as a transcriptional activator, and the activation domain is located in the C-terminal region. The expression level of AfAP2-1 is higher in juvenile plants than in adult plants, and the AfAP2-1 transcript level was rapidly and transiently reduced in plants treated with exogenous ethylene. Overexpression of AfAP2-1 in Arabidopsis thaliana results in an extremely delayed flowering phenotype. These results suggested that AfAP2-1 responds to ethylene and is a putative age-dependent flowering regulator in A. fasciata.

Origination, Expansion, Evolutionary Trajectory, and Expression Bias of AP2/ERF Superfamily in Brassica napus

The AP2/ERF superfamily, one of the most important transcription factor families, plays crucial roles in response to biotic and abiotic stresses. So far, a comprehensive evolutionary inference of its origination and expansion has not been available. Here, we identified 515 AP2/ERF genes in B. napus, a neo-tetraploid forming ~7500 years ago, and found that 82.14% of them were duplicated in the tetraploidization. A prominent subgenome bias was revealed in gene expression, tissue-specific, and gene conversion. Moreover, a large-scale analysis across plants and alga suggested that this superfamily could have been originated from AP2 family, expanding to form other families (ERF, and RAV). This process was accompanied by duplicating and/or alternative deleting AP2 domain, intragenic domain sequence conversion, and/or by acquiring other domains, resulting in copy number variations, alternatively contributing to functional innovation. We found that significant positive selection occurred at certain critical nodes during the evolution of land plants, possibly responding to changing environment. In conclusion, the present research revealed origination, functional innovation, and evolutionary trajectory of the AP2/ERF superfamily, contributing to understanding their roles in plant stress tolerance.

Expression of Brassica napus TTG2, a regulator of trichome development, increases plant sensitivity to salt stress by suppressing the expression of auxin biosynthesis genes

WRKY transcription factors (TFs) are plant specific and play important roles in regulating diverse biological processes. To identify TFs with broad-spectrum effects on various stress responses in Brassica napus, an important oil crop grown across diverse ecological regions worldwide, we functionally characterized Bna.TTG2 genes, which are homologous to the Arabidopsis AtTTG2 (WRKY44) gene. Four Bna.TTG2 genes were capable of rescuing the trichome phenotypes of Arabidopsis ttg2 mutants. Overexpressing one Bna.TTG2 family member, BnaA.TTG2.a.1, remarkably increased trichome numbers in Arabidopsis and B. napus plants. Interestingly, the BnaA.TTG2.a.1-overexpressing plants of both species exhibited increased sensitivity to salt stress. In BnaA.TTG2.a.1-overexpressing Arabidopsis under salt stress, the endogenous indole-3-acetic acid (IAA) content was reduced, and the expression of two auxin biosynthesis genes, TRYPTOPHAN BIOSYNTHESIS 5 (TRP5) and YUCCA2 (YUC2), was downregulated. The results from yeast one-hybrid, electrophoretic mobility shift, and dual-luciferase reporter assays revealed that BnaA.TTG2.a.1 is able to bind to the promoters of TRP5 and YUC2. These data indicated that BnaA.TTG2.a.1 confers salt sensitivity to overexpressing plants by suppressing the expression of IAA synthesis genes and thus lowering IAA levels. Transgenic Arabidopsis plants with an N-terminus-deleted BnaA.TTG2.a.1 no longer showed hypersensitivity to salt stress, suggesting that the N terminus of BnaA.TTG2.a.1 plays a critical role in salt stress responses. Therefore, in addition to its classical function in trichome development, our study reveals a novel role for Bna.TTG2 genes in salt stress responses.

Gene-to-metabolite network for biosynthesis of lignans in MeJA-elicited Isatis indigotica hairy root cultures

Root and leaf tissue of Isatis indigotica shows notable anti-viral efficacy, and are widely used as "Banlangen" and "Daqingye" in traditional Chinese medicine. The plants' pharmacological activity is attributed to phenylpropanoids, especially a group of lignan metabolites. However, the biosynthesis of lignans in I. indigotica remains opaque. This study describes the discovery and analysis of biosynthetic genes and AP2/ERF-type transcription factors involved in lignan biosynthesis in I. indigotica. MeJA treatment revealed differential expression of three genes involved in phenylpropanoid backbone biosynthesis (IiPAL, IiC4H, Ii4CL), five genes involved in lignan biosynthesis (IiCAD, IiC3H, IiCCR, IiDIR, and IiPLR), and 112 putative AP2/ERF transcription factors. In addition, four intermediates of lariciresinol biosynthesis were found to be induced. Based on these results, a canonical correlation analysis using Pearson's correlation coefficient was performed to construct gene-to-metabolite networks and identify putative key genes and rate-limiting reactions in lignan biosynthesis. Over-expression of IiC3H, identified as a key pathway gene, was used for metabolic engineering of I. indigotica hairy roots, and resulted in an increase in lariciresinol production. These findings illustrate the utility of canonical correlation analysis for the discovery and metabolic engineering of key metabolic genes in plants.

RNA-sequencing reveals early, dynamic transcriptome changes in the corollas of pollinated petunias

BACKGROUND

Pollination reduces flower longevity in many angiosperms by accelerating corolla senescence. This response requires hormone signaling between the floral organs and results in the degradation of macromolecules and organelles within the petals to allow for nutrient remobilization to developing seeds. To investigate early pollination-induced changes in petal gene expression, we utilized high-throughput sequencing to identify transcripts that were differentially expressed between corollas of pollinated Petunia × hybrida flowers and their unpollinated controls at 12, 18, and 24 hours after opening.

RESULTS

In total, close to 0.5 billion Illumina 101 bp reads were generated, de novo assembled, and annotated, resulting in an EST library of approximately 33 K genes. Over 4,700 unique, differentially expressed genes were identified using comparisons between the pollinated and unpollinated libraries followed by pairwise comparisons of pollinated libraries to unpollinated libraries from the same time point (i.e. 12-P/U, 18-P/U, and 24-P/U) in the Bioconductor R package DESeq2. Over 500 gene ontology terms were enriched. The response to auxin stimulus and response to 1-aminocyclopropane-1-carboxylic acid terms were enriched by 12 hours after pollination (hap). Using weighted gene correlation network analysis (WGCNA), three pollination-specific modules were identified. Module I had increased expression across pollinated corollas at 12, 18, and 24 h, and modules II and III had a peak of expression in pollinated corollas at 18 h. A total of 15 enriched KEGG pathways were identified. Many of the genes from these pathways were involved in metabolic processes or signaling. More than 300 differentially expressed transcription factors were identified.

CONCLUSIONS

Gene expression changes in corollas were detected within 12 hap, well before fertilization and corolla wilting or ethylene evolution. Significant changes in gene expression occurred at 18 hap, including the up-regulation of autophagy and down-regulation of ribosomal genes and genes involved in carbon fixation. This transcriptomic database will greatly expand the genetic resources available in petunia. Additionally, it will guide future research aimed at identifying the best targets for increasing flower longevity by delaying corolla senescence.

Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis)

Background

Chinese cabbage (Brassica rapa ssp. pekinensis) is a member of one of the most important leaf vegetables grown worldwide, which has experienced thousands of years in cultivation and artificial selection. The entire Chinese cabbage genome sequence, and more than forty thousand proteins have been obtained to date. The genome has undergone triplication events since its divergence from Arabidopsis thaliana (13 to 17 Mya), however a high degree of sequence similarity and conserved genome structure remain between the two species. Arabidopsis is therefore a viable reference species for comparative genomics studies. Variation in the number of members in gene families due to genome triplication may contribute to the broad range of phenotypic plasticity, and increased tolerance to environmental extremes observed in Brassica species. Transcription factors are important regulators involved in plant developmental and physiological processes. The AP2/ERF proteins, one of the most important families of transcriptional regulators, play a crucial role in plant growth, and in response to biotic and abiotic stressors. Our analysis will provide resources for understanding the tolerance mechanisms in Brassica rapa ssp. pekinensis.

Results

In the present study, 291 putative AP2/ERF transcription factor proteins were identified from the Chinese cabbage genome database, and compared with proteins from 15 additional species. The Chinese cabbage AP2/ERF superfamily was classified into four families, including AP2, ERF, RAV, and Soloist. The ERF family was further divided into DREB and ERF subfamilies. The AP2/ERF superfamily was subsequently divided into 15 groups. The identification, classification, phylogenetic reconstruction, conserved motifs, chromosome distribution, functional annotation, expression patterns, and interaction networks of the AP2/ERF transcription factor superfamily were predicted and analyzed. Distribution mapping results showed AP2/ERF superfamily genes were localized on the 10 Chinese cabbage chromosomes. AP2/ERF transcription factor expression levels exhibited differences among six tissue types based on expressed sequence tags (ESTs). In the AP2/ERF superfamily, 214 orthologous genes were identified between Chinese cabbage and Arabidopsis. Orthologous gene interaction networks were constructed, and included seven CBF and four AP2 genes, primarily involved in cold regulatory pathways and ovule development, respectively.

Conclusions

The evolution of the AP2/ERF transcription factor superfamily in Chinese cabbage resulted from genome triplication and tandem duplications. A comprehensive analysis of the physiological functions and biological roles of AP2/ERF superfamily genes in Chinese cabbage is required to fully elucidate AP2/ERF, which provides us with rich resources and opportunities to understand crop stress tolerance mechanisms.