Search

An enhanced requirement for extracellular nutrients is a hallmark property of cancer cells. Here, we optimized anin vivogenetic screening strategy for evaluating dependencies in acute myeloid leukemia (AML), which led to the identification of the myo-inositol transporter SLC5A3 as a unique vulnerability in this disease. In accord with this transport function, we demonstrate that the SLC5A3 dependency reflects a myo-inositol auxotrophy in AML. Importantly, the commonality among SLC5A3-dependent AML lines is the transcriptional silencing ofISYNA1, which encodes the rate limiting enzyme for myoinositol biosynthesis, inositol-3-phosphate synthase 1. We used gain- and loss-of-function experiments to demonstrate a synthetic lethal genetic interaction betweenISYNA1andSLC5A3in AML, which function redundantly to sustain intracellular myo-inositol. Transcriptional silencing and DNA hypermethylation ofISYNA1occur in a recurrent manner in human AML patient samples, in association with the presence ofIDH1/IDH2andCEBPAmutations. Collectively, our findings reveal myo-inositol auxotrophy as a novel form of metabolic dysregulation in AML, which is caused by the aberrant silencing of a biosynthetic enzyme.Statement of significanceHere, we show how epigenetic silencing can provoke a nutrient dependency in AML by exploiting a synthetic lethality relationship between biosynthesis and transport of myo-inositol. Blocking the function of this solute carrier may have therapeutic potential in an epigenetically-defined subset of AML.

MaizeCODE is a project aimed at identifying and analyzing functional elements in the maize genome. In its initial phase, MaizeCODE assayed up to five tissues from four maize strains (B73, NC350, W22, TIL11) by RNA-Seq, Chip-Seq, RAMPAGE, and small RNA sequencing. To facilitate reproducible science and provide both human and machine access to the MaizeCODE data, we enhanced SciApps, a cloud-based portal, for analysis and distribution of both raw data and analysis results. Based on the SciApps workflow platform, we generated new components to support the complete cycle of MaizeCODE data management. These include publicly accessible scientific workflows for the reproducible and shareable analysis of various functional data, a RESTful API for batch processing and distribution of data and metadata, a searchable data page that lists each MaizeCODE experiment as a reproducible workflow, and integrated JBrowse genome browser tracks linked with workflows and metadata. The SciApps portal is a flexible platform that allows the integration of new analysis tools, workflows, and genomic data from multiple projects. Through metadata and a ready-to-compute cloud-based platform, the portal experience improves access to the MaizeCODE data and facilitates its analysis.

MotivationExome sequencing is a powerful technique for the identification of disease-causing genes. A number of Mendelian inherited disease genes have been identified through this method. However, it remains a challenge to leverage exome sequencing for the study of complex disorders, such as schizophrenia and bipolar disorder, due to the genetic and phenotypic heterogeneity of these disorders. Although not feasible for many studies, sequencing large sample sizes (>10,000) may improve statistical power to associate more variants, while the aggregation of distinct rare variants associated with a given disease can make the identification of causal genes statistically challenging. Therefore, new methods for rare variant association are imperative to identify causative genes of complex disorders.ResultsHere we propose a method to predict causative rare variants using a popular probabilistic problem: The Birthday Model, which estimates the probability that multiple individuals in a group share the same birthday. We consider the probability and coincidence of samples sharing a variant akin to the chance of individuals sharing the same birthday. We investigated the parameter effects of our model, providing guidelines for its use and interpretation of the results. Using published data on autism spectrum disorder, hypertriglyceridemia in addition to a current case-control study on bipolar disorder, we evaluated this probabilistic method to identify potential causative variants. Several genes in the top results of the case-control study were associated with autism spectrum and bipolar disorder. Given that the core probability based on the birthday model is very sensitive to low recurrence, the method successfully tests the association of rare variants, which generally do not provide enough signal in commonly used statistical tests. Importantly, the simplicity of the model allows quick interpretation of genomic data, enabling users to select gene candidates for further biological validation of specific mutations and downstream functional or other studies.Availabilityhttps://github.com/yberstein/Birthday-Alqorithmhttp://labshare.cshl.edu/shares/mccombielab/www-data/Birthday-Algorithm/Birthday-Alqorithm.htmlContactyberstei@cshl.edu (or yaelberstein@gmail.com)Supplementary informationSupplementary data are available online.

Mice lacking DIX domain containing-1 (DIXDC1), an intracellular Wnt/β-catenin signal pathway protein, have abnormal measures of anxiety, depression and social behavior. Pyramidal neurons in these animals' brains have reduced dendritic spines and glutamatergic synapses. Treatment with lithium or a glycogen synthase kinase-3 (GSK3) inhibitor corrects behavioral and neurodevelopmental phenotypes in these animals. Analysis of DIXDC1 in over 9000 cases of autism, bipolar disorder and schizophrenia reveals higher rates of rare inherited sequence-disrupting single-nucleotide variants (SNVs) in these individuals compared with psychiatrically unaffected controls. Many of these SNVs alter Wnt/β-catenin signaling activity of the neurally predominant DIXDC1 isoform; a subset that hyperactivate this pathway cause dominant neurodevelopmental effects. We propose that rare missense SNVs in DIXDC1 contribute to psychiatric pathogenesis by reducing spine and glutamatergic synapse density downstream of GSK3 in the Wnt/β-catenin pathway.

Chromatin accessibility provides an important window into the regulation of gene expression. Recently, the Assay of Transposase Accessible Chromatin with sequencing (ATAC-seq) was developed to profile genome-wide chromatin accessibility. In this work we analyzed 193 ATAC-seq samples from 12 studies to determine the robustness of peaks, as well as their specificity across studies. We employ a read-downsampling approach and find that even samples which have high average read-depth yield poorly powered peaks using conventional pipelines. Finding that many peaks are sensitive to slight perturbations in read-sampling, we develop an approach to improve the robustness of peak signals and apply this novel method to detect differential accessibility between the amygdala and cortex. Most peak calling results heavily and unreasonably depend on the absolute number of reads, easily yielding artifacts when differential peak evaluations are performed. These problems can be ameliorated by resampling reads to determine the robust aggregate result.

