Genomes of the Banyan Tree and Pollinator Wasp Provide Insights into Fig-Wasp Coevolution

Haplotype-resolved genome assembly of the upas tree (Antiaris toxicaria)

The upas tree (Antiaris toxicaria Lesch.) is a medically important plant that contains various specialized metabolites with significant bioactivity. The lack of a reference genome hinders the in-depth study as well as rational exploitation and conservation of this plant. Here, we present the first holotype-resolved chromosome-scale genome of the upas tree. The assembled genome consisted of 26 chromosomes that contain 1.34 Gb of sequencing data with a contig N50 length of 60 Mb. Genome annotation identified 43,500 protein-coding genes in the upas tree genome, of which 98.75% were functionally annotated. This high-quality reference genome will lay the foundation for further studies on the evolution and functional genomics of the upas tree.

Morphological and Transcriptomic Analyses Reveal the Involvement of Key Metabolic Pathways in Male Sterility in Chimonanthus praecox (L.) Genotypes

Chimonanthus praecox (Calycanthaceae family) is a unique ornamental and economic flowering tree in China, and after thousands of years of cultivation, it has produced several varieties and varietal types. Notably, male sterility is common in flowering plants and is an important tool for the genetic improvement in plants and optimization using hybrid plant technology; however, there have been no reports on male-sterile material or related studies on C. praecox. To our knowledge, this is the first time that C. praecox male sterility is dissected unveiling the involvement of key metabolic pathways. Notably, male sterility in C. praecox was observed during the budding period and likely occurred during the premature stage of pollen cell maturation. Additionally, differentially expressed genes in the starch and sucrose metabolism pathway and the plant hormone signal transduction pathway showed regular expression trends. This study reports on significant genetic differences that contribute to male sterility in C. praecox and provides a basis for further research and breeding strategies.

Comparative transcriptome analysis of differentially expressed genes and pathways in male and female flowers of Fraxinus mandshurica

Fraxinus mandshurica Rupr. (F. mandshurica) is a dioecious tree species with important ecological and application values. To delve deeper into the regulatory pathways and genes responsible for male and female flowers in F. mandshurica, we conducted transcriptome sequencing on male and female flowers at four distinct stages. The analysis revealed that the female database generated 38,319,967 reads while the male database generated 43,320,907 reads, resulting in 2930 differentially expressed genes with 1441 were up-regulated and 1489 down-regulated in males compared to females. Following an analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), four distinct pathways (hormone signal transduction, energy metabolism, flavonoid biosynthesis, and photoperiod) linked to female and male flowers were identified. Subsequently, qRT-PCR verification revealed that FmAUX/IAA, FmEIN3, and FmA-ARR genes in hormone signal transduction pathway are related to female flower development. Meanwhile, FmABF genes in hormone signal transduction pathway, FmGS and FmGDH genes in energy metabolism pathway, FmFLS genes in flavonoid biosynthesis pathway, and FmCaM, FmCRY, and FmPKA genes in photoperiod pathway are related to male flower development. This study was the first to analyze the transcriptome of male and female flowers of F. mandshurica, providing a reference for the developmental pathways and gene expression levels of male and female plants.

Putting the puzzle together: the relationship between floral characters and pollinator morphology determines pollination mode in the fig-fig wasp mutualism

The diversification of angiosperms has largely been attributed to adaptive radiation of their pollination and mating systems, which are relevant drivers of the macroevolution processes. The fig (Ficus, Moraceae) and fig wasp (Agaonidae, Hymenoptera) interaction is an example of obligate mutualism. Passive and active pollination modes have been associated with morphological traits in both partners. However, more information is required to assess the relationship between floral traits and pollination modes, particularly in Neotropical Ficus species. This study evaluates the morphological traits of figs and fig wasps regarding pollination modes in species belonging to Neotropical Ficus sections (three species each of Americanae and Pharmacosycea). Pollination mode was identified by floral morphology, anther/ovule ratio, and specialized structures fig wasps use for pollen transport (pollen pocket and coxal combs). Fig species in sect. Americanae are actively pollinated because pistillate flowers form a synstigma, present anther/ovule ratios <0.11, and their pollinator Pegoscapus fig wasps have pollen pockets and coxal combs. In contrast, species in sect. Pharmacosycea have free pistillate flowers, with anther/ovule ratios >0.27; they are pollinated by Tetrapus wasps, which lack specialized structures to carry pollen. Each species of Ficus was associated with a single morphospecies of fig wasp. The results support previous contributions that consider reciprocal morphological traits between fig species and their pollinating wasps as evidence of a close co-evolutionary history.

Phenotypic, Metabolic and Genetic Adaptations of the Ficus Species to Abiotic Stress Response: A Comprehensive Review

The Ficus genus, having radiated from the tropics and subtropics to the temperate zone worldwide, is the largest genus among woody plants, comprising over 800 species. Evolution of the Ficus species results in genetic diversity, global radiation and geographical differentiations, suggesting adaption to diverse environments and coping with stresses. Apart from familiar physiological changes, such as stomatal closure and alteration in plant hormone levels, the Ficus species exhibit a unique mechanism in response to abiotic stress, such as regulation of leaf temperature and retention of drought memory. The stress-resistance genes harbored by Ficus result in effective responses to abiotic stress. Understanding the stress-resistance mechanisms in Ficus provides insights into the genetic breeding toward stress-tolerant crop cultivars. Following upon these issues, we comprehensively reviewed recent progress concerning the Ficus genes and relevant mechanisms that play important roles in the abiotic stress responses. These highlight prospectively important application potentials of the stress-resistance genes in Ficus.

Molecular mechanism of vivipary as revealed by the genomes of viviparous mangroves and non-viviparous relatives

Vivipary is a prominent feature of mangroves, allowing seeds to complete germination while attached to the mother plant, and equips propagules to endure and flourish in challenging coastal intertidal wetlands. However, vivipary-associated genetic mechanisms remain largely elusive. Genomes of two viviparous mangrove species and a non-viviparous inland relative were sequenced and assembled at the chromosome level. Comparative genomic analyses between viviparous and non-viviparous genomes revealed that DELAY OF GERMINATION 1 (DOG1) family genes (DFGs), the proteins from which are crucial for seed dormancy, germination, and reserve accumulation, are either lost or dysfunctional in the entire lineage of true viviparous mangroves but are present and functional in their inland, non-viviparous relatives. Transcriptome dynamics at key stages of vivipary further highlighted the roles of phytohormonal homeostasis, proteins stored in mature seeds, and proanthocyanidins in vivipary under conditions lacking DFGs. Population genomic analyses elucidate dynamics of syntenic regions surrounding the missing DFGs. Our findings demonstrated the genetic foundation of constitutive vivipary in Rhizophoraceae mangroves.

The miR159a-DUO1 module regulates pollen development by modulating auxin biosynthesis and starch metabolism in citrus

