Large-scale concatenation cDNA sequencing

Mechanistic study of METTL3 inducing ferroptosis to promote cervical cancer progression through mediating m6A modification of COTE-1

Cervical Cancer (CC) is one of the leading causes of tumor-related deaths among women worldwide, and the mechanisms underlying the anti-ferroptosis of CC cells are still unclear. Methyltransferase like 3 (METTL3) is widely expressed various types of tissues and plays a crucial role in tumorigenesis in part by mediating cell death. However, its regulatory function in CC progression and especially the underlying mechanisms have not been fully elucidated. This study aims to explore the role of METTL3 in the ferroptosis of CC cells. Mechanistically, by MeRIP-seq, we identified COTE-1 as a target of METTL3 mediated m6A modification, and revealed that METTL3-mediated COTE-1 expression was dependent on the m6A reader-dependent manner. Functionally, in vitro and in vivo experiments that METTL3 promotes proliferation and metastasis of CC cells by regulating COTE-1 expression. In addition, the study verified the effect of the METTL3/COTE-1 axis on autophagy-dependent ferroptosis. In summary, METTL3 influences CC progression by mediating COTE-1 to influence autophagy-dependent ferroptosis, representing a potential therapeutic approach for treating CC.

Mitochondrial fusion and cristae reorganization facilitate acquisition of cardiomyocyte identity during reprogramming of murine fibroblasts

Cardiomyocytes (CMs) rely on mitochondrial energy produced in highly interconnected mitochondrial networks. Direct reprogramming of cardiac fibroblasts (CFs) into induced CMs (iCMs) shows promise for treating cardiac injury, but little work has investigated mitochondrial energetics and morphology during the conversion of CFs to iCMs. We characterized mitochondria during direct cardiac reprogramming of murine neonatal CFs (mnCFs). Reprogramming increased mitochondrial respiration and interconnectivity but not to the levels of native CMs. We therefore investigated whether perturbations to mitochondrial dynamics impacted reprogramming. Mitochondrial fusion (joining) was essential for iCM generation, while various fission (dividing) genes were reprogramming barriers. In particular, the loss of mitochondrial fission regulator 1 like (Mtfr1l) significantly increased the yield of functionally mature iCMs and induced mitochondrial fusion and respiration. These changes were countered by the concomitant loss of fusion effector optical atrophy protein 1 (Opa1). The present study advances our understanding of mitochondrial barriers to and mechanisms of direct cardiac reprogramming.

Highly expressed FAM189B predicts poor prognosis in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most malignant tumors with persistently high morbidity and mortality. However, the expression, prognostic and clinical significance of FAM189 family genes in HCC remain largely unknown. In this study, the expression levels of FAM189 family genes in HCC were analyzed through TCGA-LIHC and ICGC-LIRI-JP cohorts, and further validated in multiple independent GEO datasets. It was found that the expression of FAM189B was significantly upregulated in HCC tumor tissues, while the expression of FAM189A1 and FAM189A2 was not significantly changed between tumor and adjacent tissues. Further analysis revealed that upregulated copy number variation contributed to increased expression of FAM189B in HCC. Survival analysis showed that highly expressed FAM189B was significantly correlated with unfavorable prognosis, including overall survival, disease-specific survival, and progression-free interval. Univariate and multivariate Cox regression analysis showed that FAM189B was a potential novel prognosis factor for HCC patients. In addition, the association between FAM189B expression and clinical and molecular characteristics was analyzed. High expression of FAM189B was associated with high AFP level, high predicted risk metastasis signature, and TP53 mutation, while there was no significant association between FAM189B expression and cancer stage or tumor grade of HCC. Gene set enrichment analysis revealed that highly expressed FAM189B was closely related with signal pathways and biological processes associated with cell proliferation and cell cycle in HCC. In conclusion, this study suggested that FAM189B was highly expressed in HCC and highly expressed FAM189B may serve as an effective prognostic indicator and a potential therapeutic target for HCC patients.

Isoforms of Base Excision Repair Enzymes Produced by Alternative Splicing

Transcripts of many enzymes involved in base excision repair (BER) undergo extensive alternative splicing, but functions of the corresponding alternative splice variants remain largely unexplored. In this review, we cover the studies describing the common alternatively spliced isoforms and disease-associated variants of DNA glycosylases, AP-endonuclease 1, and DNA polymerase beta. We also discuss the roles of alternative splicing in the regulation of their expression, catalytic activities, and intracellular transport.

Transcriptional response in normal mouse tissues after i.v. (211)At administration - response related to absorbed dose, dose rate, and time

Background

In cancer radiotherapy, knowledge of normal tissue responses and toxicity risks is essential in order to deliver the highest possible absorbed dose to the tumor while maintaining normal tissue exposure at non-critical levels. However, few studies have investigated normal tissue responses in vivo after ²¹¹At administration. In order to identify molecular biomarkers of ionizing radiation exposure, we investigated genome-wide transcriptional responses to (very) low mean absorbed doses from ²¹¹At in normal mouse tissues.

Methods

Female BALB/c nude mice were intravenously injected with 1.7 kBq ²¹¹At and killed after 1 h, 6 h, or 7 days or injected with 105 or 7.5 kBq and killed after 1 and 6 h, respectively. Controls were mock-treated. Total RNA was extracted from tissue samples of kidney cortex and medulla, liver, lungs, and spleen and subjected to microarray analysis. Enriched biological processes were categorized after cellular function based on Gene Ontology terms.

Results

Responses were tissue-specific with regard to the number of significantly regulated transcripts and associated cellular function. Dose rate effects on transcript regulation were observed with both direct and inverse trends. In several tissues, Angptl4, Per1 and Per2, and Tsc22d3 showed consistent transcript regulation at all exposure conditions.

Conclusions

This study demonstrated tissue-specific transcriptional responses and distinct dose rate effects after ²¹¹At administration. Transcript regulation of individual genes, as well as cellular responses inferred from enriched transcript data, may serve as biomarkers in vivo. These findings expand the knowledge base on normal tissue responses and may help to evaluate and limit side effects of radionuclide therapy.

