Optimització d'una aplicació bioinformàtica d'aliniament de seqüències executada en processadors many-core (GPUs)

Las herramientas de análisis de secuencias genómicas permiten a los biólogos identificar y entender regiones fundamentales que tienen implicación en enfermedades genéticas. Actualmente existe una necesidad de dotar al ámbito científico de herramientas de análisis eficientes. Este proyecto lleva a cabo una caracterización y análisis del rendimiento de algoritmos utilizados en la comparación de secuencias genómicas completas, y ejecutadas en arquitecturas MultiCore y ManyCore. A partir del análisis se evalúa la idoneidad de este tipo de arquitecturas para resolver el problema de comparar secuencias genómicas. Finalmente se propone una serie de modificaciones en las implementaciones de estos algoritmos con el objetivo de mejorar el rendimiento.

New insights into the long-term evolutionary history of the avian MHC

SummaryGene duplication and neofunctidnalization are important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family with a central role in vertebrates' adaptive immune system. Multigene families, which evolved in large part through duplication, represent promising systems to study the still strongly depbated relative roles of neutral and adaptive processes in the evolution of phenotypic complexity. Detailed knowledge on ecological function and a well-characterized evolutionary history place the mammals' MHC amongst ideal study systems. However mammalian MHCs usually encompass several million base pairs and hold a large number of functional and non-functional duplicate genes, which makes their study complex. Avian MHCs on the other hand are usually way more compact, but the reconstruction of. their evolutionary history has proven notoriously difficult. However, no focused attempt has been undertaken so far to study the avian MHC evolutionary history in a broad phylogenetic context and using adequate gene regions.In the present PhD, we were able to make important contributions to the understanding of the long-term evolution of the avian MHC class II Β (MHCI1B). First, we isolated and characterized MHCIIB genes in barn owl (Tyto alba?, Strigiformes, Tytonidae), a species from an avian lineage in which MHC has not been studied so far. Our results revealed that with only two functional MHCIIB genes the MHC organization of barn owl may be similar to the 'minimal essential' MHC of chicken (Gallus gallus), indicating that simple MHC organization may be ancestral to birds. Taking advantage of the sequence information from barn owl, we studied the evolution of MHCIIB genes in 13 additional species of 'typical' owls (Strigiformes, Strigidae). Phylogenetic analyses revealed that according to their function, in owls the peptide-binding region (PBR) encoding exon 2 and the non-PBR encoding exon 3 evolve by different patterns. Exon 2 exhibited an evolutionary history of positive selection and recombination, while exon 3 traced duplication history and revealed two paralogs evolving divergently from each other in owls, and in a shorebird, the great snipe {Gallinago media). The results from exon 3 were the first ever from birds to demonstrate gene orthology in species that diverged tens of millions of years ago, and strongly questioned whether the taxa studied before provided an adequate picture of avian MHC evolution. In a follow-up study, we aimed at explaining a striking pattern revealed by phylogenetic trees analyzing the owl sequences along with MHCIIB sequences from other birds: One owl paralog (termed DAB1) grouped with sequences of passerines and falcons, while the other (DAB2) grouped with wildfowl, penguins and birds of prey. This could be explained by either a duplication event preceding the evolution of these bird orders, or by convergent evolution of similar sequences in a number of orders. With extensive phylogenetic analyses we were able to show, that indeed a duplication event preceeded the major avian radiation -100 my ago, and that following this duplication, the paralogs evolved under positive selection. Furthermore, we showed that the divergently evolving amino acid residues in the MHCIIB-encoded β-chain potentially interact with the MHCI I α-chain, and that molecular coevolution of the interacting residues may have been involved in the divergent evolution of the MHCIIB paralogs.The findings of this PhD are of particular interest to the understanding of the evolutionary history of the avian MHC and, by providing essential information on long-term gene history in the avian MHC, open promising perspectives for advances in the understanding of the evolution of multigene families in general, and for avian MHC organization in particular. Amongst others I discuss the importance of including protein structure in the phylogenetic study of multigene families, and the roles of ecological versus molecular selection pressures. I conclude by providing a population genomic perspective on avian MHC, which may serve as a basis for future research to investigate the relative roles of neutral processes involving effective population size effects and of adaptation in the evolution of avian MHC diversity and organization.RésuméLa duplication de gènes et leur néo-fonctionnalisation sont des processus importants dans l'évolution de la complexité phénotypique. Ils sont impliqués dans l'apparition d'importantes nouveautés évolutives favorisant l'adaptation écologique, comme c'est le cas pour le complexe majeur d'histocompatibilité

