Search

Background: Gene sequencing is the 'gold standard' for resistance testing in the clinic. We present evidence that clinically misleading results arise with mixed virai populations due to intrasample variability. Methods: An HW+ patient taking nelfinavir presented with mixed populations of virus at codons 30, 63 and 90 of protease. RNA was re-extracted from plasma and 20 aliquots of cDNA prepared. Each aliquot, including the original cDNA, was sequenced twice, using different ABI310 Genetic Analysers and two primers. Electropherograms were compared for the 42 sequences. Results: All aliquots were sequenced successfully, but produced inconsistent results at mixed codons. A codon is uninterpretable ('N') if no peak dominates at one position. Codons 36, 71, 77 and 84 gave identical results with homozygous bases in all 42 samples. Codon 30 was wild-type (WT) in 29/42, mutant (MT) in 3/42 and 'N' in 10/42. Codon 46 was WT in 37/42 and 'N' in 5/42. Codon 63 was WT in 1/42, MT in 33/42 and 'N' in 8/42. Codon 90 was WT in 4/42, MT in 29/42 and 'N' in 9/21. Conclusions: The 42 sequences should have been identical, but in the presence of a mixed population of virus there is marked variability. Over-amplification of a minor WT population may result in masking of the true dominant MT codon. A false WT could result in a failing regimen being continued or an inactive drug being prescribed. [ABSTRACT FROM AUTHOR]

When species disperse into previously unoccupied habitats, new populations encounter unfamiliar species interactions such as altered parasite loads. Theory predicts that newly founded populations should exhibit destabilized eco-evolutionary fluctuations in infection rates and immune traits. However, to understand founder effects biologists typically rely on retrospective studies of range expansions, missing early-generation infection dynamics. To remedy this, we experimentally founded whole-lake populations of threespine stickleback. Infection rates were temporally stable in native source lakes. In contrast, newly founded populations exhibit destabilized host-parasite dynamics: high starting infection rates led to increases in a heritable immune trait (peritoneal fibrosis), suppressing infection rates. The resulting temporal auto-correlation between infection and immunity suggest that newly founded populations can exhibit rapid host-parasite eco-evolutionary dynamics., Competing Interests: Competing interests: Authors declare that they have no competing interests.

