Suite of tools for statistical N-gram language modeling for pattern mining in whole genome sequences

Genome sequences contain a number of patterns that have biomedical significance. Repetitive sequences of various kinds are a primary component of most of the genomic sequence patterns. We extended the suffix-array based Biological Language Modeling Toolkit to compute n-gram frequencies as well as n-gram language-model based perplexity in windows over the whole genome sequence to find biologically relevant patterns. We present the suite of tools and their application for analysis on whole human genome sequence.

A comprehensive analysis of the chorion locus in silkmoth

Despite more than 40 years of intense study, essential features of the silkmoth chorion (eggshell) are still not fully understood. To determine the precise structure of the chorion locus, we performed extensive EST analysis, constructed a bacterial artificial chromosome (BAC) contig, and obtained a continuous genomic sequence of 871,711 base pairs. We annotated 127 chorion genes in two segments interrupted by a 164 kb region with 5 non-chorion genes, orthologs of which were on chorion bearing scaffolds in 4 ditrysian families. Detailed transcriptome analysis revealed expression throughout choriogenesis of most chorion genes originally categorized as ``middle'', and evidence for diverse regulatory mechanisms including cis-elements, alternative splicing and promoter utilization, and antisense RNA. Phylogenetic analysis revealed multigene family associations and faster evolution of early chorion genes and transcriptionally active pseudogenes. Proteomics analysis identified 99 chorion proteins in the eggshell and micropyle localization of 1 early and 6 Hc chorion proteins.

Expression regulation and functional characterization of a novel interferon inducible gene Gig2 and its promoter

Grass carp hemorrhagic virus (GCHV)-induced gene 2 (Gig2) is a novel gene previously identified from UV-inactivated GCHV-treated Carassius auratus blastulae embryonic (CAB) cells, suggesting that it should play a pivotal role in the interferon (IFN) antiviral response. In this study, a polyclonal anti-Gig2 antiserum was generated and used to study the inductive expression pattern by Western blot analysis, showing no basal expression in normal CAB cells but a significant up-regulation upon UV-inactivated GCHV, polyinosinic:polycytidylic acid (Poly I:Q and recombinant IFN (rIFN). However, constitutive expression of Gig2 is observed in all tested tissues from grass carp (Ctenopharyngodon idellus), and Poly I:C injection increases the relative amount of Gig2 protein in skin, spleen, trunk kidney, gill, hindgut and thymus. Moreover, the genomic sequence covering the whole Gig2 ORF and the upstream promoter region were amplified by genomic walking. Significantly, the Gig2 promoter contains three IFN-stimulated response elements (ISREs), nine GAAA/TfTC motifs and five gamma-IFN activating sites (GAS), which are the characteristics of genes responsive to both type I IFN and type 11 IFN. Subsequently, the complete Gig2 promoter sequence was cloned into pGL3-Basic vector, and its activity was measured by luciferase assays in the transfected CAB cells. The Gig2 promoter-driven construct is highly induced in CAB cells after treatment with Poly I:C or rIFN, and the functional capability is dependent on IFN regulatory factor 7 (IRF7), because its activity can be stimulated by IRF7. Collectively, the data provide strong evidence that Gig2 is indeed a novel IFN inducible gene and its expression is likely dependent on IRF7 upon Poly I:C or IFN. (C) 2009 Elsevier Ltd. All rights reserved.

Metallothionein-2 gene from the mandarin fish Siniperca chuatsi: cDNA cloning, tissue expression, and immunohistochemical localization

The metallothionein-2 (MT-2) gene was isolated from the mandarin fish, one of the most important industrial aquatic animals in China, by using rapid amplification of cDNA ends (RACE). The deduced amino acid sequence of MT-2 comprised 60 amino acids and showed approximately 62.3% identity to human metallothionein. Its promoter region was amplified by thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR). The MT-2 gene consists of 3 exons and 2 introns, extending approximately 900 bp of genomic sequence. Phylogenetic analysis clearly demonstrated that MT-2 formed a clade with fish metallothionein. The promoter region contained 5 putative metal-regulatory elements (MREs) and 1 TATA box. Real-time quantitative RT-PCR analysis revealed that MT-2 transcripts were significantly increased in the brain and gills and were stable in the muscles, liver, and trunk kidney in Cd2+-stimulated fish. Western blotting analysis demonstrated that the protein of the MT-2 gene was expressed mainly in the gills, liver, heart, trunk kidney, muscle, and intestine; it was weakly detected in the brain and head kidney. Moreover, the MT-2 protein was immunohistochemically detected in the cytoplasm in the liver and trunk kidney. All the above results revealed that the mandarin fish MT-2 would be a useful biomarker for metal pollution. (C) 2008 Published by Elsevier Inc.

