A New Borrelia Species Defined by Multilocus Sequence Analysis of Housekeeping Genes

Analysis of Lyme borreliosis (LB) spirochetes, using a novel multilocus sequence analysis scheme, revealed that OspA serotype 4 strains (a rodent-associated ecotype) of Borrelia garinii were sufficiently genetically distinct from bird-associated B. garinii strains to deserve species status. We suggest that OspA serotype 4 strains be raised to species status and named Borrelia bavariensis sp. nov. The rooted phylogenetic trees provide novel insights into the evolutionary history of LB spirochetes.Multilocus sequence typing (MLST) and multilocus sequence analysis (MLSA) have been shown to be powerful and pragmatic molecular methods for typing large numbers of microbial strains for population genetics studies, delineation of species, and assignment of strains to defined bacterial species (4, 13, 27, 40, 44). To date, MLST/MLSA schemes have been applied only to a few vector-borne microbial populations (1, 6, 30, 37, 40, 41, 47).Lyme borreliosis (LB) spirochetes comprise a diverse group of zoonotic bacteria which are transmitted among vertebrate hosts by ixodid (hard) ticks. The most common agents of human LB are Borrelia burgdorferi (sensu stricto), Borrelia afzelii, Borrelia garinii, Borrelia lusitaniae, and Borrelia spielmanii (7, 8, 12, 35). To date, 15 species have been named within the group of LB spirochetes (6, 31, 32, 37, 38, 41). While several of these LB species have been delineated using whole DNA-DNA hybridization (3, 20, 33), most ecological or epidemiological studies have been using single loci (5, 9-11, 29, 34, 36, 38, 42, 51, 53). Although some of these loci have been convenient for species assignment of strains or to address particular epidemiological questions, they may be unsuitable to resolve evolutionary relationships among LB species, because it is not possible to define any outgroup. For example, both the 5S-23S intergenic spacer (5S-23S IGS) and the gene encoding the outer surface protein A (ospA) are present only in LB spirochete genomes (36, 43). The advantage of using appropriate housekeeping genes of LB group spirochetes is that phylogenetic trees can be rooted with sequences of relapsing fever spirochetes. This renders the data amenable to detailed evolutionary studies of LB spirochetes.LB group spirochetes differ remarkably in their patterns and levels of host association, which are likely to affect their population structures (22, 24, 46, 48). Of the three main Eurasian Borrelia species, B. afzelii is adapted to rodents, whereas B. valaisiana and most strains of B. garinii are maintained by birds (12, 15, 16, 23, 26, 45). However, B. garinii OspA serotype 4 strains in Europe have been shown to be transmitted by rodents (17, 18) and, therefore, constitute a distinct ecotype within B. garinii. These strains have also been associated with high pathogenicity in humans, and their finer-scale geographical distribution seems highly focal (10, 34, 52, 53).In this study, we analyzed the intra- and interspecific phylogenetic relationships of B. burgdorferi, B. afzelii, B. garinii, B. valaisiana, B. lusitaniae, B. bissettii, and B. spielmanii by means of a novel MLSA scheme based on chromosomal housekeeping genes (30, 48).     

XRCC2 and XRCC3 Regulate the Balance between Short- and Long-Tract Gene Conversions between Sister Chromatids

