Sequence Variation in the Gene Encoding the Nonstructural 3 Protein of Hepatitis C Virus: Evidence for Immune Selection

To determine whether the persistent nature of hepatitis C infection is related to the emergence of antigenic variants driven by immune selection, we examined the sequence heterogeneity in a portion of the hepatitis C virus (HCV) nonstructural 3 (NS3) gene of a patient infected over the course of more than 2 years. By PCR amplification, cloning, and sequencing, we observed several variable and conserved regions in the NS3 segment of the HCV genome. All variable regions had higher ratios of nonsynonymous/synonymous mutations and encompassed immunodominant epitopes, and their locations were not essential to maintain the known function of HCV RNA helicase. In contrast, the regions that are critical for HCV RNA helicase activity were found to be conserved with lower heterogeneity or lower ratios of nonsynonymous/synonymous mutations, and none except one of these regions was encoded within immunodominant epitopes. Our results are consistent with immune selection of viral variants at the epitope and molecular levels that may enable HCV to evade host defenses over time. Plotting the relatedness of sequence variants revealed a star topology suggesting that a wild-type HCV sequence is maintained, unlike HIV. Received: 2 November 2000 / Accepted: 1 October 2001     

Genomic repeats, genome plasticity and the dynamics of Mycoplasma evolution

Mycoplasmas evolved by a drastic reduction in genome size, but their genomes contain numerous repeated sequences with important roles in their evolution. We have established a bioinformatic strategy to detect the major recombination hot-spots in the genomes of Mycoplasma pneumoniae, Mycoplasma genitalium, Ureaplasma urealyticum and Mycoplasma pulmonis. This allowed the identification of large numbers of potentially variable regions, as well as a comparison of the relative recombination potentials of different genomic regions. Different trends are perceptible among mycoplasmas, probably due to different functional and structural constraints. The largest potential for illegitimate recombination in M.pulmonis is found at the vsa locus and its comparison in two different strains reveals numerous changes since divergence. On the other hand, the main M.pneumoniae and M.genitalium adhesins rely on large distant repeats and, hence, homologous recombination for variation. However, the relation between the existence of repeats and antigenic variation is not necessarily straightforward, since repeats of P1 adhesin were found to be anti-correlated with epitopes recognized by patient antibodies. These different strategies have important consequences for the structures of genomes, since large distant repeats correlate well with the major chromosomal rearrangements. Probably to avoid such events, mycoplasmas strongly avoid inverse repeats, in comparison to co-oriented repeats.     

Evolutionary Adaptation of the Amino Acid and Codon Usage of the Mosquito Sodium Channel following Insecticide Selection in the Field Mosquitoes

Target site insensitivity resulting from point mutations within the voltage-gated sodium channel of the insect nervous system is known to be of primary importance in the development of resistance to pyrethroid insecticides. This study shifts current research paradigms by conducting, for the first time, a global analysis of all the naturally occurring mutations, both nonsynonymous and synonymous mutations, as well as mutation combinations in the entire mosquito sodium channel of Culex quinquefasciatus and analyzing their evolutionary and heritable feature and roles in insecticide resistance. Through a systematic analysis of comparing nucleotide polymorphisms in the entire sodium channel cDNAs of individuals between susceptible and resistant mosquito strains, between field parental mosquitoes and their permethrin selected offspring, and among different mosquito groups categorized by their levels of tolerance to specific permethrin concentrations within and among the mosquito strains of the field parental strains and their permethrin selected offspring, 3 nonsynonymous (A¹⁰⁹S, L⁹⁸²F, and W¹⁵⁷³R) and 6 synonymous (L⁸⁵², G⁸⁹¹, A¹²⁴¹, D¹²⁴⁵, P¹²⁴⁹, and G¹⁷³³) mutations were identified. The co-existence of all 9 mutations, both nonsynonymous and synonymous, and their homozygousity were found to be important factors for high levels of resistance. Our study, for the first time, provide a strong case demonstrating the co-existence of both nonsynonymous and synonymous mutations in the sodium channel of resistant mosquitoes in response to insecticide resistance and the inheritance of these mutations in the offspring of field mosquito strains following insecticide selection.     

