Relative metabolomic flux: Predicting changes in the environmental metabolome from the metagenome

Background: The world's oceans are home to a diverse array of microbial life whose metabolic activity helps to drive the earth's biogeochemical cycles. Metagenomic analysis has revolutionized our access to these communities, providing a system-scale perspective of microbial community interactions. However, while metagenome sequencing can provide useful estimates of the relative change in abundance of specific genes and taxa between environments or over time, this does not investigate the relative changes in the production or consumption of different metabolites.
Results: We propose a methodology, Predicted Relative Metabolic Turnover (PRMT) that defines and enables exploration of metabolite-space inferred from the metagenome. Our analysis of metagenomic data from a time-series study in the Western English Channel demonstrated considerable correlations between predicted relative metabolic turnover and seasonal changes in abundance of measured environmental parameters as well as with observed seasonal changes in bacterial population structure.
Conclusions: The PRMT method was successfully applied to metagenomic data to explore the Western English Channel microbial metabalome to generate specific, biologically testable hypotheses. Generated hypotheses linked organic phosphate utilization to Gammaproteobactaria, Plantcomycetes, and Betaproteobacteria, chitin degradation to Actinomycetes, and potential small molecule biosynthesis pathways for Lentisphaerae, Chlamydiae, and Crenarchaeota. The PRMT method can be applied as a general tool for the analysis of additional metagenomic or transcriptomic datasets.

Marine reserve designation, trophic cascades and altered community dynamics

Marine ecosystems and their associated populations are increasingly at risk from the cumulative impacts of many anthropogenic threats that increase the likelihood of species extinction and altered community dynamics. In response, marine reserves can be used to protect exploited species and conserve biodiversity. The increased abundance of predatory species in marine reserves may cause indirect effects along chains of multi-trophic interactions. These trophic cascades can arise through direct predation, density-mediated indirect interactions (DMIIs), or indirect behavioural effects, termed trait-mediated indirect interactions (TMIIs). The extent of algal cover and the abundance of 4 primary consumers were determined in Lough Hyne, which was designated Europe's first marine nature reserve in 1981. The primary consumers were the sea urchin Paracentrotus lividus, the topshell Gibbula cineraria, the oyster Anomia ephippium, and the scallop Chlamys varia. The abundances of 3 starfish species (Marthasterias glacialis, Asterias rubens, and Asterina gibbosa) were also determined, as were 2 potential crustacean predators, Necora puber and Carcinus maenas. These data were compared with historical data from a 1962 (prey) and a 1963 (predator) survey to determine the nature of community interactions over adjacent trophic levels. The present study reveals a breakdown in population structure of the 4 surveyed prey species. Marine reserve designation has led to an increase in predatory crabs and M. glacialis, a subsequent decrease in primary consumers, especially the herbivore P. lividus, and an increase in macroalgal cover which is indicative of a trophic cascade. The study shows that establishing a Marine Reserve does not guarantee that conservation benefits will be distributed equally.

Invited review - genetic and epigenetic factors influencing CKD

Chronic kidney disease (CKD) has become a serious public health problem because of its associated morbidity, premature mortality and attendant healthcare costs. The rising number of persons with CKD is linked with ageing population structure and an increased prevalence of diabetes, hypertension and obesity. There is an inherited risk associated with developing CKD as evidenced by familial clustering and differing prevalence rates across ethnic groups. Earlier studies to determine the inherited risk factors for CKD rarely identified genetic variants that were robustly replicated. However, improvements in genotyping technologies and analytical methods are now helping to identify promising genetic loci aided by international collaboration and multi-consortia efforts. More recently, epigenetic modifications have been proposed to play a role in both the inherited susceptibility to CKD and, importantly, to explain how the environment dynamically interacts with the genome to alter an individual's disease risk. Genome-wide, epigenome-wide and whole transcriptome studies have been performed and optimal approaches for integrative analysis are being developed. This review summarises recent research and the current status of genetic and epigenetic risk factors influencing CKD using population-based information.

Pseudomonas aeruginosa exhibits frequent recombination, but only a limited association between genotype and ecological setting

Pseudomonas aeruginosa is an opportunistic pathogen and an important cause of infection, particularly amongst cystic fibrosis (CF) patients. While specific strains capable of patient-to-patient transmission are known, many infections appear to be caused by unique and unrelated strains. There is a need to understand the relationship between strains capable of colonising the CF lung and the broader set of P. aeruginosa isolates found in natural environments. Here we report the results of a multilocus sequence typing (MLST)-based study designed to understand the genetic diversity and population structure of an extensive regional sample of P. aeruginosa isolates from South East Queensland, Australia. The analysis is based on 501 P. aeruginosa isolates obtained from environmental, animal and human (CF and non-CF) sources with particular emphasis on isolates from the Lower Brisbane River and isolates from CF patients obtained from the same geographical region. Overall, MLST identified 274 different sequence types, of which 53 were shared between one or more ecological settings. Our analysis revealed a limited association between genotype and environment and evidence of frequent recombination. We also found that genetic diversity of P. aeruginosa in Queensland, Australia was indistinguishable from that of the global P. aeruginosa population. Several CF strains were encountered frequently in multiple ecological settings; however, the most frequently encountered CF strains were confined to CF patients. Overall, our data confirm a non-clonal epidemic structure and indicate that most CF strains are a random sample of the broader P. aeruginosa population. The increased abundance of some CF strains in different geographical regions is a likely product of chance colonisation events followed by adaptation to the CF lung and horizontal transmission among patients.

