The nuclear and mitochondrial genome assemblies of Tetragonisca angustula (Apidae: Meliponini), a tiny yet remarkable pollinator in the Neotropics.

BACKGROUND: The field of bee genomics has considerably advanced in recent years, however, the most diverse group of honey producers on the planet, the stingless bees, are still largely neglected. In fact, only eleven of the ~ 600 described stingless bee species have been sequenced, and only three using a long-read (LR) sequencing technology. Here, we sequenced the nuclear and mitochondrial genomes of the most common, widespread and broadly reared stingless bee in Brazil and other neotropical countries-Tetragonisca angustula (popularly known in Brazil as jataí). RESULTS: A total of 48.01 Gb of DNA data were generated, including 2.31 Gb of Pacific Bioscience HiFi reads and 45.70 Gb of Illumina short reads (SRs). Our preferred assembly comprised 683 contigs encompassing 284.49 Mb, 62.84 Mb of which (22.09%) corresponded to 445,793 repetitive elements. N50, L50 and complete BUSCOs reached 1.02 Mb, 91 contigs and 97.1%, respectively. We predicted that the genome of T. angustula comprises 17,459 protein-coding genes and 4,108 non-coding RNAs. The mitogenome consisted of 17,410 bp, and all 37 genes were found to be on the positive strand, an unusual feature among bees. A phylogenomic analysis of 26 hymenopteran species revealed that six odorant receptor orthogroups of T. angustula were found to be experiencing rapid evolution, four of them undergoing significant contractions. CONCLUSIONS: Here, we provided the first nuclear and mitochondrial genome assemblies for the ecologically and economically important T. angustula, the fourth stingless bee species to be sequenced with LR technology thus far. We demonstrated that even relatively small amounts of LR data in combination with sufficient SR data can yield high-quality genome assemblies for bees.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process.

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

Functional and phylogenetic relationships link predators to plant diversity via trophic and non-trophic pathways.

Human-induced biodiversity loss negatively affects ecosystem function, but the interactive effects of biodiversity change across trophic levels remain insufficiently understood. We sampled arboreal spiders and lepidopteran larvae across seasons in 2 years in a subtropical tree diversity experiment, and then disentangled the links between tree diversity and arthropod predator diversity by deconstructing the pathways among multiple components of diversity (taxonomic, phylogenetic and functional) with structural equation models. We found that herbivores were major mediators of plant species richness effects on abundance, species richness, functional and phylogenetic diversity of predators, while phylogenetic, functional and structural diversity of trees were also important mediators of this process. However, the strength and direction differed between functional, structural and phylogenetic diversity effects, indicating different underlying mechanisms for predator community assembly. Abundance and multiple diversity components of predators were consistently affected by tree functional diversity, indicating that the variation in structure and environment caused by plant functional composition might play key roles in predator community assembly. Our study highlights the importance of an integrated approach based on multiple biodiversity components in understanding the consequences of biodiversity loss in multitrophic communities.

Tree dissimilarity determines multi-dimensional beta-diversity of herbivores and carnivores via bottom-up effects.

Global biodiversity decline and its cascading effects through trophic interactions pose a severe threat to human society. Establishing the impacts of biodiversity decline requires a more thorough understanding of multi-trophic interactions and, more specifically, the effects that loss of diversity in primary producers has on multi-trophic community assembly. Within a synthetic conceptual framework for multi-trophic beta-diversity, we tested a series of hypotheses on neutral and niche-based bottom-up processes in assembling herbivore and carnivore communities in a subtropical forest using linear models, hieratical variance partitioning based on linear mixed-effects models (LMMs) and simulation. We found that the observed taxonomic, phylogenetic and functional beta-diversity of both herbivorous caterpillars and carnivorous spiders were significantly and positively related to tree dissimilarity. Linear models and variance partitioning for LMMs jointly suggested that as a result of bottom-up effects, producer dissimilarities were predominant in structuring consumer dissimilarity, the strength of which highly depended on the trophic dependencies on producers, the diversity facet examined, and data quality. Importantly, linear models for standardized beta-diversities against producer dissimilarities implied a transition between niche-based processes such as environmental filtering and competitive exclusion, which supports the role of bottom-up effect in determining consumer community assembly. These findings enrich our mechanistic understanding of the 'Diversity Begets Diversity' hypothesis and the complexity of higher-trophic community assembly, which is fundamental for sustainable biodiversity conservation and ecosystem management.