Achieving seedlessness in citrus varieties is one of the important objectives of citrus breeding. Male sterility associated with abnormal pollen development is an important factor in seedlessness. However, our understanding of the regulatory mechanism underlying the seedlessness phenotype in citrus is still limited. Here, we determined that the miR159a-DUO1 module played an important role in regulating pollen development in citrus, which further indirectly modulated seed development and fruit size. Both the overexpression of csi-miR159a and the knocking out of DUO1 in Hong Kong kumquat (Fortunella hindsii) resulted in small and seedless fruit phenotypes. Moreover, pollen was severely aborted in both transgenic lines, with arrested pollen mitotic I and abnormal pollen starch metabolism. Through additional cross-pollination experiments, DUO1 was proven to be the key target gene for miR159a to regulate male sterility in citrus. Based on DNA affinity purification sequencing (DAP-seq), RNA-seq, and verified interaction assays, YUC2/YUC6, SS4 and STP8 were identified as downstream target genes of DUO1, those were all positively regulated by DUO1. In transgenic F. hindsii lines, the miR159a-DUO1 module down-regulated the expression of YUC2/YUC6, which decreased indoleacetic acid (IAA) levels and modulated auxin signaling to repress pollen mitotic I. The miR159a-DUO1 module reduced the expression of the starch synthesis gene SS4 and sugar transport gene STP8 to disrupt starch metabolism in pollen. Overall, this work reveals a new mechanism by which the miR159a-DUO1 module regulates pollen development and elucidates the molecular regulatory network underlying male sterility in citrus.

Plant regeneration in the new era: from molecular mechanisms to biotechnology applications

Plants or tissues can be regenerated through various pathways. Like animal regeneration, cell totipotency and pluripotency are the molecular basis of plant regeneration. Detailed systematic studies on Arabidopsis thaliana gradually unravel the fundamental mechanisms and principles underlying plant regeneration. Specifically, plant hormones, cell division, epigenetic remodeling, and transcription factors play crucial roles in reprogramming somatic cells and reestablishing meristematic cells. Recent research on basal non-vascular plants and monocot crops has revealed that plant regeneration differs among species, with various plant species using distinct mechanisms and displaying significant differences in regenerative capacity. Conducting multi-omics studies at the single-cell level, tracking plant regeneration processes in real-time, and deciphering the natural variation in regenerative capacity will ultimately help understand the essence of plant regeneration, improve crop regeneration efficiency, and contribute to future crop design.

Genome assembly and resequencing shed light on evolution, population selection, and sex identification in Vernicia montana

Vernicia montana is a dioecious plant widely cultivated for high-quality tung oil production and ornamental purposes in the Euphorbiaceae family. The lack of genomic information has severely hindered molecular breeding for genetic improvement and early sex identification in V. montana. Here, we present a chromosome-level reference genome of a male V. montana with a total size of 1.29 Gb and a contig N50 of 3.69 Mb. Genome analysis revealed that different repeat lineages drove the expansion of genome size. The model of chromosome evolution in the Euphorbiaceae family suggests that polyploidization-induced genomic structural variation reshaped the chromosome structure, giving rise to the diverse modern chromosomes. Based on whole-genome resequencing data and analyses of selective sweep and genetic diversity, several genes associated with stress resistance and flavonoid synthesis such as CYP450 genes and members of the LRR-RLK family, were identified and presumed to have been selected during the evolutionary process. Genome-wide association studies were conducted and a putative sex-linked insertion and deletion (InDel) (Chr 2: 102 799 917-102 799 933 bp) was identified and developed as a polymorphic molecular marker capable of effectively detecting the gender of V. montana. This InDel is located in the second intron of VmBASS4, suggesting a possible role of VmBASS4 in sex determination in V. montana. This study sheds light on the genome evolution and sex identification of V. montana, which will facilitate research on the development of agronomically important traits and genomics-assisted breeding.

A high-quality chromosome-level genome assembly of Ficus hirta

Ficus species (Moraceae) play pivotal roles in tropical and subtropical ecosystems. Thriving across diverse habitats, from rainforests to deserts, they harbor a multitude of mutualistic and antagonistic interactions with insects, nematodes, and pathogens. Despite their ecological significance, knowledge about the genomic background of Ficus remains limited. In this study, we report a chromosome-level reference genome of F. hirta, with a total size of 297.27 Mb, containing 28,625 protein-coding genes and 44.67% repeat sequences. These findings illuminate the genetic basis of Ficus responses to environmental challenges, offering valuable genomic resources for understanding genome size, adaptive evolution, and co-evolution with natural enemies and mutualists within the genus.

Characterisation of the fig-fig wasp holobiont

Plants and animals have long been considered distinct kingdoms, yet here a 'plant-animal' is described. An extraordinary symbiosis in which neither organism can reproduce without the other, the fig tree (Ficus) provides the habitat for its exclusive pollinator: the fig wasp (Agaonidae). Characterising the 'fig-fig wasp holobiont' acknowledges, for the first time, 'plant-animal symbiogenesis'.

Odorant-Binding and Chemosensory Proteins in Fig Wasps: Evolutionary Insights From Comparative Studies

Fig wasps (Agaonidae; Hymenoptera) are the only pollinating insects of fig trees (Ficus; Moraceae), forming the most closely and highly specific mutualism with the host. We used transcriptome sequences of 25 fig wasps from six genera to explore the evolution of key molecular components of fig wasp chemosensory genes: odorant-binding proteins (OBPs) and chemosensory proteins (CSPs). We identified a total 321 OBPs and 240 CSPs, with each species recording from 6 to 27 OBP genes and 6-19 CSP genes. 318 OBP genes are clustered into 17 orthologous groups and can be divided into two groups: PBP sensitive to pheromone and GOBP sensitive to general odor molecules, such as alcohols, esters, acids, ketones, and terpenoids. 240 CSP genes are clustered into 12 orthologous groups, which can be divided into three major groups and have functions, such as olfactory, tissue formation and/or regeneration, developmental, and some specific and unknown function. The gene sequences of most orthologous groups vary greatly among species and are consistent with the phylogenetic relationships between fig wasps. Strong purifying selection of both OBP and CSP genes was detected, as shown by low ω values. A positive selection was detected in one locus in CSP1. In conclusion, the evolution of chemosensory proteins OBPs and CSPs in fig wasps is relatively conservative, and they play an indispensable role in the life activities of fig wasps. Our results provide a starting point for understanding the molecular basis of the chemosensory systems of fig wasps.

A telomere-to-telomere reference genome of ficus (Ficus hispida) provides new insights into sex determination

A high-quality reference genome is indispensable for resolving biologically essential traits. Ficus hispida is a dioecious plant. A complete Ficus reference genome will be crucial for understanding their sex evolution and important biological characteristics, such as aerial roots, mutualistic symbiosis with ficus-wasps, and fruiting from old stems. Here, we generated a telomere-to-telomere (T2T) genome for F. hispida using PacBio HiFi and Oxford Nanopore Ultra-long sequencing technologies. The genome contiguity and completeness has shown improvement compared with the previously released genome, with the annotation of six centromeres and 28 telomeres. We have refined our previously reported 2-Mb male-specific region into a 7.2-Mb genomic region containing 51 newly predicted genes and candidate sex-determination genes AG2 and AG3. Many of these genes showed extremely low expression, likely attributed to hypermethylation in the gene body and promoter regions. Gene regulatory networks (GRNs) revealed that AG2 and AG3 are related to the regulation of stamen development in male flowers, while the AG1 gene is responsible for regulating female flowers' defense responses and secondary metabolite processes. Comparative analysis of GRNs showed that the NAC, WRKY, and MYB transcription factor families dominate the female GRN, whereas the MADS and MYB transcription factor families are prevalent in the male GRN.

Advancing understanding of Ficus carica: a comprehensive genomic analysis reveals evolutionary patterns and metabolic pathway insights

