Rare disruptive variants in the DISC1 Interactome and Regulome: association with cognitive ability and schizophrenia

Schizophrenia (SCZ), bipolar disorder (BD) and recurrent major depressive disorder (rMDD) are common psychiatric illnesses. All have been associated with lower cognitive ability, and show evidence of genetic overlap and substantial evidence of pleiotropy with cognitive function and neuroticism. Disrupted in schizophrenia 1 (DISC1) protein directly interacts with a large set of proteins (DISC1 Interactome) that are involved in brain development and signaling. Modulation of DISC1 expression alters the expression of a circumscribed set of genes (DISC1 Regulome) that are also implicated in brain biology and disorder. Here we report targeted sequencing of 59 DISC1 Interactome genes and 154 Regulome genes in 654 psychiatric patients and 889 cognitively-phenotyped control subjects, on whom we previously reported evidence for trait association from complete sequencing of the DISC1 locus. Burden analyses of rare and singleton variants predicted to be damaging were performed for psychiatric disorders, cognitive variables and personality traits. The DISC1 Interactome and Regulome showed differential association across the phenotypes tested. After family-wise error correction across all traits (FWERacross), an increased burden of singleton disruptive variants in the Regulome was associated with SCZ (FWERacross P=0.0339). The burden of singleton disruptive variants in the DISC1 Interactome was associated with low cognitive ability at age 11 (FWERacross P=0.0043). These results identify altered regulation of schizophrenia candidate genes by DISC1 and its core Interactome as an alternate pathway for schizophrenia risk, consistent with the emerging effects of rare copy number variants associated with intellectual disability.

708 Common and 2010 rare DISC1 locus variants identified in 1542 subjects: analysis for association with psychiatric disorder and cognitive traits

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 <1%. Only 38% of these variants are reported in the 1000 Genomes Project European subset. This suggests that many DISC1 SNPs remain undiscovered and are essentially private. Rare coding variants identified exclusively in patients were found in likely functional protein domains. Significant region-wide association was observed between rs16856199 and rMDD (P=0.026, unadjusted P=6.3 × 10(-5), OR=3.48). This was not replicated in additional recurrent major depression samples (replication P=0.11). Combined analysis of both the original and replication set supported the original association (P=0.0058, OR=1.46). Evidence for segregation of this variant with disease in families was limited to those of rMDD individuals referred from primary care. Burden analysis for coding and non-coding variants gave nominal associations with diagnosis and measures of mood and cognition. Together, these observations are likely to generalise to other candidate genes for major mental illness and may thus provide guidelines for the design of future studies.

Next-generation sequencing of endoscopic biopsies identifies ARID1A as a tumor-suppressor gene in Barrett's esophagus

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.

Distinct p53 genomic binding patterns in normal and cancer-derived human cells

We report here 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 ChIPseq peaks to the recently published IMR90 methylome 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.

Genomic analysis of circulating tumor cells to evaluate predictive biomarkers

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]

EST analysis in Ginkgo biloba: an assessment of conserved developmental regulators and gymnosperm specific genes

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.

EST analysis in Ginkgo biloba: an assessment of conserved developmental regulators and gymnosperm specific genes

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.

Evolution and microsynteny of the apyrase gene family in three legume genomes

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).

A survey of canine expressed sequence tags and a display of their annotations through a flexible web-based interface

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.

The genome sequence of Schizosaccharomyces pombe

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.

Initial sequencing and analysis of the human genome

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

Identification of mesoderm development (mesd) candidate genes by comparative mapping and genome sequence analysis

