RNA-templated DNA ligation for transcript analysis

Ligase-mediated gene detection has proven valuable for detection and precise distinction of DNA sequence variants. We have recently shown that T4 DNA ligase can also be used to distinguish single nucleotide variants of RNA sequences. Here we describe parameters that influence RNA-templated DNA ligation by T4 DNA ligase. The reaction proceeds much more slowly, requiring more enzyme, compared to ligation of the same oligonucleotides hybridized to the corresponding DNA sequence. The reaction is inhibited at high concentrations of ATP and NaCl and both magnesium and manganese ions can support the reaction. We define reaction conditions where 80% of RNA target molecules can template a diagnostic ligation reaction. Ligase-mediated RNA detection should provide a useful mechanism for sensitive and accurate detection and distinction of RNA sequence variants.

Phylogenetic Analysis Of The Microbial Community In Hypersaline Petroleum Produced Water From The Campos Basin.

In this work the archaea and eubacteria community of a hypersaline produced water from the Campos Basin that had been transported and discharged to an onshore storage facility was evaluated by 16S recombinant RNA (rRNA) gene sequence analysis. The produced water had a hypersaline salt content of 10 (w/v), had a carbon oxygen demand (COD) of 4,300 mg/l and contains phenol and other aromatic compounds. The high salt and COD content and the presence of toxic phenolic compounds present a problem for conventional discharge to open seawater. In previous studies, we demonstrated that the COD and phenolic content could be largely removed under aerobic conditions, without dilution, by either addition of phenol degrading Haloarchaea or the addition of nutrients alone. In this study our goal was to characterize the microbial community to gain further insight into the persistence of reservoir community members in the produced water and the potential for bioremediation of COD and toxic contaminants. Members of the archaea community were consistent with previously identified communities from mesothermic reservoirs. All identified archaea were located within the phylum Euryarchaeota, with 98 % being identified as methanogens while 2 % could not be affiliated with any known genus. Of the identified archaea, 37 % were identified as members of the strictly carbon-dioxide-reducing genus Methanoplanus and 59 % as members of the acetoclastic genus Methanosaeta. No Haloarchaea were detected, consistent with the need to add these organisms for COD and aromatic removal. Marinobacter and Halomonas dominated the eubacterial community. The presence of these genera is consistent with the ability to stimulate COD and aromatic removal with nutrient addition. In addition, anaerobic members of the phyla Thermotogae, Firmicutes, and unclassified eubacteria were identified and may represent reservoir organisms associated with the conversion hydrocarbons to methane.

Expression analysis of putative vitellogenin and lipophorin receptors in honey bee (Apis mellifera L.) queens and workers

Two members of the low density lipoprotein receptor (LDLR) family were identified as putative orthologs for a vitellogenin receptor (Amvgr) and a lipophorin receptor (Amlpr) in the Apis mellifera genome. Both receptor sequences have the structural motifs characteristic of LDLR family members and show a high degree of similarity with sequences of other insects. RT-PCR analysis of Amvgr and Amlpr expression detected the presence of both transcripts in different tissues of adult female (ovary, fat body, midgut, head and specifically hypopharyngeal gland), as well as in embryos. In the head RNA samples we found two variant forms of AmLpR: a full length one and a shorter one lacking 29 amino acids in the O-linked sugar domain. In ovaries the expression levels of the two honey bee LDLR members showed opposing trends: whereas Amvgr expression was upregulated as the ovaries became activated, Amlpr transcript levels gradually declined. In situ hybridization analysis performed on ovaries detected Amvgr mRNA exclusively in germ line cells and corroborated the qPCR results showing an increase in Amvgr gene expression concomitant with follicle growth. (C) 2008 Elsevier Ltd. All rights reserved.

Are Ichthyosporea animals or fungi? Bayesian phylogenetic analysis of elongation factor 1α of Ichthyophonus irregularis

