DSAVE: Detection of misclassified cells in single-cell RNA-Seq data

Single-cell RNA sequencing has become a valuable tool for investigating cell types in complex tissues, where clustering of cells enables the identification and comparison of cell populations. Although many studies have sought to develop and compare different clustering approaches, a deeper investigation into the properties of the resulting populations is lacking. Specifically, the presence of misclassified cells can influence downstream analyses, highlighting the need to assess subpopulation purity and to detect such cells. We developed DSAVE (Down-SAmpling based Variation Estimation), a method to evaluate the purity of single-cell transcriptome clusters and to identify misclassified cells. The method utilizes down-sampling to eliminate differences in sampling noise and uses a log-likelihood based metric to help identify misclassified cells. In addition, DSAVE estimates the number of cells needed in a population to achieve a stable average gene expression profile within a certain gene expression range. We show that DSAVE can be used to find potentially misclassified cells that are not detectable by similar tools and reveal the cause of their divergence from the other cells, such as differing cell state or cell type. With the growing use of single-cell RNA-seq, we foresee that DSAVE will be an increasingly useful tool for comparing and purifying subpopulations in single-cell RNA-Seq datasets

Do contaminants originating from state-of-the-art treated wastewater impact the ecological quality of surface waters?

Since the 1980s, advances in wastewater treatment technology have led to considerably improved surface water quality in the urban areas of many high income countries. However, trace concentrations of organic wastewater-associated contaminants may still pose a key environmental hazard impairing the ecological quality of surface waters. To identify key impact factors, we analyzed the effects of a wide range of anthropogenic and environmental variables on the aquatic macroinvertebrate community. We assessed ecological water quality at 26 sampling sites in four urban German lowland river systems with a 0–100% load of state-of-the-art biological activated sludge treated wastewater. The chemical analysis suite comprised 12 organic contaminants (five phosphor organic flame retardants, two musk fragrances, bisphenol A, nonylphenol, octylphenol, diethyltoluamide, terbutryn), 16 polycyclic aromatic hydrocarbons, and 12 heavy metals. Non-metric multidimensional scaling identified organic contaminants that are mainly wastewater-associated (i.e., phosphor organic flame retardants, musk fragrances, and diethyltoluamide) as a major impact variable on macroinvertebrate species composition. The structural degradation of streams was also identified as a significant factor. Multiple linear regression models revealed a significant impact of organic contaminants on invertebrate populations, in particular on Ephemeroptera, Plecoptera, and Trichoptera species. Spearman rank correlation analyses confirmed wastewater-associated organic contaminants as the most significant variable negatively impacting the biodiversity of sensitive macroinvertebrate species. In addition to increased aquatic pollution with organic contaminants, a greater wastewater fraction was accompanied by a slight decrease in oxygen concentration and an increase in salinity. This study highlights the importance of reducing the wastewater-associated impact on surface waters. For aquatic ecosystems in urban areas this would lead to: (i) improvement of the ecological integrity, (ii) reduction of biodiversity loss, and (iii) faster achievement of objectives of legislative requirements, e.g., the European Water Framework Directive

Comparison of Argentinean Saint Louis Encephalitis Virus Non-Epidemic and Epidemic Strain Infections in an Avian Model

St. Louis encephalitis virus (SLEV, Flavivirus, Flaviviridae) is an emerging mosquito-borne pathogen in South America, with human SLEV encephalitis cases reported in Argentina and Brazil. Genotype III strains of SLEV were isolated from Culex quinquefasciatus mosquitoes in Cordoba, Argentina in 2005, during the largest SLEV outbreak ever reported in South America. The present study tested the hypothesis that the recent, epidemic SLEV strain exhibits greater virulence in birds as compared with a non-epidemic genotype III strain isolated from mosquitoes in Santa Fe Province 27 years earlier. The observed differences in infection parameters between adult House sparrows (Passer domesticus) that were needle-inoculated with either the epidemic or historic SLEV strain were not statistically significant. However, only the House sparrows that were infected with the epidemic strain achieved infectious-level viremia titers sufficient to infect Cx. spp. mosquitoes vectors. Furthermore, the vertebrate reservoir competence index values indicated an approximately 3-fold increase in amplification potential of House sparrows infected with the epidemic strain when pre-existing flavivirus-reactive antibodies were present, suggesting the possibility that antibody-dependent enhancement may increase the risk of avian-amplified transmission of SLEV in South America

Ras Conformational Switching: Simulating Nucleotide-Dependent Conformational Transitions with Accelerated Molecular Dynamics

