Transcriptomic basis for an antiserum against Micrurus corallinus (coral snake) venom

Background: Micrurus corallinus (coral snake) is a tropical forest snake belonging to the family Elapidae. Its venom shows a high neurotoxicity associated with pre- and post-synaptic toxins, causing diaphragm paralysis, which may result in death. In spite of a relatively small incidence of accidents, serum therapy is crucial for those bitten. However, the adequate production of antiserum is hampered by the difficulty in obtaining sufficient amounts of venom from a small snake with demanding breeding conditions. In order to elucidate the molecular basis of this venom and to uncover possible immunogens for an antiserum, we generated expressed sequences tags (ESTs) from its venom glands and analyzed the transcriptomic profile. In addition, their immunogenicity was tested using DNA immunization. Results: A total of 1438 ESTs were generated and grouped into 611 clusters. Toxin transcripts represented 46% of the total ESTs. The two main toxin classes consisted of three-finger toxins (3FTx) (24%) and phospholipases A(2) (PLA(2)s) (15%). However, 8 other classes of toxins were present, including C-type lectins, natriuretic peptide precursors and even high-molecular mass components such as metalloproteases and L-amino acid oxidases. Each class included an assortment of isoforms, some showing evidence of alternative splicing and domain deletions. Five antigenic candidates were selected (four 3FTx and one PLA(2)) and used for a preliminary study of DNA immunization. The immunological response showed that the sera from the immunized animals were able to recognize the recombinant antigens. Conclusion: Besides an improvement in our knowledge of the composition of coral snake venoms, which are very poorly known when compared to Old World elapids, the expression profile suggests abundant and diversified components that may be used in future antiserum formulation. As recombinant production of venom antigens frequently fails due to complex disulfide arrangements, DNA immunization may be a viable alternative. In fact, the selected candidates provided an initial evidence of the feasibility of this approach, which is less costly and not dependent on the availability of the venom.

Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

Normal and defective mariner-like elements in Rhynchosciara species (Sciaridae, Diptera)

Mariner-like elements are widely present in diverse organisms. These elements constitute a large fraction of the eukaryotic genome; they transpose by a ""cut-and-paste"" mechanism with their own transposase protein. We found two groups of mobile elements in the genus Rhynchosciara. PCR using primers designed from R. americana transposons (Ramar1 and Ramar2) were the starting point for this comparative study. Genomic DNA templates of four species: R. hollaenderi, R. millerii, R. baschanti, and Rhynchosciara sp were used and genomic sequences were amplified, sequenced and the molecular structures of the elements characterized as being putative mariner-like elements. The first group included the putative full-length elements. The second group was composed of defective mariner elements that contain a deletion overlapping most of the internal region of the transposase open reading frame. They were named Rmar1 (type 1) and Rmar2 (type 2), respectively. Many conserved amino acid blocks were identified, as well as a specific D,D(34) D signature motif that was defective in some elements. Based on predicted transposase sequences, these elements encode truncated proteins and are phylogenetically very close to mariner-like elements of the mauritiana subfamily. The inverted terminal repeat sequences that flanked the mariner-like elements are responsible for their mobility. These inverted terminal repeat sequences were identified by inverse PCR.

Genes up- and down-regulated by dermcidin in breast cancer: a microarray analysis

Dermcidin (DCD) is a human gene mapped to chromosome 12q13 region, which is co-amplified with multiple oncogenes with a well-established role in the growth, survival and progression of breast cancers. Here, we present a summary of a DNA microarray-based study that identified the genes that are up- and down-regulated in a human MDA-361 pLKO control clone and three clones expressing short hairpin RNA against three different regions of DCD mRNA. A list of 235 genes was differentially expressed among independent clones (> 3-fold change and P < 0.005). The gene expression of 208 was reduced and of 27 was increased in the three DCD-RNAi clones compared to pLKO control clone. The expression of 77 genes (37%) encoding for enzymes involved in amino acid metabolism, glucose metabolism and oxidoreductase activity and several genes required for cell survival and DNA repair were decreased. The expression of EGFR/ErbB-1 gene, an important predictor of outcome in breast cancer, was reduced together with the genes for betacellulin and amphiregulin, two known ligands of EGFR/ErbB receptors. Many of the 27 genes up-regulated by DCD-RNAi expression have not yet been fully characterized; among those with known function, we identified the calcium-calmodulin-dependent protein kinase-II delta and calcineurin A alpha. We compared 132 up-regulated and 12 down-regulated genes in our dataset with those genes up- and down-regulated by inhibitors targeting various signaling pathway components. The analysis showed that the genes in the DCD pathway are aligned with those functionally influenced by the drugs sirolimus, LY-294002 and wortmannin. Therefore, DCD may exert its function by activating the PI3K/AKT/mTOR signaling pathway. Together, these bioinformatic approaches suggest the involvement of DCD in the regulation of genes for breast cancer cell metabolism, proliferation and survival.

Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus)microplus in response to infection with Anaplasma marginale

