Jack-of-all-trades: phenotypic plasticity facilitates the invasion of an alien slug species

Invasive alien species might benefit from phenotypic plasticity by being able to (i) maintain fitness in stressful environments (‘robust’), (ii) increase fitness in favourable environments (‘opportunistic’), or (iii) combine both abilities (‘robust and opportunistic’). Here, we applied this framework, for the first time, to an animal, the invasive slug, Arion lusitanicus, and tested (i) whether it has a more adaptive phenotypic plasticity compared with a congeneric native slug, Arion fuscus, and (ii) whether it is robust, opportunistic or both. During one year, we exposed specimens of both species to a range of temperatures along an altitudinal gradient (700–2400 m a.s.l.) and to high and low food levels, and we compared the responsiveness of two fitness traits: survival and egg production. During summer, the invasive species had a more adaptive phenotypic plasticity, and at high temperatures and low food levels, it survived better and produced more eggs than A. fuscus, representing the robust phenotype. During winter, A. lusitanicus displayed a less adaptive phenotype than A. fuscus. We show that the framework developed for plants is also very useful for a better mechanistic understanding of animal invasions. Warmer summers and milder winters might lead to an expansion of this invasive species to higher altitudes and enhance its spread in the lowlands, supporting the concern that global climate change will increase biological invasions.     

Dynamic information in uncertain and changing worlds

A general theory for information processing by organisms living in uncertain and changing worlds is developed. The three fundamental properties of the theory are: (i) the use of a memory parameter that allows the organism to forget the more distant past, (ii) a succinct representation of encounters and information and (iii) flexibility in the estimates of parameters by including the uncertainty in these estimates in a consistent manner. The theory is developed using Bayesian methods (but can also be applied to maximum likelihood estimation) and is applied to the encounter models standardly used in ecology (Poisson, binomial, and negative binomial). Two applications are discussed: (i) patch selection and the matching rule and (ii) superparasitism by a parasitoid.     

Genetic conflict,kin and the origins of novel genetic systems

Genetic conflict may have played an important role in the evolution of novel genetic systems. The ancestral system of eumendelian genetics is highly symmetrical. Those derived from it (e.g. thelytokous parthenogenesis, haplodiploidy and parent-specific allele expression) are more asymmetrical in the genetic role played by maternal versus paternal alleles. These asymmetries may have arisen from maternal–paternal genetic conflict, or cytonuclear conflict, or from an interaction between them. Asymmetric genetic systems are much more common in terrestrial and freshwater taxa than in marine taxa. We suggest three reasons for this, based on the relative inhospitability of terrestrial environments to three types of organism: (i) pathogens—departure from the marine realm meant escape from many pathogens and parasites, reducing the need for sexual reproduction; (ii) symbionts—symbionts are no more important in the terrestrial realm than the marine realm but are more likely to be obligately intracellular and vertically transmitted, making them more likely to disrupt their host''s genetic systems; (iii) Gametes and embryos—because neither gametes nor embryos can be shed into air as easily as into seawater, the mother''s body is a more important environment for both types of organisms in the terrestrial realm than in the marine realm. This environment of asymmetric kinship (with neighbours more closely related by maternal alleles than by paternal alleles) may have helped to drive asymmetries in expression and transmission.     

Hormonal correlates of individual quality in a long-lived bird: a test of the ‘corticosterone–fitness hypothesis’

Measuring individual quality in vertebrates is difficult. Focusing on allostasis mechanisms may be useful because they are functionally involved in the ability of an individual to survive and reproduce in its environment. Thus, a rise in stress hormones levels (corticosterone) occurs when an organism has to cope with challenging environmental conditions. This has recently led to the proposal of the ‘cort–fitness hypothesis’, which suggests that elevated baseline corticosterone levels should be found in individuals of poor quality that have difficulty coping with their environment. We tested this hypothesis by comparing an integrative measure of individual quality to baseline corticosterone in black-browed albatrosses (Thalassarche melanophrys). We found that individual baseline corticosterone levels were related to individual quality and highly repeatable from one breeding season to the next. Importantly, this relationship was found in males, but not in females. Therefore, we suggest that the relationship between quality and baseline corticosterone levels may depend on the environmental and energetic constraints that individuals have to cope with.     

Benchmarking of a Bayesian single cell RNAseq differential gene expression test for dose–response study designs

