Genome wide analysis of MADS-box gene family in Brassica oleracea reveals conservation and variation in flower development

Evolutionary analysis of MADS-box genes in buckwheat species and functional study of FdMADS28 in flavonoid metabolism

The MADS-box gene family is a transcription factor family that is widely expressed in plants. It controls secondary metabolic processes in plants and encourages the development of tissues like roots and flowers. However, the phylogenetic analysis and evolutionary model of MADS-box genes in Fagopyrum species has not been reported yet. This study identified the MADS-box genes of three buckwheat species at the whole genome level, and conducted systematic evolution and physicochemical analysis. The results showed that these genes can be divided into four subfamilies, with fragment duplication being the main way for the gene family expansion. During the domestication process from golden buckwheat to tartary buckwheat and the common buckwheat, the Ka/Ks ratio indicated that most members of the family experienced strong purification selection pressure, and with individual gene pairs experiencing positive selection. In addition, we combined the expression profile data of the MADS genes, mGWAS data, and WGCNA data to mine genes FdMADS28/48/50 that may be related to flavonoid metabolism. The results also showed that overexpression of FdMADS28 could increase rutin content by decreasing Kaempferol pathway content in hairy roots, and increase the resistance and growth of hairy roots to PEG and NaCl. This study systematically analyzed the evolutionary relationship of MADS-box genes in the buckwheat species, and elaborated on the expression patterns of MADS genes in different tissues under biotic and abiotic stresses, laying an important theoretical foundation for further elucidating their role in flavonoid metabolism.

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Hydrangea arborescens 'Annabelle' Flower Formation and Flowering in the Current Year

The perennial woody plant Hydrangea arborescens 'Annabelle' is of great research value due to its unique mechanism of flower development that occurs in the current year, resulting in decorative flowers that can be enjoyed for a relatively long period of time. However, the mechanisms underlying the regulation of current-year flower development in H. arborescens 'Annabelle' are still not fully understood. In this study, we conducted an associated analysis to explore the core regulating network in H. arborescens 'Annabelle' by combining phenological observations, physiological assays, and transcriptome comparisons across seven flower developmental stages. Through this analysis, we constructed a gene co-expression network (GCN) based on the highest reciprocal rank (HRR), using 509 differentially expressed genes (DEGs) identified from seven flowering-related pathways, as well as the biosynthesis of eight flowering-related phytohormones and signal transduction in the transcriptomic analysis. According to the analysis of the GCN, we identified 14 key genes with the highest functional connectivity that played critical roles in specific development stages. We confirmed that 135 transcription factors (AP2/ERF, bHLH, CO-like, GRAS, MIKC, SBP, WRKY) were highly co-expressed with the 14 key genes, indicating their close associations with the development of current-year flowers. We further proposed a hypothetical model of a gene regulatory network for the development of the whole flower. This model suggested that the photoperiod, aging, and gibberellin pathways, along with the phytohormones abscisic acid (ABA), gibberellin (GA), brassinosteroid (BR), and jasmonic acid (JA), work synergistically to promote the floral transition. Additionally, auxin, GA, JA, ABA, and salicylic acid (SA) regulated the blooming process by involving the circadian clock. Cytokinin (CTK), ethylene (ETH), and SA were key regulators that affected flower senescence. Additionally, several floral integrators (HaLFY, HaSOC1-2, HaAP1, HaFULL, HaAGL24, HaFLC, etc.) were dominant contributors to the development of H. arborescens flowers. Overall, this research provides a comprehensive understanding of the dynamic mechanism underlying the entire process of current-year flower development, thereby offering valuable insights for further studies on the flower development of H. arborescens 'Annabelle'.

Genome-Wide Identification of AhMDHs and Analysis of Gene Expression under Manganese Toxicity Stress in Arachis hypogaea