The proximal albino deletions identify several functional regions on mouse Chromosome 7 critical for differentiation of mesoderm (mesd), development of the hypothalamus neuroendocrine lineage (nelg), and function of the liver (hsdr1). Using comparative mapping and genomic sequence analysis, we have identified four novel genes and Il16 in the mesd deletion interval. Two of the novel genes, mesdc1 and mesdc2, are located within the mesd critical region defined by BAC transgenic rescue. We have investigated the fetal role of genes located outside the mesd critical region using BAC transgenic complementation of the mesd early embryonic lethality. Using human radiation hybrid mapping and BAC contig construction, we have identified a conserved region of human chromosome 15 homologous to the mesd, nelg, and hsdr1 functional regions. Three human diseases cosegregate with microsatellite markers used in construction of the human BAC/YAC physical map, including autosomal dominant nocturnal frontal lobe epilepsy (ENFL2; also known as ADNFLE), a syndrome of mental retardation, spasticity, and tapetoretinal degeneration (MRST); and a pyogenic arthritis, pyoderma gangrenosum, and acne syndrome (PAPA).

Sequence and analysis of chromosome 5 of the plant Arabidopsis thaliana

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.

Genetic definition and sequence analysis of Arabidopsis centromeres

High-precision genetic mapping was used to define the regions that contain centromere functions on each natural chromosome in Arabidopsis thaliana. These regions exhibited dramatic recombinational repression and contained complex DNA surrounding large arrays of 180-base pair repeats. Unexpectedly, the DNA within the centromeres was not merely structural but also encoded several expressed genes. The regions flanking the centromeres were densely populated by repetitive elements yet experienced normal levels of recombination. The genetically defined centromeres were well conserved among Arabidopsis ecotypes but displayed limited sequence homology between different chromosomes, excluding repetitive DNA. This investigation provides a platform for dissecting the role of individual sequences in centromeres in higher eukaryotes.

Sequence and analysis of chromosome 4 of the plant Arabidopsis thaliana

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.

Differential methylation of genes and retrotransposons facilitates shotgun sequencing of the maize genome

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.

The Arabidopsis SKP1-LIKE1 gene is essential for male meiosis and may control homologue separation

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.

zA map for sequence analysis of the Arabidopsis thaliana genome

Arabidopsis thaliana has emerged as a model system for studies of plant genetics and development, and its genome has been targeted for sequencing 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 libraries, 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.

Kaleidaseq: a Web-based tool to monitor data flow in a high throughput sequencing facility

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.

Molecular cloning and functional expression of a serotonin receptor from Caenorhabditis elegans

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.

Fluorescence-based sequencing of double-stranded DNA by hexamer string priming

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.

High-quality automated DNA sequencing primed with hexamer strings

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 tagging of PROLIFERA, an essential MCM2-3-5-like gene in Arabidopsis

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.

Automated DNA sequencing using 4-color fluorescent detection of reactions primed with hexamer strings

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.

Iatrogenic Creutzfeldt-Jakob disease: an example of the interplay between ancient genes and modern medicine

We tested DNA from 15 centrally infected cases of iatrogenic Creutzfeldt-Jakob disease (CJD) (dura mater or corneal homografts and stereotactic EEG electrodes), 11 peripherally infected cases (native human growth hormone or gonadotrophin), and 110 control individuals for the presence of mutations in the chromosome 20 amyloid gene. No patient or control had any of the known pathogenic point or insert mutations found in familial disease, but allelic homozygosity at polymorphic codon 129 was present in all but two (92%) of the 26 patients, compared with 54 (50%) of the 110 controls (p < 0.001). Pooled data from all identified and tested cases of iatrogenic disease yielded a worldwide total of 56 patients, of whom all but four were homozygous at codon 129 (p < 0.001). These findings support the thesis that homozygosity at codon 129 enhances susceptibility to iatrogenic infections of both central and peripheral origin, with evident implications for the population of dura mater homograft and pituitary hormone recipients whose lives have been complicated by the possibility of exposure to the infectious agent of CJD.

Rapid selection using G418 of high copy number transformants of Pichia pastoris for high-level foreign gene expression

