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.

Mapping Protein-DNA Interactions Using ChIP-exo and Illumina-Based Sequencing

Chromatin immunoprecipitation (ChIP) provides a means of enriching DNA associated with transcription factors, histone modifications, and indeed any other proteins for which suitably characterized antibodies are available. Over the years, sequence detection has progressed from quantitative real-time PCR and Southern blotting to microarrays (ChIP-chip) and now high-throughput sequencing (ChIP-seq). This progression has vastly increased the sequence coverage and data volumes generated. This in turn has enabled informaticians to predict the identity of multi-protein complexes on DNA based on the overrepresentation of sequence motifs in DNA enriched by ChIP with a single antibody against a single protein. In the course of the development of high-throughput sequencing, little has changed in the ChIP methodology until recently. In the last three years, a number of modifications have been made to the ChIP protocol with the goal of enhancing the sensitivity of the method and further reducing the levels of nonspecific background sequences in ChIPped samples. In this chapter, we provide a brief commentary on these methodological changes and describe a detailed ChIP-exo method able to generate narrower peaks and greater peak coverage from ChIPped material.

Inter-Laboratory Evaluation of a Next-Generation Sequencing Panel for Acute Myeloid Leukemia

INTRODUCTION: Acute myeloid leukemia (AML) is a heterogeneous clonal disorder often associated with dismal overall survival. The clinical diversity of AML is reflected in the range of recurrent somatic mutations in several genes, many of which have a prognostic and therapeutic value. Targeted next-generation sequencing (NGS) of these genes has the potential for translation into clinical practice. In order to assess this potential, an inter-laboratory evaluation of a commercially available AML gene panel across three diagnostic centres in the UK and Ireland was performed.

METHODS: DNA from six AML patient samples was distributed to each centre and processed using a standardised workflow, including a common sequencing platform, sequencing chips and bioinformatics pipeline. A duplicate sample in each centre was run to assess inter- and intra-laboratory performance.

RESULTS: An average sample read depth of 2725X (range 629-5600) was achieved using six samples per chip, with some variability observed in the depth of coverage generated for individual samples and between centres. A total of 16 somatic mutations were detected in the six AML samples, with a mean of 2.7 mutations per sample (range 1-4) representing nine genes on the panel. 15/16 mutations were identified by all three centres. Allelic frequencies of the mutations ranged from 5.6 to 53.3 % (median 44.4 %), with a high level of concordance of these frequencies between centres, for mutations detected.

CONCLUSION: In this inter-laboratory comparison, a high concordance, reproducibility and robustness was demonstrated using a commercially available NGS AML gene panel and platform.

Understanding next generation sequencing in oncology: A guide for oncologists

DNA sequencing is now faster and cheaper than ever before, due to the development of next generation sequencing (NGS) technologies. NGS is now widely used in the research setting and is becoming increasingly utilised in clinical practice. However, due to evolving clinical commitments, increased workload and lack of training opportunities, many oncologists may be unfamiliar with the terminology and technology involved. This can lead to oncologists feeling daunted by issues such as how to interpret the vast amounts of data generated by NGS and the differences between sequencing platforms. This review article explains common concepts and terminology, summarises the process of DNA sequencing (including data analysis) and discusses the main factors to consider when deciding on a sequencing method. This article aims to improve oncologists' understanding of the most commonly used sequencing platforms and the ongoing challenges faced in expanding the use of NGS into routine clinical practice.

Genetics and Prognostication in Splenic Marginal Zone Lymphoma: Revelations from Deep Sequencing

Purpose: Mounting evidence supports the clinical significance of gene mutations and immunogenetic features in common mature B-cell malignancies.

Experimental Design: We undertook a detailed characterization of the genetic background of splenic marginal zone lymphoma (SMZL), using targeted resequencing and explored potential clinical implications in a multinational cohort of 175 patients with SMZL.

Results: We identified recurrent mutations in TP53 (16%), KLF2 (12%), NOTCH2 (10%), TNFAIP3 (7%), MLL2 (11%), MYD88 (7%), and ARID1A (6%), all genes known to be targeted by somatic mutation in SMZL. KLF2 mutations were early, clonal events, enriched in patients with del(7q) and IGHV1-2*04 B-cell receptor immunoglobulins, and were associated with a short median time to first treatment (0.12 vs. 1.11 years; P = 0.01). In multivariate analysis, mutations in NOTCH2 [HR, 2.12; 95% confidence interval (CI), 1.02–4.4; P = 0.044] and 100% germline IGHV gene identity (HR, 2.19; 95% CI, 1.05–4.55; P = 0.036) were independent markers of short time to first treatment, whereas TP53 mutations were an independent marker of short overall survival (HR, 2.36; 95 % CI, 1.08–5.2; P = 0.03).