Malate dehydrogenase (MDH) is one kind of oxidation-reduction enzyme that catalyzes the reversible conversion of oxaloacetic acid to malic acid. It has vital functions in plant development, photosynthesis, abiotic stress responses, and so on. However, there are no reports on the genome-wide identification and gene expression of the MDH gene family in Arachis hypogaea. In this study, the MDH gene family of A. hypogaea was comprehensively analyzed for the first time, and 15 AhMDH sequences were identified. According to the phylogenetic tree analysis, AhMDHs are mainly separated into three subfamilies with similar gene structures. Based on previously reported transcriptome sequencing results, the AhMDH expression quantity of roots and leaves exposed to manganese (Mn) toxicity were explored in A. hypogaea. Results revealed that many AhMDHs were upregulated when exposed to Mn toxicity, suggesting that those AhMDHs might play an important regulatory role in A. hypogaea's response to Mn toxicity stress. This study lays foundations for the functional study of AhMDHs and further reveals the mechanism of the A. hypogaea signaling pathway responding to high Mn stress.

Genetic dissection of morphological variation between cauliflower and a rapid cycling Brassica oleracea line

To improve resolution to small genomic regions and sensitivity to small-effect loci in the identification of genetic factors conferring the enlarged inflorescence and other traits of cauliflower while also expediting further genetic dissection, 104 near-isogenic introgression lines (NIILs) covering 78.56% of the cauliflower genome, were selected from an advanced backcross population using cauliflower [Brassica oleracea var. botrytis L., mutant for Orange gene (ORG)] as the donor parent and a rapid cycling line (TO1434) as recurrent parent. Subsets of the advanced backcross population and NIILs were planted in the field for 8 seasons, finding 141 marker-trait associations for 15 leaf-, stem-, and flower-traits. Exemplifying the usefulness of these lines, we delineated the previously known flower color gene to a 4.5 MB interval on C3; a gene for small plant size to a 3.4 MB region on C8; and a gene for large plant size and flowering time to a 6.1 MB region on C9. This approach unmasked closely linked QTL alleles with opposing effects (on chr. 8) and revealed both alleles with expected phenotypic effects and effects opposite the parental phenotypes. Selected B. oleracea NIILs with short generation time add new value to widely used research and teaching materials.

Genome-Wide Identification and Expression Analysis of the MADS Gene Family in Tulips (Tulipa gesneriana)

To investigate the cold response mechanism and low temperature regulation of flowering in tulips, this study identified 32 MADS-box transcription factor family members in tulips based on full-length transcriptome sequencing, named TgMADS1-TgMADS32. Phylogenetic analysis revealed that these genes can be divided into two classes: type I and type II. Structural analysis showed that TgMADS genes from different subfamilies have a similar distribution of conserved motifs. Quantitative real-time PCR results demonstrated that some TgMADS genes (e.g., TgMADS3, TgMADS15, TgMADS16, and TgMADS19) were significantly upregulated in buds and stems under cold conditions, implying their potential involvement in the cold response of tulips. In summary, this study systematically identified MADS family members in tulips and elucidated their evolutionary relationships, gene structures, and cold-responsive expression patterns, laying the foundation for further elucidating the roles of these transcription factors in flowering and the cold adaptability of tulips.

Characterization and Potential Function Analysis of the SRS Gene Family in Brassica napus

SRS (SHI-related sequence) transcription factors play a crucial role in plant growth, development, and abiotic stress response. Although Brassica napus (B. napus) is one of the most important oil crops in the world, the role of SRS genes in B. napus (BnSRS) has not been well investigated. Therefore, we employed a bioinformatics approach to identify BnSRS genes from genomic data and investigated their characteristics, functions, and expression patterns, to gain a better understanding of how this gene family is involved in plant development and growth. The results revealed that there were 34 BnSRS gene family members in the genomic sequence of B. napus, unevenly distributed throughout the sequence. Based on the phylogenetic analysis, these BnSRS genes could be divided into four subgroups, with each group sharing comparable conserved motifs and gene structure. Analysis of the upstream promoter region showed that BnSRS genes may regulate hormone responses, biotic and abiotic stress response, growth, and development in B. napus. The protein-protein interaction analysis revealed the involvement of BnSRS genes in various biological processes and metabolic pathways. Our analysis of BnSRS gene expression showed that 23 BnSRS genes in the callus tissue exhibited a dominant expression pattern, suggesting their critical involvement in cell dedifferentiation, cell division, and tissue development. In addition, association analysis between genotype and agronomic traits revealed that BnSRS genes may be linked to some important agronomic traits in B. napus, suggesting that BnSRS genes were widely involved in the regulation of important agronomic traits (including C16.0, C18.0, C18.1, C18.2 C18.3, C20.1, C22.1, GLU, protein, TSW, and FFT). In this study, we predicted the evolutionary relationships and potential functions of BnSRS gene family members, providing a basis for the development of BnSRS gene functions which could facilitate targeted functional studies and genetic improvement for elite breeding in B. napus.

