The Transcriptome of Cunninghamia lanceolata male/female cone reveal the association between MIKC MADS-box genes and reproductive organs development

Molecular regulatory mechanisms of staminate strobilus development and dehiscence in Torreya grandis

Gymnosperms are mostly dioecious, and their staminate strobili undergo a longer developmental period than those of angiosperms. However, the underlying molecular mechanisms remain unclear. This study aimed to identify key genes and pathways involved in staminate strobilus development and dehiscence in Torreya grandis. Through weighted gene co-expression network analysis (WGCNA), we identified fast elongation-related genes enriched in carbon metabolism and auxin signal transduction, whereas dehiscence-related genes were abundant in alpha-linolenic acid metabolism and the phenylpropanoid pathway. Based on WGCNA, we also identified PHYTOCHROME-INTERACTING FACTOR4 (TgPIF4) as a potential regulator for fast elongation of staminate strobilus and 2 WRKY proteins (TgWRKY3 and TgWRKY31) as potential regulators for staminate strobilus dehiscence. Multiple protein-DNA interaction analyses showed that TgPIF4 directly activates the expression of TRANSPORT INHIBITOR RESPONSE2 (TgTIR2) and NADP-MALIC ENZYME (TgNADP-ME). Overexpression of TgPIF4 significantly promoted staminate strobilus elongation by elevating auxin signal transduction and pyruvate content. TgWRKY3 and TgWRKY31 bind to the promoters of the lignin biosynthesis gene PHENYLALANINE AMMONIA-LYASE (TgPAL) and jasmonic acid metabolism gene JASMONATE O-METHYLTRANSFERASE (TgJMT), respectively, and directly activate their transcription. Overexpression of TgWRKY3 and TgWRKY31 in the staminate strobilus led to early dehiscence, accompanied by increased lignin and methyl jasmonate levels, respectively. Collectively, our findings offer a perspective for understanding the growth of staminate strobili in gymnosperms.

Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

Comparative transcriptome analysis reveals candidate genes related to the sex differentiation of Schisandra chinensis

Schisandra chinensis is a monoecious plant with unisex flowers. The fruit of S. chinensis is of high medical with economic value. The yield of S. chinensis fruit is related to the ratio of its female and male flowers. However, there is little research on its floral development and sex differentiation. To elucidate the possible mechanism for the sex differentiation of S. chinensis, we collected 18 samples of female and male flowers from three developmental stages and performed a comparative RNA-seq analysis aimed at identifying differentially expressed genes (DEGs) that may be related to sex differentiation. The results showed 936, 7179, and 6890 differentially expressed genes between female and male flowers at three developmental stages, respectively, and 466 candidate genes may play roles in sex differentiation. KEGG analysis showed genes involved in the flavonoid biosynthesis pathway and DNA replication pathway were essential for the development of female flowers. 51 MADS-box genes and 10 YABBY genes were identified in S. chinensis. The DEGs analysis indicated that MADS-box and YABBY genes were strongly related to the sex determination of S. chinensis. RT-qPCR confirmed the RNA-seq results of 20 differentially expressed genes, including three male-biased genes and 17 female-biased genes. A possible regulatory model of sex differentiation in S. chinensis was proposed according to our results. This study helps reveal the sex-differentiation mechanism of S. chinensis and lays the foundation for regulating the male-female ratio of S. chinensis in the future.

Identification and Expression of the MADS-box Gene Family in Different Versions of the Ginkgo biloba Genome

MADS-box transcription factors play important roles in many organisms. These transcription factors are involved in processes such as the formation of the flower organ structure and the seed development of plants. Ginkgo biloba has two genome versions (version 2019 and version 2021), and there is no analysis or comparison of the MADS-box gene family in these two genomes. In this study, 26 and 20 MADS-box genes were identified from the two genomes of Ginkgo, of which 12 pairs of genes reached more than 80% similarity. According to our phylogenetic analysis results, we divided these genes into type I (Mα and Mγ subfamilies) and type II (MIKC and Mδ subfamilies) members. We found that both sets of genomes lacked the Mβ gene, while the MIKC gene was the most numerous. Further analysis of the gene structure showed that the MIKC genes in the two genomes had extralong introns (≥20 kb); these introns had different splicing patterns, and their expression might be more abundant. The gene expression analysis proved that GbMADS genes were expressed to varying degrees in eight Ginkgo biological tissues. Type II GbMADS genes not only were found to be related to female flower bud differentiation and development but also are important in seed development. Therefore, MADS-box genes may play important roles in the development of Ginkgo reproductive organs, which may suggest a genetic role in sexual differentiation. This study further contributes to the research on MADS-box genes and provides new insights into sex determination in Ginkgo.

Genomic and Transcriptomic Insights into the Evolution and Divergence of MIKC-Type MADS-Box Genes in Carica papaya

