Application of microarrays (phylochips) for analysis of community diversity by species identification

Abstract Molecular probe-based methods (Fluorescent in-situ hybridisation or FISH, Next Generation Sequencing or NGS) have proved successful in improving both the efficiency and accuracy of the identification of microorganisms, especially those that lack distinct morphological features, such as picoplankton. However, FISH methods have the major drawback that they can only identify one or just a few species at a time because of the reduced number of available fluorochromes that can be added to the probe. Although the length of sequence that can be obtained is continually improving, NGS still requires a great deal of handling time, its analysis time is still months and with a PCR step it will always be sensitive to natural enzyme inhibitors. With the use of DNA microarrays, it is possible to identify large numbers of taxa on a single-glass slide, the so-called phylochip, which can be semi-quantitative. This review details the major steps in probe design, design and production of a phylochip and validation of the array. Finally, major microarray studies in the phytoplankton community are reviewed to demonstrate the scope of the method.

Application of microarrays (phylochips) for analysis of community diversity by species identification

Abstract Molecular probe-based methods (Fluorescent in-situ hybridisation or FISH, Next Generation Sequencing or NGS) have proved successful in improving both the efficiency and accuracy of the identification of microorganisms, especially those that lack distinct morphological features, such as picoplankton. However, FISH methods have the major drawback that they can only identify one or just a few species at a time because of the reduced number of available fluorochromes that can be added to the probe. Although the length of sequence that can be obtained is continually improving, NGS still requires a great deal of handling time, its analysis time is still months and with a PCR step it will always be sensitive to natural enzyme inhibitors. With the use of DNA microarrays, it is possible to identify large numbers of taxa on a single-glass slide, the so-called phylochip, which can be semi-quantitative. This review details the major steps in probe design, design and production of a phylochip and validation of the array. Finally, major microarray studies in the phytoplankton community are reviewed to demonstrate the scope of the method.

Superinfection exclusion and the long-term survival of honey bees in Varroa-infested colonies

Over the past 50 years, many millions of European honey bee (Apis mellifera) colonies have died as the ectoparasitic mite, Varroa destructor, has spread around the world. Subsequent studies have indicated that the mite’s association with a group of RNA viral pathogens (Deformed Wing Virus, DWV) correlates with colony death. Here, we propose a phenomenon known as superinfection exclusion that provides an explanation of how certain A. mellifera populations have survived, despite Varroa infestation and high DWV loads. Next-generation sequencing has shown that a non-lethal DWV variant ‘type B’ has become established in these colonies and that the lethal ‘type A’ DWV variant fails to persist in the bee population. We propose that this novel stable host-pathogen relationship prevents the accumulation of lethal variants, suggesting that this interaction could be exploited for the development of an effective treatment that minimises colony losses in the future.

Superinfection exclusion and the long-term survival of honey bees in Varroa-infested colonies

Over the past 50 years, many millions of European honey bee (Apis mellifera) colonies have died as the ectoparasitic mite, Varroa destructor, has spread around the world. Subsequent studies have indicated that the mite’s association with a group of RNA viral pathogens (Deformed Wing Virus, DWV) correlates with colony death. Here, we propose a phenomenon known as superinfection exclusion that provides an explanation of how certain A. mellifera populations have survived, despite Varroa infestation and high DWV loads. Next-generation sequencing has shown that a non-lethal DWV variant ‘type B’ has become established in these colonies and that the lethal ‘type A’ DWV variant fails to persist in the bee population. We propose that this novel stable host-pathogen relationship prevents the accumulation of lethal variants, suggesting that this interaction could be exploited for the development of an effective treatment that minimises colony losses in the future.

Diversity in a honey bee pathogen: first report of a third master variant of the Deformed Wing Virus quasispecies

Treatment of emerging RNA viruses is hampered by the high mutation and replication rates that enable these viruses to operate as a quasispecies. Declining honey bee populations have been attributed to the ectoparasitic mite Varroa destructor and its affiliation with Deformed Wing Virus (DWV). In the current study we use next-generation sequencing to investigate the DWV quasispecies in an apiary known to suffer from overwintering colony losses. We show that the DWV species complex is made up of three master variants. Our results indicate that a new DWV Type C variant is distinct from the previously described types A and B, but together they form a distinct clade compared with other members of the Iflaviridae. The molecular clock estimation predicts that Type C diverged from the other variants ~319 years ago. The discovery of a new master variant of DWV has important implications for the positive identification of the true pathogen within global honey bee populations.

