Validation of Next Generation Sequencing Technologies in Comparison to Current Diagnostic Gold Standards for BRAF, EGFR and KRAS Mutational Analysis

Next Generation Sequencing (NGS) has the potential of becoming an important tool in clinical diagnosis and therapeutic decision-making in oncology owing to its enhanced sensitivity in DNA mutation detection, fast-turnaround of samples in comparison to current gold standard methods and the potential to sequence a large number of cancer-driving genes at the one time. We aim to test the diagnostic accuracy of current NGS technology in the analysis of mutations that represent current standard-of-care, and its reliability to generate concomitant information on other key genes in human oncogenesis. Thirteen clinical samples (8 lung adenocarcinomas, 3 colon carcinomas and 2 malignant melanomas) already genotyped for EGFR, KRAS and BRAF mutations by current standard-of-care methods (Sanger Sequencing and q-PCR), were analysed for detection of mutations in the same three genes using two NGS platforms and an additional 43 genes with one of these platforms. The results were analysed using closed platform-specific proprietary bioinformatics software as well as open third party applications. Our results indicate that the existing format of the NGS technology performed well in detecting the clinically relevant mutations stated above but may not be reliable for a broader unsupervised analysis of the wider genome in its current design. Our study represents a diagnostically lead validation of the major strengths and weaknesses of this technology before consideration for diagnostic use.

Mutational analysis of active-site residues in the Mycobacterium leprae RecA intein, a LAGLIDADG homing endonuclease: Asp(122) and Asp(193) are crucial to the double-stranded DNA cleavage activity whereas Asp(218) is not

Mycobacterium leprae recA harbors an in-frame insertion sequence that encodes an intein homing endonuclease (PI-MleI). Most inteins (intein endonucleases) possess two conserved LAGLIDADG (DOD) motifs at their ctive center. A common feature of LAGLIDADG-type homing endonucleases is that they recognize and cleave the same or very similar DNA sequences. However, PI-MleI is distinctive from other members of the family of LAGLIDADG-type HEases for its modular structure with functionally separable domains for DNA-binding and cleavage, each with distinct sequence preferences. Sequence alignment analyses of PI-MleI revealed three putative LAGLIDADG motifs; however, there is conflicting bioinformatics data in regard to their identity and specific location within the intein polypeptide. To resolve this conflict and to determine the active-site residues essential for DNA target site recognition and double-stranded DNA cleavage, we performed site-directed mutagenesis of presumptive catalytic residues in the LAGLIDADG motifs. Analysis of target DNA recognition and kinetic parameters of the wild-type PI-MleI and its variants disclosed that the two amino acid residues, Asp(122) (in Block C) and Asp(193) (in functional Block E), are crucial to the double-stranded DNA endonuclease activity, whereas Asp(218) (in pseudo-Block E) is not. However, despite the reduced catalytic activity, the PI-MleI variants, like the wild-type PI-MleI, generated a footprint of the same length around the insertion site. The D122T variant showed significantly reduced catalytic activity, and D122A and D193A mutations although failed to affect their DNA-binding affinities, but abolished the double-stranded DNA cleavage activity. On the other hand, D122C variant showed approximately twofold higher double-stranded DNA cleavage activity, compared with the wild-type PI-MleI. These results provide compelling evidence that Asp(122) and Asp(193) in DOD motif I and II, respectively, are bona fide active-site residues essential for DNA cleavage activity. The implications of these results are discussed in this report.

Mutational analysis of genes implicated in LPS and capsular polysaccharide biosynthesis in the opportunistic pathogen Bacteroides fragilis

The obligate anaerobe Bacteroides fragilis is a normal resident of the human gastrointestinal tract. The clinically derived B. fragilis strain NCTC 9343 produces an extensive array of extracellular polysaccharides (EPS), including antigenically distinct large, small and micro- capsules. The genome of NCTC 9343 encodes multiple gene clusters potentially involved in the biosynthesis of EPS, eight of which are implicated in production of the antigenically variable micro-capsule. We have developed a rapid and robust method for generating marked and markerless deletions, together with efficient electroporation using unmodified plasmid DNA to enable complementation of mutations. We show that deletion of a putative wzz homologue prevents production of high-molecular-mass polysaccharides (HMMPS), which form the micro-capsule. This observation suggests that micro-capsule HMMPS constitute the distal component of LPS in B. fragilis. The long chain length of this polysaccharide is strikingly different from classical enteric O-antigen, which consists of short-chain polysaccharides. We also demonstrate that deletion of a putative wbaP homologue prevents expression of the phase-variable large capsule and that expression can be restored by complementation. This suggests that synthesis of the large capsule is mechanistically equivalent to production of Escherichia coli group 1 and 4 capsules.