Molecular evolution, diversification, and expression assessment of MADS gene family in Setaria italica, Setaria viridis, and Panicum virgatum

KEY MESSAGE

This paper sheds light on the evolution and expression patterns of MADS genes in Setaria and Panicum virgatum. SiMADS51 and SiMADS64 maybe involved in the ABA-dependent pathway of drought response. The MADS gene family is a key regulatory factor family that controls growth, reproduction, and response to abiotic stress in plants. However, the molecular evolution of this family is rarely reported. Here, a total of 265 MADS genes were identified in Setaria italica (foxtail millet), Setaria viridis (green millet), and Panicum virgatum (switchgrass) and analyzed by bioinformatics, including physicochemical characteristics, subcellular localization, chromosomal position and duplicate, motif distribution, genetic structure, genetic evolvement, and expression patterns. Phylogenetic analysis was used to categorize these genes into M and MIKC types. The distribution of motifs and gene structure were similar for the corresponding types. According to a collinearity study, the MADS genes have been mostly conserved during evolution. The principal cause of their expansion is segmental duplication. However, the MADS gene family tends to shrink in foxtail millet, green millet, and switchgrass. The MADS genes were subjected to purifying selection, but several positive selection sites were also identified in three species. And most of the promoters of MADS genes contain cis-elements related to stress and hormonal response. RNA-seq and quantitative Real-time PCR (qRT-PCR) analysis also were examined. SiMADS genes expression levels are considerably changed in reaction to various treatments, following qRT-PCR analysis. This sheds fresh light on the evolution and expansion of the MADS family in foxtail millet, green millet, and switchgrass, and lays the foundation for further research on their functions.

Molecular Markers for Detecting Inflorescence Size of Brassica oleracea L. Crops and B. oleracea Complex Species (n = 9) Useful for Breeding of Broccoli (B. oleracea var. italica) and Cauliflower (B. oleracea var. botrytis)

The gene flow from Brassica oleracea L. wild relatives to B. oleracea vegetable crops have occurred and continue to occur ordinarily in several Mediterranean countries, such as Sicily, representing an important hot spot of diversity for some of them, such as broccoli, cauliflower and kale. For detecting and for exploiting the forgotten alleles lost during the domestication processes of the B. oleracea crops, attention has been pointed to the individuation of specific markers for individuating genotypes characterized by hypertrophic inflorescence traits by the marker assisted selection (MAS) during the first plant growing phases after the crosses between broccoli (B. oleracea var. italica)/cauliflower (B. oleracea var. botrytis) with B. oleracea wild relatives (n = 9), reducing the cultivation and evaluation costs. The desired traits often found in several B. oleracea wild relatives are mainly addressed to improve the plant resistance to biotic and abiotic stresses and to increase the organoleptic, nutritive and nutraceutical traits of the products. One of the targeted traits for broccoli and cauliflower breeding is represented by the inflorescences size as is documented by the domestication processes of these two crops. Based on the previous results achieved, the numerical matrix, obtained utilizing five simple sequence repeats (SSRs), was analyzed to assess the relationship among the main inflorescence characteristics and the allelic variation of the SSRs loci analyzed (BoABI1, BoAP1, BoPLD1, BoTHL1 and PBCGSSRBo39), both for the Brassica oleracea and B. oleracea wild relatives (n = 9) accessions set. The main inflorescence morphometric characteristics, such as weight, height, diameter, shape, inflorescence curvature angle and its stem diameter, were registered before the flower anthesis. We analyzed the correlations among the allelic variation of the SSRs primers utilized and the inflorescence morphometric characteristics to individuate genomic regions stimulating the hypertrophy of the reproductive organ. The relationships found explain the diversity among B. oleracea crops and the B. oleracea complex species (n = 9) for the inflorescence size and structure. The individuated markers allow important time reduction during the breeding programs after crossing wild species for transferring useful biotic and abiotic resistances and organoleptic and nutraceutical traits to the B. oleracea crops by MAS.

