CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution

Parvoviruses are rapidly evolving viruses that infect a wide range of hosts, including vertebrates and invertebrates. Extensive methylation of the parvovirus genome has been recently demonstrated. A global pattern of methylation of CpG dinucleotides is seen in vertebrate genomes, compared to “fractional” methylation patterns in invertebrate genomes. It remains unknown if the loss of CpG dinucleotides occurs in all viruses of a given DNA virus family that infect host species spanning across vertebrates and invertebrates. We investigated the link between the extent of CpG dinucleotide depletion among autonomous parvoviruses and the evolutionary lineage of the infected host. We demonstrate major differences in the relative abundance of CpG dinucleotides among autonomous parvoviruses which share similar genome organization and common ancestry, depending on the infected host species. Parvoviruses infecting vertebrate hosts had significantly lower relative abundance of CpG dinucleotides than parvoviruses infecting invertebrate hosts. The strong correlation of CpG dinucleotide depletion with the gain in TpG/CpA dinucleotides and the loss of TpA dinucleotides among parvoviruses suggests a major role for CpG methylation in the evolution of parvoviruses. Our data present evidence that links the relative abundance of CpG dinucleotides in parvoviruses to the methylation capabilities of the infected host. In sum, our findings support a novel perspective of host-driven evolution among autonomous parvoviruses.     

Systematic CpT (ApG) Depletion and CpG Excess Are Unique Genomic Signatures of Large DNA Viruses Infecting Invertebrates

Differences in the relative abundance of dinucleotides, if any may provide important clues on host-driven evolution of viruses. We studied dinucleotide frequencies of large DNA viruses infecting vertebrates (n = 105; viruses infecting mammals = 99; viruses infecting aves = 6; viruses infecting reptiles = 1) and invertebrates (n = 88; viruses infecting insects = 84; viruses infecting crustaceans = 4). We have identified systematic depletion of CpT(ApG) dinucleotides and over-representation of CpG dinucleotides as the unique genomic signature of large DNA viruses infecting invertebrates. Detailed investigation of this unique genomic signature suggests the existence of invertebrate host-induced pressures specifically targeting CpT(ApG) and CpG dinucleotides. The depletion of CpT dinucleotides among large DNA viruses infecting invertebrates is at least in part, explained by non-canonical DNA methylation by the infected host. Our findings highlight the role of invertebrate host-related factors in shaping virus evolution and they also provide the necessary framework for future studies on evolution, epigenetics and molecular biology of viruses infecting this group of hosts.     

Performance of Different Analytical Software Packages in Quantification of DNA Methylation by Pyrosequencing

Background

Pyrosequencing has emerged as an alternative method of nucleic acid sequencing, well suited for many applications which aim to characterize single nucleotide polymorphisms, mutations, microbial types and CpG methylation in the target DNA. The commercially available pyrosequencing systems can harbor two different types of software which allow analysis in AQ or CpG mode, respectively, both widely employed for DNA methylation analysis.

Objective

Aim of the study was to assess the performance for DNA methylation analysis at CpG sites of the two pyrosequencing software which allow analysis in AQ or CpG mode, respectively. Despite CpG mode having been specifically generated for CpG methylation quantification, many investigations on this topic have been carried out with AQ mode. As proof of equivalent performance of the two software for this type of analysis is not available, the focus of this paper was to evaluate if the two modes currently used for CpG methylation assessment by pyrosequencing may give overlapping results.

Methods

We compared the performance of the two software in quantifying DNA methylation in the promoter of selected genes (GSTP1, MGMT, LINE-1) by testing two case series which include DNA from paraffin embedded prostate cancer tissues (PC study, N = 36) and DNA from blood fractions of healthy people (DD study, N = 28), respectively.

Results

We found discrepancy in the two pyrosequencing software-based quality assignment of DNA methylation assays. Compared to the software for analysis in the AQ mode, less permissive criteria are supported by the Pyro Q-CpG software, which enables analysis in CpG mode. CpG mode warns the operators about potential unsatisfactory performance of the assay and ensures a more accurate quantitative evaluation of DNA methylation at CpG sites.

Conclusion

The implementation of CpG mode is strongly advisable in order to improve the reliability of the methylation analysis results achievable by pyrosequencing.     

