Factors influencing success of clinical genome sequencing across a broad spectrum of disorders

To assess factors influencing the success of whole-genome sequencing for mainstream clinical diagnosis, we sequenced 217 individuals from 156 independent cases or families across a broad spectrum of disorders in whom previous screening had identified no pathogenic variants. We quantified the number of candidate variants identified using different strategies for variant calling, filtering, annotation and prioritization. We found that jointly calling variants across samples, filtering against both local and external databases, deploying multiple annotation tools and using familial transmission above biological plausibility contributed to accuracy. Overall, we identified disease-causing variants in 21% of cases, with the proportion increasing to 34% (23/68) for mendelian disorders and 57% (8/14) in family trios. We also discovered 32 potentially clinically actionable variants in 18 genes unrelated to the referral disorder, although only 4 were ultimately considered reportable. Our results demonstrate the value of genome sequencing for routine clinical diagnosis but also highlight many outstanding challenges.

Deep Salience: Visual Salience Modeling via Deep Belief Propagation

Visual salience is an intriguing phenomenon observed in biological neural systems. Numerous attempts have been made to model visual salience mathematically using various feature contrasts, either locally or globally. However, these algorithmic models tend to ignore the problem’s biological solutions, in which visual salience appears to arise during the propagation of visual stimuli along the visual cortex. In this paper, inspired by the conjecture that salience arises from deep propagation along the visual cortex, we present a Deep Salience model where a multi-layer model based on successive Markov random fields (sMRF) is proposed to analyze the input image successively through its deep belief propagation. As a result, the foreground object can be automatically separated from the background in a fully unsupervised way. Experimental evaluation on the benchmark dataset validated that our Deep Salience model can consistently outperform eleven state-of-the-art salience models, yielding the higher rates in the precision-recall tests and attaining the best F-measure and mean-square error in the experiments.

Exome sequencing of prostate cancer supports the hypothesis of independent tumour origins

BACKGROUND: Prostate cancer (PCa) is a clinically and pathologically heterogeneous disease. The rapid development of sequencing technology has the potential to deliver new biomarkers with emphasis on aggressive disease and to revolutionise personalised cancer treatment. However, a prostate harbouring cancer commonly contains multiple separate tumour foci, with the potential to aggravate tumour sampling. The level of intraprostatic tumour heterogeneity remains to be determined.

OBJECTIVE: To determine the level of intraprostatic tumour heterogeneity through genome-wide, high-resolution profiling of multiple tumour samples from the same individual.

DESIGN, SETTINGS, AND PARTICIPANTS: Multiple tumour samples were obtained from four individuals following radical prostatectomy. One individual (SWE-1) contained >70% cancer cells in all tumour samples, whereas the other three (SWE-2 to SWE-4) required the use of laser capture microdissection for tumour cell enrichment. Subsequently, DNA was extracted from all tissue samples, and exome sequencing was performed. All tumour foci of SWE-1 were also profiled using a high-resolution array for the identification of copy number alterations (CNA).

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Shared somatic high-frequency single nucleotide variants (SNV) and CNAs were used to infer the level of intraprostatic tumour heterogeneity.

RESULTS AND LIMITATIONS: No high-frequency mutations, common for the three tumour samples of SWE-1, were identified. Ten randomly chosen positions were validated with Sanger sequencing in all foci, which verified the exome data. The high level of intraprostatic heterogeneity was consistent in all individuals. In total, three out of four individuals harboured tumours without an apparent common somatic denominator. Although we cannot exclude the presence of common structural rearrangements, a high-density array was used for the detection of deletions and amplifications in SWE-1, which agreed with the exome data.

CONCLUSIONS: We present evidence for the presence of somatically independent tumours within the same prostate. This finding will have implications for personalised cancer treatment and biomarker discovery.

Mapping protein-DNA interactions using ChIP-sequencing

Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.

Whole Genome Sequencing Reveals Local Transmission Patterns of Mycobacterium bovis in Sympatric Cattle and Badger Populations

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled 'reservoir' host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers. © 2012 Biek et al.

Deep separation of benzene from cyclohexane by liquid extraction using ionic liquids as the solvent

