The genomic architecture of the sex-determining region and sex-related metabolic variation in Ginkgobiloba

Sexually dimorphic morphology, feeding behavior and gene expression profiles in cotton aphid Aphis gossypii

BACKGROUND

Sexual dimorphism exists in most insects; however, less is known about sexual dimorphism in aphids. In this study, we identified sexually dimorphic differences in morphology, feeding behavior and gene expression between sexual females and males of the cotton aphid through electron microscopy, electrical penetration graph techniques and RNA sequencing.

RESULTS

All males were alate with a slender reddish-yellow body and abdominal yellow-black stripes, whereas all sexual females were apterous with a pudgy green body. Sensillum types on the antennae were identical between the two sexes, although males had more sensilla, possibly because the antennae are significantly longer in males compared with sexual females. In terms of feeding behavior, males spent more time probing mesophyll cells and the phloem sieve, and salivating into the phloem sieve. By contrast, sexual females spent more time ingesting xylem sap. In total, 510 and 724 genes were specifically expressed in sexual females and males, respectively, and were significantly enriched in signaling pathways related to reproduction for sexual females (e.g. ovarian steroidogenesis, oxytocin signaling pathway) and energy and flight for males (e.g. thermogenesis, insulin signaling pathway). Moreover, 8551 differentially expressed genes were identified between the two sexes, of which the 3720 upregulated genes in sexual females were mostly enriched in signaling pathways of metabolism and energy, such as thermogenesis and the citrate cycle.

CONCLUSION

This study provides insight into sexual dimorphism in aphids and lays a foundation for revealing the molecular mechanism underlying differences between the two sexes in cotton aphid. © 2023 Society of Chemical Industry.

A rapid method for assembly of single chromosome and identification of sex determination region based on single-chromosome sequencing

The sex-determining-region (SDR) may offer the best prospects for studying sex-determining gene, recombination suppression, and chromosome heteromorphism. However, current progress of SDR identification and cloning showed following shortcomings: large near-isogenic lines need to be constructed, and a relatively large population is needed; the cost of whole-genome sequencing and assembly is high. Herein, the X/Y chromosomes of Spinacia oleracea L. subsp. turkestanica were successfully microdissected and assembled using single-chromosome sequencing. The assembly length of X and Y chromosome is c. 192.1 and 195.2 Mb, respectively. Three large inversions existed between X and Y chromosome. The SDR size of X and Y chromosome is c. 13.2 and 24.1 Mb, respectively. MSY region and six male-biased genes were identified. A Y-chromosome-specific marker in SDR was constructed and used to verify the chromosome assembly quality at cytological level via fluorescence in situ hybridization. Meanwhile, it was observed that the SDR located on long arm of Y chromosome and near the centromere. Overall, a technical system was successfully established for rapid cloning the SDR and it is also applicable to rapid assembly of specific chromosome in other plants. Furthermore, this study laid a foundation for studying the molecular mechanism of sex chromosome evolution in spinach.

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.

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.

Regulatory mechanism of MeGI on sexuality in Diospyros oleifera

Dioecy system is an important strategy for maintaining genetic diversity. The transcription factor MeGI, contributes to dioecy by promoting gynoecium development in Diospyros lotus and D. kaki. However, the function of MeGI in D. oleifera has not been identified. In this study, we confirmed that MeGI, cloned from D. oleifera, repressed the androecium development in Arabidopsis thaliana. Subsequently, chromatin immunoprecipitation-sequencing (ChIP-seq), DNA affinity purification-sequencing (DAP-seq), and RNA-seq were used to uncover the gene expression response to MeGI. The results showed that the genes upregulated and downregulated in response to MeGI were mainly enriched in the circadian rhythm-related and flavonoid biosynthetic pathways, respectively. Additionally, the WRKY DNA-binding protein 28 (WRKY28) gene, which was detected by ChIP-seq, DAP-seq, and RNA-seq, was emphasized. WRKY28 has been reported to inhibit salicylic acid (SA) biosynthesis and was upregulated in MeGI-overexpressing A. thaliana flowers, suggesting that MeGI represses the SA level by increasing the expression level of WRKY28. This was confirmed that SA level was lower in D. oleifera female floral buds than male. Overall, our findings indicate that the MeGI mediates its sex control function in D. oleifera mainly by regulating genes in the circadian rhythm, SA biosynthetic, and flavonoid biosynthetic pathways.

