Genome-wide DNA polymorphism analyses using VariScan

BACKGROUND: DNA sequence polymorphisms analysis can provide valuable information on the evolutionary forces shaping nucleotide variation, and provides an insight into the functional significance of genomic regions. The recent ongoing genome projects will radically improve our capabilities to detect specific genomic regions shaped by natural selection. Current available methods and software, however, are unsatisfactory for such genome-wide analysis. RESULTS: We have developed methods for the analysis of DNA sequence polymorphisms at the genome-wide scale. These methods, which have been tested on a coalescent-simulated and actual data files from mouse and human, have been implemented in the VariScan software package version 2.0. Additionally, we have also incorporated a graphical-user interface. The main features of this software are: i) exhaustive population-genetic analyses including those based on the coalescent theory; ii) analysis adapted to the shallow data generated by the high-throughput genome projects; iii) use of genome annotations to conduct a comprehensive analyses separately for different functional regions; iv) identification of relevant genomic regions by the sliding-window and wavelet-multiresolution approaches; v) visualization of the results integrated with current genome annotations in commonly available genome browsers. CONCLUSION: VariScan is a powerful and flexible suite of software for the analysis of DNA polymorphisms. The current version implements new algorithms, methods, and capabilities, providing an important tool for an exhaustive exploratory analysis of genome-wide DNA polymorphism data.

DnaSP v5: A software for comprehensive analysis of DNA polymorphism data

DnaSP is a software package for a comprehensive analysis of DNA polymorphism data. Version 5 implements a number of new features and analytical methods allowing extensive DNA polymorphism analyses on large datasets. Among other features, the newly implemented methods allow for: (i) analyses on multiple data files; (ii) haplotype phasing; (iii) analyses on insertion/deletion polymorphism data; (iv) visualizing sliding window results integrated with available genome annotations in the UCSC browser.

Design of new tools to improve the efficacy of DNA-damaging drugs used in cancer therapy

Résumé Les tumeurs sont diverses et hétérogènes, mais toutes partagent la capacité de proliférer sans contrôle. Une prolifération dérégulée de cellules couplée à une insensibilité à une réponse apoptotique constitue une condition minimale pour que l'évolution d'une tumeur se produise. Un des traitements les plus utilisés pour traité le cancer à l'heure actuelle sont les chimiothérapies, qui sont fréquemment des composés chimiques qui induisent des dommages dans l'ADN. Les agents anticancéreux sont efficaces seulement quand les cellules tumorales sont plus aisément tuées que le tissu normal environnant. L'efficacité de ces agents est en partie déterminée par leur capacité à induire l'apoptose. Nous avons récemment démontré que la protéine RasGAP est un substrat non conventionnel des caspases parce elle peut induire à la fois des signaux anti et pro-apoptotiques, selon l'ampleur de son clivage par les caspases. A un faible niveau d'activité des caspases, RasGAP est clivé, générant deux fragments (le fragment N et le fragment C). Le fragment N semble être un inhibiteur général de l'apoptose en aval de l'activation des caspases. À des niveaux plus élevés d'activité des caspases, la capacité du fragment N de contrecarrer l'apoptose est supprimée quand il est clivé à nouveau par les caspases. Ce dernier clivage produit deux nouveaux fragments, N 1 et N2, qui contrairement au fragment N sensibilisent efficacement des cellules cancéreuses envers des agents chimiothérapeutiques. Dans cette étude nous avons prouvé qu'un peptide, appelé par la suite TAT-RasGAP317-326, qui est dérivé du fragment N2 de RasGAP et est rendu perméable aux cellules, sensibilise spécifiquement des cellules cancéreuses à trois génotoxines différentes utilisées couramment dans des traitements anticancéreux, et cela dans des modèles in vitro et in vivo. Il est important de noté que ce peptide semble ne pas avoir d'effet sur des cellules non cancéreuses. Nous avons également commencé à caractériser les mécanismes moléculaires expliquant les fonctions de sensibilisation de TAT-RasGAP317-326. Nous avons démontré que le facteur de transcription p53 et une protéine sous son activité transcriptionelle, nommée Puma, sont indispensables pour l'activité de TAT-RasGAP317-326. Nous avons également prouvé que TAT-RasGAP317-326 exige la présence d'une protéine appelée G3BP1, une protéine se liant a RasGAP, pour potentialisé les effets d'agents anticancéreux. Les données obtenues dans cette étude montrent qu'il pourrait être possible d'augmenter l'efficacité des chimiothérapies à l'aide d'un composé capable d'augmenter la sensibilité des tumeurs aux génotoxines et ainsi pourrait permettre de traiter de manière plus efficace des patients sous traitement chimiothérapeutiques. Summary Tumors are diverse and heterogeneous, but all share the ability to proliferate without control. Deregulated cell proliferation coupled with suppressed apoptotic sensitivity constitutes a minimal requirement upon which tumor evolution occurs. One of the most commonly used treatments is chemotherapy, which frequently uses chemical compounds that induce DNA damages. Anticancer agents are effective only when tumors cells are more readily killed than the surrounding normal tissue. The efficacy of these agents is partly determined by their ability to induce apoptosis. We have recently demonstrated that the protein RasGAP is an unconventional caspase substrate because it can induce both anti- and pro-apoptotic signals, depending on the extent of its cleavage by caspases. At low levels of caspase activity, RasGAP is cleaved, generating an N-terminal fragment (fragment N) and a C-terminal fragment (fragment C). Fragment N appears to be a general Mocker of apoptosis downstream of caspase activation. At higher levels of caspase activity, the ability of fragment N to counteract apoptosis is suppressed when it is further cleaved. This latter cleavage event generates two fragments, N1 and N2, which in contrast to fragment N potently sensitizes cancer cells toward DNA-damaging agents induced apoptosis. In the present study we show that a cell permeable peptide derived from the N2 fragment of RasGAP, thereafter called TAT-RasGAP317-326, specifically sensitizes cancer cells to three different genotoxins commonly used in chemotherapy in vitro and in vivo models. Importantly this peptide seems not to have any effect on non cancer cells. We have also started to characterize the molecular mechanisms underlying the sensitizing functions of TAT-RasGAP317-326. We have demonstrated that the p53 transcription factor and a protein under its transcriptional activity, called Puma, are required for the activity of TATRasGAP317-326. We have also showed that TAT-RasGAP317-326 requires the presence of a protein called G3BP1, which have been shown to interact with RasGAP, to increase the effect of the DNA-damaging drug cisplatin. The data obtained in this study showed that it is possible to increase the efficacy of current used chemotherapies with a compound able to increase the efficacy of genotoxins which could be beneficial for patients subjected to chemotherapy.

