Efficient Bayesian inference under the multispecies coalescent with migration.

Analyses of genome sequence data have revealed pervasive interspecific gene flow and enriched our understanding of the role of gene flow in speciation and adaptation. Inference of gene flow using genomic data requires powerful statistical methods. Yet current likelihood-based methods involve heavy computation and are feasible for small datasets only. Here, we implement the multispecies-coalescent-with-migration model in the Bayesian program bpp, which can be used to test for gene flow and estimate migration rates, as well as species divergence times and population sizes. We develop Markov chain Monte Carlo algorithms for efficient sampling from the posterior, enabling the analysis of genome-scale datasets with thousands of loci. Implementation of both introgression and migration models in the same program allows us to test whether gene flow occurred continuously over time or in pulses. Analyses of genomic data from Anopheles mosquitoes demonstrate rich information in typical genomic datasets about the mode and rate of gene flow.

Uncovering ghost introgression through genomic analysis of a distinct eastern Asian hickory species.

Ghost introgression, or the transfer of genetic material from extinct or unsampled lineages to sampled species, has attracted much attention. However, conclusive evidence for ghost introgression, especially in plant species, remains scarce. Here, we newly assembled chromosome-level genomes for both Carya sinensis and Carya cathayensis, and additionally re-sequenced the whole genomes of 43 C. sinensis individuals as well as 11 individuals representing 11 diploid hickory species. These genomic datasets were used to investigate the reticulation and bifurcation patterns within the genus Carya (Juglandaceae), with a particular focus on the beaked hickory C. sinensis. By combining the D-statistic and BPP methods, we obtained compelling evidence that supports the occurrence of ghost introgression in C. sinensis from an extinct ancestral hickory lineage. This conclusion was reinforced through the phylogenetic network analysis and a genome scan method VolcanoFinder, the latter of which can detect signatures of adaptive introgression from unknown donors. Our results not only dispel certain misconceptions about the phylogenetic history of C. sinensis but also further refine our understanding of Carya's biogeography via divergence estimates. Moreover, the successful integration of the D-statistic and BPP methods demonstrates their efficacy in facilitating a more precise identification of introgression types.

Bayesian Inference Under the Multispecies Coalescent with Ancient DNA Sequences.

Ancient DNA (aDNA) is increasingly being used to investigate questions such as the phylogenetic relationships and divergence times of extant and extinct species. If aDNA samples are sufficiently old, expected branch lengths (in units of nucleotide substitutions) are reduced relative to contemporary samples. This can be accounted for by incorporating sample ages into phylogenetic analyses. Existing methods that use tip (sample) dates infer gene trees rather than species trees, which can lead to incorrect or biased inferences of the species tree. Methods using a multispecies coalescent (MSC) model overcome these issues. We developed an MSC model with tip dates and implemented it in the program bpp. The method performed well for a range of biologically realistic scenarios, estimating calibrated divergence times and mutation rates precisely. Simulations suggest that estimation precision can be best improved by prioritizing sampling of many loci and more ancient samples. Incorrectly treating ancient samples as contemporary in analyzing simulated data, mimicking a common practice of empirical analyses, led to large systematic biases in model parameters, including divergence times. Two genomic datasets of mammoths and elephants were analyzed, demonstrating the method's empirical utility.

Detection of Ghost Introgression Requires Exploiting Topological and Branch Length Information.

In recent years, the study of hybridization and introgression has made significant progress, with ghost introgression-the transfer of genetic material from extinct or unsampled lineages to extant species-emerging as a key area for research. Accurately identifying ghost introgression, however, presents a challenge. To address this issue, we focused on simple cases involving 3 species with a known phylogenetic tree. Using mathematical analyses and simulations, we evaluated the performance of popular phylogenetic methods, including HyDe and PhyloNet/MPL, and the full-likelihood method, Bayesian Phylogenetics and Phylogeography (BPP), in detecting ghost introgression. Our findings suggest that heuristic approaches relying on site-pattern counts or gene-tree topologies struggle to differentiate ghost introgression from introgression between sampled non-sister species, frequently leading to incorrect identification of donor and recipient species. The full-likelihood method BPP uses multilocus sequence alignments directly-hence taking into account both gene-tree topologies and branch lengths, by contrast, is capable of detecting ghost introgression in phylogenomic datasets. We analyzed a real-world phylogenomic dataset of 14 species of Jaltomata (Solanaceae) to showcase the potential of full-likelihood methods for accurate inference of introgression.

