Functional analysis of the two Brassica AP3 genes involved in apetalous and stamen carpelloid phenotypes

Genome-wide identification and expression pattern analysis of MIKC-Type MADS-box genes in Chionanthus retusus, an androdioecy plant

BACKGROUND

The MADS-box gene family is widely distributed in the plant kingdom, and its members typically encoding transcription factors to regulate various aspects of plant growth and development. In particular, the MIKC-type MADS-box genes play a crucial role in the determination of floral organ development and identity recognition. As a type of androdioecy plant, Chionanthus retusus have unique gender differentiation. Manifested as male individuals with only male flowers and female individuals with only bisexual flowers. However, due to the lack of reference genome information, the characteristics of MIKC-type MADS-box genes in C. retusus and its role in gender differentiation of C. retusus remain largely unknown. Therefore, it is necessary to identify and characterize the MADS-box gene family within the genome of the C. retusus.

RESULTS

In this study, we performed a genome-wide identification and analysis of MIKC-type MADS-box genes in C. retusus (2n = 2x = 46), utilizing the latest reference genome, and studied its expression pattern in individuals of different genders. As a result, we identified a total of 61 MIKC-type MADS-box genes in C. retusus. 61 MIKC-type MADS-box genes can be divided into 12 subfamilies and distributed on 18 chromosomes. Genome collinearity analysis revealed their conservation in evolution, while gene structure, domains and motif analysis indicated their conservation in structure. Finally, based on their expression patterns in floral organs of different sexes, we have identified that CrMADS45 and CrMADS60 may potentially be involved in the gender differentiation of C. retusus.

CONCLUSIONS

Our studies have provided a general understanding of the conservation and characteristics of the MIKC-type MADS-box genes family in C. retusus. And it has been demonstrated that members of the AG subfamily, CrMADS45 and CrMADS60, may play important roles in the gender differentiation of C. retusus. This provides a reference for future breeding efforts to improve flower types in C. retusus and further investigate the role of MIKC-type MADS-box genes in gender differentiation.

Combined Transcriptome and Proteome Analysis Provides Insights into Petaloidy in Pomegranate

Petaloidy leads to a plump floral pattern and increases the landscape value of ornamental pomegranates; however, research on the mechanism of petaloidy in ornamental pomegranates is limited. In this study, we aimed to screen candidate genes related to petaloidy. We performed transcriptomic and proteomic sequencing of the stamens and petals of single-petal and double-petal flowers of ornamental pomegranates. Briefly, 24,567 genes and 5865 proteins were identified, of which 5721 genes were quantified at both transcriptional and translational levels. In the petal and stamen comparison groups, the association between differentially abundant proteins (DAPs) and differentially expressed genes (DEGs) was higher than that between all genes and all proteins, indicating that petaloidy impacts the correlation between genes and proteins. The enrichment results of transcriptome, proteome, and correlation analyses showed that cell wall metabolism, jasmonic acid signal transduction, redox balance, and transmembrane transport affected petaloidy. Nine hormone-related DEGs/DAPs were selected, among which ARF, ILR1, LAX2, and JAR1 may promote petal doubling. Sixteen transcription factor DEGs/DAPs were selected, among which EREBP, LOB, MEF2, MYB, C3H, and trihelix may promote petal doubling. Our results provide transcriptomic and proteomic data on the formation mechanism of petaloidy and a theoretical basis for breeding new ornamental pomegranate varieties.

miR319-Regulated TCP3 Modulates Silique Development Associated with Seed Shattering in Brassicaceae

Seed shattering is an undesirable trait that leads to crop yield loss. Improving silique resistance to shattering is critical for grain and oil crops. In this study, we found that miR319-targeted TEOSINTE BRANCHED 1, CYCLOIDEA, and PROLIFERATING CELL NUCLEAR ANTIGEN BINDING FACTOR (TCPs) inhibited the process of post-fertilized fruits (silique) elongation and dehiscence via regulation of FRUITFULL (FUL) expression in Arabidopsis thaliana and Brassica napus. AtMIR319a activation resulted in a longer silique with thickened and lignified replum, whereas overexpression of an miR319a-resistant version of AtTCP3 (mTCP3) led to a short silique with narrow and less lignified replum. Further genetic and expressional analysis suggested that FUL acted downstream of TCP3 to negatively regulate silique development. Moreover, hyper-activation of BnTCP3.A8, a B. napus homolog of AtTCP3, in rapeseed resulted in an enhanced silique resistance to shattering due to attenuated replum development. Taken together, our findings advance our knowledge of TCP-regulated silique development and provide a potential target for genetic manipulation to reduce silique shattering in Brassica crops.

