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.

The impact of heavy metal pollution on wild bee communities in smallholder farmlands.

Wild bees provide important pollination services, but they face numerous stressors that threaten them and their ecosystem services. Wild bees can be exposed to heavy metal pollution through the consumption of nectar, pollen, and water, which might cause bee decline. While some studies have measured heavy metal concentrations in honeybees, few studies have monitored heavy metal concentrations in wild bees or explored their potential effects on wild bee communities. To investigate the impact of heavy metal pollution on wild bee communities, heavy metal concentrations, including vanadium (V), chromium (Cr), nickel (Ni), cadmium (Cd), Zinc (Zn) and lead (Pb) in multiple wild bee species were measured. Multiple wild bee species, including: Xylocopa tranquabaroroum, Eucera floralia, Apis cerana, and small bee mixtures (representing multiple small wild bee species) were sampled from 18 sites in Quzhou, Zhejiang Province, China. The findings demonstrated that there were significant differences in heavy metal concentrations among different bee species. The concentrations of V, Zn, Cd, and Pb in X. tranquabaroroum, the largest bee species in this study, were lower than that in the other three sample groups. Furthermore, there were significant negative correlations between heavy metal pollution and wild bee diversity and species richness, but not with abundance. Particularly, there was no significant relationship between heavy metal pollution and the abundance of small bees. Given these worrying findings, monitoring multiple heavy metals in wild bees should be conducted for protecting wild bee diversity and securing their pollination services.

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.

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.

The molecular clock and evolutionary timescales.

The molecular clock provides a valuable means of estimating evolutionary timescales from genetic and biochemical data. Proposed in the early 1960s, it was first applied to amino acid sequences and immunological measures of genetic distances between species. The molecular clock has undergone considerable development over the years, and it retains profound relevance in the genomic era. In this mini-review, we describe the history of the molecular clock, its impact on evolutionary theory, the challenges brought by evidence of evolutionary rate variation among species, and the statistical models that have been developed to account for these heterogeneous rates of genetic change. We explain how the molecular clock can be used to infer rates and timescales of evolution, and we list some of the key findings that have been obtained when molecular clocks have been applied to genomic data. Despite the numerous challenges that it has faced over the decades, the molecular clock continues to offer the most effective method of resolving the details of the evolutionary timescale of the Tree of Life.

The complete mitochondrial genome of Zaomma eriococci (hymenoptera: encyrtidae).

The complete mitochondrial genome of the Zaomma eriococci (Ferrière, 1955) (Hymenoptera: Encyrtidae) was obtained through next-generation sequencing, making the first reported complete mitochondrial genome of the genus Zaomma. The mitochondrial genome is 15,648 bp in length and includes 37 classical eukaryotic mitochondrial genes along with an A + T rich region. All 13 protein-coding genes (PCGs) initiate with typical ATN codons. Of these, 10 PCG genes terminate with TAA, while three terminate with TAG. Additionally, there are 22 tRNA genes, ranging in size from 62 to 70 bp. The maximum likelihood phylogenetic tree was constructed based on 13 PCGs, indicates that Z. eriococci is closely related to Tassonia gloriae. This mitochondrial genome will serve as a valuable molecular resource for species identification, genetic analysis, and comparative genomic studies of Z. eriococci, contributing to the growing collection of mitochondrial genomes within the family Encyrtidae.

Cryptic diversity begets challenges and opportunities in biodiversity research.

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

A Chromosome-level Genome Assembly and Evolution Analysis of Andrena camellia (Hymenoptera: Andrenidae).

Andrena camellia, an effective pollinator of the economicallysignificant crop Camellia oleifera, can withstand the toxic pollen of C. oleifera, making A. camellia a crucial for resource conservation and cultivation of C. oleifera. In this study, the whole genome of A. camellia was sequenced on the Oxford Nanopore platform. The assembled genome size was 340.73 Mb including 50 scaffolds (N50=47.435 Mb) and 131 contigs (N50=17.2 Mb). A total of 11, 258 protein-coding genes were annotated, in addition, 1,104 non-coding RNAs were identified. Further analysis that some chromosomes of A. camellia have a high level of synteny with those of Apis mellifera, Osmia bicornis and Andrena minutula. Thus, our reported genome of A. camellia serves as a valuable resource for studying species evolution, behavioral biology, and adaption to toxic pollen of C. oleifera.

The Whole Genome of Colletes collaris (Hymenoptera: Colletidae): An Important Step in Comparative Genomics of Cellophane Bees.

Gaining knowledge on bees is of the utmost importance due to the paramount role that they play in angiosperm pollination. Herein, we provide the first genome assembly of Colletes collaris, a pan-Eurasian cellophane bee. We sequenced 50.53 Gbp of long-read data plus 57.36 Gbp of short-read data in Oxford Nanopore Technologies and Illumina platforms, respectively. The genome assembly consisted of 374.75 Mbp distributed across 374 contigs, with L50 and N50 of 9 and 8.96 Mbp, respectively. We predicted the genome to comprise 20,399 protein-coding genes, 467,947 repeats, and 4,315 non-coding RNA genes. The transcriptome and mitochondrial genome of the species were also assembled. Gene family analysis with 15 insect species identified 14,417 families, 9,517 of them found in C. collaris. A dated phylogenomic analysis revealed high numbers of orthogroups experiencing rapid evolution within Colletes.

