Large multi-gene phylogenetic trees of the grasses (Poaceae): progress towards complete tribal and generic level sampling.

In this paper we included a very broad representation of grass family diversity (84% of tribes and 42% of genera). Phylogenetic inference was based on three plastid DNA regions rbcL, matK and trnL-F, using maximum parsimony and Bayesian methods. Our results resolved most of the subfamily relationships within the major clades (BEP and PACCMAD), which had previously been unclear, such as, among others the: (i) BEP and PACCMAD sister relationship, (ii) composition of clades and the sister-relationship of Ehrhartoideae and Bambusoideae + Pooideae, (iii) paraphyly of tribe Bambuseae, (iv) position of Gynerium as sister to Panicoideae, (v) phylogenetic position of Micrairoideae. With the presence of a relatively large amount of missing data, we were able to increase taxon sampling substantially in our analyses from 107 to 295 taxa. However, bootstrap support and to a lesser extent Bayesian inference posterior probabilities were generally lower in analyses involving missing data than those not including them. We produced a fully resolved phylogenetic summary tree for the grass family at subfamily level and indicated the most likely relationships of all included tribes in our analysis.

Update of the Gene Discovery Program in Schistosoma mansoni with the Expressed Sequence Tag Approach

Continuing the Schistosoma mansoni Genome Project 363 new templates were sequenced generating 205 more ESTs corresponding to 91 genes. Seventy four of these genes (81%) had not previously been described in S. mansoni. Among the newly discovered genes there are several of significant biological interest such as synaptophysin, NIFs-like and rho-GDP dissociation inhibitor

Towards the Physical Map of the Trypanosoma cruzi Nuclear Genome: Construction of YAC and BAC Libraries of the Reference Clone T. cruzi CL-Brener

Strategies to construct the physical map of the Trypanosoma cruzi nuclear genome have to capitalize on three main advantages of the parasite genome, namely (a) its small size, (b) the fact that all chromosomes can be defined, and many of them can be isolated by pulse field gel electrophoresis, and (c) the fact that simple Southern blots of electrophoretic karyotypes can be used to map sequence tagged sites and expressed sequence tags to chromosomal bands. A major drawback to cope with is the complexity of T. cruzi genetics, that hinders the construction of a comprehensive genetic map. As a first step towards physical mapping, we report the construction and partial characterization of a T. cruzi CL-Brener genomic library in yeast artificial chromosomes (YACs) that consists of 2,770 individual YACs with a mean insert size of 365 kb encompassing around 10 genomic equivalents. Two libraries in bacterial artificial chromosomes (BACs) have been constructed, BACI and BACII. Both libraries represent about three genome equivalents. A third BAC library (BAC III) is being constructed. YACs and BACs are invaluable tools for physical mapping. More generally, they have to be considered as a common resource for research in Chagas disease

Identification of Transcribed Sequences (ESTs) in the Trypanosoma cruzi Genome Project

Random single pass sequencing of cDNA fragments, also known as generation of Expressed Sequence Tags (ESTs), has been highly successful in the study of the gene content of higher organisms, and forms an integral part of most genome projects, with the objective to identify new genes and targets for disease control and prevention and to generate mapping probes. In the Trypanosoma cruzi genome project, EST sequencing has also been a starting point, and here we report data on the first 797 sequences obtained, partly from a CL Brener epimastigote non-normalized library, partly on a normalized library. Only around 30% of the sequences obtained showed similarity with Genbank and dbEST databases, half of which with sequences already reported for T. cruzi.

Characterization of constitutive and putative differentially expressed mRNAs by means of expressed sequence tags, differential display reverse transcriptase-PCR and randomly amplified polymorphic DNA-PCR from the sand fly vector Lutzomyia longipalpis

Molecular studies of insect disease vectors are of paramount importance for understanding parasite-vector relationship. Advances in this area have led to important findings regarding changes in vectors' physiology upon blood feeding and parasite infection. Mechanisms for interfering with the vectorial capacity of insects responsible for the transmission of diseases such as malaria, Chagas disease and dengue fever are being devised with the ultimate goal of developing transgenic insects. A primary necessity for this goal is information on gene expression and control in the target insect. Our group is investigating molecular aspects of the interaction between Leishmania parasites and Lutzomyia sand flies. As an initial step in our studies we have used random sequencing of cDNA clones from two expression libraries made from head/thorax and abdomen of sugar fed L. longipalpis for the identification of expressed sequence tags (EST). We applied differential display reverse transcriptase-PCR and randomly amplified polymorphic DNA-PCR to characterize differentially expressed mRNA from sugar and blood fed insects, and, in one case, from a L. (V.) braziliensis-infected L. longipalpis. We identified 37 cDNAs that have shown homology to known sequences from GeneBank. Of these, 32 cDNAs code for constitutive proteins such as zinc finger protein, glutamine synthetase, G binding protein, ubiquitin conjugating enzyme. Three are putative differentially expressed cDNAs from blood fed and Leishmania-infected midgut, a chitinase, a V-ATPase and a MAP kinase. Finally, two sequences are homologous to Drosophila melanogaster gene products recently discovered through the Drosophila genome initiative.

