Niche theory‐based modeling of assembly processes of viral communities in bats

Understanding the assembly processes of symbiont communities, including viromes and microbiomes, is important for improving predictions on symbionts’ biogeography and disease ecology. Here, we use phylogenetic, functional, and geographic filters to predict the similarity between symbiont communities, using as a test case the assembly process in viral communities of Mexican bats. We construct generalized linear models to predict viral community similarity, as measured by the Jaccard index, as a function of differences in host phylogeny, host functionality, and spatial co‐occurrence, evaluating the models using the Akaike information criterion. Two model classes are constructed: a “known” model, where virus–host relationships are based only on data reported in Mexico, and a “potential” model, where viral reports of all the Americas are used, but then applied only to bat species that are distributed in Mexico. Although the “known” model shows only weak dependence on any of the filters, the “potential” model highlights the importance of all three filter types—phylogeny, functional traits, and co‐occurrence—in the assemblage of viral communities. The differences between the “known” and “potential” models highlight the utility of modeling at different “scales” so as to compare and contrast known information at one scale to another one, where, for example, virus information associated with bats is much scarcer.     

DNA Barcoding and Species Boundary Delimitation of Selected Species of Chinese Acridoidea (Orthoptera: Caelifera)

We tested the performance of DNA barcoding in Acridoidea and attempted to solve species boundary delimitation problems in selected groups using COI barcodes. Three analysis methods were applied to reconstruct the phylogeny. K2P distances were used to assess the overlap range between intraspecific variation and interspecific divergence. “Best match (BM)”, “best close match (BCM)”, “all species barcodes (ASB)” and “back-propagation neural networks (BP-based method)” were utilized to test the success rate of species identification. Phylogenetic species concept and network analysis were employed to delimitate the species boundary in eight selected species groups. The results demonstrated that the COI barcode region performed better in phylogenetic reconstruction at genus and species levels than at higher-levels, but showed a little improvement in resolving the higher-level relationships when the third base data or both first and third base data were excluded. Most overlaps and incorrect identifications may be due to imperfect taxonomy, indicating the critical role of taxonomic revision in DNA barcoding study. Species boundary delimitation confirmed the presence of oversplitting in six species groups and suggested that each group should be treated as a single species.     

Darwin's Manufactory Hypothesis Is Confirmed and Predicts the Extinction Risk of Extant Birds

In the Origin of Species Darwin hypothesized that the “manufactory” of species operates at different rates in different lineages and that the richness of taxonomic units is autocorrelated across levels of the taxonomic hierarchy. We confirm the manufactory hypothesis using a database of all the world''s extant avian subspecies, species and genera. The hypothesis is confirmed both in correlations across all genera and in paired comparisons controlling for phylogeny. We also find that the modern risk of extinction, as measured by “Red List” classifications, differs across the different categories of genera identified by Darwin. Specifically, species in “manufactory” genera are less likely to be threatened, endangered or recently extinct than are “weak manufactory” genera. Therefore, although Darwin used his hypothesis to investigate past evolutionary processes, we find that the hypothesis also foreshadows future changes to the evolutionary tree.     

Evolutionary Relationships of Three New Species of Enterobacteriaceae Living as Symbionts of Aphids and Other Insects

Ecological studies on three bacterial lineages symbiotic in aphids have shown that they impose a variety of effects on their hosts, including resistance to parasitoids and tolerance to heat stress. Phylogenetic analyses of partial sequences of gyrB and recA are consistent with previous analyses limited to 16S rRNA gene sequences and yield improved confidence of the evolutionary relationships of these symbionts. All three symbionts are in the Enterobacteriaceae. One of the symbionts, here given the provisional designation “Candidatus Serratia symbiotica,” is a Serratia species that has acquired a symbiotic lifestyle. The other two symbionts, here designated “Candidatus Hamiltonella defensa” and “Candidatus Regiella insecticola,” are sister groups to one another and together show a relationship to species of Photorhabdus.     

Bayesian Analysis of Congruence of Core Genes in Prochlorococcus and Synechococcus and Implications on Horizontal Gene Transfer

