Analysis of hepatitis C virus resistance to silibinin in vitro and in vivo points to a novel mechanism involving nonstructural protein 4B.

Intravenous silibinin (SIL) is an approved therapeutic that has recently been applied to patients with chronic hepatitis C, successfully clearing hepatitis C virus (HCV) infection in some patients even in monotherapy. Previous studies suggested multiple antiviral mechanisms of SIL; however, the dominant mode of action has not been determined. We first analyzed the impact of SIL on replication of subgenomic replicons from different HCV genotypes in vitro and found a strong inhibition of RNA replication for genotype 1a and genotype 1b. In contrast, RNA replication and infection of genotype 2a were minimally affected by SIL. To identify the viral target of SIL we analyzed resistance to SIL in vitro and in vivo. Selection for drug resistance in cell culture identified a mutation in HCV nonstructural protein (NS) 4B conferring partial resistance to SIL. This was corroborated by sequence analyses of HCV from a liver transplant recipient experiencing viral breakthrough under SIL monotherapy. Again, we identified distinct mutations affecting highly conserved amino acid residues within NS4B, which mediated phenotypic SIL resistance also in vitro. Analyses of chimeric viral genomes suggest that SIL might target an interaction between NS4B and NS3/4A. Ultrastructural studies revealed changes in the morphology of viral membrane alterations upon SIL treatment of a susceptible genotype 1b isolate, but not of a resistant NS4B mutant or genotype 2a, indicating that SIL might interfere with the formation of HCV replication sites. CONCLUSION: Mutations conferring partial resistance to SIL treatment in vivo and in cell culture argue for a mechanism involving NS4B. This novel mode of action renders SIL an attractive candidate for combination therapies with other directly acting antiviral drugs, particularly in difficult-to-treat patient cohorts.

Analysis of HIV-1 expression level and sense of transcription by high-throughput sequencing of the infected cell.

Next-generation sequencing offers an unprecedented opportunity to jointly analyze cellular and viral transcriptional activity without prerequisite knowledge of the nature of the transcripts. SupT1 cells were infected with a vesicular stomatitis virus G envelope protein (VSV-G)-pseudotyped HIV vector. At 24 h postinfection, both cellular and viral transcriptomes were analyzed by serial analysis of gene expression followed by high-throughput sequencing (SAGE-Seq). Read mapping resulted in 33 to 44 million tags aligning with the human transcriptome and 0.23 to 0.25 million tags aligning with the genome of the HIV-1 vector. Thus, at peak infection, 1 transcript in 143 is of viral origin (0.7%), including a small component of antisense viral transcription. Of the detected cellular transcripts, 826 (2.3%) were differentially expressed between mock- and HIV-infected samples. The approach also assessed whether HIV-1 infection modulates the expression of repetitive elements or endogenous retroviruses. We observed very active transcription of these elements, with 1 transcript in 237 being of such origin, corresponding on average to 123,123 reads in mock-infected samples (0.40%) and 129,149 reads in HIV-1-infected samples (0.45%) mapping to the genomic Repbase repository. This analysis highlights key details in the generation and interpretation of high-throughput data in the setting of HIV-1 cellular infection.

Comprehensive analysis of RET common and rare variant patientsin a series of Spanish Hirschsprung patients confirms a synergistic effect of both kinds of events

Background. RET is the major gene associated to Hirschsprung disease (HSCR) with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In the present study, we have performed a comprehensive study of our HSCR series evaluating the involvement of both RET rare variants (RVs) and common variants (CVs) in the context of the disease. Methods. RET mutational screening was performed by dHPLC and direct sequencing for the identification of RVs. In addition Taqman technology was applied for the genotyping of 3 RET CVs previously associated to HSCR, including a variant lying in an enhancer domain within RET intron 1 (rs2435357). Statistical analyses were performed using the SPSS v.17.0 to analyze the distribution of the variants. Results. Our results confirm the strongest association to HSCR for the "enhancer" variant, and demonstrate a significantly higher impact of it in male versus female patients. Integration of the RET RVs and CVs analysis showed that in 91.66% of cases with both kinds of mutational events, the enhancer allele is in trans with the allele bearing the RET RV. Conclusions. A gender effect exists on both the transmission and distribution of rare coding and common HSCR causing mutations. In addition, these RET CVs and RVs seem to act in a synergistic way leading to HSCR phenotype.

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.

