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.

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.

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.

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.

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.

The inflammatory macrophage response to MCMV in mice with a retroviral immunodeficiency syndrome (MAIDS).

We have shown that normal C57BL/6J mice are moderately resistant to infection with murine cytomegalovirus (MCMV) and that this resistance is impaired by prior infection with LP-BM5 MuLV, which causes a disease (MAIDS) similar to early HIV-induced disease. This study investigates macrophage function in MAIDS+ mice challenged with MCMV. MAIDS reduces the influx of cells into the peritoneal cavity seen in normal C57BL/6J mice 6 days after MCMV infection. The infiltrates contained cells that resembled activated macrophages, as they took up colloidal gold, expressed the macrophage marker Mac-1, had high levels of acid phosphatase activity, and were lymphocytostatic when co-cultured with activated T cells. MAIDS+ mice had a higher percentage of cells able to take up colloidal gold and higher acid phosphatase activity per cell. The cells were also more lymphocytostatic and produced higher levels of interleukin-1 and tumor necrosis factor-alpha on days 4 and 6 after MCMV infection. Hence, MAIDS enhances baseline and induced macrophage activity, but depresses infiltration into the site of inflammation.

Sporadic primary pulmonary hypertension is associated with germline mutations of the gene encoding BMPR-II, a receptor member of the TGF-beta family.

BACKGROUND: Primary pulmonary hypertension (PPH), resulting from occlusion of small pulmonary arteries, is a devastating condition. Mutations of the bone morphogenetic protein receptor type II gene (BMPR2), a component of the transforming growth factor beta (TGF-beta) family which plays a key role in cell growth, have recently been identified as causing familial PPH. We have searched for BMPR2 gene mutations in sporadic PPH patients to determine whether the same genetic defect underlies the more common form of the disorder. METHODS: We investigated 50 unrelated patients, with a clinical diagnosis of PPH and no identifiable family history of pulmonary hypertension, by direct sequencing of the entire coding region and intron/exon boundaries of the BMPR2 gene. DNA from available parent pairs (n=5) was used to assess the occurrence of spontaneous (de novo) mutations contributing to sporadic PPH. RESULTS: We found a total of 11 different heterozygous germline mutations of the BMPR2 gene in 13 of the 50 PPH patients studied, including missense (n=3), nonsense (n=3), and frameshift (n=5) mutations each predicted to alter the cell signalling response to specific ligands. Parental analysis showed three occurrences of paternal transmission and two of de novo mutation of the BMPR2 gene in sporadic PPH. CONCLUSION: The sporadic form of PPH is associated with germline mutations of the gene encoding the receptor protein BMPR-II in at least 26% of cases. A molecular classification of PPH, based upon the presence or absence of BMPR2 mutations, has important implications for patient management and screening of relatives.

Localization of central MHC genes influencing type I diabetes.

The contribution of MHC class II haplotypes to susceptibility to type I diabetes has been clearly established, and interest has now focused on the effects of additional genes in the MHC region. We have investigated the central MHC alleles on 8.1 ancestral haplotype (HLA-A1, B8, DR3, DQ2), as it is well conserved in Caucasian populations. The HLA-DR3-DQ2 genotype is a recognized risk factor for type I diabetes. Single nucleotide polymorphisms and microsatellites in the MHC were used to map segments of the 8.1 ancestral haplotype carried by type I diabetic and control subjects expressing either HLA-B8 or DR3, but not both these markers. In this way we controlled for the diabetogenic effect of carriage of DR3. Alleles of the 8.1 ancestral haplotype between TNFA-308/D6STNFa and HLA-B were carried with significantly greater frequency in B8(-), DR3(+) type I diabetic patients compared with B8(-), DR3(+) controls. This interval was marked by a BAT1 gene polymorphism and a MIB microsatellite allele.

Susceptibility to multiple sclerosis mediated by HLA-DRB1 is influenced by a second gene telomeric of the TNF cluster.