Structure, organization and expression of common carp (Cyprinus carpio L.) SLP-76 gene

SLP-76 is an important member of the SLP-76 family of adapters, and it plays a key role in TCR signaling and T cell function. Partial cDNA sequence of SLP-76 of common carp (Cyprinus carpio L.) was isolated from thymus cDNA Library by the method of suppression subtractive hybridization (SSH). Subsequently, the full length cDNA of carp SLP-76 was obtained by means of 3' RACE and 5' RACE, respectively. The full Length cDNA of carp SLP-76 was 2007 bp, consisting of a T-terminal untranslated region (UTR) of 285 bp, a T-terminal. UTR of 240 bp, and an open reading frame of 1482 bp. Sequence comparison showed that the deduced amino acid sequence of carp SLP-76 had an overall similarity of 34-73% to that of other species homotogues, and it was composed of an NH2-terminal domain, a central proline-rich domain, and a C-terminal SH2 domain. Amino acid sequence analysis indicated the existence of a Gads binding site R-X-X-K, a 10-aa-long sequence which binds to the SH3 domain of LCK in vitro, and three conserved tyrosine-containing sequence in the NH2-terminal domain. Then we used PCR to obtain a genomic DNA which covers the entire coding region of carp SLP-76. In the 9.2 k-long genomic sequence, twenty one exons and twenty introns were identified. RT-PCR results showed that carp SLP-76 was expressed predominantly in hematopoietic tissues, and was upregulated in thymus tissue of four-month carp compared to one-year old carp. RT-PCR and virtual northern hybridization results showed that carp SLP-76 was also upregulated in thymus tissue of GH transgenic carp at the age of four-months. These results suggest that the expression level of SLP-76 gene may be related to thymocyte development in teleosts. (c) 2007 Published by Elsevier Ltd.

TRAIL in the mandarin fish Siniperca chuatsi: Gene and its apoptotic effect in HeLa cells

Tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL) is one of the TNF superfamily members, participating in many biological processes including cell proliferation and apoptotic death. In this study, a TRAIL gene was cloned from a perciform fish, the mandarin fish Siniperca chuatsi, a major cultured fish in China's aquaculture, and is named as SCTRAIL for S. chuatsi TRAIL. The full-length cDNA of SCTRAIL is 1359 bp, encoding a 283-amino-acid protein. This deduced protein contains the CYS231, a 23-mer fragment of transmembrane region, a glycosylation site and a TNF family signature, all of which are conserved among TRAIL members. SCTRAIL gene consists of six exons, with five intervening introns, spaced over approximately 9 kb of genomic sequence. Southern blotting demonstrated that the SCTRAIL gene is present as a single copy in mandarin fish genome. A 620 bp promoter region obtained by genome walking contains a number of putative transcription factor binding sites, such as Oct-1, Sp-1, NF-1, RAP-1, C/EBPaLp, NF-kappa B and AP-1. The SCTRAIL is constitutively expressed in all the analyzed tissues, as revealed by RT-PCR, which is confirmed by Western blotting analysis using polyclonal antibody against bacteria-derived recombinant SCTRAIL protein. As an apoptosis-inducing ligand, the overexpression of SCTRAIL but not the mutant SCTRAIL-C203S in HeLa cells induced changes characteristic of apoptosis, including chromatin condensation, nucleus fragmentation, DNA ladder, and increase of sub-G0/G1 cells in FACS analysis. (c) 2007 Elsevier Ltd. All rights reserved.