Sister chromatid recombination (SCR) is a potentially error-free pathway for the repair of DNA lesions associated with replication and is thought to be important for suppressing genomic instability. The mechanisms regulating the initiation and termination of SCR in mammalian cells are poorly understood. Previous work has implicated all the Rad51 paralogs in the initiation of gene conversion and the Rad51C/XRCC3 complex in its termination. Here, we show that hamster cells deficient in the Rad51 paralog XRCC2, a component of the Rad51B/Rad51C/Rad51D/XRCC2 complex, reveal a bias in favor of long-tract gene conversion (LTGC) during SCR. This defect is corrected by expression of wild-type XRCC2 and also by XRCC2 mutants defective in ATP binding and hydrolysis. In contrast, XRCC3-mediated homologous recombination and suppression of LTGC are dependent on ATP binding and hydrolysis. These results reveal an unexpectedly general role for Rad51 paralogs in the control of the termination of gene conversion between sister chromatids.DNA double-strand breaks (DSBs) are potentially dangerous lesions, since their misrepair may cause chromosomal translocations, gene amplifications, loss of heterozygosity (LOH), and other types of genomic instability characteristic of human cancers (7, 9, 21, 40, 76, 79). DSBs are repaired predominantly by nonhomologous end joining or homologous recombination (HR), two evolutionarily conserved DSB repair mechanisms (8, 12, 16, 33, 48, 60, 71). DSBs generated during the S or G₂ phase of the cell cycle may be repaired preferentially by HR, using the intact sister chromatid as a template for repair (12, 26, 29, 32, 71). Sister chromatid recombination (SCR) is a potentially error-free pathway for the repair of DSBs, which has led to the proposal that SCR protects against genomic instability, cancer, and aging. Indeed, a number of human cancer predisposition genes are implicated in SCR control (10, 24, 45, 57, 75).HR entails an initial processing of the DSB to generate a free 3′ single-stranded DNA (ssDNA) overhang (25, 48, 56). This is coupled to the loading of Rad51, the eukaryotic homolog of Escherichia coli RecA, which polymerizes to form an ssDNA-Rad51 “presynaptic” nucleoprotein filament. Formation of the presynaptic filament is tightly regulated and requires the concerted action of a large number of gene products (55, 66, 68). Rad51-coated ssDNA engages in a homology search by invading homologous duplex DNA. If sufficient homology exists between the invading and invaded strands, a triple-stranded synapse (D-loop) forms, and the 3′ end of the invading (nascent) strand is extended, using the donor as a template for gene conversion. This recombination intermediate is thought to be channeled into one of the following two major subpathways: classical gap repair or synthesis-dependent strand annealing (SDSA) (48). Gap repair entails the formation of a double Holliday junction, which may resolve into either crossover or noncrossover products. Although this is a major pathway in meiotic recombination, crossing-over is highly suppressed in somatic eukaryotic cells (26, 44, 48). Indeed, the donor DNA molecule is seldom rearranged during somatic HR, suggesting that SDSA is the major pathway for the repair of somatic DSBs (26, 44, 49, 69). SDSA terminates when the nascent strand is displaced from the D-loop and pairs with the second end of the DSB to form a noncrossover product. The mechanisms underlying displacement of the nascent strand are not well understood. However, failure to displace the nascent strand might be expected to result in the production of longer gene conversion tracts during HR (36, 44, 48, 63).Gene conversion triggered in response to a Saccharomyces cerevisiae or mammalian chromosomal DSB generally results in the copying of a short (50- to 300-bp) stretch of information from the donor (short-tract gene conversion [STGC]) (14, 47, 48, 67, 69). A minority of gene conversions in mammalian cells entail more-extensive copying, generating gene conversion tracts that are up to several kilobases in length (long-tract gene conversion [LTGC]) (26, 44, 51, 54, 64). In yeast, very long gene conversions can result from break-induced replication (BIR), a highly processive form of gene conversion in which a bona fide replication fork is thought to be established at the recombination synapse (11, 36, 37, 39, 61, 63). In contrast, SDSA does not require lagging-strand polymerases and appears to be much less processive than a conventional replication fork (37, 42, 78). BIR in yeast has been proposed to play a role in LOH in aging yeast, telomere maintenance, and palindromic gene amplification (5, 41, 52). It is unclear to what extent a BIR-like mechanism operates in mammalian cells, although BIR has been invoked to explain telomere elongation in tumors lacking telomerase (13). It is currently unknown whether LTGC and STGC in somatic mammalian cells are products of mechanistically distinct pathways or whether they represent alternative outcomes of a common SDSA pathway.Vertebrate cells contain five Rad51 paralogs—polypeptides with limited sequence homology to Rad51—Rad51B, Rad51C, Rad51D, XRCC2, and XRCC3 (74). The Rad51 paralogs form the following two major complexes: Rad51B/Rad51C/Rad51D/XRCC2 (BCDX2) and Rad51C/XRCC3 (CX3) (38, 73). Genetic deletion of any one of the rad51 paralogs in the mouse germ line produces early embryonic lethality, and mouse or chicken cells lacking any of the rad51 paralogs reveal hypersensitivity to DNA-damaging agents, reduced frequencies of HR and of sister chromatid exchanges, increased chromatid-type errors, and defective sister chromatid cohesion (18, 72, 73, 82). Collectively, these data implicate the Rad51 paralogs in SCR regulation. The purified Rad51B/Rad51C complex has been shown to assist Rad51-mediated strand exchange (62). XRCC3 null or Rad51C null hamster cells reveal a bias toward production of longer gene conversion tracts, suggesting a role for the CX3 complex in late stages of SDSA (6, 44). Rad51C copurifies with branch migration and Holliday junction resolution activities in mammalian cell extracts (35), and XRCC3, but not XRCC2, facilitates telomere shortening by reciprocal crossing-over in telomeric T loops (77). These data, taken together with the meiotic defects observed in Rad51C hypomorphic mice, suggest a specialized role for CX3, but not for BCDX2, in resolving Holliday junction structures (31, 58).To further address the roles of Rad51 paralogs in late stages of recombination, we have studied the balance between long-tract (>1-kb) and short-tract (<1-kb) SCR in XRCC2 mutant hamster cells. We found that DSB-induced gene conversion in both XRCC2 and XRCC3 mutant cells is biased in favor of LTGC. These defects were suppressed by expression of wild-type (wt) XRCC2 or XRCC3, respectively, although the dependence upon ATP binding and hydrolysis differed between the two Rad51 paralogs. These results indicate that Rad51 paralogs play a more general role in determining the balance between STGC and LTGC than was previously appreciated and suggest roles for both the BCDX2 and CX3 complexes in influencing the termination of gene conversion in mammals.     