Diversity in a honey bee pathogen: first report of a third master variant of the Deformed Wing Virus quasispecies

Treatment of emerging RNA viruses is hampered by the high mutation and replication rates that enable these viruses to operate as a quasispecies. Declining honey bee populations have been attributed to the ectoparasitic mite Varroa destructor and its affiliation with Deformed Wing Virus (DWV). In the current study we use next-generation sequencing to investigate the DWV quasispecies in an apiary known to suffer from overwintering colony losses. We show that the DWV species complex is made up of three master variants. Our results indicate that a new DWV Type C variant is distinct from the previously described types A and B, but together they form a distinct clade compared with other members of the Iflaviridae. The molecular clock estimation predicts that Type C diverged from the other variants ~319 years ago. The discovery of a new master variant of DWV has important implications for the positive identification of the true pathogen within global honey bee populations.

A metagenomic comparison of endemic viruses from broiler chickens with runting stunting syndrome and from normal birds

Runting-stunting syndrome (RSS) in broiler chickens is an enteric disease that causes significant economic losses to poultry producers worldwide due to elevated feed conversion ratios, decreased body weight during growth, and excessive culling. Of specific interest are the viral agents associated with RSS which have been difficult to fully characterise to date. Past research into the aetiology of RSS has implicated a wide variety of RNA and DNA viruses however, to date, no individual virus has been identified as the main agent of RSS and the current opinion is that it may be caused by a community of viruses, collectively known as the virome. This paper attempts to characterise the viral pathogens associated with 2 – 3 week old RSS-affected and unaffected broiler chickens using next-generation sequencing and comparative metagenomics. Analysis of the viromes identified a total of 20 DNA & RNA viral families, along with 2 unidentified categories, comprised of 31 distinct viral genera and 7 unclassified genera. The most abundant viral families identified in this study were the Astroviridae, Caliciviridae, Picornaviridae, Parvoviridae, Coronaviridae, Siphoviridae, and Myoviridae. This study has identified historically significant viruses associated with the disease such as chicken astrovirus, avian nephritis virus, chicken parvovirus, and chicken calicivirus along with relatively novel viruses such as chicken megrivirus and sicinivirus 1 and will help expand the knowledge related to enteric disease in broiler chickens, provide insights into the viral constituents of a healthy avian gut, and identify a variety of enteric viruses and viral communities appropriate for further study.

Prior knowledge transfer across transcriptional data sets and technologies using compositional statistics yields new mislabelled ovarian cell line

Here, we describe gene expression compositional assignment (GECA), a powerful, yet simple method based on compositional statistics that can validate the transfer of prior knowledge, such as gene lists, into independent data sets, platforms and technologies. Transcriptional profiling has been used to derive gene lists that stratify patients into prognostic molecular subgroups and assess biomarker performance in the pre-clinical setting. Archived public data sets are an invaluable resource for subsequent in silico validation, though their use can lead to data integration issues. We show that GECA can be used without the need for normalising expression levels between data sets and can outperform rank-based correlation methods. To validate GECA, we demonstrate its success in the cross-platform transfer of gene lists in different domains including: bladder cancer staging, tumour site of origin and mislabelled cell lines. We also show its effectiveness in transferring an epithelial ovarian cancer prognostic gene signature across technologies, from a microarray to a next-generation sequencing setting. In a final case study, we predict the tumour site of origin and histopathology of epithelial ovarian cancer cell lines. In particular, we identify and validate the commonly-used cell line OVCAR-5 as non-ovarian, being gastrointestinal in origin. GECA is available as an open-source R package.

Efficient Genotyping of KRAS Mutant Non-Small Cell Lung Cancer Using a Multiplexed Droplet Digital PCR Approach.