Schizophrenia is a serious psychiatric disorder with a broadly undiscovered genetic etiology. Recent studies of de novo mutations (DNMs) in schizophrenia and autism have reinforced the hypothesis that rare genetic variation contributes to risk. We carried out exome sequencing on 57 trios with sporadic or familial schizophrenia. In sporadic trios, we observed a ~3.5-fold increase in the proportion of nonsense DNMs (0.101 vs 0.031, empirical P=0.01, Benjamini-Hochberg-corrected P=0.044). These mutations were significantly more likely to occur in genes with highly ranked probabilities of haploinsufficiency (P=0.0029, corrected P=0.006). DNMs of potential functional consequence were also found to occur in genes predicted to be less tolerant to rare variation (P=2.01 × 10(-)(5), corrected P=2.1 × 10(-)(3)). Genes with DNMs overlapped with genes implicated in autism (for example, AUTS2, CHD8 and MECP2) and intellectual disability (for example, HUWE1 and TRAPPC9), supporting a shared genetic etiology between these disorders. Functionally CHD8, MECP2 and HUWE1 converge on epigenetic regulation of transcription suggesting that this may be an important risk mechanism. Our results were consistent in an analysis of additional exome-based sequencing studies of other neurodevelopmental disorders. These findings suggest that perturbations in genes, which function in the epigenetic regulation of brain development and cognition, could have a central role in the susceptibility to, pathogenesis and treatment of mental disorders.

In this study, the enzymatic in situ cutting of linearized DNA molecules at approximately 11 kbp intervals is demonstrated using a soft lithography technique. The ultimate goal is to provide a general ordered cutting method to greatly simplify the assembly process. DNA was stretched onto PMMA (Poly methyl methacrylate) coated silicon by withdrawing the substrate from a DNA solution (a process termed “combing”). The stretched lambda DNA could be linearly cut with a soft lithography stamp used to selectively apply DNase I. After cutting the DNA on the substrate, the DNA fragments are removed from the surface by incubating PMMA in the commercial NEBuffer 3.1 at 75°C. The recovered fragments desorbed into the buffer and were sequenced using the PacBio RS II sequencer without an amplification step. The mean coverage was 2870X for the approximately 11 kbp fragmented sample and 100% of the lambda genome was sequenced. Methods to extend of the technique to ordered fragmentation are discussed.

SummaryLong read sequencing platforms, which include the widely used Pacific Biosciences (PacBio) platform and the emerging Oxford Nanopore platform, aim to produce sequence fragments in excess of 15-20 kilobases, and have proved advantageous in the identification of structural variants and easing genome assembly. However, long read sequencing remains relatively expensive and error prone, and failed sequencing runs represent a significant problem for genomics core facilities. To quantitatively assess the underlying mechanics of sequencing failure, it is essential to have highly reproducible and controllable reference data sets to which sequencing results can be compared. Here, we present SiLiCO, the first in silico simulation tool to generate standardized sequencing results from both of the leading long read sequencing platforms.AvailabilitySiLiCO is an open source package written in Python. It is freely available at https://www.github.com/ethanagbaker/SiLiCO under the GNU GPL 3.0 license.ContactSupplementary informationSupplementary data are available at Bioinformatics online.

One of the primary challenges in the field of psychiatric genetics is the lack of anin vivomodel system in which to functionally validate candidateneuropsychiatricriskgenes (NRGs) in a rapid and cost-effective manner1−3. To overcome this obstacle, we performed a candidate-based RNAi screen in whichC. elegansorthologs of human NRGs were assayed for dendritic arborization and cell specification defects usingC. elegansPVD neurons. Of 66 NRGs, identified via exome sequencing of autism (ASD)4or schizophrenia (SCZ)5−9probands and whose mutations arede novoand predicted to result in a complete or partial loss of protein function, theC. elegansorthologs of 7 NRGs were found to be required for proper neuronal development and represent a variety of functional classes, including transcriptional regulators and chromatin remodelers, molecular chaperones, and cytoskeleton-related proteins. Notably, the positive hit rate, when selectively assayingC. elegansorthologs of ASD and SCZ NRGs, is enriched >14-fold as compared to unbiased RNAi screening10. Furthermore, we find that RNAi phenotypes associated with the depletion of NRG orthologs is recapitulated in genetic mutant animals, and, via genetic interaction studies, we show that the NRG ortholog of ANK2,unc-44, is required for SAX-7/MNR-1/DMA-1 signaling. Collectively, our studies demonstrate thatC. elegansPVD neurons are a tractable model in which to discover and dissect the fundamental molecular mechanisms underlying neuropsychiatric disease pathogenesis.

DISC1 as a genetic risk factor for schizophrenia and related major mental illness: response to Sullivan

Apyrases have been suggested to play important roles in plant nutrition, photomorphogenesis, and nodulation. To help trace the evolution of these genes in the legumes--and possibly, the acquisition of new functions for nodulation--apyrase-containing BACs were sequenced from three legume genomes. Genomic sequences from Medicago truncatula, Glycine max and Lotus japonicus were compared to one another and to corresponding regions in Arabidopsis thaliana. A phylogenetic analysis of apyrase homologs from these regions and sequences from other legume species, as well as other plant families, identified a potentially legume-specific clade that contains a well-characterized soybean ( G. soja) apyrase, Gs52, as well as homologs from Dolichos, Lotus, Medicago and Pisum. Sister clades contain homologs from members of Brassicaceae, Solanaceae, Poaceae and Fabaceae. Comparisons of rates of change at synonymous and nonsynonymous sites in the Gs52 and sister clades show rapid evolution in the potentially legume-specific Gs52 clade. The genomic organization of the apyrase-containing BACs shows evidence of gene duplication, genomic rearrangement, and gene conversion among Gs52 homologs. Taken together, these results suggest a scenario of local apyrase gene duplication in an ancestor of the legumes, followed by functional diversification and increased rates of change in the new genes, and further duplications in the Galegae (which include the genera Medicago and Pisum). The study also provides a detailed comparison of genomic regions between two model genomes which are now being sequenced ( M. truncatulaand L. japonicus), and a genome from an economically important legume species ( G. max).

