The auxin response factor gene family in allopolyploid Brassica napus

Identification of the High-Affinity Potassium Transporter Gene Family (HKT) in Brassica U-Triangle Species and Its Potential Roles in Abiotic Stress in Brassica napus L

Members of the high-affinity potassium transporter (HKT) protein family regulate the uptake and homeostasis of sodium and potassium ions, but little research describes their roles in response to abiotic stresses in rapeseed (Brassica napus L.). In this study, we identified and characterized a total of 36 HKT genes from the species comprising the triangle of U model (U-triangle species): B. rapa, B. nigra, B. oleracea, B. juncea, B. napus, and B. carinata. We analyzed the phylogenetic relationships, gene structures, motif compositions, and chromosomal distributions of the HKT family members of rapeseed. Based on their phylogenetic relationships and assemblage of functional domains, we classified the HKT members into four subgroups, HKT1;1 to HKT1;4. Analysis of the nonsynonymous substitutions (Ka), synonymous substitutions (Ks), and the Ka/Ks ratios of HKT gene pairs suggested that these genes have experienced strong purifying selective pressure after duplication, with their evolutionary relationships supporting the U-triangle theory. Furthermore, the expression profiles of BnaHKT genes varies among potassium, phytohormone and heavy-metal treatment. Their repression provides resistance to heavy-metal stress, possibly by limiting uptake. Our results systematically reveal the characteristics of HKT family proteins and their encoding genes in six Brassica species and lay a foundation for further exploration of the role of HKT family genes in heavy-metal tolerance.

Mdm-miR160-MdARF17-MdWRKY33 module mediates freezing tolerance in apple

Apple (Malus domestica) trees are vulnerable to freezing temperatures. Cold resistance in woody perennial plants can be improved through biotechnological approaches. However, genetic engineering requires a thorough understanding of the molecular mechanisms of the tree's response to cold. In this study, we demonstrated that the Mdm-miR160-MdARF17-MdWRKY33 module is crucial for apple freezing tolerance. Mdm-miR160 plays a negative role in apple freezing tolerance, whereas MdARF17, one of the targets of Mdm-miR160, is a positive regulator of apple freezing tolerance. RNA sequencing analysis revealed that in apple, MdARF17 mediates the cold response by influencing the expression of cold-responsive genes. EMSA and ChIP-qPCR assays demonstrated that MdARF17 can bind to the promoter of MdWRKY33 and promotes its expression. Overexpression of MdWRKY33 enhanced the cold tolerance of the apple calli. In addition, we found that the Mdm-miR160-MdARF17-MdWRKY33 module regulates cold tolerance in apple by regulating reactive oxygen species (ROS) scavenging, as revealed by (i) increased H₂ O₂ levels and decreased peroxidase (POD) and catalase (CAT) activities in Mdm-miR160e OE plants and MdARF17 RNAi plants and (ii) decreased H₂ O₂ levels and increased POD and CAT activities in MdmARF17 OE plants and MdWRKY33 OE calli. Taken together, our study uncovered the molecular roles of the Mdm-miR160-MdARF17-MdWRKY33 module in freezing tolerance in apple, thus providing support for breeding of cold-tolerant apple cultivars.

Genome-Wide Analysis of the Molecular Functions of B3 Superfamily in Oil Biosynthesis in Olive (Olea europaea L.)

The plant B3 gene superfamily contains a large number of transcription factors playing a vital role in both vegetative growth and reproductive development in plants. Although several B3 genes have been well studied, molecular functions of the B3 genes in olive are largely unknown. In our study, a total of 200 B3 genes were identified in olive genome based on RNA-seq and comparative genomic analyses and further classified into five groups (i.e., REM, RAV, LAV, HSI, and ARF) based on phylogenetic analysis. Results of gene structure and motif composition analyses revealed diversified functions among these five groups of B3 genes. Results of genomic duplication and syntenic analyses indicated the gene expansion in the B3 genes. Results of gene expression based on both transcriptomics and relative expression revealed the tissue-biased expression patterns in B3 genes. The results of the comparative expression analysis of B3 genes between two olive cultivars with high and low oil contents identified several potential REM genes which may be involved in oil biosynthesis in olive. Based on the comprehensive characterization of the molecular structures and functions of B3 genes in olive genome, our study provided novel insights into the potential roles of B3 transcription factors in oil biosynthesis in olive and lays the groundwork for the functional explorations into this research field.