Genome-wide identification and expression analysis of MADS-box transcription factors reveal their involvement in sex determination of hardy rubber tree (Eucommia ulmoides oliv.)

Eucommia ulmoides is a famous rubber-producing and medicinal tree species that produces unisexual flowers on separate individuals from the earliest stage of stamen/pistil primordium formation. To explore the genetic regulation pathway of sex in E. ulmoides, comprehensive genome-wide analyses and tissue-/sex-specific transcriptome comparisons of MADS-box transcription factors were performed for the first time in this work. Quantitative real-time PCR technique was employed to further validate the expression of genes that are assigned to floral organ ABCDE model. A total of 66 non-redundant E. ulmoides MADS-box (EuMADS) genes were identified, they were classified into Type I (M-type, 17 genes) and Type II (MIKC, 49 genes). Complex protein-motif composition, exon-intron structure and phytohormone-response cis-elements were detected in MIKC-EuMADS genes. Furthermore, 24 differentially-expressed EuMADS genes (DEGs) between male and female flowers, and two DEGs between male and female leaves were revealed. Amongst the 14 floral organ ABCDE model-related genes, there were 6 (A/B/C/E-class) and 5 (A/D/E-class) genes displayed male- and female-biased expression respectively. In particular, one B-class gene EuMADS39 and one A-class gene EuMADS65 were almost exclusively expressed in male trees, no matter in flower or leaf tissues. Collectively, these results suggested a critical role of MADS-box transcription factors in sex determination of E. ulmoides, which is conducive to decoding the molecular regulation mechanism of sex in E. ulmoides.

Genome-wide identification, interaction of the MADS-box proteins in Zanthoxylum armatum and functional characterization of ZaMADS80 in floral development

As a typical dioecious species, Zanthoxylum armatum establishes apomictic reproduction, hence only female trees are cultivated. However, male and hermaphrodite flowers have recently appeared in female plants, resulting in a dramatic yield reduction. To date, the genetic basis underlying sex determination and apomixis in Z. armatum has been largely unknown. Here, we observed abortion of the stamen or carpel prior to primordium initiation, thus corroborating the potential regulation of MADS-box in sex determination. In Z. armatum, a total of 105 MADS-box genes were identified, harboring 86 MIKC-type MADSs with lack of FLC orthologues. Transcriptome analysis revealed candidate MADSs involved in floral organ identity, including ten male-biased MADSs, represented by ZaMADS92/81/75(AP3/PI-like), and twenty-six female-specified, represented by ZaMADS80/49 (STK/AGL11-like) and ZaMADS42 (AG-like). Overexpressing ZaMADS92 resulted in earlier flowering, while ZaMADS80 overexpression triggered precocious fruit set and parthenocarpy as well as dramatic modifications in floral organs. To characterize their regulatory mechanisms, a comprehensive protein-protein interaction network of the represented MADSs was constructed based on yeast two-hybrid and bimolecular fluorescence complementation assays. Compared with model plants, the protein interaction patterns in Z. armatum exhibited both conservation and divergence. ZaMADS70 (SEP3-like) interacted with ZaMADS42 and ZaMADS48 (AP3-like) but not ZaMADS40 (AP1-like), facilitating the loss of petals in Z. armatum. The ZaMADS92/ZaMADS40 heterodimer could be responsible for accelerating flowering in ZaMADS92-OX lines. Moreover, the interactions between ZaMADS80 and ZaMADS67(AGL32-like) might contribute to apomixis. This work provides new insight into the molecular mechanisms of MADS-boxes in sex organ identity in Z. armatum.

A convergent mechanism of sex determination in dioecious plants: Distinct sex-determining genes display converged regulation on floral B-class genes

Sex determination in dioecious plants has been broadly and progressively studied with the blooming of genome sequencing and editing techniques. This provides us with a great opportunity to explore the evolution and genetic mechanisms underlining the sex-determining system in dioecious plants. In this study, comprehensively reviewing advances in sex-chromosomes, sex-determining genes, and floral MADS-box genes in dioecious plants, we proposed a convergent model that governs plant dioecy across divergent species using a cascade regulation pathway connecting sex-determining genes and MADS-box genes e.g., B-class genes. We believe that this convergent mechanism of sex determination in dioecious plants will shed light on our understanding of gene regulation and evolution of plant dioecy. Perspectives concerning the evolutionary pathway of plant dioecy are also suggested.