Virus Entry via the Alternative Coreceptors CCR3 and FPRL1 Differs by Human Immunodeficiency Virus Type 1 Subtype

Human immunodeficiency virus type 1 (HIV-1) infects target cells by binding to CD4 and a chemokine receptor, most commonly CCR5. CXCR4 is a frequent alternative coreceptor (CoR) in subtype B and D HIV-1 infection, but the importance of many other alternative CoRs remains elusive. We have analyzed HIV-1 envelope (Env) proteins from 66 individuals infected with the major subtypes of HIV-1 to determine if virus entry into highly permissive NP-2 cell lines expressing most known alternative CoRs differed by HIV-1 subtype. We also performed linear regression analysis to determine if virus entry via the major CoR CCR5 correlated with use of any alternative CoR and if this correlation differed by subtype. Virus pseudotyped with subtype B Env showed robust entry via CCR3 that was highly correlated with CCR5 entry efficiency. By contrast, viruses pseudotyped with subtype A and C Env proteins were able to use the recently described alternative CoR FPRL1 more efficiently than CCR3, and use of FPRL1 was correlated with CCR5 entry. Subtype D Env was unable to use either CCR3 or FPRL1 efficiently, a unique pattern of alternative CoR use. These results suggest that each subtype of circulating HIV-1 may be subject to somewhat different selective pressures for Env-mediated entry into target cells and suggest that CCR3 may be used as a surrogate CoR by subtype B while FPRL1 may be used as a surrogate CoR by subtypes A and C. These data may provide insight into development of resistance to CCR5-targeted entry inhibitors and alternative entry pathways for each HIV-1 subtype.Human immunodeficiency virus type 1 (HIV-1) infects target cells by binding first to CD4 and then to a coreceptor (CoR), of which C-C chemokine receptor 5 (CCR5) is the most common (6, 53). CXCR4 is an additional CoR for up to 50% of subtype B and D HIV-1 isolates at very late stages of disease (4, 7, 28, 35). Many other seven-membrane-spanning G-protein-coupled receptors (GPCRs) have been identified as alternative CoRs when expressed on various target cell lines in vitro, including CCR1 (76, 79), CCR2b (24), CCR3 (3, 5, 17, 32, 60), CCR8 (18, 34, 38), GPR1 (27, 65), GPR15/BOB (22), CXCR5 (39), CXCR6/Bonzo/STRL33/TYMSTR (9, 22, 25, 45, 46), APJ (26), CMKLR1/ChemR23 (49, 62), FPLR1 (67, 68), RDC1 (66), and D6 (55). HIV-2 and simian immunodeficiency virus SIVmac isolates more frequently show expanded use of these alternative CoRs than HIV-1 isolates (12, 30, 51, 74), and evidence that alternative CoRs other than CXCR4 mediate infection of primary target cells by HIV-1 isolates is sparse (18, 30, 53, 81). Genetic deficiency in CCR5 expression is highly protective against HIV-1 transmission (21, 36), establishing CCR5 as the primary CoR. The importance of alternative CoRs other than CXCR4 has remained elusive despite many studies (1, 30, 70, 81). Expansion of CoR use from CCR5 to include CXCR4 is frequently associated with the ability to use additional alternative CoRs for viral entry (8, 16, 20, 63, 79) in most but not all studies (29, 33, 40, 77, 78). This finding suggests that the sequence changes in HIV-1 env required for use of CXCR4 as an additional or alternative CoR (14, 15, 31, 37, 41, 57) are likely to increase the potential to use other alternative CoRs.We have used the highly permissive NP-2/CD4 human glioma cell line developed by Soda et al. (69) to classify virus entry via the alternative CoRs CCR1, CCR3, CCR8, GPR1, CXCR6, APJ, CMKLR1/ChemR23, FPRL1, and CXCR4. Full-length molecular clones of 66 env genes from most prevalent HIV-1 subtypes were used to generate infectious virus pseudotypes expressing a luciferase reporter construct (19, 57). Two types of analysis were performed: the level of virus entry mediated by each alternative CoR and linear regression of entry mediated by CCR5 versus all other alternative CoRs. We thus were able to identify patterns of alternative CoR use that were subtype specific and to determine if use of any alternative CoR was correlated or independent of CCR5-mediated entry. The results obtained have implications for the evolution of env function, and the analyses revealed important differences between subtype B Env function and all other HIV-1 subtypes.     

Microtubule-Nucleus Interactions in Dictyostelium discoideum Mediated by Central Motor Kinesins

