Of clades and clans: terms for phylogenetic relationships in unrooted trees

The Complex GH32 Enzyme Orchestra from Priestia megaterium Holds the Key to Better Discriminate Sucrose-6-phosphate Hydrolases from Other β-Fructofuranosidases in Bacteria

Inulin is widely used as a prebiotic and emerging as a priming compound to counteract plant diseases. We isolated inulin-degrading strains from the lettuce phyllosphere, identified as Bacillus subtilis and Priestia megaterium, species hosting well-known biocontrol organisms. To better understand their varying inulin degradation strategies, three intracellular β-fructofuranosidases from P. megaterium NBRC15308 were characterized after expression in Escherichia coli: a predicted sucrose-6-phosphate (Suc6P) hydrolase (SacAP1, supported by molecular docking), an exofructanase (SacAP2), and an invertase (SacAP3). Based on protein multiple sequence and structure alignments of bacterial glycoside hydrolase family 32 enzymes, we identified conserved residues predicted to be involved in binding phosphorylated (Suc6P hydrolases) or nonphosphorylated substrates (invertases and fructanases). Suc6P hydrolases feature positively charged residues near the structural catalytic pocket (histidine, arginine, or lysine), whereas other β-fructofuranosidases contain tryptophans. This correlates with our phylogenetic tree, grouping all predicted Suc6P hydrolases in a clan associated with genomic regions coding for transporters involved in substrate phosphorylation. These results will help to discriminate between Suc6P hydrolases and other β-fructofuranosidases in future studies and to better understand the interaction of B. subtilis and P. megaterium endophytes with sucrose and/or fructans, sugars naturally present in plants or exogenously applied in the context of defense priming.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

Improving Orthologous Signal and Model Fit in Datasets Addressing the Root of the Animal Phylogeny

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference.

Filtering artifactual signal increases support for Xenacoelomorpha and Ambulacraria sister relationship in the animal tree of life

Conflicting studies place a group of bilaterian invertebrates containing xenoturbellids and acoelomorphs, the Xenacoelomorpha, as either the primary emerging bilaterian phylum^1,2,3,4,5,6 or within Deuterostomia, sister to Ambulacraria.^7,8,9,10,11 Although their placement as sister to the rest of Bilateria supports relatively simple morphology in the ancestral bilaterian, their alternative placement within Deuterostomia suggests a morphologically complex ancestral bilaterian along with extensive loss of major phenotypic traits in the Xenacoelomorpha. Recent studies have questioned whether Deuterostomia should be considered monophyletic at all.^10,12,13 Hidden paralogy and poor phylogenetic signal present a major challenge for reconstructing species phylogenies.^{14,15,16,17,18} Here, we assess whether these issues have contributed to the conflict over the placement of Xenacoelomorpha. We reanalyzed published datasets, enriching for orthogroups whose gene trees support well-resolved clans elsewhere in the animal tree.¹⁶ We find that most genes in previously published datasets violate incontestable clans, suggesting that hidden paralogy and low phylogenetic signal affect the ability to reconstruct branching patterns at deep nodes in the animal tree. We demonstrate that removing orthogroups that cannot recapitulate incontestable relationships alters the final topology that is inferred, while simultaneously improving the fit of the model to the data. We discover increased, but ultimately not conclusive, support for the existence of Xenambulacraria in our set of filtered orthogroups. At a time when we are progressing toward sequencing all life on the planet, we argue that long-standing contentious issues in the tree of life will be resolved using smaller amounts of better quality data that can be modeled adequately.¹⁹.

HUG Domain Is Responsible for Active Dimer Stabilization in an NrdJd Ribonucleotide Reductase

Ribonucleotide reductases (RNRs) catalyze the reduction of ribonucleotides to the corresponding deoxyribonucleotides. The catalytic activity of most RNRs depends on the formation of a dimer of the catalytic subunits. The active site is located at the interface, and part of the substrate binding site and regulatory mechanisms work across the subunit in the dimer. In this study, we describe and characterize a novel domain responsible for forming the catalytic dimer in several class II RNRs. The 3D structure of the class II RNR from Rhodobacter sphaeroides reveals a so far undescribed α-helical domain in the dimer interface, which is embracing the other subunit. Genetic removal of this HUG domain leads to a severe reduction of activity paired with reduced dimerization capability. In comparison with other described RNRs, the enzyme with this domain is less dependent on the presence of nucleotides to act as allosteric effectors in the formation of dimers. The HUG domain appears to serve as an interlock to keep the dimer intact and functional even at low enzyme and/or effector concentrations.