Chromosome-level Genomes Reveal the Genetic Basis of Descending Dysploidy and Sex Determination in Morus Plants

Multiple plant lineages have independently evolved sex chromosomes and variable karyotypes to maintain their sessile lifestyles through constant biological innovation. Morus notabilis, a dioecious mulberry species, has the fewest chromosomes among Morus spp., but the genetic basis of sex determination and karyotype evolution in this species has not been identified. In this study, three high-quality genome assemblies were generated for Morus spp. [including dioecious M. notabilis (male and female) and Morus yunnanensis (female)] with genome sizes of 301-329 Mb and were grouped into six pseudochromosomes. Using a combination of genomic approaches, we found that the putative ancestral karyotype of Morus species was close to 14 protochromosomes, and that several chromosome fusion events resulted in descending dysploidy (2n = 2x = 12). We also characterized a ∼ 6.2-Mb sex-determining region on chromosome 3. Four potential male-specific genes, a partially duplicatedDNA helicase gene (named MSDH) and three Ty3_Gypsy long terminal repeat retrotransposons (named MSTG1/2/3), were identified in the Y-linked area and considered to be strong candidate genes for sex determination or differentiation. Population genomic analysis showed that Guangdong accessions in China were genetically similar to Japanese accessions of mulberry. In addition, genomic areas containing selective sweeps that distinguish domesticated mulberry from wild populations in terms of flowering and disease resistance were identified. Our findings provide an important genetic resource for sex identification research and molecular breeding in mulberry.

Evolution of the sex-determining region in Ginkgo biloba

Sex chromosomes or sex-determining regions (SDR) have been discovered in many dioecious plant species, including the iconic ‘living fossil' Ginkgo biloba, though the location and size of the SDR in G. biloba remain contradictory. Here we resolve these controversies and analyse the evolution of the SDR in this species. Based on transcriptome sequencing data from four genetic crosses we reconstruct male- and female-specific genetic maps and locate the SDR to the middle of chromosome 2. Integration of the genetic maps with the genome sequence reveals that recombination in and around the SDR is suppressed in a region of about 50 Mb in both males and females. However, occasional recombination does occur except a small, less than 5 Mb long region that does not recombine in males. Based on synonymous divergence between homologous X- and Y-linked genes in this region, we infer that the Ginkgo SDR is fairly old—at least of Cretaceous origin. The analysis of substitution rates and gene expression reveals only slight Y-degeneration. These results are consistent with findings in other dioecious plants with homomorphic sex chromosomes, where the SDR is typically small and evolves in a region with pre-existing reduced recombination, surrounded by long actively recombining pseudoautosomal regions.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

Integrative Analyses of Transcriptomics and Metabolomics in Sex Differentiation of Mulberry Flowers

Sex determination and sex differentiation of plants are important physiological processes of plant development. Mulberry (Morus indica L.) is an important economic tree being cultivated in sericulture countries, and mulberry leaf is commonly used for sericulture. The transcriptomic and metabolomic differences between the staminate flowers (SFs) and pistillate flowers (PFs) of mulberry were investigated by RNA sequencing and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Overall, we uncovered 4,230 genes and 209 metabolites are significantly differentially expressed between the SFs and PFs of mulberry. The combined transcriptomic and metabolomic analysis revealed these differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) are involved in flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism, and these detected DEGs and DEMs may be associated with sex differentiation of mulberry through the regulation of the enrichment pathways, such as the MAPK pathway, flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism. This study will provide a rich source for the analysis of the molecular mechanism of mulberry sex differentiation processes.

The Cycas genome and the early evolution of seed plants

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgo form a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycas genome contains four homologues of the fitD gene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensis contains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensis genome provides an important new resource of broad utility for biologists.

Role of female-predominant MYB39-bHLH13 complex in sexually dimorphic accumulation of taxol in Taxus media