Ichthyosporea is a recently recognized group of morphologically simple eukaryotes, many of which cause disease in aquatic organisms. Ribosomal RNA sequence analyses place Ichthyosporea near the divergence of the animal and fungal lineages, but do not allow resolution of its exact phylogenetic position. Some of the best evidence for a specific grouping of animals and fungi (Opisthokonta) has come from elongation factor 1alpha, not only phylogenetic analysis of sequences but also the presence or absence of short insertions and deletions. We sequenced the EF-1alpha gene from the ichthyosporean parasite Ichthyophonus irregularis and determined its phylogenetic position using neighbor-joining, parsimony and Bayesian methods. We also sequenced EF-1alpha genes from four chytrids to provide broader representation within fungi. Sequence analyses and the presence of a characteristic 12 amino acid insertion strongly indicate that I. irregularis is a member of Opisthokonta, but do not resolve whether I. irregularis is a specific relative of animals or of fungi. However, the EF-1alpha of I. irregularis exhibits a two amino acid deletion heretofore reported only among fungi. (C) 2003 Elsevier Science (USA). All rights reserved.

MOLECULAR CHARACTERIZATION AND SEQUENCE PHYLOGENETIC ANALYSIS OF SURFACE ANTIGEN 3 (SAG3) GENE OF LOCAL INDIAN ISOLATES (CHENNAI AND IZATNAGAR) OF Toxoplasma gondii

Context and objective:The molecular characterization of local isolates of Toxoplasma gondii is considered significant so as to assess the homologous variations between the different loci of various strains of parasites.Design and setting:The present communication deals with the molecular cloning and sequence analysis of the 1158 bp entire open reading frame (ORF) of surface antigen 3 (SAG3) of two Indian T. gondii isolates (Chennai and Izatnagar) being maintained as cryostock at the IVRI.Method:The surface antigen 3 (SAG3) of two local Indian isolates were cloned and sequenced before being compared with the available published sequences.Results:The sequence comparison analysis revealed 99.9% homology with the standard published RH strain sequence of T. gondii. The strains were also compared with other established published sequences and found to be most related to the P-Br strain and CEP strain (both 99.3%), and least with PRU strain (98.4%). However, the two Indian isolates had 100% homology between them.Conclusion:Finally, it was concluded that the Indian isolates were closer to the RH strain than to the P-Br strain (Brazilian strain), the CEP strain and the PRU strains (USA), with respect to nucleotide homology. The two Indian isolates used in the present study are known to vary between themselves, as far as homologies related to other genes are concerned, but they were found to be 100% homologous as far as SAG3 locus is concerned. This could be attributed to the fact that this SAG3 might be a conserved locus and thereby, further detailed studies are thereby warranted to exploit the use of this particular molecule in diagnostics and immunoprophylactics. The findings are important from the point of view of molecular phylogeny.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita.

Plant-parasitic nematodes are major agricultural pests worldwide and novel approaches to control them are sorely needed. We report the draft genome sequence of the root-knot nematode Meloidogyne incognita, a biotrophic parasite of many crops, including tomato, cotton and coffee. Most of the assembled sequence of this asexually reproducing nematode, totaling 86 Mb, exists in pairs of homologous but divergent segments. This suggests that ancient allelic regions in M. incognita are evolving toward effective haploidy, permitting new mechanisms of adaptation. The number and diversity of plant cell wall-degrading enzymes in M. incognita is unprecedented in any animal for which a genome sequence is available, and may derive from multiple horizontal gene transfers from bacterial sources. Our results provide insights into the adaptations required by metazoans to successfully parasitize immunocompetent plants, and open the way for discovering new antiparasitic strategies.

Widespread occurrence of non-canonical transcription termination by human RNA polymerase III.