Mutational Analysis of the Rhodopsin Gene in Sector Retinitis Pigmentosa

Background: To determine the role of rhodopsin (RHO) gene mutations in patients with sector retinitis pigmentosa (RP) from Northern Ireland.

Design: A case series of sector RP in a tertiary ocular genetics clinic.

Participants: Four patients with sector RP were recruited from the Royal Victoria Hospital (Belfast, Northern Ireland) and Altnagelvin Hospital (Londonderry, Northern Ireland) following informed consent.

Methods: The diagnosis of sector RP was based on clinical examination, International Society for Clinical Electrophysiology of Vision (ISCEV) standard electrophysiology, and visual field analysis. DNA was extracted from peripheral blood leucocytes and the coding regions and adjacent flanking intronic sequences of the RHO gene were polymerase chain reaction (PCR) amplified and cycle sequenced.

Main Outcome Measure: Rhodopsin mutational status.

Results: A heterozygous missense mutation in RHO (c.173C > T) resulting in a non-conservative substitution of threonine to methionine (p. Thr58Met) was identified in one patient and was absent from 360 control individuals. This non-conservative substitution (p.Thr58Met) replaces a highly evolutionary conserved polar hydrophilic threonine residue with a non-polar hydrophobic methionine residue at position 58 near the cytoplasmic border of helix A of RHO.

Conclusions: The study identified a RHO gene mutation (p.Thr58Met) not previously reported in RP in a patient with sector RP. These findings outline the phenotypic variability associated with RHO mutations. It has been proposed that the regional effects of RHO mutations are likely to result from interplay between mutant alleles and other genetic, epigenetic and environmental factors.

Mutational analysis of NM23-H2/NDP kinase identifies the structural domains critical to recognition of a c-myc regulatory element.

NM23-H2, a presumed regulator of tumor metastasis in humans, is a hexameric protein with both enzymatic (NDP kinase) and regulatory (transcriptional activation) activity. While the structure and catalytic mechanisms have been well characterized, the mode of DNA binding is not known. We examined this latter function in a site-directed mutational study and identified residues and domains essential for the recognition of a c-myc regulatory sequence. Three amino acids, Arg-34, Asn-69, and Lys-135, were found among 30 possibilities to be critical for DNA binding. Two of these, Asn-69 and Lys-135, are not conserved between NM23 variants differing in DNA-binding potential, suggesting that DNA recognition resides partly in nonconserved amino acids. All three DNA-binding defective mutant proteins are active enzymatically and appear to be stable hexamers, suggesting that they perform at the level of DNA recognition and that separate functional domains exist for enzyme catalysis and DNA binding. In the context of the known crystal structure of NM23-H2, the DNA-binding residues are located within distinct structural motifs in the monomer, which are exposed to the surface near the 2-fold axis of adjacent subunits in the hexamer. These findings are explained by a model in which NM23-H2 binds DNA with a combinatorial surface consisting of the "outer" face of the dimer. Chemical crosslinking data support a dimeric DNA-binding mode by NM23-H2.

Mutational analysis of the conserved region 2 site of adenovirus E1A and its effect on binding to the retinoblastoma gene product: use of the "double-tagging" assay.