Taxus trees are major natural sources for the extraction of taxol, an anti-cancer agent used worldwide. Taxus media is a dioecious woody tree with high taxol yield. However, the sexually dimorphic accumulation of taxoids in T. media is largely unknown. Our study revealed high accumulation of taxoids in female T. media trees using a UPLC-MS/MS method. Thereafter, many differential metabolites and genes between female and male T. media trees were identified using metabolomic and transcriptomic analyses, respectively. Most of the taxol-related genes were predominantly expressed in female trees. A female-specific R2R3-MYB transcription factor gene, TmMYB39, was identified. Furthermore, bimolecular fluorescence complementation and yeast two-hybrid assays suggested the potential interaction between TmMYB39 and TmbHLH13. Several taxol biosynthesis-related promoter sequences were isolated and used for the screening of MYB recognition elements. The electrophoretic mobility shift assay indicated that TmMYB39 could bind to the promoters of the GGPPS, T10OH, T13OH, and TBT genes. Interaction between TmMYB39 and TmbHLH13 transactivated the expression of the GGPPS and T10OH genes. TmMYB39 might function in the transcriptional regulation of taxol biosynthesis through an MYB-bHLH module. Our results give a potential explanation for the sexually dimorphic biosynthesis of taxol in T. media.

Early sex determination of Ginkgo biloba based on the differences in the electrocatalytic performance of extracted peroxidase

Ginkgo biloba is a dioecious plant. Male ginkgoes are mainly used in landscaping, while females are mainly used for fruit production. However, sex identification of ginkgo is a difficult task, especially at the seedling stage. In this work, we present for the first time the use of electrochemical techniques for the identification of ginkgo sex based on the differences in peroxides within male and female ginkgos. Graphene was used to concentrate peroxides in ginkgo extract, thereby improving electrochemical signal sensitivity. The electrochemical reduction of hydrogen peroxide catalyzed by peroxidase was used as a prob for sex determination in ginkgo. This electrochemical identification technique can be used not only for the analysis of adult ginkgo, but also successfully for the analysis of tissue culture seedlings and live seedlings. This electrochemical sensor has excellent discrimination ability due to the difference in peroxidase content in the leaves and petiole of ginkgo of different sexes. This electrochemical sensor allows for a rapid identification of the sex of ginkgo and has a very strong potential for field analysis.

The Diversity of Plant Sex Chromosomes Highlighted through Advances in Genome Sequencing

For centuries, scientists have been intrigued by the origin of dioecy in plants, characterizing sex-specific development, uncovering cytological differences between the sexes, and developing theoretical models. Through the invention and continued improvements in genomic technologies, we have truly begun to unlock the genetic basis of dioecy in many species. Here we broadly review the advances in research on dioecy and sex chromosomes. We start by first discussing the early works that built the foundation for current studies and the advances in genome sequencing that have facilitated more-recent findings. We next discuss the analyses of sex chromosomes and sex-determination genes uncovered by genome sequencing. We synthesize these results to find some patterns are emerging, such as the role of duplications, the involvement of hormones in sex-determination, and support for the two-locus model for the origin of dioecy. Though across systems, there are also many novel insights into how sex chromosomes evolve, including different sex-determining genes and routes to suppressed recombination. We propose the future of research in plant sex chromosomes should involve interdisciplinary approaches, combining cutting-edge technologies with the classics to unravel the patterns that can be found across the hundreds of independent origins.

Metabolome and Transcriptome Association Analysis Reveals Regulation of Flavonoid Biosynthesis by Overexpression of LaMIR166a in Larix kaempferi (Lamb.) Carr

Somatic embryogenesis is an ideal model process for studying early plant development. Embryonic cell lines of Larix kaempferi (Lamb.) Carr overexpressing LaMIR166a were obtained in our previous study. Here, a combination of de novo transcriptomics and extensively targeted metabolomics was used to study the transcriptional profiles and metabolic changes in wild-type and LaMIR166a-overexpressed embryonic cell lines. A total of 459 metabolites were found in the wild-type and transgenic cell lines. Compared to those in the wild-type cell lines, transcripts and metabolites were significantly altered in the LaMIR166a-overexpressed cell lines. Among differentially expressed genes (DEGs), phenylalanine and flavonoid synthesis genes were significantly enriched, and among differentially accumulated metabolites (DAMs), phenolic acids and flavonoids accumulated in particularly high amounts. Thus, the flavonoid biosynthetic pathway seems to be the most abundant pathway in response to LaMIR166a overexpression. Based on the Kyoto Encyclopedia of Genes and Genomes database, the association analysis of metabolome and transcriptome data showed that flavonoid biosynthesis and plant hormone signal transduction processes were significantly changed in miR166a-overexpression lines, suggesting that miR166 might be involved in these processes. The present study identified a number of potential metabolites associated with LaMIR166a overexpression, providing a significant foundation for a better understanding of the regulatory mechanisms underlying miR166.