A boundary preserving numerical algorithm for the Wright-Fisher model with mutation

The Wright-Fisher model is an Itô stochastic differential equation that was originally introduced to model genetic drift within finite populations and has recently been used as an approximation to ion channel dynamics within cardiac and neuronal cells. While analytic solutions to this equation remain within the interval [0,1], current numerical methods are unable to preserve such boundaries in the approximation. We present a new numerical method that guarantees approximations to a form of Wright-Fisher model, which includes mutation, remain within [0,1] for all time with probability one. Strong convergence of the method is proved and numerical experiments suggest that this new scheme converges with strong order 1/2. Extending this method to a multidimensional case, numerical tests suggest that the algorithm still converges strongly with order 1/2. Finally, numerical solutions obtained using this new method are compared to those obtained using the Euler-Maruyama method where the Wiener increment is resampled to ensure solutions remain within [0,1].

A three-nucleotide mutation altering the Maize streak virus Rep pRBR-interaction motif reduces symptom severity in maize and partially reverts at high frequency without restoring pRBR-Rep binding

Geminivirus infectivity is thought to depend on interactions between the virus replication-associated proteins Rep or RepA and host retinoblastoma-related proteins (pRBR), which control cell-cycle progression. It was determined that the substitution of two amino acids in the Maize streak virus (MSV) RepA pRBR-interaction motif (LLCNE to LLCLK) abolished detectable RepA-pRBR interaction in yeast without abolishing infectivity in maize. Although the mutant virus was infectious in maize, it induced less severe symptoms than the wild-type virus. Sequence analysis of progeny viral DNA isolated from infected maize enabled detection of a high-frequency single-nucleotide reversion of C(601)A in the 3 nt mutated sequence of the Rep gene. Although it did not restore RepA-pRBR interaction in yeast, sequence-specific PCR showed that, in five out of eight plants, the C(601)A reversion appeared by day 10 post-inoculation. In all plants, the C(601)A revertant eventually completely replaced the original mutant population, indicating a high selection pressure for the single-nucleotide reversion. Apart from potentially revealing an alternative or possibly additional function for the stretch of DNA that encodes the apparently non-essential pRBR-interaction motif of MSV Rep, the consistent emergence and eventual dominance of the C(601)A revertant population might provide a useful tool for investigating aspects of MSV biology, such as replication, mutation and evolution rates, and complex population phenomena, such as competition between quasispecies and population turnover. © 2005 SGM.

A deletion and point mutation study of the human papillomavirus type 16 major capsid gene

Recombinant human papillomavirus (HPV) virus-like particles (VLPs) made from the major capsid protein L1 are promising vaccine candidates for use as vaccines against genital and other HPV infections, and particularly against HPV-16. However, HPV-16 genotype variants have different binding affinities for neutralising mouse Mabs raised against HPV-16 L1 VLPs. This paper analyses, using a panel of well-characterised Mabs, the effects on the antigenicity of various C- and N-terminal deletants of HPV-16 L1 made in insect cells via recombinant baculovirus, of an A → T mutation at residue 266 (A266T), and of a C → G mutation at conserved position 428 (C428G). The effects of these changes on assembly of the variant L1s were studied by electron microscopy. Binding of Mab H16:E70 to A266T was reduced by almost half in comparison to wild type L1. Retention of the C-terminal region 428-483 was critical for the binding of conformation-specific Mabs (H16:V5, H16:E70, H16:U4 and H16:9A) whereas deletion of the nuclear localisation signal (NLS) or the C428G mutation or an N-terminal deletion (residues 2-9) did not affect the antigenicity. The N-terminal deletion resulted in a mixed population of 30 and 55 nm VLPs, which differs from the same construct expressed in Escherichia coli, whereas pentamer aggregates resulted from deletion of the 428-465 region or the C428G mutation. The results have implications both for considering use of single-genotype HPV vaccines, and for design of novel second-generation vaccines. © 2006 Elsevier B.V. All rights reserved.

Increased expression of cyclin-dependent kinase inhibitor 2 (CDKN2A) gene product P16INK4A in ovarian cancer is associated with progression and unfavourable prognosis

