榕-传粉小蜂化学通讯机制研究进展

榕属植物与其传粉小蜂组成了高度专一的专性共生关系(榕-蜂共生系统),如此高度紧密的互作关系被认为是驱动两者多样化的关键因素。榕-蜂共生系统主要依靠化学通讯完成相互识别,但目前仍不清楚化学通讯是如何维系现有共生关系并促进物种形成的。结合已有研究,系统梳理了榕-蜂共生系统化学通讯的基础与两者特异性识别的机制,阐述化学通讯在物种和种群层次对维持这一专性传粉关系的重要贡献,进而探讨化学通讯如何在协同成种和宿主转移成种两种模式中介导物种形成。最后,结合生理与多组学等技术展望榕-蜂共生系统的未来研究方向,为深入解析植物与昆虫协同进化的机制以及全球变化下物种的潜在响应模式提供重要参考。     

Expression and evolutionary divergence of the non-conventional olfactory receptor in four species of fig wasp associated with one species of fig

Background  

The interactions of fig wasps and their host figs provide a model for investigating co-evolution. Fig wasps have specialized morphological characters and lifestyles thought to be adaptations to living in the fig's syconium. Although these aspects of natural history are well documented, the genetic mechanism(s) underlying these changes remain(s) unknown. Fig wasp olfaction is the key to host-specificity. The Or83b gene class, an unusual member of olfactory receptor family, plays a critical role in enabling the function of conventional olfactory receptors. Four Or83b orthologous genes from one pollinator (PFW) (Ceratosolen solmsi) and three non-pollinator fig wasps (NPFWs) (Apocrypta bakeri, Philotrypesis pilosa and Philotrypesis sp.) associated with one species of fig (Ficus hispida) can be used to better understand the molecular mechanism underlying the fig wasp's adaptation to its host. We made a comparison of spatial tissue-specific expression patterns and substitution rates of one orthologous gene in these fig wasps and sought evidence for selection pressures.     

High-quality,chromosome-scale genome assemblies: comparisons of three Diaphorina citri (Asian citrus psyllid) geographic populations

The Asian citrus psyllid, Diaphorina citri, is the insect vector of the causal agent of huanglongbing (HLB), a devastating bacterial disease of commercial citrus. Presently, few genomic resources exist for D. citri. In this study, we utilized PacBio HiFi and chromatin confirmation contact (Hi-C) sequencing to sequence, assemble, and compare three high-quality, chromosome-scale genome assemblies of D. citri collected from California, Taiwan, and Uruguay. Our assemblies had final sizes of 282.67 Mb (California), 282.89 Mb (Taiwan), and 266.67 Mb (Uruguay) assembled into 13 pseudomolecules—a reduction in assembly size of 41–45% compared with previous assemblies which we validated using flow cytometry. We identified the X chromosome in D. citri and annotated each assembly for repetitive elements, protein-coding genes, transfer RNAs, ribosomal RNAs, piwi-interacting RNA clusters, and endogenous viral elements. Between 19,083 and 20,357 protein-coding genes were predicted. Repetitive DNA accounts for 36.87–38.26% of each assembly. Comparative analyses and mitochondrial haplotype networks suggest that Taiwan and Uruguay D. citri are more closely related, while California D. citri are closely related to Florida D. citri. These high-quality, chromosome-scale assemblies provide new genomic resources to researchers to further D. citri and HLB research.     

Dynamics and Impacts of Transposable Element Proliferation in the Drosophila nasuta Species Group Radiation