Kinesins are a diverse superfamily of motor proteins that drive organelles and other microtubule-based movements in eukaryotic cells. These motors play important roles in multiple events during both interphase and cell division. Dictyostelium discoideum contains 13 kinesin motors, 12 of which are grouped into nine families, plus one orphan. Functions for 11 of the 13 motors have been previously investigated; we address here the activities of the two remaining kinesins, both isoforms with central motor domains. Kif6 (of the kinesin-13 family) appears to be essential for cell viability. The partial knockdown of Kif6 with RNA interference generates mitotic defects (lagging chromosomes and aberrant spindle assemblies) that are consistent with kinesin-13 disruptions in other organisms. However, the orphan motor Kif9 participates in a completely novel kinesin activity, one that maintains a connection between the microtubule-organizing center (MTOC) and nucleus during interphase. kif9 null cell growth is impaired, and the MTOC appears to disconnect from its normally tight nuclear linkage. Mitotic spindles elongate in a normal fashion in kif9⁻ cells, but we hypothesize that this kinesin is important for positioning the MTOC into the nuclear envelope during prophase. This function would be significant for the early steps of cell division and also may play a role in regulating centrosome replication.Directed cell migration, organelle transport, and cell division involve fundamental motilities that are necessary for eukaryotic cell viability and function. Much of the force required for these motilities is generated through the cyclical interactions of motor proteins with the cell cytoskeleton. Microtubules (MTs) and actin filaments provide structural support and directional guides, and all eukaryotic organisms have diverse, often extensive families of motors that carry out different tasks. Functional studies have revealed that many of the motors work in combination with others, and that the individual deletion of a single motor activity often is insufficient to produce a defect that substantially impairs cell growth or function. The latter phenomenon is particularly evident in some organisms with simple motor families (14, 42). By contrasting homologous motor functions between simple and complex systems, we hope to learn the details of how each motor is custom-tuned for specific tasks.Dictyostelium discoideum is a compact amoeba that exhibits robust forms of motility common to nearly all animal cells, with speeds that frequently exceed corresponding rates in vertebrate cell models (25, 33, 54). Since Dictyostelium possesses a relatively small number of motor proteins (13 kinesin, 1 dynein, and 13 myosin isoforms [23, 24, 26]), it combines advantages of terrific cytology with straightforward molecular genetics and thus represents an excellent model to investigate individual and combined motor protein actions. To date, 11 of the 13 kinesin motors have been analyzed functionally (5, 17, 18, 30, 42, 46, 51, 60). Only 1 of these 11 motors, Kif3, a member of the kinesin-1 family of organelle transporters, appears to be essential for organism viability (51). Individual disruptions of three kinesin genes (kif1, kif4, and kif12) produce distinctive defects in cell growth or organelle transport (30, 42, 46). Analyses of six of the seven other kinesins reveal important phenotypes but only when combined with other motor disruptions or cell stresses. We address here the roles of the remaining two Dictyostelium MT-based motors.kif6 and kif9 encode two central motor kinesins in the Dictyostelium genome (24). The best-studied isoforms of this motor type are represented by the kinesin-13 family, and they largely function to regulate MT length during cell division (13, 16, 40, 41). In some organisms, kinesin-13 motors also have been shown to operate during interphase and to mediate MT and flagellar length control (3, 4, 15) and perhaps even organelle transport (32, 43, 56). kif6 encodes the kinesin-13 family member in Dictyostelium. We demonstrate that Kif6 activity is essential for viability, and that it plays a primary, conserved role in chromosome segregation during cell division.The second of the central motor kinesins, Kif9, does not group with an existing family (24, 38). The gene disruption of this motor reveals a completely novel function for a kinesin in maintaining a connection between the MT-organizing center (MTOC) and the nucleus. By electron microscopy (EM), the MTOC of Dictyostelium appears as a cytoplasmic cube-shaped structure surrounded by amorphous dense material (39, 44). EM, biochemical analyses, antibody labeling, and live-cell imaging studies have demonstrated that during interphase, the cytoplasmic MTOC is firmly and closely attached to the nucleus (28, 29, 44, 48, 49, 63). Upon entry into mitosis, the MTOC duplicates during prophase and is brought to or into a fenestration of the nuclear envelope, and then it establishes an intranuclear bipolar spindle for division (39, 53, 64). While MTOCs can be purified from Dictyostelium, the methods rely heavily on reagents that actively disrupt the attached nuclei (10, 59). A recent study has identified at least one component of this connection, the nuclear envelope protein Sun-1 (67). The perturbation of Sun-1 affects nuclear shape and results in centrosome detachment, hyperamplification, and aneuploidy. We demonstrate in the current work that the disruption of the Kif9 kinesin also perturbs the MTOC-nucleus linkage. Our results suggest that an MT-mediated mechanism plays a significant role in maintaining an MTOC-nucleus connection during interphase, and we discuss how this connection could be important to regulate centrosome replication and ensure proper chromosome segregation during cell division.     

Roles of Small,Acid-Soluble Spore Proteins and Core Water Content in Survival of Bacillus subtilis Spores Exposed to Environmental Solar UV Radiation