Pichia pastoris is a methylotrophic yeast increasingly important in the production of therapeutic proteins. Expression vectors are based on the methanol-inducible AOX1 promoter and are integrated into the host chromosome. In most cases high copy number integration has been shown to be important for high-level expression. Since this occurs at low frequency during transformation, we previously used DNA dot blot screens to identify suitable clones. In this paper we report the use of vectors containing the Tn903 kanr gene conferring G418-resistance. Initial experiments demonstrated that copy number showed a tight correlation with drug-resistance. Using a G418 growth inhibition screen, we readily isolated a series of transformants, containing progressively increasing numbers (1 to 12) of a vector expressing HIV-1 ENV, which we used to examine the relationship between copy number and foreign mRNA levels. Northern blot analysis indicated that ENV mRNA levels from a single-copy clone were nearly as high as AOX1 mRNA, and increased progressively with increasing copy number so as to greatly exceed AOX1 mRNA. We have also developed protocols for the selection, using G418, of high copy number transformants following spheroplast transformation or electroporation. We anticipate that these protocols will simplify the use of Pichia as a biotechnological tool.

An insert mutation in the chromosome 20 amyloid precursor gene in a Gerstmann-Sträussler-Scheinker family

We report the finding of an insert mutation in the chromosome 20 amyloid precursor gene in a family with neuropathologically-verified, experimentally-transmitted Gerstmann-Sträussler-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.

Expressed genes, Alu repeats and polymorphisms in cosmids sequenced from chromosome 4p16.3

The sequences of three cosmids (90 kilobases) from the Huntington's disease region in chromosome 4p16.3 have been determined. A 30,837 base overlap of DNA sequenced from two individuals was found to contain 72 DNA sequence polymorphisms, an average of 2.3 polymorphisms per kilobase (kb). The assembled 58 kb contig contains 62 Alu repeats, and eleven predicted exons representing at least three expressed genes that encode previously unidentified proteins. Each of these genes is associated with a CpG island. The structure of one of the new genes, hda1-1, has been determined by characterizing cDNAs from a placental library. This gene is expressed in a variety of tissues and may encode a novel housekeeping gene.

Caenorhabditis elegans expressed sequence tags identify gene families and potential disease gene homologues

A database containing mapped partial cDNA sequences from Caenorhabditis 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 and analysis of 106 kilobases from human chromosome 19q13.3

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.

Creutzfeldt-Jakob disease cosegregates with the codon 178AsnPRNP mutation in families of European origin

We recently discovered an amino acid-altering heterozygous mutation in codon 178 of the PRNP amyloid precursor gene in patients with familial Creutzfeldt-Jakob disease. This mutation is now shown to be associated with the occurrence of disease in 7 unrelated families of Western European origin, among which a total of 65 members are known to have died from Creutzfeldt-Jakob disease. The mutation was detected in each of 17 tested patients, including at least 1 affected member of each family, and in 16 of 36 of their first-degree relatives, but not in affected families with other mutations, patients with the nonfamilial form of the disease, or 83 healthy control individuals. Linkage analysis in two informative families yielded a lod score of 5.30, which, because no recombinants were found, strongly suggests that codon 178Asn is the actual disease mutation.

Atypical Creutzfeldt-Jakob disease in an American family with an insert mutation in the PRNP amyloid precursor gene

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.

Sequencing and analysis of genomic fragments from the NF1 locus

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.

Rapid and reliable fluorescent cycle sequencing of double-stranded templates

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.

An adipose tissue-specific beta-adrenergic receptor. Molecular cloning and down-regulation in obesity

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.

Transmissible familial Creutzfeldt-Jakob disease associated with five, seven, and eight extra octapeptide coding repeats in the PRNP 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.

Cloning, localization, and permanent expression of a Drosophila octopamine receptor

A cDNA for a member of the G protein-coupled receptor family was isolated from Drosophila using a probe derived from a human beta 2-adrenergic receptor cDNA. This Drosophila receptor gene is localized at 99A10-B1 on the right arm of chromosome 3 and is preferentially expressed in Drosophila heads. The insect octopamine receptor has been permanently expressed in mammalian cells, where it mediates the attenuation of adenylate cyclase activity and exhibits a pharmacological profile consistent with an octopamine type 1 receptor. Sequence and pharmacological comparisons indicate that the octopamine receptor is unique but closely related to mammalian adrenergic receptors, perhaps as an evolutionary precursor.