Tree communities and functional traits determine herbivore compositional turnover.

There are many factors known to drive species turnover, although the mechanisms by which these operate are less clear. Based on comprehensive datasets from the largest tree diversity experiment worldwide (BEF-China), we used shared herbivore species (zeta diversity) and multi-site generalized dissimilarity modelling to investigate the patterns and determinants of species turnover of Lepidoptera herbivores among study plots across a gradient in tree species richness. We found that zeta diversity declined sharply with an increasing number of study plots, with complete changes in caterpillar species composition observed even at the fine spatial scale of our study. Plant community characteristics rather than abiotic factors were found to play key roles in driving caterpillar compositional turnover, although these effects varied with an increasing number of study plots considered, due to the varying contributions of rare and common species to compositional turnover. Our study reveals details of the impact of phylogeny- and trait-mediated processes of trees on herbivore compositional turnover, which has implications for forest management and conservation and shows potential avenues for maintenance of heterogeneity in herbivore communities.

A Simulation-Based Evaluation of Tip-Dating Under the Fossilized Birth-Death Process.

Bayesian molecular dating is widely used to study evolutionary timescales. This procedure usually involves phylogenetic analysis of nucleotide sequence data, with fossil-based calibrations applied as age constraints on internal nodes of the tree. An alternative approach is tip-dating, which explicitly includes fossil data in the analysis. This can be done, for example, through the joint analysis of molecular data from present-day taxa and morphological data from both extant and fossil taxa. In the context of tip-dating, an important development has been the fossilized birth-death process, which allows non-contemporaneous tips and sampled ancestors while providing a model of lineage diversification for the prior on the tree topology and internal node times. However, tip-dating with fossils faces a number of considerable challenges, especially, those associated with fossil sampling and evolutionary models for morphological characters. We conducted a simulation study to evaluate the performance of tip-dating using the fossilized birth-death model. We simulated fossil occurrences and the evolution of nucleotide sequences and morphological characters under a wide range of conditions. Our analyses of these data show that the number and the maximum age of fossil occurrences have a greater influence than the degree of among-lineage rate variation or the number of morphological characters on estimates of node times and the tree topology. Tip-dating with the fossilized birth-death model generally performs well in recovering the relationships among extant taxa but has difficulties in correctly placing fossil taxa in the tree and identifying the number of sampled ancestors. The method yields accurate estimates of the ages of the root and crown group, although the precision of these estimates varies with the probability of fossil occurrence. The exclusion of morphological characters results in a slight overestimation of node times, whereas the exclusion of nucleotide sequences has a negative impact on inference of the tree topology. Our results provide an overview of the performance of tip-dating using the fossilized birth-death model, which will inform further development of the method and its application to key questions in evolutionary biology.

Multi-trophic communities re-establish with canopy cover and microclimate in a subtropical forest biodiversity experiment.

Plant diversity affects multi-trophic communities, but in young regrowth forests, where forest insects are in the process of re-establishment, other biotic and also abiotic factors might be more important. We studied cavity-nesting bees, wasps and their natural enemies along an experimental tree diversity gradient in subtropical South-East China. We compared insect communities of experimental young forests with communities of established natural forests nearby the experiment and tested for direct and indirect effects of tree diversity, tree basal area (a proxy of tree biomass), canopy cover and microclimate on bee and wasp community composition, abundance and species richness. Finally, we tested if the trophic levels of bees, herbivore-hunting wasps, spider-hunting wasps and their natural enemies respond similarly. Forest bee and wasp community composition re-established towards communities of the natural forest with increasing tree biomass and canopy cover. These factors directly and indirectly, via microclimatic conditions, increased the abundance of bees, wasps and their natural enemies. While bee and wasp species richness increased with abundance and both were not related to tree diversity, abundance increased directly with canopy cover, mediated by tree biomass. Abundance of natural enemies increased with host (bee and wasp) abundance irrespective of their trophic position. In conclusion, although maximizing tree diversity is an important goal of reforestation and forest conservation, rapid closure of canopies is also important for re-establishing communities of forest bees, wasps and their natural enemies.