Molecular cloning, promoter analysis and induced expression of the complement component C9 gene in the grass carp Ctenopharyngodon idella

Complement-mediated killing of pathogens through lytic pathway is an important effector mechanism of innate immune response. C9 is the ninth member of complement components, creating the membrane attack complex (MAC). In the present study, a putative cDNA sequence encoding the 650 amino acids of C9 and its genomic organization were identified in grass carp Ctenopharyngodon idella. The deduced amino acid sequence of grass carp C9 (gcC9) showed 48% and 38.5% identity to Japanese flounder and human C9, respectively. Domain search revealed that gcC9 contains a LDL receptor domain, an EGF precursor domain, a MACPF domain and two TSP domain located in the N-terminal and C-terminal, respectively. Phylogenetic analysis demonstrated that gcC9 is clustered in a same clade with Japanese flounder, pufferfish and rainbow trout C9. The gcC9 gene consists of 11 exons with 10 introns, spacing over approximately 7 kb of genomic sequence. Analysis of gcC9 promoter region revealed the presence of a TATA box and some putative transcription factor such as C/EBP, HSF, NF-AT, CHOP-C, HNF-3B, GATA-2, IK-2, EVI- 1, AP-1, CP2 and OCT-1 binding sites. The first intron region contains C/EBPb, HFH-1 and Oct-1 binding sites. RT-PCR and Western blotting analysis demonstrated that the mRNA and protein of gcC9 gene have similar expression patterns, being constitutively expressed in all organs examined of healthy fish, with the highest level in hepatopancreas. By real-time quantitative RT-PCR analysis, gcC9 transcripts were significantly up-regulated in head kidney, spleen, hepatopancreas and down-regulated in intestine from inactivated fish bacterial pathogen Flavobacterium columnare-stimulated fish, demonstrating the role of C9 in immune response. (c) 2007 Elsevier B.V. All rights reserved.

Intelectin gene from the grass carp Ctenopharyngodon idella: cDNA cloning, tissue expression, and immunohistochemical localization

The cDNA encoding grass carp intelectin was isolated from a head kidney cDNA library, and termed gcIntL. The deduced amino acid sequence of gcIntL consists of 318 amino acids, and about 55% identical and 74% similar to human intelectin, which is a new type of lectin recognizing galactofuranose, and plays a role in the recognition of bacteria-specific components in animal hosts. The gcIntL gene consists of seven exons and six introns, spacing over approximately 3 kb of genomic sequence. Phylogenetic analysis clearly demonstrated that the gcIntL formed a clade with Danio rerio intelectin and 35 kDa serum lectin. By real-time quantitative RT-PCR analysis, gcIntL transcripts were significantly induced in head kidney, trunk kidney, spleen, and intestine from LPS-stimulated fish. RT-PCR and Western blotting analysis demonstrated that the mRNA and protein of gcIntL gene have the same expression pattern, and both were detected in brain, gill, intestine, head kidney, trunk kidney, spleen, and heart. Furthermore, gcIntL protein could be detected in gill, intestine, trunk kidney, head kidney, spleen, heart, and brain including medulla oblongata and optic lobe, as determined by immunohistochemistry. This is the first report of intelectin expression pattern in fish, and of recombinant gcIntL and polyclonal antibody against gcIntL. (c) 2006 Elsevier Ltd. All rights reserved.

Cloning, characterization and promoter analysis of common carp hairy/Enhancer-of-split-related gene, her6

Some members of hairy/Enhancer-of-split-related gene (HES) family have important effects on axial mesoderm segmentation and the establishment and maintenance of the somite fringe. In fishes. the her6 gene, a member of the HES family, is the homologue Of heS1 in mammals and chicken. In this study, the her6 gene and its full-length cDNA from the common carp (Cyprinus carpio) were isolated and characterized. The genomic sequence of common carp her6 is approximately 1.7 kb. with four exons and three introns, and the full-length cDNA of 1314 bp encodes a Putative polypeptide of 271 amino acids. To analyse the promoter sequence of common carp her6. sequences of various lengths upstream from the transcription initiation site of her6 were fused to enhanced green fluorescent. protein gene (eGFP) and introduced into zebrafish embryos by microinjection to generate transgenic embryos. Our results show that the upstream sequence of 500 bp can direct highly efficient and tissue-specific expression of eGFP in zebrafish embryos. whereas a fragment of 200 bp containing the TATA box and a partial suppressor of hairless paired site sequence (SPS) is not sufficient to drive eGFP expression in zebrafish embryos.