Spores of Bacillus subtilis contain a number of small, acid-soluble spore proteins (SASP) which comprise up to 20% of total spore core protein. The multiple α/β-type SASP have been shown to confer resistance to UV radiation, heat, peroxides, and other sporicidal treatments. In this study, SASP-defective mutants of B. subtilis and spores deficient in dacB, a mutation leading to an increased core water content, were used to study the relative contributions of SASP and increased core water content to spore resistance to germicidal 254-nm and simulated environmental UV exposure (280 to 400 nm, 290 to 400 nm, and 320 to 400 nm). Spores of strains carrying mutations in sspA, sspB, and both sspA and sspB (lacking the major SASP-α and/or SASP-β) were significantly more sensitive to 254-nm and all polychromatic UV exposures, whereas the UV resistance of spores of the sspE strain (lacking SASP-γ) was essentially identical to that of the wild type. Spores of the dacB-defective strain were as resistant to 254-nm UV-C radiation as wild-type spores. However, spores of the dacB strain were significantly more sensitive than wild-type spores to environmental UV treatments of >280 nm. Air-dried spores of the dacB mutant strain had a significantly higher water content than air-dried wild-type spores. Our results indicate that α/β-type SASP and decreased spore core water content play an essential role in spore resistance to environmentally relevant UV wavelengths whereas SASP-γ does not.Spores of Bacillus spp. are highly resistant to inactivation by different physical stresses, such as toxic chemicals and biocidal agents, desiccation, pressure and temperature extremes, and high fluences of UV or ionizing radiation (reviewed in references 33, 34, and 48). Under stressful environmental conditions, cells of Bacillus spp. produce endospores that can stay dormant for extended periods. The reason for the high resistance of bacterial spores to environmental extremes lies in the structure of the spore. Spores possess thick layers of highly cross-linked coat proteins, a modified peptidoglycan spore cortex, a low core water content, and abundant intracellular constituents, such as the calcium chelate of dipicolinic acid and α/β-type small, acid-soluble spore proteins (α/β-type SASP), the last two of which protect spore DNA (6, 42, 46, 48, 52). DNA damage accumulated during spore dormancy is also efficiently repaired during spore germination (33, 47, 48). UV-induced DNA photoproducts are repaired by spore photoproduct lyase and nucleotide excision repair, DNA double-strand breaks (DSB) by nonhomologous end joining, and oxidative stress-induced apurinic/apyrimidinic (AP) sites by AP endonucleases and base excision repair (15, 26-29, 34, 43, 53, 57).Monochromatic 254-nm UV radiation has been used as an efficient and cost-effective means of disinfecting surfaces, building air, and drinking water supplies (31). Commonly used test organisms for inactivation studies are bacterial spores, usually spores of Bacillus subtilis, due to their high degree of resistance to various sporicidal treatments, reproducible inactivation response, and safety (1, 8, 19, 31, 48). Depending on the Bacillus species analyzed, spores are 10 to 50 times more resistant than growing cells to 254-nm UV radiation. In addition, most of the laboratory studies of spore inactivation and radiation biology have been performed using monochromatic 254-nm UV radiation (33, 34). Although 254-nm UV-C radiation is a convenient germicidal treatment and relevant to disinfection procedures, results obtained by using 254-nm UV-C are not truly representative of results obtained using UV wavelengths that endospores encounter in their natural environments (34, 42, 50, 51, 59). However, sunlight reaching the Earth''s surface is not monochromatic 254-nm radiation but a mixture of UV, visible, and infrared radiation, with the UV portion spanning approximately 290 to 400 nm (33, 34, 36). Thus, our knowledge of spore UV resistance has been constructed largely using a wavelength of UV radiation not normally reaching the Earth''s surface, even though ample evidence exists that both DNA photochemistry and microbial responses to UV are strongly wavelength dependent (2, 30, 33, 36).Of recent interest in our laboratories has been the exploration of factors that confer on B. subtilis spores resistance to environmentally relevant extreme conditions, particularly solar UV radiation and extreme desiccation (23, 28, 30, 34 36, 48, 52). It has been reported that α/β-type SASP but not SASP-γ play a major role in spore resistance to 254-nm UV-C radiation (20, 21) and to wet heat, dry heat, and oxidizing agents (48). In contrast, increased spore water content was reported to affect B. subtilis spore resistance to moist heat and hydrogen peroxide but not to 254-nm UV-C (12, 40, 48). However, the possible roles of SASP-α, -β, and -γ and core water content in spore resistance to environmentally relevant solar UV wavelengths have not been explored. Therefore, in this study, we have used B. subtilis strains carrying mutations in the sspA, sspB, sspE, sspA and sspB, or dacB gene to investigate the contributions of SASP and increased core water content to the resistance of B. subtilis spores to 254-nm UV-C and environmentally relevant polychromatic UV radiation encountered on Earth''s surface.     

Widely Conserved Recombination Patterns among Single-Stranded DNA Viruses