Ras mediates signaling pathways controlling cell proliferation and development by cycling between GTP- and GDP-bound active and inactive conformational states. Understanding the complete reaction path of this conformational change and its intermediary structures is critical to understanding Ras signaling. We characterize nucleotide-dependent conformational transition using multiple-barrier-crossing accelerated molecular dynamics (aMD) simulations. These transitions, achieved for the first time for wild-type Ras, are impossible to observe with classical molecular dynamics (cMD) simulations due to the large energetic barrier between end states. Mapping the reaction path onto a conformer plot describing the distribution of the crystallographic structures enabled identification of highly populated intermediate structures. These structures have unique switch orientations (residues 25–40 and 57–75) intermediate between GTP and GDP states, or distinct loop3 (46–49), loop7 (105–110), and α5 C-terminus (159–166) conformations distal from the nucleotide-binding site. In addition, these barrier-crossing trajectories predict novel nucleotide-dependent correlated motions, including correlations of α2 (residues 66–74) with α3-loop7 (93–110), loop2 (26–37) with loop10 (145–151), and loop3 (46–49) with α5 (152–167). The interconversion between newly identified Ras conformations revealed by this study advances our mechanistic understanding of Ras function. In addition, the pattern of correlated motions provides new evidence for a dynamic linkage between the nucleotide-binding site and the membrane interacting C-terminus critical for the signaling function of Ras. Furthermore, normal mode analysis indicates that the dominant collective motion that occurs during nucleotide-dependent conformational exchange, and captured in aMD (but absent in cMD) simulations, is a low-frequency motion intrinsic to the structure

Morphological characteristics of motor neurons do not determine their relative susceptibility to degeneration in a mouse model of severe spinal muscular atrophy

Spinal muscular atrophy (SMA) is a leading genetic cause of infant mortality, resulting primarily from the degeneration and loss of lower motor neurons. Studies using mouse models of SMA have revealed widespread heterogeneity in the susceptibility of individual motor neurons to neurodegeneration, but the underlying reasons remain unclear. Data from related motor neuron diseases, such as amyotrophic lateral sclerosis (ALS), suggest that morphological properties of motor neurons may regulate susceptibility: in ALS larger motor units innervating fast-twitch muscles degenerate first. We therefore set out to determine whether intrinsic morphological characteristics of motor neurons influenced their relative vulnerability to SMA. Motor neuron vulnerability was mapped across 10 muscle groups in SMA mice. Neither the position of the muscle in the body, nor the fibre type of the muscle innervated, influenced susceptibility. Morphological properties of vulnerable and disease-resistant motor neurons were then determined from single motor units reconstructed in Thy.1-YFP-H mice. None of the parameters we investigated in healthy young adult mice - including motor unit size, motor unit arbor length, branching patterns, motor endplate size, developmental pruning and numbers of terminal Schwann cells at neuromuscular junctions - correlated with vulnerability. We conclude that morphological characteristics of motor neurons are not a major determinant of disease-susceptibility in SMA, in stark contrast to related forms of motor neuron disease such as ALS. This suggests that subtle molecular differences between motor neurons, or extrinsic factors arising from other cell types, are more likely to determine relative susceptibility in SMA

FAK acts as a suppressor of RTK-MAP kinase signalling in Drosophila melanogaster epithelia and human cancer cells

Receptor Tyrosine Kinases (RTKs) and Focal Adhesion Kinase (FAK) regulate multiple signalling pathways, including mitogen-activated protein (MAP) kinase pathway. FAK interacts with several RTKs but little is known about how FAK regulates their downstream signalling. Here we investigated how FAK regulates signalling resulting from the overexpression of the RTKs RET and EGFR. FAK suppressed RTKs signalling in Drosophila melanogaster epithelia by impairing MAPK pathway. This regulation was also observed in MDA-MB-231 human breast cancer cells, suggesting it is a conserved phenomenon in humans. Mechanistically, FAK reduced receptor recycling into the plasma membrane, which resulted in lower MAPK activation. Conversely, increasing the membrane pool of the receptor increased MAPK pathway signalling. FAK is widely considered as a therapeutic target in cancer biology; however, it also has tumour suppressor properties in some contexts. Therefore, the FAK-mediated negative regulation of RTK/MAPK signalling described here may have potential implications in the designing of therapy strategies for RTK-driven tumours

The Distinct Conformational Dynamics of K-Ras and H-Ras A59G

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species

Computational Cancer Biology: An Evolutionary Perspective

ISSN:1553-734XISSN:1553-735

The commonness of rarity: Global and future distribution of rarity across land plants

A key feature of life’s diversity is that some species are common but many more are rare. Nonetheless, at global scales, we do not know what fraction of biodiversity consists of rare species. Here, we present the largest compilation of global plant diversity to quantify the fraction of Earth’s plant biodiversity that are rare. A large fraction, ~36.5% of Earth’s ~435,000 plant species, are exceedingly rare. Sampling biases and prominent models, such as neutral theory and the k-niche model, cannot account for the observed prevalence of rarity. Our results indicate that (i) climatically more stable regions have harbored rare species and hence a large fraction of Earth’s plant species via reduced extinction risk but that (ii) climate change and human land use are now disproportionately impacting rare species. Estimates of global species abundance distributions have important implications for risk assessments and conservation planning in this era of rapid global change