Measuring the immune system of the three‐spined stickleback – investigating natural variation by quantifying immune expression in the laboratory and the wild

Current understanding of the immune system comes primarily from laboratory‐based studies. There has been substantial interest in examining how it functions in the wild, but studies have been limited by a lack of appropriate assays and study species. The three‐spined stickleback (Gasterosteus aculeatus L.) provides an ideal system in which to advance the study of wild immunology, but requires the development of suitable immune assays. We demonstrate that meaningful variation in the immune response of stickleback can be measured using real‐time PCR to quantify the expression of eight genes, representing the innate response and Th1‐, Th2‐ and Treg‐type adaptive responses. Assays are validated by comparing the immune expression profiles of wild and laboratory‐raised stickleback, and by examining variation across populations on North Uist, Scotland. We also compare the immune response potential of laboratory‐raised individuals from two Icelandic populations by stimulating cells in culture. Immune profiles of wild fish differed from laboratory‐raised fish from the same parental population, with immune expression patterns in the wild converging relative to those in the laboratory. Innate measures differed between wild populations, whilst the adaptive response was associated with variation in age, relative size of fish, reproductive status and S. solidus infection levels. Laboratory‐raised individuals from different populations showed markedly different innate immune response potential. The ability to combine studies in the laboratory and in the wild underlines the potential of this toolkit to advance our understanding of the ecological and evolutionary relevance of immune system variation in a natural setting.     

Quantification of mesocosm fish and amphibian species diversity via environmental DNA metabarcoding

Freshwater fauna are particularly sensitive to environmental change and disturbance. Management agencies frequently use fish and amphibian biodiversity as indicators of ecosystem health and a way to prioritize and assess management strategies. Traditional aquatic bioassessment that relies on capture of organisms via nets, traps and electrofishing gear typically has low detection probabilities for rare species and can injure individuals of protected species. Our objective was to determine whether environmental DNA (eDNA) sampling and metabarcoding analysis can be used to accurately measure species diversity in aquatic assemblages with differing structures. We manipulated the density and relative abundance of eight fish and one amphibian species in replicated 206‐L mesocosms. Environmental DNA was filtered from water samples, and six mitochondrial gene fragments were Illumina‐sequenced to measure species diversity in each mesocosm. Metabarcoding detected all nine species in all treatment replicates. Additionally, we found a modest, but positive relationship between species abundance and sequencing read abundance. Our results illustrate the potential for eDNA sampling and metabarcoding approaches to improve quantification of aquatic species diversity in natural environments and point the way towards using eDNA metabarcoding as an index of macrofaunal species abundance.     

Potential impacts on ecosystem services of land use transitions to second‐generation bioenergy crops in GB

We present the first assessment of the impact of land use change (LUC) to second‐generation (2G) bioenergy crops on ecosystem services (ES) resolved spatially for Great Britain (GB). A systematic approach was used to assess available evidence on the impacts of LUC from arable, semi‐improved grassland or woodland/forest, to 2G bioenergy crops, for which a quantitative ‘threat matrix’ was developed. The threat matrix was used to estimate potential impacts of transitions to either Miscanthus, short‐rotation coppice (SRC, willow and poplar) or short‐rotation forestry (SRF). The ES effects were found to be largely dependent on previous land uses rather than the choice of 2G crop when assessing the technical potential of available biomass with a transition from arable crops resulting in the most positive effect on ES. Combining these data with constraint masks and available land for SRC and Miscanthus (SRF omitted from this stage due to lack of data), south‐west and north‐west England were identified as areas where Miscanthus and SRC could be grown, respectively, with favourable combinations of economic viability, carbon sequestration, high yield and positive ES benefits. This study also suggests that not all prospective planting of Miscanthus and SRC can be allocated to agricultural land class (ALC) ALC 3 and ALC 4 and suitable areas of ALC 5 are only minimally available. Beneficial impacts were found on 146 583 and 71 890 ha when planting Miscanthus or SRC, respectively, under baseline planting conditions rising to 293 247 and 91 318 ha, respectively, under 2020 planting scenarios. The results provide an insight into the interplay between land availability, original land uses, bioenergy crop type and yield in determining overall positive or negative impacts of bioenergy cropping on ecosystems services and go some way towards developing a framework for quantifying wider ES impacts of this important LUC.     