The complete mitochondrial genome of the rice moth, Corcyra cephalonica.

The complete mitochondrial genome (mitogenome) of the rice moth, Corcyra cephalonica Stainton (Lepidoptera: Pyralidae) was determined as a circular molecular of 15,273 bp in size. The mitogenome composition (37 genes) and gene order are the same as the other lepidopterans. Nucleotide composition of the C. cephalonica mitogenome is highly A+T biased (80.43%) like other insects. Twelve protein-coding genes start with a typical ATN codon, with the exception of coxl gene, which uses CGA as the initial codon. Nine protein-coding genes have the common stop codon TAA, and the nad2, cox1, cox2, and nad4 have single T as the incomplete stop codon. 22 tRNA genes demonstrated cloverleaf secondary structure. The mitogenome has several large intergenic spacer regions, the spacer1 between trnQ gene and nad2 gene, which is common in Lepidoptera. The spacer 3 between trnE and trnF includes microsatellite-like repeat regions (AT)18 and (TTAT)(3). The spacer 4 (16 bp) between trnS2 gene and nad1 gene has a motif ATACTAT; another species, Sesamia inferens encodes ATCATAT at the same position, while other lepidopteran insects encode a similar ATACTAA motif. The spacer 6 is A+T rich region, include motif ATAGA and a 20-bp poly(T) stretch and two microsatellite (AT)(9), (AT)(8) elements.

Chromosome-Level Genome Assembly of Anthidium xuezhongi Niu & Zhu, 2020 (Hymenoptera: Apoidea: Megachilidae: Anthidiini).

Anthidiini, a large bee tribe characterized by light-colored maculations, represents nearly 1,000 pollinator species, but no genomes are yet available for this tribe. Here, we report a chromosome-level genome assembly of Anthidium xuezhongi collected from the Tibetan Plateau. Using PacBio long reads and Hi-C data, we assembled a genome of 189.14 Mb with 99.94% of the assembly located in 16 chromosomes. Our assembly contains 23 scaffolds, with the scaffold N50 length of 12.53 Mb, and BUSCO completeness of 98.70% (n = 1,367). We masked 25.98 Mb (13.74%) of the assembly as repetitive elements, identified 385 noncoding RNAs, and predicted 10,820 protein-coding genes (99.20% BUSCO completeness). Gene family evolution analyses identified 9,251 gene families, of which 31 gene families experienced rapid evolution. Interspecific chromosomal variation among A. xuezhongi, Bombus terrestris, and Apis mellifera showed strong chromosomal syntenic relationships. This high-quality genome assembly is a valuable resource for evolutionary and comparative genomic analyses of bees.

Testing the systematic status of Homalictus and Rostrohalictus with weakened cross-vein groups within Halictini (Hymenoptera: Halictidae) using low-coverage whole-genome sequencing.

The halictid genus Lasioglossum, as one of the most species-rich bee groups with persistently contentious subgeneric boundaries, is one of the most challenging bee groups from a systematic standpoint. An enduring question is the relationship of Lasioglossum and Homalictus, whether all halictine bees with weakened distal wing venation comprise one or multiple genera. Here, we analyzed the phylogenetic relationships among the subgroups within Lasioglossum s.l. based on thousands of single-copy orthologs and ultraconserved elements, which were extracted from 23 newly sequenced low-coverage whole genomes alongside a published genome (22 ingroups plus 2 outgroups). Both marker sets provided consistent results across maximum likelihood and coalescent-based species tree approaches. The phylogenetic and topology test results show that the Lasioglossum and Hemihalictus series are reciprocally monophyletic and Homalictus and Rostrohalictus are valid subgenera of Lasioglossum. Consequently, we lower Homalictus to subgenus status within Lasioglossum again, and we also raise Rostrohalictus to subgenus status from its prior synonymy with subgenus Hemihalictus. Lasioglossum przewalskyi is also transferred to the subgenus Hemihalictus. Ultimately, we redefine Lasioglossum to include all halictine bees with weakened distal wing venation.

Potential efficacy of mitochondrial genes for animal DNA barcoding: a case study using eutherian mammals.

BACKGROUND: A well-informed choice of genetic locus is central to the efficacy of DNA barcoding. Current DNA barcoding in animals involves the use of the 5' half of the mitochondrial cytochrome oxidase 1 gene (CO1) to diagnose and delimit species. However, there is no compelling a priori reason for the exclusive focus on this region, and it has been shown that it performs poorly for certain animal groups. To explore alternative mitochondrial barcoding regions, we compared the efficacy of the universal CO1 barcoding region with the other mitochondrial protein-coding genes in eutherian mammals. Four criteria were used for this comparison: the number of recovered species, sequence variability within and between species, resolution to taxonomic levels above that of species, and the degree of mutational saturation. RESULTS: Based on 1,179 mitochondrial genomes of eutherians, we found that the universal CO1 barcoding region is a good representative of mitochondrial genes as a whole because the high species-recovery rate (> 90%) was similar to that of other mitochondrial genes, and there were no significant differences in intra- or interspecific variability among genes. However, an overlap between intra- and interspecific variability was still problematic for all mitochondrial genes. Our results also demonstrated that any choice of mitochondrial gene for DNA barcoding failed to offer significant resolution at higher taxonomic levels. CONCLUSIONS: We suggest that the CO1 barcoding region, the universal DNA barcode, is preferred among the mitochondrial protein-coding genes as a molecular diagnostic at least for eutherian species identification. Nevertheless, DNA barcoding with this marker may still be problematic for certain eutherian taxa and our approach can be used to test potential barcoding loci for such groups.