Ficus carica L. (dioecious), the most significant commercial species in the genus Ficus, which has been cultivated for more than 11,000 years and was one of the first species to be domesticated. Herein, we reported the most comprehensive F. carica genome currently. The contig N50 of the Orphan fig was 9.78 Mb, and genome size was 366.34 Mb with 13 chromosomes. Based on the high-quality genome, we discovered that F. carica diverged from Ficus microcarpa ~34 MYA, and a WGD event took place about 2─3 MYA. Throughout the evolutionary history of F. carica, chromosomes 2, 8, and 10 had experienced chromosome recombination, while chromosome 3 saw a fusion and fission. It is worth proposing that the chromosome 9 experienced both inversion and translocation, which facilitated the emergence of the F. carica as a new species. And the selections of F. carica for the genes of recombination chromosomal fragment are compatible with their goal of domestication. In addition, we found that the F. carica has the FhAG2 gene, but there are structural deletions and positional jumps. This gene is thought to replace the one needed for female common type F. carica to be pollinated. Subsequently, we conducted genomic, transcriptomic, and metabolomic analysis to demonstrate significant differences in the expression of CHS among different varieties of F. carica. The CHS playing an important role in the anthocyanin metabolism pathway of F. carica. Moreover, the CHS gene of F. carica has a different evolutionary trend compared to other Ficus species. These high-quality genome assembly, transcriptomic, and metabolomic resources further enrich F. carica genomics and provide insights for studying the chromosomes evolution, sexual system, and color characteristics of Ficus.

Comparative Genomic Analysis of Asian Cultivated Rice and Its Wild Progenitor (Oryza rufipogon) Has Revealed Evolutionary Innovation of the Pentatricopeptide Repeat Gene Family through Gene Duplication

The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in land plants. However, current knowledge about the evolution of the PPR gene family remains largely limited. In this study, we performed a comparative genomic analysis of the PPR gene family in O. sativa and its wild progenitor, O. rufipogon, and outlined a comprehensive landscape of gene duplications. Our findings suggest that the majority of PPR genes originated from dispersed duplications. Although segmental duplications have only expanded approximately 11.30% and 13.57% of the PPR gene families in the O. sativa and O. rufipogon genomes, we interestingly obtained evidence that segmental duplication promotes the structural diversity of PPR genes through incomplete gene duplications. In the O. sativa and O. rufipogon genomes, 10 (~33.33%) and 22 pairs of gene duplications (~45.83%) had non-PPR paralogous genes through incomplete gene duplication. Segmental duplications leading to incomplete gene duplications might result in the acquisition of domains, thus promoting functional innovation and structural diversification of PPR genes. This study offers a unique perspective on the evolution of PPR gene structures and underscores the potential role of segmental duplications in PPR gene structural diversity.

The genetic structuring in pollinating wasps of Ficus hispida in continental Asia

The interaction between figs and fig wasps provides a striking example of obligate brood site pollination mutualism. Monoecious figs, constituting independent radiations in each tropical biome, are present in significant proportions worldwide, but in continental Asia, dioecious figs have diverged into various niches, making the region's assemblage remarkably diverse. However, the reproductive success of figs and fig wasps largely depends on the fig wasp dispersal process. Monoecious fig pollinators in continental Asian tropical rain forests exhibit high gene flow of the plant, while many dioecious fig pollinators have a more restricted gene flow. However, there are limited studies on the genetic structure of dioecious Ficus pollinators that pollinate figs with intermediate gene flow. Here, we used molecular methods to investigate the genetic structure of pollinating wasps of the widely distributed dioecious Ficus hispida in China and Southeast Asia. Sequence data from two gene regions were used: the mitochondrial protein-coding gene cytochrome c oxidase subunit I (COI) and the nuclear 28S genes. Both molecular and morphological results support two fig wasp species at our sampling sites. Our findings suggest that for widely sympatric Ficus species in continental Asia, monoecious figs presenting long gene glow have the fewest fig wasp species, followed by dioecious figs presenting intermediate gene flow, and dioecious figs presenting local gene flow have the most fig wasp species.

Specific volatiles of tea plants determine the host preference behavior of Empoasca onukii

Empoasca onukii is a major pest that attacks tea plants. To seek effective and sustainable methods to control the pest, it is necessary to assess its host preference among different species of tea and understand the critical factors behind this behavior. In this study, the behavioral preference of E. onukii for volatile organic compounds (VOCs) of three potted tea species was evaluated. The VOCs released by the three tea species were analyzed using gas chromatography-mass spectrometry, and the major components were used to test the pest's preference. Transcriptome analysis was used to infer the key genes that affect the biosyntheses of the VOCs. The results showed that the tendency of E. onukii toward the VOCs of the three tea species was the strongest in green tea, followed by white tea, and the weakest in red tea. This behavioral preference was significantly and positively correlated with the relative levels of hexanol, linalool, and geraniol in tea volatiles. Relative hexanol was significantly and positively correlated with the expression of genes TEA009423 (LOX2.1), TEA009596 (LOX1.5), TEA008699 (HPL), TEA018669 (CYPADH), and TEA015686 (ADHIII). Relative linalool was significantly and positively correlated with the expression of genes TEA001435 (CAD) and Camellia_sinensis_newGene_22126 (TPS). Relative geraniol was significantly and positively correlated with the expression of genes TEA001435 (CAD), TEA002658 (CYP76B6), TEA025455 (CYP76T24), and Camellia_sinensis_newGene_22126 (TPS). The above findings suggested that three volatiles (hexanol, linalool, and geraniol) determined the behavioral preference of E. onukii toward tea plants, and their biosynthesis was mainly affected by nine genes (TEA009423, TEA009596, TEA008699, TEA018669, TEA015686, TEA001435, TEA002658, TEA025455, and Camellia_sinensis_newGene_22126).

The identification and expression pattern of the sex determination genes and their sex-specific variants in the egg parasitoid Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae)

Introduction: Trichogramma wasps are egg parasitoids of agricultural lepidopteran pests. The sex of Trichogramma is determined by its ploidy as well as certain sex ratio distorters, such as the endosymbiotic bacteria Wolbachia spp. and the paternal sex ratio (PSR) chromosome. The sex determination systems of hymenopterans, such as Trichogramma spp., involve cascades of the genes transformer (tra), transformer-2 (tra2), and doublesex (dsx) and are associated with sex-specific tra and dsx splicing. First, these genes and their sex-specific variants must be identified to elucidate the interactions between the sex ratio disorders and the sex determination mechanism of Trichogramma. Methods: Here, we characterized the sex determination genes tra, tra2, and dsx in Trichogramma dendrolimi. Sex-specific tra and dsx variants were detected in cDNA samples obtained from both male and female Trichogramma wasps. They were observed in the early embryos (1-10 h), late embryos (12-20 h), larvae (32 h and 48 h), pre-pupae (96 h), and pupae (144 h, 168 h, 192 h, and 216 h) of both male and female T. dendrolimi offspring. Results: We detected female-specific tra variants throughout the entire early female offspring stage. The male-specific variant began to express at 9-10 h as the egg was not fertilized. However, we did not find any maternally derived, female-specific tra variant in the early male embryo. This observation suggests that the female-specific tra variant expressed in the female embryo at 1-9 h may not have originated from the maternal female wasp. Discussion: The present study might be the first to identify the sex determination genes and sex-specific gene splicing in Trichogramma wasps. The findings of this study lay the foundation for investigating the sex determination mechanisms of Trichogramma and other wasps. They also facilitate sex identification in immature T. dendrolimi and the application of this important egg parasitoid in biological insect pest control programs.

Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Echoes of ancient introgression punctuate stable genomic lineages in the evolution of figs