Cloning, sequence analysis and chromosome localization of a Drosophila muscarinic acetylcholine receptor

Two cDNA clones (3.7 kb and 4.8 kb) encoding a Drosophila muscarinic acetylcholine receptor were isolated from a Drosophila head cDNA library and characterized by automated DNA sequence analysis. The Drosophila muscarinic receptor contains 788 amino acids with a calculated Mr of 84,807 and displays greater than 60% homology with mammalian muscarinic receptors. The muscarinic receptor maps to the tip of the right arm of the second chromosome of the Drosophila genome.

Cloning, sequence analysis, and permanent expression of a human alpha 2-adrenergic receptor in Chinese hamster ovary cells. Evidence for independent pathways of receptor coupling to adenylate cyclase attenuation and activation

The gene encoding a human alpha 2-adrenergic receptor was isolated from a human genomic DNA library using a 367-base pair fragment of Drosophila genomic DNA that exhibited 54% identity with the human beta 2-adrenergic receptor and 57% identity with the human alpha 2-adrenergic receptor. The nucleotide sequence of a fragment containing the human alpha 2-receptor gene and 2.076 kilobases of untranslated 5' sequence was determined, and potential upstream regulatory regions were identified. This gene encodes a protein of 450 amino acids and was identified as an alpha 2-adrenergic receptor by homology with published sequences and by pharmacological characterization of the protein expressed in cultured cells. Permanent expression of the alpha 2-receptor was achieved by transfecting Chinese hamster ovary (CHO) cells which lack adrenergic receptors with a 1.5-kilobase NcoI-HindIII fragment of the genomic clone containing the coding region of the gene. The alpha 2-receptor expressed in CHO cells displayed pharmacology characteristic of an alpha 2 A-receptor subtype with a high affinity for yohimbine (Ki = 1 nM) and a low affinity for prazosin (Ki = 10,000 nM). Agonists displayed a rank order of potency in radioligand binding assays of para-aminoclonidine greater than or equal to UK-14304 greater than (-)-epinephrine greater than (-)-norepinephrine greater than (-)-isoproterenol, consistent with the identification of this protein as an alpha 2-receptor. The role of the alpha 2-receptor in modulating intracellular cyclic AMP concentrations was investigated in three transfected cell lines expressing 50, 200, and 1200 fmol of receptor/mg membrane protein. At low concentrations (1-100 nM), (-)-epinephrine attenuated forskolin-stimulated cyclic AMP accumulation by up to 60% in a receptor density-dependent manner. At epinephrine concentrations above 100 nM, cyclic AMP levels were increased up to 140% of the forskolin-stimulated level. Pertussis toxin pretreatment of cells eliminated alpha 2-receptor-mediated attenuation of forskolin-stimulated cyclic AMP levels and enhanced the receptor density-dependent potentiation of forskolin-stimulated cyclic AMP concentrations from 3 to 8-fold. Potentiation of forskolin-stimulated cyclic AMP levels was also elicited by the alpha 2-adrenergic agonists, UK-14304 and para-aminoclonidine, and blocked by the alpha 2-adrenergic antagonist yohimbine, but not by the alpha 1-adrenergic antagonist prazosin or the beta-adrenergic antagonist propranolol.(ABSTRACT TRUNCATED AT 400 WORDS)

Toluene degradation by Pseudomonas putida F1: genetic organization of the tod operon

Pseudomonas putida PpF1 degrades toluene through cis-toluene dihydrodiol to 3-methylcatechol. The latter compound is metabolized through the well-established meta pathway for catechol degradation. The first four steps in the pathway involve the sequential action of toluene dioxygenase (todABC1C2), cis-toluene dihydrodiol dehydrogenase (todD), 3-methylcatechol 2,3-dioxygenase (todE), and 2-hydroxy-6-oxo-2,4-heptadienoate hydrolase (todF). The genes for these enzymes form part of the tod operon which is responsible for the degradation of toluene by this organism. A combination of transposon mutagenesis of the PpF1 chromosome, as well as analysis of cloned chromosomal fragments, was used to determine the physical order of the genes in the tod operon. The genes were determined to be transcribed in the order todF, todC1, todC2, todB, todA, todD, todE.

