Ion torrent high throughput mitochondrial genome sequencing (HTMGS).

The implementation and popularity of next generation sequencing (NGS) has led to the development of various rapid whole mitochondrial genome sequencing techniques. We summarise an efficient and cost-effective NGS approach for mitochondrial genomic DNA in humans using the Ion Torrent platform, and further discuss our bioinformatics pipeline for streamlined variant calling. Ion 316 chips were utilised with the Ion Torrent semi-conductor platform Personal Genome Machine (PGM) to perform tandem sequencing of mitochondrial genomes from the core pedigree (n = 315) of the Norfolk Island Health Study. Key improvements from commercial methods focus on the initial PCR step, which currently requires extensive optimisation to ensure the accurate and reproducible elongation of each section of the complete mitochondrial genome. Dual-platform barcodes were incorporated into our protocol thereby extending its potential application onto Illumina-based systems. Our bioinformatics pipeline consists of a modified version of GATK best practices tailored for mitochondrial genomic data. When compared with current commercial methods, our method, termed high throughput mitochondrial genome sequencing (HTMGS), allows high multiplexing of samples and the use of alternate library preparation reagents at a lower cost per sample (~1.7 times) when compared to current commercial methodologies. Our HTMGS methodology also provides robust mitochondrial sequencing data (>450X average coverage) that can be applied and modified to suit various study designs. On average, we were able to identify ~30 variants per sample with 572 variants observed across 315 samples. We have developed a high throughput sequencing and analysis method targeting complete mitochondrial genomes; with the potential to be platform agnostic with analysis options that adhere to current best practices.

Analysis of pooled genome sequences from Djallonke and Sahelian sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment.

BACKGROUND: The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually. RESULTS: Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis. CONCLUSIONS: Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.

Genetic characterisation of molecular targets in carcinoma of unknown primary.

BACKGROUND: Carcinoma of unknown primary (CUP) is a metastatic epithelial malignancy in the absence of an identifiable primary tumour. Prognosis for patients with CUP is poor because treatment options are generally limited to broad spectrum chemotherapy. A shift towards personalised cancer management based on mutation profiling offers the possibility of new treatment paradigms. This study has explored whether actionable, oncogenic driver mutations are present in CUP that have potential to better inform treatment decisions. METHODS: Carcinoma of unknown primary cases (n = 21) were selected and DNA was isolated from formalin-fixed paraffin embedded sections prior to amplification and sequencing. Two distinct yet complementary targeted gene panels were used to assess variants in up to 76 known cancer-related genes for the identification of biologically relevant and actionable mutations. RESULTS: Variants were detected in 17/21 cases (81%) of which 11 (52%) were potentially actionable with drugs currently approved for use in known primary cancer types or undergoing clinical trials. The most common variants detected were in TP53 (47%), KRAS (12%), MET (12%) and MYC (12%). Differences at the molecular level were seen between common CUP histological subtypes. CUP adenocarcinomas and poorly differentiated carcinomas harboured the highest frequency of variants in genes involved in signal transduction pathways (e.g. MET, EGFR, HRAS, KRAS, and BRAF). In contrast, squamous cell carcinoma exhibited a higher frequency of variants in cell cycle control and DNA repair genes (e.g. TP53, CDKN2A and MLH1). CONCLUSION: Taken together, mutations in biologically relevant genes were detected in the vast majority of CUP tumours, of which half provided a potentially novel treatment option not generally considered in CUP.

16(th) IHIW : review of HLA typing by NGS.