Genome-Wide Identification and Characterization of Auxin Response Factor (ARF) Gene Family Involved in Wood Formation and Response to Exogenous Hormone Treatment in Populus trichocarpa

Auxin is a key regulator that virtually controls almost every aspect of plant growth and development throughout its life cycle. As the major components of auxin signaling, auxin response factors (ARFs) play crucial roles in various processes of plant growth and development. In this study, a total of 35 PtrARF genes were identified, and their phylogenetic relationships, chromosomal locations, synteny relationships, exon/intron structures, cis-elements, conserved motifs, and protein characteristics were systemically investigated. We also analyzed the expression patterns of these PtrARF genes and revealed that 16 of them, including PtrARF1, 3, 7, 11, 13-17, 21, 23, 26, 27, 29, 31, and 33, were preferentially expressed in primary stems, while 15 of them, including PtrARF2, 4, 6, 9, 10, 12, 18-20, 22, 24, 25, 28, 32, and 35, participated in different phases of wood formation. In addition, some PtrARF genes, with at least one cis-element related to indole-3-acetic acid (IAA) or abscisic acid (ABA) response, responded differently to exogenous IAA and ABA treatment, respectively. Three PtrARF proteins, namely PtrARF18, PtrARF23, and PtrARF29, selected from three classes, were characterized, and only PtrARF18 was a transcriptional self-activator localized in the nucleus. Moreover, Y2H and bimolecular fluorescence complementation (BiFC) assay demonstrated that PtrARF23 interacted with PtrIAA10 and PtrIAA28 in the nucleus, while PtrARF29 interacted with PtrIAA28 in the nucleus. Our results provided comprehensive information regarding the PtrARF gene family, which will lay some foundation for future research about PtrARF genes in tree development and growth, especially the wood formation, in response to cellular signaling and environmental cues.

Cross-kingdom transcriptomic trends in the evolution of hybrid gene expression

Hybridization is a route to speciation that occurs widely across the eukaryote tree of life. The success of allopolyploids (hybrid species with increased ploidy) and homoploid hybrids (with unchanged ploidy) is well documented. However, their formation and establishment is not straightforward, with a suite of near‐instantaneous and longer term biological repercussions faced by the new species. Central to these challenges is the rewiring of gene regulatory networks following the merger of distinct genomes inherited from both parental species. Research on the evolution of hybrid gene expression has largely involved studies on a single hybrid species or a few gene families. Here, we present the first standardized transcriptome‐wide study exploring the fates of genes following hybridization across three kingdoms: animals, plants and fungi. Within each kingdom, we pair an allopolyploid system with a closely related homoploid hybrid to decouple the influence of increased ploidy from genome merger. Genome merger, not changes in ploidy, has the greatest effect on posthybridization expression patterns across all study systems. Strikingly, we find that differentially expressed genes in parent species preferentially switch to more similar expression in hybrids across all kingdoms, likely as a consequence of regulatory trans‐acting cross‐talk within the hybrid nucleus. We also highlight the prevalence of gene loss or silencing among extremely differentially expressed genes in hybrid species across all kingdoms. These shared patterns suggest that the evolutionary process of hybridization leads to common high‐level expression outcomes, regardless of the particular species or kingdom.

Genome-Wide Characterization of High-Affinity Nitrate Transporter 2 (NRT2) Gene Family in Brassica napus

