The Signaling Pathway That cGAMP Riboswitches Found: Analysis and Application of Riboswitches to Study cGAMP Signaling in Geobacter sulfurreducens.

The Hypr cGAMP signaling pathway was discovered via the function of the riboswitch. In this study, we show the development of a method for affinity capture followed by sequencing to identify non-coding RNA regions that bind nucleotide signals such as cGAMP. The RNAseq of affinity-captured cGAMP riboswitches from the Geobacter sulfurreducens transcriptome highlights general challenges that remain for this technique. Furthermore, by applying riboswitch reporters in vivo, we identify new growth conditions and transposon mutations that affect cGAMP levels in G. sulfurreducens. This work reveals an extensive regulatory network and supports a second functional cGAMP synthase gene in G. sulfurreducens. The activity of the second synthase was validated using riboswitch-based fluorescent biosensors, and is the first known example of an active enzyme with a variant GGDDF motif.

Clines Arc through Multivariate Morphospace.

Evolutionary biologists typically represent clines as spatial gradients in a univariate character (or a principal-component axis) whose mean changes as a function of location along a transect spanning an environmental gradient or ecotone. This univariate approach may obscure the multivariate nature of phenotypic evolution across a landscape. Clines might instead be plotted as a series of vectors in multidimensional morphospace, connecting sequential geographic sites. We present a model showing that clines may trace nonlinear paths that arc through morphospace rather than elongating along a single major trajectory. Arcing clines arise because different characters diverge at different rates or locations along a geographic transect. We empirically confirm that some clines arc through morphospace, using morphological data from threespine stickleback sampled along eight independent transects from lakes down their respective outlet streams. In all eight clines, successive vectors of lake-stream divergence fluctuate in direction and magnitude in trait space, rather than pointing along a single phenotypic axis. Most clines exhibit surprisingly irregular directions of divergence as one moves downstream, although a few clines exhibit more directional arcs through morphospace. Our results highlight the multivariate complexity of clines that cannot be captured with the traditional graphical framework. We discuss hypotheses regarding the causes, and implications, of such arcing multivariate clines.

Gene expression stasis and plasticity following migration into a foreign environment.

Selection against migrants is key to maintaining genetic differences between populations linked by dispersal. However, migrants may mitigate fitness costs by proactively choosing among available habitats, or by phenotypic plasticity. We previously reported that a reciprocal transplant of lake and stream stickleback (Gasterosteus aculeatus) found little support for divergent selection. Here, we revisit that experiment to test whether phenotypic plasticity in gene expression may have helped migrants adjust to unfamiliar habitats. We measured gene expression profiles in stickleback via TagSeq and tested whether migrants between lake and stream habitats exhibited a plastic response to their new environment that allowed them to converge on the expression profile of adapted natives. We report extensive gene expression differences between genetically divergent lake and stream stickleback, despite gene flow. But for many genes, expression was highly plastic. Fish transplanted into the adjoining habitat partially converged on the expression profile typical of natives from their new habitat. This suggests that expression plasticity may soften the impact of migration. Nonetheless, lake and stream fish differed in survival rates and parasite infection rates in our study, implying that expression plasticity is not fast or extensive enough to fully homogenize fish performance.

The Caudal Skeleton of the Zebrafish, Danio rerio, from a Phylogenetic Perspective: A Polyural Interpretation of Homologous Structures.

The structure of the caudal skeleton of extant teleost fishes has been interpreted in two different ways. In a diural interpretation, a caudal skeleton is composed of two centra articulated with one to six hypurals. Most subsequent authors have followed this interpretation. In contrast, a polyural interpretation considers the teleost fin to be derived from a fully metameristic ancestral bauplan originally composed of a one-to-one relationship between neural arches, centra (when present), and hypurals. Three different interpretations of the identity and homology of skeletal components of the caudal skeleton of the teleost fish Danio rerio have been proposed, two from a diural perspective and one from a polyural perspective. We examine each caudal skeletal component of Danio rerio from both a developmental and phylogenetic perspective. We propose that a polyural interpretation of structures is consistent with the current interpretation of the basal neopterygian caudal fin for this model organism rather than the older diural interpretation that does not take into account the metamerism observed in caudal structures during development. The polyural interpretation suggests several shared evolutionary innovations of major clades that would remain undiscovered under the older diural naming paradigm and makes the terminology of the parts of the caudal fin of Danio rerio strictly comparable to more basal fishes.