Third generation single molecule sequencing technology is poised to revolutionize genomics by enabling the sequencing of long, individual molecules of DNA and RNA. These technologies now routinely produce reads exceeding 5,000 basepairs, and can achieve reads as long as 50,000 basepairs. Here we evaluate the limits of single molecule sequencing by assessing the impact of long read sequencing in the assembly of the human genome and 25 other important genomes across the tree of life. From this, we develop a new data-driven model using support vector regression that can accurately predict assembly performance. We also present a novel hybrid error correction algorithm for long PacBio sequencing reads that uses pre-assembled Illumina sequences for the error correction. We apply it several prokaryotic and eukaryotic genomes, and show it can achieve near-perfect assemblies of small genomes (< 100Mbp) and substantially improved assemblies of larger ones. All source code and the assembly model are available open-source.

A balanced t(1;11) translocation that transects the Disrupted in schizophrenia 1 (DISC1) gene shows genome-wide significant linkage for schizophrenia and recurrent major depressive disorder (rMDD) in a single large Scottish family, but genome-wide and exome sequencing-based association studies have not supported a role for DISC1 in psychiatric illness. To explore DISC1 in more detail, we sequenced 528 kb of the DISC1 locus in 653 cases and 889 controls. We report 2718 validated single-nucleotide polymorphisms (SNPs) of which 2010 have a minor allele frequency of

The genome project has resulted in an avalanche of data, as well as an avalanche of philosophical revisionism to justify the lack of a hypothesis behind most genomics research. While sometimes heralded as a new age in biological reasoning, little has changed in the way conclusions are reached. What has changed is the richness of the information space within which experiments can be designed. Plants were the first organisms subjected to genetic analysis, and some of the first for which genetic maps assigned functions to chromosomal regions known as genes. The resolution provided by the genome project catapults genetic analysis into the computer age. As in Physics, it is only when significant algorithmic power is brought to bear that Popper's revolution can occur. As at the dawn of genetics, plants have properties particularly favorable to the study of higher eukaryotic genomes. Perhaps the real Green Revolution is only just beginning.Lin Kaul Rounsley Shea Benito Town Fujii Mason Bowman Barnstead 1999xLin, X., Kaul, S., Rounsley, S., Shea, T.P., Benito, M.I., Town, C.D., Fujii, C.Y., Mason, T., Bowman, C.L., Barnstead, M. et al. Nature. 1999; 402: 761–768Crossref | PubMed | Scopus (490)See all References, Wang Stec Hey Lukens Doebley 1999xWang, R.L., Stec, A., Hey, J., Lukens, L., and Doebley, J. Nature. 1999; 18: 236–239See all References

Heterochromatin, constitutively condensed chromosomal material, is widespread among eukaryotes but incompletely characterized at the nucleotide level. We have sequenced and analyzed 2.1 megabases (Mb) of Arabidopsis thaliana chromosome 4 that includes 0.5–0.7 Mb of isolated heterochromatin that resembles the chromosomal knobs described by Barbara McClintock in maize. This isolated region has a low density of expressed genes, low levels of recombination and a low incidence of genetrap insertion. Satellite repeats were absent, but tandem arrays of long repeats and many transposons were found. Methylation of these sequences was dependent on chromatin remodeling. Clustered repeats were associated with condensed chromosomal domains elsewhere. The complete sequence of a heterochromatic island provides an opportunity to study sequence determinants of chromosome condensation.

The higher plant Arabidopsis thaliana (Arabidopsis) is an important model for identifying plant genes and determining their function. To assist biological investigations and to define chromosome structure, a coordinated effort to sequence the Arabidopsis genome was initiated in late 1996. Here we report one of the first milestones of this project, the sequence of chromosome 4. Analysis of 17.38 megabases of unique sequence, representing about 17% of the genome, reveals 3,744 protein coding genes, 81 transfer RNAs and numerous repeat elements. Heterochromatic regions surrounding the putative centromere, which has not yet been completely sequenced, are characterized by an increased frequency of a variety of repeats, new repeats, reduced recombination, lowered gene density and lowered gene expression. Roughly 60% of the predicted protein-coding genes have been functionally characterized on the basis of their homology to known genes. Many genes encode predicted proteins that are homologous to human and Caenorhabditis elegans proteins.

The yeast and human SKP1 genes regulate the mitotic cell cycle but are not yet known to be required for meiosis. Nine Arabidopsis SKP1 homologues have been uncovered and are named ASK1 through ASK9 . Here, we report the isolation and characterization of a male sterile Arabidopsis mutant and show that the mutant defect was caused by a Ds transposon insertion into the ASK1 gene. In the ask1-1 mutant, abnormal microspores exhibit a range of sizes. Furthermore, during mutant male meiosis, although homologous chromosome pairing appeared normal at metaphase I, chromosome segregation at anaphase I is unequal, and some chromosomes are abnormally extended. Therefore, in ask1-1 , at least some homologues remain associated after metaphase I. In addition, immunofluorescence microscopy indicates that the mutant spindle morphology at both metaphase I and early anaphase I is normal; thus, the abnormal chromosome segregation is not likely caused by a spindle defect. Because the yeast Skp1p is required for targeting specific proteins for ubiquitin-mediated proteolysis, we propose that ASK1 controls homologue separation by degrading or otherwise removing a protein that is required directly or indirectly for homologue association before anaphase I.