Studies investigating the evolution of flowering plants have long focused on isolating mechanisms such as pollinator specificity. Some recent studies have proposed a role for introgressive hybridization between species, recognizing that isolating processes such as pollinator specialization may not be complete barriers to hybridization. Occasional hybridization may therefore lead to distinct yet reproductively connected lineages. We investigate the balance between introgression and reproductive isolation in a diverse clade using a densely sampled phylogenomic study of fig trees (Ficus, Moraceae). Codiversification with specialized pollinating wasps (Agaonidae) is recognized as a major engine of fig diversity, leading to about 850 species. Nevertheless, some studies have focused on the importance of hybridization in Ficus, highlighting the consequences of pollinator sharing. Here, we employ dense taxon sampling (520 species) throughout Moraceae and 1,751 loci to investigate phylogenetic relationships and the prevalence of introgression among species throughout the history of Ficus. We present a well-resolved phylogenomic backbone for Ficus, providing a solid foundation for an updated classification. Our results paint a picture of phylogenetically stable evolution within lineages punctuated by occasional local introgression events likely mediated by local pollinator sharing, illustrated by clear cases of cytoplasmic introgression that have been nearly drowned out of the nuclear genome through subsequent lineage fidelity. The phylogenetic history of figs thus highlights that while hybridization is an important process in plant evolution, the mere ability of species to hybridize locally does not necessarily translate into ongoing introgression between distant lineages, particularly in the presence of obligate plant-pollinator relationships.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

Whole-genome Duplication Reshaped Adaptive Evolution in A Relict Plant Species, Cyclocarya paliurus

Cyclocarya paliurus is a relict plant species that survived the last glacial period and shows a population expansion recently. Its leaves have been traditionally used to treat obesity and diabetes with the well-known active ingredient cyclocaric acid B. Here, we presented three C. paliurus genomes from two diploids with different flower morphs and one haplotype-resolved tetraploid assembly. Comparative genomic analysis revealed two rounds of recent whole-genome duplication events and identified 691 genes with dosage effects that likely contribute to adaptive evolution through enhanced photosynthesis and increased accumulation of triterpenoids. Resequencing analysis of 45 C. paliurus individuals uncovered two bottlenecks, consistent with the known events of environmental changes, and many selectively swept genes involved in critical biological functions, including plant defense and secondary metabolite biosynthesis. We also proposed the biosynthesis pathway of cyclocaric acid B based on multi-omics data and identified key genes, in particular gibberellin-related genes, associated with the heterodichogamy in C. paliurus species. Our study sheds light on evolutionary history of C. paliurus and provides genomic resources to study the medicinal herbs.

Parasitism Features of a Fig Wasp of Genus Apocrypta (Pteromalidae: Pteromalinae) Associated with a Host Belonging to Ficus Subgenus Ficus

Non-pollinating fig wasps (NPFWs), particularly long-ovipositored Sycoryctina wasps, exhibit a high species specificity and exert complex ecological effects on the obligate mutualism between the plant genus Ficus and pollinating fig wasps. Apocrypta is a genus of NPFWs that mostly interacts with the Ficus species under the subgenus Sycomorus, and the symbiosis case between Apocrypta and F. pedunculosa var. mearnsii, a Ficus species under subgenus Ficus, is unique. As fig's internal environments and the wasp communities are distinct between the two subgenera, we addressed the following two questions: (1) Are the parasitism features of the Apocrypta wasp associated with F. pedunculosa var. mearnsii different from those of other congeneric species? (2) Is this Apocrypta species an efficient wasp that lives in its unique host? Our observation revealed that this wasp is an endoparasitic idiobiont parasitoid, as most congeneric species are, but developed a relatively long ovipositor. Furthermore, the relationships of the parasitism rate versus the pollinator number, the fig wall, and the sex ratio of the pollinator, respectively, showed that it possessed a higher parasitism ability than that of other congeners. However, its parasitism rate was low, and thus it was not an efficient wasp in its habitat. This difference between parasitism ability and parasitism rate might be a consequence of its oviposition strategy and the severe habitat conditions. These findings may also provide insights into the mechanism to maintain the interaction between the fig tree and the fig wasp community.

Genomic and transcriptomic analysis of sacred fig (Ficus religiosa)

BACKGROUND

Peepal/Bodhi tree (Ficus religiosa L.) is an important, long-lived keystone ecological species. Communities on the Indian subcontinent have extensively employed the plant in Ayurveda, traditional medicine, and spiritual practices. The Peepal tree is often thought to produce oxygen both during the day and at night by Indian folks. The goal of our research was to produce molecular resources using whole-genome and transcriptome sequencing techniques.

RESULTS

The complete genome of the Peepal tree was sequenced using two next-generation sequencers Illumina HiSeq1000 and MGISEQ-2000. We assembled the draft genome of 406 Mb, using a hybrid assembly workflow. The genome annotation resulted in 35,093 protein-coding genes; 53% of its genome consists of repetitive sequences. To understand the physiological pathways in leaf tissues, we analyzed photosynthetically distinct conditions: bright sunny days and nights. The RNA-seq analysis supported the expression of 26,479 unigenes. The leaf transcriptomic analysis of the diurnal and nocturnal periods revealed the expression of the significant number of genes involved in the carbon-fixation pathway.

CONCLUSIONS

This study presents a draft hybrid genome assembly for F. religiosa and its functional annotated genes. The genomic and transcriptomic data-derived pathways have been analyzed for future studies on the Peepal tree.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Olfactory and gustatory receptor genes in fig wasps: Evolutionary insights from comparative studies

The mechanisms of chemoreception in fig wasps (Hymenoptera, Agaonidae) are of primary importance in their co-evolutionary relationship with the fig trees they pollinate. We used transcriptome sequences of 25 fig wasps in six genera that allowed a comparative approach to the evolution of key molecular components of fig wasp chemoreception: their odorant (OR) and gustatory (GR) receptor genes. In total, we identified 311 ORs and 47 GRs, with each species recording from 5 to 30 OR genes and 1-4 GR genes. 304 OR genes clustered into 18 orthologous groups known to be sensitive to cuticular hydrocarbons (CHC), pheromones, acids, alcohols and a variety of floral scents such as cineole, Linalool, and Heptanone. 45 GR genes clustered into 4 orthologous groups that contain sweet, bitter, CO₂ and undocumented receptors. Gene sequences in most orthologous groups varied greatly among species, except for ORco (60.0% conserved) and sweet receptors (30.7% conserved). Strong purifying selection of both odorant and gustatory genes was detected, as shown by low ω values. Signatures of positive selection were detected in loci from both OR and GR orthologous groups. Fig wasps have relatively few olfactory and especially gustatory receptors, reflecting the natural history of the system. Amino acid sequences nonetheless vary significantly between species and are consistent with the phylogenetic relationships among fig wasps. The differences in ORs within some orthologous groups from the same species, but different hosts and from closely related species from one host can reach as low as 49.3% and 9.8% respectively, implying the ORs of fig wasps can evolve rapidly to novel ecological environments. Our results provide a starting point for understanding the molecular basis of the chemosensory systems of fig wasps.

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.