Nitrate transporter 2 (NRT2) plays an essential role in Nitrogen (N) uptake, transport, utilization, and stress resistance. In this study, the NRT2 gene family in two sequenced Brassica napus ecotypes were identified, including 31 genes in ‘Zhongshuang11’ (BnaZSNRT2s) and 19 in ‘Darmor-bzh’ (BnaDarNRT2s). The candidate genes were divided into three groups (Group I−III) based on phylogenetic analyses, supported by a conserved intron-exon structure in each group. Collinearity analysis revealed that the large expansion of BnaZSNRT2s attributed to allopolyploidization of ancestors Brassica rapa and Brassica oleracea, and small-scale duplication events in B. napus. Transcription factor (TF) binding site prediction, cis-element analysis, and microRNA prediction suggested that the expressions of BnaZSNRT2s are regulated by multiple factors, and the regulatory pattern is relatively conserved in each group and is tightly connected between groups. Expression assay showed the diverse and differentiated spatial-temporal expression profiles of BnaZSNRT2s in Group I, but conserved patterns were observed in Group II/III; and the low nitrogen (LN) stress up-regulated expression profiles were presented in Group I−III, based on RNA-seq data. RT-qPCR analyses confirmed that BnaZSNRT2.5A-1 and BnaZSNRT2.5C-1 in Group II were highly up-regulated under LN stress in B. napus roots. Our results offer valid information and candidates for further functional BnaZSNRT2s studies.

Genome-Wide Identification of ARF Gene Family Suggests a Functional Expression Pattern during Fruitlet Abscission in Prunus avium L

Auxin response factors (ARFs) play a vital role in plant growth and development. In the current study, 16 ARF members have been identified in the sweet cherry (Prunus avium L.) genome. These genes are all located in the nucleus. Sequence analysis showed that genes in the same subgroup have similar exon-intron structures. A phylogenetic tree has been divided into five groups. The promoter sequence includes six kinds of plant hormone-related elements, as well as abiotic stress response elements such as low temperature or drought. The expression patterns of PavARF in different tissues, fruitlet abscission, cold and drought treatment were comprehensively analyzed. PavARF10/13 was up-regulated and PavARF4/7/11/12/15 was down-regulated in fruitlet abscising. These genes may be involved in the regulation of fruit drop in sweet cherry fruits. This study comprehensively analyzed the bioinformatics and expression pattern of PavARF, which can lay the foundation for further understanding the PavARF family in plant growth development and fruit abscission.

Genome-wide identification and expression analysis of auxin response factors in peanut (Arachis hypogaea L.)

Auxin response factors (ARFs) are transcription factors that regulate the expression of auxin response genes, and have important functions in plant growth and development. In this study, available genome data for peanut (Arachis hypogaea L.) were used to identify AhARF genes. In total, 61 AhARFs and 23 AtARFs were divided into six groups (I-VI). Molecular structural analysis revealed that the protein members of AhARF contain at least two domains, the B3 domain and the Auxin-resp domain, and that some have a C-terminal dimerisation domain. Screening of the transcriptome data of 22 tissues of A. hypogaea cv. Tifrunner in a public database showed high expression levels of AhARF2 and AhARF6. AhARF6 was expressed more highly in the stem and branch than in the root and leaf of the wild species Arachis monticola (A. mon) and cultivated species H103. After treatment with exogenous auxin (NAA), the expression of AhARF6 was inhibited, and this inhibition was greater in A. mon than in H103. The transcriptome map revealed that the expression of AhARF6 was higher in the larger pods of H8107 and ZP06 than in the medium pods of H103 and small pods of A. mon. Moreover, AhARF6-5 was proven to be localised in the nucleus, consistent with the location of AtARF6. These results suggest that AhARF6 may play an important role in pod development in peanut.

Genome-wide analysis of ARF transcription factors reveals HcARF5 expression profile associated with the biosynthesis of β-ocimene synthase in Hedychium coronarium

