Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

Alveolar macrophages and lung dendritic cells sense RNA and drive mucosal IgA responses.

The mechanisms regulating systemic and mucosal IgA responses in the respiratory tract are incompletely understood. Using virus-like particles loaded with single-stranded RNA as a ligand for TLR7, we found that systemic vs mucosal IgA responses in mice were differently regulated. Systemic IgA responses following s.c. immunization were T cell independent and did not require TACI or TGFbeta, whereas mucosal IgA production was dependent on Th cells, TACI, and TGFbeta. Strikingly, both responses required TLR7 signaling, but systemic IgA depended upon TLR7 signaling directly to B cells whereas mucosal IgA required TLR7 signaling to lung dendritic cells and alveolar macrophages. Our data show that IgA switching is controlled differently according to the cell type receiving TLR signals. This knowledge should facilitate the development of IgA-inducing vaccines.

High throughput sequencing and proteomics to identify immunogenic proteins of a new pathogen: the dirty genome approach.

With the availability of new generation sequencing technologies, bacterial genome projects have undergone a major boost. Still, chromosome completion needs a costly and time-consuming gap closure, especially when containing highly repetitive elements. However, incomplete genome data may be sufficiently informative to derive the pursued information. For emerging pathogens, i.e. newly identified pathogens, lack of release of genome data during gap closure stage is clearly medically counterproductive. We thus investigated the feasibility of a dirty genome approach, i.e. the release of unfinished genome sequences to develop serological diagnostic tools. We showed that almost the whole genome sequence of the emerging pathogen Parachlamydia acanthamoebae was retrieved even with relatively short reads from Genome Sequencer 20 and Solexa. The bacterial proteome was analyzed to select immunogenic proteins, which were then expressed and used to elaborate the first steps of an ELISA. This work constitutes the proof of principle for a dirty genome approach, i.e. the use of unfinished genome sequences of pathogenic bacteria, coupled with proteomics to rapidly identify new immunogenic proteins useful to develop in the future specific diagnostic tests such as ELISA, immunohistochemistry and direct antigen detection. Although applied here to an emerging pathogen, this combined dirty genome sequencing/proteomic approach may be used for any pathogen for which better diagnostics are needed. These genome sequences may also be very useful to develop DNA based diagnostic tests. All these diagnostic tools will allow further evaluations of the pathogenic potential of this obligate intracellular bacterium.

Alteration of glucose metabolism in cultured astrocytes after AQP9-small interference RNA application.

Aquaglyceroporin-9 (AQP9) facilitates diffusion of water and energy substrates such as glycerol and monocarboxylates. AQP9 is present in plasma membrane and mitochondria of astrocytes and catecholaminergic neurons, suggesting that it plays a role in the energetic status of these cells. Using specific small interference RNA directed against AQP9 in astrocyte cultures, we showed that glycerol uptake is decreased which is associated with an increase in glucose uptake and oxidative metabolism. Our results not only confirm the presence of AQP9 in astrocytes but also suggest that changes in AQP9 expression alter glial energy metabolism.

RNA-based gene duplication: mechanistic and evolutionary insights.

Gene copies that stem from the mRNAs of parental source genes have long been viewed as evolutionary dead-ends with little biological relevance. Here we review a range of recent studies that have unveiled a significant number of functional retroposed gene copies in both mammalian and some non-mammalian genomes. These studies have not only revealed previously unknown mechanisms for the emergence of new genes and their functions but have also provided fascinating general insights into molecular and evolutionary processes that have shaped genomes. For example, analyses of chromosomal gene movement patterns via RNA-based gene duplication have shed fresh light on the evolutionary origin and biology of our sex chromosomes.

Exome sequencing identifies INPPL1 mutations as a cause of opsismodysplasia.