The application of single-cell RNA sequencing (scRNAseq) for the evaluation of chemicals, drugs, and food contaminants presents the opportunity to consider cellular heterogeneity in pharmacological and toxicological responses. Current differential gene expression analysis (DGEA) methods focus primarily on two group comparisons, not multi-group dose–response study designs used in safety assessments. To benchmark DGEA methods for dose–response scRNAseq experiments, we proposed a multiplicity corrected Bayesian testing approach and compare it against 8 other methods including two frequentist fit-for-purpose tests using simulated and experimental data. Our Bayesian test method outperformed all other tests for a broad range of accuracy metrics including control of false positive error rates. Most notable, the fit-for-purpose and standard multiple group DGEA methods were superior to the two group scRNAseq methods for dose–response study designs. Collectively, our benchmarking of DGEA methods demonstrates the importance in considering study design when determining the most appropriate test methods.     

Trait–fitness associations do not predict within-species phenotypic evolution over 2 million years

Long-term patterns of phenotypic change are the cumulative results of tens of thousands to millions of years of evolution. Yet, empirical and theoretical studies of phenotypic selection are largely based on contemporary populations. The challenges in studying phenotypic evolution, in particular trait–fitness associations in the deep past, are barriers to linking micro- and macroevolution. Here, we capitalize on the unique opportunity offered by a marine colonial organism commonly preserved in the fossil record to investigate trait–fitness associations over 2 Myr. We use the density of female polymorphs in colonies of Antartothoa tongima as a proxy for fecundity, a fitness component, and investigate multivariate signals of trait–fitness associations in six time intervals on the backdrop of Pleistocene climatic shifts. We detect negative trait–fitness associations for feeding polymorph (autozooid) sizes, positive associations for autozooid shape but no particular relationship between fecundity and brood chamber size. In addition, we demonstrate that long-term trait patterns are explained by palaeoclimate (as approximated by ∂¹⁸O), and to a lesser extent by ecological interactions (i.e. overgrowth competition and substrate crowding). Our analyses show that macroevolutionary outcomes of trait evolution are not a simple scaling-up from the trait–fitness associations.     

BASiCS: Bayesian Analysis of Single-Cell Sequencing Data

Single-cell mRNA sequencing can uncover novel cell-to-cell heterogeneity in gene expression levels in seemingly homogeneous populations of cells. However, these experiments are prone to high levels of unexplained technical noise, creating new challenges for identifying genes that show genuine heterogeneous expression within the population of cells under study. BASiCS (Bayesian Analysis of Single-Cell Sequencing data) is an integrated Bayesian hierarchical model where: (i) cell-specific normalisation constants are estimated as part of the model parameters, (ii) technical variability is quantified based on spike-in genes that are artificially introduced to each analysed cell’s lysate and (iii) the total variability of the expression counts is decomposed into technical and biological components. BASiCS also provides an intuitive detection criterion for highly (or lowly) variable genes within the population of cells under study. This is formalised by means of tail posterior probabilities associated to high (or low) biological cell-to-cell variance contributions, quantities that can be easily interpreted by users. We demonstrate our method using gene expression measurements from mouse Embryonic Stem Cells. Cross-validation and meaningful enrichment of gene ontology categories within genes classified as highly (or lowly) variable supports the efficacy of our approach.     

Post-GWAS: where next? More samples,more SNPs or more biology?

The power of genome-wide association studies (GWAS) rests on several foundations: (i) there is a significant amount of additive genetic variation, (ii) individual causal polymorphisms often have sizable effects and (iii) they segregate at moderate-to-intermediate frequencies, or will be effectively ‘tagged'' by polymorphisms that do. Each of these assumptions has recently been questioned. (i) Why should genetic variation appear additive given that the underlying molecular networks are highly nonlinear? (ii) A new generation of relatedness-based analyses directs us back to the nearly infinitesimal model for effect sizes that quantitative genetics was long based upon. (iii) Larger effect causal polymorphisms are often low frequency, as selection might lead us to expect. Here, we review these issues and other findings that appear to question many of the foundations of the optimism GWAS prompted. We then present a roadmap emerging as one possible future for quantitative genetics. We argue that in future GWAS should move beyond purely statistical grounds. One promising approach is to build upon the combination of population genetic models and molecular biological knowledge. This combined treatment, however, requires fitting experimental data to models that are very complex, as well as accurate capturing of the uncertainty of resulting inference. This problem can be resolved through Bayesian analysis and tools such as approximate Bayesian computation—a method growing in popularity in population genetic analysis. We show a case example of anterior–posterior segmentation in Drosophila, and argue that similar approaches will be helpful as a GWAS augmentation, in human and agricultural research.     

