Recombination in hepatitis C virus genotype 1 evaluated by phylogenetic and population-genetic methods.

Although hepatitis C virus (HCV) is a major cause of viral hepatitis and hepatocellular carcinoma, many aspects of its evolution remain poorly understood. Relevant to its evolution and the development of antiviral drug resistance is the role of recombination in HCV, which has not been resolved using phylogenetic tests. In line with previous studies, we found no strong support for a role of recombination in the dominant subtypes 1A and 1B using phylogenetic approaches. In contrast, signatures of gene conversion were abundant if a population recombination model, which takes into account diversity within and between groups, was used (9676 gene conversion signatures between the genomes of subtypes 1A and 1B and 170 between the NS5A regions of subtypes 1A and 1B and the minor subtypes 1c-1g). The gene conversion signatures coincided with a striking lack of diagnostically informative sites between subtypes and a large number of shared mutations between complete subtype 1A and 1B genomes (0.76 and 62.2 % of nucleotide sites, respectively). Maximum-likelihood trees revealed significant topological incongruence among conserved PFAM domains and genome regions targeted by diagnostic assays, which underpins a major role for recombination. The same results were obtained with smaller numbers of genomes and with only synonymous sites. Topological concordance increased only marginally if larger genome regions were compared. The level of recombination in HCV subtype 1, which is probably significantly higher than can currently be measured, also illustrates the complexity of designing diagnostic assays based on the unusual patterns of genomic diversity of HCV.

Population genetic tests suggest that the epidemiologies of JCV and BKV are strikingly different.

The JCV and BKV viruses have been used as markers for the study of human evolution by assuming that these viruses coevolved with their host. However, it is currently unclear whether the details of the population expansion of these viruses and humans agree. To study this in more detail, large numbers of complete genomes were used for population genetic tests to detect evidence for population expansion. Relative to the neutral expectation of no selective forces and no demographic changes, the JCV data set contained a striking excess of synonymous and non-synonymous mutations that occur only once in the data set. The same was found for non-synonymous mutations of BKV, but not at all for synonymous mutations of BKV. The different frequency spectra of mutations in JCV and BKV do not result from the inclusion of patients with clinical symptoms associated with BKV and JCV, such as nephropathy or progressive multifocal leucoencefalopathy, nor from the different numbers of genomes available for JCV and BKV. Instead, the distribution of unique mutations and population genetic models that use older mutation classes indicate a striking difference of the historical demographies of JCV and BKV with only the former virus exhibiting the evidence of demographic expansion. Our analyses expand on recent population genetic analyses that document a global population expansion of JCV by taking into account the impact of deleterious mutations and by comparing both human viruses. The striking difference between the demographics of BKV and JCV suggests that important aspects of their epidemiology remain to be discovered.

The effect of primer-template mismatches on the detection and quantification of nucleic acids using the 5' nuclease assay.

Real-time polymerase chain reaction (PCR) is the current method of choice for detection and quantification of nucleic acids, especially for molecular diagnostics. Complementarity between primers and template is often crucial for PCR applications, as mismatches can severely reduce priming efficiency. However, little quantitative data on the effect of these mismatches is available. We quantitatively investigated the effects of primer-template mismatches within the 3'-end primer region on real-time PCR using the 5'-nuclease assay. Our results show that single mismatches instigate a broad variety of effects, ranging from minor (<1.5 cycle threshold, eg, A-C, C-A, T-G, G-T) to severe impact (>7.0 cycle threshold, eg, A-A, G-A, A-G, C-C) on PCR amplification. A clear relationship between specific mismatch types, position, and impact was found, which remained consistent for DNA versus RNA amplifications and Taq/Moloney murine leukemia virus versus rTth based amplifications. The overall size of the impact among the various master mixes used differed substantially (up to sevenfold), and for certain master mixes a reverse or forward primer-specific impact was observed, emphasizing the importance of the experimental conditions used. Taken together these data suggest that mismatch impact follows a consistent pattern and enabled us to formulate several guidelines for predicting primer-template mismatch behavior when using specific 5-nuclease assay master mixes. Our study provides novel insight into mismatch behavior and should allow for more optimized development of real-time PCR assays involving primer-template mismatches.

Mapping QTL influencing gastrointestinal nematode burden in Dutch Holstein-Friesian dairy cattle.