Loss of single HLA class I allospecificities in melanoma cells due to selective genomic abbreviations.

Expression of human leucocyte antigen (HLA) Class I molecules is essential for the recognition of malignant melanoma (MM) cells by CD8(+) T lymphocytes. A complete or partial loss of HLA Class I molecules is a potent strategy for MM cells to escape from immunosurveillance. In 2 out of 55 melanoma cell cultures we identified a complete phenotypic loss of HLA allospecificities. Both patients have been treated unsuccessfully with HLA-A2 peptides. To identify the reasons underlying the loss of single HLA-A allospecificities, we searched for genomic alterations at the locus for HLA Class I alpha-chain on chromosome 6 in melanoma cell cultures established from 2 selected patients with MM in advanced stage. This deficiency was associated with alterations of HLA-A2 gene sequences as determined by polymerase chain reaction-sequence specific primers (PCR-SSP). Karyotyping revealed a chromosomal loss in Patient 1, whereas melanoma cell cultures established from Patient 2 displayed 2 copies of chromosome 6. Loss of heterozygosity (LOH) using markers located around position 6p21 was detected in both cases. By applying group-specific primer-mixes spanning the 5'-flanking region of the HLA-A2 gene locus the relevant region was amplified by PCR and subsequent sequencing allowed alignment with the known HLA Class I reference sequences. Functional assays using HLA-A2-restricted cytotoxic T-cell clones were performed in HLA-A2 deficient MM cultures and revealed a drastically reduced susceptibility to CTL lysis in HLA-A2 negative cells. We could document the occurrence of selective HLA-A2 deficiencies in cultured advanced-stage melanoma metastases and identify their molecular causes as genomic alterations within the HLA-A gene locus.

Electrocardiographic findings in Mexican chagasic subjects living in high and low endemic regions of Trypanosoma cruzi infection

In México the first human chronic chagasic case was recognized in 1940. In spite of an increasing number of cases detected since that time, Chagas disease in México has been poorly documented. In the present work we studied 617 volunteers subjects living in high and low endemic regions of Trypanosoma cruzi infection with seroprevalence of 22% and 4% respectively. Hemoculture performed in those seropositive subjects failed to demonstrate circulating parasites, however polymerase chain reaction identified up to 60% of them as positives. A higher level of anti-T. cruzi antibodies was observed in seropositive residents in high endemic region, in spite of similar parasite persistence (p < 0.05). On standard 12 leads electrocardiogram (ECG) 20% to 22% seropositive individuals from either region showed right bundle branch block or ventricular extrasystoles which were more prevalent in seropositive than in seronegative individuals (p < 0.05). In conclusion, the frequency or type of ECG abnormality was influenced by serologic status but not by endemicity or parasite persistence. Furthermore, Mexican indeterminate patients have a similar ECG pattern to those reported in South America.

Presence of Wolbachia endosymbionts in microfilariae of Wuchereria bancrofti (Spirurida: Onchocercidae) from different geographical regions in India

In view of the recent discovery of rickettsial endosymbionts, Wolbachia in lymphatic filarial parasites, Wuchereria bancrofti and Brugia malayi and subsequently of their vital role in the survival and development of the latter, antibiotics such as tetracycline are being suggested for the treatment of lymphatic filariasis, by way of eliminating the endosymbiont. But, it is essential to assess their presence in parasites from areas endemic for lymphatic filariasis before such a new control tool is employed. In the present communication, we report the detection of Wolbachia endosymbionts in microfilariae of W. bancrofti parasites collected from geographically distant locations of India, such as Pondicherry (Union Territory), Calicut (Kerala), Jagadalpur (Madhya Pradesh), Thirukoilur (TamilNadu), Chinnanergunam (TamilNadu), Rajahmundry (Andhra Pradesh), and Varanasi (Uttar Pradesh), using Wolbachia specific 16S rDNA polymerase chain reaction.

DIP-STR: Highly Sensitive Markers for the Analysis of Unbalanced Genomic Mixtures.

