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.

The complete mitochondrial genome of Parotis chlorochroalis (Hampson, 1912) (Lepidoptera: Crambidae).

The complete mitochondrial genome of the Parotis chlorochroalis was sequenced, revaeling a length of 15239 bp with 37 genes and an A + T-rich region. All c13 PCGs begin with typical ATN codons, except COI gene, which starts with CGA. Eleven genes terminate with TAA, two with T-. All 22 tRNA genes exhibit typical cloverleaf structure except for trnS1 P. chlorochroalis has two relatively conserved intergenic regions and two relatively conserved overlapping regions. Phylogenetic analysis support P. chlorochroalis belongs to subfamily Spilomelinae, the topologies of Crambidae are highly congruent with previous studies. This newly sequences mitochondrial genome provides valuable resources for taxonomic inference and evolutionary studies of genus Parotis.

Intraspecific and interspecific resource partitioning between bumblebee workers and males related to nectar quantity and quality.

Bumblebees are important pollinators for many natural and agricultural systems in temperate regions. Interspecific and intraspecific variation in floral resource preferences have been proposed to influence bumblebee community structure. In particular, sexual dimorphism is a major source of intraspecific niche variation. Although interspecific resource partitioning is well studied, few studies have explored the intraspecific dynamics between workers and males. Here, we report a study on a total of 11 528 workers and 2220 males of 14 bumblebee species recorded over 5 years in the Hengduan Mountains of Southwest China. We first compared the potential for interspecific and intraspecific competition between workers and males using visitation records and resource partitioning indices (overlap index). We then evaluated the influence of nectar traits on flower preference, including nectar volume and the levels of hexose, sucrose and 10 essential amino acids (EAAs). We found that the niche overlap between intraspecific workers and males was higher than that between different species, and temporal overlap alone did not strongly determine diet overlap. Males of most species preferred flowers with high levels of EAAs and hexose, whereas workers of some species preferred flowers with high nectar volume and sucrose levels. This study suggests that there is floral resource partitioning among bumblebee species, and between workers and males, which may play a key role in alleviating interspecific and intraspecific competition. These findings also provide a useful guide for which kinds of plants might be most valuable for bumblebees, especially the understudied males, in this biodiversity hotspot.

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.

Plant diversity enhances ecosystem multifunctionality via multitrophic diversity.

Ecosystem functioning depends on biodiversity at multiple trophic levels, yet relationships between multitrophic diversity and ecosystem multifunctionality have been poorly explored, with studies often focusing on individual trophic levels and functions and on specific ecosystem types. Here, we show that plant diversity can affect ecosystem functioning both directly and by affecting other trophic levels. Using data on 13 trophic groups and 13 ecosystem functions from two large biodiversity experiments-one representing temperate grasslands and the other subtropical forests-we found that plant diversity increases multifunctionality through elevated multitrophic diversity. Across both experiments, the association between multitrophic diversity and multifunctionality was stronger than the relationship between the diversity of individual trophic groups and multifunctionality. Our results also suggest that the role of multitrophic diversity is greater in forests than in grasslands. These findings imply that, to promote sustained ecosystem multifunctionality, conservation planning must consider the diversity of both plants and higher trophic levels.

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.