Eco-evolutionary experiments are typically conducted in semi-unnatural controlled settings, such as mesocosms; yet inferences about how evolution and ecology interact in the real world would surely benefit from experiments in natural uncontrolled settings. Opportunities for such experiments are rare but do arise in the context of restoration ecology-where different "types" of a given species can be introduced into different "replicate" locations. Designing such experiments requires wrestling with consequential questions. (Q1) Which specific "types" of a focal species should be introduced to the restoration location? (Q2) How many sources of each type should be used-and should they be mixed together? (Q3) Which specific source populations should be used? (Q4) Which type(s) or population(s) should be introduced into which restoration sites? We recently grappled with these questions when designing an eco-evolutionary experiment with threespine stickleback ( Gasterosteus aculeatus ) introduced into nine small lakes and ponds on the Kenai Peninsula in Alaska that required restoration. After considering the options at length, we decided to use benthic versus limnetic ecotypes (Q1) to create a mixed group of colonists from four source populations of each ecotype (Q2), where ecotypes were identified based on trophic morphology (Q3), and were then introduced into nine restoration lakes scaled by lake size (Q4). We hope that outlining the alternatives and resulting choices will make the rationales clear for future studies leveraging our experiment, while also proving useful for investigators considering similar experiments in the future., Competing Interests: The authors have no conflict of interest to declare., (© 2024 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

This review article deals with the pathways of cellular and global molybdate distribution in plants, especially with a full overview for the model plant Arabidopsis thaliana . In its oxidized state as bioavailable molybdate, molybdenum can be absorbed from the environment. Especially in higher plants, molybdenum is indispensable as part of the molybdenum cofactor (Moco), which is responsible for functionality as a prosthetic group in a variety of essential enzymes like nitrate reductase and sulfite oxidase. Therefore, plants need mechanisms for molybdate import and transport within the organism, which are accomplished via high-affinity molybdate transporter (MOT) localized in different cells and membranes. Two different MOT families were identified. Legumes like Glycine max or Medicago truncatula have an especially increased number of MOT1 family members for supplying their symbionts with molybdate for nitrogenase activity. In Arabidopsis thaliana especially, the complete pathway followed by molybdate through the plant is traceable. Not only the uptake from soil by MOT1.1 and its distribution to leaves, flowers, and seeds by MOT2-family members was identified, but also that inside the cell. the transport trough the cytoplasm and the vacuolar storage mechanisms depending on glutathione were described. Finally, supplying the Moco biosynthesis complex by MOT1.2 and MOT2.1 was demonstrated.

Parasites are ubiquitous, yet their effects on hosts are difficult to quantify and generalize across ecosystems. One promising metric of parasitic impact uses the metabolic theory of ecology (MTE) to calculate energy flux, an estimate of energy lost to parasites. We investigated the feasibility of using metabolic scaling rules to compare the energetic burden of parasitism among individuals. Specifically, we found substantial sensitivity of energy flux estimates to input parameters used in the MTE equation when using available data from a model host-parasite system ( Gasterosteus aculeatus and Schistocephalus solidus ). Using literature values, size data from parasitized wild fish, and a respirometry experiment, we estimate that a single S. solidus tapeworm may extract up to 32% of its stickleback host's baseline metabolic energy requirement, and that parasites in multiple infections may collectively extract up to 46%. The amount of energy siphoned from stickleback to tapeworms is large but did not instigate an increase in respiration rate in the current study. This emphasizes the importance of future work focusing on how parasites influence ecosystem energetics. The approach of using the MTE to calculate energy flux provides great promise as a quantitative foundation for such estimates and provides a more concrete metric of parasite impact on hosts than parasite abundance alone., Competing Interests: The authors declare no conflicts of interest., (© 2023 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

What limits the size of nature's most extreme structures? For weapons like beetle horns, one possibility is a tradeoff associated with mechanical levers: as the output arm of the lever system-the beetle horn-gets longer, it also gets weaker. This "paradox of the weakening combatant" could offset reproductive advantages of additional increases in weapon size. However, in contemporary populations of most heavily weaponed species, males with the longest weapons also tend to be the strongest, presumably because selection drove the evolution of compensatory changes to these lever systems that ameliorated the force reductions of increased weapon size. Therefore, we test for biomechanical limits by reconstructing the stages of weapon evolution, exploring whether initial increases in weapon length first led to reductions in weapon force generation that were later ameliorated through the evolution of mechanisms of mechanical compensation. We describe phylogeographic relationships among populations of a rhinoceros beetle and show that the "pitchfork" shaped head horn likely increased in length independently in the northern and southern radiations of beetles. Both increases in horn length were associated with dramatic reductions to horn lifting strength-compelling evidence for the paradox of the weakening combatant-and these initial reductions to horn strength were later ameliorated in some populations through reductions to horn length or through increases in head height (the input arm for the horn lever system). Our results reveal an exciting geographic mosaic of weapon size, weapon force, and mechanical compensation, shedding light on larger questions pertaining to the evolution of extreme structures., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Molybdenum (Mo) as essential micronutrient for plants, acts as active component of molybdenum cofactor (Moco). Core metabolic processes like nitrate assimilation or abscisic-acid biosynthesis rely on Moco-dependent enzymes. Although a family of molybdate transport proteins (MOT1) is known to date in Arabidopsis, molybdate homeostasis remained unclear. Here we report a second family of molybdate transporters (MOT2) playing key roles in molybdate distribution and usage. KO phenotype-analyses, cellular and organ-specific localization, and connection to Moco-biosynthesis enzymes via protein-protein interaction suggest involvement in cellular import of molybdate in leaves and reproductive organs. Furthermore, we detected a glutathione-molybdate complex, which reveals how vacuolar storage is maintained. A putative Golgi S-adenosyl-methionine transport function was reported recently for the MOT2-family. Here, we propose a moonlighting function, since clear evidence of molybdate transport was found in a yeast-system. Our characterization of the MOT2-family and the detection of a glutathione-molybdate complex unveil the plant-wide way of molybdate., (© 2023. The Author(s).)

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism., (© 2023. The Author(s).)

What are the relative costs and benefits of mounting immune responses? Practitioners of ecoimmunology have grappled with this central question since the field's inception with the main tension being how to make tractable methodological choices that maintain the ecological relevance of induced and measured immune costs. Here, we point out two methodological approaches that we feel are underrepresented in the field, describe risks associated with neglecting these methods, and suggest modern techniques that maximize both the diversity and ecological relevance of collected data. First, it is commonly assumed that frequently used and experimentally convenient immune stimulants will induce ecologically relevant immune responses in study organisms. This can be a dangerous assumption. Even if a stimulant's general immune response properties are well characterized, it is critical to also measure the type and scale of immune responses induced by live pathogens. Second, patterns of immune defenses evolve like other traits, thus a comparative approach is essential to understand what forces shape immune variation. Finally, we describe modern genetic and immunological approaches that will soon become essential tools for ecoimmunologists, and present case studies that exemplify the utility of our recommendations., Competing Interests: Conflicts of interest: None

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.

Animals often adjust their behavior according to social context, but the capacity for such behavioral flexibility can vary among species. Here, we test for interspecific variation in behavioral flexibility by comparing burrowing behavior across three species of deer mice (genus Peromyscus ) with divergent social systems, ranging from promiscuous ( Peromyscus leucopus and Peromyscus maniculatus ) to monogamous ( Peromyscus polionotus ). First, we compared the burrows built by individual mice to those built by pairs of mice in all three species. Although burrow length did not differ in P. leucopus or P. maniculatus , we found that P. polionotus pairs cooperatively constructed burrows that were nearly twice as long as those built by individuals and that opposite-sex pairs dug longer burrows than same-sex pairs. Second, to directly observe cooperative digging behavior in P. polionotus , we designed a burrowing assay in which we could video-record active digging in narrow, transparent enclosures. Using this novel assay, we found, unexpectedly, that neither males nor females spent more time digging with an opposite-sex partner. Rather, we demonstrate that opposite-sex pairs are more socially cohesive and thus more efficient digging partners than same-sex pairs. Together, our study demonstrates how social context can modulate innate behavior and offers insight into how differences in behavioral flexibility may evolve among closely related species., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Evolution Letters published by Wiley Periodicals LLC on behalf of Society for the Study of Evolution (SSE) and European Society for Evolutionary Biology (ESEB).)