WOX11: the founder of plant organ regeneration

De novo organ regeneration is the process in which adventitious roots or shoots regenerate from detached or wounded organs. De novo organ regeneration can occur either in natural conditions, e.g. adventitious root regeneration from the wounded sites of detached leaves or stems, or in in-vitro tissue culture, e.g. organ regeneration from callus. In this review, we summarize recent advances in research on the molecular mechanism of de novo organ regeneration, focusing on the role of the WUSCHEL-RELATED HOMEOBOX11 (WOX11) gene in the model plant Arabidopsis thaliana. WOX11 is a direct target of the auxin signaling pathway, and it is expressed in, and regulates the establishment of, the founder cell during de novo root regeneration and callus formation. WOX11 activates the expression of its target genes to initiate root and callus primordia. Therefore, WOX11 links upstream auxin signaling to downstream cell fate transition during regeneration. We also discuss the role of WOX11 in diverse species and its evolution in plants.

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.

Comparison of chloroplast genomes and phylogenomics in the Ficus sarmentosa complex (Moraceae)

Due to maternal inheritance and minimal rearrangement, the chloroplast genome is an important genetic resource for evolutionary studies. However, the evolutionary dynamics and phylogenetic performance of chloroplast genomes in closely related species are poorly characterized, particularly in taxonomically complex and species-rich groups. The taxonomically unresolved Ficus sarmentosa species complex (Moraceae) comprises approximately 20 taxa with unclear genetic background. In this study, we explored the evolutionary dynamics, hotspot loci, and phylogenetic performance of thirteen chloroplast genomes (including eleven newly obtained and two downloaded from NCBI) representing the F. sarmentosa complex. Their sequence lengths, IR boundaries, repeat sequences, and codon usage were compared. Both sequence length and IR boundaries were found to be highly conserved. All four categories of long repeat sequences were found across all 13 chloroplast genomes, with palindromic and forward sequences being the most common. The number of simple sequence repeat (SSR) loci varied from 175 (F. dinganensis and F. howii) to 190 (F. polynervis), with the dinucleotide motif appearing the most frequently. Relative synonymous codon usage (RSCU) analysis indicated that codons ending with A/T were prior to those ending with C/T. The majority of coding sequence regions were found to have undergone negative selection with the exception of ten genes (accD, clpP, ndhK, rbcL, rpl20, rpl22, rpl23, rpoC1, rps15, and rps4) which exhibited potential positive selective signatures. Five hypervariable genic regions (rps15, ycf1, rpoA, ndhF, and rpl22) and five hypervariable intergenic regions (trnH-GUG-psbA, rpl32-trnL-UAG, psbZ-trnG-GCC, trnK-UUU-rps16 and ndhF-rpl32) were identified. Overall, phylogenomic analysis based on 123 Ficus chloroplast genomes showed promise for studying the evolutionary relationships in Ficus, despite cyto-nuclear discordance. Furthermore, based on the phylogenetic performance of the F. sarmentosa complex and F. auriculata complex, the chloroplast genome also exhibited a promising phylogenetic resolution in closely related species.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Jack of all trades: Genome assembly of Wild Jack and comparative genomics of Artocarpus

Artocarpus (Moraceae), known as breadfruits for their diverse nutritious fruits, is prized for its high-quality timber, medicinal value, and economic importance. Breadfruits are native to Southeast Asia but have been introduced to other continents. The most commonly cultivated species are Artocarpus heterophyllus (Jackfruit) and Artocarpus altilis (Breadfruit). With numerous smaller but nutritionally comparable fruits on a larger tree, Artocarpus hirsutus, also called "Wild Jack" or "Ayani", is an elusive forest species endemic to Indian Western Ghats. In this study, we sequenced and assembled the whole genome of Artocarpus hirsutus sampled from the sacred groves of Coorg, India. To decipher demographic and evolutionary history, we compared our Wild Jack genome with previously published Jackfruit and Breadfruit genomes. Demographic history reconstruction indicates a stronger effect of habitat rather than phylogeny on the population histories of these plants. Repetitive genomic regions, especially LTR Copia, strongly affected the demographic trajectory of A. heterophyllus. Upon further investigation, we found a recent lineage-specific accumulation of LTR Copia in A. heterophyllus, which had a major contribution to its larger genome size. Several genes from starch, sucrose metabolism, and plant hormone signal transduction pathways, in Artocarpus species had signatures of selection and gene family evolution. Our comparative genomic framework provides important insights by incorporating endemic species such as the Wild Jack.

Diversity and metabolic potentials of microbial communities associated with pollinator and cheater fig wasps in fig-fig wasp mutualism system

Microbial symbionts can influence a myriad of insect behavioral and physiological traits. However, how microbial communities may shape or be shaped by insect interactions with plants and neighboring species remains underexplored. The fig-fig wasp mutualism system offers a unique model to study the roles of microbiome in the interactions between the plants and co-habiting insects because a confined fig environment is shared by two fig wasp species, the pollinator wasp (Eupristina altissima and Eupristina verticillata) and the cheater wasp (Eupristina sp1 and Eupristina sp2). Here, we performed whole genome resequencing (WGS) on 48 individual fig wasps (Eupristina spp.) from Yunnan, China, to reveal the phylogenetic relationship and genetic divergence between pollinator and congeneric cheater wasps associated with the Ficus trees. We then extracted metagenomic sequences to explore the compositions, network structures, and functional capabilities of microbial communities associated with these wasps. We found that the cheaters and pollinators from the same fig species are sister species, which are highly genetically divergent. Fig wasps harbor diverse but stable microbial communities. Fig species dominate over the fig wasp genotype in shaping the bacterial and fungal communities. Variation in microbial communities may be partially explained by the filtering effect from fig and phylogeny of fig wasps. It is worth noting that cheaters have similar microbial communities to their sister pollinators, which may allow cheaters to coexist and gain resources from the same fig species. In terms of metabolic capabilities, some bacteria such as Desulfovibrio and Lachnospiraceae are candidates involved in the nutritional uptake of fig wasps. Our results provide novel insights into how microbiome community and metabolic functions may couple with the fig-wasp mutualistic systems.

Multi-omics analysis the differences of VOCs terpenoid synthesis pathway in maintaining obligate mutualism between Ficus hirta Vahl and its pollinators

INRODUCTION

Volatile organic compounds (VOCs) emitted by the receptive syconia of Ficus species is a key trait to attract their obligate pollinating fig wasps. Ficus hirta Vahl is a dioecious shrub, which is pollinated by a highly specialized symbiotic pollinator in southern China. Terpenoids are the main components of VOCs in F. hirta and play ecological roles in pollinator attraction, allelopathy, and plant defense. However, it remains unclear that what molecular mechanism difference in terpenoid synthesis pathways between pre-receptive stage (A-phase) and receptive stage (B-phase) of F. hirta syconia.

METHODS

Transcriptome, proteome and Gas Chromatography-Mass Spectrometer (GC-MS) were applied here to analyze these difference.

RESULTS AND DISCUSSION

