Henry M. Wood

Sequence differentiation in regions identified by a genome scan for local adaptation.

Genome scans using large numbers of randomly selected markers have revealed a small proportion of loci that deviate from neutral expectations and so may mark genomic regions that contribute to local adaptation. Measurements of sequence differentiation and identification of genes in these regions is important but difficult, especially in organisms with limited genetic information available. We have followed up a genome scan in the marine gastropod, Littorina saxatilis, by searching a bacterial artificial chromosome library with differentiated and undifferentiated markers, sequencing four bacterial artificial chromosomes and then analysing sequence variation in population samples for fragments at, and close to the original marker polymorphisms. We show that sequence differentiation follows the patterns expected from the original marker frequencies, that differentiated markers identify independent and highly localized sites and that these sites fall outside coding regions. Two differentiated loci are characterized by insertions of putative transposable elements that appear to have increased in frequency recently and which might influence expression of downstream genes. These results provide strong candidate loci for the study of local adaptation in Littorina. They demonstrate an approach that can be applied to follow up genome scans in other taxa and they show that the genome scan approach can lead rapidly to candidate genes in nonmodel organisms.

Using next-generation sequencing for high resolution multiplex analysis of copy number variation from nanogram quantities of DNA from formalin-fixed paraffin-embedded specimens

The use of next-generation sequencing technologies to produce genomic copy number data has recently been described. Most approaches, however, reply on optimal starting DNA, and are therefore unsuitable for the analysis of formalin-fixed paraffin-embedded (FFPE) samples, which largely precludes the analysis of many tumour series. We have sought to challenge the limits of this technique with regards to quality and quantity of starting material and the depth of sequencing required. We confirm that the technique can be used to interrogate DNA from cell lines, fresh frozen material and FFPE samples to assess copy number variation. We show that as little as 5 ng of DNA is needed to generate a copy number karyogram, and follow this up with data from a series of FFPE biopsies and surgical samples. We have used various levels of sample multiplexing to demonstrate the adjustable resolution of the methodology, depending on the number of samples and available resources. We also demonstrate reproducibility by use of replicate samples and comparison with microarray-based comparative genomic hybridization (aCGH) and digital PCR. This technique can be valuable in both the analysis of routine diagnostic samples and in examining large repositories of fixed archival material.

Next-generation sequencing for simultaneous determination of human papillomavirus load, subtype, and associated genomic copy number changes in tumors.

Human papillomavirus (HPV) infection in cases of squamous cell carcinoma of the oropharynx is a powerful predictive and prognostic biomarker. We describe how the use of next-generation sequencing can provide a novel method for the detection of HPV in DNA isolated from formalin-fixed paraffin-embedded tissues. Using this methodology in a cohort of 44 head and neck tumors, we identified the samples that contained HPV sequences, the viral subtype involved, and a direct readout of viral load. Specificity of HPV detection by sequencing compared to traditional detection methods using either PCR or p16 immunohistochemistry was 100%. Sensitivity was 50% when either compared to PCR [confidence interval (CI) = 29% to 71%] or 75% when compared to p16 (CI = 47% to 91%). In addition, we demonstrate the ability of next-generation sequencing to detect other HPV subtypes that would not have been detected by traditional methods, and we demonstrated the ability to apply this method to any tumor and any virus in a panel of eight human cancer cell lines. This methodology also provides a tumor genomic copy number karyogram, and in the samples analyzed here, a lower level of chromosome instability was detected in HPV-positive tumors compared to HPV-negative tumors, as observed in previous studies. Thus, the use of next-generation sequencing for the detection of HPV provides a multiplicity of data with clinical significance in a single test.

Diagnosis of copy number variation by Illumina next generation sequencing is comparable in performance to oligonucleotide array comparative genomic hybridisation.

Array comparative genomic hybridisation (aCGH) profiling is currently the gold standard for genetic diagnosis of copy number. Next generation sequencing technologies provide an alternative and adaptable method of detecting copy number by comparing the number of sequence reads in non-overlapping windows between patient and control samples. Detection of copy number using the BlueGnome 8×60k oligonucleotide aCGH platform was compared with low resolution next generation sequencing using the Illumina GAIIx on 39 patients with developmental delay and/or learning difficulties who were referred to the Leeds Clinical Cytogenetics Laboratory. Sensitivity and workflow of the two platforms were compared. Customised copy number algorithms assessed sequence counts and detected changes in copy number. Imbalances detected on both platforms were compared. Of the thirty-nine patients analysed, all eleven imbalances detected by array CGH and confirmed by FISH or Q-PCR were also detected by CNV-seq. In addition, CNV-seq reported one purported pathogenic copy number variant that was not detected by array CGH. Non-pathogenic, unconfirmed copy number calls were detected by both platforms; however few were concordant between the two. CNV-seq offers an alternative to array CGH for copy number analysis with resolution and future costs comparable to conventional array CGH platforms and with less stringent sample requirements.