Inferring the Direction of Introgression Using Genomic Sequence Data.

Genomic data are informative about the history of species divergence and interspecific gene flow, including the direction, timing, and strength of gene flow. However, gene flow in opposite directions generates similar patterns in multilocus sequence data, such as reduced sequence divergence between the hybridizing species. As a result, inference of the direction of gene flow is challenging. Here, we investigate the information about the direction of gene flow present in genomic sequence data using likelihood-based methods under the multispecies-coalescent-with-introgression model. We analyze the case of two species, and use simulation to examine cases with three or four species. We find that it is easier to infer gene flow from a small population to a large one than in the opposite direction, and easier to infer inflow (gene flow from outgroup species to an ingroup species) than outflow (gene flow from an ingroup species to an outgroup species). It is also easier to infer gene flow if there is a longer time of separate evolution between the initial divergence and subsequent introgression. When introgression is assumed to occur in the wrong direction, the time of introgression tends to be correctly estimated and the Bayesian test of gene flow is often significant, while estimates of introgression probability can be even greater than the true probability. We analyze genomic sequences from Heliconius butterflies to demonstrate that typical genomic datasets are informative about the direction of interspecific gene flow, as well as its timing and strength.

BPP0974 is a Bordetella parapertussis adhesin expressed in the avirulent phase, implicated in biofilm formation and intracellular survival.

B. parapertussis is a bacterium that causes whooping cough, a severe respiratory infection disease, that has shown an increased incidence in the population. Upon transmission through aerosol droplets, the initial steps of host colonization critically depend on the bacterial adhesins. We here described BPP0974, a B. parapertussis protein that exhibits the typical domain architecture of the large repetitive RTX adhesin family. BPP0974 was found to be retained in the bacterial membrane and secreted into the culture medium. This protein was found overexpressed in the avirulent phase of B. parapertussis, the phenotype proposed for initial host colonization. Interestingly, BPP0974 was found relevant for the biofilm formation as well as involved in the bacterial attachment to and survival within the respiratory epithelial cells. Taken together, our results suggest a role for BPP0974 in the early host colonization and pathogenesis of B. parapertussis.

The Danish Aging and Cognition (DanACo) cohort.

BACKGROUND: With aging populations worldwide, identification of predictors of age-related cognitive decline is becoming increasingly important. The Danish Aging and Cognition Cohort (DanACo) including more than 5000 Danish men was established to investigate predictors of age-related cognitive decline from young adulthood to late mid-life. CONSTRUCTION AND CONTENT: The DanACo cohort was established through two separate data collections with identical designs involving a follow-up examination in late mid-life of men for whom intelligence test scores were available from their mandatory conscription board examination. The cohort consists of 5,183 men born from 1949 through 1961, with a mean age of 20.4 years at baseline and a mean age of 64.4 years at follow-up. The baseline measures consisted of height, weight, intelligence test score and educational level collected at the conscription board examination. The follow-up assessment consisted of a re-administration of the same intelligence test and a comprehensive questionnaire covering socio-demographic factors, lifestyle, and health-related factors. The data were collected in test sessions with up to 24 participants per session. Using the unique personal identification number assigned to all Danes, the cohort has been linked to data from national administrative and health registers for prospectively collected data on socioeconomic and health-related factors. UTILITY AND DISCUSSION: The DanACo cohort has some major strengths compared to existing cognitive aging cohorts such as a large sample size (n = 5,183 men), a validated global measure of cognitive ability, a long retest interval (mean 44.0 years) and the availability of prospectively collected data from registries as well as comprehensive questionnaire data. The main weakness is the low participation rate (14.3%) and that the cohort consists of men only. CONCLUSION: Cognitive decline is a result of a summary of factors across the life-course. The DanACo cohort is characterized by a long retest interval and contains data on a wealth of factors across adult life which is essential to establish evidence on predictors of cognitive decline. Moreover, the size of the cohort ensures sufficient statistical power to identify even relatively weak predictors of cognitive decline.

Estimation of Cross-Species Introgression Rates Using Genomic Data Despite Model Unidentifiability.