We have explored the feasibility of using a "double-tagging" assay for assessing which amino acids of a protein are responsible for its binding to another protein. We have chosen the adenovirus E1A-retinoblastoma gene product (pRB) proteins for a model system, and we focused on the high-affinity conserved region 2 of adenovirus E1A (CR2). We used site-specific mutagenesis to generate a mutant E1A gene with a lysine instead of an aspartic acid at position 121 within the CR2 site. We demonstrated that this mutant exhibited little binding to pRB by the double-tagging assay. We also have shown that this lack of binding is not due to any significant decrease in the level of expression of the beta-galactosidase-E1A fusion protein. We then created a "library" of phage expressing beta-galactosidase-E1A fusion proteins with a variety of different mutations within CR2. This library of E1A mutations was used in a double-tagging screening to identify mutant clones that bound to pRB. Three classes of phage were identified: the vast majority of clones were negative and exhibited no binding to pRB. Approximately 1 in 10,000 bound to pRB but not to E1A ("true positives"). A variable number of clones appeared to bind equally well to both pRB and E1A ("false positives"). The DNA sequence of 10 true positive clones yielded the following consensus sequence: DLTCXEX, where X = any amino acid. The recovery of positive clones with only one of several allowed amino acids at each position suggests that most, if not all, of the conserved residues play an important role in binding to pRB. On the other hand, the DNA sequence of the negative clones appeared random. These results are consistent with those obtained from other sources. These data suggest that a double-tagging assay can be employed for determining which amino acids of a protein are important for specifying its interaction with another protein if the complex forms within bacteria. This assay is rapid and up to 1 x 10(6) mutations can be screened at one time.

Analysis of the CDKN2A, CDKN2B and CDK4 genes in 48 Australian melanoma kindreds

Germline mutations within the cyclin-dependent kinase inhibitor 2A (CDKN2A) gene and one of its targets, the cyclin dependent kinase 4 (CDK4) gene, have been identified in a proportion of melanoma kindreds. In the case of CDK4, only one specific mutation, resulting in the substitution of a cysteine for an arginine at codon 24 (R24C), has been found to be associated with melanoma. We have previously reported the identification of germline CDKN2A mutations in 7/18 Australian melanoma kindreds and the absence of the R24C CDK4 mutation in 21 families lacking evidence of a CDKN2A mutation. The current study represents an expansion of these efforts and includes a total of 48 melanoma families from Australia. All of these families have now been screened for mutations within CDKN2A and CDK4, as well as for mutations within the CDKN2A homolog and 9p21 neighbor, the CDKN2B gene, and the alternative exon 1 (E1beta) of CDKN2A. Families lacking CDKN2A mutations, but positive for a polymorphism(s) within this gene, were further evaluated to determine if their disease was associated with transcriptional silencing of one CDKN2A allele. Overall, CDKN2A mutations were detected in 3/30 (10%) of the new kindreds. Two of these mutations have been observed previously: a 24 bp duplication at the 5' end of the gene and a G to C transversion in exon 2 resulting in an M531 substitution. A novel G to A transition in exon 2, resulting in a D108N substitution was also detected. Combined with our previous findings, we have now detected germline CDKN2A mutations in 10/48 (21%) of our melanoma kindreds. In none of the 'CDKN2A-negative' families was melanoma found to segregate with either an untranscribed CDKN2A allele, an R24C CDK4 mutation, a CDKN2B mutation, or an E1beta mutation. The last three observations suggest that these other cell cycle control genes (or alternative gene products) are either not involved at all, or to any great extent, in melanoma predisposition.

Mutation analysis of the CDKN2A promoter in Australian melanoma families

Approximately 50% of all melanoma families worldwide show linkage to 9p21-22, but only about half of these have been shown to contain germ line CDKN2A mutations. It has been hypothesized that a proportion of these families carry mutations in the noncoding regions of CDKN2A. Several Canadian families have been reported to carry a mutation in the 5' UTR, at position -34 relative to the start site, which gives rise to a novel AUG translation initiation codon that markedly decreases translation from the wild-type AUG (Liu et al., 1999). Haplotype sharing in these Canadian families suggested that this mutation is of British origin. We sequenced 1,327 base pairs (bp) of CDKN2A, making up 1,116 bp of the 5' UTR and promoter, all of exon 1, and 61 bp of intron 1, in at least one melanoma case from 110 Australian families with three or more affected members known not to carry mutations within the p16 coding region. In addition, 431 bp upstream of the start codon was sequenced in an additional 253 affected probands from two-case melanoma families for which the CDKN2A mutation status was unknown. Several known polymorphisms at positions -33, -191, -493, and -735 were detected, in addition to four novel variants at positions 120, -252, -347, and -981 relative to the start codon. One of the probands from a two-case family was found to have the previously reported Q50R mutation. No family member was found to carry the mutation at position -34 or any other disease-associated mutation. For further investigation of noncoding CDKN2A mutations that may affect transcription, allele-specific expression analysis was carried out in 31 of the families with at least three affected members who showed either complete or "indeterminate" 9p haplotype sharing without CDKN2A exonic mutations. Reverse transcription polymerase chain reaction and automated sequencing showed expression of both CDKN2A alleles in all family members tested. The lack of CDKN2A promoter mutations and the absence of transcriptional silencing in the germ line of this cohort of families suggest that mutations in the promoter and 5' UTR play a very limited role in melanoma predisposition.