Transposable element (TE) mobilization is a constant threat to genome integrity. Eukaryotic organisms have evolved robust defensive mechanisms to suppress their activity, yet TEs can escape suppression and proliferate, creating strong selective pressure for host defense to adapt. This genomic conflict fuels a never-ending arms race that drives the rapid evolution of TEs and recurrent positive selection of genes involved in host defense; the latter has been shown to contribute to postzygotic hybrid incompatibility. However, how TE proliferation impacts genome and regulatory divergence remains poorly understood. Here, we report the highly complete and contiguous (N50 = 33.8–38.0 Mb) genome assemblies of seven closely related Drosophila species that belong to the nasuta species group—a poorly studied group of flies that radiated in the last 2 My. We constructed a high-quality de novo TE library and gathered germline RNA-seq data, which allowed us to comprehensively annotate and compare TE insertion patterns between the species, and infer the evolutionary forces controlling their spread. We find a strong negative association between TE insertion frequency and expression of genes nearby; this likely reflects survivor bias from reduced fitness impact of TEs inserting near lowly expressed, nonessential genes, with limited TE-induced epigenetic silencing. Phylogenetic analyses of insertions of 147 TE families reveal that 53% of them show recent amplification in at least one species. The most highly amplified TE is a nonautonomous DNA element (Drosophila INterspersed Element; DINE) which has gone through multiple bouts of expansions with thousands of full-length copies littered throughout each genome. Across all TEs, we find that TEs expansions are significantly associated with high expression in the expanded species consistent with suppression escape. Thus, whereas horizontal transfer followed by the invasion of a naïve genome has been highlighted to explain the long-term survival of TEs, our analysis suggests that evasion of host suppression of resident TEs is a major strategy to persist over evolutionary times. Altogether, our results shed light on the heterogenous and context-dependent nature in which TEs affect gene regulation and the dynamics of rampant TE proliferation amidst a recently radiated species group.     

When parasitic wasps hijacked viruses: genomic and functional evolution of polydnaviruses

The Polydnaviridae (PDV), including the Bracovirus (BV) and Ichnovirus genera, originated from the integration of unrelated viruses in the genomes of two parasitoid wasp lineages, in a remarkable example of convergent evolution. Functionally active PDVs represent the most compelling evolutionary success among endogenous viral elements (EVEs). BV evolved from the domestication by braconid wasps of a nudivirus 100 Ma. The nudivirus genome has become an EVE involved in BV particle production but is not encapsidated. Instead, BV genomes have co-opted virulence genes, used by the wasps to control the immunity and development of their hosts. Gene transfers and duplications have shaped BV genomes, now encoding hundreds of genes. Phylogenomic studies suggest that BVs contribute largely to wasp diversification and adaptation to their hosts. A genome evolution model explains how multidirectional wasp adaptation to different host species could have fostered PDV genome extension. Integrative studies linking ecological data on the wasp to genomic analyses should provide new insights into the adaptive role of particular BV genes. Forthcoming genomic advances should also indicate if the associations between endoparasitoid wasps and symbiotic viruses evolved because of their particularly intimate interactions with their hosts, or if similar domesticated EVEs could be uncovered in other parasites.     

Chromosome-level genome assembly and characterization of Sophora Japonica

Sophora japonica is a medium-size deciduous tree belonging to Leguminosae family and famous for its high ecological, economic and medicinal value. Here, we reveal a draft genome of S. japonica, which was ∼511.49 Mb long (contig N50 size of 17.34 Mb) based on Illumina, Nanopore and Hi-C data. We reliably assembled 110 contigs into 14 chromosomes, representing 91.62% of the total genome, with an improved N50 size of 31.32 Mb based on Hi-C data. Further investigation identified 271.76 Mb (53.13%) of repetitive sequences and 31,000 protein-coding genes, of which 30,721 (99.1%) were functionally annotated. Phylogenetic analysis indicates that S. japonica separated from Arabidopsis thaliana and Glycine max ∼107.53 and 61.24 million years ago, respectively. We detected evidence of species-specific and common-legume whole-genome duplication events in S. japonica. We further found that multiple TF families (e.g. BBX and PAL) have expanded in S. japonica, which might have led to its enhanced tolerance to abiotic stress. In addition, S. japonica harbours more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed no obvious differentiation among geographical groups and the effective population size continuously declined since 2 Ma. Our genomic data provide a powerful comparative framework to study the adaptation, evolution and active ingredients biosynthesis in S. japonica. More importantly, our high-quality S. japonica genome is important for elucidating the biosynthesis of its main bioactive components, and improving its production and/or processing.     