Opsismodysplasia (OPS) is a severe autosomal-recessive chondrodysplasia characterized by pre- and postnatal micromelia with extremely short hands and feet. The main radiological features are severe platyspondyly, squared metacarpals, delayed skeletal ossification, and metaphyseal cupping. In order to identify mutations causing OPS, a total of 16 cases (7 terminated pregnancies and 9 postnatal cases) from 10 unrelated families were included in this study. We performed exome sequencing in three cases from three unrelated families and only one gene was found to harbor mutations in all three cases: inositol polyphosphate phosphatase-like 1 (INPPL1). Screening INPPL1 in the remaining cases identified a total of 12 distinct INPPL1 mutations in the 10 families, present at the homozygote state in 7 consanguinous families and at the compound heterozygote state in the 3 remaining families. Most mutations (6/12) resulted in premature stop codons, 2/12 were splice site, and 4/12 were missense mutations located in the catalytic domain, 5-phosphatase. INPPL1 belongs to the inositol-1,4,5-trisphosphate 5-phosphatase family, a family of signal-modulating enzymes that govern a plethora of cellular functions by regulating the levels of specific phosphoinositides. Our finding of INPPL1 mutations in OPS, a severe spondylodysplastic dysplasia with major growth plate disorganization, supports a key and specific role of this enzyme in endochondral ossification.

TLR3 and Rig-like receptor on myeloid dendritic cells and Rig-like receptor on human NK cells are both mandatory for production of IFN-gamma in response to double-stranded RNA.

Cross-talk between NK cells and dendritic cells (DCs) is critical for the potent therapeutic response to dsRNA, but the receptors involved remained controversial. We show in this paper that two dsRNAs, polyadenylic-polyuridylic acid and polyinosinic-polycytidylic acid [poly(I:C)], similarly engaged human TLR3, whereas only poly(I:C) triggered human RIG-I and MDA5. Both dsRNA enhanced NK cell activation within PBMCs but only poly(I:C) induced IFN-gamma. Although myeloid DCs (mDCs) were required for NK cell activation, induction of cytolytic potential and IFN-gamma production did not require contact with mDCs but was dependent on type I IFN and IL-12, respectively. Poly(I:C) but not polyadenylic-polyuridylic acid synergized with mDC-derived IL-12 for IFN-gamma production by acting directly on NK cells. Finally, the requirement of both TLR3 and Rig-like receptor (RLR) on mDCs and RLRs but not TLR3 on NK cells for IFN-gamma production was demonstrated using TLR3- and Cardif-deficient mice and human RIG-I-specific activator. Thus, we report the requirement of cotriggering TLR3 and RLR on mDCs and RLRs on NK cells for a pathogen product to induce potent innate cell activation.

Genomic and transcriptomic analysis of the streptomycin-dependent Mycobacterium tuberculosis strain 18b.

The ability of Mycobacterium tuberculosis to establish a latent infection (LTBI) in humans confounds the treatment of tuberculosis. Consequently, there is a need to discover new therapeutic agents that can kill M. tuberculosis both during active disease and LTBI. The streptomycin-dependent strain of M. tuberculosis, 18b, provides a useful tool for this purpose since upon removal of streptomycin (STR) it enters a non-replicating state that mimics latency both in vitro and in animal models. The 4.41 Mb genome sequence of M. tuberculosis 18b was determined and this revealed the strain to belong to clade 3 of the ancient ancestral lineage of the Beijing family. STR-dependence was attributable to insertion of a single cytosine in the 530 loop of the 16S rRNA and to a single amino acid insertion in the N-terminal domain of initiation factor 3. RNA-seq was used to understand the genetic programme activated upon STR-withdrawal and hence to gain insight into LTBI. This revealed reconfiguration of gene expression and metabolic pathways showing strong similarities between non-replicating 18b and M. tuberculosis residing within macrophages, and with the core stationary phase and microaerophilic responses. The findings of this investigation confirm the validity of 18b as a model for LTBI, and provide insight into both the evolution of tubercle bacilli and the functioning of the ribosome.

Cloning, sequencing and characterization of a chitin synthase gene fragment from the fungus Phascolomyces articulosus /