BACKGROUND: Parasitic gastroenteritis caused by nematodes is only second to mastitis in terms of health costs to dairy farmers in developed countries. Sustainable control strategies complementing anthelmintics are desired, including selective breeding for enhanced resistance. RESULTS AND CONCLUSION: To quantify and characterize the genetic contribution to variation in resistance to gastro-intestinal parasites, we measured the heritability of faecal egg and larval counts in the Dutch Holstein-Friesian dairy cattle population. The heritability of faecal egg counts ranged from 7 to 21% and was generally higher than for larval counts. We performed a whole genome scan in 12 paternal half-daughter groups for a total of 768 cows, corresponding to the approximately 10% most and least infected daughters within each family (selective genotyping). Two genome-wide significant QTL were identified in an across-family analysis, respectively on chromosomes 9 and 19, coinciding with previous findings in orthologous chromosomal regions in sheep. We identified six more suggestive QTL by within-family analysis. An additional 73 informative SNPs were genotyped on chromosome 19 and the ensuing high density map used in a variance component approach to simultaneously exploit linkage and linkage disequilibrium in an initial inconclusive attempt to refine the QTL map position.

An improved approach to identify epidemiological and phylogenetic transmission pairs of source and contact tracing of hepatitis B.

The transmission of infectious diseases can be traced using epidemiological and molecular information. In the current study, the congruence was assessed between sequence data of the hepatitis B virus (HBV) and epidemiological information resulting from source and contact tracing of patients seen at the Municipal Public Health Service in Rotterdam between 2002 and 2005. HBV genotypes A-G were present in 62 acute and 334 chronic HBV patients. At the sequence level, the identical sequences of members of epidemiological transmission pairs and the rarity of such pairs provided strong support for correctness of the hypothesized transmission routes. The molecular support for epidemiological transmission pairs derived from source and contact tracing was further assessed by using topological constraints in parsimony analyses in agreement with epidemiological information, and by taking the presence of polymorphic sites of HBV within patients into account. This, in principle, allows mutations in epidemiological clusters. Of 22 epidemiological clusters, six could be refuted, four clusters received support from the molecular analysis, and support for the remaining twelve clusters was ambiguous. Two of the four epidemiological pairs that received molecular support had diverged (by 3 and 15 mutations). These results show that levels of divergence cannot be used simply as an indicator of the likelihood that groups of sequences constitute transmission pairs. Instead, to confirm or refute transmission pairs, it is necessary to assess the likelihood of a common origin of HBV variants in epidemiologically defined transmission groups relative to the HBV diversity in the local community.

Microbial diversity--insights from population genetics.

Although many environmental microbial populations are large and genetically diverse, both the level of diversity and the extent to which it is ecologically relevant remain enigmatic. Because the effective (or long-term) population size, N(e), is one of the parameters that determines population genetic diversity, tests and simulations that assume selectively neutral mutations may help to identify the processes that have shaped microbial diversity. Using ecologically important genes, tests of selective neutrality suggest that adaptive as well as non-adaptive types of selection act and that departure from neutrality may be widespread or restricted to small groups of genotypes. Population genetic simulations using population sizes between 10(3) and 10(7) suggest extremely high levels of microbial diversity in environments that sustain large populations. However, census and effective population sizes may differ considerably, and because we know nothing of the evolutionary history of environmental microbial populations, we also have no idea what N(e) of environmental populations is. On the one hand, this reflects our ignorance of the microbial world. On the other hand, the tests and simulations illustrate interactions between microbial diversity and microbial population genetics that should inform our thinking in microbial ecology. Because of the different views on microbial diversity across these disciplines, such interactions are crucial if we are to understand the role of genes in microbial communities.

Positively selected codons in immune-exposed loops of the vaccine candidate OMP-P1 of Haemophilus influenzae.