Systematic analysis of MADS-box gene family in the U's triangle species and targeted mutagenesis of BnaAG homologs to explore its role in floral organ identity in Brassica napus

MADS-box transcription factors play an important role in regulating floral organ development and participate in environmental responses. To date, the MADS-box gene family has been widely identified in Brassica rapa (B. rapa), Brassica oleracea (B. oleracea), and Brassica napus (B. napus); however, there are no analogous reports in Brassica nigra (B. nigra), Brassica juncea (B. juncea), and Brassica carinata (B. carinata). In this study, a whole-genome survey of the MADS-box gene family was performed for the first time in the triangle of U species, and a total of 1430 MADS-box genes were identified. Based on the phylogenetic relationship and classification of MADS-box genes in Arabidopsis thaliana (A. thaliana), 1430 MADS-box genes were categorized as M-type subfamily (627 genes), further divided into Mα, Mβ, Mγ, and Mδ subclades, and MIKC-type subfamily (803 genes), further classified into 35 subclades. Gene structure and conserved protein motifs of MIKC-type MADS-box exhibit diversity and specificity among different subclades. Comparative analysis of gene duplication events and syngenic gene pairs among different species indicated that polyploidy is beneficial for MIKC-type gene expansion. Analysis of transcriptome data within diverse tissues and stresses in B. napus showed tissue-specific expression of MIKC-type genes and a broad response to various abiotic stresses, particularly dehydration stress. In addition, four representative floral organ mutants (wtl, feml, aglf-2, and aglf-1) in the T0 generation were generated by editing four AGAMOUS (BnaAG) homoeologs in B. napus that enriched the floral organ variant phenotype. In brief, this study provides useful information for investigating the function of MADS-box genes and contributes to revealing the regulatory mechanisms of floral organ development in the genetic improvement of new varieties.

Genome-wide identification and expression analysis of the MADS-box gene family during female and male flower development in Juglans mandshurica

The MADS-box gene family plays a crucial role in multiple developmental processes of plants, especially in floral organ specification and the regulation of fruit development and ripening. Juglans mandshurica is a precious fruit material whose quality and yield are determined by floral organ development. The molecular mechanism of J. mandshurica female and male flower development depending on MADS-box genes remains unclear. In our study, 67 JmMADS genes were identified and unevenly distributed on 15 of 16 J. mandshurica chromosomes. These genes were divided into two types [type I (Mα, Mγ, Mδ) and type II (MIKC)]. The gene structure and motif analyses showed that most genes belonging to the same type had similar gene structures and conserved motifs. The analysis of syntenic relationships showed that MADS-box genes in J. mandshurica, J. sigillata, and J. regia exhibited the highest homology and great collinearity. Analysis of cis-acting elements showed that JmMADS gene promoter regions contained light, stress and hormone response cis-acting elements. The gene expression patterns demonstrated that 30 and 26 JmMADS genes were specifically expressed in the female and male flowers, respectively. In addition, 12 selected genes common to J. mandshurica female and male flowers were significantly upregulated at the mature stage and were used to validate the reliability of the transcriptome data using quantitative real-time PCR. This comprehensive and systematic analysis of J. mandshurica MADS-box genes lays a foundation for future studies on MADS-box gene family functions.

Genome-wide analysis of MADS-box families and their expressions in flower organs development of pineapple (Ananas comosus (L.) Merr.)