Droplet digital PCR (ddPCR) can be used to detect low frequency mutations in oncogene-driven lung cancer. The range of KRAS point mutations observed in NSCLC necessitates a multiplex approach to efficient mutation detection in circulating DNA. Here we report the design and optimisation of three discriminatory ddPCR multiplex assays investigating nine different KRAS mutations using PrimePCR™ ddPCR™ Mutation Assays and the Bio-Rad QX100 system. Together these mutations account for 95% of the nucleotide changes found in KRAS in human cancer. Multiplex reactions were optimised on genomic DNA extracted from KRAS mutant cell lines and tested on DNA extracted from fixed tumour tissue from a cohort of lung cancer patients without prior knowledge of the specific KRAS genotype. The multiplex ddPCR assays had a limit of detection of better than 1 mutant KRAS molecule in 2,000 wild-type KRAS molecules, which compared favourably with a limit of detection of 1 in 50 for next generation sequencing and 1 in 10 for Sanger sequencing. Multiplex ddPCR assays thus provide a highly efficient methodology to identify KRAS mutations in lung adenocarcinoma.

Challenges in molecular testing in non-small-cell lung cancer patients with advanced disease

Lung cancer diagnostics have progressed greatly in the previous decade. Development of molecular testing to identify an increasing number of potentially clinically actionable genetic variants, using smaller samples obtained via minimally invasive techniques, is a huge challenge. Tumour heterogeneity and cancer evolution in response to therapy means that repeat biopsies or circulating biomarkers are likely to be increasingly useful to adapt treatment as resistance develops. We highlight some of the current challenges faced in clinical practice for molecular testing of EGFR, ALK, and new biomarkers such as PDL1. Implementation of next generation sequencing platforms for molecular diagnostics in non-small-cell lung cancer is increasingly common, allowing testing of multiple genetic variants from a single sample. The use of next generation sequencing to recruit for molecularly stratified clinical trials is discussed in the context of the UK Stratified Medicine Programme and The UK National Lung Matrix Trial.

Microbial phylogenetic and functional responses within acidified wastewater communities exhibiting enhanced phosphate uptake

Acid stimulated accumulation of insoluble phosphorus within microbial cells is highly beneficial to wastewater treatment but remains largely unexplored. Using single cell analyses and next generation sequencing, the response of active polyphosphate accumulating microbial communities under conditions of enhanced phosphorus uptake under both acidic and aerobic conditions was characterised. Phosphorus accumulation activities were highest under acidic conditions (pH 5.5 > 8.5), where a significant positive effect on bioaccumulation was observed at pH 5.5 when compared to pH 8.5. In contrast to the Betaproteobacteria and Actinobacteria dominated enhanced biological phosphorus removal process, the functionally active polyP accumulators at pH 5.5 belonged to the Gammaproteobacteria, with key accumulators identified as members of the families Aeromonadaceae and Enterobacteriaceae. This study demonstrated a significant enrichment of key polyphosphate kinase and exopolyphosphatase genes within the community metagenome after acidification, concomitant with an increase in P accumulation kinetics.

Genomics and virulence factors of Fusobacterium necrophorum

Fusobacterium necrophorum, a Gram negative, anaerobic bacterium, is a common cause of acute pharyngitis and tonsillitis and a rare cause of more severe infections of the head and neck. At the beginning of the project, there was no available genome sequence for F. necrophorum. The aim of this project was to sequence the F. necrophorum genome and identify and study its putative virulence factors contained using in silico and in vitro analysis. Type strains JCM 3718 and JCM 3724,F. necrophorum subspecies necrophorum (Fnn) and funduliforme (Fnf), respectively, and strain ARU 01 (Fnf), isolated from a patient with LS, were commercially sequenced by Roche 454 GS-FLX+ next generation sequencing and assembled into contigs using Roche GS Assembler. Sequence data was annotated semi-automatically, using the xBASE pipeline, BLASTp and Pfam. The F. necrophorum genome was determined to be approximately 2.1 – 2.3 Mb in size, with an estimated 1,950 ORFs and includes genes for a leukotoxin, ecotin, haemolysin, haemagglutinin, haemin receptor, adhesin and type Vb and Vc secretion systems. The prevalence of the leukotoxin gene was investigated in strains JCM 3718, JCM 3724 and ARU 01, as well as a clinical collection of 25 Fnf strains, identified using biochemical and molecular tests. The leukotoxin operon was found to be universal within the strain collection by PCR. HL-60 cells subjected to aliquots of concentrated high molecular weight culture supernatant, predicted to contain the secreted leukotoxins of strains JCM 3718, JCM 3724 and ARU 01, were killed in a dose-dependent manner. The cytotoxic effect of the leukotoxin against human donor white blood cells was also tested to validate the HL-60 assay. The differences in the results between the two assays were not statistically significant. Ecotin, a serine protease inhibitor, was found to be present in 100 % of the strain collection and had a highly conserved sequence with primary and secondary binding sites exposed on opposing sides of the protein. During enzyme inhibition studies, a purified recombinant F. necrophorum ecotin protein inhibited human neutrophil elastase, a protease that degrades bacteria at inflammation sites, and human plasma kallikrein, a component of the host clotting cascade. The recombinant ecotin also prolonged human plasma clotting times by up to 7-fold for the extrinsic pathway, and up to 40-fold for the intrinsic pathway. The genome sequence data provides important information about F. necrophorum type strains and enables comparative study between strains and subspecies. Results from the leukotoxin and ecotin assays can be used to build up an understanding of how the organism behaves during infection.

