Rooting Trees, Methods for

Torpor energetics are related to the interaction between body mass and climate in bats of the family Vespertilionidae

Torpor is an adaptive strategy allowing heterothermic animals to cope with energy limitations. In birds and mammals, intrinsic and extrinsic factors, such as body mass and ambient temperature, are the main variables influencing torpor use. A theoretical model of the relationship between metabolic rate during torpor and ambient temperature has been proposed. Nevertheless, no empirical attempts have been made to assess the model predictions under different climates. Using open-flow respirometry, we evaluated the ambient temperature at which bats entered torpor and when torpid metabolic rate reached its minimum, the reduction in metabolic rate below basal values, and minimum torpid metabolic rate in 11 bat species of the family Vespertilionidae with different body mass from warm and cold climates. We included data on the minimum torpid metabolic rate of five species we retrieved from the literature. We tested the effects using mixed-effect phylogenetic models. All models showed a significant interaction between body mass and climate. Smaller bats went into torpor and reached minimum torpid metabolic rates at warmer temperatures, showed a higher reduction in the metabolic rate below basal values, and presented lower torpid metabolic rates than larger ones. The slopes of the models were different for bats from different climates. These results are likely explained by differences in body mass and the metabolic rate of bats, which may favor larger bats expressing torpor in colder sites and smaller bats in the warmer ones. Further studies to assess torpor use in bats from different climates are proposed.

Beyond stomatal development: SMF transcription factors as versatile toolkits for land plant evolution

As master transcription factors of stomatal development, SPEECHLESS, MUTE, and FAMA, collectively termed SMFs, are primary targets of molecular genetic analyses in the model plant Arabidopsis thaliana. Studies in other model systems identified SMF orthologs as key players in evolutionary developmental biology studies on stomata. However, recent studies on the astomatous liverwort Marchantia polymorpha revealed that the functions of these genes are not limited to the stomatal development, but extend to other types of tissues, namely sporophytic setal and gametophytic epidermal tissues. These studies provide insightful examples of gene-regulatory network co-opting, and highlight SMFs and related transcription factors as general toolkits for novel trait evolution in land plant lineages. Here, we critically review recent literature on the SMF-like gene in M. polymorpha and discuss their implications for plant evolutionary biology.

Investigation of CaMV-host co-evolution through synonymous codon pattern

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

Community ecology and functional potential of bacteria, archaea, eukarya and viruses in Guerrero Negro microbial mat

In this study, the microbial ecology, potential environmental adaptive mechanisms, and the potential evolutionary interlinking of genes between bacterial, archaeal and viral lineages in Guerrero Negro (GN) microbial mat were investigated using metagenomic sequencing across a vertical transect at millimeter scale. The community composition based on unique genes comprised bacteria (98.01%), archaea (1.81%), eukarya (0.07%) and viruses (0.11%). A gene-focused analysis of bacteria archaea, eukarya and viruses showed a vertical partition of the community. The greatest coverages of genes of bacteria and eukarya were detected in first layers, while the highest coverages of genes of archaea and viruses were found in deeper layers. Many genes potentially related to adaptation to the local environment were detected, such as UV radiation, multidrug resistance, oxidative stress, heavy metals, salinity and desiccation. Those genes were found in bacterial, archaeal and viral lineages with 6477, 44, and 1 genes, respectively. The evolutionary histories of those genes were studied using phylogenetic analysis, showing an interlinking between domains in GN mat.

The impact of the cardiovascular component and somatic mutations on ageing

Mechanistic insight into ageing may empower prolonging the lifespan of humans; however, a complete understanding of this process is still lacking despite a plethora of ageing theories. In order to address this, we investigated the association of lifespan with eight phenotypic traits, that is, litter size, body mass, female and male sexual maturity, somatic mutation, heart, respiratory, and metabolic rate. In support of the somatic mutation theory, we analysed 15 mammalian species and their whole-genome sequencing deriving somatic mutation rate, which displayed the strongest negative correlation with lifespan. All remaining phenotypic traits showed almost equivalent strong associations across this mammalian cohort, however, resting heart rate explained additional variance in lifespan. Integrating somatic mutation and resting heart rate boosted the prediction of lifespan, thus highlighting that resting heart rate may either directly influence lifespan, or represents an epiphenomenon for additional lower-level mechanisms, for example, metabolic rate, that are associated with lifespan.