Herein, 37 ARF genes were identified and analyzed in Hedychium coronarium and HcARF5 showed a potential role in the regulation of HcTPS3. Auxin is an important plant hormone, implicated in various aspects of plant growth and development processes especially in the biosynthesis of various secondary metabolites. Auxin response factors (ARF) belong to the transcription factors (TFs) gene family and play a crucial role in transcriptional activation/repression of auxin-responsive genes by directly binding to their promoter region. Nevertheless, whether ARF genes are involved in the regulatory mechanism of volatile compounds in flowering plants is largely unknown. β-ocimene is a key floral volatile compound synthesized by terpene synthase 3 (HcTPS3) in Hedychium coronarium. A comprehensive analysis of H. coronarium genome reveals 37 candidate ARF genes in the whole genome. Tissue-specific expression patterns of HcARFs family members were assessed using available transcriptome data. Among them, HcARF5 showed a higher expression level in flowers, and significantly correlated with the key structural β-ocimene synthesis gene (HcTPS3). Furthermore, transcript levels of both genes were associated with the flower development. Under hormone treatments, the response of HcARF5 and HcTPS3, and the emission level of β-ocimene contents were evaluated. Subcellular and transcriptional activity assay showed that HcARF5 localizes to the nucleus and possesses transcriptional activity. Yeast one-hybrid (Y1H) and dual-luciferase assays revealed that HcARF5 directly regulates the transcriptional activity of HcTPS3. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that HcARF5 interacts with scent-related HcIAA4, HcIAA6, and HcMYB1 in vivo. Overall, these results indicate that HcARF5 is potentially involved in the regulation of β-ocimene synthesis in H. coronarium.

Genome-wide analysis of ARF transcription factors reveals HcARF5 expression profile associated with the biosynthesis of β-ocimene synthase in Hedychium coronarium

Herein, 37 ARF genes were identified and analyzed in Hedychium coronarium and HcARF5 showed a potential role in the regulation of HcTPS3. Auxin is an important plant hormone, implicated in various aspects of plant growth and development processes especially in the biosynthesis of various secondary metabolites. Auxin response factors (ARF) belong to the transcription factors (TFs) gene family and play a crucial role in transcriptional activation/repression of auxin-responsive genes by directly binding to their promoter region. Nevertheless, whether ARF genes are involved in the regulatory mechanism of volatile compounds in flowering plants is largely unknown. β-ocimene is a key floral volatile compound synthesized by terpene synthase 3 (HcTPS3) in Hedychium coronarium. A comprehensive analysis of H. coronarium genome reveals 37 candidate ARF genes in the whole genome. Tissue-specific expression patterns of HcARFs family members were assessed using available transcriptome data. Among them, HcARF5 showed a higher expression level in flowers, and significantly correlated with the key structural β-ocimene synthesis gene (HcTPS3). Furthermore, transcript levels of both genes were associated with the flower development. Under hormone treatments, the response of HcARF5 and HcTPS3, and the emission level of β-ocimene contents were evaluated. Subcellular and transcriptional activity assay showed that HcARF5 localizes to the nucleus and possesses transcriptional activity. Yeast one-hybrid (Y1H) and dual-luciferase assays revealed that HcARF5 directly regulates the transcriptional activity of HcTPS3. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that HcARF5 interacts with scent-related HcIAA4, HcIAA6, and HcMYB1 in vivo. Overall, these results indicate that HcARF5 is potentially involved in the regulation of β-ocimene synthesis in H. coronarium.

Genome-Wide Identification and Characterization of the IGT Gene Family in Allotetraploid Rapeseed (Brassica napus L.)