Quantifying Selection against Synonymous Mutations in HIV-1 env Evolution

Intrapatient evolution of human immunodeficiency virus type 1 (HIV-1) is driven by the adaptive immune system resulting in rapid change of HIV-1 proteins. When cytotoxic CD8⁺ T cells or neutralizing antibodies target a new epitope, the virus often escapes via nonsynonymous mutations that impair recognition. Synonymous mutations do not affect this interplay and are often assumed to be neutral. We test this assumption by tracking synonymous mutations in longitudinal intrapatient data from the C2-V5 part of the env gene. We find that most synonymous variants are lost even though they often reach high frequencies in the viral population, suggesting a cost to the virus. Using published data from SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) assays, we find that synonymous mutations that disrupt base pairs in RNA stems flanking the variable loops of gp120 are more likely to be lost than other synonymous changes: these RNA hairpins might be important for HIV-1. Computational modeling indicates that, to be consistent with the data, a large fraction of synonymous mutations in this genomic region need to be deleterious with a cost on the order of 0.002 per day. This weak selection against synonymous substitutions does not result in a strong pattern of conservation in cross-sectional data but slows down the rate of evolution considerably. Our findings are consistent with the notion that large-scale patterns of RNA structure are functionally relevant, whereas the precise base pairing pattern is not.     

Persistence of Mycoplasma genitalium following azithromycin therapy

Background

To determine clinical outcomes and cure rates for M.genitalium genital infection in men and women following azithromycin 1 g.

Methodology

Patients attending Melbourne Sexual Health Centre between March 2005 and November 2007 with urethritis/epididymitis, cervicitis/pelvic inflammatory disease and sexual contacts of M.genitalium were tested for M.genitalium by polymerase chain reaction (PCR). M.genitalium-infection was treated with 1 g of azithromycin and a test-of-cure (toc) was performed one month post-azithromycin. Response to azithromycin, and response to moxifloxacin (400 mg daily for 10 days) in individuals with persistent infection post-azithromycin, was determined.

Principal Findings

Of 1538 males and 313 females tested, 161 males (11%) and 30 females (10%) were infected with M.genitalium. A toc was available on 131 (69%) infected individuals (median = 36 days [range 12-373]). Of 120 individuals prescribed azithromycin only pre-toc, M.genitalium was eradicated in 101 (84%, 95% confidence intervals [CI]: 77–90%) and persisted in 19 (16%, 95% CI: 10–23%). Eleven individuals with persistent infection (9%, 95% CI: 5–15%) had no risk of reinfection from untreated-partners, while eight (7%, 95% CI: 3–12%) may have been at risk of reinfection from doxycycline-treated or untreated-partners. Moxifloxacin was effective in eradicating persistent infection in all cases not responding to azithromycin. Patients with persistent-M.genitalium were more likely to experience persistent symptoms (91%), compared to patients in whom M.genitalium was eradicated (17%), p<0.0001.

Conclusion

Use of azithromycin 1 g in M.genitalium-infected patients was associated with unacceptable rates of persistent infection, which was eradicated with moxifloxacin. These findings highlight the importance of follow-up in M.genitalium-infected patients prescribed azithromycin, and the need to monitor for the development of resistance. Research to determine optimal first and second-line therapeutic agents for M.genitalium is needed.     

Novel point mutations in the ERG11 gene in clinical isolates of azole resistant Candida species