Spectral top-down recovery of latent tree models

Modeling the distribution of high-dimensional data by a latent tree graphical model is a prevalent approach in multiple scientific domains. A common task is to infer the underlying tree structure, given only observations of its terminal nodes. Many algorithms for tree recovery are computationally intensive, which limits their applicability to trees of moderate size. For large trees, a common approach, termed divide-and-conquer, is to recover the tree structure in two steps. First, separately recover the structure of multiple, possibly random subsets of the terminal nodes. Second, merge the resulting subtrees to form a full tree. Here, we develop spectral top-down recovery (STDR), a deterministic divide-and-conquer approach to infer large latent tree models. Unlike previous methods, STDR partitions the terminal nodes in a non random way, based on the Fiedler vector of a suitable Laplacian matrix related to the observed nodes. We prove that under certain conditions, this partitioning is consistent with the tree structure. This, in turn, leads to a significantly simpler merging procedure of the small subtrees. We prove that STDR is statistically consistent and bound the number of samples required to accurately recover the tree with high probability. Using simulated data from several common tree models in phylogenetics, we demonstrate that STDR has a significant advantage in terms of runtime, with improved or similar accuracy.

The structural characterization and bioactivity assessment of nonspecific lipid transfer protein 1 (nsLTP1) from caraway (Carum carvi) seeds

BACKGROUND

Carum carvi (caraway) of the Apiaceae family has been used in many cultures as a cooking spice and part of the folk medicine. Previous reports primarily focus on the medicinal properties of caraway seed essential oil and the whole seeds extract. However, no effort has been made to study caraway proteins and their potential pharmacological properties, including nonspecific lipid transfer protein (nsLTP), necessitating further research. The current study aimed to characterize nonspecific lipid transfer protein 1 (nsLTP1) from caraway seed, determine its three-dimensional structure, and analyze protein-ligand complex interactions through docking studies. We also evaluated nsLTP1 in vitro cytotoxic effect and antioxidant capacity. Additionally, nsLTP1 thermal- and pH- stability were investigated.

METHODS

Caraway nsLTP1 was purified using two-dimensional chromatography. The complete amino acid sequence of nsLTP1 was achieved by intact protein sequence for the first 20 residues and the overlapping digested peptides. The three-dimensional structure was predicted using MODELLER. Autodock Vina software was employed for docking fatty acids against caraway nsLTP1. Assessment of nsLTP1 cytotoxic activity was achieved by MTS assay, and the Trolox equivalent antioxidant capacity (TAC) was determined. Thermal and pH stability of the nsLTP1 was examined by circular dichroism (CD) spectroscopy.

RESULTS

Caraway nsLTP1 is composed of 91 residues and weighs 9652 Da. The three-dimensional structure of caraway nsLTP1 sequence was constructed based on searching known structures in the PDB. We chose nsLTP of Solanum melongena (PDB ID: 5TVI) as the modeling template with the highest identity among all other homologous proteins. Docking linolenic acid with caraway protein showed a maximum binding score of -3.6 kcal/mol. A preliminary screening of caraway nsLTP1 suppressed the proliferation of human breast cancer cell lines MDA-MB-231 and MCF-7 in a dose‑dependent manner with an IC₅₀ value of 52.93 and 44.76 μM, respectively. Also, nsLTP1 (41.4 μM) showed TAC up to 750.4 μM Trolox equivalent. Assessment of nsLTP1 demonstrated high thermal/pH stability.

CONCLUSION

To the best of our knowledge, this is the first study carried out on nsLTP1 from caraway seeds. We hereby report the sequence of nsLTP1 from caraway seeds and its possible interaction with respective fatty acids using in silico approach. Our data indicated that the protein had anticancer and antioxidant activities and was thermally stable.

Recombination between heterologous human acrocentric chromosomes