Inferring the Genetic Basis of Sex Determination from the Genome of a Dioecious Nightshade

Dissecting the genetic mechanisms underlying dioecy (i.e., separate female and male individuals) is critical for understanding the evolution of this pervasive reproductive strategy. Nonetheless, the genetic basis of sex determination remains unclear in many cases, especially in systems where dioecy has arisen recently. Within the economically important plant genus Solanum (∼2,000 species), dioecy is thought to have evolved independently at least 4 times across roughly 20 species. Here, we generate the first genome sequence of a dioecious Solanum and use it to ascertain the genetic basis of sex determination in this species. We de novo assembled and annotated the genome of Solanum appendiculatum (assembly size: ∼750 Mb scaffold N50: 0.92 Mb; ∼35,000 genes), identified sex-specific sequences and their locations in the genome, and inferred that males in this species are the heterogametic sex. We also analyzed gene expression patterns in floral tissues of males and females, finding approximately 100 genes that are differentially expressed between the sexes. These analyses, together with observed patterns of gene-family evolution specific to S. appendiculatum, consistently implicate a suite of genes from the regulatory network controlling pectin degradation and modification in the expression of sex. Furthermore, the genome of a species with a relatively young sex-determination system provides the foundational resources for future studies on the independent evolution of dioecy in this clade.     

A spinach genome assembly with remarkable completeness,and its use for rapid identification of candidate genes for agronomic traits

Spinach (Spinacia oleracea) is grown as a nutritious leafy vegetable worldwide. To accelerate spinach breeding efficiency, a high-quality reference genome sequence with great completeness and continuity is needed as a basic infrastructure. Here, we used long-read and linked-read technologies to construct a de novo spinach genome assembly, designated SOL_r1.1, which was comprised of 287 scaffolds (total size: 935.7 Mb; N₅₀ = 11.3 Mb) with a low proportion of undetermined nucleotides (Ns = 0.34%) and with high gene completeness (BUSCO complete 96.9%). A genome-wide survey of resistance gene analogues identified 695 genes encoding nucleotide-binding site domains, receptor-like protein kinases, receptor-like proteins and transmembrane-coiled coil domains. Based on a high-density double-digest restriction-site associated DNA sequencing-based linkage map, the genome assembly was anchored to six pseudomolecules representing ∼73.5% of the whole genome assembly. In addition, we used SOL_r1.1 to identify quantitative trait loci for bolting timing and fruit/seed shape, which harbour biologically plausible candidate genes, such as homologues of the FLOWERING LOCUS T and EPIDERMAL PATTERNING FACTOR-LIKE genes. The new genome assembly, SOL_r1.1, will serve as a useful resource for identifying loci associated with important agronomic traits and for developing molecular markers for spinach breeding/selection programs.     

Genome assembly of Scorias spongiosa and comparative genomics provide insights into ecological adaptation of honeydew-dependent sooty mould fungi

Sooty moulds are fungi of economic importance and with unique lifestyle mainly growing on insect honeydew. However, the lack of genomic data hinders investigation of genetic mechanisms underlying their ecological adaptation. With long-read sequencing technology, we generated the genome of Scorias spongiosa, an extraordinary sooty mould fungus associated with honeydew of colony aphids and producing large fruiting bodies. A 24.21 Mb high-quality genome assembly with a N50 length of 3.37 Mb was obtained. The genome contained 7758 protein coding genes, 97.13% of which were homologous to known genes, and approximately 0.29 Mb repeat sequences. Comparative genomics showed S. spongiosa lost relatively more gene families and contained fewer species-specific genes and gene families, with many CAZyme families and sugar transporters reduced or absent. This study not only promotes understanding of the ecological adaptation of sooty moulds, but also provides valuable genomic data resource for future comparative genomic and genetic studies.     

Divergence of a genomic island leads to the evolution of melanization in a halophyte root fungus