Tracking data flow in high throughput sequencing is important in maintaining a consistent number of successfully sequenced samples, making decisions on scheduling the flow of sequencing steps, resolving problems at various steps and tracking the status of different projects. This is especially critical when the laboratory is handling a multitude of projects. We have built a Web-based data flow tracking package, called Kaleidaseq, which allows us to monitor the flow and quality of sequencing samples through the steps of preparation of library plates, plaque-picking, preparation of templates, conducting sequencing reactions, loading of samples on gels, base-calling the traces, and calculating the quality of the sequenced samples. Kaleidaseq’s suite of displays allows for outstanding monitoring of the production sequencing process. The online display of current information that Kaleidaseq provides on both project status and process queues sorted by project enables accurate real-time assessment of the necessary samples that must be processed to complete the project. This information allows the process manager to allocate future resources optimally and schedule tasks according to scientific priorities. Quality of the sequenced samples can be tracked on a daily basis, which allows the sequencing laboratory to maintain a steady performance level and quickly resolve dips in quality. Kaleidaseq has a simple easy-to-use interface that allows access to all major functions and process queues from one Web page. This software package is modular and designed to allow additional processing steps and new monitoring variables to be added and tracked with ease. Access to the underlying relational database is through the Perl DBI interface, which allows for the use of different relational databases. Kaleidaseq is available for free use by the academic community from http://www.cshl.org/kaleidaseq.

A cDNA encoding a serotonin receptor has been isolated from a Caenorhabditis elegans mixed stage cDNA library. The nematode serotonin receptor, designated 5HT-Ce, was permanently expressed in murine Ltk-cells, where it mediates adenylate cyclase attenuation. Sequence analysis and the pharmacological profiles demonstrate its relatedness not only to Drosophila and Lymnae 5HT receptors but also to mammalian 5HT1a receptors. The 5HT-Ce-gene does not map close to the position of any known serotonergic mutations.

Genomic instability is one of the hallmarks of cancer, leading to widespread copy number variations, chromosomal fusions, and other structural variations in many cancers. The breast cancer cell line SK-BR-3 is an important model for HER2+ breast cancers, which are among the most aggressive forms of the disease and affect one in five cases. Through short read sequencing, copy number arrays, and other technologies, the genome of SK-BR-3 is known to be highly rearranged with many copy number variations, including an approximately twenty-fold amplification of the HER2 oncogene, along with numerous other amplifications and deletions. However, these technologies cannot precisely characterize the nature and context of the identified genomic events and other important mutations may be missed altogether because of repeats, multi-mapping reads, and the failure to reliably anchor alignments to both sides of a variation. To address these challenges, we have sequenced SK-BR-3 using PacBio long read technology. Using the new P6-C4 chemistry, we generated more than 70X coverage of the genome with average read lengths of 9-13kb (max: 71kb). Using Lumpy for split-read alignment analysis, as well as our novel assembly-based algorithms for finding complex variants, we have developed a detailed map of structural variations in this cell line. Taking advantage of the newly identified breakpoints and combining these with copy number assignments, we have developed an algorithm to reconstruct the mutational history of this cancer genome. From this we have characterized the amplifications of the HER2 region, discovering a complex series of nested duplications and translocations between chr17 and chr8, two of the most frequent translocation partners in primary breast cancers. We have also carried out full-length transcriptome sequencing using PacBio's Iso-Seq technology, which has revealed a number of previously unrecognized gene fusions and isoforms. Combining long-read genome and transcriptome sequencing technologies enables an in-depth analysis of how changes in the genome affect the transcriptome, including how gene fusions are created across multiple chromosomes. This analysis has established the most complete cancer reference genome available to date, and is already opening the door to applying long-read sequencing to patient samples with complex genome structures. Citation Format: Maria Nattestad, Karen Ng, Sara Goodwin, Timour Baslan, Fritz Sedlazeck, James Gurtowski, Elizabeth Hutton, Yogi Sundaravadanam, Tyler Garvin, Marley Alford, Elizabeth Tseng, Philipp Rescheneder, Jason Chin, Timothy Beck, Melissa Kramer, John McPherson, James Hicks, Michael C. Schatz, William R. McCombie. Comprehensive genome and transcriptome structural analysis of a breast cancer cell line using single molecule sequencing. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 850.

The incidence of Barrett’s esophagus (BE)-associated esophageal adenocarcinoma (EAC) is increasing. Next-generation sequencing (NGS) provides an unprecedented opportunity to uncover genomic alterations during BE pathogenesis and progression to EAC, but treatment-naive surgical specimens are scarce. The objective of this study was to establish the feasibility of using widely available endoscopic mucosal biopsies for successful NGS, using samples obtained from a BE ‘progressor’. Paired-end whole-genome NGS was performed on the Illumina platform using libraries generated from mucosal biopsies of normal squamous epithelium (NSE), BE and EAC obtained from a patient who progressed to adenocarcinoma during endoscopic surveillance. Selective validation studies, including Sanger sequencing, immunohistochemistry and functional assays, were performed to confirm the NGS findings. NGS identified somatic nonsense mutations of AT-rich interactive domain 1A (SWI like) (ARID1A) and PPIE and an additional 37 missense mutations in BE and/or EAC, which were confirmed by Sanger sequencing. ARID1A mutations were detected in 15% (3/20) high-grade dysplasia (HGD)/EAC patients. Immunohistochemistry performed on an independent archival cohort demonstrated ARID1A protein loss in 0% (0/76), 4.9% (2/40), 14.3% (4/28), 16.0% (8/50) and 12.2% (12/98) of NSE, BE, low-grade dysplasia, HGD and EAC tissues, respectively, and was inversely associated with nuclear p53 accumulation (P = 0.028). Enhanced cell growth, proliferation and invasion were observed on ARID1A knockdown in EAC cells. In addition, genes downstream of ARID1A that potentially contribute to the ARID1A knockdown phenotype were identified. Our studies establish the feasibility of using mucosal biopsies for NGS, which should enable the comparative analysis of larger ‘progressor’ versus ‘non-progressor’ cohorts. Further, we identify ARID1A as a novel tumor-suppressor gene in BE pathogenesis, reiterating the importance of aberrant chromatin in the metaplasia– dysplasia sequence.