Samples containing highly unbalanced DNA mixtures from two individuals commonly occur both in forensic mixed stains and in peripheral blood DNA microchimerism induced by pregnancy or following organ transplant. Because of PCR amplification bias, the genetic identification of a DNA that contributes trace amounts to a mixed sample represents a tremendous challenge. This means that standard genetic markers, namely microsatellites, also referred as short tandem repeats (STR), and single-nucleotide polymorphism (SNP) have limited power in addressing common questions of forensic and medical genetics. To address this issue, we developed a molecular marker, named DIP-STR that relies on pairing deletion-insertion polymorphisms (DIP) with STR. This novel analytical approach allows for the unambiguous genotyping of a minor component in the presence of a major component, where DIP-STR genotypes of the minor were successfully procured at ratios up to 1:1,000. The compound nature of this marker generates a high level of polymorphism that is suitable for identity testing. Here, we demonstrate the power of the DIP-STR approach on an initial set of nine markers surveyed in a Swiss population. Finally, we discuss the limitations and potential applications of our new system including preliminary tests on clinical samples and estimates of their performance on simulated DNA mixtures.

Genomic characterization of a Brazilian TT Virus isolate closely related to SEN virus-F

SEN virus (SENV) is a circular, single stranded DNA virus that has been first characterized in the serum of a human immunodeficiency virus type 1 (HIV-1)-infected patient. Eight genotypes of SENV (A-H) have been identified and further recognized as variants of TT virus (TTV) in the family Circoviridae. Here we describe the first genomic characterization of a SENV isolate (5-A) from South America. Using 'universal' primers, able to amplify most, if not all, TTV/SENV genotypes, a segment of > 3 kb was amplified by polymerase chain reaction from the serum of an HIV-1 infected patient. The amplicon was cloned and a 3087-nucleotide sequence was determined, that showed a high (85%) homology with the sequence of the Italian isolate SENV-F. Proteins encoded by open reading frames (ORFs) 1 to 4 consisted of 758, 129, 276, and 267 amino acids, respectively. By phylogenetic analysis, isolate 5-A was classified into TTV genotype 19 (phylogenetic group 3), together with SENV-F and TTV isolate SAa-38.

Detection of Brazilian hantavirus by reverse transcription polymerase chain reaction amplification of N gene in patients with hantavirus cardiopulmonary syndrome

We report a nested reverse transcription-polymerase chain reaction (RT-PCR) assay for hantavirus using primers selected to match high homology regions of hantavirus genomes detected from the whole blood of hantavirus cardiopulmonary syndrome (HCPS) patients from Brazil, also including the N gene nucleotide sequence of Araraquara virus. Hantavirus genomes were detected in eight out of nine blood samples from the HCPS patients by RT-PCR (88.9% positivity) and in all 9 blood samples (100% positivity) by nested-PCR. The eight amplicons obtained by RT-PCR (P1, P3-P9), including one obtained by nested-PCR (P-2) and not obtained by RT-PCR, were sequenced and showed high homology (94.8% to 99.1%) with the N gene of Araraquara hantavirus. Although the serologic method ELISA is the most appropriate test for HCPS diagnosis, the use of nested RT-PCR for hantavirus in Brazil would contribute to the diagnosis of acute hantavirus disease detecting viral genomes in patient specimens as well as initial genomic characterization of circulating hantaviruses.

Testing advances in molecular discrimination among Chinook salmon life histories: evidence from a blind test.

The application of DNA-based markers toward the task of discriminating among alternate salmon runs has evolved in accordance with ongoing genomic developments and increasingly has enabled resolution of which genetic markers associate with important life-history differences. Accurate and efficient identification of the most likely origin for salmon encountered during ocean fisheries, or at salvage from fresh water diversion and monitoring facilities, has far-reaching consequences for improving measures for management, restoration and conservation. Near-real-time provision of high-resolution identity information enables prompt response to changes in encounter rates. We thus continue to develop new tools to provide the greatest statistical power for run identification. As a proof of concept for genetic identification improvements, we conducted simulation and blind tests for 623 known-origin Chinook salmon (Oncorhynchus tshawytscha) to compare and contrast the accuracy of different population sampling baselines and microsatellite loci panels. This test included 35 microsatellite loci (1266 alleles), some known to be associated with specific coding regions of functional significance, such as the circadian rhythm cryptochrome genes, and others not known to be associated with any functional importance. The identification of fall run with unprecedented accuracy was demonstrated. Overall, the top performing panel and baseline (HMSC21) were predicted to have a success rate of 98%, but the blind-test success rate was 84%. Findings for bias or non-bias are discussed to target primary areas for further research and resolution.

An evaluation of the inappropriate prescribing in older residents in long term care facilities in the greater Cork and Northern Ireland regions using the STOPP and Beers' criteria.