Separation of benzene and cyclohexane is one of the most important and difficult processes in the petrochemical industry, especially for low benzene concentration. In this work, three ionic liquids (ILs), [Bmim][BF ₄], [Bpy][BF ₄], and [Bmim][SCN], were investigated as the solvent in the extraction of benzene from cyclohexane. The corresponding ternary liquid-liquid equilibria (LLE) were experimentally determined at T = 298.15 K and atmospheric pressure. The LLE data were correlated with the nonrandom two-liquid model, and the parameters were fitted. The separation capabilities of the ILs were evaluated in terms of the benzene distribution coefficient and solvent selectivity. The effect of the IL structure on the separation was explained based on a well-founded physical model, COSMO-RS. Finally, the extraction processes were defined, and the operation parameters were analyzed. It shows that the ILs studied are suitable solvents for the extractive separation of benzene and cyclohexane, and their separation efficiency can be generally ranked as [Bmim][BF ₄] > [Bpy][BF ₄] > [Bmim][SCN]. The extraction process for a feed with 15 mol % benzene was optimized. High product purity (cyclohexane 0.997) and high recovery efficiency (cyclohexane 96.9% and benzene 98.1%) can be reached. © 2012 American Chemical Society.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

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. 

Seeking the sun in deep, dark places: mesopelagic sightings of ocean sunfishes (Molidae)

Evidence is presented from publicly available remotely operated vehicle (ROV) footage that suggests deep-water ranging in ocean sunfishes (family Molidae) is more common than typically thought, including a new maximum depth recorded for the southern sunfish Mola ramsayi.

Diagnostic value of increased sensitivity by massively parallel sequencing

Routine molecular diagnostics modalities are unable to confidently detect low frequency mutations (<5-15%) that may indicate response to targeted therapies. We confirm the presence of a low frequency NRAS mutation in a rectal cancer patient using massively parallel sequencing when previous Sanger sequencing results proved negative and Q-PCR testing inconclusive. There is increasing evidence that these low frequency mutations may confer resistance to anti-EGFR therapy. In view of negative/inconclusive Sanger sequencing and Q-PCR results for NRAS mutations in a KRAS wt rectal case, the diagnostic biopsy and 4 distinct subpopulations of cells in the resection specimen after conventional chemo/radiotherapy were massively parallel sequenced using the Ion Torrent PGM. DNA was derived from FFPE rectal cancer tissue and amplicons produced using the Cancer Hotspot Panel V2 and sequenced using semiconductor technology. NRAS mutations were observed at varying frequencies in the patient biopsy (12.2%) and all four subpopulations of cells in the resection with an average frequency of 7.3% (lowest 2.6%). The results of the NGS also provided the mutational status of 49 other genes that may have prognostic or predictive value, including KRAS and PIK3CA. NGS technology has been postulated in diagnostics because of its capability to generate results in large panels of clinically meaningful genes in a cost-effective manner. This case illustrates another potential advantage of this technology: its use for detecting low frequency mutations that may influence therapeutic decisions in cancer treatment.

Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis

Mycobacterium bovis is the causal agent of bovine tuberculosis, one of the most important diseases currently facing the UK cattle industry. Here, we use high-density whole genome sequencing (WGS) in a defined sub-population of M. bovis in 145 cattle across 66 herd breakdowns to gain insights into local spread and persistence. We show that despite low divergence among isolates, WGS can in principle expose contributions of under-sampled host populations to M. bovis transmission. However, we demonstrate that in our data such a signal is due to molecular type switching, which had been previously undocumented for M. bovis. Isolates from farms with a known history of direct cattle movement between them did not show a statistical signal of higher genetic similarity. Despite an overall signal of genetic isolation by distance, genetic distances also showed no apparent relationship with spatial distance among affected farms over distances <5 km. Using simulations, we find that even over the brief evolutionary timescale covered by our data, Bayesian phylogeographic approaches are feasible. Applying such approaches showed that M. bovis dispersal in this system is heterogeneous but slow overall, averaging 2 km/year. These results confirm that widespread application of WGS to M. bovis will bring novel and important insights into the dynamics of M. bovis spread and persistence, but that the current questions most pertinent to control will be best addressed using approaches that more directly integrate WGS with additional epidemiological data.

Moytirra:Discovery of the first known deep-sea hydrothermal vent field on the slow-spreading Mid-Atlantic Ridge north of the Azores