Background: Bovine anaplasmosis, caused by the rickettsial tick-borne pathogen Anaplasma marginale (Rickettsiales: Anaplasmataceae), is vectored by Rhipicephalus (Boophilus) microplus in many tropical and subtropical regions of the world. A. marginale undergoes a complex developmental cycle in ticks which results in infection of salivary glands from where the pathogen is transmitted to cattle. In previous studies, we reported modification of gene expression in Dermacentor variabilis and cultured Ixodes scapularis tick cells in response to infection with A. marginale. In these studies, we extended these findings by use of a functional genomics approach to identify genes differentially expressed in R. microplus male salivary glands in response to A. marginale infection. Additionally, a R. microplus-derived cell line, BME26, was used for the first time to also study tick cell gene expression in response to A. marginale infection. Results: Suppression subtractive hybridization libraries were constructed from infected and uninfected ticks and used to identify genes differentially expressed in male R. microplus salivary glands infected with A. marginale. A total of 279 ESTs were identified as candidate differentially expressed genes. Of these, five genes encoding for putative histamine-binding protein (22Hbp), von Willebrand factor (94Will), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz) and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were confirmed by real-time RT-PCR to be down-regulated in tick salivary glands infected with A. marginale. The impact of selected tick genes on A. marginale infections in tick salivary glands and BME26 cells was characterized by RNA interference. Silencing of the gene encoding for putative flagelliform silk protein (100Silk) resulted in reduced A. marginale infection in both tick salivary glands and cultured BME26 cells, while silencing of the gene encoding for subolesin (4D8) significantly reduced infection only in cultured BME26 cells. The knockdown of the gene encoding for putative metallothionein (93 Meth), significantly up-regulated in infected cultured BME26 cells, resulted in higher A. marginale infection levels in tick cells. Conclusions: Characterization of differential gene expression in salivary glands of R. microplus in response to A. marginale infection expands our understanding of the molecular mechanisms at the tick-pathogen interface. Functional studies suggested that differentially expressed genes encoding for subolesin, putative von Willebrand factor and flagelliform silk protein could play a role in A. marginale infection and multiplication in ticks. These tick genes found to be functionally relevant for tick-pathogen interactions will likely be candidates for development of vaccines designed for control of both ticks and tick-borne pathogens.

Leishmania-Specific Surface Antigens Show Sub-Genus Sequence Variation and Immune Recognition

Background: A family of hydrophilic acylated surface (HASP) proteins, containing extensive and variant amino acid repeats, is expressed at the plasma membrane in infective extracellular (metacyclic) and intracellular (amastigote) stages of Old World Leishmania species. While HASPs are antigenic in the host and can induce protective immune responses, the biological functions of these Leishmania-specific proteins remain unresolved. Previous genome analysis has suggested that parasites of the sub-genus Leishmania (Viannia) have lost HASP genes from their genomes. Methods/Principal Findings: We have used molecular and cellular methods to analyse HASP expression in New World Leishmania mexicana complex species and show that, unlike in L. major, these proteins are expressed predominantly following differentiation into amastigotes within macrophages. Further genome analysis has revealed that the L. (Viannia) species, L. (V.) braziliensis, does express HASP-like proteins of low amino acid similarity but with similar biochemical characteristics, from genes present on a region of chromosome 23 that is syntenic with the HASP/SHERP locus in Old World Leishmania species and the L. (L.) mexicana complex. A related gene is also present in Leptomonas seymouri and this may represent the ancestral copy of these Leishmania-genus specific sequences. The L. braziliensis HASP-like proteins (named the orthologous (o) HASPs) are predominantly expressed on the plasma membrane in amastigotes and are recognised by immune sera taken from 4 out of 6 leishmaniasis patients tested in an endemic region of Brazil. Analysis of the repetitive domains of the oHASPs has shown considerable genetic variation in parasite isolates taken from the same patients, suggesting that antigenic change may play a role in immune recognition of this protein family. Conclusions/Significance: These findings confirm that antigenic hydrophilic acylated proteins are expressed from genes in the same chromosomal region in species across the genus Leishmania. These proteins are surface-exposed on amastigotes (although L. (L.) major parasites also express HASPB on the metacyclic plasma membrane). The central repetitive domains of the HASPs are highly variant in their amino acid sequences, both within and between species, consistent with a role in immune recognition in the host.