MADS-box genes play crucial roles in plant vegetative and reproductive growth, better development of inflorescences, flower, and fruit. Pineapple is a typical collective fruit, and a comprehensive analysis of the MADS-box gene family in the development of floral organs of pineapple is still lacking. In this study, the whole-genome survey and expression profiling of the MADS-box family in pineapple were introduced. Forty-four AcMADS genes were identified in pineapple, 39 of them were located on 18 chromosomes and five genes were distributed in five scaffolds. Twenty-two AcMADS genes were defined as 15 pairs of segmental duplication events. Most members of the type II subfamily of AcMADS genes had higher expression levels in floral organs compared with type I subfamily, thereby suggesting that AcMADS of type II may play more crucial roles in the development of floral organs of pineapple. Six AcMADS genes have significant tissue-specificity expression, thereby suggesting that they may participate in the formation of one or more floral organs. This study provides valuable insights into the role of MADS-box gene family in the floral organ development of pineapple.

Genome-Wide Identification and Expression Analysis of the MADS-Box Gene Family in Sweet Potato [Ipomoea batatas (L.) Lam]

MADS-box gene, one of the largest transcription factor families in plants, is a class of transcription factors widely present in eukaryotes. It plays an important role in plant growth and development and participates in the growth and development of flowers and fruits. Sweet potato is the seventh most important food crop in the world. Its tuberous roots, stems, and leaves contain a large number of proteins, lipids, carotenoids, anthocyanins, conjugated phenolic acids, and minerals, which have high edible, forage, and medicinal value, and is also an important energy crop. At present, MADS-box genes in sweet potato have rarely been reported, and there has been no study on the genome-wide identification and classification of MADS-box genes in Ipomoea batatas. This study provided the first comprehensive analysis of sweet potato MADS-box genes. We identified 95 MADS-box genes, analyzed the structure and protein of sweet potato MADS-box genes, and categorized them based on phylogenetic analysis with Arabidopsis MADS-box proteins. Chromosomal localization indicated an unequal number of MADS-box genes in all 14 chromosomes except LG3, with more than 10 MADS-box genes located on chromosomes LG7, LG11, and LG15. The MADS domain and core motifs of the sweet potato MADS-box genes were identified by motif analysis. We identified 19 MADS-box genes with collinear relationships and analyzed duplication events. Cis-acting elements, such as light-responsive, auxin-responsive, drought-inducible, and MeJA-responsive elements, were found in the promoter region of the MADS-box genes in sweet potato, which further indicates the basis of MADS-box gene regulation in response to environmental changes and hormones. RNA-seq suggested that sweet potato MADS-box genes exhibit tissue-specific expression patterns, with 34 genes highly expressed in sweet potato flowers and fruits, and 19 genes highly expressed in the tuberous root, pencil root, or fibrous root. qRT-PCR again validated the expression levels of the 10 genes and found that IbMADS1, IbMADS18, IbMADS19, IbMADS79, and IbMADS90 were highly expressed in the tuberous root or fibrous root, and IbMADS18, IbMADS31, and IbMADS83 were highly expressed in the fruit. In this study, the molecular basis of MADS-box genes of sweet potato was analyzed from various angles. The effects of MADS-box genes on the growth and development of sweet potato were investigated, which may provide a certain theoretical basis for molecular breeding of sweet potato.

The Melastoma dodecandrum genome and the evolution of Myrtales

Melastomataceae have abundant morphological diversity with high economic and ornamental merit in Myrtales. The phylogenetic position of Myrtales is still contested. Here, we report the first chromosome-level genome assembly of Melastoma dodecandrum in Melastomataceae. The assembled genome size was 299.81 Mb with a contig N50 value of 3.00 Mb. Genome evolution analysis indicated that M. dodecandrum, Eucalyptus grandis and Punica granatum were clustered into a clade of Myrtales and formed a sister group with the ancestor of fabids and malvids. We found that M. dodecandrum experienced four whole-genome polyploidization events: the ancient event was shared with most eudicots, one event was shared with Myrtales, and the other two events were unique to M. dodecandrum. Moreover, we identified MADS-box genes and found that the AP1-like genes expanded, and AP3-like genes might have undergone subfunctionalization. We found that the SUAR63-like genes and AG-like genes showed different expression patterns in stamens, which may be associated with heteranthery. In addition, we found that LAZY1-like genes were involved in the negative regulation of stem branching development, which may be related to its creeping features. Our study sheds new light on the evolution of Melastomataceae and Myrtales, which provides a comprehensive genetic resource for future research.