IGT family genes function critically to regulate lateral organ orientation in plants. However, little information is available about this family of genes in Brassica napus. In this study, 27 BnIGT genes were identified on 16 chromosomes and divided into seven clades, namely LAZY1∼LAZY6 and TAC1 (Tiller Angle Control 1), based on their phylogenetic relationships. Duplication analysis revealed that 91.1% of the gene pairs were derived from whole-genome duplication. Most BnIGT genes had a similar structural pattern with one or two very short exons followed by a long and a shorter exon. Common and specific motifs were identified among the seven clades, and motif 1, containing the family-specific GφL(A/T)IGT sequence, was observed in all clades except LAZY5. Three types of cis-elements pertinent to transcription factor binding, light responses, and hormone signaling were detected in the BnIGT promoters. Intriguingly, more than half of the BnIGT genes exhibited no or very low expression in various tissues, and the LAZY1 and TAC1 clade members showed distinct tissue expression preferences. Coexpression analysis revealed that the LAZY1 members had strong associations with cell wall biosynthesis genes. This analysis provides a deeper understanding of the BnIGT gene family and will facilitate further deduction of their role in regulating plant architecture in B. napus.

Genome-Wide Identification and Analysis of the Valine-Glutamine Motif-Containing Gene Family in Brassica napus and Functional Characterization of BnMKS1 in Response to Leptosphaeria maculans

Proteins containing valine-glutamine (VQ) motifs play important roles in plant growth and development as well as in defense responses to both abiotic and biotic stresses. Blackleg disease, which is caused by Leptosphaeria maculans, is the most important disease in canola (Brassica napus) worldwide; however, the identification of Brassica napus VQs and their functions in response to blackleg disease have not yet been reported. In this study, we conducted a genome-wide identification and characterization of the VQ gene family in Brassica napus, including chromosome location, phylogenetic relations, gene structure, motif domain, synteny analysis, and cis-elements categorization of their promoter regions. To understand Brassica napus VQ gene function in response to blackleg disease, we overexpressed BnVQ7 (BnaA01g36880D, also known as the mitogen-activated protein kinase 4 substrate 1 [MKS1] gene) in a blackleg-susceptible canola variety, Westar. Overexpression of BnMKS1 in canola did not improve its resistance to blackleg disease at the seedling stage; however, transgenic canola plants overexpressing BnMKS1 displayed an enhanced resistance to L. maculans infection at the adult plant stage. Expression levels of downstream and defense marker genes in cotyledons increased significantly at the necrotrophic stage of L. maculans infection in the overexpression line of BnMKS1, suggesting that the salicylic acid- and jasmonic acid-mediated signaling pathways were both involved in the defense responses. Together, these results suggest that BnMKS1 might play an important role in defense against L. maculans.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Identification of high-copy number long terminal repeat retrotransposons and their expansion in Phalaenopsis orchids

BACKGROUND

Transposable elements (TEs) are fragments of DNA that can insert into new chromosomal locations. They represent a great proportion of eukaryotic genomes. The identification and characterization of TEs facilitates understanding the transpositional activity of TEs with their effects on the orchid genome structure.

RESULTS

We combined the draft whole-genome sequences of Phalaenopsis equestris with BAC end sequences, Roche 454, and Illumina/Solexa, and identified long terminal repeat (LTR) retrotransposons in these genome sequences by using LTRfinder and classified by using Gepard software. Among the 10 families Gypsy-like retrotransposons, three families Gypsy1, Gypsy2, and Gypsy3, contained the most copies among these predicted elements. In addition, six high-copy retrotransposons were identified according to their reads in the sequenced raw data. The 12-kb Orchid-rt1 contains 18,000 copies representing 220 Mbp of the P. equestris genome. Southern blot and slot blot assays showed that these four retrotransposons Gypsy1, Gypsy2, Gypsy3, and Orchid-rt1 contained high copies in the large-genome-size/large-chromosome species P. violacea and P. bellina. Both Orchid-rt1 and Gypsy1 displayed various ratios of copy number for the LTR sequences versus coding sequences among four Phalaenopsis species, including P. violacea and P. bellina and small-genome-size/small-chromosome P. equestris and P. ahprodite subsp. formosana, which suggests that Orchid-rt1 and Gypsy1 have been through various mutations and homologous recombination events. FISH results showed amplification of Orchid-rt1 in the euchromatin regions among the four Phalaenopsis species. The expression levels of Peq018599 encoding copper transporter 1 is highly upregulated with the insertion of Orchid-rt1, while it is down regulated for Peq009948 and Peq014239 encoding for a 26S proteasome non-ATP regulatory subunit 4 homolog and auxin-responsive factor AUX/IAA-related. In addition, insertion of Orchid-rt1 in these three genes are all in their intron regions.