The short arms of the human acrocentric chromosomes 13, 14, 15, 21 and 22 (SAACs) share large homologous regions, including ribosomal DNA repeats and extended segmental duplications1,2. Although the resolution of these regions in the first complete assembly of a human genome-the Telomere-to-Telomere Consortium's CHM13 assembly (T2T-CHM13)-provided a model of their homology3, it remained unclear whether these patterns were ancestral or maintained by ongoing recombination exchange. Here we show that acrocentric chromosomes contain pseudo-homologous regions (PHRs) indicative of recombination between non-homologous sequences. Utilizing an all-to-all comparison of the human pangenome from the Human Pangenome Reference Consortium4 (HPRC), we find that contigs from all of the SAACs form a community. A variation graph5 constructed from centromere-spanning acrocentric contigs indicates the presence of regions in which most contigs appear nearly identical between heterologous acrocentric chromosomes in T2T-CHM13. Except on chromosome 15, we observe faster decay of linkage disequilibrium in the pseudo-homologous regions than in the corresponding short and long arms, indicating higher rates of recombination6,7. The pseudo-homologous regions include sequences that have previously been shown to lie at the breakpoint of Robertsonian translocations8, and their arrangement is compatible with crossover in inverted duplications on chromosomes 13, 14 and 21. The ubiquity of signals of recombination between heterologous acrocentric chromosomes seen in the HPRC draft pangenome suggests that these shared sequences form the basis for recurrent Robertsonian translocations, providing sequence and population-based confirmation of hypotheses first developed from cytogenetic studies 50 years ago9.

Roadmap to the study of gene and protein phylogeny and evolution-A practical guide

Developments in sequencing technologies and the sequencing of an ever-increasing number of genomes have revolutionised studies of biodiversity and organismal evolution. This accumulation of data has been paralleled by the creation of numerous public biological databases through which the scientific community can mine the sequences and annotations of genomes, transcriptomes, and proteomes of multiple species. However, to find the appropriate databases and bioinformatic tools for respective inquiries and aims can be challenging. Here, we present a compilation of DNA and protein databases, as well as bioinformatic tools for phylogenetic reconstruction and a wide range of studies on molecular evolution. We provide a protocol for information extraction from biological databases and simple phylogenetic reconstruction using probabilistic and distance methods, facilitating the study of biodiversity and evolution at the molecular level for the broad scientific community.

Trait biases in microbial reference genomes

Common culturing techniques and priorities bias our discovery towards specific traits that may not be representative of microbial diversity in nature. So far, these biases have not been systematically examined. To address this gap, here we use 116,884 publicly available metagenome-assembled genomes (MAGs, completeness ≥80%) from 203 surveys worldwide as a culture-independent sample of bacterial and archaeal diversity, and compare these MAGs to the popular RefSeq genome database, which heavily relies on cultures. We compare the distribution of 12,454 KEGG gene orthologs (used as trait proxies) in the MAGs and RefSeq genomes, while controlling for environment type (ocean, soil, lake, bioreactor, human, and other animals). Using statistical modeling, we then determine the conditional probabilities that a species is represented in RefSeq depending on its genetic repertoire. We find that the majority of examined genes are significantly biased for or against in RefSeq. Our systematic estimates of gene prevalences across bacteria and archaea in nature and gene-specific biases in reference genomes constitutes a resource for addressing these issues in the future.

Topiary: Pruning the manual labor from ancestral sequence reconstruction

Ancestral sequence reconstruction (ASR) is a powerful tool to study the evolution of proteins and thus gain deep insight into the relationships among protein sequence, structure, and function. A major barrier to its broad use is the complexity of the task: it requires multiple software packages, complex file manipulations, and expert phylogenetic knowledge. Here we introduce topiary, a software pipeline that aims to overcome this barrier. To use topiary, users prepare a spreadsheet with a handful of sequences. Topiary then: (1) Infers the taxonomic scope for the ASR study and finds relevant sequences by BLAST; (2) Does taxonomically informed sequence quality control and redundancy reduction; (3) Constructs a multiple sequence alignment; (4) Generates a maximum-likelihood gene tree; (5) Reconciles the gene tree to the species tree; (6) Reconstructs ancestral amino acid sequences; and (7) Determines branch supports. The pipeline returns annotated evolutionary trees, spreadsheets with sequences, and graphical summaries of ancestor quality. This is achieved by integrating modern phylogenetics software (Muscle5, RAxML-NG, GeneRax, and PastML) with online databases (NCBI and the Open Tree of Life). In this paper, we introduce non-expert readers to the steps required for ASR, describe the specific design choices made in topiary, provide a detailed protocol for users, and then validate the pipeline using datasets from a broad collection of protein families. Topiary is freely available for download: https://github.com/harmslab/topiary.