Integrated analysis of mRNA-seq and miRNA-seq reveals the advantage of polyploid Solidago canadensis in sexual reproduction

BACKGROUND

The invasion of Solidago canadensis probably related to polyploidy, which may promotes its potential of sexual reproductive. S. canadensis as an invasive species which rapidly widespread through yield huge numbers of seed, but the mechanism remains unknown. To better understand the advantages of sexual reproduction in hexaploid S. canadensis, transcriptome and small RNA sequencing of diploid and hexaploid cytotypes in flower bud and fruit development stages were performed in this study.

RESULTS

The transcriptome analysis showed that in the flower bud stage, 29 DEGs were MADS-box related genes with 14 up-regulated and 15 down-regulated in hexaploid S. canadensis; 12 SPL genes were detected differentially expressed with 5 up-regulated and 7 down-regulated. In the fruit development stage, 26 MADS-box related genes with 20 up-regulated and 6 down-regulated in hexaploid S. canadensis; 5 SPL genes were all up-regulated; 28 seed storage protein related genes with 18 were up-regulated and 10 down-regulated. The weighted gene co-expression network analysis (WGCNA) identified 19 modules which consisted of co-expressed DEGs with functions such as sexual reproduction, secondary metabolism and transcription factors. Furthermore, we discovered 326 miRNAs with 67 known miRNAs and 259 novel miRNAs. Some of miRNAs, such as miR156, miR156a and miR156f, which target the sexual reproduction related genes.

CONCLUSION

Our study provides a global view of the advantages of sexual reproduction in hexaploid S. canadensis based on the molecular mechanisms, which may promote hexaploid S. canadensis owing higher yield and fruit quality in the process of sexual reproduction and higher germination rate of seeds, and finally conductive to diffusion, faster propagation process and enhanced invasiveness.

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.

Genome-wide identification of MIKC-type genes related to stamen and gynoecium development in Liriodendron

The organogenesis and development of reproductive organs, i.e., stamen and gynoecium, are important floral characteristics that are closely related to pollinators and reproductive fitness. As a genus from Magnoliaceae, Liriodendron has only two relict species: L. chinense and L. tulipifera. Despite the similar flower shapes of these species, their natural seed-setting rates differ significantly, implying interspecies difference in floral organogenesis and development. MADS-box genes, which participate in floral organogenesis and development, remain unexplored in Liriodendron. Here, to explore the interspecies difference in floral organogenesis and development and identify MADS-box genes in Liriodendron, we examined the stamen and gynoecium primordia of the two Liriodendron species by scanning electron microscopy combined with paraffin sectioning, and then collected two types of primordia for RNA-seq. A total of 12 libraries were constructed and 42,268 genes were identified, including 35,269 reference genes and 6,999 new genes. Monoterpenoid biosynthesis was enriched in L. tulipifera. Genome-wide analysis of 32 MADS-box genes was conducted, including phylogenetic trees, exon/intron structures, and conserved motif distributions. Twenty-six genes were anchored on 17 scaffolds, and six new genes had no location information. The expression profiles of MIKC-type genes via RT-qPCR acrossing six stamen and gynoecium developmental stages indicates that the PI-like, AG/STK-like, SEP-like, and SVP-like genes may contribute to the species-specific differentiation of the organogenesis and development of reproductive organs in Liriodendron. Our findings laid the groundwork for the future exploration of the mechanism underlying on the interspecific differences in reproductive organ development and fitness in Liriodendron.

Genome-wide identification and expression analysis of the MADS-box transcription factor family in Camellia sinensis