Characterisation of a novel panel of polymorphic microsatellite loci for the liver fluke, Fasciola hepatica, using a next generation sequencing approach

The liver fluke, Fasciola hepatica is an economically important pathogen of sheep and cattle and has been described by the WHO as a re-emerging zoonosis. Control is heavily reliant on the use of drugs, particularly triclabendazole and as a result resistance has now emerged. The population structure of F. hepatica is not well known, yet it can impact on host-parasite interactions and parasite control with drugs, particularly regarding the spread of triclabendazole resistance. We have identified 2448 potential microsatellites from 83Mb of F. hepatica genome sequence using msatfinder. Thirty-five loci were developed and optimised for microsatellite PCR, resulting in a panel of 15 polymorphic loci, with a range of three to 15 alleles. This panel was validated on genomic DNA from 46 adult F. hepatica; 38 liver flukes sourced from a Northwest abattoir, UK and 8 liver flukes from an established isolate (Shrewsbury; Ridgeway Research). Evidence for null alleles was found at four loci (Fh_1, Fh_8, Fh_13 and Fh_14), which showed markedly higher levels of homozygosity than the remaining 11 loci. Of the 38 liver flukes isolated from cattle livers (n=10) at the abattoir, 37 genotypes were identified. Using a multiplex approach all 15 loci could be amplified from several life cycle stages that typically yield low amounts of DNA, including metacercariae, the infective life cycle stage present on pasture, highlighting the utility of this multiplex microsatellite panel. This study reports the largest panel of microsatellite markers available to date for population studies of F. hepatica and the first multiplex panel of microsatellite markers that can be used for several life cycle stages.

Using particle dispersal models to assist in the conservation and recovery of the overexploited native oyster (Ostrea edulis) in an enclosed sea lough

Oyster populations around the world have seen catastrophic decline which has been largely attributed to overexploitation, disease and pollution. While considerable effort and resources have been implemented into restoring these important environmental engineers, the success of oyster populations is often limited by poor understanding of site-specific dispersal patterns of propagules. Water-borne transport is a key factor controlling or regulating the dispersal of the larval stage of benthic marine invertebrates which have limited mobility. The distribution of the native oyster Ostrea edulis in Strangford Lough, Northern Ireland, together with their densities and population structure at subtidal and intertidal sites has been documented at irregular intervals between 1997 and 2013. This paper revisits this historical data and considers whether different prevailing environmental conditions can be used to explain the distribution, densities and population structure of O. edulis in Strangford Lough. The approach adopted involved comparing predictive 2D hydrodynamic models coupled with particle tracking to simulate the dispersal of oyster larvae with historical and recent field records of the distribution of both subtidal and intertidal, populations since 1995. Results from the models support the hypothesis that commercial stocks of O. edulis introduced into Strangford Lough in the 1990s resulted in the re-establishment of wild populations of oysters in the Northern Basin which in turn provided a potential source of propagules for subtidal populations. These results highlight that strategic site selection (while inadvertent in the case of the introduced population in 1995) for the re-introduction of important shellfish species can significantly accelerate their recovery and restoration.

Outlier SNP markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome-wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large- and fine-scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F(CT) = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F(CT) range 0.275-0.705) and fine-scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.

Gene-associated markers can assign origin in a weakly structured fish, Atlantic herring

Regulations on the exploitation of populations of commercially important fish species and the ensuing consumer interest in sustainable products have increased the need to accurately identify the population of origin of fish and fish products. Although genomics-based tools have proven highly useful, there are relatively few examples in marine fish displaying accurate origin assignment. We synthesize data for 156 single-nucleotide polymorphisms typed in 1039 herring, Clupea harengus L., spanning the Northeast Atlantic to develop a tool that allows assignment of individual herring to their regional origin. We show the method's suitability to address specific biological questions, as well as management applications. We analyse temporally replicated collections from two areas, the Skagerrak (n = 81, 84, 66) and the western Baltic (n = 52, 52). Both areas harbour heavily fished mixed-origin stocks, complicating management issues. We report novel genetic evidence that herring from the Baltic Sea contribute to catches in the North Sea, and find support that western Baltic feeding aggregations mainly constitute herring from the western Baltic with contributions from the Eastern Baltic. Our study describes a general approach and outlines a database allowing individual assignment and traceability of herring across a large part of its East Atlantic distribution.