CONCLUSION

Orchid-rt1 and Gypsy1-3 have amplified within Phalaenopsis orchids concomitant with the expanded genome sizes, and Orchid-rt1 and Gypsy1 may have gone through various mutations and homologous recombination events. Insertion of Orchid-rt1 is in the introns and affects gene expression levels.

Identification of Candidate Auxin Response Factors Involved in Pomegranate Seed Coat Development

Auxin response factors (ARFs) are transcription factors, regulating the auxin signaling pathways involved in plant development and related processes. In this study, we performed the genome-wide identification and characterization of ARFs in pomegranate and compared them with ARFs from three other species. Seventeen PgrARFs were identified and clustered into four groups, according to their phylogenetic relationship with the remaining 59 ARFs. A recent whole-genome duplication event in pomegranate may have contributed to the expansion and diversification of PgrARFs. Genomic truncation and variant splicing mechanisms contributed to the divergence of PgrARFs, a conclusion that was supported by different exon-intron structures of genes and incomplete conserved domains of PgrARFs in a specific phylogenetic group (group III). Interestingly, the absence of motifs from certain PgrARF genes corresponded to their low transcription levels, which contrasted to the highly expressed PgrARFs with intact motifs. Specifically, PgrARF1 and PgrARF2 highly expressed in both inner and outer seed coat, and phylogenetically related to Arabidopsis orthologs which mediates cell divisions in seed coat. We infer these two PgrARFs might involve in seed coat development through cell divisions in response to auxin regulation. These findings provided information on the characteristics and evolutionary relationships of PgrARFs, but also shed lights on their potential roles during seed coat development in pomegranate.

Genome-wide identification and functional analysis of ARF transcription factors in Brassica juncea var. tumida

Auxin signalling is vital for plant growth and development, from embryogenesis to senescence. Recent studies have shown that auxin regulates biological processes by mediating gene expression through a family of functionally original DNA-binding auxin response factors, which exist in a large multi-gene family in plants. However, to date, no information has been available about characteristics of the ARF gene family in Brassica juncea var. tumida. In this study, 65 B. juncea genes that encode ARF proteins were identified in the B. juncea whole-genome, classified into three phylogenetical groups and found to be widely and randomly distributed in the A-and B-genome. Highly conserved proteins were also found within each ortholog based on gene structure and conserved motifs, as well as clustering level. Furthermore, promoter cis-element analysis of BjARFs demonstrated that these genes affect the levels of plant hormones, such as auxin, salicylic, gibberellin acid, MeJA, abscisic acid, and ethylene. Expression analysis showed that differentially expressed BjARF genes were detected during the seedling stage, tumor stem development and the flowering period of B. juncea. Interestingly, we found that BjARF2b_A, BjARF3b_A, BjARF6b_A, and BjARF17a_B were significantly expressed in tumor stem, and an exogenous auxin assay indicated that these genes were sensitive to auxin and IAA signaling. Moreover, eight of the nine BjARF10/16/17 genes and all of the BjARF6/8 genes were involved in post-transcriptional regulation, targeted by Bj-miR160 and Bj-miR167c, respectively. This analysis provides deeper insight of diversification for ARFs and will facilitate further dissection of ARF gene function in B. juncea.

Comprehensive analyses of the annexin (ANN) gene family in Brassica rapa, Brassica oleracea and Brassica napus reveals their roles in stress response