Understanding how organisms adapt to extreme living conditions is central to evolutionary biology. Dark septate endophytes (DSEs) constitute an important component of the root mycobiome and they are often able to alleviate host abiotic stresses. Here, we investigated the molecular mechanisms underlying the beneficial association between the DSE Laburnicola rhizohalophila and its host, the native halophyte Suaeda salsa, using population genomics. Based on genome-wide Fst (pairwise fixation index) and Vst analyses, which compared the variance in allele frequencies of single-nucleotide polymorphisms (SNPs) and copy number variants (CNVs), respectively, we found a high level of genetic differentiation between two populations. CNV patterns revealed population-specific expansions and contractions. Interestingly, we identified a ~20 kbp genomic island of high divergence with a strong sign of positive selection. This region contains a melanin-biosynthetic polyketide synthase gene cluster linked to six additional genes likely involved in biosynthesis, membrane trafficking, regulation, and localization of melanin. Differences in growth yield and melanin biosynthesis between the two populations grown under 2% NaCl stress suggested that this genomic island contributes to the observed differences in melanin accumulation. Our findings provide a better understanding of the genetic and evolutionary mechanisms underlying the adaptation to saline conditions of the L. rhizohalophila–S. salsa symbiosis.Subject terms: Population dynamics, Fungal ecology, Population genetics     

Chromosome-level genome of Tibetan naked carp (Gymnocypris przewalskii) provides insights into Tibetan highland adaptation

Gymnocypris przewalskii, a cyprinid fish endemic to the Qinghai-Tibetan Plateau, has evolved unique morphological, physiological and genetic characteristics to adapt to the highland environment. Herein, we assembled a high-quality G. przewalskii tetraploid genome with a size of 2.03 Gb and scaffold N50 of 44.93 Mb, which was anchored onto 46 chromosomes. The comparative analysis suggested that gene families related to highland adaptation were significantly expanded in G. przewalskii. According to the G. przewalskii genome, we evaluated the phylogenetic relationship of 13 schizothoracine fishes, and inferred that the demographic history of G. przewalskii was strongly associated with geographic and eco-environmental alterations. We noticed that G. przewalskii experienced whole-genome duplication, and genes preserved post duplication were functionally associated with adaptation to high salinity and alkalinity. In conclusion, a chromosome-scale G. przewalskii genome provides an important genomic resource for teleost fish, and will particularly promote our understanding of the molecular evolution and speciation of fish in the highland environment.     

Population Genomics Reveals Incipient Speciation,Introgression, and Adaptation in the African Mona Monkey (Cercopithecus mona)

Guenons (tribe Cercopithecini) are the most widely distributed nonhuman primate in the tropical forest belt of Africa and show considerable phenotypic, taxonomic, and ecological diversity. However, genomic information for most species within this group is still lacking. Here, we present a high-quality de novo genome (total 2.90 Gb, contig N50 equal to 22.7 Mb) of the mona monkey (Cercopithecus mona), together with genome resequencing data of 13 individuals sampled across Nigeria. Our results showed differentiation between populations from East and West of the Niger River ∼84 ka and potential ancient introgression in the East population from other mona group species. The PTPRK, FRAS1, BNC2, and EDN3 genes related to pigmentation displayed signals of introgression in the East population. Genomic scans suggest that immunity genes such as AKT3 and IL13 (possibly involved in simian immunodeficiency virus defense), and G6PD, a gene involved in malaria resistance, are under positive natural selection. Our study gives insights into differentiation, natural selection, and introgression in guenons.     

Draft Genome Assembly and Annotation for Cutaneotrichosporon dermatis NICC30027, an Oleaginous Yeast Capable of Simultaneous Glucose and Xylose Assimilation