BnA.JAZ5 Attenuates Drought Tolerance in Rapeseed through Mediation of ABA–JA Crosstalk

Drought stress reduces water availability in plant cells and influences rapeseed yield. Currently, key genetic regulators that contribute to rapeseed response to drought remain largely unexplored, which limits breeding of drought-resistant rapeseed. In this study, we found that Brassica napus JASMONATE ZIM-DOMAIN 5 (BnA.JAZ5), one of the transcriptional repressors functioning in the jasmonate (JA) signaling pathway, was triggered by drought treatment in rapeseed, and drought-susceptibility increased in BnA.JAZ5-overexpressing rapeseed plants as compared to wild-type plants, resulting in a lower survival rate after recovery from dehydration. After recovery for 3 days, 22–40% of p35S::BnA.JAZ5 transgenic plants survived, while approximately 61% of wild-type plants survived. Additionally, seed germination of BnA.JAZ5-overexpressing rapeseed was hyposensitive to abscisic acid (ABA). The germination rate of five transgenic lines was 32~42% under 9 µM ABA treatment, while the germination rate of wild-type plants was 14%. We also found that the average stomatal density of five overexpressing lines was 371~446/mm2, which is higher than that of wild-type (232/mm2) plants under normal conditions. These results indicate that BnA.JAZ5 regulated drought response in an ABA-dependent manner, possibly by affecting stomatal density. Interestingly, methyl jasmonate (MeJA) treatment rescued the ABA-hyposensitive seed germination, revealing crosstalk between JAZ5-meidated JA and the ABA signaling pathway. Taken together, our results suggest that BnA.JAZ5 attenuated drought resistance through the ABA-dependent pathway, which could represent important genetic loci for drought-resistant rapeseed breeding.

Elucidating Mitochondrial DNA Markers of Ogura-Based CMS Lines in Indian Cauliflowers (Brassica oleracea var. botrytis L.) and Their Floral Abnormalities Due to Diversity in Cytonuclear Interactions

Mitochondrial markers can be used to differentiate diverse mitotypes as well as cytoplasms in angiosperms. In cauliflower, cultivation of hybrids is pivotal in remunerative agriculture and cytoplasmic male sterile lines constitute an important component of the hybrid breeding. In diversifying the source of male sterility, it is essential to appropriately differentiate among the available male sterile cytoplasms in cauliflower. PCR polymorphism at the key mitochondrial genes associated with male sterility will be instrumental in analyzing, molecular characterization, and development of mitotype-specific markers for differentiation of different cytoplasmic sources. Presence of auto- and alloplasmic cytonuclear combinations result in complex floral abnormalities. In this context, the present investigation highlighted the utility of organelle genome-based markers in distinguishing cytoplasm types in Indian cauliflowers and unveils the epistatic effects of the cytonuclear interactions influencing floral phenotypes. In PCR-based analysis using a set of primers targeted to orf-138, 76 Indian cauliflower lines depicted the presence of Ogura cytoplasm albeit the amplicons generated exhibited polymorphism within the ofr-138 sequence. The polymorphic fragments were found to be spanning over 200-280 bp and 410-470 bp genomic regions of BnTR4 and orf125, respectively. Sequence analysis revealed that such cytoplasmic genetic variations could be attributed to single nucleotide polymorphisms and insertion or deletions of 31/51 nucleotides. The cytoplasmic effects on varying nuclear-genetic backgrounds rendered an array of floral abnormalities like reduction in flower size, fused flowers, splitted style with the exposed ovule, absence of nonfunctional stamens, and petaloid stamens. These floral malformations caused dysplasia of flower structure affecting female fertility with inefficient nectar production. The finding provides an important reference to ameliorate understanding of mechanism of cytonuclear interactions in floral organ development in Brassicas. The study paves the way for unraveling developmental biology of CMS phenotypes in eukaryotic organisms and intergenomic conflict in plant speciation.

CRISPR/Cas9-Mediated Multiplex Genome Editing of JAGGED Gene in Brassica napus L