Human leucocyte antigen (HLA) genes play an important role in the success of organ transplantation and are associated with autoimmune and infectious diseases. Current DNA-based genotyping methods, including Sanger sequence-based typing (SSBT), have identified a high degree of polymorphism. This level of polymorphism makes high-resolution HLA genotyping challenging, resulting in ambiguous typing results due to an inability to resolve phase and/or defining polymorphisms lying outside the region amplified. Next-generation sequencing (NGS) may resolve the issue through the combination of clonal amplification, which provides phase information, and the ability to sequence larger regions of genes, including introns, without the additional effort or cost associated with current methods. The NGS HLA sequencing project of the 16IHIW aimed to discuss the different approaches to (i) template preparation including short- and long-range PCR amplicons, exome capture and whole genome; (ii) sequencing platforms, including GS 454 FLX, Ion Torrent PGM, Illumina MiSeq/HiSeq and Pacific Biosciences SMRT; (iii) data analysis, specifically allele-calling software. The pilot studies presented at the workshop demonstrated that although individual sequencers have very different performance characteristics, all produced sequence data suitable for the resolution of HLA genotyping ambiguities. The developments presented at this workshop clearly highlight the potential benefits of NGS in the HLA laboratory.

A patient with a metastatic gastroenteropancreatic endocrine carcinoma causing hyperinsulinaemic hypoglycaemia and the carcinoid syndrome.

We present the case of a 57-year-old patient who initially presented with a constellation of symptoms including intense pruritis, flushing and diarrhoea. Following several months clinical deterioration, the patient was investigated radiologically, where multiple hepatic tumours were identified. Liver biopsy confirmed the presence of a well-differentiated metastatic gastroenteropancreatic endocrine carcinoma with biochemical evidence of serotonin secretion. Over a period of six months, the clinical course of the patient's disease progressed whereby severe hypoglycaemia became the major manifestation. Subsequent biochemical investigations confirmed the diagnosis of an insulinoma. Extensive radiological investigation revealed a solitary primary pancreatic tumour, indicating the presence of a metastatic pancreatic endocrine tumour (PET) secreting both insulin and serotonin. The patient was treated with a chemotherapy regimen consisting of 12 cycles of 5-fluorouracil/oxaliplatin, responding clinically - improved World Health Organization performance score from 3 to 1, biochemically - significantly reduced plasma chromogranin A and cancer antigen 19-9 concentrations and improved liver function tests, and radiologically - reduced pancreatic and hepatic tumour size. This is the first report of a primary PET secreting insulin and serotonin. Due to the association of serotonin-secreting gastroenteropancreatic endocrine tumours (GEP-ETs) with multiple endocrine neoplasia type-1 (MEN1) and biochemical evidence of an insulinoma, MEN1 should also be considered in such cases. The case provides further evidence for the biological heterogeneity of GEP-ETs and the myriad secretory humoral products and resultant clinical syndromes arising from such tumours.

Matrix metalloproteinase-2 gene variants and abdominal aortic aneurysm.

OBJECTIVE: To investigate associations between two polymorphisms of the matrix metalloproteinase-2 gene (MMP2) and the incidence and progression of abdominal aortic aneurysm (AAA). METHODS: Cases and controls were recruited from a trial of screening for AAAs. The association between two variants of MMP2 (-1360C>T, and +649C>T) in men with AAA (n=678) and in controls (n=659) was examined using multivariate analyses. The association with AAA expansion (n=638) was also assessed. RESULTS: In multivariate analyses with adjustments for multiple testing, no association between either SNP and AAA presence or expansion was detected. CONCLUSION: MMP2 -1360C>T and +649C>T variants are not risk factors for AAA.

The evolution and diversity of TNF block haplotypes in European, Asian and Australian Aboriginal populations.

The region spanning the tumour necrosis factor (TNF) cluster in the human major histocompatibility complex is implicated in susceptibility to immunopathological disease, but ethnic differences and linkage disequilibrium have hampered identification of critical polymorphisms. Here, we investigate Europeans, Asians (Bidayuh, Chinese, Indian, Jehai, Malay, Temuan) and Australian Aborigines to provide a framework for disease-association studies. DNA from 999 unrelated healthy donors was genotyped at 38 loci, primarily in coding and promoter regions over a 60-kb region spanning seven genes near TNF. The PHASE algorithm was used to statistically infer TNF block haplotypes and estimate their frequencies in each population. The TNF block is carried as 31 haplotypes in all populations combined, with <19 in any single population. Only six haplotypes have a unique tag single nucleotide polymorphism (SNP) valid for all populations, but seven haplotypes could be tagged with individual SNPs in selected populations. Four to eight TNF block haplotypes exist across all ethnicities, and hence must pre-date the divergence of these populations from a common ancestor >160,000 years ago. Some haplotypes are unique to isolated populations, but they do not contain unique SNP. Hence, they reflect restricted migration and/or extinction of some families rather than de novo mutation.

