KRAS G13D sensitivity to neurofibromin-mediated GTP hydrolysis.

KRAS mutations occur in ∼35% of colorectal cancers and promote tumor growth by constitutively activating the mitogen-activated protein kinase (MAPK) pathway. KRAS mutations at codons 12, 13, or 61 are thought to prevent GAP protein-stimulated GTP hydrolysis and render KRAS-mutated colorectal cancers unresponsive to epidermal growth factor receptor (EGFR) inhibitors. We report here that KRAS G13-mutated cancer cells are frequently comutated with NF1 GAP but NF1 is rarely mutated in cancers with KRAS codon 12 or 61 mutations. Neurofibromin protein (encoded by the NF1 gene) hydrolyzes GTP directly in complex with KRAS G13D, and KRAS G13D-mutated cells can respond to EGFR inhibitors in a neurofibromin-dependent manner. Structures of the wild type and G13D mutant of KRAS in complex with neurofibromin (RasGAP domain) provide the structural basis for neurofibromin-mediated GTP hydrolysis. These results reveal that KRAS G13D is responsive to neurofibromin-stimulated hydrolysis and suggest that a subset of KRAS G13-mutated colorectal cancers that are neurofibromin-competent may respond to EGFR therapies.

Performance comparison of SNP detection tools with illumina exome sequencing data--an assessment using both family pedigree information and sample-matched SNP array data.

To apply exome-seq-derived variants in the clinical setting, there is an urgent need to identify the best variant caller(s) from a large collection of available options. We have used an Illumina exome-seq dataset as a benchmark, with two validation scenarios--family pedigree information and SNP array data for the same samples, permitting global high-throughput cross-validation, to evaluate the quality of SNP calls derived from several popular variant discovery tools from both the open-source and commercial communities using a set of designated quality metrics. To the best of our knowledge, this is the first large-scale performance comparison of exome-seq variant discovery tools using high-throughput validation with both Mendelian inheritance checking and SNP array data, which allows us to gain insights into the accuracy of SNP calling through such high-throughput validation in an unprecedented way, whereas the previously reported comparison studies have only assessed concordance of these tools without directly assessing the quality of the derived SNPs. More importantly, the main purpose of our study was to establish a reusable procedure that applies high-throughput validation to compare the quality of SNP discovery tools with a focus on exome-seq, which can be used to compare any forthcoming tool(s) of interest.

Guanine holes are prominent targets for mutation in cancer and inherited disease.

Single base substitutions constitute the most frequent type of human gene mutation and are a leading cause of cancer and inherited disease. These alterations occur non-randomly in DNA, being strongly influenced by the local nucleotide sequence context. However, the molecular mechanisms underlying such sequence context-dependent mutagenesis are not fully understood. Using bioinformatics, computational and molecular modeling analyses, we have determined the frequencies of mutation at G • C bp in the context of all 64 5'-NGNN-3' motifs that contain the mutation at the second position. Twenty-four datasets were employed, comprising >530,000 somatic single base substitutions from 21 cancer genomes, >77,000 germline single-base substitutions causing or associated with human inherited disease and 16.7 million benign germline single-nucleotide variants. In several cancer types, the number of mutated motifs correlated both with the free energies of base stacking and the energies required for abstracting an electron from the target guanines (ionization potentials). Similar correlations were also evident for the pathological missense and nonsense germline mutations, but only when the target guanines were located on the non-transcribed DNA strand. Likewise, pathogenic splicing mutations predominantly affected positions in which a purine was located on the non-transcribed DNA strand. Novel candidate driver mutations and tissue-specific mutational patterns were also identified in the cancer datasets. We conclude that electron transfer reactions within the DNA molecule contribute to sequence context-dependent mutagenesis, involving both somatic driver and passenger mutations in cancer, as well as germline alterations causing or associated with inherited disease.

Genetic loci associated with delayed clearance of Plasmodium falciparum following artemisinin treatment in Southeast Asia.