SpeciesRax: A tool for maximum likelihood species tree inference from gene family trees under duplication, transfer, and loss

Species tree inference from gene family trees is becoming increasingly popular because it can account for discordance between the species tree and the corresponding gene family trees. In particular, methods that can account for multiple-copy gene families exhibit potential to leverage paralogy as informative signal. At present, there does not exist any widely adopted inference method for this purpose. Here, we present SpeciesRax, the first maximum likelihood method that can infer a rooted species tree from a set of gene family trees and can account for gene duplication, loss, and transfer events. By explicitly modeling events by which gene trees can depart from the species tree, SpeciesRax leverages the phylogenetic rooting signal in gene trees. SpeciesRax infers species tree branch lengths in units of expected substitutions per site and branch support values via paralogy-aware quartets extracted from the gene family trees. Using both empirical and simulated data sets we show that SpeciesRax is at least as accurate as the best competing methods while being one order of magnitude faster on large data sets at the same time. We used SpeciesRax to infer a biologically plausible rooted phylogeny of the vertebrates comprising 188 species from 31,612 gene families in 1 h using 40 cores. SpeciesRax is available under GNU GPL at https://github.com/BenoitMorel/GeneRax and on BioConda.

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.

A rooted phylogeny resolves early bacterial evolution

A rooted bacterial tree is necessary to understand early evolution, but the position of the root is contested. Here, we model the evolution of 11,272 gene families to identify the root, extent of horizontal gene transfer (HGT), and the nature of the last bacterial common ancestor (LBCA). Our analyses root the tree between the major clades Terrabacteria and Gracilicutes and suggest that LBCA was a free-living flagellated, rod-shaped double-membraned organism. Contrary to recent proposals, our analyses reject a basal placement of the Candidate Phyla Radiation, which instead branches sister to Chloroflexota within Terrabacteria. While most gene families (92%) have evidence of HGT, overall, two-thirds of gene transmissions have been vertical, suggesting that a rooted tree provides a meaningful frame of reference for interpreting bacterial evolution.

Inferring the deep past from molecular data

There is an expectation that analyses of molecular sequences might be able to distinguish between alternative hypotheses for ancient relationships, but the phylogenetic methods used and types of data analyzed are of critical importance in any attempt to recover historical signal. Here, we discuss some common issues that can influence the topology of trees obtained when using overly simple models to analyze molecular data that often display complicated patterns of sequence heterogeneity. To illustrate our discussion, we have used three examples of inferred relationships which have changed radically as models and methods of analysis have improved. In two of these examples, the sister-group relationship between thermophilic Thermus and mesophilic Deinococcus, and the position of long-branch Microsporidia among eukaryotes, we show that recovering what is now generally considered to be the correct tree is critically dependent on the fit between model and data. In the third example, the position of eukaryotes in the tree of life, the hypothesis that is currently supported by the best available methods is fundamentally different from the classical view of relationships between major cellular domains. Since heterogeneity appears to be pervasive and varied among all molecular sequence data, and even the best available models can still struggle to deal with some problems, the issues we discuss are generally relevant to phylogenetic analyses. It remains essential to maintain a critical attitude to all trees as hypotheses of relationship that may change with more data and better methods.

Atypical Neurogenesis in Induced Pluripotent Stem Cells From Autistic Individuals

BACKGROUND

Autism is a heterogeneous collection of disorders with a complex molecular underpinning. Evidence from postmortem brain studies have indicated that early prenatal development may be altered in autism. Induced pluripotent stem cells (iPSCs) generated from individuals with autism with macrocephaly also indicate prenatal development as a critical period for this condition. But little is known about early altered cellular events during prenatal stages in autism.

METHODS

iPSCs were generated from 9 unrelated individuals with autism without macrocephaly and with heterogeneous genetic backgrounds, and 6 typically developing control individuals. iPSCs were differentiated toward either cortical or midbrain fates. Gene expression and high throughput cellular phenotyping was used to characterize iPSCs at different stages of differentiation.