MicroRNAs and head and neck cancer: Reviewing the first decade of research

MicroRNAs are a class of non-coding RNA which regulate gene expression. Their discovery in humans in 2000 has led to an explosion in research in this area in terms of their role as a biomarker, therapeutic target as well as trying to elucidate their function. This review aims to summarise the function of microRNAs as well as to examine how dysregulation at any step in their biogenesis or functional pathway can play a role in the development of cancer. We review which microRNAs are implicated as oncogenic or tumour suppressor in head and neck cancer as well as the data available on the use of microRNAs as diagnostic and prognostic marker. We also discuss routes for future research.

HER2 overexpression and amplification as a potential therapeutic target in colorectal cancer: analysis of 3256 patients enrolled in the QUASAR, FOCUS and PICCOLO colorectal cancer trials

HER2 overexpression/amplification is linked to trastuzumab response in breast/gastric cancers. One suggested anti‐EGFR resistance mechanism in colorectal cancer (CRC) is aberrant MEK–AKT pathway activation through HER2 up‐regulation. We assessed HER2‐amplification/overexpression in stage II–III and IV CRC patients, assessing relationships to KRAS/BRAF and outcome. Pathological material was obtained from 1914 patients in the QUASAR stage II–III trial and 1342 patients in stage IV trials (FOCUS and PICCOLO). Tissue microarrays were created for HER2 immunohistochemistry. HER2‐amplification was assessed using FISH and copy number variation. KRAS/BRAF mutation status was assessed by pyrosequencing. Progression‐free survival (PFS) and overall survival (OS) data were obtained for FOCUS/PICCOLO and recurrence and mortality for QUASAR; 29/1342 (2.2%) stage IV and 25/1914 (1.3%) stage II–III tumours showed HER2 protein overexpression. Of the HER2‐overexpressing cases, 27/28 (96.4%) stage IV tumours and 20/24 (83.3%) stage II–III tumours demonstrated HER2 amplification by FISH; 41/47 (87.2%) also showed copy number gains. HER2‐overexpression was associated with KRAS/BRAF wild‐type (WT) status at all stages: in 5.2% WT versus 1.0% mutated tumours (p < 0.0001) in stage IV and 2.1% versus 0.2% in stage II–III tumours (p = 0.01), respectively. HER2 was not associated with OS or PFS. At stage II–III, there was no significant correlation between HER2 overexpression and 5FU/FA response. A higher proportion of HER2‐overexpressing cases experienced recurrence, but the difference was not significant. HER2‐amplification/overexpression is identifiable by immunohistochemistry, occurring infrequently in stage II–III CRC, rising in stage IV and further in KRAS/BRAF WT tumours. The value of HER2‐targeted therapy in patients with HER2‐amplified CRC must be tested in a clinical trial.

A randomised trial of the effect of omega-3 polyunsaturated fatty acid supplements on the human intestinal microbiota

Objective Omega-3 polyunsaturated fatty acids (PUFAs) have anticolorectal cancer (CRC) activity. The intestinal microbiota has been implicated in colorectal carcinogenesis. Dietary omega-3 PUFAs alter the mouse intestinal microbiome compatible with antineoplastic activity. Therefore, we investigated the effect of omega-3 PUFA supplements on the faecal microbiome in middle-aged, healthy volunteers (n=22). Design A randomised, open-label, cross-over trial of 8 weeks’ treatment with 4 g mixed eicosapentaenoic acid/docosahexaenoic acid in two formulations (soft-gel capsules and Smartfish drinks), separated by a 12-week ‘washout’ period. Faecal samples were collected at five time-points for microbiome analysis by 16S ribosomal RNA PCR and Illumina MiSeq sequencing. Red blood cell (RBC) fatty acid analysis was performed by liquid chromatography tandem mass spectrometry. Results Both omega-3 PUFA formulations induced similar changes in RBC fatty acid content, except that drinks were associated with a larger, and more prolonged, decrease in omega-6 PUFA arachidonic acid than the capsule intervention (p=0.02). There were no significant changes in α or β diversity, or phyla composition, associated with omega-3 PUFA supplementation. However, a reversible increased abundance of several genera, including Bifidobacterium, Roseburia and Lactobacillus was observed with one or both omega-3 PUFA interventions. Microbiome changes did not correlate with RBC omega-3 PUFA incorporation or development of omega-3 PUFA-induced diarrhoea. There were no treatment order effects. Conclusion Omega-3 PUFA supplementation induces a reversible increase in several short-chain fatty acid-producing bacteria, independently of the method of administration. There is no simple relationship between the intestinal microbiome and systemic omega-3 PUFA exposure. Trial registration number ISRCTN18662143.