The recent emergence of artemisinin-resistant Plasmodium falciparum malaria in western Cambodia could threaten prospects for malaria elimination. Identification of the genetic basis of resistance would provide tools for molecular surveillance, aiding efforts to contain resistance. Clinical trials of artesunate efficacy were conducted in Bangladesh, in northwestern Thailand near the Myanmar border, and at two sites in western Cambodia. Parasites collected from trial participants were genotyped at 8,079 single nucleotide polymorphisms (SNPs) using a P. falciparum-specific SNP array. Parasite genotypes were examined for signatures of recent positive selection and association with parasite clearance phenotypes to identify regions of the genome associated with artemisinin resistance. Four SNPs on chromosomes 10 (one), 13 (two), and 14 (one) were significantly associated with delayed parasite clearance. The two SNPs on chromosome 13 are in a region of the genome that appears to be under strong recent positive selection in Cambodia. The SNPs on chromosomes 10 and 13 lie in or near genes involved in postreplication repair, a DNA damage-tolerance pathway. Replication and validation studies are needed to refine the location of loci responsible for artemisinin resistance and to understand the mechanism behind it; however, two SNPs on chromosomes 10 and 13 may be useful markers of delayed parasite clearance in surveillance for artemisinin resistance in Southeast Asia.

The somatic autosomal mutation matrix in cancer genomes.

DNA damage in somatic cells originates from both environmental and endogenous sources, giving rise to mutations through multiple mechanisms. When these mutations affect the function of critical genes, cancer may ensue. Although identifying genomic subsets of mutated genes may inform therapeutic options, a systematic survey of tumor mutational spectra is required to improve our understanding of the underlying mechanisms of mutagenesis involved in cancer etiology. Recent studies have presented genome-wide sets of somatic mutations as a 96-element vector, a procedure that only captures the immediate neighbors of the mutated nucleotide. Herein, we present a 32 × 12 mutation matrix that captures the nucleotide pattern two nucleotides upstream and downstream of the mutation. A somatic autosomal mutation matrix (SAMM) was constructed from tumor-specific mutations derived from each of 909 individual cancer genomes harboring a total of 10,681,843 single-base substitutions. In addition, mechanistic template mutation matrices (MTMMs) representing oxidative DNA damage, ultraviolet-induced DNA damage, (5m)CpG deamination, and APOBEC-mediated cytosine mutation, are presented. MTMMs were mapped to the individual tumor SAMMs to determine the maximum contribution of each mutational mechanism to the overall mutation pattern. A Manhattan distance across all SAMM elements between any two tumor genomes was used to determine their relative distance. Employing this metric, 89.5% of all tumor genomes were found to have a nearest neighbor from the same tissue of origin. When a distance-dependent 6-nearest neighbor classifier was used, 10.4% of the SAMMs had an Undetermined tissue of origin, and 92.2% of the remaining SAMMs were assigned to the correct tissue of origin. [corrected]. Thus, although tumors from different tissues may have similar mutation patterns, their SAMMs often display signatures that are characteristic of specific tissues.

Mouse endogenous retroviruses can trigger premature transcriptional termination at a distance.

Endogenous retrotransposons have caused extensive genomic variation within mammalian species, but the functional implications of such mobilization are mostly unknown. We mapped thousands of endogenous retrovirus (ERV) germline integrants in highly divergent, previously unsequenced mouse lineages, facilitating a comparison of gene expression in the presence or absence of local insertions. Polymorphic ERVs occur relatively infrequently in gene introns and are particularly depleted from genes involved in embryogenesis or that are highly expressed in embryonic stem cells. Their genomic distribution implies ongoing negative selection due to deleterious effects on gene expression and function. A polymorphic, intronic ERV at Slc15a2 triggers up to 49-fold increases in premature transcriptional termination and up to 39-fold reductions in full-length transcripts in adult mouse tissues, thereby disrupting protein expression and functional activity. Prematurely truncated transcripts also occur at Polr1a, Spon1, and up to ∼5% of other genes when intronic ERV polymorphisms are present. Analysis of expression quantitative trait loci (eQTLs) in recombinant BxD mouse strains demonstrated very strong genetic associations between the polymorphic ERV in cis and disrupted transcript levels. Premature polyadenylation is triggered at genomic distances up to >12.5 kb upstream of the ERV, both in cis and between alleles. The parent of origin of the ERV is associated with variable expression of nonterminated transcripts and differential DNA methylation at its 5'-long terminal repeat. This study defines an unexpectedly strong functional impact of ERVs in disrupting gene transcription at a distance and demonstrates that ongoing retrotransposition can contribute significantly to natural phenotypic diversity.