Here, we report genome-wide analysis of the tumor suppressor p53 binding sites in normal human cells. 743 high-confidence ChIP-seq peaks representing putative genomic binding sites were identified in normal IMR90 fibroblasts using a reference chromatin sample. More than 40% were located within 2 kb of a transcription start site (TSS), a distribution similar to that documented for individually studied, functional p53 binding sites and, to date, not observed by previous p53 genome-wide studies. Nearly half of the high-confidence binding sites in the IMR90 cells reside in CpG islands in marked contrast to sites reported in cancer-derived cells. The distinct genomic features of the IMR90 binding sites do not reflect a distinct preference for specific sequences, since the de novo developed p53 motif based on our study is similar to those reported by genome-wide studies of cancer cells. More likely, the different chromatin landscape in normal, compared with cancer-derived cells, influences p53 binding via modulating availability of the sites. We compared the IMR90 ChIP-seq peaks to the recently published IMR90 methylome1 and demonstrated that they are enriched at hypomethylated DNA. Our study represents the first genome-wide, de novo mapping of p53 binding sites in normal human cells and reveals that p53 binding sites reside in distinct genomic landscapes in normal and cancer-derived human cells.

MicroRNAs (miRNAs) are key regulators of gene expression and contribute to a variety of biological processes. Abnormal miRNA expression has been reported in various diseases including pathophysiology of breast cancer, where they regulate protumorigenic processes including vascular invasiveness, estrogen receptor status, chemotherapy resistance, invasion and metastasis. The miRBase sequence database, a public repository for newly discovered miRNAs, has grown rapidly with approximately >10,000 entries to date. Despite this rapid growth, many miRNAs have not yet been validated, and several others are yet to be identified. A lack of a full complement of miRNAs has imposed limitations on recognizing their important roles in cancer, including breast cancer. Using deep sequencing technology, we have identified 189 candidate novel microRNAs in human breast cancer cell lines with diverse tumorigenic potential. We further show that analysis of 500-nucleotide pri-microRNA secondary structure constitutes a reliable method to predict bona fide miRNAs as judged by experimental validation. Candidate novel breast cancer miRNAs with stem lengths of greater than 30 bp resulted in the generation of precursor and mature sequences in vivo. On the other hand, candidates with stem length less than 30 bp were less efficient in producing mature miRNA. This approach may be used to predict which candidate novel miRNA would qualify as bona fide miRNAs from deep sequencing data with approximately 90% accuracy.

A database containing mapped partial cDNA sequences from Caenorhabdhitis elegans will provide a ready starting point for identifying nematode homologues of important human genes and determining their functions in C. elegans. A total of 720 expressed sequence tags (ESTs) have been generated from 585 clones randomly selected from a mixed-stage C. elegans cDNA library. Comparison of these ESTs with sequence databases identified 422 new C. elegans genes, of which 317 are not similar to any sequences in the database. Twenty-six new genes have been mapped by YAC clone hybridization. Members of several gene families, including cuticle collagens, GTP-binding proteins, and RNA helicases were discovered. Many of the new genes are similar to known or potential human disease genes, including CFTR and the LDL receptor.

Automated DNA sequencing is an extremely valuable technique which requires very high quality DNA templates to be carried out successfully. While it has been possible to readily produce large numbers of such templates from M13 or other single-stranded vectors for several years, the sequencing of double-stranded DNA templates using the ABI 373 DNA Sequencer has had a considerably lower success rate. We describe how the combination of a new fluorescent, dideoxy sequencing method, called cycle-sequencing, coupled with modifications to template isolation procedures based on Qiagen columns, makes fluorescent sequencing of double-stranded templates a reliable procedure. From a single five milliliter culture enough DNA can be isolated (up to 20 micrograms) to do 4-8 sequencing reactions, each of which yields 400-500 bases of high quality sequence data. These procedures make the routine use of double-stranded DNA templates a viable strategy in automated DNA sequencing projects.

The sequence of five non-contiguous genomic fragments encompassing 14.4 kilobases from the NF1 locus have been determined by fluorescence-based automated DNA sequence analysis. These fragments included one kilobase of the NF1 coding region, which resulted in the identification of the intron/exon boundaries of five exons. Based on these sequences, five new NF1 exon-PCR assays have been developed, that could be useful for detecting new NF1 mutations. The genomic sequences were analyzed for the presence of Alu repetitive elements and their classification is described. This analysis may provide some insight into the characterization of genetic rearrangements resulting in disruption of the NF1 gene.