Annexins (ANN) are a multigene, evolutionarily conserved family of calcium-dependent and phospholipid-binding proteins that play important roles in plant development and stress resistance. However, a systematic comprehensive analysis of ANN genes of Brassicaceae species (Brassica rapa, Brassica oleracea, and Brassica napus) has not yet been reported. In this study, we identified 13, 12, and 26 ANN genes in B. rapa, B. oleracea, and B. napus, respectively. About half of these genes were clustered on various chromosomes. Molecular evolutionary analysis showed that the ANN genes were highly conserved in Brassicaceae species. Transcriptome analysis showed that different group ANN members exhibited varied expression patterns in different tissues and under different (abiotic stress and hormones) treatments. Meanwhile, same group members from Arabidopsis thaliana, B. rapa, B. oleracea, and B. napus demonstrated conserved expression patterns in different tissues. The weighted gene coexpression network analysis (WGCNA) showed that BnaANN genes were induced by methyl jasmonate (MeJA) treatment and played important roles in jasmonate (JA) signaling and multiple stress response in B. napus.

Engineering Multiple Abiotic Stress Tolerance in Canola, Brassica napus

Impacts of climate change like global warming, drought, flooding, and other extreme events are posing severe challenges to global crop production. Contribution of Brassica napus towards the oilseed industry makes it an essential component of international trade and agroeconomics. Consequences from increasing occurrences of multiple abiotic stresses on this crop are leading to agroeconomic losses making it vital to endow B. napus crop with an ability to survive and maintain yield when faced with simultaneous exposure to multiple abiotic stresses. For an improved understanding of the stress sensing machinery, there is a need for analyzing regulatory pathways of multiple stress-responsive genes and other regulatory elements such as non-coding RNAs. However, our understanding of these pathways and their interactions in B. napus is far from complete. This review outlines the current knowledge of stress-responsive genes and their role in imparting multiple stress tolerance in B. napus. Analysis of network cross-talk through omics data mining is now making it possible to unravel the underlying complexity required for stress sensing and signaling in plants. Novel biotechnological approaches such as transgene-free genome editing and utilization of nanoparticles as gene delivery tools are also discussed. These can contribute to providing solutions for developing climate change resilient B. napus varieties with reduced regulatory limitations. The potential ability of synthetic biology to engineer and modify networks through fine-tuning of stress regulatory elements for plant responses to stress adaption is also highlighted.

Morphological Characterization of Flower Buds Development and Related Gene Expression Profiling at Bud Break Stage in Heterodichogamous Cyclocarya paliurus (Batal.) lljinskaja

Cyclocarya paliurus (Batal.) Iljinskaja, a unique species growing in southern China, is a multi-function tree species with medicinal, healthcare, material, and ornamental values. So far, sexual reproduction is the main method for extensive cultivation of C. paliurus plantations, but this is limited by low seed plumpness resulted from the character of heterodichogamy. Phenological observations have revealed the asynchronism of flower development in this species. However, its molecular mechanism remains largely unknown. To reveal molecular mechanism of heterodichogamy in C. paliurus, transcriptome of female (F) and male (M) buds from two mating types (protandry, PA; protogyny, PG) at bud break stage were sequenced using Illumina Hiseq 4000 platform. The expression patterns of both 32 genes related to flowering and 58 differentially expressed transcription factors (DETFs) selected from 6 families were divided four groups (PG-F, PG-M, PA-F, and PA-M) into two categories: first flowers (PG-F and PA-M) and later flowers (PA-F and PG-M). The results indicated that genes related to plant hormones (IAA, ABA, and GA) synthesis and response, glucose metabolism, and transcription factors (especially in MIKC family) played significant roles in regulating asynchronism of male and female flowers in the same mating type. The expression of DETFs showed two patterns. One contained DETFs up-regulated in first flowers in comparison to later flowers, and the other was the reverse. Nine genes related to flowering were selected for qRT-PCR to confirm the accuracy of RNA-seq, and generally, the RPKM values of these genes were consistent with the result of qRT-PCR. The results of this work could improve our understanding in asynchronism of floral development within one mating type in C. paliurus at transcriptional level, as well as lay a foundation for further study in heterodichogamous plants.