Characterization and application of monoclonal antibodies against turbot (Scophthalmus maximus) rhabdovirus

Five monoclonal antibodies (mAbs), 1G8, 1H9, 2D2, 2D3, and 2F5, against Scophthalmus maximus rhabdovirus (SMRV) were prepared. Characterization of the mAbs included indirect enzyme-linked immunosorbent assay, isotyping, viral inhibition assay, immunofluorescence staining of virus-infected cell cultures, and Western blot analysis. Isotyping revealed that 1G8 and 1H9 were of the IgG2b subclass and that the other three were IgM. 2D2, 2D3, and 2F5 partially inhibited SMRV infection in epithelioma. papulosum cyprinid (EPC) cell culture. Western blotting showed that all five mAbs could react with two SMRV proteins with molecular masses of approximately 30 kDa (P) and 26 kDa (M). These two proteins were localized within the cytoplasm of SMRV-infected EPC cells by immunofluorescence assay. Also, progressive foci of viral replication in cell cultures were monitored from 6 to 24 h, using mAb 2D3 as the primary antibody. A flow cytometry procedure was used to detect and quantify SMRV-infected (0.01 PFU/cell) EPC cells with mAb 2D3, and 10.8% of cells could be distinguished as infected 36 h postinfection. Moreover, mAb 2D3 was successfully applied for the detection of viral antigen in cryosections from flounder tissues by immunohistochemistry tests.

The first non-mammalian CXCR3 in a teleost fish: Gene and expression in blood cells and central nervous system in the grass carp (Ctenopharyngodon idella)

A novel fish chemokine receptor gene, chemokine (C-X-C motif) receptor 3 (CXCR3)-like was isolated from the grass carp Ctenopharyngodon idella , with its full-length genomic sequence. The cDNA of grass carp CXCR3-like (gcCXCR3-like) consists of 1261 bp with a 49bp 5'-UTR and a 189 bp 3'-UTR. An open reading frame of 1023 bp encodes a 341-amino acid peptide, with seven transmembrane helices. The deduced amino acid sequence showed the same sequence identities (37.8%) with its counterparts in goat and human. The gcCXCR3-like gene consists of two exons, with one intervening intron, spaced over approximately 2 kb of genomic sequence. Phylogenetic analyses clearly demonstrated that the gcCXCR3-like resembles the CXCR3s of other vertebrates. Real-time PCR analysis showed that gcCXCR3-like was expressed in all tested organs except heart and the expression level of gcCXCR3-like was highest in brain. Flow cytometric analyses showed the positive rate of labelled leukocytes from the healthy grass carp was 17.3%, and the labelled leukocytes were divided into three types by cell sorting. Immunohistochemical localization revealed that gcCXCR3-like expressed in whole brain regions including cerebel, diencephalon, medulla oblongata, optic lobe, and rhinencephalon, and that the labelled leukocytes are actually populations of monocyte and/or phagocyte, lymphocyte and the granulocyte. It is considered that fish CXCR expression and their function may need to be investigated in both nervous and immune systems. (c) 2006 Elsevier Ltd. All rights reserved.

Characterization and expression analysis of TNF-related apoptosis inducing ligand (TRAIL) in grass carp Ctenopharyngodon idella