It is often suggested that horizontal gene transfer is so ubiquitous in microbes that the concept of a phylogenetic tree representing the pattern of vertical inheritance is oversimplified or even positively misleading. “Universal proteins” have been used to infer the organismal phylogeny, but have been criticized as being only the “tree of one percent.” Currently, few options exist for those wishing to rigorously assess how well a universal protein phylogeny, based on a relative handful of well-conserved genes, represents the phylogenetic histories of hundreds of genes. Here, we address this problem by proposing a visualization method and a statistical test within a Bayesian framework. We use the genomes of marine cyanobacteria, a group thought to exhibit substantial amounts of HGT, as a test case. We take 379 orthologous gene families from 28 cyanobacteria genomes and estimate the Bayesian posterior distributions of trees – a “treecloud” – for each, as well as for a concatenated dataset based on putative “universal proteins.” We then calculate the average distance between trees within and between all treeclouds on various metrics and visualize this high-dimensional space with non-metric multidimensional scaling (NMMDS). We show that the tree space is strongly clustered and that the universal protein treecloud is statistically significantly closer to the center of this tree space than any individual gene treecloud. We apply several commonly-used tests for incongruence/HGT and show that they agree HGT is rare in this dataset, but make different choices about which genes were subject to HGT. Our results show that the question of the representativeness of the “tree of one percent” is a quantitative empirical question, and that the phylogenetic central tendency is a meaningful observation even if many individual genes disagree due to the various sources of incongruence.     

Stable-Isotope Probing of Microorganisms Thriving at Thermodynamic Limits: Syntrophic Propionate Oxidation in Flooded Soil

Propionate is an important intermediate of the degradation of organic matter in many anoxic environments. In methanogenic environments, due to thermodynamic constraints, the oxidation of propionate requires syntrophic cooperation of propionate-fermenting proton-reducing bacteria and H₂-consuming methanogens. We have identified here microorganisms that were active in syntrophic propionate oxidation in anoxic paddy soil by rRNA-based stable-isotope probing (SIP). After 7 weeks of incubation with [¹³C]propionate (<10 mM) and the oxidation of ~30 μmol of ¹³C-labeled substrate per g dry weight of soil, we found that archaeal nucleic acids were ¹³C labeled to a larger extent than those of the bacterial partners. Nevertheless, both terminal restriction fragment length polymorphism and cloning analyses revealed Syntrophobacter spp., Smithella spp., and the novel Pelotomaculum spp. to predominate in “heavy” ¹³C-labeled bacterial rRNA, clearly showing that these were active in situ in syntrophic propionate oxidation. Among the Archaea, mostly Methanobacterium and Methanosarcina spp. and also members of the yet-uncultured “rice cluster I” lineage had incorporated substantial amounts of ¹³C label, suggesting that these methanogens were directly involved in syntrophic associations and/or thriving on the [¹³C]acetate released by the syntrophs. With this first application of SIP in an anoxic soil environment, we were able to clearly demonstrate that even guilds of microorganisms growing under thermodynamic constraints, as well as phylogenetically diverse syntrophic associations, can be identified by using SIP. This approach holds great promise for determining the structure and function relationships of further syntrophic or other nutritional associations in natural environments and for defining metabolic functions of yet-uncultivated microorganisms.     

Classification of Domain Movements in Proteins Using Dynamic Contact Graphs

A new method for the classification of domain movements in proteins is described and applied to 1822 pairs of structures from the Protein Data Bank that represent a domain movement in two-domain proteins. The method is based on changes in contacts between residues from the two domains in moving from one conformation to the other. We argue that there are five types of elemental contact changes and that these relate to five model domain movements called: “free”, “open-closed”, “anchored”, “sliding-twist”, and “see-saw.” A directed graph is introduced called the “Dynamic Contact Graph” which represents the contact changes in a domain movement. In many cases a graph, or part of a graph, provides a clear visual metaphor for the movement it represents and is a motif that can be easily recognised. The Dynamic Contact Graphs are often comprised of disconnected subgraphs indicating independent regions which may play different roles in the domain movement. The Dynamic Contact Graph for each domain movement is decomposed into elemental Dynamic Contact Graphs, those that represent elemental contact changes, allowing us to count the number of instances of each type of elemental contact change in the domain movement. This naturally leads to sixteen classes into which the 1822 domain movements are classified.     

No need to replace an “anomalous” primate (Primates) with an “anomalous” bear (Carnivora,Ursidae)