Human RNA polymerase (Pol) III-transcribed genes are thought to share a simple termination signal constituted by four or more consecutive thymidine residues in the coding DNA strand, just downstream of the RNA 3'-end sequence. We found that a large set of human tRNA genes (tDNAs) do not display any T(≥4) stretch within 50 bp of 3'-flanking region. In vitro analysis of tDNAs with a distanced T(≥4) revealed the existence of non-canonical terminators resembling degenerate T(≥5) elements, which ensure significant termination but at the same time allow for the production of Pol III read-through pre-tRNAs with unusually long 3' trailers. A panel of such non-canonical signals was found to direct transcription termination of unusual Pol III-synthesized viral pre-miRNA transcripts in gammaherpesvirus 68-infected cells. Genome-wide location analysis revealed that human Pol III tends to trespass into the 3'-flanking regions of tDNAs, as expected from extensive terminator read-through. The widespread occurrence of partial termination suggests that the Pol III primary transcriptome in mammals is unexpectedly enriched in 3'-trailer sequences with the potential to contribute novel functional ncRNAs.

Mutation analysis of genes that control the G1/S cell cycle in melanoma: TP53, CDKN1A, CDKN2A, and CDKN2B.

BACKGROUND The role of genes involved in the control of progression from the G1 to the S phase of the cell cycle in melanoma tumors in not fully known. The aim of our study was to analyse mutations in TP53, CDKN1A, CDKN2A, and CDKN2B genes in melanoma tumors and melanoma cell lines METHODS We analysed 39 primary and metastatic melanomas and 9 melanoma cell lines by single-stranded conformational polymorphism (SSCP). RESULTS The single-stranded technique showed heterozygous defects in the TP53 gene in 8 of 39 (20.5%) melanoma tumors: three new single point mutations in intronic sequences (introns 1 and 2) and exon 10, and three new single nucleotide polymorphisms located in introns 1 and 2 (C to T transition at position 11701 in intron 1; C insertion at position 11818 in intron 2; and C insertion at position 11875 in intron 2). One melanoma tumor exhibited two heterozygous alterations in the CDKN2A exon 1 one of which was novel (stop codon, and missense mutation). No defects were found in the remaining genes. CONCLUSION These results suggest that these genes are involved in melanoma tumorigenesis, although they may be not the major targets. Other suppressor genes that may be informative of the mechanism of tumorigenesis in skin melanomas should be studied.

Comparative analyses of classical phenotypic method and ribosomal RNA gene sequencing for identification of medically relevant Candida species

As the distribution of Candida species and their susceptibility to antifungal agents have changed, a new means of accurately and rapidly identifying these species is necessary for the successful early resolution of infection and the subsequent reduction of morbidity and mortality. The current work aimed to evaluate ribosomal RNA gene sequencing for the identification of medically relevant Candida species in comparison with a standard phenotypic method. Eighteen reference strains (RSs), 69 phenotypically identified isolates and 20 inconclusively identified isolates were examined. Internal transcribed spaces (ITSs) and D1/D2 of the 26S ribosomal RNA gene regions were used as targets for sequencing. Additionally, the sequences of the ITS regions were used to establish evolutionary relationships. The sequencing of the ITS regions was successful for 88% (94/107) of the RS and isolates, whereas 100% of the remaining 12% (13/107) of the samples were successfully analysed by sequencing the D1/D2 region. Similarly, genotypic analysis identified all of the RS and isolates, including the 20 isolates that were not phenotypically identified. Phenotypic analysis, however, misidentified 10% (7/69) of the isolates. Phylogenetic analysis allowed the confirmation of the relationships between evolutionarily close species. Currently, the use of genotypic methods is necessary for the correct identification of Candida species.

Analysis of epidermis- and mesophyll-specific transcript accumulation in powdery mildew-inoculated wheat leaves.

Powdery mildew is an important disease of wheat caused by the obligate biotrophic fungus Blumeria graminis f. sp. tritici. This pathogen invades exclusively epidermal cells after penetrating directly through the cell wall. Because powdery mildew colonizes exclusively epidermal cells, it is of importance not only to identify genes which are activated, but also to monitor tissue specificity of gene activation. Acquired resistance of wheat to powdery mildew can be induced by a previous inoculation with the non-host pathogen B. graminis f. sp. hordei, the causal agent of barley powdery mildew. The establishment of the resistant state is accompanied by the activation of genes. Here we report the tissue-specific cDNA-AFLP analysis and cloning of transcripts accumulating 6 and 24 h after the resistance-inducing inoculation with B. graminis f. sp. hordei. A total of 25,000 fragments estimated to represent about 17,000 transcripts were displayed. Out of these, 141 transcripts, were found to accumulate after Bgh inoculation using microarray hybridization analysis. Forty-four accumulated predominantly in the epidermis whereas 76 transcripts accumulated mostly in mesophyll tissue.

Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.