Host functional and phylogenetic composition rather than host diversity structure plant-herbivore networks.

Declining plant diversity alters ecological networks, such as plant-herbivore interactions. However, our knowledge of the potential mechanisms underlying effects of plant species loss on plant-herbivore network structure is still limited. We used DNA barcoding to identify herbivore-host plant associations along declining levels of tree diversity in a large-scale, subtropical biodiversity experiment. We tested for effects of tree species richness, host functional and phylogenetic diversity, and host functional (leaf trait) and phylogenetic composition on species, phylogenetic and network composition of herbivore communities. We found that phylogenetic host composition and related palatability/defence traits but not tree species richness significantly affected herbivore communities and interaction network complexity at both the species and community levels. Our study indicates that evolutionary dependencies and functional traits of host plants determine the composition of higher trophic levels and corresponding interaction networks in species-rich ecosystems. Our findings highlight that characteristics of the species lost have effects on ecosystem structure and functioning across trophic levels that cannot be predicted from mere reductions in species richness.

Comparison of Methods for Molecular Species Delimitation Across a Range of Speciation Scenarios.

Species are fundamental units in biological research and can be defined on the basis of various operational criteria. There has been growing use of molecular approaches for species delimitation. Among the most widely used methods, the generalized mixed Yule-coalescent (GMYC) and Poisson tree processes (PTP) were designed for the analysis of single-locus data but are often applied to concatenations of multilocus data. In contrast, the Bayesian multispecies coalescent approach in the software Bayesian Phylogenetics and Phylogeography (BPP) explicitly models the evolution of multilocus data. In this study, we compare the performance of GMYC, PTP, and BPP using synthetic data generated by simulation under various speciation scenarios. We show that in the absence of gene flow, the main factor influencing the performance of these methods is the ratio of population size to divergence time, while number of loci and sample size per species have smaller effects. Given appropriate priors and correct guide trees, BPP shows lower rates of species overestimation and underestimation, and is generally robust to various potential confounding factors except high levels of gene flow. The single-threshold GMYC and the best strategy that we identified in PTP generally perform well for scenarios involving more than a single putative species when gene flow is absent, but PTP outperforms GMYC when fewer species are involved. Both methods are more sensitive than BPP to the effects of gene flow and potential confounding factors. Case studies of bears and bees further validate some of the findings from our simulation study, and reveal the importance of using an informed starting point for molecular species delimitation. Our results highlight the key factors affecting the performance of molecular species delimitation, with potential benefits for using these methods within an integrative taxonomic framework.

Tree phylogenetic diversity promotes host-parasitoid interactions.

Evidence from grassland experiments suggests that a plant community's phylogenetic diversity (PD) is a strong predictor of ecosystem processes, even stronger than species richness per se This has, however, never been extended to species-rich forests and host-parasitoid interactions. We used cavity-nesting Hymenoptera and their parasitoids collected in a subtropical forest as a model system to test whether hosts, parasitoids, and their interactions are influenced by tree PD and a comprehensive set of environmental variables, including tree species richness. Parasitism rate and parasitoid abundance were positively correlated with tree PD. All variables describing parasitoids decreased with elevation, and were, except parasitism rate, dependent on host abundance. Quantitative descriptors of host-parasitoid networks were independent of the environment. Our study indicates that host-parasitoid interactions in species-rich forests are related to the PD of the tree community, which influences parasitism rates through parasitoid abundance. We show that effects of tree community PD are much stronger than effects of tree species richness, can cascade to high trophic levels, and promote trophic interactions. As during habitat modification phylogenetic information is usually lost non-randomly, even species-rich habitats may not be able to continuously provide the ecosystem process parasitism if the evolutionarily most distinct plant lineages vanish.

A protocol for species delineation of public DNA databases, applied to the Insecta.