Pathogen- and Host-Directed Antileishmanial Effects Mediated by Polyhexanide (PHMB)

Background

Cutaneous leishmaniasis (CL) is a neglected tropical disease caused by protozoan parasites of the genus Leishmania. CL causes enormous suffering in many countries worldwide. There is no licensed vaccine against CL, and the chemotherapy options show limited efficacy and high toxicity. Localization of the parasites inside host cells is a barrier to most standard chemo- and immune-based interventions. Hence, novel drugs, which are safe, effective and readily accessible to third-world countries and/or drug delivery technologies for effective CL treatments are desperately needed.

Methodology/Principal Findings

Here we evaluated the antileishmanial properties and delivery potential of polyhexamethylene biguanide (PHMB; polyhexanide), a widely used antimicrobial and wound antiseptic, in the Leishmania model. PHMB showed an inherent antileishmanial activity at submicromolar concentrations. Our data revealed that PHMB kills Leishmania major (L. major) via a dual mechanism involving disruption of membrane integrity and selective chromosome condensation and damage. PHMB’s DNA binding and host cell entry properties were further exploited to improve the delivery and immunomodulatory activities of unmethylated cytosine-phosphate-guanine oligodeoxynucleotides (CpG ODN). PHMB spontaneously bound CpG ODN, forming stable nanopolyplexes that enhanced uptake of CpG ODN, potentiated antimicrobial killing and reduced host cell toxicity of PHMB.

Conclusions

Given its low cost and long history of safe topical use, PHMB holds promise as a drug for CL therapy and delivery vehicle for nucleic acid immunomodulators.     

Non-Steroidal Anti-Inflammatory Drug Use and Genomic DNA Methylation in Blood

Background

Non-steroidal anti-inflammatory drug (NSAID) use is associated with decreased risk of some cancers. NSAID use modulates the epigenetic profile of normal colonic epithelium and may reduce risk of colon cancer through this pathway; however, the effect of NSAID use on the DNA methylation profile of other tissues including whole blood has not yet been examined.

Findings

Using the Sister Study cohort, we examined the association between NSAID usage and whole genome methylation patterns in blood DNA. Blood DNA methylation status across 27,589 CpG sites was evaluated for 871 women using the Illumina Infinium HumanMethylation27 Beadchip, and in a non-overlapping replication sample of 187 women at 485,512 CpG sites using the Infinium HumanMethylation450 Beadchip. We identified a number of CpG sites that were differentially methylated in regular, long-term users of NSAIDs in the discovery group, but none of these sites were statistically significant in our replication group.

Conclusions

We found no replicable methylation differences in blood related to NSAID usage. If NSAID use does effect blood DNA methylation patterns, differences are likely small.     

A Hybrid Approach for CpG Island Detection in the Human Genome

Background

CpG islands have been demonstrated to influence local chromatin structures and simplify the regulation of gene activity. However, the accurate and rapid determination of CpG islands for whole DNA sequences remains experimentally and computationally challenging.

Methodology/Principal Findings

A novel procedure is proposed to detect CpG islands by combining clustering technology with the sliding-window method (PSO-based). Clustering technology is used to detect the locations of all possible CpG islands and process the data, thus effectively obviating the need for the extensive and unnecessary processing of DNA fragments, and thus improving the efficiency of sliding-window based particle swarm optimization (PSO) search. This proposed approach, named ClusterPSO, provides versatile and highly-sensitive detection of CpG islands in the human genome. In addition, the detection efficiency of ClusterPSO is compared with eight CpG island detection methods in the human genome. Comparison of the detection efficiency for the CpG islands in human genome, including sensitivity, specificity, accuracy, performance coefficient (PC), and correlation coefficient (CC), ClusterPSO revealed superior detection ability among all of the test methods. Moreover, the combination of clustering technology and PSO method can successfully overcome their respective drawbacks while maintaining their advantages. Thus, clustering technology could be hybridized with the optimization algorithm method to optimize CpG island detection.

Conclusion/Significance

The prediction accuracy of ClusterPSO was quite high, indicating the combination of CpGcluster and PSO has several advantages over CpGcluster and PSO alone. In addition, ClusterPSO significantly reduced implementation time.     

A Case Control Study Reveals that Polyomaviruria Is Significantly Associated with Interstitial Cystitis and Vesical Ulceration