Full-likelihood implementations of the multispecies coalescent with introgression (MSci) model treat genealogical fluctuations across the genome as a major source of information to infer the history of species divergence and gene flow using multilocus sequence data. However, MSci models are known to have unidentifiability issues, whereby different models or parameters make the same predictions about the data and cannot be distinguished by the data. Previous studies of unidentifiability have focused on heuristic methods based on gene trees and do not make an efficient use of the information in the data. Here we study the unidentifiability of MSci models under the full-likelihood methods. We characterize the unidentifiability of the bidirectional introgression (BDI) model, which assumes that gene flow occurs in both directions. We derive simple rules for arbitrary BDI models, which create unidentifiability of the label-switching type. In general, an MSci model with k BDI events has 2k unidentifiable modes or towers in the posterior, with each BDI event between sister species creating within-model parameter unidentifiability and each BDI event between nonsister species creating between-model unidentifiability. We develop novel algorithms for processing Markov chain Monte Carlo samples to remove label-switching problems and implement them in the bpp program. We analyze real and synthetic data to illustrate the utility of the BDI models and the new algorithms. We discuss the unidentifiability of heuristic methods and provide guidelines for the use of MSci models to infer gene flow using genomic data.

Bayesian Phylogenetic Inference using Relaxed-clocks and the Multispecies Coalescent.

The multispecies coalescent (MSC) model accommodates both species divergences and within-species coalescent and provides a natural framework for phylogenetic analysis of genomic data when the gene trees vary across the genome. The MSC model implemented in the program bpp assumes a molecular clock and the Jukes-Cantor model, and is suitable for analyzing genomic data from closely related species. Here we extend our implementation to more general substitution models and relaxed clocks to allow the rate to vary among species. The MSC-with-relaxed-clock model allows the estimation of species divergence times and ancestral population sizes using genomic sequences sampled from contemporary species when the strict clock assumption is violated, and provides a simulation framework for evaluating species tree estimation methods. We conducted simulations and analyzed two real datasets to evaluate the utility of the new models. We confirm that the clock-JC model is adequate for inference of shallow trees with closely related species, but it is important to account for clock violation for distant species. Our simulation suggests that there is valuable phylogenetic information in the gene-tree branch lengths even if the molecular clock assumption is seriously violated, and the relaxed-clock models implemented in bpp are able to extract such information. Our Markov chain Monte Carlo algorithms suffer from mixing problems when used for species tree estimation under the relaxed clock and we discuss possible improvements. We conclude that the new models are currently most effective for estimating population parameters such as species divergence times when the species tree is fixed.

Inference of Gene Flow between Species under Misspecified Models.

Genomic sequence data provide a rich source of information about the history of species divergence and interspecific hybridization or introgression. Despite recent advances in genomics and statistical methods, it remains challenging to infer gene flow, and as a result, one may have to estimate introgression rates and times under misspecified models. Here we use mathematical analysis and computer simulation to examine estimation bias and issues of interpretation when the model of gene flow is misspecified in analysis of genomic datasets, for example, if introgression is assigned to the wrong lineages. In the case of two species, we establish a correspondence between the migration rate in the continuous migration model and the introgression probability in the introgression model. When gene flow occurs continuously through time but in the analysis is assumed to occur at a fixed time point, common evolutionary parameters such as species divergence times are surprisingly well estimated. However, the time of introgression tends to be estimated towards the recent end of the period of continuous gene flow. When introgression events are assigned incorrectly to the parental or daughter lineages, introgression times tend to collapse onto species divergence times, with introgression probabilities underestimated. Overall, our analyses suggest that the simple introgression model is useful for extracting information concerning between-specific gene flow and divergence even when the model may be misspecified. However, for reliable inference of gene flow it is important to include multiple samples per species, in particular, from hybridizing species.

Purification and Structure Characterization of the Crude Polysaccharide from the Fruiting Bodies of Butyriboletus pseudospeciosus and Its Modulation Effects on Gut Microbiota.