AVIA: an interactive web-server for annotation, visualization and impact analysis of genomic variations.

MOTIVATION: The plethora of information that emerges from large-scale genome characterization studies has triggered the development of computational frameworks and tools for efficient analysis, interpretation and visualization of genomic data. Functional annotation of genomic variations and the ability to visualize the data in the context of whole genome and/or multiple genomes has remained a challenging task. We have developed an interactive web-based tool, AVIA (Annotation, Visualization and Impact Analysis), to explore and interpret large sets of genomic variations (single nucleotide variations and insertion/deletions) and to help guide and summarize genomic experiments. The annotation, summary plots and tables are packaged and can be downloaded by the user from the email link provided. AVAILABILITY AND IMPLEMENTATION: http://avia.abcc.ncifcrf.gov.

Non-B DB v2.0: a database of predicted non-B DNA-forming motifs and its associated tools.

The non-B DB, available at http://nonb.abcc.ncifcrf.gov, catalogs predicted non-B DNA-forming sequence motifs, including Z-DNA, G-quadruplex, A-phased repeats, inverted repeats, mirror repeats, direct repeats and their corresponding subsets: cruciforms, triplexes and slipped structures, in several genomes. Version 2.0 of the database revises and re-implements the motif discovery algorithms to better align with accepted definitions and thresholds for motifs, expands the non-B DNA-forming motifs coverage by including short tandem repeats and adds key visualization tools to compare motif locations relative to other genomic annotations. Non-B DB v2.0 extends the ability for comparative genomics by including re-annotation of the five organisms reported in non-B DB v1.0, human, chimpanzee, dog, macaque and mouse, and adds seven additional organisms: orangutan, rat, cow, pig, horse, platypus and Arabidopsis thaliana. Additionally, the non-B DB v2.0 provides an overall improved graphical user interface and faster query performance.

Unique microRNA molecular profiles in lung cancer diagnosis and prognosis.

MicroRNA (miRNA) expression profiles for lung cancers were examined to investigate miRNA's involvement in lung carcinogenesis. miRNA microarray analysis identified statistical unique profiles, which could discriminate lung cancers from noncancerous lung tissues as well as molecular signatures that differ in tumor histology. miRNA expression profiles correlated with survival of lung adenocarcinomas, including those classified as disease stage I. High hsa-mir-155 and low hsa-let-7a-2 expression correlated with poor survival by univariate analysis as well as multivariate analysis for hsa-mir-155. The miRNA expression signature on outcome was confirmed by real-time RT-PCR analysis of precursor miRNAs and cross-validated with an independent set of adenocarcinomas. These results indicate that miRNA expression profiles are diagnostic and prognostic markers of lung cancer.

MicroRNA-1 is a candidate tumor suppressor and prognostic marker in human prostate cancer.

We previously reported that miR-1 is among the most consistently down-regulated miRs in primary human prostate tumors. In this follow-up study, we further corroborated this finding in an independent data set and made the novel observation that miR-1 expression is further reduced in distant metastasis and is a candidate predictor of disease recurrence. Moreover, we performed in vitro experiments to explore the tumor suppressor function of miR-1. Cell-based assays showed that miR-1 is epigenetically silenced in human prostate cancer. Overexpression of miR-1 in these cells led to growth inhibition and down-regulation of genes in pathways regulating cell cycle progression, mitosis, DNA replication/repair and actin dynamics. This observation was further corroborated with protein expression analysis and 3'-UTR-based reporter assays, indicating that genes in these pathways are either direct or indirect targets of miR-1. A gene set enrichment analysis revealed that the miR-1-mediated tumor suppressor effects are globally similar to those of histone deacetylase inhibitors. Lastly, we obtained preliminary evidence that miR-1 alters the cellular organization of F-actin and inhibits tumor cell invasion and filipodia formation. In conclusion, our findings indicate that miR-1 acts as a tumor suppressor in prostate cancer by influencing multiple cancer-related processes and by inhibiting cell proliferation and motility.