Clinical and genomic analysis of a randomised phase II study evaluating anastrozole and fulvestrant in postmenopausal patients treated for large operable or locally advanced hormone-receptor-positive breast cancer

Background:The aim of this study was to assess the efficacy of neoadjuvant anastrozole and fulvestrant treatment of large operable or locally advanced hormone-receptor-positive breast cancer not eligible for initial breast-conserving surgery, and to identify genomic changes occurring after treatment.Methods:One hundred and twenty post-menopausal patients were randomised to receive 1 mg anastrozole (61 patients) or 500 mg fulvestrant (59 patients) for 6 months. Genomic DNA copy number profiles were generated for a subgroup of 20 patients before and after treatment.Results:A total of 108 patients were evaluable for efficacy and 118 for toxicity. The objective response rate determined by clinical palpation was 58.9% (95% CI=45.0–71.9) in the anastrozole arm and 53.8% (95% CI=39.5–67.8) in the fulvestrant arm. The breast-conserving surgery rate was 58.9% (95% CI=45.0–71.9) in the anastrozole arm and 50.0% (95% CI=35.8–64.2) in the fulvestrant arm. Pathological responses >50% occurred in 24 patients (42.9%) in the anastrozole arm and 13 (25.0%) in the fulvestrant arm. The Ki-67 score fell after treatment but there was no significant difference between the reduction in the two arms (anastrozole 16.7% (95% CI=13.3–21.0) before, 3.2% (95% CI=1.9–5.5) after, n=43; fulvestrant 17.1% (95%CI=13.1–22.5) before, 3.2% (95% CI=1.8–5.7) after, n=38) or between the reduction in Ki-67 in clinical responders and non-responders. Genomic analysis appeared to show a reduction of clonal diversity following treatment with selection of some clones with simpler copy number profiles.Conclusions:Both anastrozole and fulvestrant were effective and well-tolerated, enabling breast-conserving surgery in over 50% of patients. Clonal changes consistent with clonal selection by the treatment were seen in a subgroup of patients.

IMSA: Integrated Metagenomic Sequence Analysis for Identification of Exogenous Reads in a Host Genomic Background

Metagenomics, the study of microbial genomes within diverse environments, is a rapidly developing field. The identification of microbial sequences within a host organism enables the study of human intestinal, respiratory, and skin microbiota, and has allowed the identification of novel viruses in diseases such as Merkel cell carcinoma. There are few publicly available tools for metagenomic high throughput sequence analysis. We present Integrated Metagenomic Sequence Analysis (IMSA), a flexible, fast, and robust computational analysis pipeline that is available for public use. IMSA takes input sequence from high throughput datasets and uses a user-defined host database to filter out host sequence. IMSA then aligns the filtered reads to a user-defined universal database to characterize exogenous reads within the host background. IMSA assigns a score to each node of the taxonomy based on read frequency, and can output this as a taxonomy report suitable for cluster analysis or as a taxonomy map (TaxMap). IMSA also outputs the specific sequence reads assigned to a taxon of interest for downstream analysis. We demonstrate the use of IMSA to detect pathogens and normal flora within sequence data from a primary human cervical cancer carrying HPV16, a primary human cutaneous squamous cell carcinoma carrying HPV 16, the CaSki cell line carrying HPV16, and the HeLa cell line carrying HPV18.

A computational index derived from whole-genome copy number analysis is a novel tool for prognosis in early stage lung squamous cell carcinoma.

Squamous cell carcinoma of the lung is remarkable for the extent to which the same chromosomal abnormalities are detected in individual tumours. We have used next generation sequencing at low coverage to produce high resolution copy number karyograms of a series of 89 non-small cell lung tumours specifically of the squamous cell subtype. Because this methodology is able to create karyograms from formalin-fixed paraffin-embedded material, we were able to use archival stored samples for which survival data were available and correlate frequently occurring copy number changes with disease outcome. No single region of genomic change showed significant correlation with survival. However, adopting a whole-genome approach, we devised an algorithm that relates to total genomic damage, specifically the relative ratios of copy number states across the genome. This algorithm generated a novel index, which is an independent prognostic indicator in early stage squamous cell carcinoma of the lung.