Polysaccharides from the species of Boletaceae (Boletales, Agaricomycetes, Basidiomycota) are economically significant to both functional foods and medicinal industries. The crude polysaccharide from Butyriboletus pseudospeciosus (BPP) was prepared, and its physicochemical properties were characterized through the use of consecutive experimental apparatus, and its impact on the gut microbiota of Kunming mice was evaluated. Analyses of the structure characteristics revealed that BPP was mainly composed of Man, Glc, and Gal, possessing the pyranose ring and ß/α-glycosidic linkages. TG analysis exhibited that BPP had great heat stability. The SEM observation performed demonstrated that BPP appeared with a rough, dense, and porous shape. Through the BPP intervention, the serum and fecal biochemical index in mice can be improved obviously (p < 0.05). The abundance of beneficial microbiota in the BPP-treated group was significantly increased, while the abundance of harmful microbiota was significantly decreased (p < 0.05). Based on the Tax4Fun, we also revealed the relationship between the species of gut microbiota and showed that the high dose of BPP has significantly changed the functional diversities compared with those in other groups (p < 0.05). The results suggest that B. pseudospeciosus could serve as potential functional food or medicine.

A simulation study to examine the impact of recombination on phylogenomic inferences under the multispecies coalescent model.

Phylogenomic analyses under the multispecies coalescent model assume no recombination within locus and free recombination among loci. Yet, in real data sets intralocus recombination causes different sites of the same locus to have different genealogical histories so that the model is misspecified. The impact of recombination on various coalescent-based phylogenomic analyses has not been systematically examined. Here, we conduct a computer simulation to examine the impact of recombination on several Bayesian analyses of multilocus sequence data, including species tree estimation, species delimitation (by Bayesian selection of delimitation models) and estimation of evolutionary parameters such as species divergence and introgression times, population sizes for modern and extinct species, and cross-species introgression probabilities. We found that recombination, at rates comparable to estimates from the human being, has little impact on coalescent-based species tree estimation, species delimitation and estimation of population parameters. At rates 10 times higher than the human rate, recombination may affect parameter estimation, causing positive biases in introgression times and ancestral population sizes, although species divergence times and cross-species introgression probabilities are estimated with little bias. Overall, the simulation suggests that phylogenomic inferences under the multispecies coalescent model are robust to realistic amounts of intralocus recombination.

Phylogenetic placement and species delimitation of the crab spider genus Phrynarachne (Araneae: Thomisidae) from China.

Evolutionary biologists have long been fascinated by the striking resemblance to bird droppings of the sit-and-wait crab spiders of the genus Phrynarachne. In doing so, species of Phrynarachne have evolved not to avoid detection, but rather, to cause predators to misidentify them as inedible and/or inanimate bird droppings. However, the lack of a phylogeny for Phrynarachne impedes our understanding of the evolution of this trait in the genus. Here we explore species boundaries in species of Phrynarachne from China using single- and multi-locus species delimitation approaches based on 30 Phrynarachne samples. All species delimitation approaches supported six species of Phrynarachne in China. We further present the first phylogenetic analysis of the genus Phrynarachne and estimate divergence times using two mitochondrial and three nuclear genes. All of our phylogenetic analyses supported the monophyly of Phrynarachne in China, with the genus still included within the higher 'Thomisus group' based on our results. Our dating analyses place the crown age of Phrynarachne in China to the middle Miocene. Taken together, our study provides a time-calibrated phylogeny of the genus Phrynarachne in China for testing hypotheses regarding the evolution of the lineage and bird dropping masquerade.

Is There a Correlation Between the Number of Brain Cells and IQ?

Our access to a unique material of postmortem brains obtained from decades of data collection enabled a stereological analysis of the neuron numbers and correlation of results with individual premorbid intelligence quotient (IQ) data. In our sample of 50 brains from men, we find that IQ does not correlate with the number of brain cells in the human neocortex and was only weakly correlated to brain weight. Our stereological examination extended to measures of several other parameters that might be of relevance to intelligence, including numbers of cerebral glial cells (astrocytes, oligodendrocytes, and microglia) and the volume of key areas in the gray and white matter and of the cerebral ventricles, also showing near-zero nonsignificant correlations to IQ.

A Bayesian Implementation of the Multispecies Coalescent Model with Introgression for Phylogenomic Analysis.