The azoles are the class of medications most commonly used to fight infections caused by Candida sp. Typically, resistance can be attributed to mutations in ERG11 gene (CYP51) which encodes the cytochrome P450 14α-demethylase, the primary target for the activity of azoles. The objective of this study was to identify mutations in the coding region of theERG11 gene in clinical isolates of Candidaspecies known to be resistant to azoles. We identified three new synonymous mutations in the ERG11 gene in the isolates of Candida glabrata (C108G, C423T and A1581G) and two new nonsynonymous mutations in the isolates of Candida krusei - A497C (Y166S) and G1570A (G524R). The functional consequence of these nonsynonymous mutations was predicted using evolutionary conservation scores. The G524R mutation did not have effect on 14α-demethylase functionality, while the Y166S mutation was found to affect the enzyme. This observation suggests a possible link between the mutation and dose-dependent sensitivity to voriconazole in the clinical isolate of C. krusei. Although the presence of the Y166S in phenotype of reduced azole sensitivity observed in isolate C. kruseidemands investigation, it might contribute to the search of new therapeutic agents against resistant Candida isolates.     

Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans,but not in chimpanzees.

To test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution, nucleotide sequences were determined for the 1041 bp of the NADH dehydrogenase subunit 2 (ND2) gene in 20 geographically diverse humans and 20 common chimpanzees. Contingency tests of neutrality were performed using four mutational categories for the ND2 molecule: synonymous and nonsynonymous mutations in the transmembrane regions, and synonymous and nonsynonymous mutations in the surface regions. The following three topological mutational categories were also used: intraspecific tips, intraspecific interiors, and interspecific fixed differences. The analyses reveal a significantly greater number of nonsynonymous polymorphisms within human transmembrane regions than expected based on interspecific comparisons, and they are inconsistent with a neutral equilibrium model. This pattern of excess nonsynonymous polymorphism is not seen within chimpanzees. Statistical tests of neutrality, such as TAJIMA''s D test, and the D and F tests proposed by FU and LI, indicate an excess of low frequency polymorphisms in the human data, but not in the chimpanzee data. This is consistent with recent directional selection, a population bottleneck or background selection of slightly deleterious mutations in human mtDNA samples. The analyses further support the idea that mitochondrial genome evolution is governed by selective forces that have the potential to affect its use as a "neutral" marker in evolutionary and population genetic studies.     

Assessing in vivo mutation frequencies and creating a high-resolution genome-wide map of fitness costs of Hepatitis C virus

Like many viruses, Hepatitis C Virus (HCV) has a high mutation rate, which helps the virus adapt quickly, but mutations come with fitness costs. Fitness costs can be studied by different approaches, such as experimental or frequency-based approaches. The frequency-based approach is particularly useful to estimate in vivo fitness costs, but this approach works best with deep sequencing data from many hosts are. In this study, we applied the frequency-based approach to a large dataset of 195 patients and estimated the fitness costs of mutations at 7957 sites along the HCV genome. We used beta regression and random forest models to better understand how different factors influenced fitness costs. Our results revealed that costs of nonsynonymous mutations were three times higher than those of synonymous mutations, and mutations at nucleotides A or T had higher costs than those at C or G. Genome location had a modest effect, with lower costs for mutations in HVR1 and higher costs for mutations in Core and NS5B. Resistance mutations were, on average, costlier than other mutations. Our results show that in vivo fitness costs of mutations can be site and virus specific, reinforcing the utility of constructing in vivo fitness cost maps of viral genomes.     

Identification of polycystin-1 and Gα12 binding regions necessary for regulation of apoptosis