TNF block haplotypes associated with conserved MHC haplotypes in European, Asian and Australian Aboriginal donors.

Associations between major histocompatibility complex (MHC) ancestral haplotypes (AHs) and immunopathological diseases are traditionally ascribed to human leukocyte antigen (HLA) class I or class II alleles. However, polymorphisms in TNF and nearby genes in the central MHC can influence risk. We have defined TNF block haplotypes in Asian, European and Australian Aboriginal donors and shown conservation of TNF block haplotypes in geographically distinct populations, consistent with a common evolutionary origin. Here we show that most TNF block haplotypes do not align with a single MHC AH and associations often vary with ethnicity. This suggests more recent recombination events between the TNF block and the HLA alleles.

Polymorphisms of the matrix metalloproteinase 9 gene and abdominal aortic aneurysm.

BACKGROUND: Increased matrix metalloproteinase (MMP) 9 activity has been implicated in the formation of abdominal aortic aneurysm (AAA). The aim was to explore the association between potentially functional variants of the MMP-9 gene and AAA. METHODS: The -1562C > T and -1811A > T variants of the MMP-9 gene were genotyped in 678 men with an AAA (at least 30 mm in diameter) and 659 control subjects (aortic diameter 19-22 mm) recruited from a population-based trial of screening for AAA. Levels of MMP-9 were measured in a random subset of 300 cases and 84 controls. The association between genetic variants (including haplotypes) and AAA was assessed by multivariable logistic regression. RESULTS: There was no association between the MMP-9-1562C > T (odds ratio (OR) 0.70 (95 per cent confidence interval (c.i.) 0.27 to 1.82)) or -1811A > T (OR 0.71 (95 per cent c.i. 0.28 to 1.85)) genotypes, or the most common haplotype (OR 0.81 (95 per cent c.i. 0.62 to 1.05)) and AAA. The serum MMP-9 concentration was higher in cases than controls, and in minor allele carriers in cases and controls, although the differences were not statistically significant. CONCLUSION: In this study, the genetic tendency to higher levels of circulating MMP-9 was not associated with AAA.

Polymorphisms of the interleukin-6 gene promoter and abdominal aortic aneurysm.

BACKGROUND: Elevated levels of circulating interleukin-6 (IL-6) have been reported in patients with abdominal aortic aneurysms (AAAs). Although this implicates inflammation as a cause of AAAs, there is also evidence that the aneurysmal aorta may secrete IL-6 into the circulation as a result of aortic proteolysis. Genetic association studies are one means of trying to clarify the role of specific mediators in the causal pathway. The aim of the present study was to examine the association between variants of the IL-6 gene and AAAs. METHODS: An association study involving 677 men with screen-detected AAAs and 656 age-matched controls was performed. Three variants in the IL-6 promoter region were analysed: IL-6-174G>C (rs1800795), IL-6-572G>C (rs1800796) and IL-6-597G>A (rs1800797). Univariate regression of SNP genotype on AAA as a binary outcome was initially performed under a range of genetic models (additive, dominant and recessive). This was followed by multivariate analyses, testing the same models but including risk factors known to be associated with AAAs. All analyses and haplotype estimation were performed under a generalized linear model framework. RESULTS: IL-6-572G>C polymorphism (frequency 1.5% in cases) was identified as an independent risk factor for AAA with an odds ratio (OR) of 6.00 (95%CI: 1.22, 29.41) when applied to the recessive model. No association was seen in the additive or dominant models. In a multivariate analysis using the most common haplotype (h.111, frequency 48.7%) as a reference, h.211 (frequency 4.4%) was an independent risk factor for AAA (OR 1.56, 95%CI: 1.02, 2.39). CONCLUSION: The IL-6 572G>C polymorphism (and h.211 haplotype) is associated with AAA, however it is too rare to be an important cause of most AAAs. This does not support the concept that the elevated level of IL-6 reported in patients with AAAs is a primary cause of the aneurysmal process.