Phascolomyces articulosus genomic DNA was isolated from 48 h old hyphae and was used for amplification of a chitin synthase fragment by the polymerase chain reaction method. The primers used in the amplification corresponded to two widely conserved amino acid regions found in chitin synthases of many fimgi. Amphfication resulted in four bands (820, 900, 1000 and 1500 bp, approximately) as visualized in a 1.2% agarose gel. The lowest band (820 bp) was selected as a candidate for chitin synthase because most amplified regions from other fimgi so far exhibited similar sizes (600-750 bp). The selected fragment was extracted from the gel and cloned in the Hinc n site of pUC19. The derived plasmid and insert were designated ^\5C\9'PaCHS and PaCHS respectively. The plasmid pUC19-PaC/fS was digested by several restriction enzymes and was found to contain BamHl and HincU sites. Sequencing of PaCHS revealed two intron sequences and a total open reading frame of 200 amino acids. The derived polypeptide was compared with other related sequences from the EMBL database (Heidelberg, Germany) and was matched to 36 other fiilly or partially sequenced fimgal chitin synthase genes. The closest resemblance was with two genes (74.5% and 73.1% identity) from Rhizopus oligosporus. Southern hybridization with the cloned fragment as a probe to the PCR reaction showed a strong signal at the fragment selected for cloning and weaker signals at the other two fragments. Southern hybridization with partially digested Phascolomyces articulosus genomic DNA showed a single band. The amino acid sequence was compared with sequences from other chitin synthase gene classes using the CLUSTALW program. The chitin synthase fragment from Phascolomyces articulosus was initially grouped in class n along with chitin synthase fragments from Rhizopus oligosporus and Phycomyces blakesleeanus which also belong to the same class, Zygomycetes. Bootstrap analysis using the neighbor-joining method available by CLUSTALW verified such classification. Comparison of PaCHS revealed conservation of intron positions that are characteristic of chitin synthase gene fragments of zygomycetous fungi.

Construction of bovine adenovirus type 3 E1 and E3, substitution plasmids and the sequencing analysis of DNA pol and pTP /

The relative ease to concentrate and purify adenoviruses, their well characterized mid-sized genome, and the ability to delete non-essential regions from their genome to accommodate foreign gene, made adenoviruses a suitable candidate for the construction of vectors. The use of adenoviral vectors in gene therapy, vaccination, and as a general vector system for expressing foreign genes have been documented for some time. In this study, the objective was to rescue a BAV3 E1 or E3 recombinant vector carrying the kanamycin resistant gene, a dominant selectable marker with useful applications in studying vectored gene expression in mammalian cells. To accomplish the objective of this study, more information about BAV3 DNA sequences was required in order to make the manipulation of the virus genome accessible. Therefore, sequencing of the BAV3 genome from 1 1 .7% to 30.8% was carried out. Analysis of the determined sequences revealed the primary structure of important viral gene products coded by E2 including BAV3 DNA pol and precursor to terminal protein. Comparative analysis of these proteins with their counterparts from human and non human adenoviruses revealed important insights as to the evolutionary lineage of BAV3. In order to insert the kanamycin resistance gene in either E1 or E3, it was necessary to delete BAV3 sequences to accommodate the foreign gene so as not to exceed the limit of the packaging capacity of the virus. To construct a recombinant BAV3 in which a foreign gene was inserted in the deleted E1 region, an E1 shuttle vector was constructed. This involved the deletion from the viral sequences a region between 1.3% to 9% and inserting the kanamycin resistance gene to replace the deletion. The E1 shuttle vector contained the left (0%- 53.9%) segment of the genome and was expected to generate BAV3 recombinants that can be grown and propagated in cells that can complement the missing E1 functions. To construct a similar shuttle vector for E3 deletion, DNA sequences extending from 78.9% to 82.5% (1281 bp) were deleted from within the E3 region that had been cloned into a plasmid vector. The deleted region corresponds to those that have been shown to be non-essential for viral replication in cell culture. The resulting plasmid was used to construct another recombinant plasmid with BAV3 DNA sequences extending from 37.1% to 100% and with a deletion of E3 sequences that were replaced by kanamycin resistance gene. This shuttle plasmid was used in cotransfections with digested viral DNA in an attempt to rescue a recombinant BAV3 carrying the kanamycin resistance gene to replace the deleted E3. In spite of repeated attempts of transfection, El or E3 recombinant BAV3 were not isolated. It seems that other approaches should be applied to make a final conclusion on BAV3 infectivity.

An investigation into differences between indoleacetic acid and fusicoccin in their influence on RNA synthesis, protein synthesis and growth in Avena coleoptile tissue