Most patients with autosomal dominant polycystic kidney disease (ADPKD) harbor mutations in PKD1, the gene for polycystin-1 (PC1), a transmembrane protein with a cytoplasmic C-terminus that interacts with numerous signaling molecules, including Gα12. The functions of PC1 and the mechanisms of cyst development leading to renal failure are complex. Recently, we reported that PC1 expression levels modulate activity of Gα12-stimulated apoptosis (Yu et al., J. Biol. Chem. 2010 285(14):10243-51). Herein, a mutational analysis of Gα12 and PC1 was undertaken to identify regions required for their interaction and ability to modulate apoptosis. A set of Gα12 mutations with systematic replacement of six amino acids with NAAIRS was tested for binding to the PC1 C-terminus in GST pulldowns. Additionally, a series of deletions within the PC1 C-terminus was examined for binding to Gα12. We identified 3 NAAIRS substitutions in Gα12 that completely abrogated binding, and identified a previously described 74 amino acid Gαi/o binding domain in the PC1 C-terminus as necessary for Gα12 interaction. The functional consequences of uncoupling PC1/Gα12 binding were studied in apoptosis assays utilizing HEK293 cells with inducible PC1 overexpression. Gα12 mutants deficient in PC1 binding were refractory to PC1 inhibition of Gα12-stimulated apoptosis. Likewise, deletion of the Gα12-interacting sequence from the PC1 cytoplasmic domain abrogated its inhibition of Gα12-stimulated apoptosis. Based on the crystal structure of Gα12, the PC1 interaction sites are likely to reside on exposed regions within the G protein helical domain. These structural details should facilitate the design of reagents to uncouple PC1/Gα12 signaling in ADPKD.     

Production and characterization of recombinant P1 adhesin essential for adhesion,gliding, and antigenic variation in the human pathogenic bacterium,Mycoplasma pneumoniae

Mycoplasma pneumoniae forms an attachment organelle at one cell pole, binds to the host cell surface, and glides via a unique mechanism. A 170-kDa protein, P1 adhesin, present on the organelle surface plays a critical role in the binding and gliding process. In this study, we obtained a recombinant P1 adhesin comprising 1476 amino acid residues, excluding the C-terminal domain of 109 amino acids that carried the transmembrane segment, that were fused to additional 17 amino acid residues carrying a hexa-histidine (6?×?His) tag using an Escherichia coli expression system. The recombinant protein showed solubility, and chirality in circular dichroism (CD). The results of analytical gel filtration, ultracentrifugation, negative-staining electron microscopy, and small-angle X-ray scattering (SAXS) showed that the recombinant protein exists in a monomeric form with a uniformly folded structure. SAXS analysis suggested the presence of a compact and ellipsoidal structure rather than random or molten globule-like conformation. Structure model based on SAXS results fitted well with the corresponding structure obtained with cryo-electron tomography from a closely related species, M. genitalium. This recombinant protein may be useful for structural and functional studies as well as for the preparation of antibodies for medical applications.     

Evaluation of accumulation of hepatitis C virus mutations in a chronically infected chimpanzee: comparison of the core, E1, HVR1, and NS5b regions

Four hepatitis C virus genome regions (the core, E1, HVR1, and NS5b) were amplified and sequenced from yearly samples obtained from a chronically infected chimpanzee over a 12-year span. Nucleotide substitutions were found to accumulate in the core, E1, and HVR1 regions during the course of chronic infection; substitutions within the NS5b region were not detected for the first 8 years and were found to be minimal during the last 4 years. The rate of accumulation of mutations in the core and E1 regions, based on a direct comparison between the first 1979 sequence and the last 1990 sequence, was 1.120 x 10(-3), while phylogenetic ancestral comparison using the 12 yearly sequences showed a rate of 0.816 x 10(-3) bases per site per year. Temporal evaluation of the sequences revealed that there appeared to be periods in which substitutions accumulated and became fixed, followed by periods with relative stasis or random substitutions that did not persist. Synonymous and nonsynonymous substitutions within the core, E1, and HVR1 regions were also analyzed. In the core and E1 regions, synonymous substitutions predominated and gradually increased over time. However, within the HVR1 region, nonsynonymous substitutions predominated but gradually decreased over time.     

Draft Genome Sequences of Four Axenic Mycoplasma genitalium Strains Isolated from Denmark,Japan, and Australia

The DNA genome of Mycoplasma genitalium currently represents the smallest of all known human bacterial pathogens. Despite their clinical importance in sexually transmitted infection and relevance as model bacterial pathogens, genomic diversity among M. genitalium strains worldwide is unknown. Herein we present the complete draft genome sequences of four geographically diverse strains of M. genitalium.     