Paraffin sections from 190 epithelial ovarian tumours, including 159 malignant and 31 benign epithelial tumours, were analysed immunohistochemically for expression of cyclin-dependent kinase inhibitor 2 (CDKN2A) gene product p16INK4A (p16). Most benign tumours showed no p16 expression in the tumour cells, whereas only 11% of malignant cancers were p16 negative. A high proportion of p16-positive tumour cells was associated with advanced stage and grade, and with poor prognosis in cancer patients. For FIGO stage 1 tumours, a high proportion of p16-positive tumour cells was associated with poorer survival, suggesting that accumulation of p16 is an early event of ovarian tumorigenesis. In contrast to tumour cells, high expression of p16 in the surrounding stromal cells was not associated with the stage and grade, but was associated with longer survival. When all parameters were combined in multivariate analysis, high p16 expression in stromal cells was not an independent predictor for survival, indicating that low p16 expression in stromal cells is associated with other markers of tumour progression. High expression of p16 survival in the stromal cells of tumours from long-term survivors suggests that tumour growth is limited to some extent by factors associated with p16 expression in the matrix.

Unique X-linked familial FSGS with co-segregating heart block disorder is associated with a mutation in the NXF5 gene

Focal segmental glomerulosclerosis (FSGS) is the consequence of a disease process that attacks the kidney's filtering system, causing serious scarring. More than half of FSGS patients develop chronic kidney failure within 10 years, ultimately requiring dialysis or renal transplantation. There are currently several genes known to cause the hereditary forms of FSGS (ACTN4, TRPC6, CD2AP, INF2, MYO1E and NPHS2). This study involves a large, unique, multigenerational Australian pedigree in which FSGS co-segregates with progressive heart block with apparent X-linked recessive inheritance. Through a classical combined approach of linkage and haplotype analysis, we identified a 21.19 cM interval implicated on the X chromosome. We then used a whole exome sequencing approach to identify two mutated genes, NXF5 and ALG13, which are located within this linkage interval. The two mutations NXF5-R113W and ALG13-T141L segregated perfectly with the disease phenotype in the pedigree and were not found in a large healthy control cohort. Analysis using bioinformatics tools predicted the R113W mutation in the NXF5 gene to be deleterious and cellular studies support a role in the stability and localization of the protein suggesting a causative role of this mutation in these co-morbid disorders. Further studies are now required to determine the functional consequence of these novel mutations to development of FSGS and heart block in this pedigree and to determine whether these mutations have implications for more common forms of these diseases in the general population.

Detection of a novel mutation in the CACNA1A gene

Familial hemiplegic migraine (FHM) is a rare autosomal dominant subtype of migraine with aura. It is divided into three subtypes FHM1, FHM2 and FHM3, which are caused by mutations in the CACNA1A, ATP1A2 and SCN1A genes respectively. As part of a regular diagnostic service, we investigated 168 patients with FHM symptoms. Samples were tested for mutations contained within the CACNA1A gene. Some tested samples (4.43%) showed an FHM1 mutation, with five of the mutations found in exon 5, one mutation in exon 16 and one in exon 17. Four polymorphisms were also detected, one of which occurred in a large percentage of samples (14.88%). The exon 16 2094G>A polymorphism, however, has been found to occur in healthy Caucasian control populations up to a frequency of 16% and is not considered to be significantly associated with FHM. A finding of significance, found in a single patient, was the detection of a novel mutation in exon 5 that results in a P225H change. The affected individual was an 8-year-old female. The exact phenotypic effect of this mutation is unknown, and further studies are needed to understand the pathophysiology of this mutation in FHM1. New information will allow for diagnostic procedures to be constantly updated, thus improving accuracy of diagnosis. It is possible that new information will also aid the development of new therapeutic agents for the treatment of FHM.

A dominant-negative mutation in the TRESK potassium channel is linked to familial migraine with aura

Migraine with aura is a common, debilitating, recurrent headache disorder associated with transient and reversible focal neurological symptoms. A role has been suggested for the two-pore domain (K2P) potassium channel, TWIK-related spinal cord potassium channel (TRESK, encoded by KCNK18), in pain pathways and general anaesthesia. We therefore examined whether TRESK is involved in migraine by screening the KCNK18 gene in subjects diagnosed with migraine. Here we report a frameshift mutation, F139WfsX24, which segregates perfectly with typical migraine with aura in a large pedigree. We also identified prominent TRESK expression in migraine-salient areas such as the trigeminal ganglion. Functional characterization of this mutation demonstrates that it causes a complete loss of TRESK function and that the mutant subunit suppresses wild-type channel function through a dominant-negative effect, thus explaining the dominant penetrance of this allele. These results therefore support a role for TRESK in the pathogenesis of typical migraine with aura and further support the role of this channel as a potential therapeutic target.