The high levels of variation in surface epitopes can be considered as an evolutionary hallmark of immune selection. New computational tools enable analysis of this variation by identifying codons that exhibit high rates of amino acid changes relative to the synonymous substitution rate. In the outer membrane protein P1 of Haemophilus influenzae, a vaccine candidate for nontypeable strains, we identified four codons with this attribute in domains that did not correspond to known or assumed B- and T-cell epitopes of OMP-P1. These codons flank hypervariable domains and do not appear to be false positives as judged from parsimony and maximum likelihood analyses. Some closely spaced positively selected codons have been previously considered part of a transmembrane domain, which would render this region unsuited for inclusion in a vaccine. Secondary structure analysis, three-dimensional structural database searches, and homology modeling using FadL of E. coli as a structural homologue, however, revealed that all positively selected codons are located in or near extracellular looping domains. The spacing and level of diversity of these positively selected and exposed codons in OMP-P1 suggest that vaccine targets based on these and conserved flanking residues may provide broad coverage in H. influenzae.

A simple, robust and semi-automated parasite egg isolation protocol.

Large-scale parasite quantification is required for improving our understanding of the epidemiology and genetics of host-parasite interactions. We describe a protocol that uses a low-density salt solution for flotation and centrifugation of nematode eggs. Subsequently, sucrose flotation and precipitation are used to obtain clear egg preparations. Most traditional quantification protocols such as the McMaster technique are unsuited for the standardized processing of large numbers of samples and the analysis of large amounts of feces per sample. Consequently, they are suited only for small-scale surveys. Our protocol, which can be used to analyze up to 6 g of feces, results in clear egg preparations that are concentrated in wells of a microtiter plate and that are suited for digital recording and automated counting. Starting from a fecal suspension in the first flotation solution to a digital recording requires approximately 40 min per 24 samples.

Positive selection on transposase genes of insertion sequences in the Crocosphaera watsonii genome.

Insertion sequences (ISs) are mobile elements that are commonly found in bacterial genomes. Here, the structural and functional diversity of these mobile elements in the genome of the cyanobacterium Crocosphaera watsonii WH8501 is analyzed. The number, distribution, and diversity of nucleotide and amino acid stretches with similarity to the transposase gene of this IS family suggested that this genome harbors many functional as well as truncated IS fragments. The selection pressure acting on full-length transposase open reading frames of these ISs suggested (i) the occurrence of positive selection and (ii) the presence of one or more positively selected codons. These results were obtained using three data sets of transposase genes from the same IS family that were collected based on the level of amino acid similarity, the presence of an inverted repeat, and the number of sequences in the data sets. Neither recombination nor ribosomal frameshifting, which may interfere with the selection analyses, appeared to be important forces in the transposase gene family. Some positively selected codons were located in a conserved domain, suggesting that these residues are functionally important. The finding that this type of selection acts on IS-carried genes is intriguing, because although ISs have been associated with the adaptation of the bacterial host to new environments, this has typically been attributed to transposition or transformation, thus involving different genomic locations. Intragenic adaptation of IS-carried genes identified here may constitute a novel mechanism associated with bacterial diversification and adaptation.

Selection on protein-coding genes of natural cyanobacterial populations.

We examined the distribution of synonymous and non-synonymous changes in 12 protein-coding genes of natural populations of cyanobacteria to infer changes in gene functionality. By comparing mutation distributions within and across species using the McDonald-Kreitman test, we found data sets to contain evidence for purifying selection (hetR of Trichodesmium, nifH of Cylindrospermopsis raceborskii and rpoC1 of Anabaena lemmermannii) and positive selection (kaiC of Microcoleus chthonoplastes and rbcX of Anabaena and Aphanizomenon sp.). Other genes from the same set of clonal isolates (petB and rbcL in M. chthonoplastes and Anabaena/Aphanizomenon, respectively) did not harbour evidence for either form of selection. The results of branch models of codon evolution agreed fully with the results of the McDonald-Kreitman test in terms of significance and absolute value of the dN/dS estimates. The high frequency of gene-specific mutation patterns and their association with branches that separate closely related cyanobacterial genera suggest that evolutionary tests are suited to uncover gene-specific selective differentiation in cyanobacterial genomes. At the same time, given the lack of information about the history of cyanobacteria, analysis of larger numbers of protein-coding genes of clonal cyanobacterial isolates will produce more detailed pictures of the effects of natural selection.

Purifying selection and demographic expansion affect sequence diversity of the ligand-binding domain of a glutamate-gated chloride channel gene of Haemonchus placei.