Conclusions: We identify key associations between gene mutations and clinical outcome, demonstrating for the first time that NOTCH2 and TP53 gene mutations are independent markers of reduced treatment-free and overall survival, respectively.

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.

Meningeal alveolar soft part sarcoma confirmed by characteristic ASPCR1-TFE3 fusion

Sarcoma metastatic to the brain is uncommon and rarely occurs as the initial manifestation of tumor. Alveolar soft part sarcoma (ASPS) is a rare but well-studied subtype of sarcoma. A 39-year-old man presented with seizures due to a left temporal meningeal-enhancing lesion with striking brain edema on MRI. The patient underwent neurosurgical resection for suspected meningioma. Histology showed large tumor cells clustering and forming small nests, in places with pseudoalveolar pattern. Diastase-resistant periodic acid-Schiff revealed very rare granular and rod-like cytoplasmic inclusions. Immunohistochemistry showed convincing positivity only with vimentin and smooth muscle actin. The histological features were strongly suggestive of ASPS. At the molecular level RT-PCR and sequencing analysis demonstrated ASPCR1-TFE3 fusion confirming the histological diagnosis of ASPS. There was no evidence of primary extracranial tumor by physical examination and on chest and abdominal CT scan 11 months after presentation. ASPS typically arise from the soft tissues of the extremities and develop multiple metastatic deposits usually with a long clinical course. This case may represent primary meningeal ASPS although metastatic deposit from an undiscovered primary site cannot be entirely excluded.

Next-Generation Sequencing of Hereditary Hemochromatosis-Related Genes: Novel Likely Pathogenic Variants Found in the Portuguese Population

Hereditary hemochromatosis (HH) is an autosomal recessive disorder characterized by excessive iron absorption resulting in pathologically increased body iron stores. It is typically associated with common HFE gene mutation (p.Cys282Tyr and p.His63Asp). However, in Southern European populations up to one third of HH patients do not carry the risk genotypes. This study aimed to explore the use of next-generation sequencing (NGS) technology to analyse a panel of iron metabolism-related genes (HFE, TFR2, HJV, HAMP, SLC40A1, and FTL) in 87 non-classic HH Portuguese patients. A total of 1241 genetic alterations were detected corresponding to 53 different variants, 13 of which were not described in the available public databases. Among them, five were predicted to be potentially pathogenic: three novel mutations in TFR2 [two missense (p.Leu750Pro and p.Ala777Val) and one intronic splicing mutation (c.967-1G>C)], one missense mutation in HFE (p.Tyr230Cys), and one mutation in the 5'-UTR of HAMP gene (c.-25G>A). The results reported here illustrate the usefulness of NGS for targeted iron metabolism-related gene panels, as a likely cost-effective approach for molecular genetics diagnosis of non-classic HH patients. Simultaneously, it has contributed to the knowledge of the pathophysiology of those rare iron metabolism-related disorders.

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).

Genetic diversity analysis of isolates belonging to the Photobacterium phosphoreum species group collected from salmon products using AFLP fingerprinting

An accurate amplified fragment length polymorphism (AFLP) method, including three primer sets for the selective amplification step, was developed to display the phylogenetic position of Photobacterium isolates collected from salmon products. This method was efficient for discriminating the three species Photobacterium phosphoreum, Photobacterium iliopiscarium and Photobacterium kishitanii, until now indistinctly gathered in the Photobacterium phosphoreum species group known to be strongly responsible for seafood spoilage. The AFLP fingerprints enabled the isolates to be separated into two main clusters that, according to the type strains, were assigned to the two species P. phosphoreum and P. iliopiscarium. P. kishitanii was not found in the collection. The accuracy of the method was validated by using gyrB-gene sequencing and luxA-gene PCR amplification, which confirmed the species delineation. Most of the isolates of each species were clonally distinct and even those that were isolated from the same source showed some diversity. Moreover, this AFLP method may be an excellent tool for genotyping isolates in bacterial communities and for clarifying our knowledge of the role of the different members of the Photobacterium species group in seafood spoilage.

Production of bioplastics from hydrolysates of paper industry by Haloferax mediterranei