An androgen receptor mutation in the MDA-MB-453 cell line model of molecular apocrine breast cancer compromises receptor activity

Recent evidence indicates that the estrogen receptor-a-negative, androgen receptor (AR)- positive molecular apocrine subtype of breast cancer is driven by AR signaling. The MDA-MB-453 cell line is the prototypical model of this breast cancer subtype; its proliferation is stimulated by androgens such as 5a-dihydrotestosterone (DHT) but inhibited by the progestin medroxyprogesterone acetate (MPA) via AR-mediated mechanisms. We report here that the AR gene in MDAMB- 453 cells contains a G-T transversion in exon 7, resulting in a receptor variant with a glutamine to histidine substitution at amino acid 865 (Q865H) in the ligand binding domain. Compared with wild-type AR, the Q865H variant exhibited reduced sensitivity to DHT and MPA in transactivation assays in MDA-MB-453 and PC-3 cells but did not respond to non-androgenic ligands or receptor antagonists. Ligand binding, molecular modeling, mammalian two-hybrid and immunoblot assays revealed effects of the Q865H mutation on ligand dissociation, AR intramolecular interactions, and receptor stability. Microarray expression profiling demonstrated that DHT and MPA regulate distinct transcriptional programs in MDA-MB-453 cells. Gene Set Enrichment Analysis revealed that DHT- but not MPA-regulated genes were associated with estrogen-responsive transcriptomes from MCF-7 cells and the Wnt signaling pathway. These findings suggest that the divergent proliferative responses of MDA-MB-453 cells to DHT and MPA result from the different genetic programs elicited by these two ligands through the AR-Q865H variant. This work highlights the necessity to characterize additional models of molecular apocrine breast cancer to determine the precise role of AR signaling in this breast cancer subtype. Endocrine-Related Cancer (2012) 19 599–613

A high-throughput panel for identifying clinically relevant mutation profiles in melanoma

Success with molecular-based targeted drugs in the treatment of cancer has ignited extensive research efforts within the field of personalized therapeutics. However, successful application of such therapies is dependent on the presence or absence of mutations within the patient's tumor that can confer clinical efficacy or drug resistance. Building on these findings, we developed a high-throughput mutation panel for the identification of frequently occurring and clinically relevant mutations in melanoma. An extensive literature search and interrogation of the Catalogue of Somatic Mutations in Cancer database identified more than 1,000 melanoma mutations. Applying a filtering strategy to focus on mutations amenable to the development of targeted drugs, we initially screened 120 known mutations in 271 samples using the Sequenom MassARRAY system. A total of 252 mutations were detected in 17 genes, the highest frequency occurred in BRAF (n = 154, 57%), NRAS (n = 55, 20%), CDK4 (n = 8, 3%), PTK2B (n = 7, 2.5%), and ERBB4 (n = 5, 2%). Based on this initial discovery screen, a total of 46 assays interrogating 39 mutations in 20 genes were designed to develop a melanoma-specific panel. These assays were distributed in multiplexes over 8 wells using strict assay design parameters optimized for sensitive mutation detection. The final melanoma-specific mutation panel is a cost effective, sensitive, high-throughput approach for identifying mutations of clinical relevance to molecular-based therapeutics for the treatment of melanoma. When used in a clinical research setting, the panel may rapidly and accurately identify potentially effective treatment strategies using novel or existing molecularly targeted drugs

Identification of TFG (TRK-fused gene) as a putative metastatic melanoma tumor suppressor gene