The combinatorial nature of genetic recombination can potentially provide organisms with immediate access to many more positions in sequence space than can be reached by mutation alone. Recombination features particularly prominently in the evolution of a diverse range of viruses. Despite rapid progress having been made in the characterization of discrete recombination events for many species, little is currently known about either gross patterns of recombination across related virus families or the underlying processes that determine genome-wide recombination breakpoint distributions observable in nature. It has been hypothesized that the networks of coevolved molecular interactions that define the epistatic architectures of virus genomes might be damaged by recombination and therefore that selection strongly influences observable recombination patterns. For recombinants to thrive in nature, it is probably important that the portions of their genomes that they have inherited from different parents work well together. Here we describe a comparative analysis of recombination breakpoint distributions within the genomes of diverse single-stranded DNA (ssDNA) virus families. We show that whereas nonrandom breakpoint distributions in ssDNA virus genomes are partially attributable to mechanistic aspects of the recombination process, there is also a significant tendency for recombination breakpoints to fall either outside or on the peripheries of genes. In particular, we found significantly fewer recombination breakpoints within structural protein genes than within other gene types. Collectively, these results imply that natural selection acting against viruses expressing recombinant proteins is a major determinant of nonrandom recombination breakpoint distributions observable in most ssDNA virus families.Genetic recombination is a ubiquitous biological process that is both central to DNA repair pathways (10, 57) and an important evolutionary mechanism. By generating novel combinations of preexisting nucleotide polymorphisms, recombination can potentially accelerate evolution by increasing the population-wide genetic diversity upon which adaptive selection relies. Recombination can paradoxically also prevent the progressive accumulation of harmful mutations within individual genomes (18, 35, 53). Whereas its ability to defend high-fitness genomes from mutational decay possibly underlies the evolutionary value of sexuality in higher organisms, in many microbial species where pseudosexual genetic exchange is permissible among even highly divergent genomes, recombination can enable access to evolutionary innovations that would otherwise be inaccessible by mutation alone.Such interspecies recombination is fairly common in many virus families (8, 17, 27, 44, 82). It is becoming clear, however, that as with mutation events, most recombination events between distantly related genomes are maladaptive (5, 13, 38, 50, 63, 80). As genetic distances between parental genomes increase, so too does the probability of fitness defects in their recombinant offspring (16, 51). The viability of recombinants is apparently largely dependent on how severely recombination disrupts coevolved intragenome interaction networks (16, 32, 51). These networks include interacting nucleotide sequences that form secondary structures, sequence-specific protein-DNA interactions, interprotein interactions, and amino acid-amino acid interactions within protein three-dimensional folds.One virus family where such interaction networks appear to have a large impact on patterns of natural interspecies recombination are the single-stranded DNA (ssDNA) geminiviruses. As with other ssDNA viruses, recombination is very common among the species of this family (62, 84). Partially conserved recombination hot and cold spots have been detected in different genera (39, 81) and are apparently caused by both differential mechanistic predispositions of genome regions to recombination and natural selection disfavoring the survival of recombinants with disrupted intragenome interaction networks (38, 51).Genome organization and rolling circle replication (RCR)—the mechanism by which geminiviruses and many other ssDNA viruses replicate (9, 67, 79; see reference 24 for a review)—seem to have a large influence on basal recombination rates in different parts of geminivirus genomes (20, 33, 39, 61, 81). To initiate RCR, virion-strand ssDNA molecules are converted by host-mediated pathways into double-stranded “replicative-form” (RF) DNAs (34, 67). Initiated by a virus-encoded replication-associated protein (Rep) at a well-defined virion-strand replication origin (v-ori), new virion strands are synthesized on the complementary strand of RF DNAs (28, 73, 74) by host DNA polymerases. Virion-strand replication is concomitant with the displacement of old virion strands, which, once complete, yields covalently closed ssDNA molecules which are either encapsidated or converted into additional RF DNAs. Genome-wide basal recombination rates in ssDNA viruses are probably strongly influenced by the specific characteristics of host DNA polymerases that enable RCR. Interruption of RCR has been implicated directly in geminivirus recombination (40) and is most likely responsible for increased basal recombination rates both within genes transcribed in the opposite direction from that of virion-strand replication (40, 71) and at the v-ori (1, 9, 20, 69, 74).Whereas most ssDNA virus families replicate via either a rolling circle mechanism (the Nanoviridae, Microviridae, and Geminiviridae) (3, 23, 24, 31, 59, 67, 74) or a related rolling hairpin mechanism (the Parvoviridae) (25, 76), among the Circoviridae only the Circovirus genus is known to use RCR (45). Although the Gyrovirus genus (the other member of the Circoviridae) and the anelloviruses (a currently unclassified ssDNA virus group) might also use RCR, it is currently unknown whether they do or not (78). Additionally, some members of the Begomovirus genus of the Geminiviridae either have a second genome component, called DNA-B, or are associated with satellite ssDNA molecules called DNA-1 and DNA-Beta, all of which also replicate by RCR (1, 47, 68).Recombination is known to occur in the parvoviruses (19, 43, 70), microviruses (66), anelloviruses (40, 46), circoviruses (11, 26, 60), nanoviruses (30), geminivirus DNA-B components, and geminivirus satellite molecules (2, 62). Given that most, if not all, of these ssDNA replicons are evolutionarily related to and share many biological features with the geminiviruses (22, 31, 36), it is of interest to determine whether conserved recombination patterns observed in the geminiviruses (61, 81) are evident in these other groups. To date, no comparative analyses have ever been performed with different ssDNA virus families to identify, for example, possible influences of genome organization on recombination breakpoint distributions found in these viruses.Here we compare recombination frequencies and recombination breakpoint distributions in most currently described ssDNA viruses and satellite molecules and identify a number of sequence exchange patterns that are broadly conserved across this entire group.     