Growth stimulation of Avena coleoptile tissue by indoleacetic acid (IAA) and fusicoccin (FC) was compared by measuring both their influence on RNA and protein synthesis during IAA or FC stimulated growth. FC stimulated growth more than IAA during the initial four hour exposure, after which the growth rate gradually declined to the control rate. FC, but not IAA, increased the uptake of 3H-Ieucine into tissue and the specific radioactivity of extracted protein. Cycloheximide inhibited the incorporation of 3H-Ieucine into protein by approximately 60% to 70% in all cases. In the presence of cycloheximide 3H-radioactivity accumulated in FC-treated tissue, whereas IAA did not seem to influence 3H-accumulation. These results suggest that FC stimulated leucine uptake into the tissue and that increased specific activity of coleoptile protein is due to increased leucine uptake, not an increased rate of protein synthesis. There was no measurable influence of IAA and/or FC on RNA and protein synthesis during the initial hours of a growth stimulation. Inhibitors of RNA and protein synthesis, actinomycin D and cycloheximide, respectively, severely inhibited IAA enhanced growth but only partially inhibited FC stimulated growth. The data are consistent with suggestions that a rapidly turning over protein participates in IAA stimulated growth, and that a continual synthesis of RNA and proteins is an absolute requirement for a long term growth response to IAA. On the contrary, FC-stimulated growth exhibited less dependency on the transcription and translation processes. The data are consistent with proposals suggesting different sites of action for FC and IAA stimulated growth. l?hen compared to CO2-free air, CO2 at 300 ppm had no significant influence on coleoptile growth and protein synthesis in the presence or absence of lAA or FC. Also, I mM malate, pH 6.0 did not influence growth of coleoptiles in the presence or absence of lAA. This result was obtained despite reports indicating that 300 ppm CO2 or I mM malate stimulates growth and protein synthesis. This lack of difference between CO2-treated and untreated tissue could indicate either that the interstitial space CO2 concentration is not actually different in the two treatments due to significant endogenous respiratory CO2 or else the data would suggest a very loose coupling between dark CO2 fixation and growth. IAA stimulated the in vivo fixation of 14c-bicarbonate (NaHI4c03) by about 25% and the addition of cycloheximide caused an inhibition of bicarbonate fixation within 30 min. Cycloheximide has also been reported to inhibit IAA-stimulated H+ excretion. These data are consistent with the acid growth theory and suggest that lAA stimulated growth involves dark CO2 fixation. The roles of dark CO2 fixation in lAA-stimulated growth are discussed.

Locating and sequencing of the E3 and neighbouring regions in bovine advenovirus type 2

Adenoviruses are nonenveloped icosahedral shaped particles. The double stranded DNA viral genome is divided into 5 major early transcription units, designated E1 A, E1 B, and E2 to E4, which are expressed in a regulated manner soon after infection. The gene products of the early region 3 (E3), shown to be nonessential for viral replication in vitro, are believed to be involved in counteracting host immunosurveillance. In order to sequence the E3 region of Bovine adenovirus type 2 (BAV2) it was necessary to determine the restriction map for the plasmid pEA48. A physical restriction endonuclease map for BamHl, Clal, Eco RI, Hindlll, Kpnl, Pstt, Sail, and Xbal was constructed. The DNA insert in pEA48 was determined to be viral in origin using Southern hybridization. A human adenovirus type 5 recombinant plasmid, containing partial DNA fragments of the two transcription units L4 and L5 that lie just outside the E3, was used to localize this region. The recombinant plasmid pEA was subcloned to facilitate sequencing. The DNA sequences between 74.8 and 90.5 map units containing the E3, the hexon associated protein (pVIII), and the fibre gene were determined. Homology comparison revealed that the genes for the hexon associated pV11I and the fibre protein are conserved. The last 70 amino acids of the BAV2 pV11I were the most conserved, showing a similarity of 87 percent with Ad2 pV1I1. A comparison between the predicted amino acid sequences of BAV2 and Ad40, Ad41 , Ad2 and AdS, revealed that they have an identical secondary structure consisting of a tail, a shaft and a knob. The shaft is composed of 22, 15 amino acid motifs, with periodic glycines and hydrophobic residues. The E3 region was found to consist of about 2.3 Kbp and to encode four proteins that were greater than 60 amino acids. However, these four open reading frames did not show significant homology to any other known adenovirus DNA or protein sequence.