This work has highlighted a number of areas of prescribing concern, for example, the long term use of both benzodiazepines and hypnotics, in older residents residing in long term care facilities. Each of these individual areas should be further investigated to determine the underlying reason(s) for the prescribing concerns in these areas and strategic methods of addressing and preventing further issues should be developed on a national level.This resource was contributed by The National Documentation Centre on Drug Use.

Evolution of a cellular immune response in Drosophila: a phenotypic and genomic comparative analysis.

Understanding the genomic basis of evolutionary adaptation requires insight into the molecular basis underlying phenotypic variation. However, even changes in molecular pathways associated with extreme variation, gains and losses of specific phenotypes, remain largely uncharacterized. Here, we investigate the large interspecific differences in the ability to survive infection by parasitoids across 11 Drosophila species and identify genomic changes associated with gains and losses of parasitoid resistance. We show that a cellular immune defense, encapsulation, and the production of a specialized blood cell, lamellocytes, are restricted to a sublineage of Drosophila, but that encapsulation is absent in one species of this sublineage, Drosophila sechellia. Our comparative analyses of hemopoiesis pathway genes and of genes differentially expressed during the encapsulation response revealed that hemopoiesis-associated genes are highly conserved and present in all species independently of their resistance. In contrast, 11 genes that are differentially expressed during the response to parasitoids are novel genes, specific to the Drosophila sublineage capable of lamellocyte-mediated encapsulation. These novel genes, which are predominantly expressed in hemocytes, arose via duplications, whereby five of them also showed signatures of positive selection, as expected if they were recruited for new functions. Three of these novel genes further showed large-scale and presumably loss-of-function sequence changes in D. sechellia, consistent with the loss of resistance in this species. In combination, these convergent lines of evidence suggest that co-option of duplicated genes in existing pathways and subsequent neofunctionalization are likely to have contributed to the evolution of the lamellocyte-mediated encapsulation in Drosophila.

Trypanosoma rangeli sialidase: lack of activity in stocks from Panama and other regions

A total of 33 crude and cloned Trypanosoma rangeli stocks found as natural infections in human from Panama and other endemic areas of Central and South America were evaluated as producers of sialidase (SA) activity through the MU-NANA fluorescence test. Negative results were observed in 6 of the isolates: Panama (4), Honduras (1), and Brazil (1). In addition, an immunoblotting analysis confirm the presence of the SA antigen in these stocks without enzymatic activity. These findings must be considered in the interpretation of the biological significance of T. rangeli SA and in the proper characterization and identification of this parasite.

The role of CTCFL/BORIS in Meiosis

