In silico single strand melting curve: a new approach to identify nucleic acid polymorphisms in Totiviridae

Background

The PCR technique and its variations have been increasingly used in the clinical laboratory and recent advances in this field generated new higher resolution techniques based on nucleic acid denaturation dynamics. The principle of these new molecular tools is based on the comparison of melting profiles, after denaturation of a DNA double strand. Until now, the secondary structure of single-stranded nucleic acids has not been exploited to develop identification systems based on PCR. To test the potential of single-strand RNA denaturation as a new alternative to detect specific nucleic acid variations, sequences from viruses of the Totiviridae family were compared using a new in silico melting curve approach. This family comprises double-stranded RNA virus, with a genome constituted by two ORFs, ORF1 and ORF2, which encodes the capsid/RNA binding proteins and an RNA-dependent RNA polymerase (RdRp), respectively.

Results

A phylogenetic tree based on RdRp amino acid sequences was constructed, and eight monophyletic groups were defined. Alignments of RdRp RNA sequences from each group were screened to identify RNA regions with conserved secondary structure. One region in the second half of ORF2 was identified and individually modeled using the RNAfold tool. Afterwards, each DNA or RNA sequence was denatured in silico using the softwares MELTSIM and RNAheat that generate melting curves considering the denaturation of a double stranded DNA and single stranded RNA, respectively. The same groups identified in the RdRp phylogenetic tree were retrieved by a clustering analysis of the melting curves data obtained from RNAheat. Moreover, the same approach was used to successfully discriminate different variants of Trichomonas vaginalis virus, which was not possible by the visual comparison of the double stranded melting curves generated by MELTSIM.

Conclusion

In silico analysis indicate that ssRNA melting curves are more informative than dsDNA melting curves. Furthermore, conserved RNA structures may be determined from analysis of individuals that are phylogenetically related, and these regions may be used to support the reconstitution of their phylogenetic groups. These findings are a robust basis for the development of in vitro systems to ssRNA melting curves detection.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-243) contains supplementary material, which is available to authorized users.     

Positive selection on apoptosis related genes

Apoptosis is a form of programmed cell death crucial for development, homeostasis, immunity, spermatogenesis, and prevention of cancer. Positive selection acting on mammalian apoptosis related genes targets protein interfaces that interact with pathogens and also elements of signaling complexes. Selection appears primarily to be driven by the immune/defense related function of these genes. Moreover, competitive interactions could be driving positive selection among sperm cells, as well as the need for protection against female anti-sperm immune responses. Trade-offs in fitness are expected out of these selective pressures, which could explain the involvement of these genes in various diseases, including cancer.     

Regulation of IFN response gene activity during infliximab treatment in rheumatoid arthritis is associated with clinical response to treatment

Introduction  

Cross-regulation between TNF and type I IFN has been postulated to play an important role in autoimmune diseases. Therefore, we determined the effect of TNF blockade in rheumatoid arthritis (RA) on the type I IFN response gene activity in relation to clinical response.     

Disease-modifying antirheumatic drugs are associated with a reduced risk for cardiovascular disease in patients with rheumatoid arthritis: a case control study

Rheumatoid arthritis (RA) is characterized by inflammation and an increased risk for cardiovascular disease (CVD). This study investigates possible associations between CVD and the use of conventional disease-modifying antirheumatic drugs (DMARDs) in RA. Using a case control design, 613 RA patients (5,649 patient-years) were studied, 72 with CVD and 541 without CVD. Data on RA, CVD and drug treatment were evaluated from time of RA diagnosis up to the first cardiovascular event or the end of the follow-up period. The dataset was categorized according to DMARD use: sulfasalazine (SSZ), hydroxychloroquine (HCQ) or methotrexate (MTX). Odds ratios (ORs) for CVD, corrected for age, gender, smoking and RA duration, were calculated per DMARD group. Patients who never used SSZ, HCQ or MTX were used as a reference group. MTX treatment was associated with a significant CVD risk reduction, with ORs (95% CI): 'MTX only', 0.16 (0.04 to 0.66); 'MTX and SSZ ever', 0.20 (0.08 to 0.51); and 'MTX, SSZ and HCQ ever', 0.20 (0.08 to 0.54). The risk reductions remained significant after additional correction for the presence of rheumatoid factor and erosions. After correction for hypertension, diabetes and hypercholesterolemia, 'MTX or SSZ ever' and 'MTX, SSZ and HCQ ever' showed significant CVD risk reduction. Rheumatoid factor positivity and erosions both increased CVD risk, with ORs of 2.04 (1.02 to 4.07) and 2.36 (0.92 to 6.08), respectively. MTX and, to a lesser extent, SSZ were associated with significantly lower CVD risk compared to RA patients who never used SSZ, HCQ or MTX. We hypothesize that DMARD use, in particular MTX use, results in powerful suppression of inflammation, thereby reducing the development of atherosclerosis and subsequently clinically overt CVD.     

Avian population consequences of climate change are most severe for long-distance migrants in seasonal habitats

One consequence of climate change is an increasing mismatch between timing of food requirements and food availability. Such a mismatch is primarily expected in avian long-distance migrants because of their complex annual cycle, and in habitats with a seasonal food peak. Here we show that insectivorous long-distance migrant species in The Netherlands declined strongly (1984–2004) in forests, a habitat characterized by a short spring food peak, but that they did not decline in less seasonal marshes. Also, within generalist long-distance migrant species, populations declined more strongly in forests than in marshes. Forest-inhabiting migrant species arriving latest in spring declined most sharply, probably because their mismatch with the peak in food supply is greatest. Residents and short-distance migrants had non-declining populations in both habitats, suggesting that habitat quality did not deteriorate. Habitat-related differences in trends were most probably caused by climate change because at a European scale, long-distance migrants in forests declined more severely in western Europe, where springs have become considerably warmer, when compared with northern Europe, where temperatures during spring arrival and breeding have increased less. Our results suggest that trophic mismatches may have become a major cause for population declines in long-distance migrants in highly seasonal habitats.     

A Model-Based Analysis of GC-Biased Gene Conversion in the Human and Chimpanzee Genomes

GC-biased gene conversion (gBGC) is a recombination-associated process that favors the fixation of G/C alleles over A/T alleles. In mammals, gBGC is hypothesized to contribute to variation in GC content, rapidly evolving sequences, and the fixation of deleterious mutations, but its prevalence and general functional consequences remain poorly understood. gBGC is difficult to incorporate into models of molecular evolution and so far has primarily been studied using summary statistics from genomic comparisons. Here, we introduce a new probabilistic model that captures the joint effects of natural selection and gBGC on nucleotide substitution patterns, while allowing for correlations along the genome in these effects. We implemented our model in a computer program, called phastBias, that can accurately detect gBGC tracts about 1 kilobase or longer in simulated sequence alignments. When applied to real primate genome sequences, phastBias predicts gBGC tracts that cover roughly 0.3% of the human and chimpanzee genomes and account for 1.2% of human-chimpanzee nucleotide differences. These tracts fall in clusters, particularly in subtelomeric regions; they are enriched for recombination hotspots and fast-evolving sequences; and they display an ongoing fixation preference for G and C alleles. They are also significantly enriched for disease-associated polymorphisms, suggesting that they contribute to the fixation of deleterious alleles. The gBGC tracts provide a unique window into historical recombination processes along the human and chimpanzee lineages. They supply additional evidence of long-term conservation of megabase-scale recombination rates accompanied by rapid turnover of hotspots. Together, these findings shed new light on the evolutionary, functional, and disease implications of gBGC. The phastBias program and our predicted tracts are freely available.