Pooled-Peptide Epitope Mapping Strategies Are Efficient and Highly Sensitive: An Evaluation of Methods for Identifying Human T Cell Epitope Specificities in Large-Scale HIV Vaccine Efficacy Trials

The interferon gamma, enzyme-linked immunospot (IFN-γ ELISpot) assay is widely used to identify viral antigen-specific T cells is frequently employed to quantify T cell responses in HIV vaccine studies. It can be used to define T cell epitope specificities using panels of peptide antigens, but with sample and cost constraints there is a critical need to improve the efficiency of epitope mapping for large and variable pathogens. We evaluated two epitope mapping strategies, based on group testing, for their ability to identify vaccine-induced T-cells from participants in the Step HIV-1 vaccine efficacy trial, and compared the findings to an approach of assaying each peptide individually. The group testing strategies reduced the number of assays required by >7-fold without significantly altering the accuracy of T-cell breadth estimates. Assays of small pools containing 7–30 peptides were highly sensitive and effective at detecting single positive peptides as well as summating responses to multiple peptides. Also, assays with a single 15-mer peptide, containing an identified epitope, did not always elicit a response providing validation that 15-mer peptides are not optimal antigens for detecting CD8+ T cells. Our findings further validate pooling-based epitope mapping strategies, which are critical for characterizing vaccine-induced T-cell responses and more broadly for informing iterative vaccine design. We also show ways to improve their application with computational peptide:MHC binding predictors that can accurately identify the optimal epitope within a 15-mer peptide and within a pool of 15-mer peptides.     

Fibroblast Growth Factor 21 (FGF21) Protects against High Fat Diet Induced Inflammation and Islet Hyperplasia in Pancreas

Fibroblast growth factor 21 (FGF21) is an important endocrine metabolic regulator expressed in multiple tissues including liver and adipose tissue. Although highest levels of expression are in pancreas, little is known about the function of FGF21 in this tissue. In order to understand the physiology of FGF21 in the pancreas, we analyzed its expression and regulation in both acinar and islet tissues. We found that acinar tissue express 20-fold higher levels than that observed in islets. We also observed that pancreatic FGF21 is nutritionally regulated; a marked reduction in FGF21 expression was noted with fasting while obesity is associated with 3–4 fold higher expression. Acinar and islet cells are targets of FGF21, which when systemically administered, leads to phosphorylation of the downstream target ERK 1/2 in about half of acinar cells and a small subset of islet cells. Chronic, systemic FGF21 infusion down-regulates its own expression in the pancreas. Mice lacking FGF21 develop significant islet hyperplasia and periductal lymphocytic inflammation when fed with a high fat obesogenic diet. Inflammatory infiltrates consist of TCRb+ Thy1+ T lymphocytes with increased levels of Foxp3+ regulatory T cells. Increased levels of inflammatory cells were coupled with elevated expression of cytokines such as TNFα, IFNγ and IL1β. We conclude that FGF21 acts to limit islet hyperplasia and may also prevent pancreatic inflammation.     

<Emphasis Type="Italic">In silico</Emphasis> coronary wave intensity analysis: application of an integrated one-dimensional and poromechanical model of cardiac perfusion

Coronary wave intensity analysis (cWIA) is a diagnostic technique based on invasive measurement of coronary pressure and velocity waveforms. The theory of WIA allows the forward- and backward-propagating coronary waves to be separated and attributed to their origin and timing, thus serving as a sensitive and specific cardiac functional indicator. In recent years, an increasing number of clinical studies have begun to establish associations between changes in specific waves and various diseases of myocardium and perfusion. These studies are, however, currently confined to a trial-and-error approach and are subject to technological limitations which may confound accurate interpretations. In this work, we have developed a biophysically based cardiac perfusion model which incorporates full ventricular–aortic–coronary coupling. This was achieved by integrating our previous work on one-dimensional modelling of vascular flow and poroelastic perfusion within an active myocardial mechanics framework. Extensive parameterisation was performed, yielding a close agreement with physiological levels of global coronary and myocardial function as well as experimentally observed cumulative wave intensity magnitudes. Results indicate a strong dependence of the backward suction wave on QRS duration and vascular resistance, the forward pushing wave on the rate of myocyte tension development, and the late forward pushing wave on the aortic valve dynamics. These findings are not only consistent with experimental observations, but offer a greater specificity to the wave-originating mechanisms, thus demonstrating the value of the integrated model as a tool for clinical investigation.     