Functional and transcriptional induction of aquaporin-1 gene by hypoxia; analysis of promoter and role of Hif-1α.

Aquaporin-1 (AQP1) is a water channel that is highly expressed in tissues with rapid O(2) transport. It has been reported that this protein contributes to gas permeation (CO(2), NO and O(2)) through the plasma membrane. We show that hypoxia increases Aqp1 mRNA and protein levels in tissues, namely mouse brain and lung, and in cultured cells, the 9L glioma cell line. Stopped-flow light-scattering experiments confirmed an increase in the water permeability of 9L cells exposed to hypoxia, supporting the view that hypoxic Aqp1 up-regulation has a functional role. To investigate the molecular mechanisms underlying this regulatory process, transcriptional regulation was studied by transient transfections of mouse endothelial cells with a 1297 bp 5' proximal Aqp1 promoter-luciferase construct. Incubation in hypoxia produced a dose- and time-dependent induction of luciferase activity that was also obtained after treatments with hypoxia mimetics (DMOG and CoCl(2)) and by overexpressing stabilized mutated forms of HIF-1α. Single mutations or full deletions of the three putative HIF binding domains present in the Aqp1 promoter partially reduced its responsiveness to hypoxia, and transfection with Hif-1α siRNA decreased the in vitro hypoxia induction of Aqp1 mRNA and protein levels. Our results indicate that HIF-1α participates in the hypoxic induction of AQP1. However, we also demonstrate that the activation of Aqp1 promoter by hypoxia is complex and multifactorial and suggest that besides HIF-1α other transcription factors might contribute to this regulatory process. These data provide a conceptual framework to support future research on the involvement of AQP1 in a range of pathophysiological conditions, including edema, tumor growth, and respiratory diseases.

Clinical pharmacogenomic testing of KRAS, BRAF and EGFR mutations by high resolution melting analysis and ultra-deep pyrosequencing

BACKGROUND: Epidermal growth factor receptor (EGFR) and its downstream factors KRAS and BRAF are mutated in several types of cancer, affecting the clinical response to EGFR inhibitors. Mutations in the EGFR kinase domain predict sensitivity to the tyrosine kinase inhibitors gefitinib and erlotinib in lung adenocarcinoma, while activating point mutations in KRAS and BRAF confer resistance to the anti-EGFR monoclonal antibody cetuximab in colorectal cancer. The development of new generation methods for systematic mutation screening of these genes will allow more appropriate therapeutic choices. METHODS: We describe a high resolution melting (HRM) assay for mutation detection in EGFR exons 19-21, KRAS codon 12/13 and BRAF V600 using formalin-fixed paraffin-embedded samples. Somatic variation of KRAS exon 2 was also analysed by massively parallel pyrosequencing of amplicons with the GS Junior 454 platform. RESULTS: We tested 120 routine diagnostic specimens from patients with colorectal or lung cancer. Mutations in KRAS, BRAF and EGFR were observed in 41.9%, 13.0% and 11.1% of the overall samples, respectively, being mutually exclusive. For KRAS, six types of substitutions were detected (17 G12D, 9 G13D, 7 G12C, 2 G12A, 2 G12V, 2 G12S), while V600E accounted for all the BRAF activating mutations. Regarding EGFR, two cases showed exon 19 deletions (delE746-A750 and delE746-T751insA) and another two substitutions in exon 21 (one showed L858R with the resistance mutation T590M in exon 20, and the other had P848L mutation). Consistent with earlier reports, our results show that KRAS and BRAF mutation frequencies in colorectal cancer were 44.3% and 13.0%, respectively, while EGFR mutations were detected in 11.1% of the lung cancer specimens. Ultra-deep amplicon pyrosequencing successfully validated the HRM results and allowed detection and quantitation of KRAS somatic mutations. CONCLUSIONS: HRM is a rapid and sensitive method for moderate-throughput cost-effective screening of oncogene mutations in clinical samples. Rather than Sanger sequence validation, next-generation sequencing technology results in more accurate quantitative results in somatic variation and can be achieved at a higher throughput scale.

A new approach for potential drug target discovery through in silico metabolic pathway analysis using Trypanosoma cruzi genome information