By means of mitochondrial 12S rRNA sequencing of putative “yeti”, “bigfoot”, and other “anomalous primate” hair samples, a recent study concluded that two samples, presented as from the Himalayas, do not belong to an “anomalous primate”, but to an unknown, anomalous type of ursid. That is, that they match 12S rRNA sequences of a fossil Polar Bear (Ursus maritimus), but neither of modern Polar Bears, nor of Brown Bears (Ursus arctos), the closest relative of Polar Bears, and one that occurs today in the Himalayas. We have undertaken direct comparison of sequences; replication of the original comparative study; inference of phylogenetic relationships of the two samples with respect to those from all extant species of Ursidae (except for the Giant Panda, Ailuropoda melanoleuca) and two extinct Pleistocene species; and application of a non-tree-based population aggregation approach for species diagnosis and identification. Our results demonstrate that the very short fragment of the 12S rRNA gene sequenced by Sykes et al. is not sufficiently informative to support the hypotheses provided by these authors with respect to the taxonomic identity of the individuals from which these sequences were obtained. We have concluded that there is no reason to believe that the two samples came from anything other than Brown Bears. These analyses afforded an opportunity to test the monophyly of morphologically defined species and to comment on both their phylogenetic relationships and future efforts necessary to advance our understanding of ursid systematics.     

Segregating the Effects of Seed Traits and Common Ancestry of Hardwood Trees on Eastern Gray Squirrel Foraging Decisions

The evolution of specific seed traits in scatter-hoarded tree species often has been attributed to granivore foraging behavior. However, the degree to which foraging investments and seed traits correlate with phylogenetic relationships among trees remains unexplored. We presented seeds of 23 different hardwood tree species (families Betulaceae, Fagaceae, Juglandaceae) to eastern gray squirrels (Sciurus carolinensis), and measured the time and distance travelled by squirrels that consumed or cached each seed. We estimated 11 physical and chemical seed traits for each species, and the phylogenetic relationships between the 23 hardwood trees. Variance partitioning revealed that considerable variation in foraging investment was attributable to seed traits alone (27–73%), and combined effects of seed traits and phylogeny of hardwood trees (5–55%). A phylogenetic PCA (pPCA) on seed traits and tree phylogeny resulted in 2 “global” axes of traits that were phylogenetically autocorrelated at the family and genus level and a third “local” axis in which traits were not phylogenetically autocorrelated. Collectively, these axes explained 30–76% of the variation in squirrel foraging investments. The first global pPCA axis, which produced large scores for seed species with thin shells, low lipid and high carbohydrate content, was negatively related to time to consume and cache seeds and travel distance to cache. The second global pPCA axis, which produced large scores for seeds with high protein, low tannin and low dormancy levels, was an important predictor of consumption time only. The local pPCA axis primarily reflected kernel mass. Although it explained only 12% of the variation in trait space and was not autocorrelated among phylogenetic clades, the local axis was related to all four squirrel foraging investments. Squirrel foraging behaviors are influenced by a combination of phylogenetically conserved and more evolutionarily labile seed traits that is consistent with a weak or more diffuse coevolutionary relationship between rodents and hardwood trees rather than a direct coevolutionary relationship.     

Macroevolutionary Dynamics and Historical Biogeography of Primate Diversification Inferred from a Species Supermatrix

Phylogenetic relationships, divergence times, and patterns of biogeographic descent among primate species are both complex and contentious. Here, we generate a robust molecular phylogeny for 70 primate genera and 367 primate species based on a concatenation of 69 nuclear gene segments and ten mitochondrial gene sequences, most of which were extracted from GenBank. Relaxed clock analyses of divergence times with 14 fossil-calibrated nodes suggest that living Primates last shared a common ancestor 71–63 Ma, and that divergences within both Strepsirrhini and Haplorhini are entirely post-Cretaceous. These results are consistent with the hypothesis that the Cretaceous-Paleogene mass extinction of non-avian dinosaurs played an important role in the diversification of placental mammals. Previous queries into primate historical biogeography have suggested Africa, Asia, Europe, or North America as the ancestral area of crown primates, but were based on methods that were coopted from phylogeny reconstruction. By contrast, we analyzed our molecular phylogeny with two methods that were developed explicitly for ancestral area reconstruction, and find support for the hypothesis that the most recent common ancestor of living Primates resided in Asia. Analyses of primate macroevolutionary dynamics provide support for a diversification rate increase in the late Miocene, possibly in response to elevated global mean temperatures, and are consistent with the fossil record. By contrast, diversification analyses failed to detect evidence for rate-shift changes near the Eocene-Oligocene boundary even though the fossil record provides clear evidence for a major turnover event (“Grande Coupure”) at this time. Our results highlight the power and limitations of inferring diversification dynamics from molecular phylogenies, as well as the sensitivity of diversification analyses to different species concepts.     