Senescence-induced immune remodeling facilitates metastatic adrenal cancer in a sex-dimorphic manner.

Aging markedly increases cancer risk, yet our mechanistic understanding of how aging influences cancer initiation is limited. Here we demonstrate that the loss of ZNRF3, an inhibitor of Wnt signaling that is frequently mutated in adrenocortical carcinoma, leads to the induction of cellular senescence that remodels the tissue microenvironment and ultimately permits metastatic adrenal cancer in old animals. The effects are sexually dimorphic, with males exhibiting earlier senescence activation and a greater innate immune response, driven in part by androgens, resulting in high myeloid cell accumulation and lower incidence of malignancy. Conversely, females present a dampened immune response and increased susceptibility to metastatic cancer. Senescence-recruited myeloid cells become depleted as tumors progress, which is recapitulated in patients in whom a low myeloid signature is associated with worse outcomes. Our study uncovers a role for myeloid cells in restraining adrenal cancer with substantial prognostic value and provides a model for interrogating pleiotropic effects of cellular senescence in cancer.

Predation of Threespine Stickleback by Dragonfly Naiads.

Threespine stickleback (Gasterosteus aculeatus) populations that have evolved pelvic girdle reduction are most commonly found in lakes with low dissolved ion concentration, a lack of piscivorous fishes, and abundant macroinvertebrate predators. Researchers have speculated that macroinvertebrates have a propensity to consume prey with pelvic spines. If this is true, perhaps macroinvertebrates use the stickleback's spines to facilitate capture and manipulation. This study tested whether dragonfly naiads differentially prey upon stickleback possessing either a complete or reduced pelvis and documented naiad hunting and capturing behavior. Results from an arena experiment suggest that naiads do not prey more heavily upon individuals with one pelvic phenotype over the other. However, results from trials where the naiads were presented with one stickleback with pelvic spines and another without suggest that naiads prey more heavily upon small stickleback with pelvic spines and large stickleback without pelvic spines and that they adjust their predatory behavior based upon the pelvic phenotype of the prey.

A phospho-tyrosine-based signaling module using SPOP, CSK, and LYN controls TLR-induced IRF activity.

Toll-like receptors (TLRs) recognize pathogen- and host-derived factors and control immune responses via the adaptor protein MyD88 and members of the interferon regulatory transcription factor (IRF) family. IRFs orchestrate key effector functions, including cytokine release, cell differentiation, and, under certain circumstances, inflammation pathology. Here, we show that IRF activity is generically controlled by the Src kinase family member LYN, which phosphorylates all TLR-induced IRFs at a conserved tyrosine residue, resulting in K48-linked polyubiquitination and proteasomal degradation of IRFs. We further show that LYN activity is controlled by the upstream kinase C-terminal Src kinase (CSK), whose activity, in turn, is controlled by the adaptor protein SPOP, which serves as molecular bridge to recruit CSK into the TLR signaling complex and to activate CSK catalytic activity. Consistently, deletion of SPOP or CSK results in increased LYN activity, LYN-directed IRF degradation, and inhibition of IRF transcriptional activity. Together, the data reveal a key regulatory mechanism for IRF family members controlling TLR biology.

Evolutionary gain and loss of a pathological immune response to parasitism.

Parasites impose fitness costs on their hosts. Biologists often assume that natural selection favors infection-resistant hosts. Yet, when the immune response itself is costly, theory suggests that selection may sometimes favor loss of resistance, which may result in alternative stable states where some populations are resistant and others are tolerant. Intraspecific variation in immune costs is rarely surveyed in a manner that tests evolutionary patterns, and there are few examples of adaptive loss of resistance. Here, we show that when marine threespine stickleback colonized freshwater lakes, they gained resistance to the freshwater-associated cestode Schistocephalus solidus. Extensive peritoneal fibrosis and inflammation are a commonly observed phenotype that contributes to suppression of cestode growth and viability but also imposes a substantial cost on fecundity. Combining genetic mapping and population genomics, we find that opposing selection generates immune system differences between tolerant and resistant populations, consistent with divergent optimization.