Molybdate uptake and molybdenum cofactor (Moco) biosynthesis were investigated in detail in the last few decades. The present study critically reviews our present knowledge about eukaryotic molybdate transporters (MOT) and focuses on the model plant Arabidopsis thaliana , complementing it with new experiments, filling missing gaps, and clarifying contradictory results in the literature. Two molybdate transporters, MOT1.1 and MOT1.2, are known in Arabidopsis , but their importance for sufficient molybdate supply to Moco biosynthesis remains unclear. For a better understanding of their physiological functions in molybdate homeostasis, we studied the impact of mot1.1 and mot1.2 knock-out mutants, including a double knock-out on molybdate uptake and Moco-dependent enzyme activity, MOT localisation, and protein-protein interactions. The outcome illustrates different physiological roles for Moco biosynthesis: MOT1.1 is plasma membrane located and its function lies in the efficient absorption of molybdate from soil and its distribution throughout the plant. However, MOT1.1 is not involved in leaf cell imports of molybdate and has no interaction with proteins of the Moco biosynthesis complex. In contrast, the tonoplast-localised transporter MOT1.2 exports molybdate stored in the vacuole and makes it available for re-localisation during senescence. It also supplies the Moco biosynthesis complex with molybdate by direct interaction with molybdenum insertase Cnx1 for controlled and safe sequestering.