Unbiased estimation of the rates of synonymous and nonsynonymous substitution

Summary The current convention in estimating the number of substitutions per synonymous site (K _S ) and per nonsynonymous site (K _A ) between two protein-coding genes is to count each twofold degenerate site as one-third synonymous and two-thirds nonsynonymous because one of the three possible changes at such a site is synonymous and the other two are nonsynonymous. This counting rule can considerably overestimate theK _S value because transitional mutations tend to occur more often than transversional mutations and because most transitional mutations at twofold degenerate sites are synonymous. A new method that gives unbiased estimates is proposed. An application of the new and the old method to 14 pairs of mouse and rat genes shows that the new method gives aK _S value very close to the number of substitutions per fourfold degenerate site whereas the old method gives a value 30% higher. Both methods give aK _A value close to the number of substitutions per nondegenerate site.     

Gene conversion generates hypervariability at the variable regions of kallikreins and their inhibitors.

The two mechanisms for generating hypervariability at the reactive center of serine proteases and their inhibitors are gene conversion followed by natural selection and natural selection for point mutation. One way to clarify the effects of these two mechanisms is to calculate separately the number of nonsynonymous substitutions and that of synonymous substitutions at the variable regions and at the conserved regions. Our data analysis shows that not only the number of nonsynonymous substitutions but also the number of synonymous substitutions at the variable regions exceed the corresponding numbers at the conserved regions. Thus gene conversion has provided needed variability at the variable regions of serine proteases and their inhibitors. Natural selection has helped perpetuate such variability.     

Defective pre-mRNA splicing in PKD1 due to presumed missense and synonymous mutations causing autosomal dominant polycystic disease

Autosomal dominant polycystic kidney disease is the most common human monogenic disorder and is caused by mutations in the PKD1 or PKD2 genes. Most patients with the disease present mutations in PKD1, and a considerable number of these alterations are single base substitutions within the coding sequence that are usually predicted to lead to missense or synonymous mutations. There is growing evidence that some of these mutations can be detrimental by affecting the pre-mRNA splicing process. The aim of our study was to test PKD1 mutations, described as missense or synonymous in the literature or databases, for their effects on exon inclusion. Bioinformatics tools were used to select mutations with a potential effect on pre-mRNA splicing. Mutations were experimentally tested using minigene assays. Exons and adjacent intronic sequences were PCR-amplified and cloned in the splicing reporter minigene, and selected mutations were introduced by site-directed mutagenesis. Minigenes were transfected into kidney derived cell lines. RNA from cultured cells was analyzed by RT-PCR and DNA sequencing. Analysis of thirty-three PKD1 exonic mutations revealed three mutations that induce splicing defects. The substitution c.11156G > A, previously predicted as missense mutation p.R3719Q, abolished the donor splice site of intron 38 and resulted in the incorporation of exon 38 with 117 bp of intron 38 and skipping of exon 39. Two synonymous variants, c.327A > T (p.G109G) and c.11257C > A (p.R3753R), generated strong donor splice sites within exons 3 and 39 respectively, resulting in incorporation of incomplete exons. These three nucleotide substitutions represent the first PKD1 exonic mutations that induce aberrant mRNAs. Our results strengthen the importance to evaluate the consequences of presumed missense and synonymous mutations at the mRNA level.     

How Common Are Intragene Windows with <Emphasis Type="Italic">K</Emphasis><Subscript>A</Subscript> > <Emphasis Type="Italic">K</Emphasis><Subscript>S</Subscript> Owing to Purifying Selection on Synonymous Mutations?