Locked nucleic acid (LNA) single nucleotide polymorphism (SNP) genotype analysis and validation using real-time PCR

With an increased emphasis on genotyping of single nucleotide polymorphisms (SNPs) in disease association studies, the genotyping platform of choice is constantly evolving. In addition, the development of more specific SNP assays and appropriate genotype validation applications is becoming increasingly critical to elucidate ambiguous genotypes. In this study, we have used SNP specific Locked Nucleic Acid (LNA) hybridization probes on a real-time PCR platform to genotype an association cohort and propose three criteria to address ambiguous genotypes. Based on the kinetic properties of PCR amplification, the three criteria address PCR amplification efficiency, the net fluorescent difference between maximal and minimal fluorescent signals and the beginning of the exponential growth phase of the reaction. Initially observed SNP allelic discrimination curves were confirmed by DNA sequencing (n = 50) and application of our three genotype criteria corroborated both sequencing and observed real-time PCR results. In addition, the tested Caucasian association cohort was in Hardy-Weinberg equilibrium and observed allele frequencies were very similar to two independently tested Caucasian association cohorts for the same tested SNP. We present here a novel approach to effectively determine ambiguous genotypes generated from a real-time PCR platform. Application of our three novel criteria provides an easy to use semi-automated genotype confirmation protocol.

Genome-wide DNA methylation analysis of human brain tissue from schizophrenia patients

Recent studies suggest that genetic and environmental factors do not account for all the schizophrenia risk and epigenetics also plays a role in disease susceptibility. DNA methylation is a heritable epigenetic modification that can regulate gene expression. Genome-Wide DNA methylation analysis was performed on post-mortem human brain tissue from 24 patients with schizophrenia and 24 unaffected controls. DNA methylation was assessed at over 485 000 CpG sites using the Illumina Infinium Human Methylation450 Bead Chip. After adjusting for age and post-mortem interval (PMI), 4 641 probes corresponding to 2 929 unique genes were found to be differentially methylated. Of those genes, 1 291 were located in a CpG island and 817 were in a promoter region. These include NOS1, AKT1, DTNBP1, DNMT1, PPP3CC and SOX10 which have previously been associated with schizophrenia. More than 100 of these genes overlap with a previous DNA methylation study of peripheral blood from schizophrenia patients in which 27 000 CpG sites were analysed. Unsupervised clustering analysis of the top 3 000 most variable probes revealed two distinct groups with significantly more people with schizophrenia in cluster one compared to controls (p = 1.74x10-4). The first cluster was composed of 88% of patients with schizophrenia and only 12% controls while the second cluster was composed of 27% of patients with schizophrenia and 73% controls. These results strongly suggest that differential DNA methylation is important in schizophrenia etiology and add support for the use of DNA methylation profiles as a future prognostic indicator of schizophrenia.

Cloning, expression, characterisation and mutational analysis of l-aspartate oxidase from Pseudomonas putida