Vertebrate immunity is a complex system consisting of a mix of constitutive and inducible defences. Furthermore, host immunity is subject to selective pressure from a range of parasites and pathogens which can produce variation in these defences across populations. As populations evolve immune responses to parasites, they may adapt via a combination of (1) constitutive differences, (2) shared inducible responses, or (3) divergent inducible responses. Here, we leverage a powerful natural host-parasite model system (Gasterosteus aculeatus and Schistochephalus solidus) to tease apart the relative contributions of these three types of adaptations to among-population divergence in response to parasites. Gene expression analyses revealed limited evidence of significant divergence in constitutive expression of immune defence, and strong signatures of conserved inducible responses to the parasite. Furthermore, our results highlight a handful of immune-related genes which show divergent inducible responses which may contribute disproportionately to functional differences in infection success or failure. In addition to investigating variation in evolutionary adaptation to parasite selection, we also leverage this unique data set to improve understanding of cellular mechanisms underlying a putative resistance phenotype (fibrosis). Combined, our results provide a case study in evolutionary immunology showing that a very small number of genes may contribute to genotype differences in infection response., (© 2021 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

An ever-increasing number of intracellular multi-protein networks have been identified in plant cells. Split-GFP-based protein-protein interaction assays combine the advantages of in vivo interaction studies in a native environment with additional visualization of protein complex localization. Because of their simple protocols, they have become some of the most frequently used methods. However, standard fluorescent proteins present several drawbacks for sophisticated microscopy. With the HaloTag system, these drawbacks can be overcome, as this reporter forms covalent irreversible bonds with synthetic photostable fluorescent ligands. Dyes can be used in adjustable concentrations and are suitable for advanced microscopy methods. Therefore, we have established the Split-HaloTag imaging assay in plants, which is based on the reconstitution of a functional HaloTag protein upon protein-protein interaction and the subsequent covalent binding of an added fluorescent ligand. Its suitability and robustness were demonstrated using a well-characterized interaction as an example of protein-protein interaction at cellular structures: the anchoring of the molybdenum cofactor biosynthesis complex to filamentous actin. In addition, a specific interaction was visualized in a more distinctive manner with subdiffractional polarization microscopy, Airyscan, and structured illumination microscopy to provide examples of sophisticated imaging. Split-GFP and Split-HaloTag can complement one another, as Split-HaloTag represents an alternative option and an addition to the large toolbox of in vivo methods. Therefore, this promising new Split-HaloTag imaging assay provides a unique and sensitive approach for more detailed characterization of protein-protein interactions using specific microscopy techniques, such as 3D imaging, single-molecule tracking, and super-resolution microscopy., (© 2021 The Authors.)

Commensal microbial communities have immense effects on their vertebrate hosts, contributing to a number of physiological functions, as well as host fitness. In particular, host immunity is strongly linked to microbiota composition through poorly understood bi-directional links. Gene expression may be a potential mediator of these links between microbial communities and host function. However, few studies have investigated connections between microbiota composition and expression of host immune genes in complex systems. Here, we leverage a large study of laboratory-raised fish from the species Gasterosteus aculeatus (three-spined stickleback) to document correlations between gene expression and microbiome composition. First, we examined correlations between microbiome alpha diversity and gene expression. Our results demonstrate robust positive associations between microbial alpha diversity and expression of host immune genes. Next, we examined correlations between host gene expression and abundance of microbial taxa. We identified 15 microbial families that were highly correlated with host gene expression. These families were all tightly correlated with host expression of immune genes and processes, falling into one of three categories-those positively correlated, negatively correlated, and neutrally related to immune processes. Furthermore, we highlight several important immune processes that are commonly associated with the abundance of these taxa, including both macrophage and B cell functions. Further functional characterization of microbial taxa will help disentangle the mechanisms of the correlations described here. In sum, our study supports prevailing hypotheses of intimate links between host immunity and gut microbiome composition. IMPORTANCE Here, we document associations between host gene expression and gut microbiome composition in a nonmammalian vertebrate species. We highlight associations between expression of immune genes and both microbiome diversity and abundance of specific microbial taxa. These findings support other findings from model systems which have suggested that gut microbiome composition and host immunity are intimately linked. Furthermore, we demonstrate that these correlations are truly systemic; the gene expression detailed here was collected from an important fish immune organ (the head kidney) that is anatomically distant from the gut. This emphasizes the systemic impact of connections between gut microbiota and host immune function. Our work is a significant advancement in the understanding of immune-microbiome links in nonmodel, natural systems., (Copyright © 2021 Fuess et al.)

Frequency and intensity of wildfire occurrences are dramatically increasing worldwide due to global climate change, having a devastating effect on the entire ecosystem including plants. Moreover, distribution of fire-smoke can influence the natural environment over very long distances, i.e. hundreds of kilometres. Dry plant matter contains 0.1-0.9% (w/w) sulphur, which is mainly released during combustion into the atmosphere as sulphur dioxide (SO 2 ) resulting in local concentrations of up to 3000 nL L -1 . SO 2 is a highly hazardous gas, which enters plants mostly via the stomata. Toxic sulphite is formed inside the leaves due to conversion of SO 2 . Plants as sessile organisms cannot escape from threats, why they evolved an impressive diversity of molecular defence mechanisms. In the present study, two recent wildfires in Germany were evaluated to analyse the effect of SO 2 released into the atmosphere on deciduous trees: the Meppen peat fire in 2018 and the forest fire close to Luebtheen in 2019. Collected leaf material from beech (Fagus sylvatica) and oak (Quercus robur) was examined with respect to detoxification of sulphur surplus due to the exposure to elevated SO 2 . An induced stress reaction in both species was indicated by a 1.5-fold increase in oxidized glutathione. In beech leaves, the enzymatic activities of the sulphite detoxification enzymes sulphite oxidase and apoplastic peroxidases were increased 5-fold and a trend of sulphate accumulation was observed. In contrast, oaks did not regulate these enzymes during smoke exposure, however, the constitutive activity is 10-fold and 3-fold higher than in beech. These results show for the first time sulphite detoxification strategies of trees in situ after natural smoke exposure. Beech and oak trees survived short-term SO 2 fumigation due to exclusion of toxic gases and different oxidative detoxification strategies. Beeches use efficient upregulation of oxidative sulphite detoxification enzymes, while oaks hold a constitutively high enzyme-pool available., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Background: Previous trials suggest lower long-term risk of mortality after invasive rather than non-invasive management of patients with non-ST elevation myocardial infarction (NSTEMI), but the trials excluded very elderly patients. We aimed to estimate the effect of invasive versus non-invasive management within 3 days of peak troponin concentration on the survival of patients aged 80 years or older with NSTEMI., Methods: Routine clinical data for this study were obtained from five collaborating hospitals hosting NIHR Biomedical Research Centres in the UK (all tertiary centres with emergency departments). Eligible patients were 80 years old or older when they underwent troponin measurements and were diagnosed with NSTEMI between 2010 (2008 for University College Hospital) and 2017. Propensity scores (patients' estimated probability of receiving invasive management) based on pretreatment variables were derived using logistic regression; patients with high probabilities of non-invasive or invasive management were excluded. Patients who died within 3 days of peak troponin concentration without receiving invasive management were assigned to the invasive or non-invasive management groups based on their propensity scores, to mitigate immortal time bias. We estimated mortality hazard ratios comparing invasive with non-invasive management, and compared the rate of hospital admissions for heart failure., Findings: Of the 1976 patients with NSTEMI, 101 died within 3 days of their peak troponin concentration and 375 were excluded because of extreme propensity scores. The remaining 1500 patients had a median age of 86 (IQR 82-89) years of whom (845 [56%] received non-invasive management. During median follow-up of 3·0 (IQR 1·2-4·8) years, 613 (41%) patients died. The adjusted cumulative 5-year mortality was 36% in the invasive management group and 55% in the non-invasive management group (adjusted hazard ratio 0·68, 95% CI 0·55-0·84). Invasive management was associated with lower incidence of hospital admissions for heart failure (adjusted rate ratio compared with non-invasive management 0·67, 95% CI 0·48-0·93)., Interpretation: The survival advantage of invasive compared with non-invasive management appears to extend to patients with NSTEMI who are aged 80 years or older., Funding: NIHR Imperial Biomedical Research Centre, as part of the NIHR Health Informatics Collaborative., (Copyright © 2020 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Background Patients presenting with atrial fibrillation (AF) often undergo a blood test to measure troponin, but interpretation of the result is impeded by uncertainty about its clinical importance. We investigated the relationship between troponin level, coronary angiography, and all-cause mortality in real-world patients presenting with AF. Methods and Results We used National Institute of Health Research Health Informatics Collaborative data to identify patients admitted between 2010 and 2017 at 5 tertiary centers in the United Kingdom with a primary diagnosis of AF. Peak troponin results were scaled as multiples of the upper limit of normal. A total of 3121 patients were included in the analysis. Over a median follow-up of 1462 (interquartile range, 929-1975) days, there were 586 deaths (18.8%). The adjusted hazard ratio for mortality associated with a positive troponin (value above upper limit of normal) was 1.20 (95% CI, 1.01-1.43; P <0.05). Higher troponin levels were associated with higher risk of mortality, reaching a maximum hazard ratio of 2.6 (95% CI, 1.9-3.4) at ≈250 multiples of the upper limit of normal. There was an exponential relationship between higher troponin levels and increased odds of coronary angiography. The mortality risk was 36% lower in patients undergoing coronary angiography than in those who did not (adjusted hazard ratio, 0.61; 95% CI, 0.42-0.89; P =0.01). Conclusions Increased troponin was associated with increased risk of mortality in patients presenting with AF. The lower hazard ratio in patients undergoing invasive management raises the possibility that the clinical importance of troponin release in AF may be mediated by coronary artery disease, which may be responsive to revascularization.

To explore landscape genomics at the range limit of an obligate mutualism, we use genotyping-by-sequencing (ddRADseq) to quantify population structure and the effect of host-symbiont interactions between the northernmost fungus-farming leafcutter ant Atta texana and its two main types of cultivated fungus. Genome-wide differentiation between ants associated with either of the two fungal types is of the same order of magnitude as differentiation associated with temperature and precipitation across the ant's entire range, suggesting that specific ant-fungus genome-genome combinations may have been favoured by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. This population-genetic structure was concordant between the ants and one cultivar type (M-fungi, concordant clines) but discordant for the other cultivar type (T-fungi). Discordance in population-genetic structures between ant hosts and a fungal symbiont is surprising because the ant farmers codisperse with their vertically transmitted fungal symbionts. Discordance implies that (a) the fungi disperse also through between-nest horizontal transfer or other unknown mechanisms, and (b) genetic drift and gene flow can differ in magnitude between each partner and between different ant-fungus combinations. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome-genome interactions between host and symbiont influence adaptive genetic differentiation of the host during range evolution in this obligate mutualism., (© 2019 John Wiley & Sons Ltd.)

Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present-day tropical biota. Using genome-wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation of Hypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence in H. lundii across the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post-Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as "isolation by instability," highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas., (© 2019 John Wiley & Sons Ltd.)

Background: Severely injured children and adolescents in clinical practice are rare. For adequate treatment of these patients, detailed knowledge of anatomical and physiological peculiarities, as well as abundant injury patterns, are indispensable. Traumatic brain injuries are known to lead to an unfavorable outcome. In addition, thoracic trauma is regarded as prognostically unfavorable., Objectives: This study depicts epidemiology and injury patterns of severely injured children and adolescents focusing on peculiarities in the severely injured with associated thoracic injuries., Materials and Methods: A retrospective analysis of underaged patients with suspicion of severe injuries who obtained emergency-room treatment in our level-one trauma center during a four-year time period was performed. The data was collected prospectively using the TraumaRegister® of the German Trauma Society as well as an extended house-internal dataset including data of daily clinical routine. The patients were divided into subgroups with (TT) and without (KT) thoracic trauma based on whether a thoracic injury was present or not. For further analysis, four age groups were established., Results: In all, 256 patients younger than 18 years were eligible. Of these, 46 patients revealed thoracic injuries. The mean age of patients with thoracic trauma (12.4 ± 4.9 years) was significantly higher than for patients without thoracic trauma (8.0 ± 5.2 years). In both subgroups, most patients were male (TT: 69.9%, KT: 64.8%). Patients with concomitant thoracic trauma showed a significantly higher injury severity score (ISS) than patients without thoracic trauma (ISS: TT: 26.7 ± 15.8 vs. KT: 8.1 ± 6.8 points). Mortality was higher for TT as well (TT: 6.9% vs. KT: 1.9%). For both groups, traffic accidents were the most common cause of injury. Of patients with thoracic injuries, 52.2% developed at least one complication during their hospital stay (KT: 12.9%)., Conclusions: Thoracic trauma is a relevant factor in children with regard to the severity of total injury and complications. Particular attention should therefore be paid to early diagnosis and treatment.

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.

Parasites can be a major cause of natural selection on hosts, which consequently evolve a variety of strategies to avoid, eliminate, or tolerate infection. When ecologically similar host populations present disparate infection loads, this natural variation can reveal immunological strategies underlying adaptation to infection and population divergence. For instance, the tapeworm Schistocephalus solidus persistently infects 0-80% of threespine stickleback ( Gasterosteus aculeatus ) in lakes on Vancouver Island. To test whether these heterogeneous infection rates result from evolved differences in immunity, we experimentally exposed laboratory-reared fish from ecologically similar high-infection and no-infection populations to controlled doses of Schistocephalus We observed heritable between-population differences in several immune traits: Fish from the naturally uninfected population initiated a stronger granulocyte response to Schistocephalus infection, and their granulocytes constitutively generate threefold more reactive oxygen species in cell culture. Despite these immunological differences, Schistocephalus was equally successful at establishing initial infections in both host populations. However, the no-infection fish dramatically suppressed tapeworm growth relative to high-infection fish, and parasite size was intermediate in F1 hybrid hosts. Our results show that stickleback recently evolved heritable variation in their capacity to suppress helminth growth by two orders of magnitude. Data from many natural populations indicate that growth suppression is widespread but not universal and, when present, is associated with reduced infection prevalence. Host suppression of helminth somatic growth may be an important immune strategy that aids in parasite clearance or in mitigating the fitness costs of persistent infection., Competing Interests: The authors declare no conflict of interest.

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.

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation-by-distance (IBD) and isolation-by-environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake-stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake-stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ∼1.9 km of IBD. Remarkably, of our three lake-stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population., (© 2016 The Author(s). Evolution © 2016 The Society for the Study of Evolution.)

Parasite infections are a product of both ecological processes affecting host-parasite encounter rates and evolutionary dynamics affecting host susceptibility. However, few studies examine natural infection variation from both ecological and evolutionary perspectives. Here, we describe the ecological and evolutionary factors generating variation in infection rates by a tapeworm (Schistocephalus solidus) in a vertebrate host, the threespine stickleback (Gasterosteus aculeatus). To explore ecological aspects of infection, we measured tapeworm prevalence in Canadian stickleback inhabiting two distinct environments: marine and freshwater. Consistent with ecological control of infection, the tapeworm is very rare in marine environments, even though marine fish are highly susceptible. Conversely, commonly infected freshwater stickleback exhibit substantial resistance in controlled laboratory trials, suggesting that high exposure risk overwhelms their recently evolved resistance. We also tested for parasite adaptation to its host by performing transcontinental reciprocal infections, using stickleback and tapeworm populations from Europe and western Canada. More infections occurred in same-continent host-parasite combinations, indicating parasite "local" adaptation, at least on the scale of continents. However, the recently evolved immunity of freshwater hosts applies to both local and foreign parasites. The pattern of adaptation described here is not wholly compatible with either of the common models of host-parasite coevolution (i.e., matching infection or targeted recognition). Instead, we propose a hybrid, eco-evolutionary model to explain the remarkable pattern of global host resistance and local parasite infectivity.

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., (© 2016 John Wiley & Sons Ltd.)

Fetal alcohol spectrum disorders range in severity depending on the amount, timing, and frequency of alcohol exposure. Regardless of severity, sensorimotor defects are commonly reported. Sensorimotor information travels through three tracts of the internal capsule: thalamocortical axons, corticothalamic axons, and corticospinal axons. Here we describe the effects of binge ethanol exposure during the first-trimester equivalent on corticothalamic neurons using Swiss Webster mice. We injected pregnant mice with ethanol (2.9 g/kg, intraperitoneal, followed by 1.45 g/kg, intraperitoneal, 2 h later) on embryonic days (E) 11.5, 12.5, and 13.5. Our paradigm resulted in a mean maternal blood ethanol content of 294.8±15.4 mg/dl on E12.5 and 258.3±22.2 mg/dl on E13.5. Control dams were injected with an equivalent volume of PBS. Bromodeoxyuridine birthdating was carried out on E11.5 to label S-phase neurons. The days of injection were chosen because they are at the onset of neurogenesis and axon extension for corticothalamic, thalamocortical, and corticospinal neurons. Ethanol-exposed pups exhibited no differences compared with controls on day of birth in litter size, body weight, or brain weight. Corticothalamic neurons labeled with bromodeoxyuridine and T-box brain 1 were located in the deep layers of the cortex and did not differ in number in both groups. These results contrast several studies demonstrating alcohol-related differences in these parameters using chronic ethanol exposure paradigms and inbred mouse strains. Therefore, our findings highlight the importance of expanding the mouse strains used to model fetal alcohol spectrum disorder to enhance our understanding of its complex etiology.

Advances in immunoglobulin (Ig) sequencing technology are leading to new perspectives on immune system dynamics. Much research in this nascent field has focused on resolving immune responses to viral infection. However, the dynamics of B-cell diversity in early HIV infection, and in response to anti-retroviral therapy, are still poorly understood. Here, we investigate these dynamics through bulk Ig sequencing of samples collected over 2 years from a group of eight HIV-1 infected patients, five of whom received anti-retroviral therapy during the first half of the study period. We applied previously published methods for visualizing and quantifying B-cell sequence diversity, including the Gini index, and compared their efficacy to alternative measures. While we found significantly greater clonal structure in HIV-infected patients versus healthy controls, within HIV patients, we observed no significant relationships between statistics of B-cell clonal expansion and clinical variables such as viral load and CD4(+) count. Although there are many potential explanations for this, we suggest that important factors include poor sampling resolution and complex B-cell dynamics that are difficult to summarize using simple summary statistics. Importantly, we find a significant association between observed Gini indices and sequencing read depth, and we conclude that more robust analytical methods and a closer integration of experimental and theoretical work is needed to further our understanding of B-cell repertoire diversity during viral infection.

We are writing in response to the population and phylogenomics meeting review by Andrews & Luikart (2014) entitled 'Recent novel approaches for population genomics data analysis'. Restriction-site-associated DNA (RAD) sequencing has become a powerful and useful approach in molecular ecology, with several different published methods now available to molecular ecologists, none of which can be considered the best option in all situations. A&L report that the original RAD protocol of Miller et al. (2007) and Baird et al. (2008) is superior to all other RAD variants because putative PCR duplicates can be identified (see Baxter et al. 2011), thereby reducing the impact of PCR artefacts on allele frequency estimates (Andrews & Luikart 2014). In response, we (i) challenge the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, (ii) present additional biases in RADseq that are at least as important as PCR artefacts in selecting a RAD protocol and (iii) highlight the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants: the original RAD protocol (mbRAD, Miller et al. 2007; Baird et al. 2008), double digest RAD (ddRAD, Peterson et al. 2012), ezRAD (Toonen et al. 2013) and 2bRAD (Wang et al. 2012). With an understanding of the strengths and weaknesses of different RAD protocols, researchers can make a more informed decision when selecting a RAD protocol., (© 2014 John Wiley & Sons Ltd.)

Background: Dynamic changes in Human Immunodeficiency Virus 1 (HIV-1) sequence diversity and divergence are associated with immune control during primary infection and progression to AIDS. Consensus sequencing or single genome amplification sequencing of the HIV-1 envelope (env) gene, in particular the variable (V) regions, is used as a marker for HIV-1 genome diversity, but population diversity is only minimally, or semi-quantitatively sampled using these methods., Results: Here we use second generation deep sequencing to determine inter-and intra-patient sequence heterogeneity and to quantify minor variants in a cohort of individuals either receiving or not receiving antiretroviral treatment following seroconversion; the SPARTAC trial. We show, through a cross-sectional study of sequence diversity of the env V3 in 30 antiretroviral-naive patients during primary infection that considerable population structure diversity exists, with some individuals exhibiting highly constrained plasma virus diversity. Diversity was independent of clinical markers (viral load, time from seroconversion, CD4 cell count) of infection. Serial sampling over 60 weeks of non-treated individuals that define three initially different diversity profiles showed that complex patterns of continuing HIV-1 sequence diversification and divergence could be readily detected. Evidence for minor sequence turnover, emergence of new variants and re-emergence of archived variants could be inferred from this analysis. Analysis of viral divergence over the same time period in patients who received short (12 weeks, ART12) or long course antiretroviral therapy (48 weeks, ART48) and a non-treated control group revealed that ART48 successfully suppressed viral divergence while ART12 did not have a significant effect., Conclusions: Deep sequencing is a sensitive and reliable method for investigating the diversity of the env V3 as an important component of HIV-1 genome diversity. Detailed insights into the complex early intra-patient dynamics of env V3 diversity and divergence were explored in antiretroviral-naïve recent seroconverters. Long course antiretroviral therapy, initiated soon after seroconversion and administered for 48 weeks, restricts HIV-1 divergence significantly. The effect of ART12 and ART48 on clinical markers of HIV infection and progression is currently investigated in the SPARTAC trial.

Relative to morphological traits, we know little about how genetics influence the evolution of complex behavioural differences in nature. It is unclear how the environment influences natural variation in heritable behaviour, and whether complex behavioural differences evolve through few genetic changes, each affecting many aspects of behaviour, or through the accumulation of several genetic changes that, when combined, give rise to behavioural complexity. Here we show that in nature, oldfield mice (Peromyscus polionotus) build complex burrows with long entrance and escape tunnels, and that burrow length is consistent across populations, although burrow depth varies with soil composition. This burrow architecture is in contrast with the small, simple burrows of its sister species, deer mice (P. maniculatus). When investigated under laboratory conditions, both species recapitulate their natural burrowing behaviour. Genetic crosses between the two species reveal that the derived burrows of oldfield mice are dominant and evolved through the addition of multiple genetic changes. In burrows built by first-generation backcross mice, entrance-tunnel length and the presence of an escape tunnel can be uncoupled, suggesting that these traits are modular. Quantitative trait locus analysis also indicates that tunnel length segregates as a complex trait, affected by at least three independent genetic regions, whereas the presence of an escape tunnel is associated with only a single locus. Together, these results suggest that complex behaviours--in this case, a classic 'extended phenotype'--can evolve through multiple genetic changes each affecting distinct behaviour modules.

The ability to efficiently and accurately determine genotypes is a keystone technology in modern genetics, crucial to studies ranging from clinical diagnostics, to genotype-phenotype association, to reconstruction of ancestry and the detection of selection. To date, high capacity, low cost genotyping has been largely achieved via "SNP chip" microarray-based platforms which require substantial prior knowledge of both genome sequence and variability, and once designed are suitable only for those targeted variable nucleotide sites. This method introduces substantial ascertainment bias and inherently precludes detection of rare or population-specific variants, a major source of information for both population history and genotype-phenotype association. Recent developments in reduced-representation genome sequencing experiments on massively parallel sequencers (commonly referred to as RAD-tag or RADseq) have brought direct sequencing to the problem of population genotyping, but increased cost and procedural and analytical complexity have limited their widespread adoption. Here, we describe a complete laboratory protocol, including a custom combinatorial indexing method, and accompanying software tools to facilitate genotyping across large numbers (hundreds or more) of individuals for a range of markers (hundreds to hundreds of thousands). Our method requires no prior genomic knowledge and achieves per-site and per-individual costs below that of current SNP chip technology, while requiring similar hands-on time investment, comparable amounts of input DNA, and downstream analysis times on the order of hours. Finally, we provide empirical results from the application of this method to both genotyping in a laboratory cross and in wild populations. Because of its flexibility, this modified RADseq approach promises to be applicable to a diversity of biological questions in a wide range of organisms.