Compared to A-phase syconia, the genes (ACAT2, HMGR3, GGPS2, HDR, GPS2, TPS2, TPS4, TPS10-4, TPS14) related to the terpenoid synthesis pathway had higher expression level in receptive syconia (B-phase) according to transcriptome sequencing. Seven differentially expressed transcription factors were screened, namely bHLH7, MYB1R1, PRE6, AIL1, RF2b, ANT, VRN1. Specifically, bHLH7 was only specifically expressed in B-phase. 235 differentially expressed proteins (DEPs) were mainly located in the cytoplasm and chloroplasts. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the DEPs were mainly enriched in the metabolic process. A total of 9 terpenoid synthesis proteins were identified in the proteome. Among them, 4 proteins in methylerythritol phosphate (MEP) pathway were all down-regulated. Results suggested the synthesis of terpenoids precursors in B-phase bracts were mainly accomplished through the mevalonic acid (MVA) pathway in cytoplasm. Correlation analysis between the transcriptome and proteome, we detected a total of 1082 transcripts/proteins, three of which are related to stress. From the VOCs analysis, the average percent of monoterpenoids emitted by A-phase and B-phase syconia were 8.29% and 37.08%, while those of sesquiterpenes were 88.43% and 55.02% respectively. Monoterpenes (camphene, myrcene, camphor, menthol) were only detected in VOCs of B-phase syconia. To attract pollinators, B-phase syconia of F. hirta need more monoterpenoids and less sesquiterpenes. We speculate that transcription factor bHLH7 may regulate the terpenoid synthesis pathway between A- and B-phase syconia. Our research provided the first global analysis of mechanism differences of terpenoid synthesis pathways between A and B phases in F. hirta syconia.

AP2/ERF genes associated with superfast fig (Ficus carica L.) fruit ripening

Fig fruits have significant health value and are culturally important. Under suitable climatic conditions, fig fruits undergo a superfast ripening process, nearly doubling in size, weight, and sugar content over three days in parallel with a sharp decrease in firmness. In this study, 119 FcAP2/ERF genes were identified in the fig genome, namely 95 ERFs, 20 AP2s, three RAVs, and one soloist. Most of the ERF subfamily members (76) contained no introns, whereas the majority of the AP2 subfamily members had at least two introns each. Three previously published transcriptome datasets were mined to discover expression patterns, encompassing the fruit peel and flesh of the 'Purple Peel' cultivar at six developmental stages; the fruit receptacle and flesh of the 'Brown Turkey' cultivar after ethephon treatment; and the receptacle and flesh of parthenocarpic and pollinated fruits of the 'Brown Turkey' cultivar. Eighty-three FcAP2/ERFs (68 ERFs, 13 AP2s, one RAV, and one soloist) were expressed in the combined transcriptome dataset. Most FcAP2/ERFs were significantly downregulated (|log₂(fold change) | ≥ 1 and p-adjust < 0.05) during both normal fruit development and ethephon-induced accelerated ripening, suggesting a repressive role of these genes in fruit ripening. Five significantly downregulated ERFs also had repression domains in the C-terminal. Seven FcAP2/ERFs were identified as differentially expressed during ripening in all three transcriptome datasets. These genes were strong candidates for future functional genetic studies to elucidate the major FcAP2/ERF regulators of the superfast fig fruit ripening process.

Genome sequencing and comparative analysis of Ficus benghalensis and Ficus religiosa species reveal evolutionary mechanisms of longevity

Ficus benghalensis and Ficus religiosa are large woody trees well known for their long lifespan, ecological and traditional significance, and medicinal properties. To understand the genomic and evolutionary aspects of these characteristics, the whole genomes of these Ficus species were sequenced using 10x Genomics linked reads and Oxford Nanopore long reads. The draft genomes of F. benghalensis and F. religiosa comprised of 392.89 Mbp and 332.97 Mbp, respectively. We established the genome-wide phylogenetic positions of the two Ficus species with respect to 50 other Angiosperm species. Comparative evolutionary analyses with other phylogenetically closer Eudicot species revealed adaptive evolution in genes involved in key cellular mechanisms associated with prolonged survival including phytohormones signaling, senescence, disease resistance, and abiotic stress tolerance, which provide genomic insights into the mechanisms conferring longevity and suggest that longevity is a multifaceted phenomenon. This study also provides clues on the existence of CAM pathway in these Ficus species.

Hemiepiphytic figs kill their host trees: acquiring phosphorus is a driving factor

Hemiepiphytic figs killing their host trees is an ecological process unique to the tropics. Yet the benefits and adaptive strategies of their special life history remain poorly understood. We compared leaf phosphorus (P) content data of figs and palms worldwide, and functional traits and substrate P content of hemiepiphytic figs (Ficus tinctoria), their host palm and nonhemiepiphytic conspecifics at different growth stages in a common garden. We found that leaf P content of hemiepiphytic figs and their host palms significantly decreased when they were competing for soil resources, but that of hemiepiphytic figs recovered after host death. P availability in the canopy humus and soil decreased significantly with the growth of hemiepiphytic figs. Functional trait trade-offs of hemiepiphytic figs enabled them to adapt to the P shortage while competing with their hosts. From the common garden to a global scale, the P competition caused by high P demand of figs may be a general phenomenon. Our results suggest that P competition is an important factor causing host death, except for mechanically damaging and shading hosts. Killing hosts benefits hemiepiphytic figs by reducing interspecific P competition and better acquiring P resources in the P-deficient tropics, thereby linking the life history strategy of hemiepiphytic figs to the widespread P shortage in tropical soils.

Comparative chloroplast genome analysis of Ficus (Moraceae): Insight into adaptive evolution and mutational hotspot regions

As the largest genus in Moraceae, Ficus is widely distributed across tropical and subtropical regions and exhibits a high degree of adaptability to different environments. At present, however, the phylogenetic relationships of this genus are not well resolved, and chloroplast evolution in Ficus remains poorly understood. Here, we sequenced, assembled, and annotated the chloroplast genomes of 10 species of Ficus, downloaded and assembled 13 additional species based on next-generation sequencing data, and compared them to 46 previously published chloroplast genomes. We found a highly conserved genomic structure across the genus, with plastid genome sizes ranging from 159,929 bp (Ficus langkokensis) to 160,657 bp (Ficus religiosa). Most chloroplasts encoded 113 unique genes, including a set of 78 protein-coding genes, 30 transfer RNA (tRNA) genes, four ribosomal RNA (rRNA) genes, and one pseudogene (infA). The number of simple sequence repeats (SSRs) ranged from 67 (Ficus sagittata) to 89 (Ficus microdictya) and generally increased linearly with plastid size. Among the plastomes, comparative analysis revealed eight intergenic spacers that were hotspot regions for divergence. Additionally, the clpP, rbcL, and ccsA genes showed evidence of positive selection. Phylogenetic analysis indicated that none of the six traditionally recognized subgenera of Ficus were monophyletic. Divergence time analysis based on the complete chloroplast genome sequences showed that Ficus species diverged rapidly during the early to middle Miocene. This research provides basic resources for further evolutionary studies of Ficus.

Genomic insights into rapid speciation within the world's largest tree genus Syzygium