Susceptibility to multiple sclerosis (MS) is clearly associated with human leukocyte antigen (HLA)-DRB1*1501, but some studies show associations with HLA-B7 and -B18. These are often co-expressed with DRB1*1501 in the ancestral haplotypes (AH) denoted 7.1 (HLA-A3, B7, tumor necrosis factor [TNF]a11b4, DRB1*1501) and 18.1 (HLA-A25, B18, TNFa10b4, DRB 1*1501). Here we present a systematic study of 218 patients and 274 controls typed at all standard class II and TNF microsatellite loci, and a novel non-synonymous polymorphism in the central major histocompatibility complex gene, inhibitor of kappa B-like protein (IKBL). The C allele at IKBL+738 is only found on the 7.1 haplotype. HLA-DRB1*1501 was associated with disease, as expected. When subjects expressing DRB 1*501 were analyzed separately, TNFa11b4 and IKBL+738C were less common in the patients and, hence, mark an allele that mediates resistance which lies telomeric of IKBL. TNFa10b4 and TNFa1b5 were more common in DRB1*1501 patients than in controls. These alleles have been associated with the 18.1 and 18.2 AH, respectively. Since no component of these haplotypes was an independent risk factor in this study, it appears likely that a gene linked to TNFa10b4 and TNFa1b5 modifies the effect of the susceptibility locus marked by HLA-DRB1*1501. Potential candidate genes telomeric of the TNF cluster are discussed.

Hemangioblastoma and renal clear cell carcinoma distinguished by means of the AgNOR method.

In view of the difficulty encountered in distinguishing between 2 degrees renal cell carcinoma (RCC) and hemangioblastoma (HBl) in the central nervous system, the AgNOR technique has been applied empirically to a series of 16 specimens of HBl, 5 primary RCC, and 6 specimens of secondary RCC in the CNS. To avoid tautology, the nature of these was confirmed by immunostaining for epithelial membrane antigen (EMA) and factor VIII-related antigen (FVII RAg). It was found that mean nuclear AgNOR counts in the stromal and endothelial cells of HBl exceeded significantly the counts in the tumor and endothelial cells of RCC, with no overlap in values. It is suggested that the AgNOR method is a useful adjunct in achieving the differential diagnosis of HBl and RCC in the nervous system.

Structure and polymorphism of two stress-activated protein kinase genes centromeric of the MHC: SAPK2a and SAPK4.

As MHC genes are potent determinants of susceptibility to immunopathological diseases, the mapping of SAPK2a (CSBP) and SAPK4 to chromosome 6p 21.2-21.3 suggested that these genes may mediate the effects of the MHC on disease. Here we describe the genomic structure and localisation of both genes approximately 2.3Mb centromeric of HLA-DP. Examination of the complete coding region and selected intronic regions of SAPK2a and SAPK4 from 22 human EBV-transformed B-cell lines of different MHC haplotypes and racial background revealed complete sequence conservation. There were no notable differences in levels of expression of SAPK2a and SAPK4 mRNA in cell lines of different MHC haplotypes or racial origin. Examination of the SAPK2a and SAPK4 sequences from two chimpanzees revealed 3 nucleotide differences between human and chimpanzee in each gene resulting in only one amino acid change in SAPK4, and 6 nucleotide substitutions plus 2 deletions in 600bp of intronic sequence from SAPK4. This highlights the selective pressure placed on these genes to maintain their protein sequence, but does not favour a role in genetic regulation of disease or provide evidence of linkage disequilibrium with the MHC.

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.

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.

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.

The central MHC gene IKBL carries a structural polymorphism that is associated with HLA-A3,B7,DR15.

Susceptibility to several disorders, including insulin-dependent diabetes mellitus and multiple sclerosis, has been associated with alleles of HLA class II genes and loci in the TNF cluster in the central major histocompatibility complex (MHC) region. As recombination within this region is rare, it is difficult to determine which genes are important. This will be facilitated by the identification of functional polymorphisms. Hence we are sequencing reverse transcription-polymerase chain reaction products derived from central MHC genes in well characterized and conserved ancestral haplotypes. Here we address the IKBL gene, which lies near the TNF cluster at the telomeric end of the central MHC. Although the IKBL cDNA sequence was conserved between most ancestral haplotypes, a synonymous nucleotide substitution, a 3' untranslated region substitution, and a single nonsynonymous substitution were identified. The latter (IKBL+738) was present in multiple examples of the 7.1 haplotype [HLA-A3, B7, DR2 (DR15)] and resulted in a cysteine to arginine substitution in a predicted protein kinase C phosphorylation site. This polymorphism did not occur in 18 other common haplotypes from the 10th International Histocompatibility Workshop and thus appears haplospecific. A role for IKBL+738 in the association between HLA-A3,B7,DR2(DR15) and susceptibility to multiple sclerosis is discussed.