High density SNP arrays can be used to identify DNA copy number changes in tumors such as homozygous deletions of tumor suppressor genes and focal amplifications of oncogenes. Illumina Human CNV370 Bead chip arrays were used to assess the genome for unbalanced chromosomal events occurring in 39 cell lines derived from stage III metastatic melanomas. A number of genes previously recognized to have an important role in the development and progression of melanoma were identified including homozygous deletions of CDKN2A (13 of 39 samples), CDKN2B (10 of 39), PTEN (3 of 39), PTPRD (3 of 39), TP53 (1 of 39), and amplifications of CCND1 (2 of 39), MITF (2 of 39), MDM2 (1 of 39), and NRAS (1 of 39). In addition, a number of focal homozygous deletions potentially targeting novel melanoma tumor suppressor genes were identified. Because of their likely functional significance for melanoma progression, FAS, CH25H, BMPR1A, ACTA2, and TFG were investigated in a larger cohort of melanomas through sequencing. Nonsynonymous mutations were identified in BMPR1A (1 of 43), ACTA2 (3 of 43), and TFG (5 of 103). A number of potentially important mutation events occurred in TFG including the identification of a mini mutation ‘‘hotspot’’ at amino acid residue 380 (P380S and P380L) and the presence of multiple mutations in two melanomas. Mutations in TFG may have important clinical relevance for current therapeutic strategies to treat metastatic melanoma.

B-Raf mutation : a key player in molecular biology of cancer

B-Raf is one of the more commonly mutated proto-oncogenes implicated in the development of cancers. In this review, we consider the mechanisms and clinical impacts of B-Raf mutations in cancer and discuss the implications for the patient in melanoma, thyroid cancer and colorectal cancer, where B-Raf mutations are particularly common.

Correlation between BRAF mutation and the clinicopathological parameters in papillary thyroid carcinoma with particular reference to follicular variant

Mutation of the BRAF gene is common in thyroid cancer. Follicular variant of papillary thyroid carcinoma is a variant of papillary thyroid carcinoma that has created continuous diagnostic controversies among pathologists. The aims of this study are to (1) investigate whether follicular variant of papillary thyroid carcinoma has a different pattern of BRAF mutation than conventional papillary thyroid carcinoma in a large cohort of patients with typical features of follicular variant of papillary thyroid carcinoma and (2) to study the relationship of clinicopathological features of papillary thyroid carcinomas with BRAF mutation. Tissue blocks from 76 patients with diagnostic features of papillary thyroid carcinomas (40 with conventional type and 36 with follicular variant) were included in the study. From these, DNA was extracted and BRAF V600E mutations were detected by polymerase chain reaction followed by restriction enzyme digestion and sequencing of exon 15. Analysis of the data indicated that BRAF V600E mutation is significantly more common in conventional papillary thyroid carcinoma (58% versus 31%, P = .022). Furthermore, the mutation was often noted in female patients (P = .017), in high-stage cancers (P = .034), and in tumors with mild lymphocytic thyroiditis (P = .006). We concluded that follicular variant of papillary thyroid carcinoma differs from conventional papillary thyroid carcinoma in the rate of BRAF mutation. The results of this study add further information indicating that mutations in BRAF play a role in thyroid cancer development and progression.

Relationship between EGFR expression, EGFR mutation status, and the efficacy of chemotherapy plus cetuximab in FLEX study patients with advanced non-small-cell lung cancer

INTRODUCTION: The phase III FLEX study (NCT00148798) in advanced non-small-cell lung cancer indicated that the survival benefit associated with the addition of cetuximab to cisplatin and vinorelbine was limited to patients whose tumors expressed high levels of epidermal growth factor receptor (EGFR) (immunohistochemistry score of >/=200; scale 0-300). We assessed whether the treatment effect was also modulated in FLEX study patients by tumor EGFR mutation status. METHODS: A tumor mutation screen of EGFR exons 18 to 21 included 971 of 1125 (86%) FLEX study patients. Treatment outcome in low and high EGFR expression groups was analyzed across efficacy endpoints according to tumor EGFR mutation status. RESULTS: Mutations in EGFR exons 18 to 21 were detected in 133 of 971 tumors (14%), 970 of which were also evaluable for EGFR expression level. The most common mutations were exon 19 deletions and L858R (124 of 133 patients; 93%). In the high EGFR expression group (immunohistochemistry score of >/=200), a survival benefit for the addition of cetuximab to chemotherapy was demonstrated in patients with EGFR wild-type (including T790M mutant) tumors. Although patient numbers were small, those in the high EGFR expression group whose tumors carried EGFR mutations may also have derived a survival benefit from the addition of cetuximab to chemotherapy. Response data suggested a cetuximab benefit in the high EGFR expression group regardless of EGFR mutation status. CONCLUSIONS: The survival benefit associated with the addition of cetuximab to first-line chemotherapy for advanced non-small-cell lung cancer expressing high levels of EGFR is not limited by EGFR mutation status.