Genome comparison of progressively drug resistant Plasmodium falciparum lines derived from drug sensitive clone

Chloroquine has been the mainstay of malaria chemotherapy for the past five decades, but resistance is now widespread. Pyrimethamine or proguanil form an important component of some alternate drug combinations being used for treatment of uncomplicated Plasmodium falciparum infections in areas of chloroquine resistance. Both pyrimethamine and proguanil are dihydrofolate reductase (DHFR) inhibitors, the proguanil acting primarily through its major metabolite cycloguanil. Resistance to these drugs arises due to specific point mutations in the dhfr gene. Cross resistance between cycloguanil and pyrimethamine is not absolute. It is, therefore, important to investigate mutation rates in P. falciparum for pyrimethamine and proguanil so that DHFR inhibitor with less mutation rate is favored in drug combinations. Hence, we have compared mutation rates in P. falciparum genome for pyrimethamine and cycloguanil. Using erythrocytic stages of P. falciparum cultures, progressively drug resistant lines were selected in vitro and comparing their RFLP profile with a repeat sequence. Our finding suggests that pyrimethamine has higher mutation rate compared to cycloguanil. It enhances the degree of genomic polymorphism leading to diversity of natural parasite population which in turn is predisposes the parasites for faster selection of resistance to some other antimalarial drugs.

Vertebrate conserved non coding DNA regions have a high persistence length and a short persistence time.

BACKGROUND: The comparison of complete genomes has revealed surprisingly large numbers of conserved non-protein-coding (CNC) DNA regions. However, the biological function of CNC remains elusive. CNC differ in two aspects from conserved protein-coding regions. They are not conserved across phylum boundaries, and they do not contain readily detectable sub-domains. Here we characterize the persistence length and time of CNC and conserved protein-coding regions in the vertebrate and insect lineages. RESULTS: The persistence length is the length of a genome region over which a certain level of sequence identity is consistently maintained. The persistence time is the evolutionary period during which a conserved region evolves under the same selective constraints.Our main findings are: (i) Insect genomes contain 1.60 times less conserved information than vertebrates; (ii) Vertebrate CNC have a higher persistence length than conserved coding regions or insect CNC; (iii) CNC have shorter persistence times as compared to conserved coding regions in both lineages. CONCLUSION: Higher persistence length of vertebrate CNC indicates that the conserved information in vertebrates and insects is organized in functional elements of different lengths. These findings might be related to the higher morphological complexity of vertebrates and give clues about the structure of active CNC elements.Shorter persistence time might explain the previously puzzling observations of highly conserved CNC within each phylum, and of a lack of conservation between phyla. It suggests that CNC divergence might be a key factor in vertebrate evolution. Further evolutionary studies will help to relate individual CNC to specific developmental processes.

Sequencing of simple sequence repeat anchored polymerase chain reaction amplification products of Biomphalaria glabrata

Simple sequence repeat anchored polymerase chain reaction amplification (SSR-PCR) is a genetic typing technique based on primers anchored at the 5' or 3' ends of microsatellites, at high primer annealing temperatures. This technique has already been used in studies of genetic variability of several organisms, using different primer designs. In order to conduct a detailed study of the SSR-PCR genomic targets, we cloned and sequenced 20 unique amplification products of two commonly used primers, CAA(CT)6 and (CA)8RY, using Biomphalaria glabrata genomic DNA as template. The sequences obtained were novel B. glabrata genomic sequences. It was observed that 15 clones contained microsatellites between priming sites. Out of 40 clones, seven contained complex sequence repetitions. One of the repeats that appeared in six of the amplified fragments generated a single band in Southern analysis, indicating that the sequence was not widespread in the genome. Most of the annealing sites for the CAA(CT)6 primer contained only the six repeats found within the primer sequence. In conclusion, SSR-PCR is a useful genotyping technique. However, the premise of the SSR-PCR technique, verified with the CAA(CT)6 primer, could not be supported since the amplification products did not result necessarily from microsatellite loci amplification.