Selective sweep sites and SNP dense regions differentiate Mycobacterium bovis isolates across scales

Mycobacterium bovis, a bacterial zoonotic pathogen responsible for the economically and agriculturally important livestock disease bovine tuberculosis (bTB), infects a broad mammalian host range worldwide. This characteristic has led to bidirectional transmission events between livestock and wildlife species as well as the formation of wildlife reservoirs, impacting the success of bTB control measures. Next Generation Sequencing (NGS) has transformed our ability to understand disease transmission events by tracking variant sites, however the genomic signatures related to host adaptation following spillover, alongside the role of other genomic factors in the M. bovis transmission process are understudied problems. We analyzed publicly available M. bovis datasets collected from 700 hosts across three countries with bTB endemic regions (United Kingdom, United States, and New Zealand) to investigate if genomic regions with high SNP density and/or selective sweep sites play a role in Mycobacterium bovis adaptation to new environments (e.g., at the host-species, geographical, and/or sub-population levels). A simulated M. bovis alignment was created to generate null distributions for defining genomic regions with high SNP counts and regions with selective sweeps evidence. Random Forest (RF) models were used to investigate evolutionary metrics within the genomic regions of interest to determine which genomic processes were the best for classifying M. bovis across ecological scales. We identified in the M. bovis genomes 14 and 132 high SNP density and selective sweep regions, respectively. Selective sweep regions were ranked as the most important in classifying M. bovis across the different scales in all RF models. SNP dense regions were found to have high importance in the badger and cattle specific RF models in classifying badger derived isolates from livestock derived ones. Additionally, the genes detected within these genomic regions harbor various pathogenic functions such as virulence and immunogenicity, membrane structure, host survival, and mycobactin production. The results of this study demonstrate how comparative genomics alongside machine learning approaches are useful to investigate further the nature of M. bovis host-pathogen interactions.

Quintet Rooting: rooting species trees under the multi-species coalescent model

Motivation

Rooted species trees are a basic model with multiple applications throughout biology, including understanding adaptation, biodiversity, phylogeography and co-evolution. Because most species tree estimation methods produce unrooted trees, methods for rooting these trees have been developed. However, most rooting methods either rely on prior biological knowledge or assume that evolution is close to clock-like, which is not usually the case. Furthermore, most prior rooting methods do not account for biological processes that create discordance between gene trees and species trees.

Results

We present Quintet Rooting (QR), a method for rooting species trees based on a proof of identifiability of the rooted species tree under the multi-species coalescent model established by Allman, Degnan and Rhodes (J. Math. Biol., 2011). We show that QR is generally more accurate than other rooting methods, except under extreme levels of gene tree estimation error.

Availability and implementation

Quintet Rooting is available in open source form at https://github.com/ytabatabaee/Quintet-Rooting. The simulated datasets used in this study are from a prior study and are available at https://www.ideals.illinois.edu/handle/2142/55319. The biological dataset used in this study is also from a prior study and is available at http://gigadb.org/dataset/101041.

Contact

warnow@illinois.edu

Supplementary information

Supplementary data are available at Bioinformatics online.

The rates of global bacterial and archaeal dispersal