Frequent Cross-Species Transmission of Parvoviruses among Diverse Carnivore Hosts

Although parvoviruses are commonly described in domestic carnivores, little is known about their biodiversity in nondomestic species. A phylogenetic analysis of VP2 gene sequences from puma, coyote, gray wolf, bobcat, raccoon, and striped skunk revealed two major groups related to either feline panleukopenia virus (“FPV-like”) or canine parvovirus (“CPV-like”). Cross-species transmission was commonplace, with multiple introductions into each host species but, with the exception of raccoons, relatively little evidence for onward transmission in nondomestic species.     

Evolution of Pentameric Ligand-Gated Ion Channels: Pro-Loop Receptors

Pentameric ligand-gated ion channels (pLGICs) are ubiquitous neurotransmitter receptors in Bilateria, with a small number of known prokaryotic homologues. Here we describe a new inventory and phylogenetic analysis of pLGIC genes across all kingdoms of life. Our main finding is a set of pLGIC genes in unicellular eukaryotes, some of which are metazoan-like Cys-loop receptors, and others devoid of Cys-loop cysteines, like their prokaryotic relatives. A number of such “Cys-less” receptors also appears in invertebrate metazoans. Together, those findings draw a new distribution of pLGICs in eukaryotes. A broader distribution of prokaryotic channels also emerges, including a major new archaeal taxon, Thaumarchaeota. More generally, pLGICs now appear nearly ubiquitous in major taxonomic groups except multicellular plants and fungi. However, pLGICs are sparsely present in unicellular taxa, suggesting a high rate of gene loss and a non-essential character, contrasting with their essential role as synaptic receptors of the bilaterian nervous system. Multiple alignments of these highly divergent sequences reveal a small number of conserved residues clustered at the interface between the extracellular and transmembrane domains. Only the “Cys-loop” proline is absolutely conserved, suggesting the more fitting name “Pro loop” for that motif, and “Pro-loop receptors” for the superfamily. The infered molecular phylogeny shows a Cys-loop and a Cys-less clade in eukaryotes, both containing metazoans and unicellular members. This suggests new hypotheses on the evolutionary history of the superfamily, such as a possible origin of the Cys-loop cysteines in an ancient unicellular eukaryote. Deeper phylogenetic relationships remain uncertain, particularly around the split between bacteria, archaea, and eukaryotes.     

Developmentally Regulated Expression of the Gene Family for Cytosolic Glutamine Synthetase in Pisum sativum

In Pisum sativum, two classes of genes encode distinct isoforms of cytosolic glutamine synthetase (GS). The first class comprises two nearly identical or “twin” GS genes (GS341 and GS132), while the second comprises a single GS gene (GS299) distinct in both coding and noncoding regions from the “twin” GS genes. Gene-specific analyses were used to monitor the individual contribution of each gene for cytosolic GS during root nodule development and in cotyledons during germination, two contexts where large amounts of ammonia must be assimilated by GS for nitrogen transport. mRNAs corresponding to all three genes for cytosolic GS were shown to accumulate coordinately during a time course of nodule development. All the GS mRNAs also accumulate to wild-type levels in mutant nodules formed by a nifD⁻ strain of Rhizobium leguminosarum indicating that induced GS expression in pea root nodules does not depend on the production of ammonia. Distinct patterns of expression for the two classes of GS genes were observed in certain mutant root nodules and most dramatically in cotyledons of germinating seedlings. The different patterns of expression between the two classes of genes for cytosolic GS suggests that their distinct gene products may serve nonoverlapping functions during pea development.     

Relative Influence of Genetics and Shared Environment on Child Mental Health Symptoms Depends on Comorbidity

Background

Comorbidity among childhood mental health symptoms is common in clinical and community samples and should be accounted for when investigating etiology. We therefore aimed to uncover latent classes of mental health symptoms in middle childhood in a community sample, and to determine the latent genetic and environmental influences on those classes.

Methods

The sample comprised representative cohorts of twins. A questionnaire-based assessment of mental health symptoms was used in latent class analyses. Data on 3223 twins (1578 boys and 1645 girls) with a mean age of 7.5 years were analyzed. The sample was predominantly non-Hispanic Caucasian (92.1%).