Single-nucleotide polymorphism, linkage disequilibrium and geographic structure in the malaria parasite Plasmodium vivax: prospects for genome-wide association studies

Background: The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized. Results: We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for similar to 40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance. Conclusion: These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.

Psychogenetics of Turner syndrome: an investigation of 28 subjects and respective controls using the Bender test and Piagetian scales

Piagetian scales and the Bender visual motor gestalt test (BT) were applied to 28 subjects with universal 45, X Turner syndrome (TS), and their respective controls, in order to investigate their cognitive performance. Dermatoglyphics were also analyzed to obtain clues concerning embryological changes that may have appeared during development of the nervous system and could be associated with cognitive performance of TS patients. Dermatoglyphic pattern distribution was similar to that reported in previous studies of TS individuals: ulnar loops in the digital patterns and finger ridge, a-b, and A'-d counts were more frequent, while arch and whorl patterns were less frequent compared to controls. However, we did not find higher frequencies of hypothenar pattern, maximum atd angle, and ulnarity index in our TS subjects, unlike other investigations. Furthermore, we found significant differences between TS and control T line index values. The BT scores were also lower in probands, as has been previously reported, revealing a neurocognitive deficit of visual motor perception in TS individuals, which could be due to an absence of, or deficiency in, cerebral hemispheric lateralization. However, TS subjects seemed to improve their performance on BT with age. Cognitive performance of the TS subjects was not significantly different from that of controls, confirming a previous study in which TS performance was found to be similar to that of the normal Brazilian population. There were significant correlations between BT scores and Piagetian scale levels with dermatoglyphic parameters. This association could be explained by changes in the common ectodermal origin of the epidermis and the central nervous system. TS subjects seem to succeed in compensating their spatial impairments in adapting their cognitive and social contacts. We concluded that genetic counseling should consider cognitive and psychosocial difficulties presented by TS subjects, providing appropriate treatment and orientation for them and their families.

Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp aurantifolii

Background: Citrus canker is a disease that has severe economic impact on the citrus industry worldwide. There are three types of canker, called A, B, and C. The three types have different phenotypes and affect different citrus species. The causative agent for type A is Xanthomonas citri subsp. citri, whose genome sequence was made available in 2002. Xanthomonas fuscans subsp. aurantifolii strain B causes canker B and Xanthomonas fuscans subsp. aurantifolii strain C causes canker C. Results: We have sequenced the genomes of strains B and C to draft status. We have compared their genomic content to X. citri subsp. citri and to other Xanthomonas genomes, with special emphasis on type III secreted effector repertoires. In addition to pthA, already known to be present in all three citrus canker strains, two additional effector genes, xopE3 and xopAI, are also present in all three strains and are both located on the same putative genomic island. These two effector genes, along with one other effector-like gene in the same region, are thus good candidates for being pathogenicity factors on citrus. Numerous gene content differences also exist between the three cankers strains, which can be correlated with their different virulence and host range. Particular attention was placed on the analysis of genes involved in biofilm formation and quorum sensing, type IV secretion, flagellum synthesis and motility, lipopolysacharide synthesis, and on the gene xacPNP, which codes for a natriuretic protein. Conclusion: We have uncovered numerous commonalities and differences in gene content between the genomes of the pathogenic agents causing citrus canker A, B, and C and other Xanthomonas genomes. Molecular genetics can now be employed to determine the role of these genes in plant-microbe interactions. The gained knowledge will be instrumental for improving citrus canker control.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

The role of exon shuffling in shaping protein-protein interaction networks

Background: Physical protein-protein interaction (PPI) is a critical phenomenon for the function of most proteins in living organisms and a significant fraction of PPIs are the result of domain-domain interactions. Exon shuffling, intron-mediated recombination of exons from existing genes, is known to have been a major mechanism of domain shuffling in metazoans. Thus, we hypothesized that exon shuffling could have a significant influence in shaping the topology of PPI networks. Results: We tested our hypothesis by compiling exon shuffling and PPI data from six eukaryotic species: Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Cryptococcus neoformans and Arabidopsis thaliana. For all four metazoan species, genes enriched in exon shuffling events presented on average higher vertex degree (number of interacting partners) in PPI networks. Furthermore, we verified that a set of protein domains that are simultaneously promiscuous (known to interact to multiple types of other domains), self-interacting (able to interact with another copy of themselves) and abundant in the genomes presents a stronger signal for exon shuffling. Conclusions: Exon shuffling appears to have been a recurrent mechanism for the emergence of new PPIs along metazoan evolution. In metazoan genomes, exon shuffling also promoted the expansion of some protein domains. We speculate that their promiscuous and self-interacting properties may have been decisive for that expansion.