Geological, biological, morphological, and hydrochemical data are presented for the newly discovered Moytirra vent field at 45^oN. This is the only high temperature hydrothermal vent known between the Azores and Iceland, in the North Atlantic and is located on a slow to ultraslow-spreading mid-ocean ridge uniquely situated on the 300 m high fault scarp of the eastern axial wall, 3.5 km from the axial volcanic ridge crest. Furthermore, the Moytirra vent field is, unusually for tectonically controlled hydrothermal vents systems, basalt hosted and perched midway up on the median valley wall and presumably heated by an off-axis magma chamber. The Moytirra vent field consists of an alignment of four sites of venting, three actively emitting "black smoke," producing a complex of chimneys and beehive diffusers. The largest chimney is 18 m tall and vigorously venting. The vent fauna described here are the only ones documented for the North Atlantic (Azores to Reykjanes Ridge) and significantly expands our knowledge of North Atlantic biodiversity. The surfaces of the vent chimneys are occupied by aggregations of gastropods (Peltospira sp.) and populations of alvinocaridid shrimp (Mirocaris sp. with Rimicaris sp. also present). Other fauna present include bythograeid crabs (Segonzacia sp.) and zoarcid fish (Pachycara sp.), but bathymodiolin mussels and actinostolid anemones were not observed in the vent field. The discovery of the Moytirra vent field therefore expands the known latitudinal distributions of several vent-endemic genera in the north Atlantic, and reveals faunal affinities with vents south of the Azores rather than north of Iceland. © 2013. American Geophysical Union. All Rights Reserved.

VentBase:Developing a consensus among stakeholders in the deep-sea regarding environmental impact assessment for deep-sea mining-A workshop report

Mining seafloor massive sulfides for metals is an emergent industry faced with environmental management challenges. These revolve largely around limits to our current understanding of biological variability in marine systems, a challenge common to all marine environmental management. VentBase was established as a forum where academic, commercial, governmental, and non-governmental stakeholders can develop a consensus regarding the management of exploitative activities in the deep-sea. Participants advocate a precautionary approach with the incorporation of lessons learned from coastal studies. This workshop report from VentBase encourages the standardization of sampling methodologies for deep-sea environmental impact assessment. VentBase stresses the need for the collation of spatial data and importance of datasets amenable to robust statistical analyses. VentBase supports the identification of set-asides to prevent the local extirpation of vent-endemic communities and for the post-extraction recolonization of mine sites. © 2013.

Innate immunity in the deep sea hydrothermal vent mussel Bathymodiolus azoricus

The interaction between microorganisms and host defense mechanisms is a decisive factor for the survival of marine bivalves. They rely on cell-mediated and humoral reactions to overcome the pathogens that naturally occur in the marine environment. In order to understand host defense reactions in animals inhabiting extreme environments we investigated some of the components from the immune system of the deep sea hydrothermal vent mussel Bathymodiolus azoricus. Cellular constituents in the hemolymph and extrapallial fluid were examined and led to the identification of three types of hemocytes revealing the granulocytes as the most abundant type of cell. To further characterize hemocyte types, the presence of cell surface carbohydrate epitopes was demonstrated with fluorescent WGA lectin, which was mostly ascribed to the granulocytes. Cellular reactions were then investigated by means of phagocytosis and by the activation of putative MAPKs using the microbial compounds zymosan, glucan, peptidoglycan and lipopolysaccharide. Two bacterial agents, Bacillus subtilis and Vibrio parahaemolyticus, were also used to stimulate hemocytes. The results showed that granulocytes were the main phagocytic cells in both hemolymph and extrapallial fluid of B. azoricus. Western blotting analyses using commercially available antibodies against ERK, p38 and JNK, suggested that these putative kinases are involved in signal transduction pathways during experimental stimulation of B. azoricus hemocytes. The fluorescent Ca²⁺ indicator Fura-2 AM was also insightful in demonstrating hemocyte stimulation in the presence of laminarin or live V. parahaemolyticus. Finally, the expression of the antibacterial gene mytilin was analyzed in gill tissues by means of RT-PCR and whole-mount in situ hybridization. Mytilin transcripts were localized in hemocytes underlying gill epithelium. Moreover, mytilin was induced by exposure of live animals to V. parahaemolyticus. These findings support the premise of a conserved innate immune system in B. azoricus. Such system is comparable to other Bivalves and involves the participation of cellular and humoral components. © 2008 Elsevier Inc. All rights reserved.