The identification of oleaginous yeast species capable of simultaneously utilizing xylose and glucose as substrates to generate value-added biological products is an area of key economic interest. We have previously demonstrated that the Cutaneotrichosporon dermatis NICC30027 yeast strain is capable of simultaneously assimilating both xylose and glucose, resulting in considerable lipid accumulation. However, as no high-quality genome sequencing data or associated annotations for this strain are available at present, it remains challenging to study the metabolic mechanisms underlying this phenotype. Herein, we report a 39,305,439 bp draft genome assembly for C. dermatis NICC30027 comprised of 37 scaffolds, with 60.15% GC content. Within this genome, we identified 524 tRNAs, 142 sRNAs, 53 miRNAs, 28 snRNAs, and eight rRNA clusters. Moreover, repeat sequences totaling 1,032,129 bp in length were identified (2.63% of the genome), as were 14,238 unigenes that were 1,789.35 bp in length on average (64.82% of the genome). The NCBI non-redundant protein sequences (NR) database was employed to successfully annotate 11,795 of these unigenes, while 3,621 and 11,902 were annotated with the Swiss-Prot and TrEMBL databases, respectively. Unigenes were additionally subjected to pathway enrichment analyses using the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster of Orthologous Groups of proteins (COG), Clusters of orthologous groups for eukaryotic complete genomes (KOG), and Non-supervised Orthologous Groups (eggNOG) databases. Together, these results provide a foundation for future studies aimed at clarifying the mechanistic basis for the ability of C. dermatis NICC30027 to simultaneously utilize glucose and xylose to synthesize lipids.     

Phylogeography and virulence structure of the powdery mildew population on its 'new' host triticale

Background

The interaction between insects and plants takes myriad forms in the generation of spectacular diversity. In this association a species host range is fundamental and often measured using an estimate of phylogenetic concordance between species. Pollinating fig wasps display extreme host species specificity, but the intraspecific variation in empirical accounts of host affiliation has previously been underestimated. In this investigation, lineage delimitation and codiversification tests are used to generate and discuss hypotheses elucidating on pollinating fig wasp associations with Ficus.

Results

Statistical parsimony and AMOVA revealed deep divergences at the COI locus within several pollinating fig wasp species that persist on the same host Ficus species. Changes in branching patterns estimated using the generalized mixed Yule coalescent test indicated lineage duplication on the same Ficus species. Conversely, Elisabethiella and Alfonsiella fig wasp species are able to reproduce on multiple, but closely related host fig species. Tree reconciliation tests indicate significant codiversification as well as significant incongruence between fig wasp and Ficus phylogenies.

Conclusions

The findings demonstrate more relaxed pollinating fig wasp host specificity than previously appreciated. Evolutionarily conservative host associations have been tempered by horizontal transfer and lineage duplication among closely related Ficus species. Independent and asynchronistic diversification of pollinating fig wasps is best explained by a combination of both sympatric and allopatric models of speciation. Pollinator host preference constraints permit reproduction on closely related Ficus species, but uncertainty of the frequency and duration of these associations requires better resolution.     

Bracoviruses recruit host integrases for their integration into caterpillar’s genome

Some DNA viruses infect host animals usually by integrating their DNAs into the host genome. However, the mechanisms for integration remain largely unknown. Here, we find that Cotesia vestalis bracovirus (CvBV), a polydnavirus of the parasitic wasp C. vestalis (Haliday), integrates its DNA circles into host Plutella xylostella (L.) genome by two distinct strategies, conservatively and randomly, through high-throughput sequencing analysis. We confirmed that the conservatively integrating circles contain an essential “8+5” nucleotides motif which is required for integration. Then we find CvBV circles are integrated into the caterpillar’s genome in three temporal patterns, the early, mid and late stage-integration. We further identify that three CvBV-encoded integrases are responsible for some, but not all of the virus circle integrations, indeed they mainly participate in the processes of early stage-integration. Strikingly, we find two P. xylostella retroviral integrases (PxIN1 and PxIN2) are highly induced upon wasp parasitism, and PxIN1 is crucial for integration of some other early-integrated CvBV circles, such as CvBV_04, CvBV_12 and CvBV_24, while PxIN2 is important for integration of a late-integrated CvBV circle, CvBV_21. Our data uncover a novel mechanism in which CvBV integrates into the infected host genome, not only by utilizing its own integrases, but also by recruiting host enzymes. These findings will strongly deepen our understanding of how bracoviruses regulate and integrate into their hosts.