L-Amino acid oxidases (LAAOs) are useful catalysts for the deracemisation of racemic amino acid sub-strates when combined with abiotic reductants. The gene nadB encoding the L-aspartate amino acid oxidase from Pseudomonas putida (PpLASPO) has been cloned and expressed in E. coli. The purified PpLASPO enzyme displayed a K M for l-aspartic acid of 2.26 mM and a k cat = 10.6 s −1 , with lower activity also displayed towards L-asparagine, for which pronounced substrate inhibition was also observed. The pH optimum of the enzyme was recorded at pH 7.4. The enzyme was stable for 60 min at up to 40 • C, but rapid losses in activity were observed at 50 • C. A mutational analysis of the enzyme, based on its sequence homology with the LASPO from E. coli of known structure, appeared to confirm roles in substrate binding or catalysis for residues His244, His351, Arg386 and Arg290 and also for Thr259 and Gln242. The high activity of the enzyme, and its promiscuous acceptance of both L-asparagine and L-glutamate as substrates, if with low activity, suggests that PpLASPO may provide a good model enzyme for evolution studies towards AAOs of altered or improved properties in the future.

Mitochondrial DNA sequence analysis from multiple gene fragments reveals genetic heterogeneity of Crassostrea ariakensis in East Asia

The native Asian oyster, Crassostrea ariakensis is one of the most common and important Crassostrea species that occur naturally along the coast of East Asia. Molecular species diagnosis is a prerequisite for population genetic analysis of wild oyster populations because oyster species cannot be discriminated reliably using external morphological characters alone due to character ambiguity. To date there have been few phylogeographic studies of natural edible oyster populations in East Asia, in particular this is true of the common species in Korea C. ariakensis. We therefore assessed the levels and patterns of molecular genetic variation in East Asian wild populations of C. ariakensis from Korea, Japan, and China using DNA sequence analysis of five concatenated mtDNA regions namely; 16S rRNA, cytochrome oxidase I, cytochrome oxidase II, cytochrome oxidase III, and cytochrome b. Two divergent C. ariakensis clades were identified between southern China and remaining sites from the northern region. In addition, hierarchical AMOVA and pairwise UST analyses showed that genetic diversity was discontinuous among wild populations of C. ariakensis in East Asia. Biogeographical and historical sea level changes are discussed as potential factors that may have influenced the genetic heterogeneity of wild C. ariakensis stocks across this region.

Identification of a discrete intermediate in the assembly disassembly of physalis mottle tymovirus through mutational analysis

Assembly intermediates of icosahedral viruses are usually transient and are difficult to identify. In the present investigation, site-specific and deletion mutants of the coat protein gene of physalis mottle tymovirus (PhMV) were used to delineate the role of specific amino acid residues in the assembly of the virus and to identify intermediates in this process. N-terminal 30, 34, 35 and 39 amino acid deletion and single C-terminal (N188) deletion mutant proteins of PhMV were expressed in Escherichia coli. Site-specific mutants H69A, C75A, W96A, D144N, D144N-T151A, K143E and N188A were also constructed and expressed. The mutant protein lacking 30 amino acid residues from the N terminus self-assembled to T = 3 particles in vivo while deletions of 34, 35 and 39 amino acid residues resulted in the mutant proteins that were insoluble. Interestingly, the coat protein (pR PhCP) expressed using pRSET B vector with an additional 41 amino acid residues at the N terminus also assembled into T = 3 particles that were more compact and had a smaller diameter. These results demonstrate that the amino-terminal segment is flexible and either the deletion or addition of amino acid residues at the N terminus does not affect T = 3 capsid assembly, in contrast, the deletion of even a single residue from the C terminus (PhN188 Delta 1) resulted in capsids that were unstable. These capsids disassembled to a discrete intermediate with a sedimentation coefficent of 19.4 S. However, the replacement of C-terminal asparagine 188 by alanine led to the formation of stable capsids. The C75A and D144N mutant proteins also assembled into capsids that were as stable as the pR PhCP, suggesting that C75A and D144 are not crucial for the T = 3 capsid assembly. pR PhW96A and pR PhD144N-T151A mutant proteins failed to form capsids and were present as heterogeneous aggregates. Interestingly, the pR PhK143E mutant protein behaved in a manner similar to the C-terminal deletion protein in forming unstable capsids. The intermediate with an s value of 19.4 S was the major assembly product of pR PhH69A mutant protein and could correspond to a 30mer. It is possible that the assembly or disassembly is arrested at a similar stage in pR PhN188 Delta 1, pR PhH69A and pR PhK143E mutant proteins.