Does a central MHC gene in linkage disequilibrium with HLA-DRB1*0401 affect susceptibility to type 1 diabetes?

Subtypes of HLA-DR4 are associated with susceptibility or protection against type 1 diabetes (T1DM). We addressed whether this reflects linkage disequilibrium with the true susceptibility locus by studying broader MHC haplotypes marked by alleles of HLA-B, IKBL (adjacent to TNFA) and complement C4. The study used a largely Caucasian cohort from Western Australia. HLA-DRB1*0401 and HLA-DRB1*0405 marked susceptibility to T1DM. In Caucasians, DRB1*0401 occurs predominantly in the 44.1 ancestral haplotype (AH; HLA-A2,B44, DRB1*0401,DQB1*0301) and the 62.1AH (HLA-A2,B15(62),DRB1*0401,DQB1*0302). HLA-B15 marked susceptibility and HLA-B44 marked with resistance to T1DM in patients and controls preselected for HLA-DRB1*0401. A gene between TNFA and HLA-B on the 8.1AH (HLA-A1,B8,;DR3,DQ2) modifies the effects of the class II alleles. Here, alleles characteristic of the 62.1AH (C4B3, IKBL+446*T and HLA-A2,B15) were screened in donors preselected for HLA-DRB1*0401. C4B3 was associated with diabetes, consistent with a diabetes gene telomeric of MHC class II. However, increases in carriage of IKBL+446*T and HLA-A2,B15 were marginal, as too few control subjects were available with the diabetogenic alleles. However, with these tools, selection of HLA-DRB1*0401, DQB1*0302 donors who are positive and negative for C4B3 will allow bidirectional mapping of diabetes genes in the central MHC.

Two major histocompatibility complex haplotypes influence susceptibility to sporadic inclusion body myositis: critical evaluation of an association with HLA-DR3.