The current drug options for the treatment of chronic Chagas disease have not been sufficient and high hopes have been placed on the use of genomic data from the human parasite Trypanosoma cruzi to identify new drug targets and develop appropriate treatments for both acute and chronic Chagas disease. However, the lack of a complete assembly of the genomic sequence and the presence of many predicted proteins with unknown or unsure functions has hampered our complete view of the parasite's metabolic pathways. Moreover, pinpointing new drug targets has proven to be more complex than anticipated and has revealed large holes in our understanding of metabolic pathways and their integrated regulation, not only for this parasite, but for many other similar pathogens. Using an in silicocomparative study on pathway annotation and searching for analogous and specific enzymes, we have been able to predict a considerable number of additional enzymatic functions in T. cruzi. Here we focus on the energetic pathways, such as glycolysis, the pentose phosphate shunt, the Krebs cycle and lipid metabolism. We point out many enzymes that are analogous to those of the human host, which could be potential new therapeutic targets.

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.

An isoform of microtubule-associated protein 2 (MAP2) containing four repeats of the tubulin-binding motif.

Microtubule-associated protein 2 (MAP2) exists in both high- and low-molecular mass isoforms, each of which has a tubulin-binding domain consisting of 3 imperfect tandem repeats of 31 amino acids containing a more highly conserved 18 amino acid 'core' sequence. We describe here a novel form of low molecular mass MAP2 (MAP2c) that contains an additional 4th repeat of this tubulin-binding motif. Like the 3 previously known repeat sequences, this 4th copy is highly conserved between MAP2 and the two other known members of the same gene family, tau and MAP4. In each of these three genes the additional 4th repeat is inserted between the 1st and 2nd repeats of the 3-repeat form of the molecule. Experiments with brain cell cultures, in which the relative proportions of neurons and glia had been manipulated by drug treatment, showed that 4-repeat MAP2c is associated with glial cells whereas 3-repeat MAP2c is expressed in neurons. Whereas 3-repeat MAP2c is expressed early in development and then declines, the level of 4-repeat MAP2c increases later in development, corresponding to the relatively late differentiation of glial cells compared to neurons. When transfected into non-neuronal cells, the 4-repeat version of MAP2c behaved indistinguishably from the 3-repeat form in stabilising and rearranging cellular microtubules. The presence of an additional 4th repeat of the tubulin-binding motif in all three members of the MAP2 gene family suggests that this variant arose prior to their differentiation from an ancestral gene.

Estimating the motion of an underwater robot from a monocular image sequence

When underwater vehicles perform navigation close to the ocean floor, computer vision techniques can be applied to obtain quite accurate motion estimates. The most crucial step in the vision-based estimation of the vehicle motion consists on detecting matchings between image pairs. Here we propose the extensive use of texture analysis as a tool to ameliorate the correspondence problem in underwater images. Once a robust set of correspondences has been found, the three-dimensional motion of the vehicle can be computed with respect to the bed of the sea. Finally, motion estimates allow the construction of a map that could aid to the navigation of the robot

Structural analysis of TCR-ligand interactions studied on H-2Kd-restricted cloned CTL specific for a photoreactive peptide derivative.

To study the interaction of the TCR with its ligand, the complex of a MHC molecule and an antigenic peptide, we modified a TCR contact residue of a H-2Kd-restricted antigenic peptide with photoreactive 4-azidobenzoic acid. The photoreactive group was a critical component of the epitope recognized by CTL clones derived from mice immunized with such a peptide derivative. The majority of these clones expressed V beta 1-encoded beta chains that were paired with J alpha TA28-encoded alpha chains. For one of these TCR, the photoaffinity labeled sites were mapped on the alpha chain as a J alpha TA28-encoded tryptophan and on the beta chain as a residue of the C' strand of V beta 1. Molecular modeling of this TCR suggested the presence of a hydrophobic pocket that harbors this tryptophan as well as a tyrosine on the C' strand of V beta 1 between which the photoreactive side chain inserts. It is concluded that this avid binding principle may account for the preferential selection of V beta 1 and J alpha TA28-encoded TCR.

Analysis of the cytolytic T lymphocyte response of melanoma patients to the naturally HLA-A*0201-associated tyrosinase peptide 368-376.