No longitudinal mitochondrial DNA sequence changes in HIV-infected individuals with and without lipoatrophy.

The potential for mitochondrial (mt) DNA mutation accumulation during antiretroviral therapy (ART), and preferential accumulation in patients with lipoatrophy compared with control participants, remains controversial. We sequenced the entire mitochondrial genome, both before ART and after ART exposure, in 29 human immunodeficiency virus (HIV)-infected Swiss HIV Cohort Study participants initiating a first-line thymidine analogue-containing ART regimen. No accumulation of mtDNA mutations or deletions was detected in 13 participants who developed lipoatrophy or in 16 control participants after significant and comparable ART exposure (median duration, 3.3 and 3.7 years, respectively). In HIV-infected persons, the development of lipoatrophy is unlikely to be associated with accumulation of mtDNA mutations detectable in peripheral blood.

Sedimentation of macroscopic rigid knots and its relation to gel electrophoretic mobility of DNA knots.

We address the general question of the extent to which the hydrodynamic behaviour of microscopic freely fluctuating objects can be reproduced by macrosopic rigid objects. In particular, we compare the sedimentation speeds of knotted DNA molecules undergoing gel electrophoresis to the sedimentation speeds of rigid stereolithographic models of ideal knots in both water and silicon oil. We find that the sedimentation speeds grow roughly linearly with the average crossing number of the ideal knot configurations, and that the correlation is stronger within classes of knots. This is consistent with previous observations with DNA knots in gel electrophoresis.

Molecular dating of caprines using ancient DNA sequences of Myotragus balearicus, an extinct endemic Balearic mammal