Electronic supplementary material

The online version of this article (doi:10.1186/s13550-014-0078-7) contains supplementary material, which is available to authorized users.

Upregulation of the putative oncogene COTE1 contributes to human hepatocarcinogenesis through modulation of WWOX signaling

Family with sequence similarity 189, also known as COTE1, has been found to be significantly upregulated in hepatocellular carcinoma (HCC) specimens and cell lines and is associated with tumor size and differentiation. Furthermore, COTE1 contributes to hepatocellular carcinogenesis. The overexpression of COTE1 enhanced in vitro cell viability and colony formation in soft agar, and in vivo tumorigenicity of HCC-derived Focus and Huh7 cells. In contrast, COTE1 knockdown via RNAi markedly suppressed these phenotypes in YY-8103 and WRL-68 HCC cell lines. Mechanistic analyses indicated that COTE1 physically associated with WW domain-containing oxidoreductase (WWOX) and modulated WWOX tyrosine phosphorylation. The ectopic overexpression of COTE1 inhibited the WWOX-p53 signaling pathway by reducing the phosphorylation of WWOX at the Tyr33 residue in Focus cells. Conversely, COTE1 silencing activated tyrosine 33 phosphorylation of WWOX and induced WWOX-p53 mediated mitochondrial apoptosis in WRL-68 cells. In addition, COTE1 upregulation in Huh7 cells blocked the WWOX-cyclin D1 pathway via dephosphorylation of WWOX Tyr287, stimulating cell cycle progression whereas phosphorylation of Tyr287 of WWOX induced by COTE1 silencing resulted in activation of WWOX-cyclin D1 signaling, leading to cell cycle arrest in YY-8103 cells. Together, our findings suggest that the cytoplasmic protein COTE1 contributes to HCC tumorigenesis by regulating cell proliferation through the modulation of WWOX signaling.

Genomics and proteomics in solving brain complexity

The human brain is extraordinarily complex, composed of billions of neurons and trillions of synaptic connections. Neurons are organized into circuit assemblies that are modulated by specific interneurons and non-neuronal cells, such as glia and astrocytes. Data on human genome sequences predicts that each of these cells in the human brain has the potential of expressing ∼20 000 protein coding genes and tens of thousands of noncoding RNAs. A major challenge in neuroscience is to determine (1) how individual neurons and circuitry utilize this potential during development and maturation of the nervous system, and for higher brain functions such as cognition, and (2) how this potential is altered in neurological and psychiatric disorders. In this review, we will discuss how recent advances in next generation sequencing, proteomics and bioinformatics have transformed our understanding of gene expression and the functions of neural circuitry, memory storage, and disorders of cognition.

Identification of transmembrane protein 134 as a novel LMP1-binding protein by using bimolecular fluorescence complementation and an enhanced retroviral mutagen

Latent membrane protein 1 (LMP1) of Epstein-Barr virus induces constitutive signaling in infected cells. LMP1 signaling requires oligomerization of LMP1 via its transmembrane domain, localization to lipid rafts in the membrane, and association of the LMP1 cytoplasmic domain to adaptor proteins, such as the tumor necrosis factor receptor-associated factors (TRAFs). Protein complementation is a novel technique to examine protein-protein interaction through the assembly of functional fluorescent proteins or enzymes from inactive fragments. A previous study in our lab demonstrated the use of bimolecular fluorescence complementation (BiFC) to study the assembly of the LMP1 signaling complexes within the plasma membrane of mammalian cells. In the present study, LMP1 was used as bait in a genome-wide BiFC screen with an enhanced retroviral mutagen to identify new LMP1-binding proteins. Our screen identified a novel LMP1-binding protein, transmembrane protein 134 (Tmem134). Tmem134 is a candidate oncogene that is amplified in breast cancer cell lines. Binding, colocalization, and cofractionation between LMP1 and Tmem134 were confirmed. Finally, Tmem134 affected LMP1-induced NF-κB induction. Together, these data suggest that BiFC is a unique and novel platform to identify proteins recruited to the LMP1-signaling complex.

Ectopic overexpression of COTE1 promotes cellular invasion of hepatocellular carcinoma

Family with sequence similarity 189, member B (FAM189B), alias COTE1, a putative oncogene selected by microarray, for the first time was here found to be significantly up-regulated in hepatocellular carcinoma (HCC) specimens and HCC cell lines. mRNA expression of COTE1 in HCC samples and cell lines was detected by reverse transcription-polymerase chain reaction (RT-PCR) and real-time PCR, while protein expression of COTE1 in HCC tissues was assessed by immunohistochemistry. In addition, invasion of HCC cells was observed after overexpressing or silencing COTE1. In the total of 48 paired HCC specimens, compared with the adjacent non-cancer tissues, the expression of COTE1 was up-regulated in 31 (p<0.01). In HCC cell lines, COTE1 expression was significantly higher than in normal human adult liver (p<0.01). Overexpression of COTE1 enhanced HCC-derived LM6 and MHCC-L cellular invasion in vitro. In contrast, COTE1 knockdown via RNAi markedly suppressed these phenotypes, as documented in LM3 and MHCC-H HCC cells. Mechanistic analyses indicated that COTE1 could physically associate with WW domain oxidoreductase (WWOX), a tumor suppressor. COTE1 may be closely correlated with invasion of hepatocellular carcinoma (HCC) cells and thus may serve as an effective target for gene therapy.

Cost-effective sequencing of full-length cDNA clones powered by a de novo-reference hybrid assembly

BACKGROUND

Sequencing full-length cDNA clones is important to determine gene structures including alternative splice forms, and provides valuable resources for experimental analyses to reveal the biological functions of coded proteins. However, previous approaches for sequencing cDNA clones were expensive or time-consuming, and therefore, a fast and efficient sequencing approach was demanded.

METHODOLOGY