Probabilistic Daily ILI Syndromic Surveillance with a Spatio-Temporal Bayesian Hierarchical Model

Background

For daily syndromic surveillance to be effective, an efficient and sensible algorithm would be expected to detect aberrations in influenza illness, and alert public health workers prior to any impending epidemic. This detection or alert surely contains uncertainty, and thus should be evaluated with a proper probabilistic measure. However, traditional monitoring mechanisms simply provide a binary alert, failing to adequately address this uncertainty.

Methods and Findings

Based on the Bayesian posterior probability of influenza-like illness (ILI) visits, the intensity of outbreak can be directly assessed. The numbers of daily emergency room ILI visits at five community hospitals in Taipei City during 2006–2007 were collected and fitted with a Bayesian hierarchical model containing meteorological factors such as temperature and vapor pressure, spatial interaction with conditional autoregressive structure, weekend and holiday effects, seasonality factors, and previous ILI visits. The proposed algorithm recommends an alert for action if the posterior probability is larger than 70%. External data from January to February of 2008 were retained for validation. The decision rule detects successfully the peak in the validation period. When comparing the posterior probability evaluation with the modified Cusum method, results show that the proposed method is able to detect the signals 1–2 days prior to the rise of ILI visits.

Conclusions

This Bayesian hierarchical model not only constitutes a dynamic surveillance system but also constructs a stochastic evaluation of the need to call for alert. The monitoring mechanism provides earlier detection as well as a complementary tool for current surveillance programs.     

Evolutionary Tuning of Protein Expression Levels of a Positively Autoregulated Two-Component System

Cellular adaptation relies on the development of proper regulatory schemes for accurate control of gene expression levels in response to environmental cues. Over- or under-expression can lead to diminished cell fitness due to increased costs or insufficient benefits. Positive autoregulation is a common regulatory scheme that controls protein expression levels and gives rise to essential features in diverse signaling systems, yet its roles in cell fitness are less understood. It remains largely unknown how much protein expression is ‘appropriate’ for optimal cell fitness under specific extracellular conditions and how the dynamic environment shapes the regulatory scheme to reach appropriate expression levels. Here, we investigate the correlation of cell fitness and output response with protein expression levels of the E. coli PhoB/PhoR two-component system (TCS). In response to phosphate (Pi)-depletion, the PhoB/PhoR system activates genes involved in phosphorus assimilation as well as genes encoding themselves, similarly to many other positively autoregulated TCSs. We developed a bacteria competition assay in continuous cultures and discovered that different Pi conditions have conflicting requirements of protein expression levels for optimal cell fitness. Pi-replete conditions favored cells with low levels of PhoB/PhoR while Pi-deplete conditions selected for cells with high levels of PhoB/PhoR. These two levels matched PhoB/PhoR concentrations achieved via positive autoregulation in wild-type cells under Pi-replete and -deplete conditions, respectively. The fitness optimum correlates with the wild-type expression level, above which the phosphorylation output saturates, thus further increase in expression presumably provides no additional benefits. Laboratory evolution experiments further indicate that cells with non-ideal protein levels can evolve toward the optimal levels with diverse mutational strategies. Our results suggest that the natural protein expression levels and feedback regulatory schemes of TCSs are evolved to match the phosphorylation output of the system, which is determined by intrinsic activities of TCS proteins.     

No evidence that relatedness or familiarity modulates male harm in Drosophila melanogaster flies from a wild population

Sexual selection frequently promotes the evolution of aggressive behaviors that help males compete against their rivals, but which may harm females and hamper their fitness. Kin selection theory predicts that optimal male–male competition levels can be reduced when competitors are more genetically related to each other than to the population average, contributing to resolve this sexual conflict. Work in Drosophila melanogaster has spearheaded empirical tests of this idea, but studies so far have been conducted in laboratory‐adapted populations in homogeneous rearing environments that may hamper kin recognition, and used highly skewed sex ratios that may fail to reflect average natural conditions. Here, we performed a fully factorial design with the aim of exploring how rearing environment (i.e., familiarity) and relatedness affect male–male aggression, male harassment, and overall male harm levels in flies from a wild population of Drosophila melanogaster, under more natural conditions. Namely, we (a) manipulated relatedness and familiarity so that larvae reared apart were raised in different environments, as is common in the wild, and (b) studied the effects of relatedness and familiarity under average levels of male–male competition in the field. We show that, contrary to previous findings, groups of unrelated‐unfamiliar males were as likely to fight with each other and harass females than related‐familiar males and that overall levels of male harm to females were similar across treatments. Our results suggest that the role of kin selection in modulating sexual conflict is yet unclear in Drosophila melanogaster, and call for further studies that focus on natural populations and realistic socio‐sexual and ecological environments.     