Twenty-two novel microsatellite loci for lumpfish (Cyclopterus lumpus)

Lumpfish (Cyclopterus lumpus L. 1758) are widely distributed on both sides of the North Atlantic. They are a commercially important species, but stock size estimates have declined since the mid-1980s in Canada, Norway and Iceland. Little is known about the biology of this species, in particular the breeding migrations and population structure which are fundamental for effective management. This paper describes the development and characterization of twenty-two polymorphic microsatellite loci using next generation sequencing. The number of alleles per locus ranged from 3 to 27 in two geographically distant North Atlantic populations, with observed and expected heterozygosities ranging between 0. 0625-0. 979 and 0. 0618-0. 946, respectively. These loci are an important resource that will allow assessment of the population genetic structure of this species, and contribute to its appropriate management.

New microsatellite loci for the longnose velvet dogfish Centroselachus crepidater (Squaliformes: Somniosidae) and other deep sea sharks

The continuing over-exploitation of traditional coastal stocks has resulted in the shift of commercial fishing towards deep-sea ecosystems in many parts of the world. The effects on target and non-target species have been dramatic; particularly for the deep-sea sharks. With the aim of providing tools that will allow the assessment of population genetic structure of Centroselachus crepidater, novel microsatellite loci have been developed for this deep-sea elasmobranch. Seven of these markers showed between 3 and 7 alleles per locus in two North Atlantic populations, with observed and expected heterozygosities between 0.18-0.95 and 0.25-0.82, respectively. Additionally, ten loci cross-amplify in other Elasmobranch species.

Novel microsatellite loci for a deep sea fish (Macrourus berglax) and their amplification in other grenadiers (Gadiformes: Macrouridae)

Over-exploitation of traditional coastal stocks and a rising demand for seafood have resulted in the shift of commercial fishing towards less-known, deep-sea species in many parts of the world. Yet, the lack of knowledge of the biology, ecology and life-history of these species represents a serious impediment for establishing sound stock management plans. With the aim of providing tools that will allow assessment of the population genetic structure of Macrourus berglax, we have isolated and characterised a suite of novel microsatellite loci for this deep sea grenadier. Eight of these markers showed between 4 and 11 alleles per locus in two distant North Atlantic populations, with observed and expected heterozygosities between 0.17-0.83 and 0.35-0.87, respectively. Importantly, eight of these loci also cross-amplify in other Macrourid species. 

Directional genetic differentiation and relative migration

Understanding the population structure and patterns of gene flow within species is of fundamental importance to the study of evolution. In the fields of population and evolutionary genetics, measures of genetic differentiation are commonly used to gather this information. One potential caveat is that these measures assume gene flow to be symmetric. However, asymmetric gene flow is common in nature, especially in systems driven by physical processes such as wind or water currents. As information about levels of asymmetric gene flow among populations is essential for the correct interpretation of the distribution of contemporary genetic diversity within species, this should not be overlooked. To obtain information on asymmetric migration patterns from genetic data, complex models based on maximum-likelihood or Bayesian approaches generally need to be employed, often at great computational cost. Here, a new simpler and more efficient approach for understanding gene flow patterns is presented. This approach allows the estimation of directional components of genetic divergence between pairs of populations at low computational effort, using any of the classical or modern measures of genetic differentiation. These directional measures of genetic differentiation can further be used to calculate directional relative migration and to detect asymmetries in gene flow patterns. This can be done in a user-friendly web application called divMigrate-online introduced in this study. Using simulated data sets with known gene flow regimes, we demonstrate that the method is capable of resolving complex migration patterns under a range of study designs.

The factor structure of the revised Illness Perceptions Questionnaire in a population of oesophageal cancer survivors

Objectives
To determine whether the proposed 7-factor structure of the Illness Perception Questionnaire-Revised (Timeline Acute/Chronic, Timeline Cyclical, Consequences, Personal Control, Treatment Control, Illness Coherence and Emotional Representations) is appropriate among a population of oesophageal cancer survivors.
Methods
Everyone registered with the Oesophageal Patients’ Association in the UK (n=2185) was mailed a questionnaire booklet which included the Illness Perception Questionnaire-Revised. Responses from 587 oesophageal cancer survivors (27%) were subjected to a confirmatory factor analysis.
Results
The proposed 7 factor structure provided a reasonable fit of the data. Modification indices suggested that a significantly better fit could be provided if one of the items on the Timeline Acute/Chronic factor loaded on the Treatment Control factor and an error covariance was added between 2 other items on the Timeline Acute/Chronic factor.
Conclusions
The model fit for the 7 factor structure proposed by Moss-Morris et al. (2002) was found to be adequate in our study. However, the structure of the timeline acute/chronic factor needs to be considered, particularly when the IPQ-R is to be used among older people with a potentially life-threatening illness or those receiving palliative care.