Four putative entomopathogenic fungi of armoured scale insects on <Emphasis Type="Italic">Citrus</Emphasis> in Australia

This study led to the discovery of four putative entomopathogenic fungi of armoured scale insects on citrus trees in coastal New South Wales. Two of these species belong in Podonectria as P. coccicola (Ellis & Everh.) Petch (syn. Tetracrium coccicola (Höhn.) Ellis & Everh.) and P. novae-zelandiae Dingley. Members of this genus are grown in culture for the first time. Formerly placed in the Pleosporales, Tubeufiaceae, or more recently in the Tubeufiales, these species are herein placed in the new family Podonectriaceae fam. nov., Pleosporales. Another species is placed in the Hypocreales, Bionectriaceae as Clonostachys coccicola (J.A. Stev.) H.T. Dao comb. nov. (basionym Tubercularia coccicola J.A. Stev., syn. Nectria tuberculariae Petch). The fourth species is Myriangium citri Henn. (Myriangiales, Myriangiaceae). Each fungal species is characterized and the phylogenetic placement confirmed by molecular analyses of the ITS and 28 s rDNA regions. In addition, their biology is noted, including location of the fungi within tree canopies.     

Consideration of multiple load cases is critical in modelling orthotropic bone adaptation in the femur

Functional adaptation of the femur has been investigated in several studies by embedding bone remodelling algorithms in finite element (FE) models, with simplifications often made to the representation of bone’s material symmetry and mechanical environment. An orthotropic strain-driven adaptation algorithm is proposed in order to predict the femur’s volumetric material property distribution and directionality of its internal structures within a continuum. The algorithm was applied to a FE model of the femur, with muscles, ligaments and joints included explicitly. Multiple load cases representing distinct frames of two activities of daily living (walking and stair climbing) were considered. It is hypothesised that low shear moduli occur in areas of bone that are simply loaded and high shear moduli in areas subjected to complex loading conditions. In addition, it is investigated whether material properties of different femoral regions are stimulated by different activities. The loading and boundary conditions were considered to provide a physiological mechanical environment. The resulting volumetric material property distribution and directionalities agreed with ex vivo imaging data for the whole femur. Regions where non-orthogonal trabecular crossing has been documented coincided with higher values of predicted shear moduli. The topological influence of the different activities modelled was analysed. The influence of stair climbing on the properties of the femoral neck region is highlighted. It is recommended that multiple load cases should be considered when modelling bone adaptation. The orthotropic model of the complete femur is released with this study.     

The dynamics of coupled populations subject to control

The dynamics of coupled populations have mostly been studied in the context of metapopulation viability with application to, for example, species at risk. However, when considering pests and pathogens, eradication, not persistence, is often the end goal. Humans may intervene to control nuisance populations, resulting in reciprocal interactions between the human and natural systems that can lead to unexpected dynamics. The incidence of these human-natural couplings has been increasing, hastening the need to better understand the emergent properties of such systems in order to predict and manage outbreaks of pests and pathogens. For example, the success of the growing aquaculture industry depends on our ability to manage pathogens and maintain a healthy environment for farmed and wild fish. We developed a model for the dynamics of connected populations subject to control, motivated by sea louse parasites that can disperse among salmon farms. The model includes exponential population growth with a forced decline when populations reach a threshold, representing control interventions. Coupling two populations with equal growth rates resulted in phase locking or synchrony in their dynamics. Populations with different growth rates had different periods of oscillation, leading to quasiperiodic dynamics when coupled. Adding small amounts of stochasticity destabilized quasiperiodic cycles to chaos, while stochasticity was damped for periodic or stable dynamics. Our analysis suggests that strict treatment thresholds, although well intended, can complicate parasite dynamics and hinder control efforts. Synchronizing populations via coordinated management among farms leads to more effective control that is required less frequently. Our model is simple and generally applicable to other systems where dispersal affects the management of pests and pathogens.