Objectives

To investigate whether polyomaviruses contribute to interstitial cystitis pathogenesis.

Subjects and Methods

A prospective study was performed with 50 interstitial cystitis cases compared with 50 age-matched, disease-free controls for the frequency of polyomaviruria. Associations between polyomaviruria and disease characteristics were analysed in cases. Polyomavirus in urine and bladder tissue was detected with species (JC virus vs. BK virus) specific, real-time PCR.

Results

Case patients were reflective of interstitial cystitis epidemiology with age range from 26–88 years (median 58) and female predominance (41/50 F). There was a significant increase in the frequency of polyomavirus shedding between cases and controls (p<0.02). Polyomavirus shedding, in particular BK viruria, was associated with vesical ulceration, a marker of disease severity, among interstitial cystitis cases after adjustment for age and sex (OR 6.8, 95% CI 1.89–24.4). There was a significant association among cases between the presence of BK viruria and response to intravesical Clorpactin therapy (OR 4.50, 95% CI 1.17–17.4).

Conclusion

The presence of polyomaviruria was found to be associated with the ulcerative form of interstitial cystitis. Clorpactin, which has anti-DNA virus activity, was more likely to improve symptoms in the presence of BK viruria. These data from this pilot study suggest associations between polyomaviruria and interstitial cystitis warranting further investigation.     

Patterns of DNA methylation in development, division of labor and hybridization in an ant with genetic caste determination

Background

DNA methylation is a common regulator of gene expression, including acting as a regulator of developmental events and behavioral changes in adults. Using the unique system of genetic caste determination in Pogonomyrmex barbatus, we were able to document changes in DNA methylation during development, and also across both ancient and contemporary hybridization events.

Methodology/Principal Findings

Sodium bisulfite sequencing demonstrated in vivo methylation of symmetric CG dinucleotides in P. barbatus. We also found methylation of non-CpG sequences. This validated two bioinformatics methods for predicting gene methylation, the bias in observed to expected ratio of CpG dinucleotides and the density of CpG/TpG single nucleotide polymorphisms (SNP). Frequencies of genomic DNA methylation were determined for different developmental stages and castes using ms-AFLP assays. The genetic caste determination system (GCD) is probably the product of an ancestral hybridization event between P. barbatus and P. rugosus. Two lineages obligately co-occur within a GCD population, and queens are derived from intra-lineage matings whereas workers are produced from inter-lineage matings. Relative DNA methylation levels of queens and workers from GCD lineages (contemporary hybrids) were not significantly different until adulthood. Virgin queens had significantly higher relative levels of DNA methylation compared to workers. Worker DNA methylation did not vary among developmental stages within each lineage, but was significantly different between the currently hybridizing lineages. Finally, workers of the two genetic caste determination lineages had half as many methylated cytosines as workers from the putative parental species, which have environmental caste determination.

Conclusions/Significance

These results suggest that DNA methylation may be a conserved regulatory mechanism moderating division of labor in both bees and ants. Current and historic hybridization appear to have altered genomic methylation levels suggesting a possible link between changes in overall DNA methylation and the origin and regulation of genetic caste determination in P. barbatus.     

Effect of the assignment of ancestral CpG state on the estimation of nucleotide substitution rates in mammals

Background  

Molecular evolutionary studies in mammals often estimate nucleotide substitution rates within and outside CpG dinucleotides separately. Frequently, in alignments of two sequences, the division of sites into CpG and non-CpG classes is based simply on the presence or absence of a CpG dinucleotide in either sequence, a procedure that we refer to as CpG/non-CpG assignment. Although it likely that this procedure is biased, it is generally assumed that the bias is negligible if species are very closely related.     

Dystroglycan Depletion Impairs Actin-Dependent Functions of Differentiated Kasumi-1 Cells

Background

Dystroglycan has recently been characterised in blood tissue cells, as part of the dystrophin glycoprotein complex involved in the differentiation process of neutrophils.

Purpose

In the present study we have investigated the role of dystroglycan in the human promyelocytic leukemic cell line Kasumi-1 differentiated to macrophage-like cells.

Methods

We characterised the pattern expression and subcellular distribution of dystroglycans in non-differentiated and differentiated Kasumi-1 cells.

Results