Validation of dot blot hybridization and denaturing high performance liquid chromatography as reliable methods for TP53 codon 72 genotyping in molecular epidemiologic studies

Background: Mutations in TP53 are common events during carcinogenesis. In addition to gene mutations, several reports have focused on TP53 polymorphisms as risk factors for malignant disease. Many studies have highlighted that the status of the TP53 codon 72 polymorphism could influence cancer susceptibility. However, the results have been inconsistent and various methodological features can contribute to departures from Hardy-Weinberg equilibrium, a condition that may influence the disease risk estimates. The most widely accepted method of detecting genotyping error is to confirm genotypes by sequencing and/or via a separate method. Results: We developed two new genotyping methods for TP53 codon 72 polymorphism detection: Denaturing High Performance Liquid Chromatography (DHPLC) and Dot Blot hybridization. These methods were compared with Restriction Fragment Length Polymorphism (RFLP) using two different restriction enzymes. We observed high agreement among all methodologies assayed. Dot-blot hybridization and DHPLC results were more highly concordant with each other than when either of these methods was compared with RFLP. Conclusions: Although variations may occur, our results indicate that DHPLC and Dot Blot hybridization can be used as reliable screening methods for TP53 codon 72 polymorphism detection, especially in molecular epidemiologic studies, where high throughput methodologies are required.

Evaluating Theobroma grandiflorum for comparative genomic studies with Theobroma cacao

The seeds of Theobroma cacao (cacao) are the source of cocoa, the raw material for the multi-billion dollar chocolate industry. Cacao`s two most important traits are its unique seed storage triglyceride (cocoa butter) and the flavor of its fermented beans (chocolate). The genome of T. cacao is being sequenced, and to expand the utility of the genome sequence to the improvement of cacao, we are evaluating Theobroma grandiflorum, the closest economically important species of Theobroma for its potential use in a comparative genomic study. T. grandiflorum differs from cacao in important agronomic traits such as flavor of the fermented beans, disease resistance to witches` broom and abscission of mature fruits. By comparing genomic sequences and analyzing viable inter-specific hybrids, we hope to identify the key genes that regulate cacao`s most important traits. We have investigated the utility in T. grandiflorum of three types of markers (microsatellite markers, single-strand conformational polymorphism markers and single nucleotide polymorphism (SNP) markers) developed in cacao. Through sequencing of amplicons of 12 diverse individuals of both cacao and T. grandiflorum, we have identified new intra- and inter-specific SNPs. Two markers which had no overlap of alleles between the species were used to genotype putative inter-specific hybrid seedlings. Sequence conservation was significant and species-specific differences numerous enough to suggest that comparative genomics of T. grandiflorum and T. cacao will be useful in elucidating the genetic differences that lead to a variety of important agronomic trait differences.

The mitochondrial genome of the phytopathogenic basidiomycete Moniliophthora perniciosa is 109 kb in size and contains a stable integrated plasmid

We present here the sequence of the mitochondrial genome of the basidiomycete phytopathogenic hemibiotrophic fungus Moniliophthora perniciosa, causal agent of the Witches` Broom Disease in Theobroma cacao. The DNA is a circular molecule of 109103 base pairs, with 31.9 % GC, and is the largest sequenced so far. This size is due essentially to the presence of numerous non-conserved hypothetical ORFs. It contains the 14 genes coding for proteins involved in the oxidative phosphorylation, the two rRNA genes, one ORF coding for a ribosomal protein (rps3), and a set of 26 tRNA genes that recognize codons for all amino acids. Seven homing endonucleases are located inside introns. Except atp8, all conserved known genes are in the same orientation. Phylogenetic analysis based on the cox genes agrees with the commonly accepted fungal taxonomy. An uncommon feature of this mitochondrial genome is the presence of a region that contains a set of four, relatively small, nested, inverted repeats enclosing two genes coding for polymerases with an invertron-type structure and three conserved hypothetical genes interpreted as the stable integration of a mitochondrial linear plasmid. The integration of this plasmid seems to be a recent evolutionary event that could have implications in fungal biology. This sequence is available under GenBank accession number AY376688. (c) 2008 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.