We developed a program, MuSICA 2, that assembles millions of short (36-nucleotide) reads collected from a single flow cell lane of Illumina Genome Analyzer to shotgun-sequence approximately 800 human full-length cDNA clones. MuSICA 2 performs a hybrid assembly in which an external de novo assembler is run first and the result is then improved by reference alignment of shotgun reads. We compared the MuSICA 2 assembly with 200 pooled full-length cDNA clones finished independently by the conventional primer-walking using Sanger sequencers. The exon-intron structure of the coding sequence was correct for more than 95% of the clones with coding sequence annotation when we excluded cDNA clones insufficiently represented in the shotgun library due to PCR failure (42 out of 200 clones excluded), and the nucleotide-level accuracy of coding sequences of those correct clones was over 99.99%. We also applied MuSICA 2 to full-length cDNA clones from Toxoplasma gondii, to confirm that its ability was competent even for non-human species.

CONCLUSIONS

The entire sequencing and shotgun assembly takes less than 1 week and the consumables cost only approximately US$3 per clone, demonstrating a significant advantage over previous approaches.

A 24-residue presequence localizes human factor B to mitochondria

We reported previously that the human factor B precursor is a 215-amino acid polypeptide, the first 40 amino acid residues of which function as a mitochondrial targeting presequence [G.I. Belogrudov, Y. Hatefi, J. Biol. Chem. 277 (2002) 6097-6103]. Confocal microscopy of live HEK293 cells, transiently transfected with factor B constructs tagged at the C-terminus with green fluorescent protein (GFP) revealed that either a 40- or 25-residue presequence localized factor B to mitochondria. Indirect immunofluorescent labeling of fixed, permeabilized HEK293 cells that were transiently transfected with a construct lacking a presequence, showed diffuse, intracellular staining that was consistent with targeting of ectopically expressed factor B to cellular compartments distinct from the mitochondria. Mutants in which either Met(-25) or both Met(-25)/Met(-24) residues of the presequence were deleted exhibited decreased or undetectable levels, respectively, of the GFP-tagged factor B. The factor B presequence alone was shown to target a reporter polypeptide GFP to mitochondria. Our studies, therefore, demonstrate that a 24-residue presequence is sufficient to localize factor B to mitochondria, and suggest that the human factor B precursor is a 199-amino acid polypeptide.

Prolactin-induced changes in protein expression in human pancreatic islets

Ex vivo islet cell culture prior to transplantation appears as an attractive alternative for treatment of type 1 diabetes. Previous results from our laboratory have demonstrated beneficial effects of human prolactin (rhPRL) treatment on human islet primary cultures. In order to probe into the molecular events involved in the intracellular action of rhPRL in these cells, we set out to identify proteins with altered expression levels upon rhPRL cell treatment, using two-dimensional (2D) gel electrophoresis and mass spectrometry (MS). An average of 300 different protein spots were detected, 14 of which were modified upon rhPRL treatment (p<0.01), of which 12 were successfully identified using MS and grouped according to their biological functions. In conclusion, our study provides, for the first time, information about proteins that could be critically involved in PRL's action on human pancreatic islets, and facilitate identification of new and specific targets involved in islet cell function and proliferation.

Gene Expression Profiling of Ovarian Tissues for Determination of Molecular Pathways Reflective of Tumorigenesis

Ovarian cancer is the fourth leading cause of gynecological cancer death among women in the United States. Early detection is a critical prerequisite to initiating effective cancer therapy. Gene microarray technology and proteomics have provided much of the biomarkers with potential use for diagnosis. However, more research is needed to fully understand disease onset and progression. To this end, we have performed microarray analysis with the goal of identifying molecular interaction networks defining tumor growth. Microarray analysis was performed on a limited set of ovarian tissues with various pathological diagnoses using Human Genome Focus Array (HGFA) for the detection of approximately 8500 human transcripts. Hierarchical clustering identified groups of ovarian tissues reflective of low malignant potential/early cancer onset and possible pre-cancerous stages involving small molecule, cytokine and/or hormone-dependent feed-back responses specific to the pelvic reproductive system and a priori initiated tumor suppression mechanisms. ANOVA followed by post hoc Scheffe confirmed our hypotheses. Moreover, we established a protein/protein interaction database associated with HGFA probe sets. This database was used to build and visualize molecular networks integrating small but significant changes in gene expression. In conclusion, we were able for the first time to delineate an intersecting genetic pattern linking ovarian tissues reflective of low potential malignancy/early cancer onset stages via long distance signaling between tissues of gynecological origin.

Pooled Genomic Indexing (PGI): analysis and design of experiments

Pooled Genomic Indexing (PGI) is a novel method for physical mapping of clones onto known sequences. PGI is carried out by pooling arrayed clones and generating shotgun sequence reads from the pools. The shotgun sequences are compared to a reference sequence. In the simplest case, clones are placed on an array and are pooled by rows and columns. If a shotgun sequence from a row pool and another shotgun sequence from a column pool match the reference sequence at a close distance, they are both assigned to the clone at the intersection of the two pools. Accordingly, the clone is mapped onto the region of the reference sequence between the two matches. A probabilistic model for PGI is developed, and several pooling designs are described and analyzed, including transversal designs and designs from linear codes. The probabilistic model and the pooling schemes are validated in simulated experiments where 625 rat bacterial artificial chromosome (BAC) clones and 207 mouse BAC clones are mapped onto homologous human sequence.

HYPERSENSITIVE TO RED AND BLUE 1, a ZZ-type zinc finger protein, regulates phytochrome B-mediated red and cryptochrome-mediated blue light responses