Recent analyses suggest that cross-species gene flow or introgression is common in nature, especially during species divergences. Genomic sequence data can be used to infer introgression events and to estimate the timing and intensity of introgression, providing an important means to advance our understanding of the role of gene flow in speciation. Here, we implement the multispecies-coalescent-with-introgression model, an extension of the multispecies-coalescent model to incorporate introgression, in our Bayesian Markov chain Monte Carlo program Bpp. The multispecies-coalescent-with-introgression model accommodates deep coalescence (or incomplete lineage sorting) and introgression and provides a natural framework for inference using genomic sequence data. Computer simulation confirms the good statistical properties of the method, although hundreds or thousands of loci are typically needed to estimate introgression probabilities reliably. Reanalysis of data sets from the purple cone spruce confirms the hypothesis of homoploid hybrid speciation. We estimated the introgression probability using the genomic sequence data from six mosquito species in the Anopheles gambiae species complex, which varies considerably across the genome, likely driven by differential selection against introgressed alleles.

A Simulation Study to Examine the Information Content in Phylogenomic Data Sets under the Multispecies Coalescent Model.

We use computer simulation to examine the information content in multilocus data sets for inference under the multispecies coalescent model. Inference problems considered include estimation of evolutionary parameters (such as species divergence times, population sizes, and cross-species introgression probabilities), species tree estimation, and species delimitation based on Bayesian comparison of delimitation models. We found that the number of loci is the most influential factor for almost all inference problems examined. Although the number of sequences per species does not appear to be important to species tree estimation, it is very influential to species delimitation. Increasing the number of sites and the per-site mutation rate both increase the mutation rate for the whole locus and these have the same effect on estimation of parameters, but the sequence length has a greater effect than the per-site mutation rate for species tree estimation. We discuss the computational costs when the data size increases and provide guidelines concerning the subsampling of genomic data to enable the application of full-likelihood methods of inference.

Bursal peptide BP-IV as a novel immunoadjuvant enhances the protective efficacy of an epitope peptide vaccine containing T and B cell epitopes of the H9N2 avian influenza virus.

Short peptide antigens covering conserved T or B cell epitopes have been investigated in influenza vaccines. Bursal pentapeptide V (BPP-V) and bursal peptide IV (BP-IV) are small molecular peptides that were isolated and identified from the bursa of Fabricius (BF) and induce a strong immune response at both the humoural and cellular levels. To explore the molecular adjuvant potential of BPP-V and BP-IV with an epitope vaccine, an epitope peptide (HA284-298, GNCVVQCQTERGGLN) rich in T and B cell epitopes for the H9N2 avian influenza virus (AIV) haemagglutinin (HA) protein was selected. BPP-V and BP-IV were coupled with the epitope peptide sequence to form BPP-V and BP-IV-epitope vaccines, respectively. The immunoefficacy of BPP-V and BP-IV-epitope peptide vaccines was evaluated. The results showed that the epitope peptide had weak immunogenicity. The BPP-V-epitope peptide vaccine promoted only the secretion of anti-HA IgG and IgG1 antibodies. The BP-IV-epitope peptide vaccine not only promoted the production of anti-HA IgG and IgG1 antibodies but also significantly induced the production of the IgG2a antibody. The BP-IV-epitope peptide vaccine significantly promoted the production of interleukin (IL-4) and interferon-γ (IFN-γ) (the BPP-V epitope peptide vaccine promoted only the production of IL-4), enhanced the cytotoxic T lymphocyte (CTL) response, and increased the proportion of CD3+ T lymphocytes. Moreover, the BP-IV-epitope peptide vaccine promoted a cell-mediated immune response similar to that of the AIV vaccine group. Furthermore, BPP-V and BP-IV-epitope peptide vaccines could also accelerate the clearance of pulmonary virus and reduce pathological damage after the challenge with H9N2 AIV. This study demonstrates the potential of BP-IV as an effective adjuvant for the epitope peptide vaccine for the H9N2 AIV.

Testing species hypotheses for Fridericia magna, an enchytraeid worm (Annelida: Clitellata) with great mitochondrial variation.