The MADS-box genes are an important class of transcription factors and play critical roles in flower development. However, the functions of these genes in the economically important drinking plant, Camellia sinensis, are still not reported. Here, an evolutionary analysis of tea MADS-box genes was performed at whole genome level. A total of 83 MADS-box genes were identified in tea, and their gene structures and expression patterns were further analyzed. The tea MADS-box genes were classified into Mα (26), Mβ (12), Mγ (9), MIKC* (7), and MIKCC (29) clade according to their phylogenetic relationship with Arabidopsis thaliana. Several cis-elements were identified in the promoter regions of the CsMADS genes that are important in regulating growth, development, light responses, and the response to several stresses. Most CsMADS genes display clear different expression patterns in different organs and different species of tea plant. The expression of CsMADS genes can be regulated by abiotic stresses and phytohormone treatment. Our results lay the foundation for future research on the function of CsMADS genes and beneficial for improving tea agricultural traits in the future.

Genome-wide identification and comparative analysis of diacylglycerol kinase (DGK) gene family and their expression profiling in Brassica napus under abiotic stress

BACKGROUND

Diacylglycerol kinases (DGKs) are signaling enzymes that play pivotal roles in response to abiotic and biotic stresses by phosphorylating diacylglycerol (DAG) to form phosphatidic acid (PA). However, no comprehensive analysis of the DGK gene family had previously been reported in B. napus and its diploid progenitors (B. rapa and B. oleracea).

RESULTS

In present study, we identified 21, 10, and 11 DGK genes from B. napus, B. rapa, and B. oleracea, respectively, which all contained conserved catalytic domain and were further divided into three clusters. Molecular evolutionary analysis showed that speciation and whole-genome triplication (WGT) was critical for the divergence of duplicated DGK genes. RNA-seq transcriptome data revealed that, with the exception of BnaDGK4 and BnaDGK6, BnaDGK genes have divergent expression patterns in most tissues. Furthermore, some DGK genes were upregulated or downregulated in response to hormone treatment and metal ion (arsenic and cadmium) stress. Quantitative real-time PCR analysis revealed that different BnaDGK genes contribute to seed oil content.

CONCLUSIONS

Together, our results indicate that DGK genes have diverse roles in plant growth and development, hormone response, and metal ion stress, and in determining seed oil content, and lay a foundation for further elucidating the roles of DGKs in Brassica species.

Chromosome-scale assembly of the Kandelia obovata genome

The mangrove Kandelia obovata (Rhizophoraceae) is an important coastal shelterbelt and landscape tree distributed in tropical and subtropical areas across East Asia and Southeast Asia. Herein, a chromosome-level reference genome of K. obovata based on PacBio, Illumina, and Hi-C data is reported. The high-quality assembled genome size is 177.99 Mb, with a contig N50 value of 5.74 Mb. A large number of contracted gene families and a small number of expanded gene families, as well as a small number of repeated sequences, may account for the small K. obovata genome. We found that K. obovata experienced two whole-genome polyploidization events: one whole-genome duplication shared with other Rhizophoreae and one shared with most eudicots (γ event). We confidently annotated 19,138 protein-coding genes in K. obovata and identified the MADS-box gene class and the RPW8 gene class, which might be related to flowering and resistance to powdery mildew in K. obovata and Rhizophora apiculata, respectively. The reference K. obovata genome described here will be very useful for further molecular elucidation of various traits, the breeding of this coastal shelterbelt species, and evolutionary studies with related taxa.

Identification of Candidate Genes Involved in Curd Riceyness in Cauliflower

“Riceyness” refers to the precocious development of flower bud initials over the curd surface of cauliflower, and it is regarded as undesirable for the market. The present study aimed to identify the candidate loci and genes responsible for the morphological difference in riceyness between a pair of cauliflower sister lines. Genetic analysis revealed that riceyness is controlled by a single dominant locus. An F₂ population derived from the cross between these sister lines was used to construct “riceyness” and “non-riceyness” bulks, and then it was subjected to BSA-seq. On the basis of the results of Δ(SNP-index) analysis, a 4.0 Mb candidate region including 22 putative SNPs was mapped on chromosome C04. Combining the RNA-seq, gene function annotation, and target sequence analysis among two parents and other breeding lines, an orthologous gene of the Arabidopsis gene SOC1, Bo4g024850 was presumed as the candidate gene, and an upstream SNP likely resulted in riceyness phenotype via influencing the expression levels of Bo4g024850. These results are helpful to understand the genetic mechanism regulating riceyness, and to facilitate the molecular improvement on cauliflower curds.