Our results demonstrated by WB and flow cytometer assays that during the differentiation process to macrophages, dystroglycans were down-regulated; these results were confirmed with qRT-PCR assays. Additionally, depletion of dystroglycan by RNAi resulted in altered morphology and reduced properties of differentiated Kasumi-1 cells, including morphology, migration and phagocytic activities although secretion of IL-1β and expression of markers of differentiation are not altered.

Conclusion

Our findings strongly implicate dystroglycan as a key membrane adhesion protein involved in actin-based structures during the differentiation process in Kasumi-1 cells.     

A Global Phylogeny of Leafmining Ectoedemia Moths (Lepidoptera: Nepticulidae): Exploring Host Plant Family Shifts and Allopatry as Drivers of Speciation

Background

Host association patterns in Ectoedemia (Lepidoptera: Nepticulidae) are also encountered in other insect groups with intimate plant relationships, including a high degree of monophagy, a preference for ecologically dominant plant families (e.g. Fagaceae, Rosaceae, Salicaceae, and Betulaceae) and a tendency for related insect species to feed on related host plant species. The evolutionary processes underlying these patterns are only partly understood, we therefore assessed the role of allopatry and host plant family shifts in speciation within Ectoedemia.

Methodology

Six nuclear and mitochondrial DNA markers with a total aligned length of 3692 base pairs were used to infer phylogenetic relationships among 92 species belonging to the subgenus Ectoedemia of the genus Ectoedemia, representing a thorough taxon sampling with a global coverage. The results support monophyletic species groups that are congruent with published findings based on morphology. We used the obtained phylogeny to explore host plant family association and geographical distribution to investigate if host shifts and allopatry have been instrumental in the speciation of these leafmining insects.

Significance

We found that, even though most species within species groups commonly feed on plants from one family, shifts to a distantly related host family have occasionally occurred throughout the phylogeny and such shifts are most commonly observed towards Betulaceae. The largest radiations have occurred within species groups that feed on Fagaceae, Rosaceae, and Salicaceae. Most species are restricted to one of the seven global biogeographic regions, but within species groups representatives are commonly found in different biogeographic regions. Although we find general patterns with regard to host use and biogeography, there are differences between clades that suggest that different drivers of speciation, and perhaps drivers that we did not examine, have shaped diversity patterns in different clades.     

Involvement of Cytochrome P450 1A1 and Glutathione S-Transferase P1 Polymorphisms and Promoter Hypermethylation in the Progression of Anti-Tuberculosis Drug-Induced Liver Injury: A Case–Control Study

Background

Anti-tuberculosis (anti-TB) drug-induced liver injury (ADLI) is one of the most common adverse effects associated with TB treatment. Cytochrome P450 (CYP) 1A1 and glutathione S-transferase (GST) P1 are important phase I/II metabolizing enzymes involved in drug metabolism and detoxification. Genetic polymorphism and CpG island methylation have been reported as factors influencing the expression of CYP1A1 and GSTP1.

Objective

This study aimed to determine the potential relationships of CYP1A1 and GSTP1 polymorphisms and CpG island methylation with ADLI risk.

Design

This was a population-based one-to-one matched case–control study.

Setting

The subjects were patients with TB receiving treatment in China from December 2010 to June 2013.

Patients

In total, 127 patients with TB and ADLI (case group) and 127 patients with TB but without liver injury (control group) were included in this study. Subjects were matched in terms of sex, age, and therapeutic regimen.

Methods

The general condition of each patient was assessed using questionnaires. The CYP1A1 MspI and GSTP1 Ile105Val polymorphisms as well as methylation status were detected by polymerase chain reaction (PCR)–restriction fragment length polymorphism and the methylation-specific PCR method.

Results

We found no significant difference in GSTP1 and CYP1A1 genotypes between the two groups, probably because the sample size was not large enough; however, patients with ADLI had significantly higher GSTP1 and CYP1A1 promoter methylation rates than control subjects [odds ratio (OR) = 2.467 and 2.000, respectively]. After adjusting for drinking, which significantly differed between the groups as per univariate analysis, we found that hypermethylation of GSTP1 and CYP1A1 promoters was associated with ADLI (OR = 2.645 and 2.090, respectively).

Conclusion