Population Structure of the Lyme Borreliosis Spirochete Borrelia burgdorferi in the Western Black-Legged Tick (Ixodes pacificus) in Northern California

Factors potentially contributing to the lower incidence of Lyme borreliosis (LB) in the far-western than in the northeastern United States include tick host-seeking behavior resulting in fewer human tick encounters, lower densities of Borrelia burgdorferi-infected vector ticks in peridomestic environments, and genetic variation among B. burgdorferi spirochetes to which humans are exposed. We determined the population structure of B. burgdorferi in over 200 infected nymphs of the primary bridging vector to humans, Ixodes pacificus, collected in Mendocino County, CA. This was accomplished by sequence typing the spirochete lipoprotein ospC and the 16S-23S rRNA intergenic spacer (IGS). Thirteen ospC alleles belonging to 12 genotypes were found in California, and the two most abundant, ospC genotypes H3 and E3, have not been detected in ticks in the Northeast. The most prevalent ospC and IGS biallelic profile in the population, found in about 22% of ticks, was a new B. burgdorferi strain defined by ospC genotype H3. Eight of the most common ospC genotypes in the northeastern United States, including genotypes I and K that are associated with disseminated human infections, were absent in Mendocino County nymphs. ospC H3 was associated with hardwood-dominated habitats where western gray squirrels, the reservoir host, are commonly infected with LB spirochetes. The differences in B. burgdorferi population structure in California ticks compared to the Northeast emphasize the need for a greater understanding of the genetic diversity of spirochetes infecting California LB patients.In the United States, Lyme borreliosis (LB) is the most commonly reported vector-borne illness and is caused by infection with the spirochete Borrelia burgdorferi (3, 9, 52). The signs and symptoms of LB can include a rash, erythema migrans, fever, fatigue, arthritis, carditis, and neurological manifestations (50, 51). The black-legged tick, Ixodes scapularis, and the western black-legged tick, Ixodes pacificus, are the primary vectors of B. burgdorferi to humans in the United States, with the former in the northeastern and north-central parts of the country and the latter in the Far West (9, 10). These ticks perpetuate enzootic transmission cycles together with a vertebrate reservoir host such as the white-footed mouse, Peromyscus leucopus, in the Northeast and Midwest (24, 35), or the western gray squirrel, Sciurus griseus, in California (31, 46).B. burgdorferi is a spirochete species with a largely clonal population structure (14, 16) comprising several different strains or lineages (8). The polymorphic ospC gene of B. burgdorferi encodes a surface lipoprotein that increases expression within the tick during blood feeding (47) and is required for initial infection of mammalian hosts (25, 55). To date, approximately 20 North American ospC genotypes have been described (40, 45, 49, 56). At least four, and possibly up to nine, of these genotypes are associated with B. burgdorferi invasiveness in humans (1, 15, 17, 49, 57). Restriction fragment length polymorphism (RFLP) and, subsequently, sequence analysis of the 16S-23S rRNA intergenic spacer (IGS) are used as molecular typing tools to investigate genotypic variation in B. burgdorferi (2, 36, 38, 44, 44, 57). The locus maintains a high level of variation between related species, and this variation reflects the heterogeneity found at the genomic level of the organism (37). The IGS and ospC loci appear to be linked (2, 8, 26, 45, 57), but the studies to date have not been representative of the full range of diversity of B. burgdorferi in North America.Previous studies in the northeastern and midwestern United States have utilized IGS and ospC genotyping to elucidate B. burgdorferi evolution, host strain specificity, vector-reservoir associations, and disease risk to humans. In California, only six ospC and five IGS genotypes have been described heretofore in samples from LB patients or I. pacificus ticks (40, 49, 56) compared to approximately 20 ospC and IGS genotypes identified in ticks, vertebrate hosts, or humans from the Northeast and Midwest (8, 40, 45, 49, 56). Here, we employ sequence analysis of both the ospC gene and IGS region to describe the population structure of B. burgdorferi in more than 200 infected I. pacificus nymphs from Mendocino County, CA, where the incidence of LB is among the highest in the state (11). Further, we compare the Mendocino County spirochete population to populations found in the Northeast.     