The human tyrosinase gene codes for two distinct antigens that are recognized by HLA-A*0201-restricted CTLs. For one of them, tyrosinase peptide 368-376, the sequence identified by mass spectrometry in melanoma cell eluates differs from the gene-encoded sequence as a result of posttranslational modification of amino acid residue 370 (asparagine to aspartic acid). Here, we used fluorescent tetrameric complexes ("tetramers") of HLA-A*0201 and tyrosinase peptide 368-376 (YMDGTMSQV) to characterize the CD8+ T-cell response to this antigen in lymphoid cell populations from HLA-A2 melanoma patients. Taking advantage of the presence of significant numbers of tetramer-positive CD8+ T cells in tumor-infiltrated lymph node cells from a melanoma patient, we derived polyclonal and monoclonal tyrosinase peptide 368-376-specific CTLs by tetramer-guided flow cytometric sorting. These CTLs efficiently and specifically lysed HLA-A*0201- and tyrosinase-positive melanoma cells. As assessed with tyrosinase peptide variants, the fine antigen specificity of the CTLs was quite diverse at the clonal level. Flow cytometric analysis of PBMCs stained with tetramers showed that tyrosinase peptide 368-376-specific CD8+ T cells were hardly detectable in peripheral blood of melanoma patients. However, significant numbers of such cells were detected after short-term stimulation of CD8+ lymphocytes with tyrosinase peptide 368-376 in 6 of 10 HLA-A2 melanoma patients. Taken together, these findings emphasize the significant contribution of the natural tyrosinase peptide 368-376 to the antigenic specificities recognized by the tumor-reactive CTLs that may develop in HLA-A2 melanoma patients.

Molecular analysis of non-melanoma skin cancer

SUMMARY:Cylindroma, trichoepithelioma and spiradenoma are benign tumors of hair follicle. They are caused by mutations and loss of heterozygosity in the CYLD gene. CYLD is a ubiquitously expressed, but the tumors are restricted to skin, suggesting that the tumorigenesis is influenced by skin-specific regulators and probably by mutations in other genes. The objectives of the thesis were to analyze the molecular mechanisms leading to the aforementioned tumors. In the first project, we have identified five new mutations in CYLD gene in tive families affected with different combinations of these skin appendage tumors. F our of these mutations caused the introduction of a premature stop codon in CYLD protein sequence, but one was a missense mutation changing aspartic acid 681 into glycine (D68lG), in patients exhibiting multiple trichoepitheliomas. CYLD is a deubiquitinase which can downregulate NF-κB and INK pathways through the deubiquitination of TRAF2, for example. We showed that the CYLD-D681G mutant was unable to remove polyubiquitin chains from TRAF2. We also proved that CYLD-D68lG could not inhibit TRAP 2- or TNFα- mediated NF-κB or INK activations in 293T cells. These results underlined the importance of the D68l residue for the enzymatic activity of CYLD. TRAP-interacting protein (TRIP), which is a E3-Ubiquitin ligase, is a partner of CYLD. In the second project of the thesis, we studied the function of TRIP in the epidermis. We found that TRIP was a nucleolar protein in cultured human primary keratinocytes (HEK) and HeLa cells, and was detected in the midbody of HeLa cells. Moreover, TRIP expression was shown to be downregulated through a PKC-dependent mechanism before induction of keratinocyte differentiation. We also proved that TRIP was upregulated in basal cell carcinomas. Furthermore, TRIP was found to be important for keratinocyte survival and proliferation through the regulation of the Gl/S transition. Our results suggest that TRIP may be involved in keratinocyte tumorigenesis.RÉSUMÉ :Les cylindromes, trichoépithéliomes et spiradénomes sont des tumeurs bénignes du follicule pileux causées par des mutations et une perte d'hétérozygotie du gène CYLD. CYLD est ubiquitaire mais les tumeurs sont limitées à la peau, suggérant que la tumorigénèse est influencée par des protéines spécifiques de la peau et par des mutations dans d'autres gènes. Les objectifs de la thèse étaient d'2malyser les mécanismes moléculaires aboutissant à la formation de ces tumeurs. Dans le premier projet, cinq nouvelles mutations du gène CYLD ont été identifiées chez cinq familles présentant différentes combinaisons des tumeurs citées ci- dessus. Quatre de ces mutations causaient I' introduction d'un codon stop prématuré dans la séquence protéique, mais une était une mutation «misser1se» changeant l'aspartate 681 en résidu glycine (D68lG) chez des patients présentant des trichoépithéliomes multiples. CYLD est une déubiquitinase qui inhibe les voies de signalisation de NF-κB et JNK, en déubiquitinant notamment TRAF2. Nous avons montré que la protéine mutante CYLD- D68lG ne pouvait pas cliver la chaîne de poly-ubiquitines liée à TRAF2. CYLD-D68lG était aussi incapable d'inhiber l'activation de NF-κB ou de JNK induite par TRAF2 ou TNF-o dans les cellules 293T. Ces résultats ont donc souligné l'impo1tance du résidu D68l pour l'activité de CYLD. «TRAF-interacting protein (TRIP)», qui est une «E3-ubiquitin-ligase», est un partenaire de CYLD. Dans le second proj et de la thèse, nous avons étudié la fonction de TRIP dans l'épidenne. Nous avons montrépque TRIP était nucléolaire dans les cellules HeLa et les kératinocytes primaires humains en culture et était détectée dans le «midbody» des cellules HeLa. Nous avons prouvé que l'ARNm de TRIP était diminué avant l'induction de la différentiation des kératinocytes, par un mécanisme dépendent de la protéine kinase C, tandis qu'il était augmenté dans les carcinomes baso-cellulaires. Nous avons aussi montré que TRIP influençait la prolifération et la survie des kératinocytes en régulant la transition G1/S, Nos résultats suggèrent que TRIP est peut-être impliquée dans la tumorigénèse des kératinocytes. 7