RESULTS

A subset of autism-iPSC cortical neurons were RNA-sequenced to reveal autism-specific signatures similar to postmortem brain studies, indicating a potential common biological mechanism. Autism-iPSCs differentiated toward a cortical fate displayed impairments in the ability to self-form into neural rosettes. In addition, autism-iPSCs demonstrated significant differences in rate of cell type assignment of cortical precursors and dorsal and ventral forebrain precursors. These cellular phenotypes occurred in the absence of alterations in cell proliferation during cortical differentiation, differing from previous studies. Acquisition of cell fate during midbrain differentiation was not different between control- and autism-iPSCs.

CONCLUSIONS

Taken together, our data indicate that autism-iPSCs diverge from control-iPSCs at a cellular level during early stage of neurodevelopment. This suggests that unique developmental differences associated with autism may be established at early prenatal stages.

Spectral neighbor joining for reconstruction of latent tree Models

A common assumption in multiple scientific applications is that the distribution of observed data can be modeled by a latent tree graphical model. An important example is phylogenetics, where the tree models the evolutionary lineages of a set of observed organisms. Given a set of independent realizations of the random variables at the leaves of the tree, a key challenge is to infer the underlying tree topology. In this work we develop Spectral Neighbor Joining (SNJ), a novel method to recover the structure of latent tree graphical models. Given a matrix that contains a measure of similarity between all pairs of observed variables, SNJ computes a spectral measure of cohesion between groups of observed variables. We prove that SNJ is consistent, and derive a sufficient condition for correct tree recovery from an estimated similarity matrix. Combining this condition with a concentration of measure result on the similarity matrix, we bound the number of samples required to recover the tree with high probability. We illustrate via extensive simulations that in comparison to several other reconstruction methods, SNJ requires fewer samples to accurately recover trees with a large number of leaves or long edges.

The Origin and Evolution of Antistasin-like Proteins in Leeches (Hirudinida, Clitellata)

Bloodfeeding is employed by many parasitic animals and requires specific innovations for efficient feeding. Some of these innovations are molecular features that are related to the inhibition of hemostasis. For example, bloodfeeding insects, bats, and leeches release proteins with anticoagulatory activity through their salivary secretions. The antistasin-like protein family, composed of serine protease inhibitors with one or more antistasin-like domains, is tightly linked to inhibition of hemostasis in leeches. However, this protein family has been recorded also in non-bloodfeeding invertebrates, such as cnidarians, mollusks, polychaetes, and oligochaetes. The present study aims to 1) root the antistasin-like gene tree and delimit the major orthologous groups, 2) identify potential independent origins of salivary proteins secreted by leeches, and 3) identify major changes in domain and/or motif structure within each orthologous group. Five clades containing leech antistasin-like proteins are distinguishable through rigorous phylogenetic analyses based on nine new transcriptomes and a diverse set of comparative data: the trypsin + leukocyte elastase inhibitors clade, the antistasin clade, the therostasin clade, and two additional, unnamed clades. The antistasin-like gene tree supports multiple origins of leech antistasin-like proteins due to the presence of both leech and non-leech sequences in one of the unnamed clades, but a single origin of factor Xa and trypsin + leukocyte elastase inhibitors. This is further supported by three sequence motifs that are exclusive to antistasins, the trypsin + leukocyte elastase inhibitor clade, and the therostasin clade, respectively. We discuss the implications of our findings for the evolution of this diverse family of leech anticoagulants.

Redefining the H-NS protein family: a diversity of specialized core and accessory forms exhibit hierarchical transcriptional network integration