High level expression of heterologous proteins in methylotrophic yeast Pichia pastoris

Expression of human tumor necrosis factor-alpha (TNF) and four different TNF analogs has been studied in Pichia pastoris by utilizing the alcohol oxidase gene promoter. TNF expression level in certain transformants accounted for as much as 36% of the soluble protein. TNF expression was stably maintained during high cell density fermentation (100 g dry cell weight/liter) resulting in a TNF production level of 6-10 g/liter. TNF contained in P. pastoris cell lysates was biologically active as determined by its cytotoxic effect on murine L-929 fibroblast cells and the bioactivity was retained for at least 6 months in the lysates stored frozen at -20 degrees C in the presence of protease inhibitor PMSF. TNF expressed in P. pastoris was recognized by monoclonal antibodies prepared against recombinant Escherichia coli derived TNF. TNF purified from P. pastoris has the expected N-terminal amino acid sequence and specific activity of 10(7) units/mg protein. TNF analogs were also expressed at levels comparable to that of native TNF. Three of the four analogs were insoluble when produced in P. pastoris.

Identification of cis-diols as intermediates in the oxidation of aromatic acids by a strain of Pseudomonas putida that contains a TOL plasmid

Pseudomonas putida BG1 was isolated from soil by enrichment with p-toluate and selection for growth with p-xylene. Other hydrocarbons that served as growth substrates were toluene, m-xylene, 3-ethyltoluene, and 1,2,4-trimethylbenzene. The enzymes responsible for growth on these substrates are encoded by a large plasmid with properties similar to those of TOL plasmids isolated from other strains of Pseudomonas. Treatment of P. putida BG1 with nitrosoguanidine led to the isolation of a mutant strain which, when grown with fructose, oxidized both p-xylene and p-toluate to (-)-cis-1,2-dihydroxy-4-methylcyclohexa-3,5-diene-1-carboxylic acid (cis-p-toluate diol). The structure of the diol was determined by conventional chemical techniques including identification of the products formed by acid-catalyzed dehydration and characterization of a methyl ester derivative. The cis-relative stereochemistry of the hydroxyl groups was determined by the isolation and characterization of an isopropylidene derivative. p-Xylene-grown cells contained an inducible NAD+-dependent dehydrogenase which formed catechols from cis-p-toluate diol and the analogous acid diols formed from the other hydrocarbon substrates listed above. The catechols were converted to meta ring fission products by an inducible catechol-2,3-dioxygenase which was partially purified from p-xylene-grown cells of P. putida BG1.

Pseudomonas streptomycin resistance transposon associated with R-plasmid mobilization

Plasmid pMG1 encodes resistance to gentamicin, streptomycin, sulfonamides, and mercuric ions and also mobilizes pRO161, a transfer-deficient plasmid derived from RP1. Upon mobilization, pRO161 acquires streptomycin resistance (Smr) and can subsequently be remobilized by pMG1 at significantly higher frequencies than pRO161 itself. Both the initial acquisition of Smr and the subsequent mobilization of the transfer-deficient plasmid are recA independent: thus, the Smr determinant appears to be located on a transposon, disignated Tn904. Tn904 transposes to a variety of other plasmids, including RP1, FP2, R388, K, pRO1600, and pBR322, and in some cases the acquisition of this transposon accompanied deletions in the target plasmid. When no deletion occurred, target plasmids gained 5.2 kilobase pairs of DNA and new restriction endonuclease cleavage sites for AvaI, BglII, PstI, SmaI, and SstI. Physical analysis of such plasmids showed that the Tn904 termini are inverted repeat DNA sequences of approximately 124 base pairs. After cloning into vector pRO1723, a single site for restriction endonuclease AvaI was identified within the Smr determinant of Tn904. In Escherichia coli, but not in Pseudomonas aeruginosa. Tn904 shows a gene dosage-dependent expression of streptomycin resistance.