Hypermethylation of CpG islands of GSTP1 and CYP1A1 promoters may thus play important roles in the development of ADLI and provide evidence of being used as novel markers for ADLI risk prediction.     

Kidney Dysfunction in Adult Offspring Exposed In Utero to Type 1 Diabetes Is Associated with Alterations in Genome-Wide DNA Methylation

Background

Fetal exposure to hyperglycemia impacts negatively kidney development and function.

Objective

Our objective was to determine whether fetal exposure to moderate hyperglycemia is associated with epigenetic alterations in DNA methylation in peripheral blood cells and whether those alterations are related to impaired kidney function in adult offspring.

Design

Twenty nine adult, non-diabetic offspring of mothers with type 1 diabetes (T1D) (case group) were matched with 28 offspring of T1D fathers (control group) for the study of their leukocyte genome-wide DNA methylation profile (27,578 CpG sites, Human Methylation 27 BeadChip, Illumina Infinium). In a subset of 19 cases and 18 controls, we assessed renal vascular development by measuring Glomerular Filtration Rate (GFR) and Effective Renal Plasma Flow (ERPF) at baseline and during vasodilatation produced by amino acid infusion.

Results

Globally, DNA was under-methylated in cases vs. controls. Among the 87 CpG sites differently methylated, 74 sites were less methylated and 13 sites more methylated in cases vs. controls. None of these CpG sites were located on a gene known to be directly involved in kidney development and/or function. However, the gene encoding DNA methyltransferase 1 (DNMT1)—a key enzyme involved in gene expression during early development–was under-methylated in cases. The average methylation of the 74 under-methylated sites differently correlated with GFR in cases and controls.

Conclusion

Alterations in methylation profile imprinted by the hyperglycemic milieu of T1D mothers during fetal development may impact kidney function in adult offspring. The involved pathways seem to be a nonspecific imprinting process rather than specific to kidney development or function.     

Evaluation of microarray-based DNA methylation measurement using technical replicates: the Atherosclerosis Risk In Communities (ARIC) Study

Background

DNA methylation is a widely studied epigenetic phenomenon; alterations in methylation patterns influence human phenotypes and risk of disease. As part of the Atherosclerosis Risk in Communities (ARIC) study, the Illumina Infinium HumanMethylation450 (HM450) BeadChip was used to measure DNA methylation in peripheral blood obtained from ~3000 African American study participants. Over 480,000 cytosine-guanine (CpG) dinucleotide sites were surveyed on the HM450 BeadChip. To evaluate the impact of technical variation, 265 technical replicates from 130 participants were included in the study.

Results

For each CpG site, we calculated the intraclass correlation coefficient (ICC) to compare variation of methylation levels within- and between-replicate pairs, ranging between 0 and 1. We modeled the distribution of ICC as a mixture of censored or truncated normal and normal distributions using an EM algorithm. The CpG sites were clustered into low- and high-reliability groups, according to the calculated posterior probabilities. We also demonstrated the performance of this clustering when applied to a study of association between methylation levels and smoking status of individuals. For the CpG sites showing genome-wide significant association with smoking status, most (~96%) were seen from sites in the high reliability cluster.

Conclusions

We suggest that CpG sites with low ICC may be excluded from subsequent association analyses, or extra caution needs to be taken for associations at such sites.

Electronic supplementary material

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

High HIV Prevalence among Asylum Seekers Who Gave Birth in the Netherlands: A Nationwide Study Based on Antenatal HIV Tests

Objectives

Asylum seekers are considered to be a particularly vulnerable group with respect to HIV. Data on the HIV prevalence among asylum seekers, however, are scarce. The aim of this study is to map the HIV prevalence among asylum seekers who gave birth in The Netherlands.

Methods

We used a nationwide electronic medical records database from the community health services for asylum seekers (MOA). The study population consisted of 4,854 women and girls who delivered in asylum reception between 2000 and 2008. A unique electronic health data base was used and case allocation was based on ICPC-codes.

Results

The number of women and girls that was HIV positive during their last pregnancy was 80, of which 79 originated from sub-Saharan Africa. The prevalence for women from this region of origin (3.4%) was high compared to women from all other regions of origin (0.04%; OR = 90.2; 95%CI 12.5–648.8). The highest HIV prevalence rates were found for women from Rwanda (17.0%) and Cameroon (13.2%). HIV prevalence rates were higher among women who arrived in reception without partner (OR = 1.82; 95%CI 0.75–4.44) and unaccompanied minors (OR = 2.59; 95%CI 0.79–8.49), compared to women who arrived in reception with partner.