Role of fibroblast growth factors and their receptors in prostate cancer

Prostate cancer (PCa) is the most common non-cutaneous malignant disease among males in the developed countries. Radical prostatectomy (RP) is an effective therapy for most PCa patients with localized or locally invaded tumors but in some cases the cancer recurs after RP. PCa is a heterogeneous disease, which is regulated by many factors, such as androgen receptor (AR), estrogen receptors and  (ER and ER), fibroblast growth factors (FGFs) and their receptors (FGFRs). In this study, the role of ERβ, FGF8, FGF13 and FGFRL1 was investigated in PCa. Previous studies have suggested that ER is protective against PCa whereas FGF8 has been shown to induce PCa in transgenic mice. FGF13 and FGFRL1 are poorly understood members of the FGF and FGFR families, respectively. Transgenic mouse models were used to investigate the ability of inactivated ERβ to facilitate FGF8-induced prostate tumorigenesis. Human PCa tissue microarrays (TMAs) were used to study the expression pattern of FGF13 and FGFRL1 in PCa and the results were correlated to corresponding patient data. The targets and biological functions of FGF13 and FGFRL1 were characterized using experimental in vivo and in vitro models. The results show that deficiency of ERβ, which had been expected to have tumor suppressing capacity, seemed to influence epithelial differentiation but did not affect FGF8-induced prostate tumorigenesis. Analysis of the TMAs showed increased expression of FGF13 in PCa. The level of cytoplasmic FGF13 was associated with the PCa biochemical recurrence (BCR), demonstrated by increasing serum PSA value, and was able to act as an independent prognostic biomarker for PCa patients after RP. Expression of FGFRL1, the most recently identified FGFR, was also elevated in PCa. Cytoplasmic and nuclear FGFRL1 was associated with high Gleason score and Ki67 level whereas the opposite was true for the cell membrane FGFRL1. Silencing of FGFRL1 in PC-3M cells led to a strongly decreased growth rate of these cells as xenografts in nude mice and the experiments with PCa cell lines showed that FGFRL1 is able to modulate the FGF2- and FGF8-induced signaling pathways. The next generation sequencing (NGS) experiments with FGFRL1-silenced PC-3M cells revealed candidates for FGFRL1 target genes. In summary, these studies provide new data on the FGF/FGFR signaling pathways in normal and malignant prostate and suggest a potential role for FGF13 and FGFRL1 as novel prognostic markers for PCa patients. Keywords: FGF8, FGF13, FGFRL1, ERβ, prostate cancer, prognostic marker

High diversity of skin-associated bacterial communities of marine fishes is promoted by their high variability among body parts, individuals and species

Animal-associated microbiotas form complex communities, which are suspected to play crucial functions for their host fitness. However, the biodiversity of these communities, including their differences between host species and individuals, has been scarcely studied, especially in case of skin-associated communities. In addition, the intraindividual variability (i.e. between body parts) has never been assessed to date. The objective of this study was to characterize skin bacterial communities of two teleostean fish species, namely the European seabass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata), using a high-throughput DNA sequencing method. In order to focus on intrinsic factors of host-associated bacterial community variability, individuals of the two species were raised in controlled conditions. Bacterial diversity was assessed using a set of four complementary indices, describing the taxonomic and phylogenetic facets of biodiversity and their respective composition (based on presence/absence data) and structure (based on species relative abundances) components. Variability of bacterial diversity was quantified at the interspecific, interindividual and intraindividual scales. We demonstrated that fish surfaces host highly diverse bacterial communities, whose composition was very different from that of surrounding bacterioplankton. This high total biodiversity of skin-associated communities was supported by the important variability, between host species, individuals and the different body parts (dorsal, anal, pectoral and caudal fins).