MIKC-type MADS-box genes, also known as type II genes, play a crucial role in regulating the formation of floral organs and reproductive development in plants. However, the genome-wide identification and characterization of type II genes as well as a transcriptomic survey of their potential roles in Carica papaya remain unresolved. Here, we identified and characterized 24 type II genes in the C. papaya genome, and investigated their evolutional scenario and potential roles with a widespread expression profile. The type II genes were divided into thirteen subclades, and gene loss events likely occurred in papaya, as evidenced by the contracted member size of most subclades. Gene duplication mainly contributed to MIKC-type gene formation in papaya, and the duplicated gene pairs displayed prevalent expression divergence, implying the evolutionary significance of gene duplication in shaping the diversity of type II genes in papaya. A large-scale transcriptome analysis of 152 samples indicated that different subclasses of these genes showed distinct expression patterns in various tissues, biotic stress response, and abiotic stress response, reflecting their divergent functions. The hub-network of male and female flowers and qRT-PCR suggested that TT16-3 and AGL8 participated in male flower development and seed germination. Overall, this study provides valuable insights into the evolution and functions of MIKC-type genes in C. papaya.

Genome-Wide Identification of the MADS-Box Gene Family during Male and Female Flower Development in Chayote (Sechium edule)

The MADS-box gene plays an important role in plant growth and development. As an important vegetable of Cucurbitaceae, chayote has great edible and medicinal value. So far, there is little molecular research on chayote, and there are no reports on the MADS-box transcription factor of chayote. In this study, the MADS-box gene family of chayote was analyzed for the first time, and a total of 70 MADS-box genes were identified, including 14 type I and 56 type II MICK MADS genes. They were randomly distributed on 13 chromosomes except for chromosome 11. The light response element, hormone response element and abiotic stress response element were found in the promoter region of 70 MADS genes, indicating that the MADS gene can regulate the growth and development of chayote, resist abiotic stress, and participate in hormone response; GO and KEGG enrichment analysis also found that SeMADS genes were mainly enriched in biological regulation and signal regulation, which further proved the important role of MADS-box gene in plant growth and development. The results of collinearity showed that segmental duplication was the main driving force of MADS gene expansion in chayote. RNA-seq showed that the expression levels of SeMADS06, SeMADS13, SeMADS26, SeMADS28, SeMADS36 and SeMADS37 gradually increased with the growth of chayote, indicating that these genes may be related to the development of root tubers of 'Tuershao'. The gene expression patterns showed that 12 SeMADS genes were specifically expressed in the male flower in 'Tuershao' and chayote. In addition, SeMADS03 and SeMADS52 may be involved in regulating the maturation of male flowers of 'Tuershao' and chayote. SeMADS21 may be the crucial gene in the development stage of the female flower of 'Tuershao'. This study laid a theoretical foundation for the further study of the function of the MADS gene in chayote in the future.

Identifying Genes Associated with Female Flower Development of Phellodendron amurense Rupr. Using a Transcriptomics Approach

Phellodendron amurense Rupr., a species of Rutaceae, is a nationally protected and valuable medicinal plant. It is generally considered to be dioecious. With the discovery of monoecious P. amurense, the phenomenon that its sex development is regulated by epigenetics has been revealed, but the way epigenetics affects the sex differentiation of P. amurense is still unclear. In this study, we investigated the effect of DNA methylation on the sexual development of P. amurense. The young inflorescences of male plants were treated with the demethylation agent 5-azaC, and the induced female flowers were obtained. The induced female flowers’ morphological functions and transcriptome levels were close to those of normally developed plants. Genes associated with the development of female flowers were studied by comparing the differences in transcriptome levels between the male and female flowers. Referring to sex-related genes reported in other plants, 188 candidate genes related to the development of female flowers were obtained, including sex-regulating genes, genes related to the formation and development of sexual organs, genes related to biochemical pathways, and hormone-related genes. RPP0W, PAL3, MCM2, MCM6, SUP, PIN1, AINTEGUMENTA, AINTEGUMENTA-LIKE6, AGL11, SEUSS, SHI-RELATED SEQUENCE 5, and ESR2 were preliminarily considered the key genes for female flower development. This study has demonstrated that epigenetics was involved in the sex regulation of P. amurense, with DNA methylation as one of its regulatory modes. Moreover, some candidate genes related to the sexual differentiation of P. amurense were obtained with analysis. These results are of great significance for further exploring the mechanism of sex differentiation of P. amurense and studying of sex differentiation of plants.

Application of a novel strong promoter from Chinese fir (Cunninghamia lanceolate) in the CRISPR/Cas mediated genome editing of its protoplasts and transgenesis of rice and poplar

Novel constitutive promoters are essential for plant biotechnology. Although in angiosperms, a number of promoters were applied in monocots or dicots genetic engineering, only a few promoters were used in gymnosperm. Here we identified two strong promoters (Cula11 and Cula08) from Chinese fir (C. lanceolate) by screening the transcriptomic data and preliminary promoter activity assays in tobacco. By using the newly established Chinese fir protoplast transient expression technology that enables in vivo molecular biology studies in its homologous system, we compared the activities of Cula11 and Cula08 with that of the commonly used promoters in genetic engineering of monocots or dicots, such as CaM35S, CmYLCV, and ZmUbi, and our results revealed that Cula11 and Cula08 promoters have stronger activities in Chinese fir protoplasts. Furthermore, the vector containing Cas gene driven by Cula11 promoter and sgRNA driven by the newly isolated CulaU6b polyIII promoters were introduced into Chinese fir protoplasts, and CRISPR/Cas mediated gene knock-out event was successfully achieved. More importantly, compared with the commonly used promoters in the genetic engineering in angiosperms, Cula11 promoter has much stronger activity than CaM35S promoter in transgenic poplar, and ZmUbi promoter in transgenic rice, respectively, indicating its potential application in poplar and rice genetic engineering. Overall, the novel putative constitutive gene promoters reported here will have great potential application in gymnosperm and angiosperm biotechnology, and the transient gene expression system established here will serve as a useful tool for the molecular and genetic analyses of Chinese fir genes.