Plant photoreceptors that regulate photomorphogenic development include red/far-red-light-absorbing phytochromes and blue/UV-A-light-absorbing cryptochromes. We have undertaken a genetic screen to identify additional components downstream of the photoreceptors in Arabidopsis thaliana. We identified a short hypocotyl mutant under red and blue light, hypersensitive to red and blue 1 (hrb1). Mutation in HRB1 also enhances the end-of-day far-red light response, inhibits leaf expansion and petiole elongation, and attenuates the expression of CAB3 and CHS. Double mutant analysis indicates that phyB is epistatic to hrb1 under red light, and cry1 cry2 is epistatic to hrb1 under blue light for both hypocotyl growth and light-regulated gene expression responses. HRB1 localizes to the nucleus and belongs to a protein family of Drought induced 19 (Di19). HRB1 and all other family members contain a ZZ-type zinc finger domain, which in other organisms is implicated in protein-protein interactions between dystrophin and calmodulin and between transcriptional adaptors and activators. HRB1 activity is also required for red and blue light-induced expression of PHYTOCHROME INTERACTING FACTOR 4 (PIF4). pif4 shows a very similar hypersensitive response as hrb1 to both red light and blue light and is epistatic to hrb1 in control of light-regulated gene expression responses. Thus, the roles of HRB1 and PIF4 together in regulating both red and blue light responses may represent points where red light signaling and blue light signaling intersect.

Identification of a novel partner of duox: EFP1, a thioredoxin-related protein

H(2)O(2) is a crucial substrate of thyroproxidase (TPO) to iodinate thyroglobulin and synthesize thyroid hormones in thyroid. ThOX proteins (thyroid oxidase) also called Duox are believed to be responsible for H(2)O(2) generation. Duoxs expressed in transfected cells do not generate an active system, nor permit their membrane localization suggesting that other proteins are required to fulfill these functions. In this study, we demonstrate interactions of Duoxs with TPO and with p22(phox) without any effect on Duox activity. By yeast two-hybrid method using EF-hand fragment of dog Duox1 as the bait we have isolated EFP1 (EF-hand binding protein 1), one partner of Duoxs that belongs to the thioredoxin-related protein family. EFP1 shares moderate similarities with other members of thioredoxin-related proteins, but the characteristic active site and the folding structures are well conserved. EFP1 can be co-immunoprecipitated with Duoxs in transfected COS cells as well as in primary cultured human thyrocytes. It interacts also with TPO but not thyroglobulin. Immunofluorescence studies show that EFP1 and Duox proteins are co-localized inside the transfected cells, suggesting that EFP1 is not sufficient to induce either the expression of Duox at the plasma membrane or to permit H(2)O(2) production. EFP1 and Duox mRNA share similar distribution in nine different tissues. These results suggest that EFP1 could be one of the partners in the assembly of the multiprotein complex constituting the thyroid H(2)O(2) generating system but is certainly not sufficient to permit H(2)O(2) generation.

Analysis of human mRNAs with the reference genome sequence reveals potential errors, polymorphisms, and RNA editing

The NCBI Reference Sequence (RefSeq) project and the NIH Mammalian Gene Collection (MGC) together define a set of approximately 30,000 nonredundant human mRNA sequences with identified coding regions representing 17,000 distinct loci. These high-quality mRNA sequences allow for the identification of transcribed regions in the human genome sequence, and many researchers accept them as the correct representation of each defined gene sequence. Computational comparison of these mRNA sequences and the recently published essentially finished human genome sequence reveals several thousand undocumented nonsynonymous substitution and frame shift discrepancies between the two resources. Additional analysis is undertaken to verify that the euchromatic human genome is sufficiently complete--containing nearly the whole mRNA collection, thus allowing for a comprehensive analysis to be undertaken. Many of the discrepancies will prove to be genuine polymorphisms in the human population, somatic cell genomic variants, or examples of RNA editing. It is observed that the genome sequence variant has significant additional support from other mRNAs and ESTs, almost four times more often than does the mRNA variant, suggesting that the genome sequence is more accurate. In approximately 15% of these cases, there is substantial support for both variants, suggestive of an undocumented polymorphism. An initial screening against a 24-individual genomic DNA diversity panel verified 60% of a small set of potential single nucleotide polymorphisms from which successful results could be obtained. We also find statistical evidence that a few of these discrepancies are due to RNA editing. Overall, these results suggest that the mRNA collections may contain a substantial number of errors. For current and future mRNA collections, it may be prudent to fully reconcile each genome sequence discrepancy, classifying each as a polymorphism, site of RNA editing or somatic cell variation, or genome sequence error.

Large-scale RT-PCR recovery of full-length cDNA clones

Pseudogenes, alternative transcripts, noncoding RNA, and polymorphisms each add extensive complexity to the mammalian transcriptome and confound estimation of the total number of genes. Despite advanced algorithms for gene prediction and several large-scale efforts to obtain cDNA clones for all human open reading frames (ORFs), no single collection is complete. To enhance this effort, we have developed a high-throughput pipeline for reverse transcription PCR (RT-PCR) gene recovery. Most importantly, novel molecular strategies for improving RT-PCR yield of transcripts that have been difficult to isolate by other means and computational strategies for clone sequence validation have been developed and optimized. This systematic gene recovery pipeline allows both rescue of predicted human and rat genes and provides insight into the complexity of the transcriptome through comparisons with existing data sets.

hIscA: a protein implicated in the biogenesis of iron–sulfur clusters

This article describes the gene called hIscA, its transcription product and protein (hIscA) and its putative function. We screened a human brain cDNA expression library with serum from a patient suffering from the autoimmune Sjögren's syndrome (S5:823/94). One cDNA of 1.6-kbp clone was isolated. This clone contains the entire coding sequence for a protein unknown in human. IscA is ubiquitously expressed and expression levels vary among tissues. The 15.5-kDa predicted protein contains a structural domain named HESB, is located in the mitochondria and is implicated in the biogenesis of iron-sulfur clusters. Since this unknown protein is related to IscA-like protein, we suggest as name for this protein hIscA. The recombinant protein is recognized by a rabbit polyclonal antiserum generated against the carboxyl extreme of the Saccharomyces cerevisiae orthologue Isa1. In this article, we demonstrate the functional homology between hIscA and Isa1 proteins using Isa1 null mutant S. cerevisiae transformed with hIscA in a yeast functional complementation test. We also describe the rat homologue to this gene.

Concatenation cDNA sequencing for transcriptome analysis