A novel reiterated family of transcribed oligo(A)-terminated, interspersed DNA elements in the genome of Trypanosoma cruzi

We report the molecular characterization of a novel reiterated family of transcribed oligo(A)-terminated, interspersed DNA elements in the genome of Trypanosoma cruzi. Steady-state level of transcripts of this sequence family appeared to be developmentally regulated, since only in the replicative forms the parasite showed expression of related sequences with a major band around 3 kb. The presence of frame shifts or premature stop codons predicts that transcripts are not translated. The sequence family also contains truncated forms of retrotransposons elements that may become potential hot spots for retroelement insertion. Sequences homologous to this family are interspersed at many chromosomes including the subtelomeric regions.

Sequence-Based Hepatitis B Virus Antiviral Resistance Testing in Switzerland

BACKGROUND: A growing number of patients with chronic hepatitis B is being treated for extended periods with nucleoside and/or nucleotide analogs. In this context, antiviral resistance represents an increasingly common and complex issue. METHODS: Mutations in the hepatitis B virus (HBV) reverse transcriptase (rt) gene and viral genotypes were determined by direct sequencing of PCR products and alignment with reference sequences deposited in GenBank. RESULTS: Plasma samples from 60 patients with chronic hepatitis B were analyzed since March 2009. The predominant mutation pattern identified in patients with virological breakthrough was rtM204V/I ± different compensatory mutations, conferring resistance to L-nucleosides (lamivudine, telbivudine, emtricitabine) and predisposing to entecavir resistance (n = 18). Complex mutation patterns with a potential for multidrug resistance were identified in 2 patients. Selection of a fully entecavir resistant strain was observed in a patient exposed to lamivudine alone. Novel mutations were identified in 1 patient. Wild-type HBV was identified in 9 patients with suspected virological breakthrough, raising concerns about treatment adherence. No preexisting resistance mutations were identified in treatment-naïve patients (n = 13). Viral genome amplification and sequencing failed in 16 patients, of which only 2 had a documented HBV DNA > 1000 IU/ml. HBV genotypes were D in 28, A in 6, B in 4, C in 3 and E in 3 patients. Results will be updated in August 2010 and therapeutic implications discussed. CONCLUSIONS: With expanding treatment options and a growing number of patients exposed to nucleoside and/or nucleotide analogs, sequence-based HBV antiviral resistance testing is expected to become a cornerstone in the management of chronic hepatitis B.

Developmental and Anatomical Constraints on Vertebrate Genome Evolution

Pendant ma thèse de doctorat, j'ai utilisé des espèces modèles, comme la souris et le poisson-zèbre, pour étudier les facteurs qui affectent l'évolution des gènes et leur expression. Plus précisément, j'ai montré que l'anatomie et le développement sont des facteurs clés à prendre en compte, car ils influencent la vitesse d'évolution de la séquence des gènes, l'impact sur eux de mutations (i.e. la délétion du gène est-elle létale ?), et leur tendance à se dupliquer. Où et quand il est exprimé impose à un gène certaines contraintes ou au contraire lui donne des opportunités d'évoluer. J'ai pu comparer ces tendances aux modèles classiques d'évolution de la morphologie, que l'on pensait auparavant refléter directement les contraintes s'appliquant sur le génome. Nous avons montré que les contraintes entre ces deux niveaux d'organisation ne peuvent pas être transférées simplement : il n'y a pas de lien direct entre la conservation du génotype et celle de phénotypes comme la morphologie. Ce travail a été possible grâce au développement d'outils bioinformatiques. Notamment, j'ai travaillé sur le développement de la base de données Bgee, qui a pour but de comparer l'expression des gènes entre différentes espèces de manière automatique et à large échelle. Cela implique une formalisation de l'anatomie, du développement et de concepts liés à l'homologie grâce à l'utilisation d'ontologies. Une intégration cohérente de données d'expression hétérogènes (puces à ADN, marqueurs de séquence exprimée, hybridations in situ) a aussi été nécessaire. Cette base de données est mise à jour régulièrement et disponible librement. Elle devrait contribuer à étendre les possibilités de comparaison de l'expression des gènes entre espèces pour des études d'évo-devo (évolution du développement) et de génomique. During my PhD, I used model species of vertebrates, such as mouse and zebrafish, to study factors affecting the evolution of genes and their expression. More precisely I have shown that anatomy and development are key factors to take into account, influencing the rate of gene sequence evolution, the impact of mutations (i.e. is the deletion of a gene lethal?), and the propensity of a gene to duplicate. Where and when genes are expressed imposes constraints, or on the contrary leaves them some opportunity to evolve. We analyzed these patterns in relation to classical models of morphological evolution in vertebrates, which were previously thought to directly reflect constraints on the genomes. We showed that the patterns of evolution at these two levels of organization do not translate smoothly: there is no direct link between the conservation of genotype and phenotypes such as morphology. This work was made possible by the development of bioinformatics tools. Notably, I worked on the development of the database Bgee, which aims at comparing gene expression between different species in an automated and large-scale way. This involves the formalization of anatomy, development, and concepts related to homology, through the use of ontologies. A coherent integration of heterogeneous expression data (microarray, expressed sequence tags, in situ hybridizations) is also required. This database is regularly updated and freely available. It should contribute to extend the possibilities for comparison of gene expression between species in evo-devo and genomics studies.