H-NS is a nucleoid structuring protein and global repressor of virulence and horizontally-acquired genes in bacteria. H-NS can interact with itself or with homologous proteins, but protein family diversity and regulatory network overlap remain poorly defined. Here, we present a comprehensive phylogenetic analysis that revealed deep-branching clades, dispelling the presumption that H-NS is the progenitor of varied molecular backups. Each clade is composed exclusively of either chromosome-encoded or plasmid-encoded proteins. On chromosomes, stpA and newly discovered hlpP are core genes in specific genera, whereas hfp and newly discovered hlpC are sporadically distributed. Six clades of H-NS plasmid proteins (Hpp) exhibit ancient and dedicated associations with plasmids, including three clades with fidelity for plasmid incompatibility groups H, F or X. A proliferation of H-NS homologs in Erwiniaceae includes the first observation of potentially co-dependent H-NS forms. Conversely, the observed diversification of oligomerization domains may facilitate stable co-existence of divergent homologs in a genome. Transcriptomic and proteomic analysis in Salmonella revealed regulatory crosstalk and hierarchical control of H-NS homologs. We also discovered that H-NS is both a repressor and activator of Salmonella Pathogenicity Island 1 gene expression, and both regulatory modes are restored by Sfh (HppH) in the absence of H-NS.

Undinarchaeota illuminate DPANN phylogeny and the impact of gene transfer on archaeal evolution

The recently discovered DPANN archaea are a potentially deep-branching, monophyletic radiation of organisms with small cells and genomes. However, the monophyly and early emergence of the various DPANN clades and their role in life's evolution are debated. Here, we reconstructed and analysed genomes of an uncharacterized archaeal phylum (Candidatus Undinarchaeota), revealing that its members have small genomes and, while potentially being able to conserve energy through fermentation, likely depend on partner organisms for the acquisition of certain metabolites. Our phylogenomic analyses robustly place Undinarchaeota as an independent lineage between two highly supported 'DPANN' clans. Further, our analyses suggest that DPANN have exchanged core genes with their hosts, adding to the difficulty of placing DPANN in the tree of life. This pattern can be sufficiently dominant to allow identifying known symbiont-host clades based on routes of gene transfer. Together, our work provides insights into the origins and evolution of DPANN and their hosts.

Draft genome of the European medicinal leech Hirudo medicinalis (Annelida, Clitellata, Hirudiniformes) with emphasis on anticoagulants

The European medicinal leech has been used for medicinal purposes for millennia, and continues to be used today in modern hospital settings. Its utility is granted by the extremely potent anticoagulation factors that the leech secretes into the incision wound during feeding and, although a handful of studies have targeted certain anticoagulants, the full range of anticoagulation factors expressed by this species remains unknown. Here, we present the first draft genome of the European medicinal leech, Hirudo medicinalis, and estimate that we have sequenced between 79–94% of the full genome. Leveraging these data, we searched for anticoagulation factors across the genome of H. medicinalis. Following orthology determination through a series of BLAST searches, as well as phylogenetic analyses, we estimate that fully 15 different known anticoagulation factors are utilized by the species, and that 17 other proteins that have been linked to antihemostasis are also present in the genome. We underscore the utility of the draft genome for comparative studies of leeches and discuss our results in an evolutionary context.

Inadvertent Paralog Inclusion Drives Artifactual Topologies and Timetree Estimates in Phylogenomics

Increasingly, large phylogenomic data sets include transcriptomic data from nonmodel organisms. This not only has allowed controversial and unexplored evolutionary relationships in the tree of life to be addressed but also increases the risk of inadvertent inclusion of paralogs in the analysis. Although this may be expected to result in decreased phylogenetic support, it is not clear if it could also drive highly supported artifactual relationships. Many groups, including the hyperdiverse Lissamphibia, are especially susceptible to these issues due to ancient gene duplication events and small numbers of sequenced genomes and because transcriptomes are increasingly applied to resolve historically conflicting taxonomic hypotheses. We tested the potential impact of paralog inclusion on the topologies and timetree estimates of the Lissamphibia using published and de novo sequencing data including 18 amphibian species, from which 2,656 single-copy gene families were identified. A novel paralog filtering approach resulted in four differently curated data sets, which were used for phylogenetic reconstructions using Bayesian inference, maximum likelihood, and quartet-based supertrees. We found that paralogs drive strongly supported conflicting hypotheses within the Lissamphibia (Batrachia and Procera) and older divergence time estimates even within groups where no variation in topology was observed. All investigated methods, except Bayesian inference with the CAT-GTR model, were found to be sensitive to paralogs, but with filtering convergence to the same answer (Batrachia) was observed. This is the first large-scale study to address the impact of orthology selection using transcriptomic data and emphasizes the importance of quality over quantity particularly for understanding relationships of poorly sampled taxa.