The phylogenetic resolution at which microorganisms display geographic endemism, the rates at which they disperse at global scales, and the role of humans on global microbial dispersal are largely unknown. Answering these questions is necessary for interpreting microbial biogeography, ecology, and macroevolution and for predicting the spread of emerging pathogenic strains. To resolve these questions, I analyzed the geographic and evolutionary relationships between 36,795 bacterial and archaeal ("prokaryotic") genomes from ∼7000 locations around the world. I find clear signs of continental-scale endemism, including strong correlations between phylogenetic divergence and geographic distance. However, the phylogenetic scale at which endemism generally occurs is extremely small, and most "species" (defined by an average nucleotide identity ≥ 95%) and even closely related strains (average nucleotide identity ≥ 99.9%) are globally distributed. Human-associated lineages display faster dispersal rates than other terrestrial lineages; the average net distance between any two human-associated cell lineages diverging 50 years ago is roughly 580 km. These results suggest that many previously reported global-scale microbial biogeographical patterns are likely the result of recent or current environmental filtering rather than geographic endemism. For human-associated lineages, estimated transition rates between Europe and North America are particularly high, and much higher than for non-human associated terrestrial lineages, highlighting the role that human movement plays in global microbial dispersal. Dispersal was slowest for hot spring- and terrestrial subsurface-associated lineages, indicating that these environments may act as "isolated islands" of microbial evolution.

Standard Candles for Dating Microbial Lineages

Molecular clock analyses are challenging for microbial phylogenies, due to a lack of fossil calibrations that can reliably provide absolute time constraints. An alternative source of temporal constraints for microbial groups is provided by the inheritance of proteins that are specific for the utilization of eukaryote-derived substrates, which have often been dispersed across the Tree of Life via horizontal gene transfer. In particular, animal, algal, and plant-derived substrates are often produced by groups with more precisely known divergence times, providing an older-bound on their availability within microbial environments. Therefore, these ages can serve as "standard candles" for dating microbial groups across the Tree of Life, expanding the reach of informative molecular clock investigations. Here, we formally develop the concept of substrate standard candles and describe how they can be propagated and applied using both microbial species trees and individual gene family phylogenies. We also provide detailed evaluations of several candidate standard candles and discuss their suitability in light of their often complex evolutionary and metabolic histories.

Phylogenetic Analysis Indicates That Evasin-Like Proteins of Ixodid Ticks Fall Into Three Distinct Classes

Chemokines are structurally related proteins that activate leucocyte migration in response to injury or infection. Tick saliva contains chemokine-binding proteins or evasins which likely neutralize host chemokine function and inflammation. Biochemical characterisation of 50 evasins from Ixodes, Amblyomma and Rhipicephalus shows that they fall into two functional classes, A and B, with exclusive binding to either CC- or CXC- chemokines, respectively. Class A evasins, EVA1 and EVA4 have a four-disulfide-bonded core, whereas the class B evasin EVA3 has a three-disulfide-bonded “knottin” structure. All 29 class B evasins have six cysteine residues conserved with EVA3, arrangement of which defines a Cys6-motif. Nineteen of 21 class A evasins have eight cysteine residues conserved with EVA1/EVA4, the arrangement of which defines a Cys8-motif. Two class A evasins from Ixodes (IRI01, IHO01) have less than eight cysteines. Many evasin-like proteins have been identified in tick salivary transcriptomes, but their phylogenetic relationship with respect to biochemically characterized evasins is not clear. Here, using BLAST searches of tick transcriptomes with biochemically characterized evasins, we identify 292 class A and 157 class B evasins and evasin-like proteins from Prostriate (Ixodes), and Metastriate (Amblyomma, Dermacentor, Hyalomma, Rhipicephalus) ticks. Phylogenetic analysis shows that class A evasins/evasin-like proteins segregate into two classes, A1 and A2. Class A1 members are exclusive to Metastriate ticks and typically have a Cys8-motif and include EVA1 and EVA4. Class A2 members are exclusive to Prostriate ticks, lack the Cys8-motif, and include IHO01 and IRI01. Class B evasins/evasin-like proteins are present in both Prostriate and Metastriate lineages, typically have a Cys6-motif, and include EVA3. Most evasins/evasin-like proteins in Metastriate ticks belong to class A1, whereas in Prostriate species they are predominantly class B. In keeping with this, the majority of biochemically characterized Metastriate evasins bind CC-chemokines, whereas the majority of Prostriate evasins bind CXC-chemokines. While the origin of the structurally dissimilar classes A1 and A2 is yet unresolved, these results suggest that class B evasin-like proteins arose before the divergence of Prostriate and Metastriate lineages and likely functioned to neutralize CXC-chemokines and support blood feeding.