One method for diagnosing the mode of sequence evolution considers the ratio of nonsynonymous substitutions per nonsynonymous site (K _A) to the corresponding figure for synonymous substitutions (K _S). A ratio (K _A/K _S) greater than unity is taken as evidence for positive selection. This, however, need not necessarily be the case. Notably, there is one instance of a high intragenic K _A/K _S peak, revealed by sliding window analysis and observed in two pairwise comparisons, better accounted for by localised purifying selection on synonymous mutations that affect splicing. Is this example exceptional? To address this we isolate intragenic domains with K _A/K _S > 1 from more than 1000 long mouse-rat orthologues. Approximately one K _A/K _S > 1 peak is found per 12–15 kb of coding sequence. Surprisingly, low synonymous substitution rates underpin more incidences than do high nonsynonymous rates. Several reasons, however, prevent us from supposing that the low synonymous rates reflect purifying selection on synonymous mutations. First, for many peaks, the null that the peak is no higher than expected given the underlying rates of evolution, cannot be rejected. Second, of 18 statistically significant incidences with unusually low K _S values, only 3 are repeatable across independent comparisons. At least two of these are within alternatively spliced exons. We conclude that repeatable statistically significant intragenic domains of low intragenic K _S are rare. As so few K _A/K _S peaks reflect increased rates of protein evolution and so few hold statistical support, we additionally conclude that sliding window analysis to infer domains of positive selection is highly error-prone.     

A survey of genome‐wide single nucleotide polymorphisms through genome resequencing in the Périgord black truffle (Tuber melanosporum Vittad.)

The Périgord black truffle (Tuber melanosporum Vittad.), considered a gastronomic delicacy worldwide, is an ectomycorrhizal filamentous fungus that is ecologically important in Mediterranean French, Italian and Spanish woodlands. In this study, we developed a novel resource of single nucleotide polymorphisms (SNPs) for T. melanosporum using Illumina high‐throughput resequencing. The genome from six T. melanosporum geographical accessions was sequenced to a depth of approximately 20×. These geographical accessions were selected from different populations within the northern and southern regions of the geographical species distribution. Approximately 80% of the reads for each of the six resequenced geographical accessions mapped against the reference T. melanosporum genome assembly, estimating the core genome size of this organism to be approximately 110 Mbp. A total of 442 326 SNPs corresponding to 3540 SNPs/Mbps were identified as being included in all seven genomes. The SNPs occurred more frequently in repeated sequences (85%), although 4501 SNPs were also identified in the coding regions of 2587 genes. Using the ratio of nonsynonymous mutations per nonsynonymous site (pN) to synonymous mutations per synonymous site (pS) and Tajima's D index scanning the whole genome, we were able to identify genomic regions and genes potentially subjected to positive or purifying selection. The SNPs identified represent a valuable resource for future population genetics and genomics studies.     

RecA mediates MgpB and MgpC phase and antigenic variation in Mycoplasma genitalium, but plays a minor role in DNA repair

Mycoplasma genitalium, a sexually transmitted human pathogen, encodes MgpB and MgpC adhesins that undergo phase and antigenic variation through recombination with archived ‘MgPar’ donor sequences. The mechanism and molecular factors required for this genetic variation are poorly understood. In this study, we estimate that sequence variation at the mgpB/C locus occurs in vitro at a frequency of > 1.25 × 10^?4 events per genome per generation using a quantitative anchored PCR assay. This rate was dramatically reduced in a recA deletion mutant and increased in a complemented strain overexpressing RecA. Similarly, the frequency of haemadsorption‐deficient phase variants was reduced in the recA mutant, but restored by complementation. Unlike Escherichia coli, inactivation of recA in M. genitalium had a minimal effect on survival after exposure to mitomycin C or UV irradiation. In contrast, a deletion mutant for the predicted nucleotide excision repair uvrC gene showed growth defects and was exquisitely sensitive to DNA damage. We conclude that M. genitalium RecA has a primary role in mgpB/C–MgPar recombination leading to antigenic and phase variation, yet plays a minor role in DNA repair. Our results also suggest that M. genitalium possesses an active nucleotide excision repair system, possibly representing the main DNA repair pathway in this minimal bacterium.