Species radiations, despite immense phenotypic variation, can be difficult to resolve phylogenetically when genetic change poorly matches the rapidity of diversification. Genomic potential furnished by palaeopolyploidy, and relative roles for adaptation, random drift and hybridisation in the apportionment of genetic variation, remain poorly understood factors. Here, we study these aspects in a model radiation, Syzygium, the most species-rich tree genus worldwide. Genomes of 182 distinct species and 58 unidentified taxa are compared against a chromosome-level reference genome of the sea apple, Syzygium grande. We show that while Syzygium shares an ancient genome doubling event with other Myrtales, little evidence exists for recent polyploidy events. Phylogenomics confirms that Syzygium originated in Australia-New Guinea and diversified in multiple migrations, eastward to the Pacific and westward to India and Africa, in bursts of speciation visible as poorly resolved branches on phylogenies. Furthermore, some sublineages demonstrate genomic clines that recapitulate cladogenetic events, suggesting that stepwise geographic speciation, a neutral process, has been important in Syzygium diversification.

Transcriptomic profiles of poplar (Populus simonii × P. nigra) cuttings during adventitious root formation

The formation of adventitious roots (ARs) is vital for the vegetative propagation of poplars. However, the relevant mechanisms remain unclear. To reveal the underlying molecular mechanism, we used RNA-seq to investigate the transcriptional alterations of poplar cuttings soaked in water for 0, 2, 4, 6, 8, and 10 d; 3,798 genes were differentially expressed at all the time points, including 2,448 upregulated and 1,350 downregulated genes. Biological processes including “cell cycle,” “photosynthesis,” “regulation of hormone levels,” and “auxin transport” were enriched in the differentially expressed genes (DEGs). KEGG results showed that the common DEGs were most enriched in the pathway of “Carbon fixation in photosynthetic organisms” and “Starch and sucrose metabolism.” We further dissected 38 DEGs related to root and auxin, including two lateral root primordium 1 (LRP1), one root meristem growth factor (RGF9), one auxin-induced in the root (AIR12), three rooting-associated genes (AUR1 and AUR3), eight auxin transcription factors (ARFs and LBDs), 10 auxin respective genes (SAURs and GH3s), nine auxin transporters (PINs, ABCs, LAX2, and AUXs), and four auxin signal genes (IAAs and TIR1). We found that the rooting abilities of poplar cuttings with and without leaves are different. By applying different concentrations of IBA and sucrose to the top of cuttings without leaves, we found that 0.2 mg/ml IBA and 2 mg/ml sucrose had the best effect on promoting AR formation. The transcriptome results indicated photosynthesis may influence AR formation in poplar cuttings with leaves and revealed a potential regulatory mechanism of leafy cuttage from poplar cuttings. In addition, we provided a new perspective to resolve rooting difficulties in recalcitrant species.

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.

Identification and characterization of auxin response factor (ARF) family members involved in fig (Ficus carica L.) fruit development

The auxin response factor (ARF) combines with AuxREs cis-acting elements in response to auxin to regulate plant development. To date, no comprehensive analysis of ARF genes expressed during fruit development has been conducted for common fig (Ficus carica L.). In this study, members of the FcARF gene family were screened, identified in the fig genome database and their features characterized using bioinformatics. Twenty FcARF genes were clustered into three classes, with almost similar highly conserved DBD (B3-like DNA binding domain), AUX/IAA (auxin/indole-3-acetic acid gene family) and MR domain structure among class members. Analysis of amino acid species in MR domain revealed 10 potential transcription activators and 10 transcription inhibitors, and 17 FcARF members were predicted to be located in the nucleus. DNA sequence analysis showed that the ARF gene family consisted of 4-25 exons, and the promoter region contained 16 cis-acting elements involved in stress response, hormone response and flavonoid biosynthesis. ARF genes were expressed in most tissues of fig, especially flower and peel. Transcriptomics analysis results showed that FcARF2, FcARF11 and FcARF12, belonging to class-Ia, were stably and highly expressed in the early development stage of flower and peel of 'Purple peel' fig. However, their expression levels decreased after maturity. Expression of class-Ic member FcARF3 conformed to the regularity of fig fruit development. These four potential transcription inhibitors may regulate fruit growth and development of 'Purple Peel' fig. This study provides comprehensive information on the fig ARF gene family, including gene structure, chromosome position, phylogenetic relationship and expression pattern. Our work provides a foundation for further research on auxin-mediated fig fruit development.

Genome-Wide Analysis of Anthocyanin Biosynthesis Regulatory WD40 Gene FcTTG1 and Related Family in Ficus carica L

WD40 proteins serve as crucial regulators in a broad spectrum of plant developmental and physiological processes, including anthocyanin biosynthesis. However, in fig (Ficus carica L.), neither the WD40 family nor any member involved in anthocyanin biosynthesis has been elucidated. In the present study, 204 WD40 genes were identified from the fig genome and phylogenetically classified into 5 clusters and 12 subfamilies. Bioinformatics analysis prediction localized 109, 69, and 26 FcWD40 proteins to the cytoplasm, nucleus and other cellular compartments, respectively. RNA-seq data mining revealed 127 FcWD40s expressed at FPKM > 10 in fig fruit. Most of these genes demonstrated higher expression in the early stages of fruit development. FcWD40-97 was recruited according to three criteria: high expression in fig fruit, predicted nuclear localization, and closest clustering with TTG1s identified in other plants. FcWD40-97, encoding 339 amino acids including 5 WD-repeat motifs, showed 88.01 and 87.94% amino acid sequence similarity to apple and peach TTG1, respectively. The gene is located on fig chromosome 4, and is composed of 1 intron and 2 exons. Promoter analysis revealed multiple light-responsive elements, one salicylic acid-responsive element, three methyl jasmonate-responsive elements, and one MYB-binding site involved in flavonoid biosynthesis gene regulation. FcWD40-97 was in the FPKM > 100 expression level group in fig fruit, and higher expression was consistently found in the peel compared to the flesh at the same development stages. Expression level did not change significantly under light deprivation, whereas in leaves and roots, its expression was relatively low. Transient expression verified FcWD40-97's localization to the nucleus. Yeast two-hybrid (Y2H) and biomolecular fluorescence complementation (BiFC) assays revealed that FcWD40-97 interacts with FcMYB114, FcMYB123, and FcbHLH42 proteins in vitro and in vivo, showing that FcWD40-97 functions as a member of the MYB-bHLH-WD40 (MBW) complex in anthocyanin-biosynthesis regulation in fig. We therefore renamed FcWD40-97 as FcTTG1. Our results provide the first systematic analysis of the FcWD40 family and identification of FcTTG1 in fig pigmentation.

Proteome Analysis of Vacuoles Isolated from Fig (Ficus carica L.) Flesh during Fruit Development

Fruit flesh cell vacuoles play a pivotal role in fruit growth and quality formation. In the present study, intact vacuoles were carefully released and collected from protoplasts isolated from flesh cells at five sampling times along fig fruit development. Label-free quantification and vacuole proteomic analysis identified 1,251 proteins, 1,137 of which were recruited as differentially abundant proteins (DAPs) by fold change ≥ 1.5, P < 0.05. DAPs were assigned to 10 functional categories; among them, 238, 186, 109, 93 and 90 were annotated as metabolism, transport proteins, membrane fusion or vesicle trafficking, protein fate and stress response proteins, respectively. Decreased numbers of DAPs were uncovered along fruit development. The overall changing pattern of DAPs revealed two major proteome landscape conversions in fig flesh cell vacuoles: the first occurred when fruit developed from late-stage I to mid-stage II, and the second occurred when the fruit started ripening. Metabolic proteins related to glycosidase, lipid and extracellular proteins contributing to carbohydrate storage and vacuole expansion, and protein-degrading proteins determining vacuolar lytic function were revealed. Key tonoplast proteins contributing to vacuole expansion, cell growth and fruit quality formation were also identified. The revealed comprehensive changes in the vacuole proteome during flesh development were compared with our previously published vacuole proteome of grape berry. The information expands our knowledge of the vacuolar proteome and the protein basis of vacuole functional evolution during fruit development and quality formation.