The salivary transcriptome of Limnobdella mexicana (Annelida: Clitellata: Praobdellidae) and orthology determination of major leech anticoagulants

Bloodfeeding requires several adaptations that allow the parasite to feed efficiently. Leeches and other hematophagous animals have developed different mechanisms to inhibit hemostasis, one of the main barriers imposed by their hosts. Limnobdella mexicana is a member of the leech family Praobdellidae, a family of host generalists known for their preference to attach on mucosal membranes of mammals, such as those in nasopharyngeal cavities, bladders and ocular orbits. Previous studies have hypothesized a positive relationship between diversity of anticoagulants and diversity of hosts in bloodfeeding leeches. However, orthology determination of putative anticoagulants and the lack of standardization of sequencing effort and method hinder comparisons between publicly available transcriptomes generated in different laboratories. In the present study, we examine the first transcriptome of a praobdellid leech and identify 15 putative anticoagulants using a phylogeny-based inference approach, amino-acid conservation, Pfam domains and BLAST searches. Our phylogenetic analyses suggest that the ancestral leech was able to inhibit factor Xa and that some hirudins that have been reported in previous studies on leech anticoagulants may not be orthologous with the archetypal hirudin.

Investigating the Origins of Membrane Phospholipid Biosynthesis Genes Using Outgroup-Free Rooting

One of the key differences between Bacteria and Archaea is their canonical membrane phospholipids, which are synthesized by distinct biosynthetic pathways with nonhomologous enzymes. This “lipid divide” has important implications for the early evolution of cells and the type of membrane phospholipids present in the last universal common ancestor. One of the main challenges in studies of membrane evolution is that the key biosynthetic genes are ancient and their evolutionary histories are poorly resolved. This poses major challenges for traditional rooting methods because the only available outgroups are distantly related. Here, we address this issue by using the best available substitution models for single-gene trees, by expanding our analyses to the diversity of uncultivated prokaryotes recently revealed by environmental genomics, and by using two complementary approaches to rooting that do not depend on outgroups. Consistent with some previous analyses, our rooted gene trees support extensive interdomain horizontal transfer of membrane phospholipid biosynthetic genes, primarily from Archaea to Bacteria. They also suggest that the capacity to make archaeal-type membrane phospholipids was already present in last universal common ancestor.

Identification and characterization of putative Aeromonas spp. T3SS effectors

The genetic determinants of bacterial pathogenicity are highly variable between species and strains. However, a factor that is commonly associated with virulent Gram-negative bacteria, including many Aeromonas spp., is the type 3 secretion system (T3SS), which is used to inject effector proteins into target eukaryotic cells. In this study, we developed a bioinformatics pipeline to identify T3SS effector proteins, applied this approach to the genomes of 105 Aeromonas strains isolated from environmental, mutualistic, or pathogenic contexts and evaluated the cytotoxicity of the identified effectors through their heterologous expression in yeast. The developed pipeline uses a two-step approach, where candidate Aeromonas gene families are initially selected using Hidden Markov Model (HMM) profile searches against the Virulence Factors DataBase (VFDB), followed by strict comparisons against positive and negative control datasets, greatly reducing the number of false positives. This approach identified 21 Aeromonas T3SS likely effector families, of which 8 represent known or characterized effectors, while the remaining 13 have not previously been described in Aeromonas. We experimentally validated our in silico findings by assessing the cytotoxicity of representative effectors in Saccharomyces cerevisiae BY4741, with 15 out of 21 assayed proteins eliciting a cytotoxic effect in yeast. The results of this study demonstrate the utility of our approach, combining a novel in silico search method with in vivo experimental validation, and will be useful in future research aimed at identifying and authenticating bacterial effector proteins from other genera.

Transcriptomics and differential gene expression in Whitmania pigra (Annelida: Clitellata: Hirudinida: Hirudinidae): Contrasting feeding and fasting modes