Fundamental evolution of all Orthocoronavirinae including three deadly lineages descendent from Chiroptera-hosted coronaviruses: SARS-CoV, MERS-CoV and SARS-CoV-2

The severe acute respiratory syndrome coronavirus (SARS-CoV) emerged in humans in 2002. Despite reports showing Chiroptera as the original animal reservoir of SARS-CoV, many argue that Carnivora-hosted viruses are the most likely origin. The emergence of the Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012 also involves Chiroptera-hosted lineages. However, factors such as the lack of comprehensive phylogenies hamper our understanding of host shifts once MERS-CoV emerged in humans and Artiodactyla. Since 2019, the origin of SARS-CoV-2, causative agent of coronavirus disease 2019 (COVID-19), added to this episodic history of zoonotic transmission events. Here we introduce a phylogenetic analysis of 2006 unique and complete genomes of different lineages of Orthocoronavirinae. We used gene annotations to align orthologous sequences for total evidence analysis under the parsimony optimality criterion. Deltacoronavirus and Gammacoronavirus were set as outgroups to understand spillovers of Alphacoronavirus and Betacoronavirus among ten orders of animals. We corroborated that Chiroptera-hosted viruses are the sister group of SARS-CoV, SARS-CoV-2 and MERS-related viruses. Other zoonotic events were qualified and quantified to provide a comprehensive picture of the risk of coronavirus emergence among humans. Finally, we used a 250 SARS-CoV-2 genomes dataset to elucidate the phylogenetic relationship between SARS-CoV-2 and Chiroptera-hosted coronaviruses.

Evolutionary and functional analysis of an NRPS condensation domain integrates β-lactam, ᴅ-amino acid, and dehydroamino acid synthesis

Nonribosomal peptide synthetases (NRPSs) are large, multidomain biosynthetic enzymes involved in the assembly-line-like synthesis of numerous peptide natural products. Among these are clinically useful antibiotics including three classes of β-lactams: the penicillins/cephalosporins, the monobactams, and the monocyclic nocardicins, as well as the vancomycin family of glycopeptides and the depsipeptide daptomycin. During NRPS synthesis, peptide bond formation is catalyzed by condensation (C) domains, which couple the nascent peptide with the next programmed amino acid of the sequence. A growing number of additional functions are linked to the activity of C domains. In the biosynthesis of the nocardicins, a specialized C domain prepares the embedded β-lactam ring from a serine residue. Here, we examine the evolutionary descent of this unique β-lactam-synthesizing C domain. Guided by its ancestry, we predict and demonstrate in vitro that this C domain alternatively performs peptide bond formation when a single stereochemical change is introduced into its peptide starting material. Remarkably, the function of the downstream thioesterase (TE) domain also changes. Natively, the TE directs C terminus epimerization prior to hydrolysis when the β-lactam is made but catalyzes immediate release of the alternative peptide. In addition, we investigate the roles of C-domain histidine residues in light of clade-specific sequence motifs, refining earlier mechanistic proposals of both β-lactam formation and canonical peptide synthesis. Finally, expanded phylogenetic analysis reveals unifying connections between β-lactam synthesis and allied C domains associated with the appearance of ᴅ-amino acid and dehydroamino acid residues in other NRPS-derived natural products.

Molecular genomic features associated with in vitro response of the NCI-60 cancer cell line panel to natural products