The cloning and reconstitution of a bovine adenovirus type 2 E3 deletion mutant and the sequencing and analysis of the early 4 region

Recombinant Adenoviruses (Ads) have been shown to have potential applications in three areas: gene therapy, high level protein expression and recombinant vaccines.' At least three different locations within the Ad genome can be deleted and subsequently used for the insertion of foreign sequences. These include the Early 3 (E3), Early 1 (E1) and Early 4 (E4) regions. Viral vectors of this type have been well studied in Human Ads 2 and 5, however one has not yet been constructed for Bovine Adenovirus Type 2 (BAV2). The E3 region is located between 76.6 and 86 m.u. on the r-strand and is transcribed in a rightward direction. The gene products of the Early 3 region (E3) have been shown to be non-essential for viral replication, in vitro, but are required for host immunosurveillance. This study represents the cloning and reconstitution of a BAV2 E3 deletion mutant. A deletion of 1800bp was made within the E3 region of BAV2 and the thymidine kinase gene was subsequently inserted in the deleted area . . The plasmid pdlE3-4tk1 (23.4Kbp) was constructed and used to to facilitate homologous recombination with the wild type BAV2 to produce a mutant. Southern Blotting and Hybridization results suggest the presence of a BAV2 E3 deletion mutant with thymidine kinase sequences present. The E4 region of Human Adenovirus types 2 and 5 is located at the extreme right end of the genome (91.3 map units - 99.1 map units) and is transcribed in a leftward direction giving rise to a complicated set of differentially spliced mRNAs. Essentially there are 7 open reading frames (ORFs) encoding for at least 7 polypeptides. The gene products encoded by the E4 region have been shown to be essential for the expression of late viral genes, host cell shutoff and normal viral growth. We have cloned and sequenced the right end segment between 90.5 map units and 100 map units of the BAV2 genome. The results show several open reading frames which encode polypeptides exhibiting homology to three polypeptides encoded by the E4 region of human adenovirus type 2. These include the 14kDa protein encoded by ORF1, the 34kDa protein encoded by ORF6 and the 13kDa protein encoded by ORF3. The nucleotide sequence, restriction enzyme map, and ORF map of the E4 region could be very useful in future molecular manipulation of this region and could possibly explain the slow growth rate of BAV2 in MDBK cells.

(A) Optimization of solid body phase synthesis of RNA using the Cpep chemistry. (B) Synthesis of BODIPY labeled oligonucleotides

(A) Solid phase synthesis of oligonucleotides are well documented and are extensively studied as the demands continue to rise with the development of antisense, anti-gene, RNA interference, and aptamers. Although synthesis of RNA sequences faces many challenges, most notably the choice of the 2' -hydroxy protecting group, modified 2' -O-Cpep protected ribonucleotides were synthesized as alternitive building blocks. Altering phosphitylation procedures to incorporate 3' -N,N-diethyl phosphoramidites enhanced the overall reactivity, thus, increased the coupling efficiency without loss of integrety. Furthermore, technical optimizations of solid phase synthesis cycles were carried out to allow for successful synthesis of a homo UIO sequences with a stepwise coupling efficiency reaching 99% and a final yield of 91 %. (B) Over the past few decades, dipyrrometheneboron difluoride (BODIPY) has gained recognition as one of the most versatile fluorophores. Currently, BODIPY labeling of oligonucleotides are carried out post-synthetically and to date, there lacks a method that allows for direct incorporation of BODIPY into oligonucleotides during solid phase synthesis. Therefore, synthesis of BODIPY derived phosphoramidites will provide an alternative method in obtaining fluorescently labelled oligonucleotides. A method for the synthesis and incorporation of the BODIPY analogues into oligonucleotides by phosphoramidite chemistry-based solid phase DNA synthesis is reported here. Using this approach, BODIPY-labeled TlO homopolymer and ISIS 5132 were successfully synthesized.