BACKGROUND: Deep mitochondrial divergences were observed in Scandinavian populations of the terrestrial to semi-aquatic annelid Fridericia magna (Clitellata: Enchytraeidae). This raised the need for testing whether the taxon is a single species or a complex of cryptic species. RESULTS: A total of 62 specimens from 38 localities were included in the study, 44 of which were used for species delimitation. First, the 44 specimens were divided into clusters using ABGD (Automatic Barcode Gap Discovery) on two datasets, consisting of sequences of the mitochondrial markers COI and 16S. For each dataset, the worms were divided into six not completely congruent clusters. When they were combined, a maximum of seven clusters, or species hypotheses, were obtained, and the seven clusters were used as input in downstream analyses. We tested these hypotheses by constructing haplowebs for two nuclear markers, H3 and ITS, and in both haplowebs the specimens appeared as a single species. Multi-locus species delimitation analyses performed with the Bayesian BPP program also mainly supported a single species. Furthermore, no apparent morphological differences were found between the clusters. Two of the clusters were partially separated from each other and the other clusters, but not strongly enough to consider them as separate species. All 62 specimens were used to visualise the Scandinavian distribution, of the species, and to compare with published COI data from other Fridericia species. CONCLUSION: We show that the morphospecies Fridericia magna is a single species, harbouring several distinct mitochondrial clusters. There is partial genetic separation between some of them, which may be interpreted as incipient speciation. The study shows the importance of rigorous species delimitation using several independent markers when deep mitochondrial divergences might give the false impression of cryptic speciation.

Assessment of branch point prediction tools to predict physiological branch points and their alteration by variants.

BACKGROUND: Branch points (BPs) map within short motifs upstream of acceptor splice sites (3'ss) and are essential for splicing of pre-mature mRNA. Several BP-dedicated bioinformatics tools, including HSF, SVM-BPfinder, BPP, Branchpointer, LaBranchoR and RNABPS were developed during the last decade. Here, we evaluated their capability to detect the position of BPs, and also to predict the impact on splicing of variants occurring upstream of 3'ss. RESULTS: We used a large set of constitutive and alternative human 3'ss collected from Ensembl (n = 264,787 3'ss) and from in-house RNAseq experiments (n = 51,986 3'ss). We also gathered an unprecedented collection of functional splicing data for 120 variants (62 unpublished) occurring in BP areas of disease-causing genes. Branchpointer showed the best performance to detect the relevant BPs upstream of constitutive and alternative 3'ss (99.48 and 65.84% accuracies, respectively). For variants occurring in a BP area, BPP emerged as having the best performance to predict effects on mRNA splicing, with an accuracy of 89.17%. CONCLUSIONS: Our investigations revealed that Branchpointer was optimal to detect BPs upstream of 3'ss, and that BPP was most relevant to predict splicing alteration due to variants in the BP area.

Cryptic species delineation in freshwater planarians of the genus Dugesia (Platyhelminthes, Tricladida): Extreme intraindividual genetic diversity, morphological stasis, and karyological variability.

The keystone of planarian taxonomy traditionally has been the anatomy of the copulatory apparatus. However, many planarian species comprise asexual fissiparous populations, with the fissiparous animals not developing a copulatory apparatus, thus precluding their morphological identification. Incorporation of molecular data into planarian systematics has been of great value, not only in the identification of fissiparous individuals but also as an additional source of information for determining species boundaries. Nevertheless, the discrepancy between morphological and molecular data has highlighted the need for extra sources of taxonomic information. Moreover, a recent study has pointed out that fissiparous reproduction may lead to high levels of intraindividual genetic diversity in planarians, which may mislead molecular analyses. In the present study we aim to test a new up-to-date integrative taxonomic procedure for planarians, including intraindividual genetic data and additional sources of taxonomic information, besides morphology and DNA, using Dugesia subtentaculata sensu lato as a model organism, a species with an intricate taxonomic history. First, we used three different methods for molecular species delimitation on single locus datasets, both with and without intraindividual information, for formulating Primary Species Hypotheses (PSHs). Subsequently, Secondary Species Hypotheses (SSHs) were formulated on the basis of three types of information: (1) a coalescent-based species delimitation method applied to multilocus data, (2) morphology of the copulatory apparatus, and (3) karyological metrics. This resulted in the delimitation of four morphologically cryptic species within the nominal species D. subtentaculata. Our results provide evidence that the analysis of intraindividual genetic data is essential for properly developing PSHs in planarians. Our study reveals also that karyological differentiation, rather than morphological differentiation, may play an important role in speciation processes in planarians, thus suggesting that the currently known diversity of the group could be highly underestimated.