Public DNA databases are composed of data from many different taxa, although the taxonomic annotation on sequences is not always complete, which impedes the utilization of mined data for species-level applications. There is much ongoing work on species identification and delineation based on the molecular data itself, although applying species clustering to whole databases requires consolidation of results from numerous undefined gene regions, and introduces significant obstacles in data organization and computational load. In the current paper, we demonstrate an approach for species delineation of a sequence database. All DNA sequences for the insects were obtained and processed. After filtration of duplicated data, delineation of the database into species or molecular operational taxonomic units (MOTUs) followed a three-step process in which (i) the genetic loci L are partitioned, (ii) the species S are delineated within each locus, then (iii) species units are matched across loci to form the matrix L × S, a set of global (multilocus) species units. Partitioning the database into a set of homologous gene fragments was achieved by Markov clustering using edge weights calculated from the amount of overlap between pairs of sequences, then delineation of species units and assignment of species names were performed for the set of genes necessary to capture most of the species diversity. The complexity of computing pairwise similarities for species clustering was substantial at the cytochrome oxidase subunit I locus in particular, but made feasible through the development of software that performs pairwise alignments within the taxonomic framework, while accounting for the different ranks at which sequences are labeled with taxonomic information. Over 24 different homologs, the unidentified sequences numbered approximately 194,000, containing 41,525 species IDs (98.7% of all found in the insect database), and were grouped into 59,173 single-locus MOTUs by hierarchical clustering under parameters optimized independently for each locus. Species units from different loci were matched using a multipartite matching algorithm to form multilocus species units with minimal incongruence between loci. After matching, the insect database as represented by these 24 loci was found to be composed of 78,091 species units in total. 38,574 of these units contained only species labeled data, 34,891 contained only unlabeled data, leaving 4,626 units composed both of labeled and unlabeled sequences. In addition to giving estimates of species diversity of sequence repositories, the protocol developed here will facilitate species-level applications of modern-day sequence data sets. In particular, the L × S matrix represents a post-taxonomic framework that can be used for species-level organization of metagenomic data, and incorporation of these methods into phylogenetic pipelines will yield matrices more representative of species diversity.

Molecular phylogeny reveals independent origins of body scales in Entomobryidae (Hexapoda: Collembola).

Entomobryidae is the largest family in Collembola but relationships within the family have never been subjected to rigorous phylogenetic analyses. Within the family, body scales are present in many species, and are fundamental in the classification at the subfamilial and tribal levels. A molecular phylogeny was reconstructed using the nuclear 18SrRNA and partial 28SrRNA and the mitochondrial 16SrRNA to examine the evolution of scales across Entomobryidae subfamilies. These datasets were analyzed separately and combined, with parsimony, likelihood and Bayesian algorithms. Monophyly of Orchesellinae was not recovered, and it was split into a scaled clade and an unscaled clade, contradicting to all recent taxonomic conceptions. The monophyly of Entomobryinae, Seirinae and Lepidocyrtinae is well supported however within Entomobryinae, the polyphyly of Entomobryini and Willowsiini implies that classification using the presence/absence of scales is not valid. Analyses of ancestral character state reconstruction in Entomobryidae indicate that the presence of body scales have evolved independently at least five times, with a loss of scales occurring independently at least twice. A revision of the family Entomobryidae on molecular and morphological basis is clearly needed.

Bees of the Colletes flavicornis-group from China with description of one new species (Hymenoptera: Apoidea: Colletidae).

Two species of the Colletes flavicornis-group from China are treated in this paper. C. vestitus sp. n. from Xinjiang is illustrated and described, and C. popovi Noskiewicz, 1936 is illustrated and redescribed. Both sexes of the two species are in addition characterized by DNA barcodes. The type specimens of the new species are deposited in the Insect Collection of Institute of Zoology, Chinese Academy of Sciences, Beijing, China.

A chromosome-level genome assembly of the forestry pest Coronaproctus castanopsis.

As an important forestry pest, Coronaproctus castanopsis (Monophlebidae) has caused serious damage to the globally valuable Gutianshan ecosystem, China. In this study, we assembled the first chromosome-level genome of the female specimen of C. castanopsis by merging BGI reads, HiFi long reads and Hi-C data. The assembled genome size is 700.81 Mb, with a scaffold N50 size of 273.84 Mb and a contig N50 size of 12.37 Mb. Hi-C scaffolding assigned 98.32% (689.03 Mb) of C. Castanopsis genome to three chromosomes. The BUSCO analysis (n = 1,367) showed a completeness of 91.2%, comprising 89.2% of single-copy BUSCOs and 2.0% of multicopy BUSCOs. The mapping ratio of BGI, second-generation RNA, third-generation RNA and HiFi reads are 97.84%, 96.15%, 97.96%, and 99.33%, respectively. We also identified 64.97% (455.3 Mb) repetitive elements, 1,373 non-coding RNAs and 10,542 protein-coding genes. This study assembled a high-quality genome of C. castanopsis, which accumulated valuable molecular data for scale insects.