Results

Latent class models delineated groups of children according to symptom profiles–not necessarily clinical groups but groups representing the general population, most with scores in the normative range. The best-fitting models suggested 9 classes for both girls and boys. Eight of the classes were very similar across sexes; these classes ranged from a “Low Symptom” class to a “Moderately Internalizing & Severely Externalizing” class. In addition, a “Moderately Anxious” class was identified for girls but not boys, and a “Severely Impulsive & Inattentive” class was identified for boys but not girls. Sex-combined analyses implicated moderate genetic influences for all classes. Shared environmental influences were moderate for the “Low Symptom” and “Moderately Internalizing & Severely Externalizing” classes, and small to zero for other classes.

Conclusions

We conclude that symptom classes are largely similar across sexes in middle childhood. Heritability was moderate for all classes, but shared environment played a greater role for classes in which no one type of symptom predominated.     

Evolution and Functional Characterisation of Melanopsins in a Deep-Sea Chimaera (Elephant Shark,Callorhinchus milii)

Non-visual photoreception in mammals is primarily mediated by two splice variants that derive from a single melanopsin (OPN4M) gene, whose expression is restricted to a subset of retinal ganglion cells. Physiologically, this sensory system regulates the photoentrainment of many biological rhythms, such as sleep via the melatonin endocrine system and pupil constriction. By contrast, melanopsin exists as two distinct lineages in non-mammals, opn4m and opn4x, and is broadly expressed in a wide range of tissue types, including the eye, brain, pineal gland and skin. Despite these findings, the evolution and function of melanopsin in early vertebrates are largely unknown. We, therefore, investigated the complement of opn4 classes present in the genome of a model deep-sea cartilaginous species, the elephant shark (Callorhinchus milii), as a representative vertebrate that resides at the base of the gnathostome (jawed vertebrate) lineage. We reveal that three melanopsin genes, opn4m1, opn4m2 and opn4x, are expressed in multiple tissues of the elephant shark. The two opn4m genes are likely to have arisen as a result of a lineage-specific duplication, whereas “long” and “short” splice variants are generated from a single opn4x gene. By using a heterologous expression system, we suggest that these genes encode functional photopigments that exhibit both “invertebrate-like” bistable and classical “vertebrate-like” monostable biochemical characteristics. We discuss the evolution and function of these melanopsin pigments within the context of the diverse photic and ecological environments inhabited by this chimaerid holocephalan, as well as the origin of non-visual sensory systems in early vertebrates.     

Network Signatures of Survival in Glioblastoma Multiforme

To determine a molecular basis for prognostic differences in glioblastoma multiforme (GBM), we employed a combinatorial network analysis framework to exhaustively search for molecular patterns in protein-protein interaction (PPI) networks. We identified a dysregulated molecular signature distinguishing short-term (survival<225 days) from long-term (survival>635 days) survivors of GBM using whole genome expression data from The Cancer Genome Atlas (TCGA). A 50-gene subnetwork signature achieved 80% prediction accuracy when tested against an independent gene expression dataset. Functional annotations for the subnetwork signature included “protein kinase cascade,” “IκB kinase/NFκB cascade,” and “regulation of programmed cell death” – all of which were not significant in signatures of existing subtypes. Finally, we used label-free proteomics to examine how our subnetwork signature predicted protein level expression differences in an independent GBM cohort of 16 patients. We found that the genes discovered using network biology had a higher probability of dysregulated protein expression than either genes exhibiting individual differential expression or genes derived from known GBM subtypes. In particular, the long-term survivor subtype was characterized by increased protein expression of DNM1 and MAPK1 and decreased expression of HSPA9, PSMD3, and CANX. Overall, we demonstrate that the combinatorial analysis of gene expression data constrained by PPIs outlines an approach for the discovery of robust and translatable molecular signatures in GBM.     

Unification of the Globally Distributed Spindle-Shaped Viruses of the Archaea

Viruses with spindle-shaped virions are abundant in diverse environments. Over the years, such viruses have been isolated from a wide range of archaeal hosts. Evolutionary relationships between them remained enigmatic, however. Here, using structural proteins as markers, we define familial ties among these “dark horses” of the virosphere and segregate all spindle-shaped viruses into two distinct evolutionary lineages, corresponding to Bicaudaviridae and Fuselloviridae. Our results illuminate the utility of structure-based virus classification and bring additional order to the virosphere.