SRAdb: query and use public next-generation sequencing data from within R.

BACKGROUND: The Sequence Read Archive (SRA) is the largest public repository of sequencing data from the next generation of sequencing platforms including Illumina (Genome Analyzer, HiSeq, MiSeq, .etc), Roche 454 GS System, Applied Biosystems SOLiD System, Helicos Heliscope, PacBio RS, and others. RESULTS: SRAdb is an attempt to make queries of the metadata associated with SRA submission, study, sample, experiment and run more robust and precise, and make access to sequencing data in the SRA easier. We have parsed all the SRA metadata into a SQLite database that is routinely updated and can be easily distributed. The SRAdb R/Bioconductor package then utilizes this SQLite database for querying and accessing metadata. Full text search functionality makes querying metadata very flexible and powerful. Fastq files associated with query results can be downloaded easily for local analysis. The package also includes an interface from R to a popular genome browser, the Integrated Genomics Viewer. CONCLUSIONS: SRAdb Bioconductor package provides a convenient and integrated framework to query and access SRA metadata quickly and powerfully from within R.

An integrated understanding of the physiological response to elevated extracellular phosphate.

Recent studies have suggested that changes in serum phosphate levels influence pathological states associated with aging such as cancer, bone metabolism, and cardiovascular function, even in individuals with normal renal function. The causes are only beginning to be elucidated but are likely a combination of endocrine, paracrine, autocrine, and cell autonomous effects. We have used an integrated quantitative biology approach, combining transcriptomics and proteomics to define a multi-phase, extracellular phosphate-induced, signaling network in pre-osteoblasts as well as primary human and mouse mesenchymal stromal cells. We identified a rapid mitogenic response stimulated by elevated phosphate that results in the induction of immediate early genes including c-fos. The mechanism of activation requires FGF receptor signaling followed by stimulation of N-Ras and activation of AP-1 and serum response elements. A distinct long-term response also requires FGF receptor signaling and results in N-Ras activation and expression of genes and secretion of proteins involved in matrix regulation, calcification, and angiogenesis. The late response is synergistically enhanced by addition of FGF23 peptide. The intermediate phase results in increased oxidative phosphorylation and ATP production and is necessary for the late response providing a functional link between the phases. Collectively, the results define elevated phosphate, as a mitogen and define specific mechanisms by which phosphate stimulates proliferation and matrix regulation. Our approach provides a comprehensive understanding of the cellular response to elevated extracellular phosphate, functionally connecting temporally coordinated signaling, transcriptional, and metabolic events with changes in long-term cell behavior.

Transcription signatures encoded by ultraconserved genomic regions in human prostate cancer.