Analysis of susceptibility to human immunodeficiency virus 1 (HIV-1) by comparative genetics

Summary [résumé français voir ci-dessous] From the beginning of the 20th century the world population has been confronted with the human immune deficiency virus 1 (HIV-1). This virus has the particularity to mutate fast, and could thus evade and adapt to the human host. Our closest evolutionary related organisms, the non-human primates, are less susceptible to HIV-1. In a broader sense, primates are differentially susceptible to various retrovirus. Species specificity may be due to genetic differences among primates. In the present study we applied evolutionary and comparative genetic techniques to characterize the evolutionary pattern of host cellular determinants of HIV-1 pathogenesis. The study of the evolution of genes coding for proteins participating to the restriction or pathogenesis of HIV-1 may help understanding the genetic basis of modern human susceptibility to infection. To perform comparative genetics analysis, we constituted a collection of primate DNA and RNA to allow generation of de novo sequence of gene orthologs. More recently, release to the public domain of two new primate complete genomes (bornean orang-utan and common marmoset) in addition of the three previously available genomes (human, chimpanzee and Rhesus monkey) help scaling up the evolutionary and comparative genome analysis. Sequence analysis used phylogenetic and statistical methods for detecting molecular adaptation. We identified different selective pressures acting on host proteins involved in HIV-1 pathogenesis. Proteins with HIV-1 restriction properties in non-human primates were under strong positive selection, in particular in regions of interaction with viral proteins. These regions carried key residues for the antiviral activity. Proteins of the innate immunity presented an evolutionary pattern of conservation (purifying selection) but with signals of relaxed constrain if we compared them to the average profile of purifying selection of the primate genomes. Large scale analysis resulted in patterns of evolutionary pressures according to molecular function, biological process and cellular distribution. The data generated by various analyses served to guide the ancestral reconstruction of TRIM5a a potent antiviral host factor. The resurrected TRIM5a from the common ancestor of Old world monkeys was effective against HIV-1 and the recent resurrected hominoid variants were more effective against other retrovirus. Thus, as the result of trade-offs in the ability to restrict different retrovirus, human might have been exposed to HIV-1 at a time when TRIM5a lacked the appropriate specific restriction activity. The application of evolutionary and comparative genetic tools should be considered for the systematical assessment of host proteins relevant in viral pathogenesis, and to guide biological and functional studies. Résumé La population mondiale est confrontée depuis le début du vingtième siècle au virus de l'immunodéficience humaine 1 (VIH-1). Ce virus a un taux de mutation particulièrement élevé, il peut donc s'évader et s'adapter très efficacement à son hôte. Les organismes évolutivement le plus proches de l'homme les primates nonhumains sont moins susceptibles au VIH-1. De façon générale, les primates répondent différemment aux rétrovirus. Cette spécificité entre espèces doit résider dans les différences génétiques entre primates. Dans cette étude nous avons appliqué des techniques d'évolution et de génétique comparative pour caractériser le modèle évolutif des déterminants cellulaires impliqués dans la pathogenèse du VIH- 1. L'étude de l'évolution des gènes, codant pour des protéines impliquées dans la restriction ou la pathogenèse du VIH-1, aidera à la compréhension des bases génétiques ayant récemment rendu l'homme susceptible. Pour les analyses de génétique comparative, nous avons constitué une collection d'ADN et d'ARN de primates dans le but d'obtenir des nouvelles séquences de gènes orthologues. Récemment deux nouveaux génomes complets ont été publiés (l'orang-outan du Bornéo et Marmoset commun) en plus des trois génomes déjà disponibles (humain, chimpanzé, macaque rhésus). Ceci a permis d'améliorer considérablement l'étendue de l'analyse. Pour détecter l'adaptation moléculaire nous avons analysé les séquences à l'aide de méthodes phylogénétiques et statistiques. Nous avons identifié différentes pressions de sélection agissant sur les protéines impliquées dans la pathogenèse du VIH-1. Des protéines avec des propriétés de restriction du VIH-1 dans les primates non-humains présentent un taux particulièrement haut de remplacement d'acides aminés (sélection positive). En particulier dans les régions d'interaction avec les protéines virales. Ces régions incluent des acides aminés clé pour l'activité de restriction. Les protéines appartenant à l'immunité inné présentent un modèle d'évolution de conservation (sélection purifiante) mais avec des traces de "relaxation" comparé au profil général de sélection purifiante du génome des primates. Une analyse à grande échelle a permis de classifier les modèles de pression évolutive selon leur fonction moléculaire, processus biologique et distribution cellulaire. Les données générées par les différentes analyses ont permis la reconstruction ancestrale de TRIM5a, un puissant facteur antiretroviral. Le TRIM5a ressuscité, correspondant à l'ancêtre commun entre les grands singes et les groupe des catarrhiniens, est efficace contre le VIH-1 moderne. Les TRIM5a ressuscités plus récents, correspondant aux ancêtres des grands singes, sont plus efficaces contre d'autres rétrovirus. Ainsi, trouver un compromis dans la capacité de restreindre différents rétrovirus, l'homme aurait été exposé au VIH-1 à une période où TRIM5a manquait d'activité de restriction spécifique contre celui-ci. L'application de techniques d'évolution et de génétique comparative devraient être considérées pour l'évaluation systématique de protéines impliquées dans la pathogenèse virale, ainsi que pour guider des études biologiques et fonctionnelles