The biorefinery concept has attracted much attention over the last decade due to increasing concerns about the use of fossil resources. In this context emerged the use of bioplastics, namely polyhydroxyalkanoates (PHA). PHA are biocompatible and biodegradable plastics that can be obtained from renewable raw materials and can constitute an alternative solution to conventional plastics. In this work, hydrolysed cellulose pulp, coming from Eucalyptus globulus wood cooking, was used as substrate to the PHA-storing bacteria Haloferax mediterranei. The hydrolysed pulp is rich in simple sugars, mainly glucose (81.96 g.L-1) and xylose (20.90 g.L-1). Tests were made in defined medium with glucose and xylose and in hydrolysate supplemented with salts and yeast extract. Different concentrations of glucose were tested, namely 10, 15, 20, 30 and 40 g.L-1. The best accumulation results (27.1 % of PHA) were obtained in hydrolysate medium with 10 g.L-1. Using this concentration, assays were performed in fed-batch and sequencing batch reactor conditions in order to determine the best feeding strategy. The strategy that led to the best results was fed-batch assay with 24.7 % of PHA. An assay without sterile conditions was performed, in which was obtained the same growth than in sterilization test. Finally it was performed an assay in a bioreactor and a fast growth (0.14 h-1) with high glucose and xylose consumption rates (0.368 g.L-1.h-1 and 0.0947 g.L-1.h-1, respectively) were obtained. However 1.50 g.L-1 of PHA, corresponding to 16.1 % (92.52 % of 3HB and 3HV of 7.48 %) of % PHA were observed. The polymer was further characterized by DSC with a glass transition temperature of -6.07 °C, a melting temperature of 156.3 °C and a melting enthalpy of 63.07 J.g-1, values that are in accordance with the literature. This work recognizes for the first time the suitability of the pulp paper hydrolysate as a substrate for PHA production by H. mediterranei.

Macrophages and CD4 T-cells in rheumatoid arthritis and their modulation by JAK inhibitors

Background: Rheumatoid arthritis (RA) is a chronic inflammatory arthritis that causes significant morbidity and mortality and has no cure. Although early treatment strategies and biologic therapies such as TNFα blocking antibodies have revolutionised treatment, there still remains considerable unmet need. JAK kinase inhibitors, which target multiple inflammatory cytokines, have shown efficacy in treating RA although their exact mechanism of action remains to be determined. Stratified medicine promises to deliver the right drug to the right patient at the right time by using predictive ‘omic biomarkers discovered using bioinformatic and “Big Data” techniques. Therefore, knowledge across the realms of clinical rheumatology, applied immunology, bioinformatics and data science is required to realise this goal. Aim: To use bioinformatic tools to analyse the transcriptome of CD14 macrophages derived from patients with inflammatory arthritis and define a JAK/STAT signature. Thereafter to investigate the role of JAK inhibition on inflammatory cytokine production in a macrophage cell contact activation assay. Finally, to investigate JAK inhibition, following RA synovial fluid stimulation of monocytes. Methods and Results: Using bioinformatic software such as limma from the Bioconductor repository, I determined that there was a JAK/STAT signature in synovial CD14 macrophages from patients with RA and this differed from psoriatic arthritis samples. JAK inhibition using a JAK1/3 inhibitor tofacitinib reduced TNFα production when macrophages were cell contact activated by cytokine stimulated CD4 T-cells. Other pro-inflammatory cytokines such as IL-6 and chemokines such as IP-10 were also reduced. RA synovial fluid failed to stimulate monocytes to phosphorylate STAT1, 3 or 6 but CD4 T-cells activated STAT3 with this stimulus. RNA sequencing of synovial fluid stimulated CD4 T-cells showed an upregulation of SOCS3, BCL6 and SBNO2, a gene associated with RA but with unknown function and tofacitinib reversed this. Conclusion: These studies demonstrate that tofacitinib is effective at reducing inflammatory mediator production in a macrophage cell contact assay and also affects soluble factor mediated stimulation of CD4 T-cells. This suggests that the effectiveness of JAK inhibition is due to inhibition of multiple cytokine pathways such as IL-6, IL-15 and interferon. RNA sequencing is a useful tool to identify non-coding RNA transcripts that are associated with synovial fluid stimulation and JAK inhibition but these require further validation. SBNO2, a gene that is associated with RA, may be biomarker of tofacitinib treatment but requires further investigation and validation in wider disease cohorts.