The PRNP gene, encoding the amyloid precursor protein that is centrally involved in Creutzfeldt-Jakob disease (CJD), has an unstable region of five variant tandem octapeptide coding repeats between codons 51 and 91. We screened a total of 535 individuals for the presence of extra repeats in this region, including patients with sporadic and familial forms of spongiform encephalopathy, members of their families, other neurological and non-neurological patients, and normal controls. We identified three CJD families (in each of which the proband's disease was neuropathologically confirmed and experimentally transmitted to primates) that were heterozygous for alleles with 10, 12, or 13 repeats, some of which had "wobble" nucleotide substitutions. We also found one individual with 9 repeats and no nucleotide substitutions who had no evidence of neurological disease. These observations, together with data on published British patients with 11 and 14 repeats, strongly suggest that the occurrence of 10 or more octapeptide repeats in the encoded amyloid precursor protein predisposes to CJD.

Background Ginkgo biloba L. is the only surviving member of one of the oldest living seed plant groups with medicinal, spiritual and horticultural importance worldwide. As an evolutionary relic, it displays many characters found in the early, extinct seed plants and extant cycads. To establish a molecular base to understand the evolution of seeds and pollen, we created a cDNA library and EST dataset from the reproductive structures of male (microsporangiate), female (megasporangiate), and vegetative organs (leaves) of Ginkgo biloba. Results RNA from newly emerged male and female reproductive organs and immature leaves was used to create three distinct cDNA libraries from which 6,434 ESTs were generated. These 6,434 ESTs from Ginkgo biloba were clustered into 3,830 unigenes. A comparison of our Ginkgo unigene set against the fully annotated genomes of rice and Arabidopsis, and all available ESTs in Genbank revealed that 256 Ginkgo unigenes match only genes among the gymnosperms and non-seed plants – many with multiple matches to genes in non-angiosperm plants. Conversely, another group of unigenes in Gingko had highly significant homology to transcription factors in angiosperms involved in development, including MADS box genes as well as post-transcriptional regulators. Several of the conserved developmental genes found in Ginkgo had top BLAST homology to cycad genes. We also note here the presence of ESTs in G. biloba similar to genes that to date have only been found in gymnosperms and an additional 22 Ginkgo genes common only to genes from cycads. Conclusion Our analysis of an EST dataset from G. biloba revealed genes potentially unique to gymnosperms. Many of these genes showed homology to fully sequenced clones from our cycad EST dataset found in common only with gymnosperms. Other Ginkgo ESTs are similar to developmental regulators in higher plants. This work sets the stage for future studies on Ginkgo to better understand seed and pollen evolution, and to resolve the ambiguous phylogenetic relationship of G. biloba among the gymnosperms.

We have initially sequenced approximately 8,000 canine expressed sequence tags (ESTs) from several complementary DNA (cDNA) libraries: testes, whole brain, and Madin-Darby canine kidney (MDCK) cells. Analysis of these sequences shows that they provide partial sequence information for about 5%-10% of the canine genes. An analysis pipeline has been created to cluster the ESTs and to map individual ESTs as well as clustered ESTs to both the human genome and the human proteome. Gene ontology (GO) terms have been assigned to the ESTs and clusters based on their top matches to the International Protein Index (IPI) set of human proteins. The data generated is stored in a MySQL relational database for analysis and display. A Web-based Perl script has been written to display the analyzed data to the scientific community.

10 páginas, 3 figuras, 8 tablas.-- et al., We have sequenced and annotated the genome of fission yeast (Schizosaccharomyces pombe), which contains the smallest number of protein-coding genes yet recorded for a eukaryote: 4,824. The centromeres are between 35 and 110 kilobases (kb) and contain related repeats including a highly conserved 1.8-kb element. Regions upstream of genes are longer than in budding yeast (Saccharomyces cerevisiae), possibly reflecting more-extended control regions. Some 43% of the genes contain introns, of which there are 4,730. Fifty genes have significant similarity with human disease genes; half of these are cancer related. We identify highly conserved genes important for eukaryotic cell organization including those required for the cytoskeleton, compartmentation, cell-cycle control, proteolysis, protein phosphorylation and RNA splicing. These genes may have originated with the appearance of eukaryotic life. Few similarly conserved genes that are important for multicellular organization were identified, suggesting that the transition from prokaryotes to eukaryotes required more new genes than did the transition from unicellular to multicellular organization., We thank the European Commission, the Wellcome Trust and Cancer Research UK for financial support.

The genome of the model plant Arabidopsis thaliana has been sequenced by an international collaboration, The Arabidopsis Genome Initiative. Here we report the complete sequence of chromosome 5. This chromosome is 26 megabases long; it is the second largest Arabidopsis chromosome and represents 21% of the sequenced regions of the genome. The sequence of chromosomes 2 and 4 have been reported previously and that of chromosomes 1 and 3, together with an analysis of the complete genome sequence, are reported in this issue. Analysis of the sequence of chromosome 5 yields further insights into centromere structure and the sequence determinants of heterochromatin condensation. The 5,874 genes encoded on chromosome 5 reveal several new functions in plants, and the patterns of gene organization provide insights into the mechanisms and extent of genome evolution in plants.

The genomes of higher plants and animals are highly differentiated, and are composed of a relatively small number of genes and a large fraction of repetitive DNA. The bulk of this repetitive DNA constitutes transposable, and especially retrotransposable, elements. It has been hypothesized that most of these elements are heavily methylated relative to genes, but the evidence for this is controversial. We show here that repeat sequences in maize are largely excluded from genomic shotgun libraries by the selection of an appropriate host strain because of their sensitivity to bacterial restriction-modification systems. In contrast, unmethylated genic regions are preserved in these genetically filtered libraries if the insert size is less than the average size of genes. The representation of unique maize sequences not found in plant reference genomes is also greatly enriched. This demonstrates that repeats, and not genes, are the primary targets of methylation in maize. The use of restrictive libraries in genome shotgun sequencing in plant genomes should allow significant representation of genes, reducing the number of reactions required.

Arabidopsis thaliana has emerged as a model system for studies of plant genetics and development, and its genome has been targeted for sequencing1 by an international consortium (the Arabidopsis Genome Initiative; http://genome-www.stanford.edu/Arabidopsis/agi.html ). To support the genome-sequencing effort, we fingerprinted more than 20,000 BACs (ref. 2) from two high-quality publicly available libraries3,4,5, generating an estimated 17-fold redundant coverage of the genome, and used the fingerprints to nucleate assembly of the data by computer. Subsequent manual revision of the assemblies resulted in the incorporation of 19,661 fingerprinted BACs into 169 ordered sets of overlapping clones ('contigs'), each containing at least 3 clones. These contigs are ideal for parallel selection of BACs for large-scale sequencing and have supported the generation of more than 5.8 Mb of finished genome sequence submitted to GenBank; analysis of the sequence has confirmed the integrity of contigs constructed using this fingerprint data. Placement of contigs onto chromosomes can now be performed, and is being pursued by groups involved in both sequencing and positional cloning studies. To our knowledge, these data provide the first example of whole-genome random BAC fingerprint analysis of a eucaryote, and have provided a model essential to efforts aimed at generating similar databases of fingerprint contigs to support sequencing of other complex genomes, including that of human.

A fluorescence-based, T7 (Sequenase) dye terminator method for sequencing double-stranded DNA using strings of three contiguous hexamers as primers and single-stranded binding protein is described. In this method, the circular, supercoiled DNA vector pUC19 is first linearized with a restriction enzyme to create a sequenceable template. Sequencing is then accomplished using three cycles of "denaturation," annealing, and extension/termination. Twenty-two of 33 hexamer strings tested in a controlled study produced acceptable sequence, with read lengths varying from the mid 300s to the low 400s and a base-calling accuracy of at least 97%. To test its potential utility in directed DNA sequencing, the protocol was then used to completely sequence both strands of pUC19. For this test project, a total of 28 hexamer strings was used with an overall successful priming rate of 75%. The current protocol appears to be sufficiently robust to be used in the finishing phase of a shotgun sequencing project and is amenable to semiautomation. Prospects for using the protocol for full-scale directed sequencing as well as for full automation are discussed.

The finishing phase of genome sequencing projects is expensive, in part, because of the cost of de novo synthesis of custom primers and the management burden associated with obtaining and using them for primer walking. One approach to reduce these high costs is the use of a presynthesized library of short oligonucleotides (8-10 bases) rather than long primers. The use of such a library eliminates the need for custom synthesis of oligonucleotides, providing the convenience of priming from any site by combining two to three short oligonucleotides to form a string with the required specificity. The first practical implementation of this strategy presented a robust protocol for using hexamer strings with radioisotopic labelling. Whereas versions of this technique have subsequently been implemented on fluorescent sequencers we felt that there was a need to develop and extensively test a protocol that consistently gave read lengths comparable to dye-terminator sequencing with longer primers. We have developed a new two-cycle fluorescent Sequenase terminator procedure for using hexamer strings. We tested this procedure using a set of 32 different 3 hexamer primer strings, each known to be functional to some degree in radioisotopic sequencing, on single-stranded M13mp18 template and ABI 373 DNA sequencers. The overall success rate of priming with these hexamer primer strings is 97% with the failure of only one string. In this case, the corresponding 18-mer primer also failed to produce usable sequence from M13mp18 template. The average read length from reactions successfully primed with the 31 different hexamer strings was 461 bases with > 99% base-calling accuracy. The current protocol is robust enough to be used in virtually any situation where primer walking on single-stranded templates is used. The success rate and read lengths make it universally applicable to the sequencing of single-stranded templates on automated sequencers. It is also amenable to automation.

Gene trap transposon mutagenesis can identify essential genes whose functions in later development are obscured by an early lethal phenotype. In higher plants, many genes are required for haploid gametophyte viability, so that the phenotypic effects of their disruption cannot be readily observed in the diploid plant body. The PROLIFERA (PRL) gene, identified by gene trap transposon mutagenesis in Arabidopsis, is required for megaga-metophyte and embryo development. Reporter gene expression patterns revealed that PRL was expressed in dividing cells throughout the plant. PRL is related to the MCM2-3-5 family of yeast genes that are required for the initiation of DNA replication.

179 Background: Although designed to retain activity against known androgen receptor (AR) mutants, it is predicted that new MDV3100-resistant mutations may predict for sensitivity to the drug in clinic. To predict for tumor sensitivity to treatment with novel AR targeted therapies, we explored the frequency of mutation detection and copy number alteration in CTC isolated from patients with CRPC enrolled on a trial with MDV3100. Methods: We used fluorescence-activated cell sorting (FACS) methodology to enrich EpCAM+, CD45−, DAPI− cells. For mutation detection and genomic copy number alteration in CTC in low number of cancer cells found in clinical samples, we optimized next-gen deep sequencing by Illumina. Results: In patients with progressive CRPC, >10 or >50 EpCAM+ events (EPE) were isolated by FACS in 88% or 58% of patients, in whom 32% and 10% had unfavorable (>5 cells/7.5 ml) CTC counts using CellSearch. EPE, expressing prostate-specific mRNAs, provide sufficient high quality DNA for genomic sequencing and copy number analysis. Adequate coverage was obtained from as few 50 EPE, with a recovery rate of 89% from FACS sorted samples. The detection threshold of a mutation was established at 1:4 alleles. To further expand genomic profiling in CTC, we optimized Nimblegen exome mutation detection by deep sequencing on HiSeq PE101. Our initial analysis established the polymorphism frequency detection thresholds in heterogeneous cell populations, and confirmed the sequencing coverage. Somatic missense mutations in AR, APC and TP53 found in CTC but not in paired WBC were confirmed by Sanger sequencing. In parallel, copy number alterations in CTC are studied. Conclusions: Somatic mutations detected in CTC isolated from patients with CRPC can serve as predictive markers of tumor sensitivity to targeted therapies. We established standard operating procedures for specimen processing, and confirmed the sequencing coverage and polymorphism detection thresholds in heterogeneous cell population. Currently we are proceeding to clinical samples to study the associations between specific molecular alterations in CTC as predictive markers of sensitivity and clinical outcomes.

The use of strings of contiguous short primers rather than a single long primer has the potential to greatly diminish the cost, time and management efforts associated with sequencing by primer walking. For maximum impact, this chemistry must be adaptable to current fluorescent automated sequencers, which, while allowing automated acquisition of data, have reduced sensitivity when compared with radioactive-based sequencing. The ability to use hexamer strings on Applied Biosystems DNA sequencers has now been demonstrated for single-stranded templates. Procedures are described that allow sequence to be obtained up to 400 bases from the priming site. Signal strength is sufficient in this region to allow single base resolution and highly accurate automatic base calling to be performed by the sequencer. While these conditions can no doubt be further optimized, these results show the feasibility of inexpensive primer walking using hexamer string primers on currently available commercial DNA sequencers. This should have a wide range of applications from genome sequencing projects to the sequencing of cDNA clones without the necessity of creating nested deletions or the necessity of spending inordinate amounts of time and money on oligonucleotide synthesis.

4540 Background: CTC are representative of the genomic makeup of a tumor and can provide the molecular snapshot of the disease at the time treatment is being considered. We explored the frequency o...

38 Background: To estimate the association between molecular biomarkers detected in circulating tumor cells (CTC) and tumor sensitivity to treatment, robust assays are needed before qualification in prospective clinical trials. Methods: To address the limitations of the current FDA cleared technology, we focused on improving our ability to isolate more purified CTC populations based on fluorescence-activated cell sorting (FACS) to capture EpCAM+, CD45−, DAPI− cells from patients with castration-resistant prostate cancer (CRPC). Androgen receptor (AR) and genes frequently mutated in CRPC have been selected from the integrative genomic profiling at MSK. We optimized the RainDance microfluidic PCR followed by targeted sequencing in low number of cancer cells, before proceeding to clinical samples. Results: On blood samples from124 patients with progressive CRPC, FACS method isolates an average 100-fold more EpCAM+ events compared to the current FDA cleared CellSearch assay. By FACS, >10 or >50 events were isolated in 88% or 58% of patients, compared to 32% or 10% of patients by CellSearch, respectively. FACS isolated cells express prostate-specific mRNAs (PSA, AR, TMPRSS2-ERG), as detected by an analytically validated multiplex RT-PCR, indicating that these EpCAM+ events are bona fide CTC. For genomic profiling, sufficient high quality DNA was obtained from as little as 50 CTC, with a recovery rate of 89% from FACS sorted samples and adequate sequencing coverage, and 1:4 detection threshold in a heterogeneous cell population. Selected missense mutations in AR, PIK3CA and TP53 found in CTC but not in WBC from same patient are further analyzed. Conclusions: Molecular alterations in CTC can potentially serve as predictive markers of sensitivity and clinical outcomes as surrogate tissue in clinical practice. We established standard operating procedures for specimen processing, and confirmed the sequencing coverage and polymorphism detection thresholds in a heterogeneous cell population. In the context of available samples collected from patients enrolled on AR-targeted therapies, we will generate data to qualify CTC as biomarkers under the Oncology Biomarker Qualification Initiative. [Table: see text]

We report the finding of an insert mutation in the chromosome 20 amyloid precursor gene in a family with neuropathologically-verified, experimentally-transmitted Gerstmann-Straussler-Scheinker syndrome (GSS). The insert consisted of 8 extra copies of a repeating octapeptide coding sequence in the region between codons 51 and 91; it was identified in the proband and a presently unaffected at-risk niece by full sequencing of the open reading frame, and was visualized electrophoretically in the proband and 6 of 12 at-risk relatives. Although affected members in this French-Breton family have shown a variety of clinical profiles, including durations of illness that ranged from 3 months to 13 years, all autopsied cases (including the patient with the shortest illness) have had the distinctive multicentric amyloid plaques that define GSS as a nosologic entity.

A total of 116,118 basepairs (bp) derived from three cosmids spanning the ERCC1 locus of human chromosome 19q13.3 have been sequenced with automated fluorescence-based sequencers and analysed by polymerase chain reaction amplification and computer methods. The assembled sequence forms two contigs totalling 105,831 bp, which contain a human fosB proto-oncogene, a gene encoding a protein phosphatase, two genes of unknown function and the previously-characterized ERCC1 DNA repair gene. This light band region has a high average density of 1.4 Alu repeats per kilobase. Human chromosome light bands could therefore contain up to 75,000 genes and 1.5 million Alu repeats.

An American family of English origin with an unusually early onset and long-duration form of Creutzfeldt-Jakob disease (CJD) had a heterozygous insert mutation in the region of repeating octapeptide coding sequences between codons 51 and 91 of the PRNP gene on chromosome 20. Affected members were 23 to 35 years old at the onset of illnesses that lasted from 4 to 13 years, yet experimental transmission of disease from the proband (11-year duration) produced a typically brief incubation period and duration of illness in each of three inoculated primates. Also, the PrP amyloid protein that accumulates in CJD brain was only barely detectable in extracted brain tissue from one case with massive spongiform change and was undetectable in another case with no spongiform change, perhaps because of epitope shielding by a configurational change in the protein induced by the mutation. Analysis of this and other families with similar inserts suggests that such mutations in the PRNP gene not only predispose to CJD, but also modify its phenotypic expression.

Clones encoding an atypical beta-adrenergic receptor were isolated from a rat brown adipose tissue cDNA library. This receptor expressed in Chinese hamster ovary (CHO) cells displays a low affinity for beta-adrenergic antagonists and a high affinity for BRL 37344, an agonist that selectively stimulates lipolysis in adipose tissue. The rank order of potency for agonist-mediated increases in intracellular cAMP in transfected cells correlates with that for agonist-mediated stimulation of lipolysis in brown adipocytes. Northern blot analysis demonstrates that this receptor subtype is expressed only in brown and white adipose tissue where it represents the predominant beta-receptor subtype. The amount of atypical beta-adrenergic receptor present in adipose tissue of obese (fa/fa) Zucker rats is reduced by up to 71% as compared with lean (Fa/Fa) control animals. These findings suggest that a change in the expression of this beta-adrenergic receptor subtype may play a role in obesity.