TRAIL (Apo2 ligand) described as a type II transmembrame protein belonging to the TNF superfamily can induce apoptotic cell death in a variety of cell types. In the present study, a putative cDNA sequence encoding the 299 amino acids of TRAIL (GC-TRAIL) and its genomic organization were identified in grass carp Ctenopharyngodon idella. The predicted GC-TRAIL sequence showed 44 and 41% identities to chicken and human TRAILs, respectively. In a domain search, a tumor necrosis factor homology domain (THD) was identified in the C-terminal portion of TRAILs. The GC-TRAIL gene consists of five exons, with four intervening introns, spaced over approximately 4 kb of genomic sequence. Analysis of GC-TRAlL promoter region revealed the presence of a number of putative transcription factor binding sites, such as Sp1, NF-kappaB, AP-1, GATA, NFAT, HNF, STAT, P53 and IRFI sequences which are important for the expression of other TNF family members. Phylogenetic analysis placed GC-TRAIL and the putative zebrafish (Danio rerio) TRAIL obtained from searching the zebrafish database into one separate cluster near mammalian TRAIL genes, but apart from the reported zebrafish TRAIL-like protein, indicating that the GC-TRAIL is an authentic fish TRAIL. Expression analysis revealed that GC-TRAIL is expressed in many tissues, such as in gills, liver, trunk kidney, head kidney, intestine and spleen. (c) 2005 Elsevier B.V. All rights reserved.

Molecular structure and expression patterns of flounder (Paralichthys olivaceus) Myf-5, a myogenic regulatory factor

Myf-5, a member of the myogenic regulatory factors (MRF), has been shown to be expressed in muscle precursors in early stage zebrafish embryos. The MRFs, including MyoD, Myf-5, Myogenin and MR-F4, belong to the basic Helix-Loop-Helix transcription factors that contain a conserved basic Helix-Loop-Helix (bHLH) domain. To better understand the role of Myf-5 in the development of fish muscles, we have isolated the Myf-5 genomic sequence and cDNA from Flounder (Paralichthys olivaceus), and analyzed its structures and patterns of expression. Promoter analysis identified several putative transcription factor binding sites such as an E-box, NF-Y sites that might confer muscle-specific expression. Myf-5 transcripts were first detected in the paraxial mesoderm that gives rise to slow muscles. During somitogenesis, Myf-5 expression was found in developing somites. Myf-5 expression decreased gradually in somites in the anterior region, but remained strong in the newly formed somites. In the hatching stage, the expression was also detected in other muscle cells such as head muscle and fin muscle. In the growing fish, RT-PCR results showed that Myf-5 was expressed in the skeletal muscle and intestine. (c) 2006 Elsevier Inc. All rights reserved.

High-coverage sequencing and annotated assemblies of the budgerigar genome.

BACKGROUND: Parrots belong to a group of behaviorally advanced vertebrates and have an advanced ability of vocal learning relative to other vocal-learning birds. They can imitate human speech, synchronize their body movements to a rhythmic beat, and understand complex concepts of referential meaning to sounds. However, little is known about the genetics of these traits. Elucidating the genetic bases would require whole genome sequencing and a robust assembly of a parrot genome. FINDINGS: We present a genomic resource for the budgerigar, an Australian Parakeet (Melopsittacus undulatus) -- the most widely studied parrot species in neuroscience and behavior. We present genomic sequence data that includes over 300× raw read coverage from multiple sequencing technologies and chromosome optical maps from a single male animal. The reads and optical maps were used to create three hybrid assemblies representing some of the largest genomic scaffolds to date for a bird; two of which were annotated based on similarities to reference sets of non-redundant human, zebra finch and chicken proteins, and budgerigar transcriptome sequence assemblies. The sequence reads for this project were in part generated and used for both the Assemblathon 2 competition and the first de novo assembly of a giga-scale vertebrate genome utilizing PacBio single-molecule sequencing. CONCLUSIONS: Across several quality metrics, these budgerigar assemblies are comparable to or better than the chicken and zebra finch genome assemblies built from traditional Sanger sequencing reads, and are sufficient to analyze regions that are difficult to sequence and assemble, including those not yet assembled in prior bird genomes, and promoter regions of genes differentially regulated in vocal learning brain regions. This work provides valuable data and material for genome technology development and for investigating the genomics of complex behavioral traits.

Copy number variation and Huntington’s disease

Tese de doutoramento, Ciências Biomédicas (Neurociências), Universidade de Lisboa, Faculdade de Medicina, 2014