Genomic Confirmation of Borrelia garinii, United States.

Lyme disease is a multisystem disorder primarily caused by Borrelia burgdorferi sensu lato. However, B. garinii, which has been identified on islands off the coast of Newfoundland and Labrador, Canada, is a cause of Lyme disease in Eurasia. We report isolation and whole-genome nucleotide sequencing of a B. garinii isolate from a cotton mouse (Peromyscus gossypinus) in South Carolina, USA. We identified a second B. garinii isolate from the same repository. Phylogenetic analysis does not associate these isolates with the previously described isolates of B. garinii from Canada.

Multipartite Genome of Lyme Disease Borrelia: Structure, Variation and Prophages.

All members of the Borrelia genus that have been examined harbour a linear chromosome that is about 900 kbp in length, as well as a plethora of both linear and circular plasmids in the 5-220 kbp size range. Genome sequences for 27 Lyme disease Borrelia isolates have been determined since the elucidation of the B. burgdorferi B31 genome sequence in 1997. The chromosomes, which carry the vast majority of the housekeeping genes, appear to be very constant in gene content and organization across all Lyme disease Borrelia species. The content of the plasmids, which carry most of the genes that encode the differentially expressed surface proteins that interact with the spirochete's arthropod and vertebrate hosts, is much more variable. Lyme disease Borrelia isolates carry between 7-21 different plasmids, ranging in size from 5-84 kbp. All strains analyzed to date harbor three plasmids, cp26, lp54 and lp17. The plasmids are unusual, as compared to most bacterial plasmids, in that they contain many paralogous sequences, a large number of pseudogenes, and, in some cases, essential genes. In addition, a number of the plasmids have features indicating that they are prophages. Numerous methods have been developed for Lyme disease Borrelia strain typing. These have proven valuable for clinical and epidemiological studies, as well as phylogenomic and population genetic analyses. Increasingly, these approaches have been displaced by whole genome sequencing techniques. Some correlations between genome content and pathogenicity have been deduced, and comparative whole genome analyses promise future progress in this arena.

Analysis of a flagellar filament cap mutant reveals that HtrA serine protease degrades unfolded flagellin protein in the periplasm of Borrelia burgdorferi.

Unlike external flagellated bacteria, spirochetes have periplasmic flagella (PF). Very little is known about how PF are assembled within the periplasm of spirochaetal cells. Herein, we report that FliD (BB0149), a flagellar cap protein (also named hook-associated protein 2), controls flagellin stability and flagellar filament assembly in the Lyme disease spirochete Borrelia burgdorferi. Deletion of fliD leads to non-motile mutant cells that are unable to assemble flagellar filaments and pentagon-shaped caps (10 nm in diameter, 12 nm in length). Interestingly, FlaB, a major flagellin protein of B. burgdorferi, is degraded in the fliD mutant but not in other flagella-deficient mutants (i.e., in the hook, rod, or MS-ring). Biochemical and genetic studies reveal that HtrA, a serine protease of B. burgdorferi, controls FlaB turnover. Specifically, HtrA degrades unfolded but not polymerized FlaB, and deletion of htrA increases the level of FlaB in the fliD mutant. Collectively, we propose that the flagellar cap protein FliD promotes flagellin polymerization and filament growth in the periplasm. Deletion of fliD abolishes this process, which leads to leakage of unfolded FlaB proteins into the periplasm where they are degraded by HtrA, a protease that prevents accumulation of toxic products in the periplasm.

Systems-level analysis of NalD mutation, a recurrent driver of rapid drug resistance in acute Pseudomonas aeruginosa infection.

Pseudomonas aeruginosa, a main cause of human infection, can gain resistance to the antibiotic aztreonam through a mutation in NalD, a transcriptional repressor of cellular efflux. Here we combine computational analysis of clinical isolates, transcriptomics, metabolic modeling and experimental validation to find a strong association between NalD mutations and resistance to aztreonam-as well as resistance to other antibiotics-across P. aeruginosa isolated from different patients. A detailed analysis of one patient's timeline shows how this mutation can emerge in vivo and drive rapid evolution of resistance while the patient received cancer treatment, a bone marrow transplantation, and antibiotics up to the point of causing the patient's death. Transcriptomics analysis confirmed the primary mechanism of NalD action-a loss-of-function mutation that caused constitutive overexpression of the MexAB-OprM efflux system-which lead to aztreonam resistance but, surprisingly, had no fitness cost in the absence of the antibiotic. We constrained a genome-scale metabolic model using the transcriptomics data to investigate changes beyond the primary mechanism of resistance, including adaptations in major metabolic pathways and membrane transport concurrent with aztreonam resistance, which may explain the lack of a fitness cost. We propose that metabolic adaptations may allow resistance mutations to endure in the absence of antibiotics and could be targeted by future therapies against antibiotic resistant pathogens.