The amphioxus genome and the evolution of the chordate karyotype.

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

Trypanosoma cruzi ribosomal protein S4: characterization of its coding locus, analysis of transcripts, and antigenicity of the protein

Two allelic genomic fragments containing ribosomal protein S4 encoding genes (rpS4) from Trypanosoma cruzi (CL-Brener strain) were isolated and characterized. One allele comprises two complete tandem repeats of a sequence encoding an rpS4 gene. In the other, only one rpS4 gene is found. Sequence comparison to the accessed data in the genome project database reveals that our two-copy allele corresponds to a variant haplotype. However, the deduced aminoacid sequence of all the gene copies is identical. The rpS4 transcripts processing sites were determined by comparison of genomic sequences with published cDNA data. The obtained sequence data demonstrates that rpS4 genes are expressed in epimastigotes, amastigotes, and trypomastigotes. A recombinant version of rpS4 was found to be an antigenic: it was recognized by 62.5% of the individuals with positive serology for T. cruzi and by 93.3% of patients with proven chronic chagasic disease.

Analysis and chromosomal mapping of Leishmania (Leishmania) amazonensis amastigote expressed sequence tags

The characterization of expressed sequence tags (ESTs) generated from a cDNA library of Leishmania (Leishmania) amazonensis amastigotes is described. The sequencing of 93 clones generated new L. (L.) amazonensis ESTs from which 32% are not related to any other sequences in database and 68% presented significant similarities to known genes. The chromosome localization of some L. (L.) amazonensis ESTs was also determined in L. (L.) amazonensis and L. (L.) major. The characterization of these ESTs is suitable for the genome physical mapping, as well as for the identification of genes encoding cysteine proteinases implicated with protective immune responses in leishmaniasis.

Expression of PROKR1 and PROKR2 in Human Enteric Neural Precursor Cells and Identification of Sequence Variants Suggest a Role in HSCR

Background. The enteric nervous system (ENS) is entirely derived from neural crest and its normal development is regulated by specific molecular pathways. Failure in complete ENS formation results in aganglionic gut conditions such as Hirschsprung's disease (HSCR). Recently, PROKR1 expression has been demonstrated in mouse enteric neural crest derived cells and Prok-1 was shown to work coordinately with GDNF in the development of the ENS. Principal Findings. In the present report, ENS progenitors were isolated and characterized from the ganglionic gut from children diagnosed with and without HSCR, and the expression of prokineticin receptors was examined. Immunocytochemical analysis of neurosphere-forming cells demonstrated that both PROKR1 and PROKR2 were present in human enteric neural crest cells. In addition, we also performed a mutational analysis of PROKR1, PROKR2, PROK1 and PROK2 genes in a cohort of HSCR patients, evaluating them for the first time as susceptibility genes for the disease. Several missense variants were detected, most of them affecting highly conserved amino acid residues of the protein and located in functional domains of both receptors, which suggests a possible deleterious effect in their biological function. Conclusions. Our results suggest that not only PROKR1, but also PROKR2 might mediate a complementary signalling to the RET/GFRα1/GDNF pathway supporting proliferation/survival and differentiation of precursor cells during ENS development. These findings, together with the detection of sequence variants in PROKR1, PROK1 and PROKR2 genes associated to HSCR and, in some cases in combination with RET or GDNF mutations, provide the first evidence to consider them as susceptibility genes for HSCR.