Estimation of Transmission Parameters of H5N1 Avian Influenza Virus in Chickens

Despite considerable research efforts, little is yet known about key epidemiological parameters of H5N1 highly pathogenic influenza viruses in their avian hosts. Here we show how these parameters can be estimated using a limited number of birds in experimental transmission studies. Our quantitative estimates, based on Bayesian methods of inference, reveal that (i) the period of latency of H5N1 influenza virus in unvaccinated chickens is short (mean: 0.24 days; 95% credible interval: 0.099–0.48 days); (ii) the infectious period of H5N1 virus in unvaccinated chickens is approximately 2 days (mean: 2.1 days; 95%CI: 1.8–2.3 days); (iii) the reproduction number of H5N1 virus in unvaccinated chickens need not be high (mean: 1.6; 95%CI: 0.90–2.5), although the virus is expected to spread rapidly because it has a short generation interval in unvaccinated chickens (mean: 1.3 days; 95%CI: 1.0–1.5 days); and (iv) vaccination with genetically and antigenically distant H5N2 vaccines can effectively halt transmission. Simulations based on the estimated parameters indicate that herd immunity may be obtained if at least 80% of chickens in a flock are vaccinated. We discuss the implications for the control of H5N1 avian influenza virus in areas where it is endemic.     

Host and brood parasite coevolutionary interactions covary with comparative patterns of the avian visual system

In coevolutionary arms-races, reciprocal ecological interactions and their fitness impacts shape the course of phenotypic evolution. The classic example of avian host–brood parasite interactions selects for host recognition and rejection of increasingly mimetic foreign eggs. An essential component of perceptual mimicry is that parasitic eggs escape detection by host sensory systems, yet there is no direct evidence that the avian visual system covaries with parasitic egg recognition or mimicry. Here, we used eye size measurements collected from preserved museum specimens as a metric of the avian visual system for species involved in host–brood parasite interactions. We discovered that (i) hosts had smaller eyes compared with non-hosts, (ii) parasites had larger eyes compared with hosts before but not after phylogenetic corrections, perhaps owing to the limited number of independent evolutionary origins of obligate brood parasitism, (iii) egg rejection in hosts with non-mimetic parasitic eggs positively correlated with eye size, and (iv) eye size was positively associated with increased avian-perceived host–parasite eggshell similarity. These results imply that both host-use by parasites and anti-parasitic responses by hosts covary with a metric of the visual system across relevant bird species, providing comparative evidence for coevolutionary patterns of host and brood parasite sensory systems.     

Social familiarity modulates personality trait in a cichlid fish

Personality traits, such as exploration–avoidance, are expected to be adaptive in a given context (e.g. low-risk environment) but to be maladaptive in others (e.g. high-risk environment). Therefore, it is expected that personality traits are flexible and respond to environmental fluctuations, given that consistency across different contexts is maintained, so that the relative individual responses in relation to others remains the same (i.e. although the magnitude of the response varies the differences between high and low responders are kept). Here, we tested the response of male cichlid fish (Oreochromis mossambicus) to a novel object (NO) in three different social contexts: (i) social isolation, (ii) in the presence of an unfamiliar conspecific, and (iii) in the presence of a familiar conspecific. Males in the familiar treatment exhibited more exploratory behaviour and less neophobia than males in either the unfamiliar or the social isolation treatments. However, there were no overall correlations in individual behaviour across the three treatments, suggesting a lack of consistency in exploration–avoidance as measured by the NO test in this species. Moreover, there were no differences in cortisol responsiveness to an acute stressor between the three treatments. Together, these results illustrate how behavioural traits usually taken as measures of personality may exhibit significant flexibility and lack the expected consistency across different social contexts.     

Rapid Renal Alpha-1 Antitrypsin Gene Induction in Experimental and Clinical Acute Kidney Injury