Background: Myotragus balearicus was an endemic bovid from the Balearic Islands (Western Mediterranean) that became extinct around 6,000-4,000 years ago. The Myotragus evolutionary lineage became isolated in the islands most probably at the end of the Messinian crisis, when the desiccation of the Mediterranean ended, in a geological date established at 5.35 Mya. Thus, the sequences of Myotragus could be very valuable for calibrating the mammalian mitochondrial DNA clock and, in particular, the tree of the Caprinae subfamily, to which Myotragus belongs. Results: We have retrieved the complete mitochondrial cytochrome b gene (1,143 base pairs), plus fragments of the mitochondrial 12S gene and the nuclear 28S rDNA multi-copy gene from a well preserved Myotragus subfossil bone. The best resolved phylogenetic trees, obtained with the cytochrome b gene, placed Myotragus in a position basal to the Ovis group. Using the calibration provided by the isolation of Balearic Islands, we calculated that the initial radiation of caprines can be dated at 6.2 ± 0.4 Mya. In addition, alpine and southern chamois, considered until recently the same species, split around 1.6 ± 0.3 Mya, indicating that the two chamois species have been separated much longer than previously thought. Conclusion: Since there are almost no extant endemic mammals in Mediterranean islands, the sequence of the extinct Balearic endemic Myotragus has been crucial for allowing us to use the Messinian crisis calibration point for dating the caprines phylogenetic tree.

Nuclear translocation and DNA recognition signals colocalized within the bZIP domain of cyclic adenosine 3',5'-monophosphate response element-binding protein CREB.

CREB is a cAMP-responsive nuclear DNA-binding protein that binds to cAMP response elements and stimulates gene transcription upon activation of the cAMP signalling pathway. The protein consists of an amino-terminal transcriptional transactivation domain and a carboxyl-terminal DNA-binding domain (bZIP domain) comprised of a basic region and a leucine zipper involved in DNA recognition and dimerization, respectively. Recently, we discovered a testis-specific transcript of CREB that contains an alternatively spliced exon encoding multiple stop codons. CREB encoded by this transcript is a truncated protein lacking the bZIP domain. We postulated that the antigen detected by CREB antiserum in the cytoplasm of germinal cells is the truncated CREB that must also lack its nuclear translocation signal (NTS). To test this hypothesis we prepared multiple expression plasmids encoding carboxyl-terminal deletions of CREB and transiently expressed them in COS-1 cells. By Western immunoblot analysis as well as immunocytochemistry of transfected cells, we show that CREB proteins truncated to amino acid 286 or shorter are sequestered in the cytoplasm, whereas a CREB of 295 amino acids is translocated into the nucleus. Chimeric CREBs containing a heterologous NTS fused to the first 248 or 261 amino acids of CREB are able to drive the translocation of the protein into the nucleus. Thus, the nine amino acids in the basic region involved in DNA recognition between positions 287 and 295 (RRKKKEYVK) of CREB contain the NTS. Further, mutation of the lysine at position 290 in CREB to an asparagine diminishes nuclear translocation of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Preclinical evaluation of the immunogenicity of C-type HIV-1-based DNA and NYVAC vaccines in the Balb/C mouse model.

As part of a European initiative (EuroVacc), we report the design, construction, and immunogenicity of two HIV-1 vaccine candidates based on a clade C virus strain (CN54) representing the current major epidemic in Asia and parts of Africa. Open reading frames encoding an artificial 160-kDa GagPolNef (GPN) polyprotein and the external glycoprotein gp120 were fully RNA and codon optimized. A DNA vaccine (DNA-GPN and DNA-gp120, referred to as DNA-C), and a replication-deficient vaccinia virus encoding both reading frames (NYVAC-C), were assessed regarding immunogenicity in Balb/C mice. The intramuscular administration of both plasmid DNA constructs, followed by two booster DNA immunizations, induced substantial T-cell responses against both antigens as well as Env-specific antibodies. Whereas low doses of NYVAC-C failed to induce specific CTL or antibodies, high doses generated cellular as well as humoral immune responses, but these did not reach the levels seen following DNA vaccination. The most potent immune responses were detectable using prime:boost protocols, regardless of whether DNA-C or NYVAC-C was used as the priming or boosting agent. These preclinical findings revealed the immunogenic response triggered by DNA-C and its enhancement by combining it with NYVAC-C, thus complementing the macaque preclinical and human phase I clinical studies of EuroVacc.

Cytoskeletal and DNA structure abnormalities result from bypass of requirement for the cdc10 start gene in the fission yeast Schizosaccharomyces pombe.