Mixed Infections of Pepino Mosaic Virus Strains Modulate the Evolutionary Dynamics of this Emergent Virus

Pepino mosaic virus (PepMV) is an emerging pathogen that causes severe economic losses in tomato crops (Solanum lycopersicum L.) in the Northern hemisphere, despite persistent attempts of control. In fact, it is considered one of the most significant viral diseases for tomato production worldwide, and it may constitute a good model for the analysis of virus emergence in crops. We have combined a population genetics approach with an analysis of in planta properties of virus strains to explain an observed epidemiological pattern. Hybridization analysis showed that PepMV populations are composed of isolates of two types (PepMV-CH2 and PepMV-EU) that cocirculate. The CH2 type isolates are predominant; however, EU isolates have not been displaced but persist mainly in mixed infections. Two molecularly cloned isolates belonging to each type have been used to examine the dynamics of in planta single infections and coinfection, revealing that the CH2 type has a higher fitness than the EU type. Coinfections expand the range of susceptible hosts, and coinfected plants remain symptomless several weeks after infection, so a potentially important problem for disease prevention and management. These results provide an explanation of the observed epidemiological pattern in terms of genetic and ecological interactions among the different viral strains. Thus, mixed infections appear to be contributing to shaping the genetic structure and dynamics of PepMV populations.Pepino mosaic virus (PepMV; genus Potexvirus, family Flexiviridae) was identified in 1974 as the agent responsible for a viral disease of pepino crops (Solanum muricatum) in Peru (30). PepMV in tomato (Solanum lycopersicum) was first reported in The Netherlands in 1999 (74) but has since spread rapidly in Europe (3, 11, 38, 48, 51, 57) and beyond (20, 35, 36, 42, 68), causing epidemics and severe economic losses (27, 29, 36, 51, 67, 69). The PepMV host range is limited mainly to the Solanaceae (59), and the virus is easily transmitted from plant to plant by contact (30), vectored by bumblebees (65), or seedborne-transmitted (37). PepMV infections in tomato are associated with a wide range of leaf symptoms: mild and severe mosaics, bubbling, laminal distortions, and stunting (26, 27, 51). Fruit symptoms occur with or without leaf symptoms, and the main impact of PepMV is on fruit quality (irregular lycopene distribution [26]) but not on yield (69). Therefore, PepMV is currently considered a dangerous pathogen and is included in the European Plant Protection Organization alert list (15) as one of the most important tomato viruses worldwide (27, 51, 57, 68, 69).The PepMV genome consists of a single, positive-sense, ∼6,400-nucleotide (nt) RNA strand containing five open reading frames (ORFs). ORF1 encodes the putative viral polymerase (RdRp) (3). ORFs 2, 3, and 4 encode the triple gene block (TGB) proteins TGBp1, TGBp2, and TGBp3, which are essential for virus movement (46, 75, 78). Potato virus X TGBp1 is a multifunctional protein that induces plasmodesmal gating, moves from cell to cell, has ATPase and RNA helicase activities, binds viral RNAs, and acts as suppressor of RNA silencing (39, 76-78). ORF5 encodes the coat protein (CP) which, in addition to its structural role, is required for cell-to-cell and long-distance movement (12). Finally, two short untranslated sequences flank the coding regions, and there is a poly(A) tail at the 3′ end of the genomic RNA (3, 11, 48).Previous studies have shown that Spanish PepMV populations sampled between 2000 and 2004 were genetically very homogeneous (∼99% nucleotide identity), most comprising isolates highly similar to the so-called European tomato strain (PepMV-EU). However, a few isolates sampled in 2004 in the Murcia region (Southeastern Spain) were distinct and highly similar to the US2 strain reported in the United States (51). U.S. isolates (US1 and US2) and a Chilean isolate from infected tomato seeds (CH2) share only 79 to 86% nucleotide identity with European (EU) isolates (36, 42). The CH2 type has been reported recently in greenhouses for tomato production in Poland (29) and Belgium (27). In this last study, CH2 was predominant in single infections and also frequent in mixed infections with isolates of the EU type (27). However, all PepMV types (EU, US1, US2, and CH2) have been found in United States, where the PepMV-EU type has been the most prevalent, and mixed infections were found in samples collected from Arizona, Colorado, and Texas (35).Several studies of plant virus populations have reported a reduced genetic diversity of populations separated in time or space (19, 40, 56) with high virus genetic stability (23). Nevertheless, how genetic and ecological factors modulate the evolutionary dynamics of viruses and determine epidemiological patterns is still poorly understood (25, 47).We have characterized the population genetic structure of PepMV in infected samples of commercial tomato crops in the Murcia region (southeastern Spain) between 2005 and 2008. Phylogenetic analysis was performed, and genetic diversity values among PepMV isolates were estimated to determine the structure of the population and the strength and direction of selection. In addition, the biological properties (host range, fitness, and virulence) of two cloned isolates of the CH2 and EU types were studied to understand the evolutionary dynamics of natural PepMV populations.