Transcriptional dynamics of the developing sweet cherry (Prunus avium L.) fruit: sequencing, annotation and expression profiling of exocarp-associated genes

The exocarp, or skin, of fleshy fruit is a specialized tissue that protects the fruit, attracts seed dispersing fruit eaters, and has large economical relevance for fruit quality. Development of the exocarp involves regulated activities of many genes. This research analyzed global gene expression in the exocarp of developing sweet cherry (Prunus avium L., 'Regina'), a fruit crop species with little public genomic resources. A catalog of transcript models (contigs) representing expressed genes was constructed from de novo assembled short complementary DNA (cDNA) sequences generated from developing fruit between flowering and maturity at 14 time points. Expression levels in each sample were estimated for 34 695 contigs from numbers of reads mapping to each contig. Contigs were annotated functionally based on BLAST, gene ontology and InterProScan analyses. Coregulated genes were detected using partitional clustering of expression patterns. The results are discussed with emphasis on genes putatively involved in cuticle deposition, cell wall metabolism and sugar transport. The high temporal resolution of the expression patterns presented here reveals finely tuned developmental specialization of individual members of gene families. Moreover, the de novo assembled sweet cherry fruit transcriptome with 7760 full-length protein coding sequences and over 20 000 other, annotated cDNA sequences together with their developmental expression patterns is expected to accelerate molecular research on this important tree fruit crop.

Next-generation sequencing of hereditary hemochromatosis-related genes: novel likely pathogenic variants found in the Portuguese population

Hereditary hemochromatosis (HH) is an autosomal recessive disorder characterized by excessive iron absorption resulting in pathologically increased body iron stores. It is typically associated with common HFE gene mutation (p.Cys282Tyr and p.His63Asp). However, in Southern European populations up to one third of HH patients do not carry the risk genotypes. This study aimed to explore the use of next-generation sequencing (NGS) technology to analyse a panel of iron metabolism-related genes (HFE, TFR2, HJV, HAMP, SLC40A1, and FTL) in 87 non-classic HH Portuguese patients. A total of 1241 genetic alterations were detected corresponding to 53 different variants, 13 of which were not described in the available public databases. Among them, five were predicted to be potentially pathogenic: three novel mutations in TFR2 [two missense (p.Leu750Pro and p.Ala777Val) and one intronic splicing mutation (c.967-1GNC)], one missense mutation in HFE (p.Tyr230Cys), and one mutation in the 5′-UTR of HAMP gene(c.-25GNA). The results reported here illustrate the usefulness of NGS for targeted iron metabolism-related gene panels, as a likely cost-effective approach for molecular genetics diagnosis of non-classic HH patients. Simultaneously, it has contributed to the knowledge of the pathophysiology of those rare iron metabolism-related disorders.

Cloning of the ω-secalin gene family in a wheat 1BL/1RS translocation line using BAC clone sequencing

Background: Wheat 1BL/1RS translocation lines are planted around the world for their disease resistance and high yield. Most of them are poor in bread making, which is partially caused by ω-secalins that are encoded by the ω-secalin gene family, which is located on the short arm of rye chromosome 1R (1RS). However, information on the structure and evolution of the ω-secalin gene family is still limited. Results: We first generated a physicalmap of the ω-secalin gene family covering 195 kb of the Sec-1 locus based on sequencing three bacterial artificial chromosome (BAC) clones of the 1BL/1RS translocation wheat cultivar Shimai 15. A BAC contig was constructed spanning 168 kb of the Sec-1 locus on 1RS. Twelve ω-secalin genes were arranged in a head-to-tail fashion, separated by 8.2–21.6 kb spacers on the contig, whereas six other ω-secalin genes were arranged head-to-tail, separated by 8.2–8.4 kb of spacers on clone BAC125. The 18 ω-secalin genes can be classified into six types among which eight ω-secalin genes were expressed during seed development. The ω-secalin genes with the 1074-bp open reading frame (ORF) represented the main population. Except for two pseudogenes, the N-terminal of the ω-secalin gene was conserved, whereas variations in the C-terminal led to a change in ORF length. The spacers can be sorted into two classes. Class-1 spacers contained conserved and non-conservative sequences. Conclusion: The ω-secalin gene family consisted of at least 18 members in the 1BL/1RS translocation line cv. Shimai 15. Eight ω-secalin genes were expressed during seed development. Eighteen members may originate from a progenitor with a 1,074-bp ORF. The spacers differed in length and sequence conservation.