Signal search analysis server.

Signal search analysis is a general method to discover and characterize sequence motifs that are positionally correlated with a functional site (e.g. a transcription or translation start site). The method has played an instrumental role in the analysis of eukaryotic promoter elements. The signal search analysis server provides access to four different computer programs as well as to a large number of precompiled functional site collections. The programs offered allow: (i) the identification of non-random sequence regions under evolutionary constraint; (ii) the detection of consensus sequence-based motifs that are over- or under-represented at a particular distance from a functional site; (iii) the analysis of the positional distribution of a consensus sequence- or weight matrix-based sequence motif around a functional site; and (iv) the optimization of a weight matrix description of a locally over-represented sequence motif. These programs can be accessed at: http://www.isrec.isb-sib.ch/ssa/.

One hundred base pairs of 5' flanking sequence of a vaccinia virus late gene are sufficient to temporally regulate late transcription.

A vaccinia virus late gene coding for a major structural polypeptide of 11 kDa was sequenced. Although the 5' flanking gene region is very A+T rich, it shows little homology either to the corresponding region of vaccinia early genes or to consensus sequences characteristic of most eukaryotic genes. Three DNA fragments (100, 200, and 500 base pairs, respectively), derived from the flanking region and including the late gene mRNA start site, were inserted into the coding sequence of the vaccinia virus thymidine kinase (TK) early gene by homologous in vivo recombination. Recombinants were selected on the basis of their TK- phenotype. Cells were infected with the recombinant viruses and RNA was isolated at 1-hr intervals. Transcripts initiating either from the TK early promoter, or from the late gene promoter at its authentic position, or from the translocated late gene promoters within the early gene were detected by nuclease S1 mapping. Early after infection, only transcripts from the TK early promoter were detected. Later in infection, however, transcripts were also initiated from the translocated late promoters. This RNA appeared at the same time and in similar quantities as the RNA from the late promoter at its authentic position. No quantitative differences in promoter efficiency between the 100-, 200-, and 500-base-pair insertions were observed. We conclude that all necessary signals for correct regulation of late-gene expression reside within only 100 base pairs of 5' flanking sequence.