The medicinal utility of leeches has been demonstrated through decades of use in modern hospital settings, mainly as relievers of venous congestion following flap or digit replantation surgery. In the present study, we sequence and annotate (through BLAST- and Gene Ontology-based approaches) the salivary transcriptome of the nonblood feeding hirudinid Whitmania pigra and assess the differential gene expression of anticoagulation factors (through both quantitative real-time PCR [qRT-PCR] and in silico-based methods) during feeding and fasting conditions. This was done in order to evince the diversity of putative anticoagulation factors, as well as estimate the levels of upregulation of genes immediately after feeding. In total, we found sequences with demonstrated orthology (via both phylogenetic analyses and BLAST-based approaches) to seven different proteins that have previously been linked to anticoagulatory capabilities-eglin C, bdellin, granulin, guamerin, hyaluronidase, destabilase I, and lipocalin. All of these were recovered from leeches both in the fasting and in the feeding conditions, but all show signs of upregulation in the feeding leeches. Interestingly, our RNA-seq effort, coupled with a hypergeometric test, indicated that the differentially expressed genes were disproportionately involved in three main immunological pathways (endocytosis, peroxisome regulation, and lysosome regulation). The results and implications of the finding of anticoagulants in this nonblood feeding leech and the putative upregulation of anticoagulation factors after feeding are briefly discussed in an evolutionary context.

The Grayling Genome Reveals Selection on Gene Expression Regulation after Whole-Genome Duplication

Whole-genome duplication (WGD) has been a major evolutionary driver of increased genomic complexity in vertebrates. One such event occurred in the salmonid family ∼80 Ma (Ss4R) giving rise to a plethora of structural and regulatory duplicate-driven divergence, making salmonids an exemplary system to investigate the evolutionary consequences of WGD. Here, we present a draft genome assembly of European grayling (Thymallus thymallus) and use this in a comparative framework to study evolution of gene regulation following WGD. Among the Ss4R duplicates identified in European grayling and Atlantic salmon (Salmo salar), one-third reflect nonneutral tissue expression evolution, with strong purifying selection, maintained over ∼50 Myr. Of these, the majority reflect conserved tissue regulation under strong selective constraints related to brain and neural-related functions, as well as higher-order protein–protein interactions. A small subset of the duplicates have evolved tissue regulatory expression divergence in a common ancestor, which have been subsequently conserved in both lineages, suggestive of adaptive divergence following WGD. These candidates for adaptive tissue expression divergence have elevated rates of protein coding- and promoter-sequence evolution and are enriched for immune- and lipid metabolism ontology terms. Lastly, lineage-specific duplicate divergence points toward underlying differences in adaptive pressures on expression regulation in the nonanadromous grayling versus the anadromous Atlantic salmon. Our findings enhance our understanding of the role of WGD in genome evolution and highlight cases of regulatory divergence of Ss4R duplicates, possibly related to a niche shift in early salmonid evolution.

Resculpting the binding pocket of APC superfamily LeuT-fold amino acid transporters

Amino acid transporters are essential components of prokaryote and eukaryote cells, possess distinct physiological functions, and differ markedly in substrate specificity. Amino acid transporters can be both drug targets and drug transporters (bioavailability, targeting) with many monogenic disorders resulting from dysfunctional membrane transport. The largest collection of amino acid transporters (including the mammalian SLC6, SLC7, SLC32, SLC36, and SLC38 families), across all kingdoms of life, is within the Amino acid-Polyamine-organoCation (APC) superfamily. The LeuT-fold is a paradigm structure for APC superfamily amino acid transporters and carriers of sugars, neurotransmitters, electrolytes, osmolytes, vitamins, micronutrients, signalling molecules, and organic and fatty acids. Each transporter is specific for a unique sub-set of solutes, specificity being determined by how well a substrate fits into each binding pocket. However, the molecular basis of substrate selectivity remains, by and large, elusive. Using an integrated computational and experimental approach, we demonstrate that a single position within the LeuT-fold can play a crucial role in determining substrate specificity in mammalian and arthropod amino acid transporters within the APC superfamily. Systematic mutation of the amino acid residue occupying the equivalent position to LeuT V104 titrates binding pocket space resulting in dramatic changes in substrate selectivity in exemplar APC amino acid transporters including PAT2 (SLC36A2) and SNAT5 (SLC38A5). Our work demonstrates how a single residue/site within an archetypal structural motif can alter substrate affinity and selectivity within this important superfamily of diverse membrane transporters.