Conclusions

We conclude that, among asylum-seeking women from sub-Saharan Africa giving birth in The Netherlands, the HIV prevalence is high compared to the host population. For women from other regions of origin, the prevalence is at the same level as in the host population. The high HIV prevalence underlines the importance of preventive interventions and voluntary HIV testing for sub-Saharan African asylum seekers as from shortly after arrival.     

Factors to preserve CpG-rich sequences in methylated CpG islands

Background

Mammalian CpG islands (CGIs) normally escape DNA methylation in all adult tissues and developmental stages. However, in our previous study we unexpectedly identified many methylated CGIs in human peripheral blood leukocytes. Methylated CpG dinucleotides convert to TpG dinucleotides through deaminization of their cytosine bases more frequently than hypomethylated CpG dinucleotides. Therefore, we wondered how methylated CGIs in germline or non-germline cells maintain their CpG-rich sequences. It is known that events such as germline hypomethylation, CpG selection, biased gene conversion (BGC), and frequent CpG fixation can contribute to the maintenance of CpG-rich sequences in methylated CGIs in germline or non-germline cells. However, it has not been investigated which of the processes maintain CpG-rich sequences of methylated CGIs in each genomic position.

Results

In this study, we comprehensively examined the contribution of the processes described above to the maintenance of CpG-rich sequences in methylated CGIs in germline and non-germline cells which were classified by genomic positions. Approximately 60–80% of CGIs with high methylation in H1 cell line (H1-HM) in all the genomic positions showed a low average CpG → TpG/CpA substitution rate. In contrast, fewer than half the numbers of CGIs with H1-HM in all the genomic positions showed a low average CpG → TpG/CpA substitution rate and low levels of methylation in sperm cells (SPM-LM). Furthermore, a small fraction of CGIs with a low average CpG → TpG/CpA substitution rate and high levels of methylation in sperm cells (SPM-HM) showed CpG selection.On the other hand, independent of the positions in genes, most CGIs with SPM-HM showed a slightly higher average TpG/CpA → CpG substitution rate compared with those with SPM-LM.

Conclusions

Relatively high numbers (approximately 60–80%) of CGIs with H1-HM in all the genomic positions preserve their CpG-rich sequences by a low CpG → TpG/CpA substitution rate caused mainly by their SPM-LM, and for those with SPM-HM partly by CpG selection and TpG/CpA → CpG fixation. BGC has little contribution to the maintenance of CpG-rich sequences of CGIs with SPM-HM which were classified by genomic positions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1286-x) contains supplementary material, which is available to authorized users.     

Motif depletion in bacteriophages infecting hosts with CRISPR systems

Background

CRISPR is a microbial immune system likely to be involved in host-parasite coevolution. It functions using target sequences encoded by the bacterial genome, which interfere with invading nucleic acids using a homology-dependent system. The system also requires protospacer associated motifs (PAMs), short motifs close to the target sequence that are required for interference in CRISPR types I and II. Here, we investigate whether PAMs are depleted in phage genomes due to selection pressure to escape recognition.

Results

To this end, we analyzed two data sets. Phages infecting all bacterial hosts were analyzed first, followed by a detailed analysis of phages infecting the genus Streptococcus, where PAMs are best understood. We use two different measures of motif underrepresentation that control for codon bias and the frequency of submotifs. We compare phages infecting species with a particular CRISPR type to those infecting species without that type. Since only known PAMs were investigated, the analysis is restricted to CRISPR types I-C and I-E and in Streptococcus to types I-C and II. We found evidence for PAM depletion in Streptococcus phages infecting hosts with CRISPR type I-C, in Vibrio phages infecting hosts with CRISPR type I-E and in Streptococcus thermopilus phages infecting hosts with type II-A, known as CRISPR3.

Conclusions

The observed motif depletion in phages with hosts having CRISPR can be attributed to selection rather than to mutational bias, as mutational bias should affect the phages of all hosts. This observation implies that the CRISPR system has been efficient in the groups discussed here.

Electronic supplementary material

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