CA88, a nuclear repetitive DNA sequence identified in Schistosoma mansoni, aids in the genotyping of nine Schistosoma species of medical and veterinary importance

CA88 is the first long nuclear repetitive DNA sequence identified in the blood fluke, Schistosoma mansoni. The assembled S. mansoni sequence, which contains the CA88 repeat, has 8,887 nucleotides and at least three repeat units of approximately 360 bp. In addition, CA88 also possesses an internal CA microsatellite, identified as SmBr18. Both PCR and BLAST analysis have been used to analyse and confirm the CA88 sequence in other S. mansoni sequences in the public database. PCR-acquired nuclear repetitive DNA sequence profiles from nine Schistosoma species were used to classify this organism into four genotypes. Included among the nine species analysed were five sequences of both African and Asian lineages that are known to infect humans. Within these genotypes, three of them refer to recognised species groups. A panel of four microsatellite loci, including SmBr18 and three previously published loci, has been used to characterise the nine Schistosoma species. Each species has been identified and classified based on its CA88 DNA fingerprint profile. Furthermore, microsatellite sequences and intra-specific variation have also been observed within the nine Schistosoma species sequences. Taken together, these results support the use of these markers in studying the population dynamics of Schistosoma isolates from endemic areas and also provide new methods for investigating the relationships between different populations of parasites. In addition, these data also indicate that Schistosoma magrebowiei is not a sister taxon to Schistosoma mattheei, prompting a new designation to a basal clade.

Phylogenetic analysis of Biomphalaria tenagophila (Orbigny, 1835) (Mollusca: Gastropoda)

Mitochondrial DNA of Biomphalaria tenagophila, a mollusc intermediate host of Schistosoma mansoni in Brazil, was sequenced and characterised. The genome size found for B. tenagophila was 13,722 bp and contained 13 messenger RNAs, 22 transfer RNAs (tRNA) and two ribosomal RNAs (rRNA). In addition to sequencing, the mitochondrial DNA (mtDNA) genome organization of B. tenagophila was analysed based on its content and localization of both coding and non-coding regions, regions of gene overlap and tRNA nucleotide sequences. Sequences of protein, rRNA 12S and rRNA 16S nucleotides as well as gene organization were compared between B. tenagophila and Biomphalaria glabrata, as the latter is the most important S. mansoni intermediate host in Brazil. Differences between such species were observed regarding rRNA composition. The complete sequence of the B. tenagophila mitochondrial genome was deposited in GenBank (accession EF433576). Furthermore, phylogenetic relationships were estimated among 28 mollusc species, which had their complete mitochondrial genome deposited in GenBank, using the neighbour-joining method, maximum parsimony and maximum likelihood bootstrap. B. tenagophila was positioned at a branch close to B. glabrata and Pulmonata molluscs, collectively comprising a paraphyletic group, contrary to Opistobranchia, which was positioned at a single branch and constituted a monophyletic group.