BpWrapper: BioPerl-based sequence and tree utilities for rapid prototyping of bioinformatics pipelines.

BACKGROUND: Automated bioinformatics workflows are more robust, easier to maintain, and results more reproducible when built with command-line utilities than with custom-coded scripts. Command-line utilities further benefit by relieving bioinformatics developers to learn the use of, or to interact directly with, biological software libraries. There is however a lack of command-line utilities that leverage popular Open Source biological software toolkits such as BioPerl ( http://bioperl.org ) to make many of the well-designed, robust, and routinely used biological classes available for a wider base of end users. RESULTS: Designed as standard utilities for UNIX-family operating systems, BpWrapper makes functionality of some of the most popular BioPerl modules readily accessible on the command line to novice as well as to experienced bioinformatics practitioners. The initial release of BpWrapper includes four utilities with concise command-line user interfaces, bioseq, bioaln, biotree, and biopop, specialized for manipulation of molecular sequences, sequence alignments, phylogenetic trees, and DNA polymorphisms, respectively. Over a hundred methods are currently available as command-line options and new methods are easily incorporated. Performance of BpWrapper utilities lags that of precompiled utilities while equivalent to that of other utilities based on BioPerl. BpWrapper has been tested on BioPerl Release 1.6, Perl versions 5.10.1 to 5.25.10, and operating systems including Apple macOS, Microsoft Windows, and GNU/Linux. Release code is available from the Comprehensive Perl Archive Network (CPAN) at https://metacpan.org/pod/Bio::BPWrapper . Source code is available on GitHub at https://github.com/bioperl/p5-bpwrapper . CONCLUSIONS: BpWrapper improves on existing sequence utilities by following the design principles of Unix text utilities such including a concise user interface, extensive command-line options, and standard input/output for serialized operations. Further, dozens of novel methods for manipulation of sequences, alignments, and phylogenetic trees, unavailable in existing utilities (e.g., EMBOSS, Newick Utilities, and FAST), are provided. Bioinformaticians should find BpWrapper useful for rapid prototyping of workflows on the command-line without creating custom scripts for comparative genomics and other bioinformatics applications.

Primordial origin and diversification of plasmids in Lyme disease agent bacteria.

BACKGROUND: With approximately one-third of their genomes consisting of linear and circular plasmids, the Lyme disease agent cluster of species has the most complex genomes among known bacteria. We report here a comparative analysis of plasmids in eleven Borreliella (also known as Borrelia burgdorferi sensu lato) species. RESULTS: We sequenced the complete genomes of two B. afzelii, two B. garinii, and individual B. spielmanii, B. bissettiae, B. valaisiana and B. finlandensis isolates. These individual isolates carry between seven and sixteen plasmids, and together harbor 99 plasmids. We report here a comparative analysis of these plasmids, along with 70 additional Borreliella plasmids available in the public sequence databases. We identify only one new putative plasmid compatibility type (the 30th) among these 169 plasmid sequences, suggesting that all or nearly all such types have now been discovered. We find that the linear plasmids in the non-B. burgdorferi species have undergone the same kinds of apparently random, chaotic rearrangements mediated by non-homologous recombination that we previously discovered in B. burgdorferi. These rearrangements occurred independently in the different species lineages, and they, along with an expanded chromosomal phylogeny reported here, allow the identification of several whole plasmid transfer events among these species. Phylogenetic analyses of the plasmid partition genes show that a majority of the plasmid compatibility types arose early, most likely before separation of the Lyme agent Borreliella and relapsing fever Borrelia clades, and this, with occasional cross species plasmid transfers, has resulted in few if any species-specific or geographic region-specific Borreliella plasmid types. CONCLUSIONS: The primordial origin and persistent maintenance of the Borreliella plasmid types support their functional indispensability as well as evolutionary roles in facilitating genome diversity. The improved resolution of Borreliella plasmid phylogeny based on conserved partition-gene clusters will lead to better determination of gene orthology which is essential for prediction of biological function, and it will provide a basis for inferring detailed evolutionary mechanisms of Borreliella genomic variability including homologous gene and plasmid exchanges as well as non-homologous rearrangements.

Genotyping and Quantifying Lyme Pathogen Strains by Deep Sequencing of the Outer Surface Protein C (ospC) Locus.

A mixed infection of a single tick or host by Lyme disease spirochetes is common and a unique challenge for the diagnosis, treatment, and surveillance of Lyme disease. Here, we describe a novel protocol for differentiating Lyme strains on the basis of deep sequencing of the hypervariable outer surface protein C locus (ospC). Improving upon the traditional DNA-DNA hybridization method, the next-generation sequencing-based protocol is high throughput, quantitative, and able to detect new pathogen strains. We applied the method to more than one hundred infected Ixodes scapularis ticks collected from New York State, USA, in 2015 and 2016. An analysis of strain distributions within individual ticks suggests an overabundance of multiple infections by five or more strains, inhibitory interactions among coinfecting strains, and the presence of a new strain closely related to Borreliella bissettiae A supporting bioinformatics pipeline has been developed. The newly designed pair of universal ospC primers target intergenic sequences conserved among all known Lyme pathogens. The protocol could be used for culture-free identification and quantification of Lyme pathogens in wildlife and potentially in clinical specimens.

Bow-tie signaling in c-di-GMP: Machine learning in a simple biochemical network.

Bacteria of many species rely on a simple molecule, the intracellular secondary messenger c-di-GMP (Bis-(3'-5')-cyclic dimeric guanosine monophosphate), to make a vital choice: whether to stay in one place and form a biofilm, or to leave it in search of better conditions. The c-di-GMP network has a bow-tie shaped architecture that integrates many signals from the outside world-the input stimuli-into intracellular c-di-GMP levels that then regulate genes for biofilm formation or for swarming motility-the output phenotypes. How does the 'uninformed' process of evolution produce a network with the right input/output association and enable bacteria to make the right choice? Inspired by new data from 28 clinical isolates of Pseudomonas aeruginosa and strains evolved in laboratory experiments we propose a mathematical model where the c-di-GMP network is analogous to a machine learning classifier. The analogy immediately suggests a mechanism for learning through evolution: adaptation though incremental changes in c-di-GMP network proteins acquires knowledge from past experiences and enables bacteria to use it to direct future behaviors. Our model clarifies the elusive function of the ubiquitous c-di-GMP network, a key regulator of bacterial social traits associated with virulence. More broadly, the link between evolution and machine learning can help explain how natural selection across fluctuating environments produces networks that enable living organisms to make sophisticated decisions.

Proteome-wide analysis of mutant p53 targets in breast cancer identifies new levels of gain-of-function that influence PARP, PCNA, and MCM4.

The gain-of-function mutant p53 (mtp53) transcriptome has been studied, but, to date, no detailed analysis of the mtp53-associated proteome has been described. We coupled cell fractionation with stable isotope labeling with amino acids in cell culture (SILAC) and inducible knockdown of endogenous mtp53 to determine the mtp53-driven proteome. Our fractionation data highlight the underappreciated biology that missense mtp53 proteins R273H, R280K, and L194F are tightly associated with chromatin. Using SILAC coupled to tandem MS, we identified that R273H mtp53 expression in MDA-MB-468 breast cancer cells up- and down-regulated multiple proteins and metabolic pathways. Here we provide the data set obtained from sequencing 73,154 peptide pairs that then corresponded to 3,010 proteins detected under reciprocal labeling conditions. Importantly, the high impact regulated targets included the previously identified transcriptionally regulated mevalonate pathway proteins but also identified two new levels of mtp53 protein regulation for nontranscriptional targets. Interestingly, mtp53 depletion profoundly influenced poly(ADP ribose) polymerase 1 (PARP1) localization, with increased cytoplasmic and decreased chromatin-associated protein. An enzymatic PARP shift occurred with high mtp53 expression, resulting in increased poly-ADP-ribosylated proteins in the nucleus. Mtp53 increased the level of proliferating cell nuclear antigen (PCNA) and minichromosome maintenance 4 (MCM4) proteins without changing the amount of pcna and mcm4 transcripts. Pathway enrichment analysis ranked the DNA replication pathway above the cholesterol biosynthesis pathway as a R273H mtp53 activated proteomic target. Knowledge of the proteome diversity driven by mtp53 suggests that DNA replication and repair pathways are major targets of mtp53 and highlights consideration of combination chemotherapeutic strategies targeting cholesterol biosynthesis and PARP inhibition.

Plasmid diversity and phylogenetic consistency in the Lyme disease agent Borrelia burgdorferi.

BACKGROUND: Bacteria from the genus Borrelia are known to harbor numerous linear and circular plasmids. We report here a comparative analysis of the nucleotide sequences of 236 plasmids present in fourteen independent isolates of the Lyme disease agent B. burgdorferi. RESULTS: We have sequenced the genomes of 14 B. burgdorferi sensu stricto isolates that carry a total of 236 plasmids. These individual isolates carry between seven and 23 plasmids. Their chromosomes, the cp26 and cp32 circular plasmids, as well as the lp54 linear plasmid, are quite evolutionarily stable; however, the remaining plasmids have undergone numerous non-homologous and often duplicative recombination events. We identify 32 different putative plasmid compatibility types among the 236 plasmids, of which 15 are (usually) circular and 17 are linear. Because of past rearrangements, any given gene, even though it might be universally present in these isolates, is often found on different linear plasmid compatibility types in different isolates. For example, the arp gene and the vls cassette region are present on plasmids of four and five different compatibility types, respectively, in different isolates. A majority of the plasmid types have more than one organizationally different subtype, and the number of such variants ranges from one to eight among the 18 linear plasmid types. In spite of this substantial organizational diversity, the plasmids are not so variable that every isolate has a novel version of every plasmid (i.e., there appears to be a limited number of extant plasmid subtypes). CONCLUSIONS: Although there have been many past recombination events, both homologous and nonhomologous, among the plasmids, particular organizational variants of these plasmids correlate with particular chromosomal genotypes, suggesting that there has not been rapid horizontal transfer of whole linear plasmids among B. burgdorferi lineages. We argue that plasmid rearrangements are essentially non-revertable and are present at a frequency of only about 0.65% that of single nucleotide changes, making rearrangement-derived novel junctions (mosaic boundaries) ideal phylogenetic markers in the study of B. burgdorferi population structure and plasmid evolution and exchange.

BorreliaBase: a phylogeny-centered browser of Borrelia genomes.

BACKGROUND: The bacterial genus Borrelia (phylum Spirochaetes) consists of two groups of pathogens represented respectively by B. burgdorferi, the agent of Lyme borreliosis, and B. hermsii, the agent of tick-borne relapsing fever. The number of publicly available Borrelia genomic sequences is growing rapidly with the discovery and sequencing of Borrelia strains worldwide. There is however a lack of dedicated online databases to facilitate comparative analyses of Borrelia genomes. DESCRIPTION: We have developed BorreliaBase, an online database for comparative browsing of Borrelia genomes. The database is currently populated with sequences from 35 genomes of eight Lyme-borreliosis (LB) group Borrelia species and 7 Relapsing-fever (RF) group Borrelia species. Distinct from genome repositories and aggregator databases, BorreliaBase serves manually curated comparative-genomic data including genome-based phylogeny, genome synteny, and sequence alignments of orthologous genes and intergenic spacers. CONCLUSIONS: With a genome phylogeny at its center, BorreliaBase allows online identification of hypervariable lipoprotein genes, potential regulatory elements, and recombination footprints by providing evolution-based expectations of sequence variability at each genomic locus. The phylo-centric design of BorreliaBase (http://borreliabase.org) is a novel model for interactive browsing and comparative analysis of bacterial genomes online.

Gene Set Enrichment Analysis (GSEA) of Toxoplasma gondii expression datasets links cell cycle progression and the bradyzoite developmental program.

BACKGROUND: Large amounts of microarray expression data have been generated for the Apicomplexan parasite Toxoplasma gondii in an effort to identify genes critical for virulence or developmental transitions. However, researchers' ability to analyze this data is limited by the large number of unannotated genes, including many that appear to be conserved hypothetical proteins restricted to Apicomplexa. Further, differential expression of individual genes is not always informative and often relies on investigators to draw big-picture inferences without the benefit of context. We hypothesized that customization of gene set enrichment analysis (GSEA) to T. gondii would enable us to rigorously test whether groups of genes serving a common biological function are co-regulated during the developmental transition to the latent bradyzoite form. RESULTS: Using publicly available T. gondii expression microarray data, we created Toxoplasma gene sets related to bradyzoite differentiation, oocyst sporulation, and the cell cycle. We supplemented these with lists of genes derived from community annotation efforts that identified contents of the parasite-specific organelles, rhoptries, micronemes, dense granules, and the apicoplast. Finally, we created gene sets based on metabolic pathways annotated in the KEGG database and Gene Ontology terms associated with gene annotations available at http://www.toxodb.org. These gene sets were used to perform GSEA analysis using two sets of published T. gondii expression data that characterized T. gondii stress response and differentiation to the latent bradyzoite form. CONCLUSIONS: GSEA provides evidence that cell cycle regulation and bradyzoite differentiation are coupled. Δgcn5A mutants unable to induce bradyzoite-associated genes in response to alkaline stress have different patterns of cell cycle and bradyzoite gene expression from stressed wild-type parasites. Extracellular tachyzoites resemble a transitional state that differs in gene expression from both replicating intracellular tachyzoites and in vitro bradyzoites by expressing genes that are enriched in bradyzoites as well as genes that are associated with the G1 phase of the cell cycle. The gene sets we have created are readily modified to reflect ongoing research and will aid researchers' ability to use a knowledge-based approach to data analysis facilitating the development of new insights into the intricate biology of Toxoplasma gondii.

Quantity of questing black-legged ticks and associated micro-scale environmental data collected from four Suburban Parks near New York City.

During 2017 and 2018, we collected the quantity of questing black-legged ticks (Ixodes scapularis), also known as deer ticks, in 124 sampling sites of 5m by 5m in four state parks-Caumsett State Historic Park, Connetquot River State Park, Rockefeller State Park, and Fire Island National Seashore-around New York City. The black-legged tick is the primary vector for the spirochete Borrelia burgdorferi, the pathogen of Lyme disease, in Northeastern United States. Using the flagging method, we collected and counted the numbers of adult and nymphal black-legged ticks at each stie. Along with these quantities, we also recorded the geographic coordinates, ambient temperature, and relative humidity at the sampling sites. Using high-resolution aerial imagery and LiDAR data, we further derived land cover composition, ecotone boundary length, normalized difference vegetation index (NDVI), elevation, solar radiation, and other environmental factors. The data could be used to conduct longitudinal analysis at the same sampling sites as well as comparison with other sites. Ecologists and environmental scientists can use the data for spatiotemporal and statistical analyses of tick ecology at the local scale.

Tumor­suppressive effects of Smad­ubiquitination regulator 2 in papillary thyroid carcinoma.

Smad-ubiquitination regulator 2 (SMURF2) functions as a homolog of E6AP carboxyl terminus-type E3 ubiquitin ligase to regulate cell cycle progression and tumor growth factor expression. SMURF2 has been revealed to function as a tumor suppressor in a number of cancers; however, its function in papillary thyroid carcinoma (PTC) remains largely unknown. Therefore, the aim of the present study was to investigate the function of SMURF2 in PTC. Reverse transcription-quantitative PCR and western blotting were used to detect cellular expression of SMURF2 in vitro. After increasing or inhibiting the expression of SMURF2, MTT was used to detect the effect on tumor cell proliferation and Transwell assays were used to detect the effect on tumor cell migration and invasion. Finally, ELISA was used to detect the effects on glucose and glutamine metabolism in tumor cells and the findings revealed that SMURF2 was downregulated in PTC tissues. Moreover, SMURF2 inhibited the proliferation, invasion and migration of PTC cells, and promoted their apoptosis. Finally, SMURF2 inhibited cell glycolysis and glutaminolysis and affected metabolism in the PTC cell line, TPC-1. Thus, the findings of the present study suggest that SMURF2 may be a potential target in the treatment of PTC.

A CANCER PERSISTENT DNA REPAIR CIRCUIT DRIVEN BY MDM2, MDM4 (MDMX), AND MUTANT P53 FOR RECRUITMENT OF MDC1 AND 53BP1 TO CHROMATIN.

The influence of the metastasis promoting proteins mutant p53 (mtp53) and MDM2 on Cancer Persistent Repair (CPR) to promote cancer cell survival is understudied. Interactions between the DNA repair choice protein 53BP1 and wild type tumor suppressor protein p53 (wtp53) regulates cell cycle control. Cancer cells often express elevated levels of transcriptionally inactive missense mutant p53 (mtp53) that interacts with MDM2 and MDM4/MDMX (herein called MDMX). The ability of mtp53 to maintain a 53BP1 interaction while in the context of interactions with MDM2 and MDMX has not been described. We asked if MDM2 regulates chromatin-based phosphorylation events in the context of mtp53 by comparing the chromatin of T47D breast cancer cells with and without MDM2 in a phospho-peptide stable isotope labeling in cell culture (SILAC) screen. We found reduced phospho-53BP1 chromatin association, which we confirmed by chromatin fractionation and immunofluorescence in multiple breast cancer cell lines. We used the Proximity Ligation Assay (PLA) in breast cancer cell lines and detected 53BP1 in close proximity to mtp53, MDM2, and the DNA repair protein MDC1. Through disruption of the mtp53-MDM2 interaction, by either Nutlin 3a or a mtp53 R273H C-terminal deletion, we uncovered that mtp53 was required for MDM2-53BP1 interaction foci. Our data suggests that mtp53 works with MDM2 and 53BP1 to promote CPR and cell survival.

Inter- and intra-specific pan-genomes of Borrelia burgdorferi sensu lato: genome stability and adaptive radiation.

BACKGROUND: Lyme disease is caused by spirochete bacteria from the Borrelia burgdorferi sensu lato (B. burgdorferi s.l.) species complex. To reconstruct the evolution of B. burgdorferi s.l. and identify the genomic basis of its human virulence, we compared the genomes of 23 B. burgdorferi s.l. isolates from Europe and the United States, including B. burgdorferi sensu stricto (B. burgdorferi s.s., 14 isolates), B. afzelii (2), B. garinii (2), B. "bavariensis" (1), B. spielmanii (1), B. valaisiana (1), B. bissettii (1), and B. "finlandensis" (1). RESULTS: Robust B. burgdorferi s.s. and B. burgdorferi s.l. phylogenies were obtained using genome-wide single-nucleotide polymorphisms, despite recombination. Phylogeny-based pan-genome analysis showed that the rate of gene acquisition was higher between species than within species, suggesting adaptive speciation. Strong positive natural selection drives the sequence evolution of lipoproteins, including chromosomally-encoded genes 0102 and 0404, cp26-encoded ospC and b08, and lp54-encoded dbpA, a07, a22, a33, a53, a65. Computer simulations predicted rapid adaptive radiation of genomic groups as population size increases. CONCLUSIONS: Intra- and inter-specific pan-genome sizes of B. burgdorferi s.l. expand linearly with phylogenetic diversity. Yet gene-acquisition rates in B. burgdorferi s.l. are among the lowest in bacterial pathogens, resulting in high genome stability and few lineage-specific genes. Genome adaptation of B. burgdorferi s.l. is driven predominantly by copy-number and sequence variations of lipoprotein genes. New genomic groups are likely to emerge if the current trend of B. burgdorferi s.l. population expansion continues.

Novelty Search Promotes Antigenic Diversity in Microbial Pathogens.

Driven by host-pathogen coevolution, cell surface antigens are often the fastest evolving parts of a microbial pathogen. The persistent evolutionary impetus for novel antigen variants suggests the utility of novelty-seeking algorithms in predicting antigen diversification in microbial pathogens. In contrast to traditional genetic algorithms maximizing variant fitness, novelty-seeking algorithms optimize variant novelty. Here, we designed and implemented three evolutionary algorithms (fitness-seeking, novelty-seeking, and hybrid) and evaluated their performances in 10 simulated and 2 empirically derived antigen fitness landscapes. The hybrid walks combining fitness- and novelty-seeking strategies overcame the limitations of each algorithm alone, and consistently reached global fitness peaks. Thus, hybrid walks provide a model for microbial pathogens escaping host immunity without compromising variant fitness. Biological processes facilitating novelty-seeking evolution in natural pathogen populations include hypermutability, recombination, wide dispersal, and immune-compromised hosts. The high efficiency of the hybrid algorithm improves the evolutionary predictability of novel antigen variants. We propose the design of escape-proof vaccines based on high-fitness variants covering a majority of the basins of attraction on the fitness landscape representing all potential variants of a microbial antigen.

Upregulation of ubiquitin carboxy­terminal hydrolase 47 (USP47) in papillary thyroid carcinoma ex vivo and reduction of tumor cell malignant behaviors after USP47 knockdown by stabilizing SATB1 expression in vitro.

Aberrant ubiquitination contributes to cancer development, including thyroid carcinoma. The present study assessed the expression of ubiquitin carboxy-terminal hydrolase 47 (USP47) and underlying molecular events in the development of papillary thyroid carcinoma (PTC). The effects of USP47 on PTC cell invasion and migration were analyzed by Transwell assays, while. the effects of USP47 and SATB1on PTC cell gene expression and changes in tumor cell metabolism were assayed by reverse transcription-quantitative PCR, western bolt, or ELISA, respectively. The expression of USP47 mRNA and protein was upregulated in PTC tissue and associated with the PTC tumor size. Knockdown of USP47 expression in PTC cell lines (TPC-1 and K1), decreased the cell proliferation mobility and invasion capacities, whereas USP47 overexpression in these cell lines showed an inverse effect and promoted cell glycolysis and glutamine metabolism. Moreover, expression of special AT-rich sequence-binding protein-1 (SATB1) was high in PTC tissue and was associated with USP47 expression. SATB1 expression promoted tumor cell glycolysis and glutamine metabolism, while USP47 protein bound to and deubiquitinated SATB1 to increase its intracellular levels, thus promoting glycolysis and glutamine metabolism. USP47 promotion of PTC development may be due to its stabilization of SATB1 protein, suggesting that targeting the USP47/SATB1 signaling axis may serve as a therapeutic intervention for PTC.

Whole-genome sequences of Borrelia bissettii, Borrelia valaisiana, and Borrelia spielmanii.

It has been known for decades that human Lyme disease is caused by the three spirochete species Borrelia burgdorferi, Borrelia afzelii, and Borrelia garinii. Recently, Borrelia valaisiana, Borrelia spielmanii, and Borrelia bissettii have been associated with Lyme disease. We report the complete genome sequences of B. valaisiana VS116, B. spielmanii A14S, and B. bissettii DN127.

Accelerated and adaptive evolution of yeast sexual adhesins.

There is a recent emergence of interest in the genes involved in gametic recognition as drivers of reproductive isolation. The recent population genomic sequencing of two species of sexually primitive yeasts (Liti G, Carter DM, Moses AM, Warringer J, Parts L, James SA, Davey RP, Roberts IN, Burt A, Koufopanou V et al. [23 co-authors]. 2009. Population genomics of domestic and wild yeasts. Nature 458:337-341.) has provided data for systematic study of the roles these genes play in the early evolution of sex and speciation. Here, we discovered that among genes encoding cell surface proteins, the sexual adhesin genes have evolved significantly more rapidly than others, both within and between Saccharomyces cerevisiae and its closest relative S. paradoxus. This result was supported by analyses using the PAML pairwise model, a modified McDonald-Kreitman test, and the PAML branch model. Moreover, using a combination of a new statistic of neutrality, an information theory-based measure of evolutionary variability, and functional characterization of amino acid changes, we found that a higher proportion of amino acid changes are fixed in the sexual adhesins than in other proteins and a greater proportion of the fixed amino acid changes either between the two species or the two subgroups of S. paradoxus are functionally dissimilar or radically different. These results suggest that the accelerated evolution of sexual adhesin genes may facilitate speciation, or incipient speciation, and promote sexual selection in general.