Previous studies of sporadic inclusion body myositis (sIBM) have shown a strong association with HLA-DR3 and other components of the 8.1 ancestral haplotype (AH) (HLA-A1, B8, DR3), where the susceptibility locus has been mapped to the central major histocompatibility complex (MHC) region between HLA-DR and C4. Here, the association with HLA-DR3 and other genes in the central MHC and class II region was further investigated in a group of 42 sIBM patients and in an ethnically similar control group (n = 214), using single-nucleotide polymorphisms and microsatellite screening. HLA-DR3 (marking DRB1*0301 in Caucasians) was associated with sIBM (Fisher's test). However, among HLA-DR3-positive patients and controls, carriage of HLA-DR3 without microsatellite and single-nucleotide polymorphism alleles of the 8.1AH (HLA-A1, B8, DRB3*0101, DRB1*0301, DQB1*0201) was marginally less common in patients. Patients showed no increase in carriage of the 18.2AH (HLA-A30, B18, DRB3*0202, DRB1*0301, DQB1*0201) or HLA-DR3 without the central MHC of the 8.1AH, further arguing against HLA-DRB1 as the direct cause of susceptibility. Genes between HLA-DRB1 and HOX12 require further investigation. BTL-II lies in this region and is expressed in muscle. Carriage of allele 2 (exon 6) was more common in patients. BTL-II(E6)*2 is characteristic of the 35.2AH (HLA-A3, B35, DRB1*01) in Caucasians and HLA-DR1, BTL-II(E6)*2, HOX12*2, RAGE*2 was carried by several patients. The 8.1AH and 35.2AH may confer susceptibility to sIBM independently or share a critical allele.

A study of the prevalence of systemic sclerosis in northeast England.

OBJECTIVES: We aimed to obtain an estimate of the prevalence and demographics of systemic sclerosis (SSc) and its subtypes at the turn of the millennium. METHODS: Case finding from multiple sources from a defined geographical area. Diagnosis confirmed by clinical examination. RESULTS: The crude prevalence of SSc in northeast England was 8.8 (95% CI: 6.8-10.8) per 100,000. The prevalence when adjusted for the entire UK is 8.2 (95% CI: 6.2-9.8) per 100,000. The ratio of women to men was 5.2:1. The median age of patients was 57.1 yr. The ratio of limited cutaneous SSc to diffuse cutaneous SSc was 4.7:1. Limited cutaneous SSc is associated with the presence of anticentromere antibodies; diffuse cutaneous SSc is associated with anti-Scl 70 antibodies, but either antibody was found in either form of SSc. CONCLUSIONS: SSc appears to be more common in northeast England than was found in the West Midlands in 1986. This may reflect changes in the diagnostic definition of SSc.

Palliation of systemic sclerosis-associated pulmonary hypertension by atrial septostomy.

The onset of pulmonary hypertension in patients with systemic sclerosis carries a poor prognosis. Atrial septostomy has been used successfully to palliate endstage primary pulmonary hypertension but has not been attempted in other forms of pulmonary vascular disease. We report substantial clinical improvement following atrial septostomy in a patient with systemic sclerosis complicated by severe, isolated pulmonary hypertension. After the procedure, exercise capacity was improved and exertional syncope abolished. We suggest that this procedure should be considered for other patients with this diagnosis.

The central MHC gene, BAT1, may encode a protein that down-regulates cytokine production.

BACKGROUND: BAT1 belongs to the DEAD-box family of RNA-binding proteins and is encoded in the central MHC. To determine whether it affects immune responses and hence diseases influenced by MHC haplotypes, U937, THP1 and Jurkat cells were stably transfected with anti-sense DNA corresponding to exons 2-5 of BAT1 using a retroviral vector. RESULTS: Anti-sense transfectants carried anti-sense DNA and expressed anti-sense mRNA. After mitogenic stimulation, they produced higher levels of TNFalpha, IL-1 and IL-6 than equivalent cells carrying the vector alone, suggesting that BAT1 may down-regulate acute phase cytokine production. Polyclonal antibodies raised against a peptide in exon 8 of BAT1 recognized approximately 50 kDa and approximately 38 kDa proteins in all cell lines tested, including the anti-sense transfectants. Expression was localized to the nucleolus in dividing fibroblasts. However the immunochemistry may be confounded by a recently described gene, DDXL, on chromosome 19, which shares a 89% amino acid identity with BAT1. RT-PCR analyses established that BAT1 and DDXL mRNA are expressed in resting U937, THP1 and Jurkat cells. BAT1 and DDXL are divergent in the exons selected for the anti-sense study. CONCLUSIONS: BAT1 is a negative regulator of inflammation. Future studies should address how its functions relate to those of DDXL.

Can MHC class II genes mediate resistance to type 1 diabetes?

Numerous studies have associated carriage of HLA-DRB1*1501, DQA1*0102 and DQB1*0602 (DR15, DQ6) with dominant resistance to type 1 diabetes and have concluded that one or more of the component HLA class II molecules mediate this effect. Mechanisms for MHC class II-mediated resistance to diabetes have been proposed from studies of transgenic mice, usually using the diabetes-prone non-obese diabetic (NOD) strain. However, these studies have not reached any consensus on a plausible mechanism. In this study we question why the role of central MHC genes in resistance to diabetes has not been addressed, as the central MHC carries markers of susceptibility to diabetes in linkage disequilibrium with several genes with known or putative immunoregulatory functions. To illustrate the type of studies required to address this issue, we selected diabetes patients and control subjects for carriage of HLA-DR15 and the C allele at position +738 in the inhibitor of kappa B-like gene (IKBL). These alleles mark the 7.1 haplotype (HLA-A3, B7, IKBL738*C, DR15, DQ6). HLA-DR15 was the most effective marker of resistance, but an effect may be evident with IKBL738*C in a larger study. Moreover, carriage of the entire haplotype was particularly rare in patients. The best explanation for this is that the critical gene lies between IKBL and HLA-DRB1, and is more closely linked to HLA-DRB1. Candidate genes at the centromeric end of the central MHC are reviewed, highlighting the need for further study.