Treasurer’s Report for Financial Year (FY) 2018Prepared by Jeffrey L. Thorne, SMBE Treasurer Audited by Dr John Archibald

Wed, 08 Apr 2020 00:00:00 GMT

Society’s net assets as of May 31, 2019:

Molecular Evolution of Auxin-Mediated Root Initiation in Plants

Sun, 05 Apr 2020 00:00:00 GMT

The root originated independently in euphyllophytes (ferns and seed plants) and lycophytes; however, the molecular evolutionary route of root initiation remains elusive. By analyses of the fern Ceratopteris richardii and seed plants, here we show that the molecular pathway involving auxin, intermediate-clade WUSCHEL-RELATED HOMEOBOX (IC-WOX) genes, and WUSCHEL-clade WOX (WC-WOX) genes could be conserved in root initiation. We propose that the “auxin>IC-WOX>WC-WOX” module in root initiation might have arisen in the common ancestor of euphyllophytes during the second origin of roots, and that this module has further developed during the evolution of different root types in ferns and seed plants.

Infectious Disease Defenses among Primary Ancient Hominid Contributions to Adaptation of Modern Humans

Thu, 27 Feb 2020 00:00:00 GMT

During the past decade, our human evolutionary tree has turned into something more resembling an unwieldy bush. Scientists have discovered swapped segments of DNA that we shared from mating between two other hominids, Neanderthals and Denisovans, which were first sequenced in 2010 and 2014, respectively.

Back from the Deep: Scientists Use “Resurrection Ecology” of 600-Year-Old Water Fleas to Better Understand How Organisms Adapt to Lake Pollution

Thu, 27 Feb 2020 00:00:00 GMT

One of the leading threats to lakes since the rise of agriculture is runoffs from fertilizer, in the form of high phosphorus levels. These can trigger devastating events like eutrophication, where deadly algal blooms thrive on phosphorus, and in the process, outcompete and choke off vital nutrients from the rest of the lake.

CellCoal: Coalescent Simulation of Single-Cell Sequencing Samples

Thu, 06 Feb 2020 00:00:00 GMT

Our capacity to study individual cells has enabled a new level of resolution for understanding complex biological systems such as multicellular organisms or microbial communities. Not surprisingly, several methods have been developed in recent years with a formidable potential to investigate the somatic evolution of single cells in both healthy and pathological tissues. However, single-cell sequencing data can be quite noisy due to different technical biases, so inferences resulting from these new methods need to be carefully contrasted. Here, I introduce CellCoal, a software tool for the coalescent simulation of single-cell sequencing genotypes. CellCoal simulates the history of single-cell samples obtained from somatic cell populations with different demographic histories and produces single-nucleotide variants under a variety of mutation models, sequencing read counts, and genotype likelihoods, considering allelic imbalance, allelic dropout, amplification, and sequencing errors, typical of this type of data. CellCoal is a flexible tool that can be used to understand the implications of different somatic evolutionary processes at the single-cell level, and to benchmark dedicated bioinformatic tools for the analysis of single-cell sequencing data. CellCoal is available at https://github.com/dapogon/cellcoal.

IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era

Mon, 03 Feb 2020 00:00:00 GMT