Genome-wide analysis of MADS-box gene family in kiwifruit (Actinidia chinensis var. chinensis) and their potential role in floral sex differentiation

Kiwifruit (Actinidia chinensis Planch.) is a functionally dioecious plant, which displays diverse morphology in male and female flowers. MADS-box is an ancient and huge gene family that plays a key role in plant floral organ differentiation. In this study, we have identified 89 MADS-box genes from A. chinensis Red 5 genome. These genes are distributed on 26 chromosomes and are classified into type I (21 genes) and type II (68 genes). Overall, type II AcMADS-box genes have more complex structures than type I with more exons, protein domains, and motifs, indicating that type II genes may have more diverse functions. Gene duplication analysis showed that most collinearity occurred in type II AcMADS-box genes, which was consistent with a large number of type II genes. Analysis of cis-acting elements in promoters showed that AcMADS-box genes are mainly associated with light and phytohormone responsiveness. The expression profile of AcMADS-box genes in different tissues showed that most genes were highly expressed in flowers. Further, the qRT-PCR analysis of the floral organ ABCDE model-related genes in male and female flowers revealed that AcMADS4, AcMADS56, and AcMADS70 were significantly expressed in female flowers. It indicated that those genes may play an important role in the sex differentiation of kiwifruit. This work provided a comprehensive analysis of the AcMADS-box genes and may help facilitate our understanding of the sex differentiation regulatory mechanism in kiwifruit.

Genome-wide identification and expression analysis of MIKCC genes in rose provide insight into their effects on flower development

The MIKC^C-type gene family plays important roles in plant growth, development, and tolerance of biotic and abiotic stress, especially during floral organ differentiation. However, there have been no studies of MIKC^C-type genes in rose, and functional differentiation of family members has not been explored. In this study, we identified 42 MIKC^C-type genes in rose, classified the genes into 12 subfamilies, and constructed a phylogenetic tree. We performed expression analysis of these genes, and found that expression patterns correlated with the predicted subfamily, indicating that the features of MIKC^C-type genes were broadly conserved during evolution. Collinear analysis of MIKC^C genes among Rosaceae species confirmed the occurrence of whole genome duplications (WGD) and revealed some species-specific MIKC^C genes. Transcriptome analysis showed that the expression of some MIKC^C-type genes responded to low temperatures (4°C, 24 h) during flower organ differentiation. These conserved, duplicated, and novel expression patterns of MIKC^C-type genes may have facilitated the adaptation of rose to various internal and external environmental changes. The results of this study provide a theoretical basis for future functional analysis of the MIKC^C genes in rose and investigation of the evolutionary pattern of the MIKC^C gene family in the Rosaceae genome.

Genome-Wide Identification and Comparative Analysis of WOX Genes in Four Euphorbiaceae Species and Their Expression Patterns in Jatropha curcas

The WUSCHEL-related homeobox (WOX) proteins are widely distributed in plants and play important regulatory roles in growth and development processes such as embryonic development and organ development. Here, series of bioinformatics methods were utilized to unravel the structural basis and genetic hierarchy of WOX genes, followed by regulation of the WOX genes in four Euphorbiaceae species. A genome-wide survey identified 59 WOX genes in Hevea brasiliensis (H. brasiliensis: 20 genes), Jatropha curcas (J. curcas: 10 genes), Manihot esculenta (M. esculenta: 18 genes), and Ricinus communis (R. communis: 11 genes). The phylogenetic analysis revealed that these WOX members could be clustered into three close proximal clades, such as namely ancient, intermediate and modern/WUS clades. In addition, gene structures and conserved motif analyses further validated that the WOX genes were conserved within each phylogenetic clade. These results suggested the relationships among WOX members in the four Euphorbiaceae species. We found that WOX genes in H. brasiliensis and M. esculenta exhibit close genetic relationship with J. curcas and R. communis. Additionally, the presence of various cis-acting regulatory elements in the promoter of J. curcas WOX genes (JcWOXs) reflected distinct functions. These speculations were further validated with the differential expression profiles of various JcWOXs in seeds, reflecting the importance of two JcWOX genes (JcWOX6 and JcWOX13) during plant growth and development. Our quantitative real-time PCR (qRT-PCR) analysis demonstrated that the JcWOX11 gene plays an indispensable role in regulating plant callus. Taken together, the present study reports the comprehensive characteristics and relationships of WOX genes in four Euphorbiaceae species, providing new insights into their characterization.