Alpha-1-antitrypsin (AAT) is a hepatic stress protein with protease inhibitor activity. Recent evidence indicates that ischemic or toxic injury can evoke selective changes within kidney that resemble a hepatic phenotype. Hence, we tested the following: i) Does acute kidney injury (AKI) up-regulate the normally renal silent AAT gene? ii) Does rapid urinary AAT excretion result? And iii) Can AAT''s anti-protease/anti-neutrophil elastase (NE) activity protect injured proximal tubule cells? CD-1 mice were subjected to ischemic or nephrotoxic (glycerol, maleate, cisplatin) AKI. Renal functional and biochemical assessments were made 4–72 hrs later. Rapidly following injury, 5–10 fold renal cortical and isolated proximal tubule AAT mRNA and protein increases occurred. These were paralleled by rapid (>100 fold) increases in urinary AAT excretion. AKI also induced marked increases in renal cortical/isolated proximal tubule NE mRNA. However, sharp NE protein levels declines resulted, which strikingly correlated (r, −0.94) with rising AAT protein levels (reflecting NE complexing by AAT/destruction). NE addition to HK-2 cells evoked ∼95% cell death. AAT completely blocked this NE toxicity, as well as Fe induced oxidant HK-2 cell attack. Translational relevance of experimental AAT gene induction was indicated by ∼100–1000 fold urinary AAT increases in 22 AKI patients (matching urine NGAL increases). We conclude: i) AKI rapidly up-regulates the renal cortical/proximal tubule AAT gene; ii) NE gene induction also results; iii) AAT can confer cytoprotection, potentially by blocking/reducing cytotoxic NE accumulation; and iv) marked increases in urinary AAT excretion in AKI patients implies clinical relevance of the AKI- AAT induction pathway.     

Transitional fossils and the origin of turtles

The origin of turtles is one of the most contentious issues in systematics with three currently viable hypotheses: turtles as the extant sister to (i) the crocodile–bird clade, (ii) the lizard–tuatara clade, or (iii) Diapsida (a clade composed of (i) and (ii)). We reanalysed a recent dataset that allied turtles with the lizard–tuatara clade and found that the inclusion of the stem turtle Proganochelys quenstedti and the ‘parareptile’ Eunotosaurus africanus results in a single overriding morphological signal, with turtles outside Diapsida. This result reflects the importance of transitional fossils when long branches separate crown clades, and highlights unexplored issues such as the role of topological congruence when using fossils to calibrate molecular clocks.     

Effects of dietary macronutrients on the hepatic transcriptome and serum metabolome in mice

Dietary macronutrient composition influences both hepatic function and aging. Previous work suggested that longevity and hepatic gene expression levels were highly responsive to dietary protein, but almost unaffected by other macronutrients. In contrast, we found expression of 4005, 4232, and 4292 genes in the livers of mice were significantly associated with changes in dietary protein (5%–30%), fat (20%–60%), and carbohydrate (10%–75%), respectively. More genes in aging‐related pathways (notably mTOR, IGF‐1, and NF‐kappaB) had significant correlations with dietary fat intake than protein and carbohydrate intake, and the pattern of gene expression changes in relation to dietary fat intake was in the opposite direction to the effect of graded levels of caloric restriction consistent with dietary fat having a negative impact on aging. We found 732, 808, and 995 serum metabolites were significantly correlated with dietary protein (5%–30%), fat (8.3%–80%), and carbohydrate (10%–80%) contents, respectively. Metabolomics pathway analysis revealed sphingosine‐1‐phosphate signaling was the significantly affected pathway by dietary fat content which has also been identified as significant changed metabolic pathway in the previous caloric restriction study. Our results suggest dietary fat has major impact on aging‐related gene and metabolic pathways compared with other macronutrients.     

Seasonal time bombs: dominant temperate viruses affect Southern Ocean microbial dynamics

Rapid warming in the highly productive western Antarctic Peninsula (WAP) region of the Southern Ocean has affected multiple trophic levels, yet viral influences on microbial processes and ecosystem function remain understudied in the Southern Ocean. Here we use cultivation-independent quantitative ecological and metagenomic assays, combined with new comparative bioinformatic techniques, to investigate double-stranded DNA viruses during the WAP spring–summer transition. This study demonstrates that (i) temperate viruses dominate this region, switching from lysogeny to lytic replication as bacterial production increases, and (ii) Southern Ocean viral assemblages are genetically distinct from lower-latitude assemblages, primarily driven by this temperate viral dominance. This new information suggests fundamentally different virus–host interactions in polar environments, where intense seasonal changes in bacterial production select for temperate viruses because of increased fitness imparted by the ability to switch replication strategies in response to resource availability. Further, temperate viral dominance may provide mechanisms (for example, bacterial mortality resulting from prophage induction) that help explain observed temporal delays between, and lower ratios of, bacterial and primary production in polar versus lower-latitude marine ecosystems. Together these results suggest that temperate virus–host interactions are critical to predicting changes in microbial dynamics brought on by warming in polar marine systems.