Pod shattering resistance is an essential component to achieving a high yield, which is a substantial objective in polyploid rapeseed cultivation. Previous studies have suggested that the Arabidopsis JAGGED (JAG) gene is a key factor implicated in the regulatory web of dehiscence fruit. However, its role in controlling pod shattering resistance in oilseed rape is still unknown. In this study, multiplex genome editing was carried out by the CRISPR/Cas9 system on five homoeologs (BnJAG.A02, BnJAG.C02, BnJAG.C06, BnJAG.A07, and BnJAG.A08) of the JAG gene. Knockout mutagenesis of all homoeologs drastically affected the development of the lateral organs in organizing pod shape and size. The cylindrical body of the pod comprised a number of undifferentiated cells like a callus, without distinctive valves, replum, septum, and valve margins. Pseudoseeds were produced, which were divided into two halves with an incomplete layer of cells (probably septum) that separated the undifferentiated cells. These mutants were not capable of generating any productive seeds for further generations. However, one mutant line was identified in which only a BnJAG.A08-NUB-Like paralog of the JAG gene was mutated. Knockout mutagenesis in BnJAG.A08-NUB gene caused significant changes in the pod dehiscence zone. The replum region of the mutant was increased to a great extent, resulting in enlarged cell size, bumpy fruit, and reduced length compared with the wild type. A higher replum-valve joint area may have increased the resistance to pod shattering by ~2-fold in JAG mutants compared with wild type. Our results offer a basis for understanding variations in Brassica napus fruit by mutating JAG genes and providing a way forward for other Brassicaceae species.

BcMAF2 activates BcTEM1 and represses flowering in Pak-choi (Brassica rapa ssp. chinensis)

BcMAF2 plays a key role in flowering regulation by controlling BcTEM1, BcSOC1 and BCSPL15 in Pak-choi. Flowering is a key event in the life cycle of plants. Flowering time shows an extensive variation from different Pak-choi (Brassica rapa ssp. chinensis) cultivars. However, the regulation mechanism of flowering in Pak-choi remains rarely known. In this study, a systematic identification and functional analysis of a Pak-choi MADS Affecting Flowering (MAF) gene, BcMAF2, was carried out. BcMAF2 encoded a protein containing a conserved MADS-box domain, which was localized in the nucleus. QPCR analysis indicated that the expression of BcMAF2 was higher in the leaves and flowers. Overexpression of BcMAF2 in Arabidopsis showed that BcMAF2 repressed flowering, which was further confirmed by silencing endogenous BcMAF2 in Pak-choi. In addition, Tempranillo 1 (TEM1) expression was up-regulated and MAF2 expression was down-regulated in the BcMAF2-overexpressing Arabidopsis. The expression of BcMAF2 and BcTEM1 was down-regulated in BcMAF2-silencing Pak-choi plants. The yeast one-hybrid, dual luciferase and qPCR results revealed that BcMAF2 protein could directly bind to BcTEM1 promoter and activate its expression, which was not reported in Arabidopsis. Meanwhile, a self-inhibition was found in BcMAF2. Taken together, this work suggested that BcMAF2 could repress flowering by directly activating BcTEM1.

Quantitative Trait Transcripts Mapping Coupled with Expression Quantitative Trait Loci Mapping Reveal the Molecular Network Regulating the Apetalous Characteristic in Brassica napus L

The apetalous trait of rapeseed (Brassica napus, AACC, 2n = 38) is important for breeding an ideal high-yield rapeseed with superior klendusity to Sclerotinia sclerotiorum. Currently, the molecular mechanism underlying the apetalous trait of rapeseed is unclear. In this study, 14 petal regulators genes were chosen as target genes (TGs), and the expression patterns of the 14 TGs in the AH population, containing 189 recombinant inbred lines derived from a cross between apetalous "APL01" and normal "Holly," were analyzed in two environments using qRT-PCR. Phenotypic data of petalous degree (PDgr) in the AH population were obtained from the two environments. Both quantitative trait transcript (QTT)-association mapping and expression QTL (eQTL) analyses of TGs expression levels were performed to reveal regulatory relationships among TGs and PDgr. QTT mapping for PDgr determined that PLURIPETALA (PLP) was the major negative QTT associated with PDgr in both environments, suggesting that PLP negatively regulates the petal development of line "APL01." The QTT mapping of PLP expression levels showed that CHROMATIN-REMODELING PROTEIN 11 (CHR11) was positively associated with PLP expression, indicating that CHR11 acts as a positive regulator of PLP expression. Similarly, QTT mapping for the remaining TGs identified 38 QTTs, associated with 13 TGs, and 31 QTTs, associated with 10 TGs, respectively, in the first and second environments. Additionally, eQTL analyses of TG expression levels showed that 12 and 11 unconditional eQTLs were detected in the first and second environment, respectively. Based on the QTTs and unconditional eQTLs detected, we presented a hypothetical molecular regulatory network in which 14 petal regulators potentially regulated the apetalous trait in "APL01" through the CHR11-PLP pathway. PLP acts directly as the terminal signal integrator negatively regulating petal development in the CHR11-PLP pathway. These findings will aid in the understanding the molecular mechanism underlying the apetalous trait of rapeseed.