IQ-TREE (http://www.iqtree.org, last accessed February 6, 2020) is a user-friendly and widely used software package for phylogenetic inference using maximum likelihood. Since the release of version 1 in 2014, we have continuously expanded IQ-TREE to integrate a plethora of new models of sequence evolution and efficient computational approaches of phylogenetic inference to deal with genomic data. Here, we describe notable features of IQ-TREE version 2 and highlight the key advantages over other software.

Evolution of the Autism-Associated Neuroligin-4 Gene Reveals Broad Erosion of Pseudoautosomal Regions in Rodents

Mon, 03 Feb 2020 00:00:00 GMT

Variants in genes encoding synaptic adhesion proteins of the neuroligin family, most notably neuroligin-4, are a significant cause of autism spectrum disorders in humans. Although human neuroligin-4 is encoded by two genes, NLGN4X and NLGN4Y, that are localized on the X-specific and male-specific regions of the two sex chromosomes, the chromosomal localization and full genomic sequence of the mouse Nlgn4 gene remain elusive. Here, we analyzed the neuroligin-4 genes of numerous rodent species by direct sequencing and bioinformatics, generated complete drafts of multiple rodent neuroligin-4 genes, and examined their evolution. Surprisingly, we find that the murine Nlgn4 gene is localized to the pseudoautosomal region (PAR) of the sex chromosomes, different from its human orthologs. We show that the sequence differences between various neuroligin-4 proteins are restricted to hotspots in which rodent neuroligin-4 proteins contain short repetitive sequence insertions compared with neuroligin-4 proteins from other species, whereas all other protein sequences are highly conserved. Evolutionarily, these sequence insertions initiate in the clade eumuroidea of the infraorder myomorpha and are additionally associated with dramatic changes in noncoding sequences and gene size. Importantly, these changes are not exclusively restricted to neuroligin-4 genes but reflect major evolutionary changes that substantially altered or even deleted genes from the PARs of both sex chromosomes. Our results show that despite the fact that the PAR in rodents and the neuroligin-4 genes within the rodent PAR underwent massive evolutionary changes, neuroligin-4 proteins maintained a highly conserved core structure, consistent with a substantial evolutionary pressure preserving its physiological function.

Independent Losses of the Hypoxia-Inducible Factor (HIF) Pathway within Crustacea

Fri, 31 Jan 2020 00:00:00 GMT

Metazoans respond to hypoxic stress via the hypoxia-inducible factor (HIF) pathway, a mechanism thought to be extremely conserved due to its importance in monitoring cellular oxygen levels and regulating responses to hypoxia. However, recent work revealed that key members of the HIF pathway have been lost in specific lineages (a tardigrade and a copepod), suggesting that this pathway is not as widespread in animals as previously assumed. Using genomic and transcriptomic data from 70 different species across 12 major crustacean groups, we assessed the degree to which the gene HIFα, the master regulator of the HIF pathway, was conserved. Mining of protein domains, followed by phylogenetic analyses of gene families, uncovered group-level losses of HIFα, including one across three orders within Cirripedia, and in three orders within Copepoda. For these groups, additional assessment showed losses of HIF repression machinery (EGLN and VHL). These results suggest the existence of alternative mechanisms for cellular response to low oxygen and highlight these taxa as models useful for probing these evolutionary outcomes.

Herpes Simplex Viruses: New Relationships between Epidemiology and History

Thu, 30 Jan 2020 00:00:00 GMT

Scientists have revised the dating of viral dispersal from Africa: it did not occur during ancient “Out of Africa” migrations, but rather, more recent events, including the transatlantic slave trade of the eighteenth century.

Evolution of a New Function by Fusion between Phage DNA and a Bacterial Gene

Tue, 21 Jan 2020 00:00:00 GMT

Mobile genetic elements, such as plasmids, phages, and transposons, are important sources for evolution of novel functions. In this study, we performed a large-scale screening of metagenomic phage libraries for their ability to suppress temperature-sensitivity in Salmonella enterica serovar Typhimurium strain LT2 mutants to examine how phage DNA could confer evolutionary novelty to bacteria. We identified an insert encoding 23 amino acids from a phage that when fused with a bacterial DNA-binding repressor protein (LacI) resulted in the formation of a chimeric protein that localized to the outer membrane. This relocalization of the chimeric protein resulted in increased membrane vesicle formation and an associated suppression of the temperature sensitivity of the bacterium. Both the host LacI protein and the extracellular 23-amino acid stretch are necessary for the generation of the novel phenotype. Furthermore, mutational analysis of the chimeric protein showed that although the native repressor function of the LacI protein is maintained in this chimeric structure, it is not necessary for the new function. Thus, our study demonstrates how a gene fusion between foreign DNA and bacterial DNA can generate novelty without compromising the native function of a given gene.

The Genomes of the Livebearing Fish Species Poeciliopsis retropinna and Poeciliopsis turrubarensis Reflect Their Different Reproductive Strategies

Tue, 21 Jan 2020 00:00:00 GMT

The evolution of a placenta is predicted to be accompanied by rapid evolution of genes involved in processes that regulate mother–offspring interactions during pregnancy, such as placenta formation, embryonic development, and nutrient transfer to offspring. However, these predictions have only been tested in mammalian species, where only a single instance of placenta evolution has occurred. In this light, the genus Poeciliopsis is a particularly interesting model for placenta evolution, because in this genus a placenta has evolved independently from the mammalian placenta. Here, we present and compare genome assemblies of two species of the livebearing fish genus Poeciliopsis (family Poeciliidae) that differ in their reproductive strategy: Poeciliopsis retropinna which has a well-developed complex placenta and P. turrubarensis which lacks a placenta. We applied different assembly strategies for each species: PacBio sequencing for P. retropinna (622-Mb assembly, scaffold N50 of 21.6 Mb) and 10× Genomics Chromium technology for P. turrubarensis (597-Mb assembly, scaffold N50 of 4.2 Mb). Using the high contiguity of these genome assemblies and near-completeness of gene annotations to our advantage, we searched for gene duplications and performed a genome-wide scan for genes evolving under positive selection. We find rapid evolution in major parts of several molecular pathways involved in parent–offspring interaction in P. retropinna, both in the form of gene duplications as well as positive selection. We conclude that the evolution of the placenta in the genus Poeciliopsis is accompanied by rapid evolution of genes involved in similar genomic pathways as found in mammals.

Evolutionary Dynamics of Abundant 7-bp Satellites in the Genome of Drosophila virilis

Tue, 21 Jan 2020 00:00:00 GMT

The factors that drive the rapid changes in abundance of tandem arrays of highly repetitive sequences, known as satellite DNA, are not well understood. Drosophila virilis has one of the highest relative amounts of simple satellites of any organism that has been studied, with an estimated >40% of its genome composed of a few related 7-bp satellites. Here, we use D. virilis as a model to understand technical biases affecting satellite sequencing and the evolutionary processes that drive satellite composition. By analyzing sequencing data from Illumina, PacBio, and Nanopore platforms, we identify platform-specific biases and suggest best practices for accurate characterization of satellites by sequencing. We use comparative genomics and cytogenetics to demonstrate that the highly abundant AAACTAC satellite family arose from a related satellite in the branch leading to the virilis phylad 4.5–11 Ma before exploding in abundance in some species of the clade. The most abundant satellite is conserved in sequence and location in the pericentromeric region but has diverged widely in abundance among species, whereas the satellites nearest the centromere are rapidly turning over in sequence composition. By analyzing multiple strains of D. virilis, we saw that the abundances of two centromere-proximal satellites are anticorrelated along a geographical gradient, which we suggest could be caused by ongoing conflicts at the centromere. In conclusion, we illuminate several key attributes of satellite evolutionary dynamics that we hypothesize to be driven by processes including selection, meiotic drive, and constraints on satellite sequence and abundance.

Evolution of Transcriptional Repressors Impacts Caenorhabditis Vulval Development

Tue, 21 Jan 2020 00:00:00 GMT

Comparative genomic sequence analysis has found that the genes for many chromatin-associated proteins are poorly conserved, but the biological consequences of these sequence changes are not understood. Here, we show that four genes identified for an Inappropriate Vulval cell Proliferation (ivp) phenotype in the nematode Caenorhabditis briggsae exhibit distinct functions and genetic interactions when compared with their orthologs in C. elegans. Specifically, we show that the four C. briggsae ivp genes encode the noncanonical histone HTZ-1/H2A.z and three nematode-specific proteins predicted to function in the nucleus. The mutants exhibit ectopic vulval precursor cell proliferation (the multivulva [Muv] phenotype) due to inappropriate expression of the lin-3/EGF gene, and RNAseq analysis suggests a broad role for these ivp genes in transcriptional repression. Importantly, although the C. briggsae phenotypes have parallels with those seen in the C. elegans synMuv system, except for the highly conserved HTZ-1/H2A.z, comparable mutations in C. elegans ivp orthologs do not exhibit synMuv gene interactions or phenotypes. These results demonstrate the evolutionary changes that can underlie conserved biological outputs and argue that proteins critical to repress inappropriate expression from the genome participate in a rapidly evolving functional landscape.

Genome-Wide Natural Selection Signatures Are Linked to Genetic Risk of Modern Phenotypes in the Japanese Population

Mon, 20 Jan 2020 00:00:00 GMT

Elucidation of natural selection signatures and relationships with phenotype spectra is important to understand adaptive evolution of modern humans. Here, we conducted a genome-wide scan of selection signatures of the Japanese population by estimating locus-specific time to the most recent common ancestor using the ascertained sequentially Markovian coalescent (ASMC), from the biobank-based large-scale genome-wide association study data of 170,882 subjects. We identified 29 genetic loci with selection signatures satisfying the genome-wide significance. The signatures were most evident at the alcohol dehydrogenase (ADH) gene cluster locus at 4q23 (P_ASMC = 2.2 × 10⁻³⁶), followed by relatively strong selection at the FAM96A (15q22), MYOF (10q23), 13q21, GRIA2 (4q32), and ASAP2 (2p25) loci (P_ASMC < 1.0 × 10⁻¹⁰). The additional analysis interrogating extended haplotypes (integrated haplotype score) showed robust concordance of the detected signatures, contributing to fine-mapping of the genes, and provided allelic directional insights into selection pressure (e.g., positive selection for ADH1B-Arg48His and HLA-DPB1*04:01). The phenome-wide selection enrichment analysis with the trait-associated variants identified a variety of the modern human phenotypes involved in the adaptation of Japanese. We observed population-specific evidence of enrichment with the alcohol-related phenotypes, anthropometric and biochemical clinical measurements, and immune-related diseases, differently from the findings in Europeans using the UK Biobank resource. Our study demonstrated population-specific features of the selection signatures in Japanese, highlighting a value of the natural selection study using the nation-wide biobank-scale genome and phenotype data.

Polygenic Patterns of Adaptive Introgression in Modern Humans Are Mainly Shaped by Response to Pathogens

Tue, 14 Jan 2020 00:00:00 GMT

Anatomically modern humans carry many introgressed variants from other hominins in their genomes. Some of them affect their phenotype and can thus be negatively or positively selected. Several individual genes have been proposed to be the subject of adaptive introgression, but the possibility of polygenic adaptive introgression has not been extensively investigated yet. In this study, we analyze archaic introgression maps with refined functional enrichment methods to find signals of polygenic adaptation of introgressed variants. We first apply a method to detect sets of connected genes (subnetworks) within biological pathways that present higher-than-expected levels of archaic introgression. We then introduce and apply a new statistical test to distinguish between epistatic and independent selection in gene sets of present-day humans. We identify several known targets of adaptive introgression, and we show that they belong to larger networks of introgressed genes. After correction for genetic linkage, we find that signals of polygenic adaptation are mostly explained by independent and potentially sequential selection episodes. However, we also find some gene sets where introgressed variants present significant signals of epistatic selection. Our results confirm that archaic introgression has facilitated local adaptation, especially in immunity related and metabolic functions and highlight its involvement in a coordinated response to pathogens out of Africa.

Disentangling Population History and Character Evolution among Hybridizing Lineages

Mon, 13 Jan 2020 00:00:00 GMT

Understanding the origin and maintenance of adaptive phenotypic novelty is a central goal of evolutionary biology. However, both hybridization and incomplete lineage sorting can lead to genealogical discordance between the regions of the genome underlying adaptive traits and the remainder of the genome, decoupling inferences about character evolution from population history. Here, to disentangle these effects, we investigated the evolutionary origins and maintenance of Batesian mimicry between North American admiral butterflies (Limenitis arthemis) and their chemically defended model (Battus philenor) using a combination of de novo genome sequencing, whole-genome resequencing, and statistical introgression mapping. Our results suggest that balancing selection, arising from geographic variation in the presence or absence of the unpalatable model, has maintained two deeply divergent color patterning haplotypes that have been repeatedly sieved among distinct mimetic and nonmimetic lineages of Limenitis via introgressive hybridization.

Genetic Basis of De Novo Appearance of Carotenoid Ornamentation in Bare Parts of Canaries

Mon, 13 Jan 2020 00:00:00 GMT

Unlike wild and domestic canaries (Serinus canaria), or any of the three dozen species of finches in genus Serinus, the domestic urucum breed of canaries exhibits bright red bills and legs. This novel trait offers a unique opportunity to understand the mechanisms of bare-part coloration in birds. To identify the mutation producing the colorful phenotype, we resequenced the genome of urucum canaries and performed a range of analyses to search for genotype-to-phenotype associations across the genome. We identified a nonsynonymous mutation in the gene BCO2 (beta-carotene oxygenase 2, also known as BCDO2), an enzyme involved in the cleavage and breakdown of full-length carotenoids into short apocarotenoids. Protein structural models and in vitro functional assays indicate that the urucum mutation abrogates the carotenoid-cleavage activity of BCO2. Consistent with the predicted loss of carotenoid-cleavage activity, urucum canaries tended to have increased levels of full-length carotenoid pigments in bill tissue and reduced levels of carotenoid-cleavage products (apocarotenoids) in retinal tissue compared with other breeds of canaries. We hypothesize that carotenoid-based bare-part coloration might be readily gained, modified, or lost through simple switches in the enzymatic activity or regulation of BCO2 and this gene may be an important mediator in the evolution of bare-part coloration among bird species.

Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9

Sat, 11 Jan 2020 00:00:00 GMT

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions—a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions—a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.

Recent Out-of-Africa Migration of Human Herpes Simplex Viruses

Thu, 09 Jan 2020 00:00:00 GMT

Herpes simplex virus types 1 and 2 (HSV-1 and HSV-2) are ubiquitous human pathogens. Both viruses evolved from simplex viruses infecting African primates and they are thus thought to have left Africa during early human migrations. We analyzed the population structure of HSV-1 and HSV-2 circulating strains. Results indicated that HSV-1 populations have limited geographic structure and the most evident clustering by geography is likely due to recent bottlenecks. For HSV-2, the only level of population structure is accounted for by the so-called “worldwide” and “African” lineages. Analysis of ancestry components and nucleotide diversity, however, did not support the view that the worldwide lineage followed early humans during out-of-Africa dispersal. Although phylogeographic analysis confirmed an African origin for both viruses, molecular dating with a method that corrects for the time-dependent rate phenomenon indicated that HSV-1 and HSV-2 migrated from Africa in relatively recent times. In particular, we estimated that the HSV-2 worldwide lineage left the continent in the 18th century, which corresponds to the height of the transatlantic slave trade, possibly explaining the high prevalence of HSV-2 in the Americas (second highest after Africa). The limited geographic clustering of HSV-1 makes it difficult to date its exit from Africa. The split between the basal clade, containing mostly African sequences, and all other strains was dated at ∼5,000 years ago. Our data do not imply that herpes simplex viruses did not infect early humans but show that the worldwide distribution of circulating strains is the result of relatively recent events.

Ancient DNA Evidence from China Reveals the Expansion of Pacific Dogs

Wed, 08 Jan 2020 00:00:00 GMT

The ancestral homeland of Australian dingoes and Pacific dogs is proposed to be in South China. However, the location and timing of their dispersal and relationship to dog domestication is unclear. Here, we sequenced 7,000- to 2,000-year-old complete mitochondrial DNA (mtDNA) genomes of 27 ancient canids (one gray wolf and 26 domestic dogs) from the Yellow River and Yangtze River basins (YYRB). These are the first complete ancient mtDNA of Chinese dogs from the cradle of early Chinese civilization. We found that most ancient dogs (18/26) belong to the haplogroup A1b lineage that is found in high frequency in present-day Australian dingoes and precolonial Pacific Island dogs but low frequency in present-day China. Particularly, a 7,000-year-old dog from the Tianluoshan site in Zhejiang province possesses a haplotype basal to the entire haplogroup A1b lineage. We propose that A1b lineage dogs were once widely distributed in the YYRB area. Following their dispersal to South China, and then into Southeast Asia, New Guinea and remote Oceania, they were largely replaced by dogs belonging to other lineages in the last 2,000 years in present-day China, especially North China.

Reduced Translational Efficiency of Eukaryotic Genes after Duplication Events

Mon, 06 Jan 2020 00:00:00 GMT

Control of gene expression has been found to be predominantly determined at the level of protein translation. However, to date, reduced expression from duplicated genes in eukaryotes for dosage maintenance has only been linked to transcriptional control involving epigenetic mechanisms. Here, we hypothesize that dosage maintenance following gene duplication also involves regulation at the protein level. To test this hypothesis, we compared transcriptome and proteome data of yeast models, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and worm models, Caenorhabditis elegans and Caenorhabditis briggsae, to investigate lineage-specifically duplicated genes. Duplicated genes in both eukaryotic models exhibited a reduced protein-to-mRNA abundance ratio. Moreover, dosage sensitive genes, represented by genes encoding protein complex subunits, reduced their protein-to-mRNA abundance ratios more significantly than the other genes after duplication events. An analysis of ribosome profiling (Ribo-Seq) data further showed that reduced translational efficiency was more prominent for dosage sensitive genes than for the other genes. Meanwhile, no difference in protein degradation rate was associated with duplication events. Translationally repressed duplicated genes were also more likely to be inhibited at the level of transcription. Taken together, these results suggest that translation-mediated dosage control is partially contributed by natural selection and it enhances transcriptional control in maintaining gene dosage after gene duplication events during eukaryotic genome evolution.

Origin and Evolution of Carboxysome Positioning Systems in Cyanobacteria

Fri, 03 Jan 2020 00:00:00 GMT

Carboxysomes are protein-based organelles that are essential for allowing cyanobacteria to fix CO₂. Previously, we identified a two-component system, McdAB, responsible for equidistantly positioning carboxysomes in the model cyanobacterium Synechococcus elongatus PCC 7942 (MacCready JS, Hakim P, Young EJ, Hu L, Liu J, Osteryoung KW, Vecchiarelli AG, Ducat DC. 2018. Protein gradients on the nucleoid position the carbon-fixing organelles of cyanobacteria. eLife 7:pii:e39723). McdA, a ParA-type ATPase, nonspecifically binds the nucleoid in the presence of ATP. McdB, a novel factor that directly binds carboxysomes, displaces McdA from the nucleoid. Removal of McdA from the nucleoid in the vicinity of carboxysomes by McdB causes a global break in McdA symmetry, and carboxysome motion occurs via a Brownian-ratchet-based mechanism toward the highest concentration of McdA. Despite the importance for cyanobacteria to properly position their carboxysomes, whether the McdAB system is widespread among cyanobacteria remains an open question. Here, we show that the McdAB system is widespread among β-cyanobacteria, often clustering with carboxysome-related components, and is absent in α-cyanobacteria. Moreover, we show that two distinct McdAB systems exist in β-cyanobacteria, with Type 2 systems being the most ancestral and abundant, and Type 1 systems, like that of S. elongatus, possibly being acquired more recently. Lastly, all McdB proteins share the sequence signatures of a protein capable of undergoing liquid–liquid phase separation. Indeed, we find that representatives of both McdB types undergo liquid–liquid phase separation in vitro, the first example of a ParA-type ATPase partner protein to exhibit this behavior. Our results have broader implications for understanding carboxysome evolution, biogenesis, homeostasis, and positioning in cyanobacteria.

Insights from Empirical Analyses and Simulations on Using Multiple Fossil Calibrations with Relaxed Clocks to Estimate Divergence Times

Fri, 03 Jan 2020 00:00:00 GMT

Relaxed clock methods account for among-branch-rate-variation when estimating divergence times by inferring different rates for individual branches. In order to infer different rates for individual branches, important assumptions are required. This is because molecular sequence data do not provide direct information about rates but instead provide direct information about the total number of substitutions along any branch, which is a product of the rate and time for that branch. Often, the assumptions required for estimating rates for individual branches depend heavily on the implementation of multiple fossil calibrations in a single phylogeny. Here, we show that the basis of these assumptions is often critically undermined. First, we highlight that the temporal distribution of the fossil record often violates key assumptions of methods that use multiple fossil calibrations with relaxed clocks. With respect to “node calibration” methods, this conclusion is based on our inference that different fossil calibrations are unlikely to reflect the relative ages of different clades. With respect to the fossilized birth–death process, this conclusion is based on our inference that the fossil recovery rate is often highly heterogeneous. We then demonstrate that methods of divergence time estimation that use multiple fossil calibrations are highly sensitive to assumptions about the fossil record and among-branch-rate-variation. Given the problems associated with these assumptions, our results highlight that using multiple fossil calibrations with relaxed clocks often does little to improve the accuracy of divergence time estimates.

Deep Residual Neural Networks Resolve Quartet Molecular Phylogenies

Mon, 23 Dec 2019 00:00:00 GMT

Phylogenetic inference is of fundamental importance to evolutionary as well as other fields of biology, and molecular sequences have emerged as the primary data for this task. Although many phylogenetic methods have been developed to explicitly take into account substitution models of sequence evolution, such methods could fail due to model misspecification or insufficiency, especially in the face of heterogeneities in substitution processes across sites and among lineages. In this study, we propose to infer topologies of four-taxon trees using deep residual neural networks, a machine learning approach needing no explicit modeling of the subject system and having a record of success in solving complex nonlinear inference problems. We train residual networks on simulated protein sequence data with extensive amino acid substitution heterogeneities. We show that the well-trained residual network predictors can outperform existing state-of-the-art inference methods such as the maximum likelihood method on diverse simulated test data, especially under extensive substitution heterogeneities. Reassuringly, residual network predictors generally agree with existing methods in the trees inferred from real phylogenetic data with known or widely believed topologies. Furthermore, when combined with the quartet puzzling algorithm, residual network predictors can be used to reconstruct trees with more than four taxa. We conclude that deep learning represents a powerful new approach to phylogenetic reconstruction, especially when sequences evolve via heterogeneous substitution processes. We present our best trained predictor in a freely available program named Phylogenetics by Deep Learning (PhyDL, https://gitlab.com/ztzou/phydl; last accessed January 3, 2020).

Horizontal Gene Transfer Phylogenetics: A Random Walk Approach

Mon, 23 Dec 2019 00:00:00 GMT

The dramatic decrease in time and cost for generating genetic sequence data has opened up vast opportunities in molecular systematics, one of which is the ability to decipher the evolutionary history of strains of a species. Under this fine systematic resolution, the standard markers are too crude to provide a phylogenetic signal. Nevertheless, among prokaryotes, genome dynamics in the form of horizontal gene transfer (HGT) between organisms and gene loss seem to provide far richer information by affecting both gene order and gene content. The “synteny index” (SI) between a pair of genomes combines these latter two factors, allowing comparison of genomes with unequal gene content, together with order considerations of their common genes. Although this approach is useful for classifying close relatives, no rigorous statistical modeling for it has been suggested. Such modeling is valuable, as it allows observed measures to be transformed into estimates of time periods during evolution, yielding the “additivity” of the measure. To the best of our knowledge, there is no other additivity proof for other gene order/content measures under HGT. Here, we provide a first statistical model and analysis for the SI measure. We model the “gene neighborhood” as a “birth–death–immigration” process affected by the HGT activity over the genome, and analytically relate the HGT rate and time to the expected SI. This model is asymptotic and thus provides accurate results, assuming infinite size genomes. Therefore, we also developed a heuristic model following an “exponential decay” function, accounting for biologically realistic values, which performed well in simulations. Applying this model to 1,133 prokaryotes partitioned to 39 clusters by the rank of genus yields that the average number of genome dynamics events per gene in the phylogenetic depth of genus is around half with significant variability between genera. This result extends and confirms similar results obtained for individual genera in different manners.

Parallel Seed Color Adaptation during Multiple Domestication Attempts of an Ancient New World Grain

Fri, 20 Dec 2019 00:00:00 GMT

Thousands of plants have been selected as crops; yet, only a few are fully domesticated. The lack of adaptation to agroecological environments of many crop plants with few characteristic domestication traits potentially has genetic causes. Here, we investigate the incomplete domestication of an ancient grain from the Americas, amaranth. Although three grain amaranth species have been cultivated as crop for millennia, all three lack key domestication traits. We sequenced 121 crop and wild individuals to investigate the genomic signature of repeated incomplete adaptation. Our analysis shows that grain amaranth has been domesticated three times from a single wild ancestor. One trait that has been selected during domestication in all three grain species is the seed color, which changed from dark seeds to white seeds. We were able to map the genetic control of the seed color adaptation to two genomic regions on chromosomes 3 and 9, employing three independent mapping populations. Within the locus on chromosome 9, we identify an MYB-like transcription factor gene, a known regulator for seed color variation in other plant species. We identify a soft selective sweep in this genomic region in one of the crop species but not in the other two species. The demographic analysis of wild and domesticated amaranths revealed a population bottleneck predating the domestication of grain amaranth. Our results indicate that a reduced level of ancestral genetic variation did not prevent the selection of traits with a simple genetic architecture but may have limited the adaptation of complex domestication traits.

Probabilities of Unranked and Ranked Anomaly Zones under Birth–Death Models

Fri, 20 Dec 2019 00:00:00 GMT

A labeled gene tree topology that is more probable than the labeled gene tree topology matching a species tree is called “anomalous.” Species trees that can generate such anomalous gene trees are said to be in the “anomaly zone.” Here, probabilities of “unranked” and “ranked” gene tree topologies under the multispecies coalescent are considered. A ranked tree depicts not only the topological relationship among gene lineages, as an unranked tree does, but also the sequence in which the lineages coalesce. In this article, we study how the parameters of a species tree simulated under a constant-rate birth–death process can affect the probability that the species tree lies in the anomaly zone. We find that with more than five taxa, it is possible for species trees to have both anomalous unranked and ranked gene trees. The probability of being in either type of anomaly zone increases with more taxa. The probability of anomalous gene trees also increases with higher speciation rates. We observe that the probabilities of unranked anomaly zones are higher and grow much faster than those of ranked anomaly zones as the speciation rate increases. Our simulation shows that the most probable ranked gene tree is likely to have the same unranked topology as the species tree. We design the software PRANC, which computes probabilities of ranked gene tree topologies given a species tree under the coalescent model.

Using Selection by Nonantibiotic Stressors to Sensitize Bacteria to Antibiotics

Wed, 18 Dec 2019 00:00:00 GMT

Evolutionary adaptation of bacteria to nonantibiotic selective forces, such as osmotic stress, has been previously associated with increased antibiotic resistance, but much less is known about potentially sensitizing effects of nonantibiotic stressors. In this study, we use laboratory evolution to investigate adaptation of Enterococcus faecalis, an opportunistic bacterial pathogen, to a broad collection of environmental agents, ranging from antibiotics and biocides to extreme pH and osmotic stress. We find that nonantibiotic selection frequently leads to increased sensitivity to other conditions, including multiple antibiotics. Using population sequencing and whole-genome sequencing of single isolates from the evolved populations, we identify multiple mutations in genes previously linked with resistance to the selecting conditions, including genes corresponding to known drug targets or multidrug efflux systems previously tied to collateral sensitivity. Finally, we hypothesized based on the measured sensitivity profiles that sequential rounds of antibiotic and nonantibiotic selection may lead to hypersensitive populations by harnessing the orthogonal collateral effects of particular pairs of selective forces. To test this hypothesis, we show experimentally that populations evolved to a sequence of linezolid (an oxazolidinone antibiotic) and sodium benzoate (a common preservative) exhibit increased sensitivity to more stressors than adaptation to either condition alone. The results demonstrate how sequential adaptation to drug and nondrug environments can be used to sensitize bacteria to antibiotics and highlight new potential strategies for exploiting shared constraints governing adaptation to diverse environmental challenges.