A review of the Colletes succinctus-group (Hymenoptera: Colletidae) from China with redescription of the male of C. gigas.

Four species of the Colletes succinctus-group known from China are reviewed and a key to species is provided. Colletes gigas Cockerell, 1918 is recognized as a member of the Colletes succinctus-group for the first time based on a redescription of the male.

Bees of the Colletes clypearis-group (Hymenoptera: Apoidea: Colletidae) from China with descriptions of seven new species.

Twenty-three known species of the Colletes clypearis speices group from China are treated in this paper. C. cinerascens Morawitz 1893, C. clypearis Morawitz 1876, C. floralis Eversmann 1852, C. impunctatus Nylander 1852, Colletes paratibeticus Kuhlmann 2002 and Colletes sodalis (Cameron 1897) are newly recorded from China. C. harrerioides sp. n., C. heilongtenensis sp. n., C. hirsutus sp. n., C. inspersus sp. n., C. xizangensis sp. n., C. xuezhongi sp. n. and C. yanruae sp. n. are described and illustrated as new species. Checklist of the known species from China in Colletes clypearis-group with distribution, floral records, and an illustrated key to all known males and females from China are provided. The type specimens of the new species are deposited in the Insect Collection of Institute of Zoology, Chinese Academy of Sciences, Beijing, China. 

First record of the bee genus Homalictus Cockerell for China with description of a new species (Hymenoptera: Halictidae: Halictini).

This paper reports the first record of the genus Homalictus from China. We describe and illustrate H. (H.) nabanensis sp. n. collected from the Naban River Watershed National Nature Reserve, Xishuangbanna, Yunnan, China. The type specimens are deposited in Institute of Zoology, Chinese Academy of Sciences, Beijing, China.

New species and records of Trypetinae (Diptera: Tephritidae) from China.

Eight new species of Trypetinae (Diptera: Tephritidae) from China are described and illustrated. The new species are Cornutrypeta chishuiensis Chen sp. nov., C. motuonia Chen sp. nov., Hoplandromyia hubeiensis Chen sp. nov., H. motuonica Chen sp. nov., Magnimyiolia disrupta Chen sp. nov., M tengchongnica Chen sp. nov., M tibetana Chen et Wang sp. nov. and M. yunnanica Chen et Wang sp. nov. Keys to the known species of the genera Cornutrypeta Han et Wang and Hoplandromyia Bezzi from the World, and the genera Magnimyiolia Shiraki from East Asia are updated and provided. C. spinifrons (Schroeder) is newly recorded from China, and the female of C. hunanica Chen et Wang is recorded, described and illustrated for the first time. Illustrations of the type of M huanana Wang and a new distribution of C. hunanica Chen et Wang are provided.

Launching insectphylo.org; a new hub facilitating construction and use of synthesis molecular phylogenies of insects.

The Holy Grail of an Insect Tree of Life can only be 'discovered' through extensive collaboration among taxon specialists, phylogeneticists and centralized frameworks such as Open Tree of Life, but insufficient effort from stakeholders has so far hampered this promising approach. The resultant unavailability of synthesis phylogenies is an unfortunate situation given the numerous practical usages of phylogenies in the near term and against the backdrop of the ongoing biodiversity crisis. To resolve this issue, we establish a new online hub that centralizes the collation of relevant phylogenetic data and provides the resultant synthesis molecular phylogenies. This is achieved through key developments in a proposed pipeline for the construction of a species-level insect phylogeny. The functionality of the framework is demonstrated through the construction of a highly supported, species-comprehensive phylogeny of Diptera, built from integrated omics data, COI DNA barcodes, and a compiled database of over 100 standardized, published Diptera phylogenies. Machine-readable forms of the phylogeny (and subsets thereof) are publicly available at insectphylo.org, a new public repository for species-comprehensive phylogenies for biological research.