Genome-wide transcriptomic analysis uncovers the molecular basis underlying early flowering and apetalous characteristic in Brassica napus L

Floral transition and petal onset, as two main aspects of flower development, are crucial to rapeseed evolutionary success and yield formation. Currently, very little is known regarding the genetic architecture that regulates flowering time and petal morphogenesis in Brassica napus. In the present study, a genome-wide transcriptomic analysis was performed with an absolutely apetalous and early flowering line, APL01, and a normally petalled line, PL01, using high-throughput RNA sequencing. In total, 13,205 differential expressed genes were detected, of which 6111 genes were significantly down-regulated, while 7094 genes were significantly up-regulated in the young inflorescences of APL01 compared with PL01. The expression levels of a vast number of genes involved in protein biosynthesis were altered in response to the early flowering and apetalous character. Based on the putative rapeseed flowering genes, an early flowering network, mainly comprised of vernalization and photoperiod pathways, was built. Additionally, 36 putative upstream genes possibly governing the apetalous character of line APL01 were identified, and six genes potentially regulating petal origination were obtained by combining with three petal-related quantitative trait loci. These findings will facilitate understanding of the molecular mechanisms underlying floral transition and petal initiation in B. napus.

Organelle Simple Sequence Repeat Markers Help to Distinguish Carpelloid Stamen and Normal Cytoplasmic Male Sterile Sources in Broccoli

We previously discovered carpelloid stamens when breeding cytoplasmic male sterile lines in broccoli (Brassica oleracea var. italica). In this study, hybrids and multiple backcrosses were produced from different cytoplasmic male sterile carpelloid stamen sources and maintainer lines. Carpelloid stamens caused dysplasia of the flower structure and led to hooked or coiled siliques with poor seed setting, which were inherited in a maternal fashion. Using four distinct carpelloid stamens and twelve distinct normal stamens from cytoplasmic male sterile sources and one maintainer, we used 21 mitochondrial simple sequence repeat (mtSSR) primers and 32 chloroplast SSR primers to identify a mitochondrial marker, mtSSR2, that can differentiate between the cytoplasm of carpelloid and normal stamens. Thereafter, mtSSR2 was used to identify another 34 broccoli accessions, with an accuracy rate of 100%. Analysis of the polymorphic sequences revealed that the mtSSR2 open reading frame of carpelloid stamen sterile sources had a deletion of 51 bases (encoding 18 amino acids) compared with normal stamen materials. The open reading frame is located in the coding region of orf125 and orf108 of the mitochondrial genomes in Brassica crops and had the highest similarity with Raphanus sativus and Brassica carinata. The current study has not only identified a useful molecular marker to detect the cytoplasm of carpelloid stamens during broccoli breeding, but it also provides evidence that the mitochondrial genome is maternally inherited and provides a basis for studying the effect of the cytoplasm on flower organ development in plants.

Comparative expression analysis of senescence gene CsNAP and B-class floral development gene CsAP3 during different stages of flower development in Saffron (Crocus sativus L.)

Crocus sativus, a monocot triploid species belonging to the Iridaceae family, is cultivated for its red stigmatic lobes of the carpel that constitute saffron. Flower development has been extensively studied in different plants. Different floral developmental pathways have been deciphered in many plants. In Crocus sativus, flower is the most important part and understanding the pathway underlying the flower development can pave the way for new avenues to improve its productivity and quality. The combination of class A genes (including APETALA1; CsAP1 and APETALA2; CsAP2), class B genes (including APETALA3; CsAP3 and PISTILLATA; CsPI) and class C genes (including AGAMOUS; CsAG) that are active in each whorl, determines the identity of the organs that will later develop in that whorl. CsAP3 is a class B homeotic gene which promotes petal and stamen formation and has a very important role in flower development. It also activates other genes playing pivotal role in flower development. It has been earlier reported that CsAP3 gene has direct role in activation of CsNAP gene which promotes senescence in plants. Present work was focused on study of relative gene expression changes of CsAP3 and CsNAP gene during different stages of flower development. CsAP3 gene expression was found maximum during late-preanthesis stages of stigma development. Expression increases from stage 5 to stage 6 of flower development and then reduces again from stage 6 to stage 7. CsNAP gene had moderate expression during stage 3 to stage 4 transition and its expression increased abruptly from stage 6 to stage 7 of flower development. There is no direct concordance in the expression of CsAP3 and CsNAP gene expression in saffron. We may conclude that some other factor(s) may be responsible for initiation of CsNAP expression and CsAP3 gene may directly/indirectly be involved in regulating the factors responsible for CsNAP activation.