Human melanocyte development and melanoma dedifferentiation at single-cell resolution.

In humans, epidermal melanocytes are responsible for skin pigmentation, defence against ultraviolet radiation and the deadliest common skin cancer, melanoma. Although there is substantial overlap in melanocyte development pathways between different model organisms, species-dependent differences are frequent and the conservation of these processes in human skin remains unresolved. Here, we used a single-cell enrichment and RNA-sequencing pipeline to study human epidermal melanocytes directly from the skin, capturing transcriptomes across different anatomical sites, developmental age, sexes and multiple skin tones. We uncovered subpopulations of melanocytes that exhibit anatomical site-specific enrichment that occurs during gestation and persists through adulthood. The transcriptional signature of the volar-enriched subpopulation is retained in acral melanomas. Furthermore, we identified human melanocyte differentiation transcriptional programs that are distinct from gene signatures generated from model systems. Finally, we used these programs to define patterns of dedifferentiation that are predictive of melanoma prognosis and response to immune checkpoint inhibitor therapy.

SV-plaudit: A cloud-based framework for manually curating thousands of structural variants.

SV-plaudit is a framework for rapidly curating structural variant (SV) predictions. For each SV, we generate an image that visualizes the coverage and alignment signals from a set of samples. Images are uploaded to our cloud framework where users assess the quality of each image using a client-side web application. Reports can then be generated as a tab-delimited file or annotated Variant Call Format (VCF) file. As a proof of principle, nine researchers collaborated for 1 hour to evaluate 1,350 SVs each. We anticipate that SV-plaudit will become a standard step in variant calling pipelines and the crowd-sourced curation of other biological results.Code available at https://github.com/jbelyeu/SV-plauditDemonstration video available at https://www.youtube.com/watch?v=ono8kHMKxDs.

Gene Expression Contributes to the Recent Evolution of Host Resistance in a Model Host Parasite System.

Heritable population differences in immune gene expression following infection can reveal mechanisms of host immune evolution. We compared gene expression in infected and uninfected threespine stickleback (Gasterosteus aculeatus) from two natural populations that differ in resistance to a native cestode parasite, Schistocephalus solidus. Genes in both the innate and adaptive immune system were differentially expressed as a function of host population, infection status, and their interaction. These genes were enriched for loci controlling immune functions known to differ between host populations or in response to infection. Coexpression network analysis identified two distinct processes contributing to resistance: parasite survival and suppression of growth. Comparing networks between populations showed resistant fish have a dynamic expression profile while susceptible fish are static. In summary, recent evolutionary divergence between two vertebrate populations has generated population-specific gene expression responses to parasite infection, affecting parasite establishment and growth.

Contrasting effects of environment and genetics generate a continuum of parallel evolution.

Parallel evolution of similar traits by independent populations in similar environments is considered strong evidence for adaptation by natural selection. Often, however, replicate populations in similar environments do not all evolve in the same way, thus deviating from any single, predominant outcome of evolution. This variation might arise from non-adaptive, population-specific effects of genetic drift, gene flow or limited genetic variation. Alternatively, these deviations from parallel evolution might also reflect predictable adaptation to cryptic environmental heterogeneity within discrete habitat categories. Here, we show that deviations from parallel evolution are the consequence of environmental variation within habitats combined with variation in gene flow. Threespine stickleback (Gasterosteus aculeatus) in adjoining lake and stream habitats (a lake-stream 'pair') diverge phenotypically, yet the direction and magnitude of this divergence is not always fully parallel among 16 replicate pairs. We found that the multivariate direction of lake-stream morphological divergence was less parallel between pairs whose environmental differences were less parallel. Thus, environmental heterogeneity among lake-stream pairs contributes to deviations from parallel evolution. Additionally, likely genomic targets of selection were more parallel between environmentally more similar pairs. In contrast, variation in the magnitude of lake-stream divergence (independent of direction) was better explained by differences in lake-stream gene flow; pairs with greater lake-stream gene flow were less morphologically diverged. Thus, both adaptive and non-adaptive processes work concurrently to generate a continuum of parallel evolution across lake-stream stickleback population pairs.