BACKGROUND: Ultraconserved regions (UCR) are genomic segments of more than 200 base pairs that are evolutionarily conserved among mammalian species. They are thought to have functions as transcriptional enhancers and regulators of alternative splicing. Recently, it was shown that numerous RNAs are transcribed from these regions. These UCR-encoded transcripts (ucRNAs) were found to be expressed in a tissue- and disease-specific manner and may interfere with the function of other RNAs through RNA: RNA interactions. We hypothesized that ucRNAs have unidentified roles in the pathogenesis of human prostate cancer. In a pilot study, we examined ucRNA expression profiles in human prostate tumors. METHODS: Using a custom microarray with 962 probesets representing sense and antisense sequences for the 481 human UCRs, we examined ucRNA expression in resected, fresh-frozen human prostate tissues (57 tumors, 7 non-cancerous prostate tissues) and in cultured prostate cancer cells treated with either epigenetic drugs (the hypomethylating agent, 5-Aza 2'deoxycytidine, and the histone deacetylase inhibitor, trichostatin A) or a synthetic androgen, R1881. Expression of selected ucRNAs was also assessed by qRT-PCR and NanoString®-based assays. Because ucRNAs may function as RNAs that target protein-coding genes through direct and inhibitory RNA: RNA interactions, computational analyses were applied to identify candidate ucRNA:mRNA binding pairs. RESULTS: We observed altered ucRNA expression in prostate cancer (e.g., uc.106+, uc.477+, uc.363 + A, uc.454 + A) and found that these ucRNAs were associated with cancer development, Gleason score, and extraprostatic extension after controlling for false discovery (false discovery rate < 5% for many of the transcripts). We also identified several ucRNAs that were responsive to treatment with either epigenetic drugs or androgen (R1881). For example, experiments with LNCaP human prostate cancer cells showed that uc.287+ is induced by R1881 (P < 0.05) whereas uc.283 + A was up-regulated following treatment with combined 5-Aza 2'deoxycytidine and trichostatin A (P < 0.05). Additional computational analyses predicted RNA loop-loop interactions of 302 different sense and antisense ucRNAs with 1058 different mRNAs, inferring possible functions of ucRNAs via direct interactions with mRNAs. CONCLUSIONS: This first study of ucRNA expression in human prostate cancer indicates an altered transcript expression in the disease.

Analysis of 454 sequencing error rate, error sources, and artifact recombination for detection of Low-frequency drug resistance mutations in HIV-1 DNA.

BACKGROUND: 454 sequencing technology is a promising approach for characterizing HIV-1 populations and for identifying low frequency mutations. The utility of 454 technology for determining allele frequencies and linkage associations in HIV infected individuals has not been extensively investigated. We evaluated the performance of 454 sequencing for characterizing HIV populations with defined allele frequencies. RESULTS: We constructed two HIV-1 RT clones. Clone A was a wild type sequence. Clone B was identical to clone A except it contained 13 introduced drug resistant mutations. The clones were mixed at ratios ranging from 1% to 50% and were amplified by standard PCR conditions and by PCR conditions aimed at reducing PCR-based recombination. The products were sequenced using 454 pyrosequencing. Sequence analysis from standard PCR amplification revealed that 14% of all sequencing reads from a sample with a 50:50 mixture of wild type and mutant DNA were recombinants. The majority of the recombinants were the result of a single crossover event which can happen during PCR when the DNA polymerase terminates synthesis prematurely. The incompletely extended template then competes for primer sites in subsequent rounds of PCR. Although less often, a spectrum of other distinct crossover patterns was also detected. In addition, we observed point mutation errors ranging from 0.01% to 1.0% per base as well as indel (insertion and deletion) errors ranging from 0.02% to nearly 50%. The point errors (single nucleotide substitution errors) were mainly introduced during PCR while indels were the result of pyrosequencing. We then used new PCR conditions designed to reduce PCR-based recombination. Using these new conditions, the frequency of recombination was reduced 27-fold. The new conditions had no effect on point mutation errors. We found that 454 pyrosequencing was capable of identifying minority HIV-1 mutations at frequencies down to 0.1% at some nucleotide positions. CONCLUSION: Standard PCR amplification results in a high frequency of PCR-introduced recombination precluding its use for linkage analysis of HIV populations using 454 pyrosequencing. We designed a new PCR protocol that resulted in a much lower recombination frequency and provided a powerful technique for linkage analysis and haplotype determination in HIV-1 populations. Our analyses of 454 sequencing results also demonstrated that at some specific HIV-1 drug resistant sites, mutations can reliably be detected at frequencies down to 0.1%.

Non-B DB: a database of predicted non-B DNA-forming motifs in mammalian genomes.

Although the capability of DNA to form a variety of non-canonical (non-B) structures has long been recognized, the overall significance of these alternate conformations in biology has only recently become accepted en masse. In order to provide access to genome-wide locations of these classes of predicted structures, we have developed non-B DB, a database integrating annotations and analysis of non-B DNA-forming sequence motifs. The database provides the most complete list of alternative DNA structure predictions available, including Z-DNA motifs, quadruplex-forming motifs, inverted repeats, mirror repeats and direct repeats and their associated subsets of cruciforms, triplex and slipped structures, respectively. The database also contains motifs predicted to form static DNA bends, short tandem repeats and homo(purine•pyrimidine) tracts that have been associated with disease. The database has been built using the latest releases of the human, chimp, dog, macaque and mouse genomes, so that the results can be compared directly with other data sources. In order to make the data interpretable in a genomic context, features such as genes, single-nucleotide polymorphisms and repetitive elements (SINE, LINE, etc.) have also been incorporated. The database is accessed through query pages that produce results with links to the UCSC browser and a GBrowse-based genomic viewer. It is freely accessible at http://nonb.abcc.ncifcrf.gov.

Long-range enhancers on 8q24 regulate c-Myc.

Recent genomewide association studies have found multiple genetic variants on chromosome 8q24 that are significantly associated with an increased susceptibility to prostate, colorectal, and breast cancer. These risk loci are located in a "gene desert," a few hundred kilobases telomeric to the Myc gene. To date, the biological mechanism(s) underlying these associations remain unclear. It has been speculated that these 8q24 genetic variant(s) might affect Myc expression by altering its regulation or amplification status. Here, we show that multiple enhancer elements are present within this region and that they can regulate transcription of Myc. We also demonstrate that one such enhancer element physically interacts with the Myc promoter via transcription factor Tcf-4 binding and acts in an allele specific manner to regulate Myc expression.

MouseIndelDB: a database integrating genomic indel polymorphisms that distinguish mouse strains.

MouseIndelDB is an integrated database resource containing thousands of previously unreported mouse genomic indel (insertion and deletion) polymorphisms ranging from approximately 100 nt to 10 Kb in size. The database currently includes polymorphisms identified from our alignment of 26 million whole-genome shotgun sequence traces from four laboratory mouse strains mapped against the reference C57BL/6J genome using GMAP. They can be queried on a local level by chromosomal coordinates, nearby gene names or other genomic feature identifiers, or in bulk format using categories including mouse strain(s), class of polymorphism(s) and chromosome number. The results of such queries are presented either as a custom track on the UCSC mouse genome browser or in tabular format. We anticipate that the MouseIndelDB database will be widely useful for research in mammalian genetics, genomics, and evolutionary biology. Access to the MouseIndelDB database is freely available at: http://variation.osu.edu/.

bioDBnet: the biological database network.

SUMMARY: bioDBnet is an online web resource that provides interconnected access to many types of biological databases. It has integrated many of the most commonly used biological databases and in its current state has 153 database identifiers (nodes) covering all aspects of biology including genes, proteins, pathways and other biological concepts. bioDBnet offers various ways to work with these databases including conversions, extensive database reports, custom navigation and has various tools to enhance the quality of the results. Importantly, the access to bioDBnet is updated regularly, providing access to the most recent releases of each individual database. AVAILABILITY: http://biodbnet.abcc.ncifcrf.gov.

Gene expression profiling of substantia nigra dopamine neurons: further insights into Parkinson's disease pathology.

Parkinson's disease is caused by a progressive loss of the midbrain dopamine (DA) neurons in the substantia nigra pars compacta. Although the main cause of Parkinson's disease remains unknown, there is increasing evidence that it is a complex disorder caused by a combination of genetic and environmental factors, which affect key signalling pathways in substantia nigra DA neurons. Insights into pathogenesis of Parkinson's disease stem from in vitro and in vivo models and from postmortem analyses. Recent technological developments have added a new dimension to this research by determining gene expression profiles using high throughput microarray assays. However, many of the studies reported to date were based on whole midbrain dissections, which included cells other than DA neurons. Here, we have used laser microdissection to isolate single DA neurons from the substantia nigra pars compacta of controls and subjects with idiopathic Parkinson's disease matched for age and postmortem interval followed by microarrays to analyse gene expression profiling. Our data confirm a dysregulation of several functional groups of genes involved in the Parkinson's disease pathogenesis. In particular, we found prominent down-regulation of members of the PARK gene family and dysregulation of multiple genes associated with programmed cell death and survival. In addition, genes for neurotransmitter and ion channel receptors were also deregulated, supporting the view that alterations in electrical activity might influence DA neuron function. Our data provide a 'molecular fingerprint identity' of late-stage Parkinson's disease DA neurons that will advance our understanding of the molecular pathology of this disease.