A chromosome-level genome assembly of Artocarpus nanchuanensis (Moraceae), an extremely endangered fruit tree

Artocarpus nanchuanensis (Moraceae), which is naturally distributed in China, is a representative and extremely endangered tree species. In this study, we obtained a high-quality chromosome-scale genome assembly and annotation information for A. nanchuanensis using integrated approaches, including Illumina, Nanopore sequencing platform, and Hi-C. A total of 128.71 Gb of raw Nanopore reads were generated from 20-kb libraries, and 123.38 Gb of clean reads were obtained after filtration with 160.34× coverage depth and a 17.48-kb average read length. The final assembled A. nanchuanensis genome was 769.44 Mb with a 2.09 Mb contig N50, and 99.62% (766.50 Mb) of the assembled data was assigned to 28 pseudochromosomes. In total, 39,596 genes (95.10%, 39,596/41,636) were successfully annotated, and 129 metabolic pathways were detected. Plants disease resistance/insect resistance genes, plant-pathogen interaction metabolic pathways, and abundant biosynthesis pathways of vitamins, flavonoid, and gingerol were detected. Unigene reveals the basis of species-specific functions, and gene family in contraction and expansion generally implies strong functional differences in the evolution. Compared with other related species, a total of 512 unigenes, 309 gene families in contraction, and 559 gene families in expansion were detected in A. nanchuanensis. This A. nanchuanensis genome information provides an important resource to expand our understanding of the unique biological processes, nutritional and medicinal benefits, and evolutionary relationship of this species. The study of gene function and metabolic pathway in A. nanchuanensis may reveal the theoretical basis of a special trait in A. nanchuanensis and promote the study and utilization of its rare medicinal value.

Unraveling coevolutionary dynamics using ecological genomics

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

Metabolic Profiling of Terpene Diversity and the Response of Prenylsynthase-Terpene Synthase Genes during Biotic and Abiotic Stresses in Dendrobium catenatum

Dendrobium catenatum is a widely cultivated Chinese orchid herb rich in abundant secondary metabolites, such as terpenes. However, terpene distribution and characterization of terpene biosynthesis-related genes remain unknown in D. catenatum. In this study, metabolic profiling was performed to analyze terpene distribution in the root, stem, leaf, and flower of D. catenatum. A total of 74 terpene compounds were identified and classified. Clustering analysis revealed that terpene compounds exhibited a tissue-specific accumulation, including monoterpenes in the flowers, sesquiterpenes in the stems, and triterpenes in the roots. Transcriptome analysis revealed that the ‘terpenoid backbone biosynthesis’ pathway was only significantly enriched in root vs. flower. The expression of terpene biosynthesis-related genes was spatiotemporal in the flowers. Prenylsynthase-terpene synthases (PS-TPSs) are the largest and core enzymes for generating terpene diversity. By systematic sequence analysis of six species, 318 PS-TPSs were classified into 10 groups and 51 DcaPS-TPSs were found in eight of them. Eighteen DcaPS-TPSs were regulated by circadian rhythm under drought stress. Most of the DcaPS-TPSs were influenced by cold stress and fungi infection. The cis-element of the majority of the DcaPS-TPS promoters was related to abiotic stress and plant development. Methyl jasmonate levels were significantly associated with DcaTPSs expression and terpene biosynthesis. These results provide insight into further functional investigation of DcaPS-TPSs and the regulation of terpene biosynthesis in Dendrobium.

A chromosome-level genome assembly of the pollinating fig wasp Valisia javana

Fig wasp has always been thought the species-specific pollinator for their host fig (Moraceae, Ficus) and constitute a model system with its host to study co-evolution and co-speciation. The availability of a high-quality genome will help to further reveal the mechanisms underlying these characteristics. Here, we present a high-quality chromosome-level genome for Valisa javana developed by a combination of PacBio long-read and Illumina short-read. The assembled genome size is 296.34 Mb from 13 contigs with a contig N50 length of 26.76 kb. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid fig wasps living in syconium of its highly specific host. Protein-coding genes associated with chemosensory, detoxification and venom genes were identified. Several differentially expressed genes in transcriptome data of V. javana between odor-stimulated samples and the controls have been identified in some olfactory signal transduction pathways, e.g. olfactory transduction, cAMP, cGMP-PKG, Calcim, Ras and Rap1. This study provides a valuable genomic resource for a fig wasp, and sheds insight into further revealing the mechanisms underlying their adaptive traits to their hosts in different places and co-speciation with their host.

SunUp and Sunset genomes revealed impact of particle bombardment mediated transformation and domestication history in papaya

Transgenic papaya is widely publicized for controlling papaya ringspot virus. However, the impact of particle bombardment on the genome remains unknown. The transgenic SunUp and its progenitor Sunset genomes were assembled into 351.5 and 350.3 Mb in nine chromosomes, respectively. We identified a 1.64 Mb insertion containing three transgenic insertions in SunUp chromosome 5, consisting of 52 nuclear-plastid, 21 nuclear-mitochondrial and 1 nuclear genomic fragments. A 591.9 kb fragment in chromosome 5 was translocated into the 1.64 Mb insertion. We assembled a gapless 9.8 Mb hermaphrodite-specific region of the Y^h chromosome and its 6.0 Mb X counterpart. Resequencing 86 genomes revealed three distinct groups, validating their geographic origin and breeding history. We identified 147 selective sweeps and defined the essential role of zeta-carotene desaturase in carotenoid accumulation during domestication. Our findings elucidated the impact of particle bombardment and improved our understanding of sex chromosomes and domestication to expedite papaya improvement.

Profiles of telomeric repeats in Insecta reveal diverse forms of telomeric motifs in Hymenopterans

Telomeres consist of highly conserved simple tandem telomeric repeat motif (TRM): (TTAGG)n in arthropods, (TTAGGG)n in vertebrates, and (TTTAGGG)n in most plants. TRM can be detected from chromosome-level assembly, which typically requires long-read sequencing data. To take advantage of short-read data, we developed an ultra-fast Telomeric Repeats Identification Pipeline and evaluated its performance on 91 species. With proven accuracy, we applied Telomeric Repeats Identification Pipeline in 129 insect species, using 7 Tbp of short-read sequences. We confirmed (TTAGG)n as the TRM in 19 orders, suggesting it is the ancestral form in insects. Systematic profiling in Hymenopterans revealed a diverse range of TRMs, including the canonical 5-bp TTAGG (bees, ants, and basal sawflies), three independent losses of tandem repeat form TRM (Ichneumonoids, hunting wasps, and gall-forming wasps), and most interestingly, a common 8-bp (TTATTGGG)n in Chalcid wasps with two 9-bp variants in the miniature wasp (TTACTTGGG) and fig wasps (TTATTGGGG). Our results identified extraordinary evolutionary fluidity of Hymenopteran TRMs, and rapid evolution of TRM and repeat abundance at all evolutionary scales, providing novel insights into telomere evolution.