The cdc10 gene of the fission yeast S. pombe is required for traverse of the start control in late G1 and commitment to the mitotic cell cycle. To increase our understanding of the events which occur at start, a pseudoreversion analysis was undertaken to identify genes whose products may interact with cdc10 or bypass the requirement for it. A single gene, sct1+ (suppressor of cdc ten), has been identified, mutation of which suppresses all conditional alleles and a null allele of cdc10. Bypass of the requirement for cdc10+ function by sct1-1 mutations leads to pleiotropic defects, including microtubule, microfilament and nuclear structural abnormalities. Our data suggest that sct1 encodes a protein that is dependent upon cdc10+ either for its normal function or expression, or is a component of a checkpoint that monitors execution of p85cdc10 function.

Length variation in HIV-1 gp120 as the product of DNA misalignment mechanism

Background: To determine whether misalignment structures such as duplications, repeats, and palindromes are associated to insertions/deletions (indels) in gp120, indicating that indels are indeed frameshift mutations generated by DNA misalignment mechanism. Methods: Cloning and sequencing of a fragment of HIV-1 gp120 spanning C2-C4 derived from plasma RNA in 12 patients with early chronic disease and naïve to antiretroviral therapy. Results: Indels in V4 involved always insertion and deletion of duplicated nucleotide segments, and AAT repeats, and were associated to the presence of palindromic sequences. No duplications were detected in V3 and C3. Palindromic sequences occurred with similar frequencies in V3, C3 and V4; the frequency of palindromes in individual genes was found to be significantly higher in structural (gp120, p ≤ 3.00E-7) and significantly lower in regulatory (Tat, p ≤ 9.00E-7) genes, as compared to the average frequency calculated over the full genome. Discussion: Indels in V4 are associated to misalignment structures (i.e. duplications repeat and palindromes) indicating DNA misalignment as the mechanism underlying length variation in V4. The finding that indels in V4 are caused by DNA misalignment has some very important implications: 1) indels in V4 are likely to occur in proviral DNA (and not in RNA), after integration of HIV into the host genome; 2) they are likely to occur as progressive modifications of the early founder virus during chronic infection, as more and more cells get infected; 3) frameshift mutations involving any number of base pairs are likely to occur evenly across gp120; however, only those mutants carrying a functional gp120 (indels as multiples of three base pairs) will be able to perpetuate the virus cycle and to keep spreading through the population.

Meta-analyses identify 13 loci associated with age at menopause and highlight DNA repair and immune pathways.

To newly identify loci for age at natural menopause, we carried out a meta-analysis of 22 genome-wide association studies (GWAS) in 38,968 women of European descent, with replication in up to 14,435 women. In addition to four known loci, we identified 13 loci newly associated with age at natural menopause (at P < 5 × 10(-8)). Candidate genes located at these newly associated loci include genes implicated in DNA repair (EXO1, HELQ, UIMC1, FAM175A, FANCI, TLK1, POLG and PRIM1) and immune function (IL11, NLRP11 and PRRC2A (also known as BAT2)). Gene-set enrichment pathway analyses using the full GWAS data set identified exoDNase, NF-κB signaling and mitochondrial dysfunction as biological processes related to timing of menopause.

Probing Rad51-DNA interactions by changing DNA twist.

In eukaryotes, Rad51 protein is responsible for the recombinational repair of double-strand DNA breaks. Rad51 monomers cooperatively assemble on exonuclease-processed broken ends forming helical nucleo-protein filaments that can pair with homologous regions of sister chromatids. Homologous pairing allows the broken ends to be reunited in a complex but error-free repair process. Rad51 protein has ATPase activity but its role is poorly understood, as homologous pairing is independent of adenosine triphosphate (ATP) hydrolysis. Here we use magnetic tweezers and electron microscopy to investigate how changes of DNA twist affect the structure of Rad51-DNA complexes and how ATP hydrolysis participates in this process. We show that Rad51 protein can bind to double-stranded DNA in two different modes depending on the enforced DNA twist. The stretching mode is observed when DNA is unwound towards a helical repeat of 18.6 bp/turn, whereas a non-stretching mode is observed when DNA molecules are not permitted to change their native helical repeat. We also show that the two forms of complexes are interconvertible and that by enforcing changes of DNA twist one can induce transitions between the two forms. Our observations permit a better understanding of the role of ATP hydrolysis in Rad51-mediated homologous pairing and strand exchange.

Topological aspects of DNA function and protein folding.

The Topological Aspects of DNA Function and Protein Folding international meeting provided an interdisciplinary forum for biological scientists, physicists and mathematicians to discuss recent developments in the application of topology to the study of DNA and protein structure. It had 111 invited participants, 48 talks and 21 posters. The present article discusses the importance of topology and introduces the articles from the meeting's speakers.