Ninety-five genomic sequences of the ligand-binding domain of glutamate-gated chloride channel genes of three populations of the parasitic nematode H. placei were evaluated for patterns of diversity, demography, and selection. These genes code for subunits of ion channels, which are involved in the mode of action of the most commonly used antiparasitic drugs, the macrocyclic lactones. An extremely high frequency of unique segregating sites in exons and introns was observed, with significantly negative neutrality tests in each population for noncoding, synonymous, and nonsynonymous sites. Several tests indicated that support for balancing selection, positive selection, and hitchhiking was lacking. McDonald-Kreitman tests using H. contortus or C. elegans as an outgroup revealed an extreme excess of replacement polymorphism, consistent with weak purifying selection. Although these tests agree that negative selection may explain the excess of replacement changes, an alternative interpretation is required for the significantly negative Fu and Li's D statistics based on silent and noncoding sites. These include homogeneous forces such as background selection and demographic expansion. The lack of population subdivision and the negative values of Tajima's D for this outbreeding parasitic nematode render background selection less likely than demographic expansion. Comparison of D statistics based on different site types using neutral coalescent simulations supported this interpretation. Although this statistic was more negative for nonsynonymous sites than for synonymous sites, most comparisons of the D statistic were not significantly different between mutation classes. A few significant site comparisons were also consistent with demographic expansion, because the observed test statistic ( D(neutral) - D(selected)) were low relative to the neutral expectations. Finally, previous mitochondrial studies also identified a demographic expansion of this parasitic nematode species, which lends further support to a scenario involving both demographic and purifying forces in the ligand-binding domain of H. placei.

Technical variability and required sample size of helminth egg isolation procedures.

Measurements of parasite load are often very variable. This implies that little confidence can be attached to single measurements of parasite numbers and egg concentrations, and that many measurements are required for the detection of differences between groups of hosts or parasites. For studies that aim to detect these differences, it is important to increase the precision (closeness of repeated measures to each other) of parasite numbers, because it determines the number of samples that is needed to find significant differences among groups. In this study, sample sizes required to detect group differences were estimated using nematode egg counts of faecal samples of dairy cattle. They were found to be much lower for a centrifugation technique than for the widely used McMaster technique in replicate samples, in spite of a generally similar mean FEC. For example, the sample size required to detect FEC differences between groups of 10, 50, and 250 eggs per gram (EPG) were 46, 25, and 27 for the McMaster technique and 8, 5, and 12 for the SSF method, respectively. Interestingly, sample sizes required for faeces with a relatively high egg concentration (approximately 1000 EPG) were also considerably lower than for the McMaster technique in spite of a higher mean EPG of the latter method. This implies that technical variation can be reduced considerably by simple methods of egg isolation. Given that the range of egg concentration is similar for a number of nematodes of livestock and human helminths, a reduction of technical error will aid studies with many group comparisons such as vaccination strategies against parasites with typically low FECs and studies of the genetics of host resistance. It may also lead to improved guidelines for measures related to public health.

Detection of genetically divergent clone mates in apomictic dandelions.

This study aims to identify genetically diverged clone mates in apomictic dandelions. Clone mates are defined as individuals that may have diverged as a result of mutation accumulation and that have undergone only clonal reproduction since their most recent common ancestor. Based on distinctive morphology and an aberrant and rare chloroplast haplotype, northwest European individuals of Taraxacum section Naevosa are well suited for the detection of clonal lineages in which mutation has occurred. In the case of strictly clonal reproduction, nuclear genetic variability was expected to be hierarchically organized. Nucleotide polymorphisms in internal transcribed spacer (ITS) sequences, however, were incompatible with a clonal structure of the Norwegian individuals, probably due to persistent ancestral polymorphisms that pre-date the origin of the Naevosa clone. This interpretation is supported by the presence of ITS variants in section Naevosa that were also found in distantly related dandelions. In contrast to the ITS sequence data, amplified fragment length polymorphisms (AFLPs), isozymes and microsatellites strongly supported the contention of prolonged clonal reproduction and mutation accumulation in Norwegian Naevosa. Because these markers are generally considered to be more variable and more rapidly evolving than ITS sequences, mutations in these markers probably evolved after the origin of the clone. Within the Norwegian clone, a surprising number of markers distinguished the clone mates. As a consequence, incorporation of mutation in the detection of clone mates is anticipated to have a big impact on estimates of size, geographical range and age of clones as well as on experimental designs of studies of clonal plants.