Evaluation of TagSeq, a reliable low-cost alternative for RNAseq.

RNAseq is a relatively new tool for ecological genetics that offers researchers insight into changes in gene expression in response to a myriad of natural or experimental conditions. However, standard RNAseq methods (e.g., Illumina TruSeq® or NEBNext® ) can be cost prohibitive, especially when study designs require large sample sizes. Consequently, RNAseq is often underused as a method, or is applied to small sample sizes that confer poor statistical power. Low cost RNAseq methods could therefore enable far greater and more powerful applications of transcriptomics in ecological genetics and beyond. Standard mRNAseq is costly partly because one sequences portions of the full length of all transcripts. Such whole-mRNA data are redundant for estimates of relative gene expression. TagSeq is an alternative method that focuses sequencing effort on mRNAs' 3' end, reducing the necessary sequencing depth per sample, and thus cost. We present a revised TagSeq library construction procedure, and compare its performance against NEBNext® , the 'gold-standard' whole mRNAseq method. We built both TagSeq and NEBNext® libraries from the same biological samples, each spiked with control RNAs. We found that TagSeq measured the control RNA distribution more accurately than NEBNext® , for a fraction of the cost per sample (~10%). The higher accuracy of TagSeq was particularly apparent for transcripts of moderate to low abundance. Technical replicates of TagSeq libraries are highly correlated, and were correlated with NEBNext® results. Overall, we show that our modified TagSeq protocol is an efficient alternative to traditional whole mRNAseq, offering researchers comparable data at greatly reduced cost.

The role of clonal interference in the evolutionary dynamics of plasmid-host adaptation.

UNLABELLED: Promiscuous plasmids replicate in a wide range of bacteria and therefore play a key role in the dissemination of various host-beneficial traits, including antibiotic resistance. Despite the medical relevance, little is known about the evolutionary dynamics through which drug resistance plasmids adapt to new hosts and thereby persist in the absence of antibiotics. We previously showed that the incompatibility group P-1 (IncP-1) minireplicon pMS0506 drastically improved its stability in novel host Shewanella oneidensis MR-1 after 1,000 generations under antibiotic selection for the plasmid. The only mutations found were those affecting the N terminus of the plasmid replication initiation protein TrfA1. Our aim in this study was to gain insight into the dynamics of plasmid evolution. Changes in stability and genotype frequencies of pMS0506 were monitored in evolving populations of MR-1 (pMS0506). Genotypes were determined by sequencing trfA1 amplicons from individual clones and by 454 pyrosequencing of whole plasmids from entire populations. Stability of pMS0506 drastically improved by generation 200. Many evolved plasmid genotypes with point mutations as well as in-frame and frameshift deletions and duplications in trfA1 were observed in all lineages with both sequencing methods. Strikingly, multiple genotypes were simultaneously present at high frequencies (>10%) in each population. Their relative abundances changed over time, but after 1,000 generations only one or two genotypes dominated the populations. This suggests that hosts with different plasmid genotypes were competing with each other, thus affecting the evolutionary trajectory. Plasmids can thus rapidly improve their stability, and clonal interference plays a significant role in plasmid-host adaptation dynamics. IMPORTANCE: Promiscuous plasmids play an important role in the spread of antibiotic resistance and many other traits between closely and distantly related bacteria. However, little is known about the dynamics by which these broad-host-range antibiotic resistance plasmids adapt to novel bacteria and thereby become more persistent, even in the absence of antibiotics. In this study, we show that after no more than 200 generations of growth in the presence of antibiotics, a plasmid that was initially poorly maintained in a novel bacterial host evolved to become drastically more persistent in the absence of antibiotics. In each of the evolving populations, an unexpectedly large number of bacterial variants arose with distinct mutations in the plasmid's replication initiation protein. Our results suggest that clonal interference, characterized by competition between variant clones in a population, plays a major role in the evolution of the persistence of drug resistance.