High-Density SNP Map Construction and QTL Identification for the Apetalous Character in Brassica napus L

The apetalous genotype is a morphological ideotype for increasing seed yield and should be of considerable agricultural use; however, only a few studies have focused on the genetic control of this trait in Brassica napus. In the present study, a recombinant inbred line, the AH population, containing 189 individuals was derived from a cross between an apetalous line 'APL01' and a normally petalled variety 'Holly'. The Brassica 60 K Infinium BeadChip Array harboring 52,157 single nucleotide polymorphism (SNP) markers was used to genotype the AH individuals. A high-density genetic linkage map was constructed based on 2,755 bins involving 11,458 SNPs and 57 simple sequence repeats, and was used to identify loci associated with petalous degree (PDgr). The linkage map covered 2,027.53 cM, with an average marker interval of 0.72 cM. The AH map had good collinearity with the B. napus reference genome, indicating its high quality and accuracy. After phenotypic analyses across five different experiments, a total of 19 identified quantitative trait loci (QTLs) distributed across chromosomes A3, A5, A6, A9 and C8 were obtained, and these QTLs were further integrated into nine consensus QTLs by a meta-analysis. Interestingly, the major QTL qPD.C8-2 was consistently detected in all five experiments, and qPD.A9-2 and qPD.C8-3 were stably expressed in four experiments. Comparative mapping between the AH map and the B. napus reference genome suggested that there were 328 genes underlying the confidence intervals of the three steady QTLs. Based on the Gene Ontology assignments of 52 genes to the regulation of floral development in published studies, 146 genes were considered as potential candidate genes for PDgr. The current study carried out a QTL analysis for PDgr using a high-density SNP map in B. napus, providing novel targets for improving seed yield. These results advanced our understanding of the genetic control of PDgr regulation in B. napus.

Expression pattern of Class B gene PAP3 in flower development of pepper

Class B gene APETALA3 (AP3) plays a key role in the development of petals and stamens. Here, we investigated the expression pattern of PAP3 gene (genbank accession number: HM104635) in the buds of cytoplasmic male sterility line 121A and its near-isogenic restorer line 121C at four developmental stages and analyzed the possible association between Class B genes and cytoplasmic male sterility of pepper. Semi-quantitative PCR and quantitative real-time RT-PCR (qRT-PCR) as well as RNA in situ hybridization showed increased expression of PAP3 at late phase of anther development and its higher expression in restorer line compared with sterility line indicating PAP3’s role at late developmental stage of anther and suppressed expression in sterility line. RNA in situ hybridization showed Class B gene features: high abundance in stamen and petal; lower expression in pistil; no expression in sepal. Results of transient expression in onion epidermal cells also showed PAP3 localized in the nucleus, which is consistent with the expression pattern of transcription factors of MADS-box gene family.

Identification of miRNAs and their targets from Brassica napus by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are endogenous regulators of a broad range of physiological processes and act by either degrading mRNA or blocking its translation. Oilseed rape (Brassica napus) is one of the most important crops in China, Europe and other Asian countries with publicly available expressed sequence tags (ESTs) and genomic survey sequence (GSS) databases, but little is known about its miRNAs and their targets. To date, only 46 miRNAs have been identified in B. napus.

RESULTS

Forty-one conserved and 62 brassica-specific candidate B. napus miRNAs, including 20 miRNA* sequences, were identified using Solexa sequencing technology. Furthermore, 33 non-redundant mRNA targets of conserved brassica miRNAs and 19 new non-redundant mRNA targets of novel brassica-specific miRNAs were identified by genome-scale sequencing of mRNA degradome.

CONCLUSIONS

This study describes large scale cloning and characterization of B. napus miRNAs and their potential targets, providing the foundation for further characterization of miRNA function in the regulation of diverse physiological processes in B. napus.