The complexity of mammalian genome organization demands a complex interplay of DNA and proteins to orchestrate proper gene regulation. CTCF, a highly conserved, ubiquitously expressed protein has been postulated as a primary organizer of genome architecture because of its roles in transcriptional activation/repression, insulation and imprinting. Diverse regulatory functions are exerted through genome wide binding via a central eleven zinc finger DNA binding domain and an array of diverse protein-protein interactions through N- and C- terminal domains. CTCFL has been identified as a paralog of CTCF expressed only in spermatogenic cells of the testis. CTCF and CTCFL have a highly homologous DNA-binding domain, while the flanking amino acid sequences exhibit no significant similarity. Genome- wide mapping of CTCF binding sites has been carried out in many cell types, but no data exist for CTCFL apart from a few identified loci. The lack of high quality antibodies prompted us to generate an endogenously flag-tagged CTCFL mouse model using BAC recombination. IHC staining using anti-flag antibodies confirmed CTCFL localization to type Β spermatogonia and preleptotene spermatocytes and a mutually exclusive pattern of expression with CTCF. ChIP followed by high-throughput sequencing identified 10,382 binding sites showing 70% overlap but representing only 20% of CTCF sites. Consensus sequence analysis identified a significantly longer binding motif with prominently less ambiguity of base calling at every position. The significant difference between CTCF and CTCFL genomic binding patterns proposes that their binding to DNA is differentially regulated. Analysis of CTCFL binding to methylated regions on a genome wide scale identified approximately 1,000 loci. Methylation-independent binding of CTCFL might be at least one of the mechanisms that ensures distinct binding patterns of CTCF and CTCFL since CTCF binding is methylation- sensitive. Co-localization of CTCF with cohesin has been well established and analysis of CTCFL and SMC3 overlap identified around 3,300 binding sites from which two related but distinct consensus sequence motifs were derived. Because virtually all data for cohesin binding originate from mitotically proliferating cells, the anticipated overlap is expected to be considerably higher in meiotic cells. Meiosis-specific cohesin subunit Rec8 is specific for spermatocytes and 6 out of the 12 identified binding sites are also bound by CTCFL. In conclusion, this was the first genome-wide mapping of CTCFL binding sites in spermatocytes, the only cell type where CTCF is not expressed. CTCFL has a unique binding site repertoire distinct from CTCF, binds to methylated sequences and shows a significant overlap with cohesin binding sites. Future efforts will be oriented towards deciphering the role CTCFL plays in conversion of chromatin structure and function from mitotic to meiotic chromosomes. - La complexité de l'organisation du génome des mammifères exige une interaction particulière entre ADN et protéines pour orchestrer une régulation appropriée de l'expression des gènes. CTCFL, une protéine ubiquitaire très conservée, serait le principal organisateur de l'architecture du génome de par son rôle dans l'activation / la répression de la transcription, la protection et la localisation des gènes. Diverses régulations sont opérées, d'une part au travers d'interactions à différents endroits du génome par le biais d'un domaine protéique central de liaison à l'ADN à onze doigts de zinc, et d'autre part par des interactions protéine-protéine variées au niveau de leur domaine N- et C-terminal. CTCFL a été identifié comme un paralogue de CTCF exprimé uniquement dans les cellules spermatiques du testicule. CTCFL et CTCF ont un domaine de liaison à l'ADN très homologue, tandis que les séquences d'acides aminés situées de part et d'autre de ce domaine ne présentent aucune similitude. Une cartographie générale des sites de liaison au CTCF a été réalisée pour de nombreux types cellulaires, mais il n'existe aucune donnée pour CTCFL à l'exception de l'identification de quelques loci. L'absence d'anticorps de bonne qualité nous a conduit à générer un modèle murin portant un CTCFL endogène taggué grâce à un procédé de recombinaison BAC. Une coloration IHC à l'aide d'anticorps anti-FLAG a confirmé la présence de CTCFL au niveau des spermatogonies de type Β et des spermatocytes au stade préleptotène, et une distribution mutuellement exclusive avec CTCF. Une méthode de Chromatine Immunoprecipitation (ChIP) suivie d'un séquençage à haut débit a permis d'identifier 10.382 sites de liaison montrant 70% d'homologie mais ne représentant que 20% des sites CTCF. L'analyse de la séquence consensus révèle un motif de fixation à l'ADN nettement plus long et qui comporte bien moins de bases aléatoires à chaque position nucléotidique. La différence significative entre les séquences génomiques des sites de liaison au CTCF et CTCFL suggère que leur fixation à l'ADN est régulée différemment. Appliquée à l'échelle du génome, l'étude de l'interaction de CTCFL avec des régions méthylées de l'ADN a permis d'identifier environ 1.000 loci. Contrairement à CTCFL, la liaison de CTCF dépend de l'état de méthylation de l'ADN ; cette modification épigénétique constitue donc au moins un des mécanismes de régulation expliquant une localisation de CTCF et CTCFL à des sites distincts du génome. La co- localisation de CTCF avec la cohésine étant établie, l'analyse de la superposition des séquences de CTCFL avec la sous-unité SMC3 identifie environ 3.300 sites de liaison parmi lesquels deux mêmes motifs consensus distincts par leur séquence sont mis en évidence. La presque quasi-totalité des données sur la cohésine ayant été établie à partir de cellules en prolifération mitotique, il est probable que la similitude au sein des séquences consensus soit encore plus grande dans le cas des cellules en méiose. La sous-unité Rec8 de la cohésine propre à l'état de méiose est spécifiquement exprimée dans les spermatocytes. Or 6 des 12 sites de liaison identifiés sont également utilisés par CTCFL. Pour conclure, ce travail constitue la première cartographie à l'échelle du génome des sites de liaison de CTCFL dans les spermatocytes, seul type cellulaire où CTCFL n'est pas exprimé. CTCFL possède un répertoire unique de sites de fixation à l'ADN distinct de CTCF, se lie à des séquences méthylées et présente un nombre important de sites de liaison communs avec la cohésine. Les perspectives futures sont d'élucider le rôle de CTCFL dans le remodelage de la structure de la chromatine et de définir sa fonction dans le processus de méiose.