Natural products remain a significant source of anticancer chemotherapeutics. The search for targeted drugs for cancer treatment includes consideration of natural products, which may provide new opportunities for antitumor cytotoxicity as single agents or in combination therapy. We examined the association of molecular genomic features in the well-characterized NCI-60 cancer cell line panel with in vitro response to treatment with 1302 small molecules which included natural products, semisynthetic natural product derivatives, and synthetic compounds based on a natural product pharmacophore from the Developmental Therapeutics Program of the US National Cancer Institute's database. These compounds were obtained from a variety of plant, marine, and microbial species. Molecular information utilized for the analysis included expression measures for 23059 annotated transcripts, lncRNAs, and miRNAs, and data on protein-changing single nucleotide variants in 211 cancer-related genes. We found associations of expression of multiple genes including SLFN11, CYP2J2, EPHX1, GPC1, ELF3, and MGMT involved in DNA damage repair, NOTCH family members, ABC and SLC transporters, and both mutations in tyrosine kinases and BRAF V600E with NCI-60 responses to specific categories of natural products. Hierarchical clustering identified groups of natural products, which correlated with a specific mechanism of action. Specifically, several natural product clusters were associated with SLFN11 gene expression, suggesting that potential action of these compounds may involve DNA damage. The associations between gene expression or genome alterations of functionally relevant genes with the response of cancer cells to natural products provide new information about potential mechanisms of action of these identified clusters of compounds with potentially similar biological effects. This information will assist in future drug discovery and in design of new targeted cancer chemotherapy agents.

FLAVi: An Enhanced Annotator for Viral Genomes of Flaviviridae