We describe a high-throughput cDNA sequencing pipeline (http://www.hgsc.bcm.tmc.edu/projects/cdna) built in response to the emerging need for rapid sequencing of large cDNA collections. Using this strategy cDNA inserts are purified and joined through concatenation into large molecules. These 'pseudo-BACs' are subjected to random shotgun sequencing whereby the majority of cDNA inserts in the pool are sequenced. Using this concatenation cDNA sequencing platform, we have contributed more than 13000 full-length cDNA sequences from human and mouse to the Mammalian Gene Collection (MGC).

Novel potential ALL low-risk markers revealed by gene expression profiling with new high-throughput SSH-CCS-PCR

The current systems of risk grouping in pediatric acute lymphoblastic leukemia (ALL) fail to predict therapeutic success in 10-35% of patients. To identify better predictive markers of clinical behavior in ALL, we have developed an integrated approach for gene expression profiling that couples suppression subtractive hybridization, concatenated cDNA sequencing, and reverse transcriptase real-time quantitative PCR. Using this approach, a total of 600 differentially expressed genes were identified between t(4;11) ALL and pre-B ALL with no determinant chromosomal translocation. The expression of 67 genes was analyzed in different cytogenetic ALL subgroups and B lymphocytes isolated from healthy donors. Three genes, BACH1, TP53BPL, and H2B/S, were consistently expressed as a significant cluster associated with the low-risk ALL subgroups. A total of 42 genes were differentially expressed in ALL vs normal B lymphocytes, with no specific association with any particular ALL subgroups. The remaining 22 genes were part of a specific expression profile associated with the hyperdiploid, t(12;21), or t(4;11) subgroups. Using an unsupervised hierarchical cluster analysis, the discriminating power of these specific expression profiles allowed the clustering of patients according to their subgroups. These genes could help to understand the difference in treatment response and become therapeutical targets to improve ALL clinical outcomes.

Cloning and expression of a novel cDNA encoding shell coat protein, cp4, from the brushtail possum (Trichosurus vulpecula)

The marsupial conceptus is enclosed by several egg coats of evolutionary significance and unknown composition, of which the shell coat in mammals occurs only in marsupials and monotremes. Intact coats are vital to marsupial embryonic development. Towards a better understanding of the marsupial coat proteins, a cDNA sequence (cp4) encoding a shell coat protein was cloned from the brushtail possum. A cDNA library of a zygote stage uterus was screened using a deduced oligonucleotide sequence based on a partial amino acid sequence of the coat protein. This study has confirmed a single copy cp4 gene encoding a unique protein of 306 amino acids, although its N-terminus shares high sequence identity with the C-terminal half of the enzyme alpha-enolase. Using Northern blots, the expression of cp4 in adult tissues showed that cp4 transcript is restricted predominantly to the uterus with two stages of expression occurring in the gravid uterus at early cleavage and bilaminar stages, suggesting an important developmental role for CP4. Using RT-PCR, cp4-specific expression as represented by the 3'-end 400 bp was present in heart, liver, oviduct, and uterus. Uterine expression reflected the increase found with Northern blot except that expression was low at unilaminar and bilaminar stages.

The gene structure and expression of human ABHD1: overlapping polyadenylation signal sequence with Sec12

BACKGROUND

Overlapping sense/antisense genes orientated in a tail-to-tail manner, often involving only the 3'UTRs, form the majority of gene pairs in mammalian genomes and can lead to the formation of double-stranded RNA that triggers the destruction of homologous mRNAs. Overlapping polyadenylation signal sequences have not been described previously.

RESULTS

An instance of gene overlap has been found involving a shared single functional polyadenylation site. The genes involved are the human alpha/beta hydrolase domain containing gene 1 (ABHD1) and Sec12 genes. The nine exon human ABHD1 gene is located on chromosome 2p23.3 and encodes a 405-residue protein containing a catalytic triad analogous to that present in serine proteases. The Sec12 protein promotes efficient guanine nucleotide exchange on the Sar1 GTPase in the ER. Their sequences overlap for 42 bp in the 3'UTR in an antisense manner. Analysis by 3' RACE identified a single functional polyadenylation site, ATTAAA, within the 3'UTR of ABHD1 and a single polyadenylation signal, AATAAA, within the 3'UTR of Sec12. These polyadenylation signals overlap, sharing three bp. They are also conserved in mouse and rat. ABHD1 was expressed in all tissues and cells examined, but levels of ABHD1 varied greatly, being high in skeletal muscle and testis and low in spleen and fibroblasts.

CONCLUSIONS

Mammalian ABHD1 and Sec12 genes contain a conserved 42 bp overlap in their 3'UTR, and share a conserved TTTATTAAA/TTTAATAAA sequence that serves as a polyadenylation signal for both genes. No inverse correlation between the respective levels of ABHD1 and Sec12 RNA was found to indicate that any RNA interference occurred.

Generation and initial analysis of more than 15,000 full-length human and mouse cDNA sequences

The National Institutes of Health Mammalian Gene Collection (MGC) Program is a multiinstitutional effort to identify and sequence a cDNA clone containing a complete ORF for each human and mouse gene. ESTs were generated from libraries enriched for full-length cDNAs and analyzed to identify candidate full-ORF clones, which then were sequenced to high accuracy. The MGC has currently sequenced and verified the full ORF for a nonredundant set of >9,000 human and >6,000 mouse genes. Candidate full-ORF clones for an additional 7,800 human and 3,500 mouse genes also have been identified. All MGC sequences and clones are available without restriction through public databases and clone distribution networks (see http:mgc.nci.nih.gov).

Biochemical characterisation and genetic analysis of aureocin A53, a new, atypical bacteriocin from Staphylococcus aureus

Aureocin A53 is produced by Staphylococcus aureus A53. It is encoded on a 10.4 kb plasmid, pRJ9, and is active against Listeria monocytogenes. Aureocin A53 is a highly cationic 51-residue peptide containing ten lysine and five tryptophan residues. Aureocin A53 was purified to homogeneity by hydrophobic-interaction, cation-exchange, and reverse-phase chromatography. MALDI-TOF mass spectrometry yielded a molecular mass of 6012.5 Da, which was 28 Da higher than predicted from the structural gene sequence of the bacteriocin. The mass increment resulted from an N-formylmethionine residue, indicating that the aureocin A53 is synthesised and secreted without a typical bacteriocin leader sequence or sec-dependent signal peptide. The structural identity of aureocin A53 was verified by Edman sequencing after de-blocking with cyanogen bromide and extensive mass spectrometry analysis of enzymatically and laser-generated fragments. The complete sequence of pRJ9 was determined and none of the open reading frames identified in the vicinity of the structural gene aucA showed similarity to genes that are typically found in bacteriocin gene clusters. Thus, neither a dedicated protease or transporter, nor modifying enzymes and regulatory elements seemed to be involved in the production of aureocin A53. Further unique features that distinguish aureocin A53 from other peptide bacteriocins include remarkable protease stability and a defined, rigid structure in aqueous solution.

Systematic sequencing of cDNA clones using the transposon Tn5

In parallel with the production of genomic sequence data, attention is being focused on the generation of comprehensive cDNA-sequence resources. Such efforts are increasingly emphasizing the production of high-accuracy sequence corresponding to the entire insert of cDNA clones, especially those presumed to reflect the full-length mRNA. The complete sequencing of cDNA clones on a large scale presents unique challenges because of the generally small, yet heterogeneous, sizes of the cloned inserts. We have developed a strategy for high-throughput sequencing of cDNA clones using the transposon Tn5. This approach has been tailored for implementation within an existing large-scale 'shotgun-style' sequencing program, although it could be readily adapted for use in virtually any sequencing environment. In addition, we have developed a modified version of our strategy that can be applied to cDNA clones with large cloning vectors, thereby overcoming a potential limitation of transposon-based approaches. Here we describe the details of our cDNA-sequencing pipeline, including a summary of the experience in sequencing more than 4200 cDNA clones to produce more than 8 million base pairs of high-accuracy cDNA sequence. These data provide both convincing evidence that the insertion of Tn5 into cDNA clones is sufficiently random for its effective use in large-scale cDNA sequencing as well as interesting insight about the sequence context preferred for insertion by Tn5.

A Drosophila full-length cDNA resource

BACKGROUND

A collection of sequenced full-length cDNAs is an important resource both for functional genomics studies and for the determination of the intron-exon structure of genes. Providing this resource to the Drosophila melanogaster research community has been a long-term goal of the Berkeley Drosophila Genome Project. We have previously described the Drosophila Gene Collection (DGC), a set of putative full-length cDNAs that was produced by generating and analyzing over 250,000 expressed sequence tags (ESTs) derived from a variety of tissues and developmental stages.

RESULTS

We have generated high-quality full-insert sequence for 8,921 clones in the DGC. We compared the sequence of these clones to the annotated Release 3 genomic sequence, and identified more than 5,300 cDNAs that contain a complete and accurate protein-coding sequence. This corresponds to at least one splice form for 40% of the predicted D. melanogaster genes. We also identified potential new cases of RNA editing.

CONCLUSIONS

We show that comparison of cDNA sequences to a high-quality annotated genomic sequence is an effective approach to identifying and eliminating defective clones from a cDNA collection and ensure its utility for experimentation. Clones were eliminated either because they carry single nucleotide discrepancies, which most probably result from reverse transcriptase errors, or because they are truncated and contain only part of the protein-coding sequence.

A clone-array pooled shotgun strategy for sequencing large genomes

A simplified strategy for sequencing large genomes is proposed. Clone-Array Pooled Shotgun Sequencing (CAPSS) is based on pooling rows and columns of arrayed genomic clones, for shotgun library construction. Random sequences are accumulated, and the data are processed by sequential comparison of rows and columns to assemble the sequence of clones at points of intersection. Compared with either a clone-by-clone approach or whole-genome shotgun sequencing, CAPSS requires relatively few library constructions and only minimal computational power for a complete genome assembly. The strategy is suitable for sequencing large genomes for which there are no sequence-ready maps, but for which relatively high resolution STS maps and highly redundant BAC libraries are available. It is immediately applicable to the sequencing of mouse, rat, zebrafish, and other important genomes, and can be managed in a cooperative fashion to take advantage of a distributed international DNA sequencing capacity.

Molecular characterisation of aureocin A70, a multi-peptide bacteriocin isolated from Staphylococcus aureus

Staphylococcus aureus A70 produces a heat-stable bacteriocin designated aureocin A70. Aureocin A70 is encoded within a mobilisable 8 kb plasmid, pRJ6, and is active against Listeria monocytogenes. Experiments of transposition mutagenesis and gene cloning had shown that aureocin A70 production and immunity were associated with the HindIII-A and B fragments of pRJ6. Therefore, a 6332 bp region of the plasmid, encompassing both these fragments, was sequenced using a concatenation DNA sequencing procedure. DNA sequence and genetic analyses revealed the presence of three transcriptional units that appear to be involved in bacteriocin activity. The first transcriptional unit contains a single gene, aurT, which encodes a protein that resembles an ATP-dependent transporter, similar to those involved in lantibiotic export. AurT is required for aureocin A70 production and it appears to be essential for mobilisation of pRJ6. The second putative operon contains two open reading frames (ORFs); the first gene, orfA, is predicted to encode a protein similar to small repressor proteins found in some Archaea, whose function remains to be elucidated. The second gene, orfB, codes for an 138 amino acid residue protein which shares a number of characteristics (high pI and hydrophobicity profile) with proteins associated with immunity, needed for self-protection against bacteriocin. Four other genes are present in the third operon, aurABCD. aurABCD encode four related peptides that are small (30-31 amino acid residues), strongly cationic (pI of 9.85 to 10.04) and highly hydrophobic. Theses peptides also have a high content of small amino acid residues like glycine and alanine, and no cysteine residue. Tn917-lac insertional mutations, which affected aureocin A70 activity, reside within operon aurABCD. Analysis of purified bacteriocin preparations by mass spectrometry demonstrated that all four peptides encoded by aurABCD operon are produced, expressed and excreted without post-translational modifications. Thus, aureocin A70 is a multi-peptide non-lantibiotic bacteriocin, which is transported without processing.

Identification of the gene altered in Berardinelli-Seip congenital lipodystrophy on chromosome 11q13

Congenital generalized lipodystrophy, or Berardinelli-Seip syndrome (BSCL), is a rare autosomal recessive disease characterized by a near-absence of adipose tissue from birth or early infancy and severe insulin resistance. Other clinical and biological features include acanthosis nigricans, hyperandrogenism, muscular hypertrophy, hepatomegaly, altered glucose tolerance or diabetes mellitus, and hypertriglyceridemia. A locus (BSCL1) has been mapped to 9q34 with evidence of heterogeneity. Here, we report a genome screen of nine BSCL families from two geographical clusters (in Lebanon and Norway). We identified a new disease locus, designated BSCL2, within the 2.5-Mb interval flanked by markers D11S4076 and D11S480 on chromosome 11q13. Analysis of 20 additional families of various ethnic origins led to the identification of 11 families in which the disease cosegregates with the 11q13 locus; the remaining families provide confirmation of linkage to 9q34. Sequence analysis of genes located in the 11q13 interval disclosed mutations in a gene homologous to the murine guanine nucleotide-binding protein (G protein), gamma3-linked gene (Gng3lg) in all BSCL2-linked families. BSCL2 is most highly expressed in brain and testis and encodes a protein (which we have called seipin) of unknown function. Most of the variants are null mutations and probably result in a severe disruption of the protein. These findings are of general importance for understanding the molecular mechanisms underlying regulation of body fat distribution and insulin resistance.

Chromosomal mapping, gene structure and characterization of the human and murine RAB27B gene

BACKGROUND

Rab GTPases are regulators of intracellular membrane traffic. The Rab27 subfamily consists of Rab27a and Rab27b. Rab27a has been recently implicated in Griscelli Disease, a disease combining partial albinism with severe immunodeficiency. Rab27a plays a key role in the function of lysosomal-like organelles such as melanosomes in melanocytes and lytic granules in cytotoxic T lymphocytes. Little is known about Rab27b.

RESULTS

The human RAB27B gene is organised in six exons, spanning about 69 kb in the chromosome 18q21.1 region. Exon 1 is non-coding and is separated from the others by 49 kb of DNA and exon 6 contains a long 3' untranslated sequence (6.4 kb). The mouse Rab27b cDNA shows 95% identity with the human cDNA at the protein level and maps to mouse chromosome 18. The mouse mRNA was detected in stomach, large intestine, spleen and eye by RT-PCR, and in heart, brain, spleen and kidney by Northern blot. Transient over-expression of EGF-Rab27b fusion protein in cultured melanocytes revealed that Rab27b is associated with melanosomes, as observed for EGF-Rab27a.

CONCLUSIONS

Our results indicate that the Rab27 subfamily of Ras-like GTPases is highly conserved in mammals. There is high degree of conservation in sequence and gene structure between RAB27A and RAB27B genes. Exogenous expression of Rab27b in melanocytes results in melanosomal association as observed for Rab27a, suggesting the two Rab27 proteins are functional homologues. As with RAB27A in Griscelli Disease, RAB27B may be also associated with human disease mapping to chromosome 18.

Survey of transcripts in the adult Drosophila brain

BACKGROUND

Classic methods of identifying genes involved in neural function include the laborious process of behavioral screening of mutagenized flies and then rescreening candidate lines for pleiotropic effects due to developmental defects. To accelerate the molecular analysis of brain function in Drosophila we constructed a cDNA library exclusively from adult brains. Our goal was to begin to develop a catalog of transcripts expressed in the brain. These transcripts are expected to contain a higher proportion of clones that are involved in neuronal function.

RESULTS

The library contains approximately 6.75 million independent clones. From our initial characterization of 271 randomly chosen clones, we expect that approximately 11% of the clones in this library will identify transcribed sequences not found in expressed sequence tag databases. Furthermore, 15% of these 271 clones are not among the 13,601 predicted Drosophila genes.

CONCLUSIONS

Our analysis of this unique Drosophila brain library suggests that the number of genes may be underestimated in this organism. This work complements the Drosophila genome project by providing information that facilitates more complete annotation of the genomic sequence. This library should be a useful resource that will help in determining how basic brain functions operate at the molecular level.

Drosophila Double parked: a conserved, essential replication protein that colocalizes with the origin recognition complex and links DNA replication with mitosis and the down-regulation of S phase transcripts

We identified a Drosophila gene, double parked(dup), that is essential for DNA replication and belongs to a new family of replication proteins conserved fromSchizosaccharomyces pombe to humans. Strong mutations indup cause embryonic lethality, preceded by a failure to undergo S phase during the postblastoderm divisions. dup is required also for DNA replication in the adult ovary, establishing thatdup is needed for DNA replication at multiple stages of development. Strikingly, DUP protein colocalizes with the origin recognition complex to specific sites in the ovarian follicle cells. This suggests that DUP plays a direct role in DNA replication. Thedup transcript is cell cycle regulated and is under the control of E2F and Cyclin E. Interestingly, dup mutant embryos fail both to downregulate S phase genes and to engage a checkpoint preventing mitosis until completion of S phase. This could be either because these events depend on progression of S phase beyond the point blocked in the dup mutants or because DUP is needed directly for these feedback mechanisms.

Drosophila double parked: a conserved, essential replication protein that colocalizes with the origin recognition complex and links DNA replication with mitosis and the down-regulation of S phase transcripts

We identified a Drosophila gene, double parked (dup), that is essential for DNA replication and belongs to a new family of replication proteins conserved from Schizosaccharomyces pombe to humans. Strong mutations in dup cause embryonic lethality, preceded by a failure to undergo S phase during the postblastoderm divisions. dup is required also for DNA replication in the adult ovary, establishing that dup is needed for DNA replication at multiple stages of development. Strikingly, DUP protein colocalizes with the origin recognition complex to specific sites in the ovarian follicle cells. This suggests that DUP plays a direct role in DNA replication. The dup transcript is cell cycle regulated and is under the control of E2F and Cyclin E. Interestingly, dup mutant embryos fail both to downregulate S phase genes and to engage a checkpoint preventing mitosis until completion of S phase. This could be either because these events depend on progression of S phase beyond the point blocked in the dup mutants or because DUP is needed directly for these feedback mechanisms.

Identification of a glialblastoma cell differentiation factor-related gene mRNA in human microvascular endothelial cells

Vascular endothelial cells (VEC) transduce mitogenic and chemoattractant signals in response to erythropoietin (Epo). An analysis of changes in gene expression in VEC would be helpful to understanding the molecular nature of mitogenic signals. An effective method for analysis of gene expression is through differential display. Using this approach, we obtained from Epo-treated human microvascular endothelial cells (HMVEC) a cDNA fragment with characteristics of the 3'end of mRNA. Using the cDNA fragment, we then isolated a full-length clone from a HMVEC cDNA library. The cDNA of interest encodes a protein consisting of 404 amino acids with a carboxy-terminal end sequence identical to glialblastoma cell differentiation factor-related protein (GBDR1). Northern blot analysis showed that GBDR1 mRNA was ubiquitously expressed in human tissues. In Southern blot analysis, GBDR1 cDNA identified a single gene on chromosome 9. Since analysis of the amino acid sequence revealed several putative phosphorylation sites for different protein kinases, the GBDR1 protein was expressed and purified from bacterial extracts and, as predicted, casein kinase II phosphorylated GBDR1 in vitro. Immunofluorescence and biochemical data revealed that the GBDR1 protein is not entirely localized in the cytosolic fraction, suggesting that it may interact with another protein(s). These findings demonstrate that GBDR1 is an intracellular signaling molecule that may play a role in the regulation of endothelial cell growth.

Octamer-primed sequencing technology: development of primer identification software

Octamer sequencing technology (OST) is a primer-directed sequencing strategy in which an individual octamer primer is selected from a pre-synthesized octamer primer library and used to sequence a DNA fragment. However, selecting candidate primers from such a library is time consuming and can be a bottleneck in the sequencing process. To accelerate the sequencing process and to obtain high quality sequencing data, a computer program, electronic OST or eOST, was developed to automatically identify candidate primers from an octamer primer library. eOST integrates the base calling software PHRED to provide a quality assessment for target sequences and identifies potential primer binding sites located within a high quality target region. To increase the sequencing success rate, eOST includes a simple dynamic folding algorithm to automatically calculate the free energy and predict the secondary structure within the template in the vicinity of the octamer-binding site. Several parameters were found to be important, including base quality threshold, the window size of the template sequence segment, and the threshold [Delta] G value. OST, coupled with the eOST software, can be used to sequence short DNA fragments or in the finishing assembly stage of large-scale sequencing of genomic DNA.

Molecular and immunochemical evidences demonstrate that endooligopeptidase A is the predominant cytosolic oligopeptidase of rabbit brain

Oligopeptidases are tissue endopeptidases that do not attack proteins and are likely to be involved in the maturation and degradation of peptide hormones and neuropeptides. The rabbit brain endooligopeptidase A and the rat testes soluble metallopeptidase (EC 3.4.24.15) are thiol-activated oligopeptidases which are able to generate enkephalin from a number of opioid peptides and to inactivate bradykinin and neurotensin by hydrolyzing the same peptide bonds. A monospecific antibody raised against the purified rabbit brain endooligopeptidase A allowed the identification of a 2. 3 kb cDNA coding for a truncated enzyme of 512 amino acids, displaying the same enzymatic features as endooligopeptidase A. In spite of all efforts, employing several strategies, the full-length cDNA could not be cloned until now. The analysis of the deduced amino acid sequence showed no similarity to the rat testes metalloendopeptidase sequence, except for the presence of the typical metalloprotease consensus sequence [HEXXH]. The antibody raised against recombinant endooligopeptidase A specifically inhibited its own activity and reduced the thiol-activated oligopeptidase activity of rabbit brain cytosol to less than 30%. Analysis of the endooligopeptidase A tissue distribution indicated that this enzyme is mainly expressed in the CNS, whereas the soluble metallo EC 3.4.24.15 is mainly expressed in peripheral tissues.

Expression patterns of two new members of the Semaphorin family in Drosophila suggest early functions during embryogenesis

We report the sequence and expression analysis of two new Drosophila members of the Semaphorin family. Both proteins show the presence of Semaphorin domains and transmembrane domains. Both genes are expressed maternally and in embryos, and reveal distinct expression patterns much earlier than the onset of neurogenesis. We also present an overview of the domain structure of all so far known semaphorins in Drosophila. Furthermore, we compared all Drosophila and C. elegans Semaphorins and discuss them in the light of their evolution.

Sequencing analysis of forty-eight human image cDNA clones similar to Drosophila mutant protein

We have sequenced 48 human IMAGE cDNA clones selected from the public EST database (dbEST) for their significant homology to Drosophila mutant genes. A dynamically updated analysis report was produced by BlastX and BlastN analysis searches in the latest databases available. This analysis led us to estimate the grade of similarity with homologous genes isolated in other species. Bottlenecks were detected in the sequencing process and here we have presented our problem-solving approach. We think the value of this full-length sequencing project is an enrichment of the sequence database information that is currently available to the human genome community.

Plexin A is a neuronal semaphorin receptor that controls axon guidance

The Semaphorins comprise a large family of secreted and transmembrane proteins, some of which function as repellents during axon guidance. Semaphorins fall into seven subclasses. Neuropilins are neuronal receptors for class III Semaphorins. In the immune system, VESPR, a member of the Plexin family, is a receptor for a viral-encoded Semaphorin. Here, we identify two Drosophila Plexins, both of which are expressed in the developing nervous system. We present evidence that Plexin A is a neuronal receptor for class I Semaphorins (Sema 1a and Sema 1b) and show that Plexin A controls motor and CNS axon guidance. Plexins, which themselves contain complete Semaphorin domains, may be both the ancestors of classical Semaphorins and binding partners for Semaphorins.