Positive selection on hemagglutinin and neuraminidase genes of H1N1 influenza viruses.

BACKGROUND: Since its emergence in March 2009, the pandemic 2009 H1N1 influenza A virus has posed a serious threat to public health. To trace the evolutionary path of these new pathogens, we performed a selection-pressure analysis of a large number of hemagglutinin (HA) and neuraminidase (NA) gene sequences of H1N1 influenza viruses from different hosts. RESULTS: Phylogenetic analysis revealed that both HA and NA genes have evolved into five distinct clusters, with further analyses indicating that the pandemic 2009 strains have experienced the strongest positive selection. We also found evidence of strong selection acting on the seasonal human H1N1 isolates. However, swine viruses from North America and Eurasia were under weak positive selection, while there was no significant evidence of positive selection acting on the avian isolates. A site-by-site analysis revealed that the positively selected sites were located in both of the cleaved products of HA (HA1 and HA2), as well as NA. In addition, the pandemic 2009 strains were subject to differential selection pressures compared to seasonal human, North American swine and Eurasian swine H1N1 viruses. CONCLUSIONS: Most of these positively and/or differentially selected sites were situated in the B-cell and/or T-cell antigenic regions, suggesting that selection at these sites might be responsible for the antigenic variation of the viruses. Moreover, some sites were also associated with glycosylation and receptor-binding ability. Thus, selection at these positions might have helped the pandemic 2009 H1N1 viruses to adapt to the new hosts after they were introduced from pigs to humans. Positive selection on position 274 of NA protein, associated with drug resistance, might account for the prevalence of drug-resistant variants of seasonal human H1N1 influenza viruses, but there was no evidence that positive selection was responsible for the spread of the drug resistance of the pandemic H1N1 strains.

The complete mitochondrial genome of Spilonota lechriaspis Meyrick (Lepidoptera: Tortricidae).

We determined the nucleotide sequence of the mitochondrial genome (mtgenome) of Spilonota lechriaspis Meyrick (Lepidoptera: Tortricidae). The entire closed circular molecule is 15,368 bp and contains 37 genes with the typical gene complement and order for lepidopteran mtgenomes. All tRNAs except tRNASer(AGN) can be folded into the typical cloverleaf secondary structures. The protein-coding genes (PCGs) have typical mitochondrial start codons, with the exception of COI, which uses the unusual CGA one as is found in all other Lepidoptera sequenced to date. In addition, six of 13 PCGs harbor the incomplete termination codons, a single T. The A+T-rich region contains some conserved structures that are similar to those found in other lepidopteran mtgenomes, including a structure combining the motif 'ATAGA', a 19-bp poly(T) stretch and three microsatellite (AT)n elements which are part of larger 122+ bp macrorepeats. This is the first report of macrorepeats in a lepidopteran mtgenome.

Review of the bee genus Pseudoanthidium Friese, 1898 (Hymenoptera: Apoidea: Megachilidae: Anthidiini) of China with descriptions of three new species.

The carder bee genus Pseudoanthidium Friese, 1898, is revised from China. Eight species are confirmed to occur in China, including three new species: Pseudoanthidium (Pseudoanthidium) yanruae Niu Zhu, sp. nov., P. (P.) kunesense Niu Zhu, sp. nov., P. (P.) chenggongense Niu Zhu, sp. nov. There is also one new generic assignment and synonymy: Anthidium kryzhanovskii Wu, 1962 is a junior synonym of P. (P.) orientale (Bingham, 1897). Pseudoanthidium (P.) campulodonta (Wu, 1990) is synonymized with P. (P.) tenellum (Mocsry, 1880). Here we provide descriptions for the three new species and an illustrated key to the known Chinese Pseudoanthidium.

The complete mitochondrial genome of Cerceris quinquefasciata (Hymenoptera: Crabronidae).

The complete mitochondrial genome of the Cerceris quinquefasciata (Rossi, 1792) (Hymenoptera: Crabronidae) was obtained via next-generation sequencing. This mitochondrial genome is 16,188 bp in length with 37 classical eukaryotic mitochondrial genes and two A + T-rich region. All the 13 PCGs begin with typical ATN codons. Among them, eleven PCG genes terminate with TAA, two with T-. All of the 22 tRNA genes, ranging from 58 to 72 bp with typical cloverleaf structure except for trnS1, whose dihydrouridine (DHU) arm forms a simple loop. Phylogenetic analysis highly supported Crabronidae is the sister group of anthophila bees.