Responding to the ongoing and severe public health threat of viruses of the family Flaviviridae, including dengue, hepatitis C, West Nile, yellow fever, and Zika, demands a greater understanding of how these viruses emerge and spread. Updated phylogenies are central to this understanding. Most cladograms of Flaviviridae focus on specific lineages and ignore outgroups, hampering the efficacy of the analysis to test ingroup monophyly and relationships. This is due to the lack of annotated Flaviviridae genomes, which has gene content variation among genera. This variation makes analysis without partitioning difficult. Therefore, we developed an annotation pipeline for the genera of Flaviviridae (Flavirirus, Hepacivirus, Pegivirus, and Pestivirus, named "Fast Loci Annotation of Viruses" (FLAVi; http://flavi-web.com/), that combines ab initio and homology-based strategies. FLAVi recovered 100% of the genes in Flavivirus and Hepacivirus genomes. In Pegivirus and Pestivirus, annotation efficiency was 100% except for one partition each. There were no false positives. The combined phylogenetic analysis of multiple genes made possible by annotation has clear impacts over the tree topology compared to phylogenies that we inferred without outgroups or data partitioning. The final tree is largely congruent with previous hypotheses and adds evidence supporting the close phylogenetic relationship between dengue and Zika.

Assessing the accuracy of phylogenetic rooting methods on prokaryotic gene families

Almost all standard phylogenetic methods for reconstructing gene trees result in unrooted trees; yet, many of the most useful applications of gene trees require that the gene trees be correctly rooted. As a result, several computational methods have been developed for inferring the root of unrooted gene trees. However, the accuracy of such methods has never been systematically evaluated on prokaryotic gene families, where horizontal gene transfer is often one of the dominant evolutionary events driving gene family evolution. In this work, we address this gap by conducting a thorough comparative evaluation of five different rooting methods using large collections of both simulated and empirical prokaryotic gene trees. Our simulation study is based on 6000 true and reconstructed gene trees on 100 species and characterizes the rooting accuracy of the four methods under 36 different evolutionary conditions and 3 levels of gene tree reconstruction error. The empirical study is based on a large, carefully designed data set of 3098 gene trees from 504 bacterial species (406 Alphaproteobacteria and 98 Cyanobacteria) and reveals insights that supplement those gleaned from the simulation study. Overall, this work provides several valuable insights into the accuracy of the considered methods that will help inform the choice of rooting methods to use when studying microbial gene family evolution. Among other findings, this study identifies parsimonious Duplication-Transfer-Loss (DTL) rooting and Minimal Ancestor Deviation (MAD) rooting as two of the most accurate gene tree rooting methods for prokaryotes and specifies the evolutionary conditions under which these methods are most accurate, demonstrates that DTL rooting is highly sensitive to high evolutionary rates and gene tree error, and that rooting methods based on branch-lengths are generally robust to gene tree reconstruction error.

From pairs of most similar sequences to phylogenetic best matches

BACKGROUND

Many of the commonly used methods for orthology detection start from mutually most similar pairs of genes (reciprocal best hits) as an approximation for evolutionary most closely related pairs of genes (reciprocal best matches). This approximation of best matches by best hits becomes exact for ultrametric dissimilarities, i.e., under the Molecular Clock Hypothesis. It fails, however, whenever there are large lineage specific rate variations among paralogous genes. In practice, this introduces a high level of noise into the input data for best-hit-based orthology detection methods.

RESULTS

If additive distances between genes are known, then evolutionary most closely related pairs can be identified by considering certain quartets of genes provided that in each quartet the outgroup relative to the remaining three genes is known. A priori knowledge of underlying species phylogeny greatly facilitates the identification of the required outgroup. Although the workflow remains a heuristic since the correct outgroup cannot be determined reliably in all cases, simulations with lineage specific biases and rate asymmetries show that nearly perfect results can be achieved. In a realistic setting, where distances data have to be estimated from sequence data and hence are noisy, it is still possible to obtain highly accurate sets of best matches.

CONCLUSION

Improvements of tree-free orthology assessment methods can be expected from a combination of the accurate inference of best matches reported here and recent mathematical advances in the understanding of (reciprocal) best match graphs and orthology relations.

AVAILABILITY

Accompanying software is available at https://github.com/david-schaller/AsymmeTree.

Global phylogenetic relationships, population structure and gene flow estimation of Trialeurodes vaporariorum (Greenhouse whitefly)

Trialeurodes vaporariorum (Westwood, 1856) (Greenhouse whitefly) is an agricultural pest of global importance. It is associated with damage to plants during feeding and subsequent virus transmission. Yet, global phylogenetic relationships, population structure, and estimation of the rates of gene flow within this whitefly species remain largely unexplored. In this study, we obtained and filtered 227 GenBank records of mitochondrial cytochrome c oxidase I (mtCOI) sequences of T. vaporariorum, across various global locations to obtain a final set of 217 GenBank records. We further amplified and sequenced a ~750 bp fragment of mtCOI from an additional 31 samples collected from Kenya in 2014. Based on a total of 248 mtCOI sequences, we identified 16 haplotypes, with extensive overlap across all countries. Population structure analysis did not suggest population differentiation. Phylogenetic analysis indicated the 2014 Kenyan collection of samples clustered with a single sequence from the Netherlands to form a well-supported clade (denoted clade 1a) nested within the total set of sequences (denoted clade 1). Pairwise distances between sequences show greater sequence divergence between clades than within clades. In addition, analysis using migrate-n gave evidence for recent gene flow between the two groups. Overall, we find that T. vaporariorum forms a single large group, with evidence of further diversification consisting primarily of Kenyan sequences and one sequence from the Netherlands forming a well-supported clade.

Phylogeography and Coevolution of Bamboo Mosaic Virus and Its Associated Satellite RNA

Bamboo mosaic virus (BaMV), a plant potexvirus, has been found only in infected bamboo species. It is frequently associated with a large, linear single-stranded satellite RNA (satBaMV) that encodes a non-structural protein. Decades of collecting across a wide geographic area in Asia have accumulated a sizable number of BaMV and satBaMV isolates. In this study, we reconstructed the BaMV phylogeny and satBaMV phylogeny with partial coat protein gene sequences and partial genomic sequences, respectively. The evolutionary relationships allowed us to infer the phylogeography of BaMV and satBaMV on the Asian continent and its outlying islands. The BaMV phylogeny suggests that the BaMV isolates from Taiwan, unsurprisingly, are most likely derived from China. Interestingly, the newly available satBaMV isolates from China were found to be most closely related to the previously established Clade III, which is found in India. The general pattern of clustering along the China/India and Taiwan divide led us to hypothesize that the Taiwan Strait has been a physical barrier to gene flow in the past evolutionary history of both BaMV and satBaMV. Lastly, cophylogeny analyses revealed a complex association pattern between BaMV and satBaMV isolates from China. In general, closely related BaMV sequences tend to carry closely related satBaMV sequences as well; but instances of mismatching with distantly related satBaMV isolates were also found. We hypothesize plausible scenarios of infection and superinfection of bamboo hosts that may be responsible for the observed association pattern. However, a more systematic sampling throughout the geographic distribution of various bamboo species is needed to unambiguously establish the origin, movement, and evolution of BaMV and satBaMV.