Quasispecies Analysis of JC Virus DNA Present in Urine of Healthy Subjects

JC virus is a human polyomavirus that infects the majority of people without apparent symptoms in healthy subjects and it is the causative agent of progressive multifocal leucoencephalopathy (PML), a disorder following lytic infection of oligodendrocytes that mainly manifests itself under immunosuppressive conditions. A hallmark for JC virus isolated from PML-brain is the presence of rearrangements in the non-coding control region (NCCR) interspersed between the early and late genes on the viral genome. Such rearrangements are believed to originate from the archetype JC virus which is shed in urine by healthy subjects and PML patients. We applied next generation sequencing to explore the non-coding control region variability in urine of healthy subjects in search for JC virus quasispecies and rearrangements reminiscent of PML. For 61 viral shedders (out of a total of 254 healthy subjects) non-coding control region DNA and VP1 (major capsid protein) coding sequences were initially obtained by Sanger sequencing. Deletions between 1 and 28 nucleotides long appeared in ∼24.5% of the NCCR sequences while insertions were only detected in ∼3.3% of the samples. 454 pyrosequencing was applied on a subset of 54 urine samples demonstrating the existence of JC virus quasispecies in four subjects (∼7.4%). Hence, our results indicate that JC virus DNA in urine is not always restricted to one unique virus variant, but can be a mixture of naturally occurring variants (quasispecies) reflecting the susceptibility of the non-coding control region for genomic rearrangements in healthy individuals. Our findings pave the way to explore the presence of viral quasispecies and the altered viral tropism that might go along with it as a potential risk factor for opportunistic secondary infections such as PML.     

Sexually-Transmitted/Founder HIV-1 Cannot Be Directly Predicted from Plasma or PBMC-Derived Viral Quasispecies in the Transmitting Partner

Objective

Characterization of HIV-1 sequences in newly infected individuals is important for elucidating the mechanisms of viral sexual transmission. We report the identification of transmitted/founder viruses in eight pairs of HIV-1 sexually-infected patients enrolled at the time of primary infection (“recipients”) and their transmitting partners (“donors”).

Methods

Using a single genome-amplification approach, we compared quasispecies in donors and recipients on the basis of 316 and 376 C2V5 env sequences amplified from plasma viral RNA and PBMC-associated DNA, respectively.

Results

Both DNA and RNA sequences indicated very homogeneous viral populations in all recipients, suggesting transmission of a single variant, even in cases of recent sexually transmitted infections (STIs) in donors (n = 2) or recipients (n = 3). In all pairs, the transmitted/founder virus was derived from an infrequent variant population within the blood of the donor. The donor variant sequences most closely related to the recipient sequences were found in plasma samples in 3/8 cases and/or in PBMC samples in 6/8 cases. Although donors were exclusively (n = 4) or predominantly (n = 4) infected by CCR5-tropic (R5) strains, two recipients were infected with highly homogeneous CXCR4/dual-mixed-tropic (X4/DM) viral populations, identified in both DNA and RNA. The proportion of X4/DM quasispecies in donors was higher in cases of X4/DM than R5 HIV transmission (16.7–22.0% versus 0–2.6%), suggesting that X4/DM transmission may be associated with a threshold population of X4/DM circulating quasispecies in donors.

Conclusions

These suggest that a severe genetic bottleneck occurs during subtype B HIV-1 heterosexual and homosexual transmission. Sexually-transmitted/founder virus cannot be directly predicted by analysis of the donor’s quasispecies in plasma and/or PBMC. Additional studies are required to fully understand the traits that confer the capacity to transmit and establish infection, and determine the role of concomitant STIs in mitigating the genetic bottleneck in mucosal HIV transmission.     

Detection of Salmonella spp. Using a Generic and Differential FRET-PCR

To facilitate the detection of Salmonella and to be able to rapidly and conveniently determine the species/subspecies present, we developed and tested a generic and differential FRET-PCR targeting their tetrathionate reductase response regulator gene. The differential pan-Salmonella FRET-PCR we developed successfully detected seven plasmids that contained partial sequences of S. bongori and the six S. enterica subspecies. The detection limit varied from ∼5 copies of target gene/per PCR reaction for S. enterica enterica to ∼200 for S. bongori. Melting curve analysis demonstrated a T _m of ∼68°C for S. enterica enterica, ∼62.5°C for S. enterica houtenae and S. enterica diarizonae, ∼57°C for S. enterica indica, and ∼54°C for S. bongori, S. enterica salamae and S. enterica arizonae. The differential pan-Salmonella FRET-PCR also detected and determined the subspecies of 4 reference strains and 47 Salmonella isolated from clinically ill birds or pigs. Finally, we found it could directly detect and differentiate Salmonella in feline (5/50 positive; 10%; one S. enterica salamae and 4 S. enterica enterica) and canine feces (15/114 positive; 13.2%; all S. enterica enterica). The differential pan-Salmonella FRET-PCR failed to react with 96 non-Salmonella bacterial strains. Our experiments show the differential pan-Salmonella FRET-PCR we developed is a rapid, sensitive and specific method to detect and differentiate Salmonella.     

The Effects of Mechanical Loading on Tendons - An In Vivo and In Vitro Model Study

Mechanical loading constantly acts on tendons, and a better understanding of its effects on the tendons is essential to gain more insights into tendon patho-physiology. This study aims to investigate tendon mechanobiological responses through the use of mouse treadmill running as an in vivo model and mechanical stretching of tendon cells as an in vitro model. In the in vivo study, mice underwent moderate treadmill running (MTR) and intensive treadmill running (ITR) regimens. Treadmill running elevated the expression of mechanical growth factors (MGF) and enhanced the proliferative potential of tendon stem cells (TSCs) in both patellar and Achilles tendons. In both tendons, MTR upregulated tenocyte-related genes: collagen type I (Coll. I ∼10 fold) and tenomodulin (∼3–4 fold), but did not affect non-tenocyte-related genes: LPL (adipocyte), Sox9 (chondrocyte), Runx2 and Osterix (both osteocyte). However, ITR upregulated both tenocyte (Coll. I ∼7–11 fold; tenomodulin ∼4–5 fold) and non-tenocyte-related genes (∼3–8 fold). In the in vitro study, TSCs and tenocytes were stretched to 4% and 8% using a custom made mechanical loading system. Low mechanical stretching (4%) of TSCs from both patellar and Achilles tendons increased the expression of only the tenocyte-related genes (Coll. I ∼5–6 fold; tenomodulin ∼6–13 fold), but high mechanical stretching (8%) increased the expression of both tenocyte (Coll. I ∼28–50 fold; tenomodulin ∼14–48 fold) and non-tenocyte-related genes (2–5-fold). However, in tenocytes, non-tenocyte related gene expression was not altered by the application of either low or high mechanical stretching. These findings indicate that appropriate mechanical loading could be beneficial to tendons because of their potential to induce anabolic changes in tendon cells. However, while excessive mechanical loading caused anabolic changes in tendons, it also induced differentiation of TSCs into non-tenocytes, which may lead to the development of degenerative tendinopathy frequently seen in clinical settings.     

Archaea Appear to Dominate the Microbiome of Inflatella pellicula Deep Sea Sponges

Microbes associated with marine sponges play significant roles in host physiology. Remarkable levels of microbial diversity have been observed in sponges worldwide through both culture-dependent and culture-independent studies. Most studies have focused on the structure of the bacterial communities in sponges and have involved sponges sampled from shallow waters. Here, we used pyrosequencing of 16S rRNA genes to compare the bacterial and archaeal communities associated with two individuals of the marine sponge Inflatella pellicula from the deep-sea, sampled from a depth of 2,900 m, a depth which far exceeds any previous sequence-based report of sponge-associated microbial communities. Sponge-microbial communities were also compared to the microbial community in the surrounding seawater. Sponge-associated microbial communities were dominated by archaeal sequencing reads with a single archaeal OTU, comprising ∼60% and ∼72% of sequences, being observed from Inflatella pellicula. Archaeal sequencing reads were less abundant in seawater (∼11% of sequences). Sponge-associated microbial communities were less diverse and less even than any other sponge-microbial community investigated to date with just 210 and 273 OTUs (97% sequence identity) identified in sponges, with 4 and 6 dominant OTUs comprising ∼88% and ∼89% of sequences, respectively. Members of the candidate phyla, SAR406, NC10 and ZB3 are reported here from sponges for the first time, increasing the number of bacterial phyla or candidate divisions associated with sponges to 43. A minor cohort from both sponge samples (∼0.2% and ∼0.3% of sequences) were not classified to phylum level. A single OTU, common to both sponge individuals, dominates these unclassified reads and shares sequence homology with a sponge associated clone which itself has no known close relative and may represent a novel taxon.     

Cuban history of CRF19 recombinant subtype of HIV-1

CRF19 is a recombinant form of HIV-1 subtypes D, A1 and G, which was first sampled in Cuba in 1999, but was already present there in 1980s. CRF19 was reported almost uniquely in Cuba, where it accounts for ∼25% of new HIV-positive patients and causes rapid progression to AIDS (∼3 years).We analyzed a large data set comprising ∼350 pol and env sequences sampled in Cuba over the last 15 years and ∼350 from Los Alamos database. This data set contained both CRF19 (∼315), and A1, D and G sequences. We performed and combined analyses for the three A1, G and D regions, using fast maximum likelihood approaches, including: (1) phylogeny reconstruction, (2) spatio-temporal analysis of the virus spread, and ancestral character reconstruction for (3) transmission mode and (4) drug resistance mutations (DRMs). We verified these results with a Bayesian approach. This allowed us to acquire new insights on the CRF19 origin and transmission patterns. We showed that CRF19 recombined between 1966 and 1977, most likely in Cuban community stationed in Congo region. We further investigated CRF19 spread on the Cuban province level, and discovered that the epidemic started in 1970s, most probably in Villa Clara, that it was at first carried by heterosexual transmissions, and then quickly spread in the 1980s within the “men having sex with men” (MSM) community, with multiple transmissions back to heterosexuals. The analysis of the transmission patterns of common DRMs found very few resistance transmission clusters.Our results show a very early introduction of CRF19 in Cuba, which could explain its local epidemiological success. Ignited by a major founder event, the epidemic then followed a similar pattern as other subtypes and CRFs in Cuba. The reason for the short time to AIDS remains to be understood and requires specific surveillance, in Cuba and elsewhere.     

Mutability Dynamics of an Emergent Single Stranded DNA Virus in a Na?ve Host

Quasispecies variants and recombination were studied longitudinally in an emergent outbreak of beak and feather disease virus (BFDV) infection in the orange-bellied parrot (Neophema chrysogaster). Detailed health monitoring and the small population size (<300 individuals) of this critically endangered bird provided an opportunity to longitudinally track viral replication and mutation events occurring in a circular, single-stranded DNA virus over a period of four years within a novel bottleneck population. Optimized PCR was used with different combinations of primers, primer walking, direct amplicon sequencing and sequencing of cloned amplicons to analyze BFDV genome variants. Analysis of complete viral genomes (n = 16) and Rep gene sequences (n = 35) revealed that the outbreak was associated with mutations in functionally important regions of the normally conserved Rep gene and immunogenic capsid (Cap) gene with a high evolutionary rate (3.41×10⁻³ subs/site/year) approaching that for RNA viruses; simultaneously we observed significant evidence of recombination hotspots between two distinct progenitor genotypes within orange-bellied parrots indicating early cross-transmission of BFDV in the population. Multiple quasispecies variants were also demonstrated with at least 13 genotypic variants identified in four different individual birds, with one containing up to seven genetic variants. Preferential PCR amplification of variants was also detected. Our findings suggest that the high degree of genetic variation within the BFDV species as a whole is reflected in evolutionary dynamics within individually infected birds as quasispecies variation, particularly when BFDV jumps from one host species to another.     

DNA assembler,an in vivo genetic method for rapid construction of biochemical pathways

The assembly of large recombinant DNA encoding a whole biochemical pathway or genome represents a significant challenge. Here, we report a new method, DNA assembler, which allows the assembly of an entire biochemical pathway in a single step via in vivo homologous recombination in Saccharomyces cerevisiae. We show that DNA assembler can rapidly assemble a functional d-xylose utilization pathway (∼9 kb DNA consisting of three genes), a functional zeaxanthin biosynthesis pathway (∼11 kb DNA consisting of five genes) and a functional combined d-xylose utilization and zeaxanthin biosynthesis pathway (∼19 kb consisting of eight genes) with high efficiencies (70–100%) either on a plasmid or on a yeast chromosome. As this new method only requires simple DNA preparation and one-step yeast transformation, it represents a powerful tool in the construction of biochemical pathways for synthetic biology, metabolic engineering and functional genomics studies.     

Unexpectedly broad target recognition of the CRISPR-mediated virus defence system in the archaeon Sulfolobus solfataricus

The hyperthermophilic archaeon Sulfolobus solfataricus carries an extensive array of clustered regularly interspaced short palindromic repeats (CRISPR) systems able to mediate DNA degradation of invading genetic elements when complementarity to the small CRISPR-derived (cr)RNAs is given. Studying virus defence in vivo with recombinant viral variants, we demonstrate here that an unexpectedly high number of mutations are tolerated between the CRISPR-derived guide RNAs (crRNAs) and their target sequences (protospacer). Up to 15 mismatches in the crRNA still led to ∼50% of DNA degradation, when these mutations were outside the ‘seed’ region. More than 15 mutations were necessary to fully abolished interference. Different from other CRISPR systems investigated in vivo, mutations outside the protospacer region indicated no need for a protospacer adjacent motif sequence to confer DNA interference. However, complementarity of only 3 nucleotides between the repeat-derived 5′ handle of the crRNA and nucleotides adjacent to the protospacer enabled self-recognition, i.e. protection of the host locus. Our findings show commonalities and differences among the various CRISPR-mediated defence systems and suggest that they should not merely be perceived as a ‘first-barrier-defence system’ but may be considered to have a broader mechanism that allows host cells to cope with viruses keeping them at reduced levels.     

Static and Dynamic Factors Limit Chromosomal Replication Complexity in Escherichia coli,Avoiding Dangers of Runaway Overreplication

We define chromosomal replication complexity (CRC) as the ratio of the copy number of the most replicated regions to that of unreplicated regions on the same chromosome. Although a typical CRC of eukaryotic or bacterial chromosomes is 2, rapidly growing Escherichia coli cells induce an extra round of replication in their chromosomes (CRC = 4). There are also E. coli mutants with stable CRC∼6. We have investigated the limits and consequences of elevated CRC in E. coli and found three limits: the “natural” CRC limit of ∼8 (cells divide more slowly); the “functional” CRC limit of ∼22 (cells divide extremely slowly); and the “tolerance” CRC limit of ∼64 (cells stop dividing). While the natural limit is likely maintained by the eclipse system spacing replication initiations, the functional limit might reflect the capacity of the chromosome segregation system, rather than dedicated mechanisms, and the tolerance limit may result from titration of limiting replication factors. Whereas recombinational repair is beneficial for cells at the natural and functional CRC limits, we show that it becomes detrimental at the tolerance CRC limit, suggesting recombinational misrepair during the runaway overreplication and giving a rationale for avoidance of the latter.     

Urine from Treated Cattle Drives Selection for Cephalosporin Resistant Escherichia coli in Soil

The U.S. Food and Drug Administration recently issued new rules for using ceftiofur in food animals in part because of an increasing prevalence of enteric bacteria that are resistant to 3^rd-generation cephalosporins. Parenteral ceftiofur treatment, however, has limited effects on enteric bacteria so we tested the hypothesis that excreted ceftiofur metabolites exert significant selection pressure for ceftiofur-resistant Escherichia coli in soil. Test matrices were prepared by mixing soil with bovine feces and adding urine containing ceftiofur metabolites (CFM) (0 ppm, ∼50 ppm and ∼100 ppm). Matrices were incubated at 23°C or 4°C for variable periods of time after which residual CFM was quantified using a bioassay. Bla _CMY-2 plasmid-bearing ceftiofur resistant (cef^R) E. coli and one-month old calves were used to study the selection effects of CFM and transmission of cef^R bacteria from the environment back to animals. Our studies showed that urinary CFM (∼13 ppm final concentration) is biologically degraded in soil within 2.7 days at 23°C, but persists up to 23.3 days at 4°C. Even short-term persistence in soil provides a >1 log₁₀ advantage to resistant E. coli populations, resulting in significantly prolonged persistence of these bacteria in the soil (∼two months). We further show that resistant strains readily colonize calves by contact with contaminated bedding and without antibiotic selection pressure. Ceftiofur metabolites in urine amplify resistant E. coli populations and, if applicable to field conditions, this effect is far more compelling than reported selection in vivo after parenteral administration of ceftiofur. Because ceftiofur degradation is temperature dependent, these compounds may accumulate during colder months and this could further enhance selection as seasonal temperatures increase. If cost-effective engineered solutions can be developed to limit ex vivo selection, this may limit proliferation for ceftiofur resistant enteric bacteria while preserving the ability to use this important antibiotic in food animal production.     

Genotypic and Functional Impact of HIV-1 Adaptation to Its Host Population during the North American Epidemic

HLA-restricted immune escape mutations that persist following HIV transmission could gradually spread through the viral population, thereby compromising host antiviral immunity as the epidemic progresses. To assess the extent and phenotypic impact of this phenomenon in an immunogenetically diverse population, we genotypically and functionally compared linked HLA and HIV (Gag/Nef) sequences from 358 historic (1979–1989) and 382 modern (2000–2011) specimens from four key cities in the North American epidemic (New York, Boston, San Francisco, Vancouver). Inferred HIV phylogenies were star-like, with approximately two-fold greater mean pairwise distances in modern versus historic sequences. The reconstructed epidemic ancestral (founder) HIV sequence was essentially identical to the North American subtype B consensus. Consistent with gradual diversification of a “consensus-like” founder virus, the median “background” frequencies of individual HLA-associated polymorphisms in HIV (in individuals lacking the restricting HLA[s]) were ∼2-fold higher in modern versus historic HIV sequences, though these remained notably low overall (e.g. in Gag, medians were 3.7% in the 2000s versus 2.0% in the 1980s). HIV polymorphisms exhibiting the greatest relative spread were those restricted by protective HLAs. Despite these increases, when HIV sequences were analyzed as a whole, their total average burden of polymorphisms that were “pre-adapted” to the average host HLA profile was only ∼2% greater in modern versus historic eras. Furthermore, HLA-associated polymorphisms identified in historic HIV sequences were consistent with those detectable today, with none identified that could explain the few HIV codons where the inferred epidemic ancestor differed from the modern consensus. Results are therefore consistent with slow HIV adaptation to HLA, but at a rate unlikely to yield imminent negative implications for cellular immunity, at least in North America. Intriguingly, temporal changes in protein activity of patient-derived Nef (though not Gag) sequences were observed, suggesting functional implications of population-level HIV evolution on certain viral proteins.     

Detection and Quantification of Flavobacterium psychrophilum-Specific Bacteriophages In Vivo in Rainbow Trout upon Oral Administration: Implications for Disease Control in Aquaculture

The use of bacteriophages in the treatment and prevention of infections by the fish pathogen Flavobacterium psychrophilum has attracted increased attention in recent years. It has been shown recently that phage delivery via the parenteral route resulted in immediate distribution of phages to the circulatory system and the different organs. However, little is known about phage dispersal and survival in vivo in rainbow trout after delivery via the oral route. Here we examined the dispersal and survival of F. psychrophilum phage FpV-9 in vivo in juvenile rainbow trout after administration by three different methods—bath, oral intubation into the stomach, and phage-coated feed—with special emphasis on the oral route of delivery. Phages could be detected in all the organs investigated (intestine, spleen, brain, and kidney) 0.5 h postadministration, reaching concentrations as high as ∼10⁵ PFU mg intestine⁻¹ and ∼10³ PFU mg spleen⁻¹ within the first 24 h following the bath and ∼10⁷ PFU mg intestine⁻¹ and ∼10⁴ PFU mg spleen⁻¹ within the first 24 h following oral intubation. The phages were most persistent in the organs for the first 24 h and then decreased exponentially; no phages were detected after 83 h in the organs investigated. Phage administration via feed resulted in the detection of phages in the intestine, spleen, and kidney 1 h after feeding. Average concentrations of ∼10⁴ PFU mg intestine⁻¹ and ∼10¹ PFU mg spleen⁻¹ were found throughout the experimental period (200 h) following continuous delivery of phages with feed. These experiments clearly demonstrate the ability of the phages to survive passage through the fish stomach and to penetrate the intestinal barrier and enter the circulatory system after oral delivery, although the quantity of phages found in the spleen was 100- to 1,000-fold lower than that in the intestine. It was also shown that phages could tolerate long periods of desiccation on the feed pellets, with 60% survival after storage at −80°C, and 10% survival after storage at 5°C, for ∼8 months. Continuous delivery of phages via coated feed pellets constitutes a promising method of treatment and especially prevention of rainbow trout fry syndrome.     

Angiopoietin-Like Protein 3 Promotes Preservation of Stemness during Ex Vivo Expansion of Murine Hematopoietic Stem Cells

Allogeneic hematopoietic stem cell (HSC) transplantations from umbilical cord blood or autologous HSCs for gene therapy purposes are hampered by limited number of stem cells. To test the ability to expand HSCs in vitro prior to transplantation, two growth factor cocktails containing stem cell factor, thrombopoietin, fms-related tyrosine kinase-3 ligand (STF) or stem cell factor, thrombopoietin, insulin-like growth factor-2, fibroblast growth factor-1 (STIF) either with or without the addition of angiopoietin-like protein-3 (Angptl3) were used. Culturing HSCs in STF and STIF media for 7 days expanded long-term repopulating stem cells content in vivo by ∼6-fold and ∼10-fold compared to freshly isolated stem cells. Addition of Angptl3 resulted in increased expansion of these populations by ∼17-fold and ∼32-fold, respectively, and was further supported by enforced expression of Angptl3 in HSCs through lentiviral transduction that also promoted HSC expansion. As expansion of highly purified lineage-negative, Sca-1⁺, c-Kit⁺ HSCs was less efficient than less pure lineage-negative HSCs, Angptl3 may have a direct effect on HCS but also an indirect effect on accessory cells that support HSC expansion. No evidence for leukemia or toxicity was found during long-term follow up of mice transplanted with ex vivo expanded HSCs or manipulated HSC populations that expressed Angptl3. We conclude that the cytokine combinations used in this study to expand HSCs ex vivo enhances the engraftment in vivo. This has important implications for allogeneic umbilical cord-blood derived HSC transplantations and autologous HSC applications including gene therapy.     

O-Linked β-N-Acetylglucosamine (O-GlcNAc) Regulates Emerin Binding to Barrier to Autointegration Factor (BAF) in a Chromatin- and Lamin B-enriched “Niche”

Emerin, a membrane component of nuclear “lamina” networks with lamins and barrier to autointegration factor (BAF), is highly O-GlcNAc-modified (“O-GlcNAcylated”) in mammalian cells. Mass spectrometry analysis revealed eight sites of O-GlcNAcylation, including Ser-53, Ser-54, Ser-87, Ser-171, and Ser-173. Emerin O-GlcNAcylation was reduced ∼50% by S53A or S54A mutation in vitro and in vivo. O-GlcNAcylation was reduced ∼66% by the triple S52A/S53A/S54A mutant, and S173A reduced O-GlcNAcylation of the S52A/S53A/S54A mutant by ∼30%, in vivo. We separated two populations of emerin, A-type lamins and BAF; one population solubilized easily, and the other required sonication and included histones and B-type lamins. Emerin and BAF associated only in histone- and lamin-B-containing fractions. The S173D mutation specifically and selectively reduced GFP-emerin association with BAF by 58% and also increased GFP-emerin hyper-phosphorylation. We conclude that β-N-acetylglucosaminyltransferase, an essential enzyme, controls two regions in emerin. The first region, defined by residues Ser-53 and Ser-54, flanks the LEM domain. O-GlcNAc modification at Ser-173, in the second region, is proposed to promote emerin association with BAF in the chromatin/lamin B “niche.” These results reveal direct control of a conserved LEM domain nuclear lamina component by β-N-acetylglucosaminyltransferase, a nutrient sensor that regulates cell stress responses, mitosis, and epigenetics.     

Deconvoluting the Composition of Low-Frequency Hepatitis C Viral Quasispecies: Comparison of Genotypes and NS3 Resistance-Associated Variants between HCV/HIV Coinfected Hemophiliacs and HCV Monoinfected Patients in Japan

Pre-existing low-frequency resistance-associated variants (RAVs) may jeopardize successful sustained virological responses (SVR) to HCV treatment with direct-acting antivirals (DAAs). However, the potential impact of low-frequency (∼0.1%) mutations, concatenated mutations (haplotypes), and their association with genotypes (Gts) on the treatment outcome has not yet been elucidated, most probably owing to the difficulty in detecting pre-existing minor haplotypes with sufficient length and accuracy. Herein, we characterize a methodological framework based on Illumina MiSeq next-generation sequencing (NGS) coupled with bioinformatics of quasispecies reconstruction (QSR) to realize highly accurate variant calling and genotype-haplotype detection. The core-to-NS3 protease coding sequences in 10 HCV monoinfected patients, 5 of whom had a history of blood transfusion, and 11 HCV/HIV coinfected patients with hemophilia, were studied. Simulation experiments showed that, for minor variants constituting more than 1%, our framework achieved a positive predictive value (PPV) of 100% and sensitivities of 91.7–100% for genotyping and 80.6% for RAV screening. Genotyping analysis indicated the prevalence of dominant Gt1a infection in coinfected patients (6/11 vs 0/10, p = 0.01). For clinical samples, minor genotype overlapping infection was prevalent in HCV/HIV coinfected hemophiliacs (10/11) and patients who experienced whole-blood transfusion (4/5) but none in patients without exposure to blood (0/5). As for RAV screening, the Q80K/R and S122K/R variants were particularly prevalent among minor RAVs observed, detected in 12/21 and 6/21 cases, respectively. Q80K was detected only in coinfected patients, whereas Q80R was predominantly detected in monoinfected patients (1/11 vs 7/10, p < 0.01). Multivariate interdependence analysis revealed the previously unrecognized prevalence of Gt1b-Q80K, in HCV/HIV coinfected hemophiliacs [Odds ratio = 13.4 (3.48–51.9), p < 0.01]. Our study revealed the distinct characteristics of viral quasispecies between the subgroups specified above and the feasibility of NGS and QSR-based genetic deconvolution of pre-existing minor Gts, RAVs, and their interrelationships.     

Deep-Sequencing of the Peach Latent Mosaic Viroid Reveals New Aspects of Population Heterogeneity

Viroids are small circular single-stranded infectious RNAs characterized by a relatively high mutation level. Knowledge of their sequence heterogeneity remains largely elusive and previous studies, using Sanger sequencing, were based on a limited number of sequences. In an attempt to address sequence heterogeneity from a population dynamics perspective, a GF305-indicator peach tree was infected with a single variant of the Avsunviroidae family member Peach latent mosaic viroid (PLMVd). Six months post-inoculation, full-length circular conformers of PLMVd were isolated and deep-sequenced. We devised an original approach to the bioinformatics refinement of our sequence libraries involving important phenotypic data, based on the systematic analysis of hammerhead self-cleavage activity. Two distinct libraries yielded a total of 3,939 different PLMVd variants. Sequence variants exhibiting up to ∼17% of mutations relative to the inoculated viroid were retrieved, clearly illustrating the high level of divergence dynamics within a unique population. While we initially assumed that most positions of the viroid sequence would mutate, we were surprised to discover that ∼50% of positions remained perfectly conserved, including several small stretches as well as a small motif reminiscent of a GNRA tetraloop which are the result of various selective pressures. Using a hierarchical clustering algorithm, the different variants harvested were subdivided into 7 clusters. We found that most sequences contained an average of 4.6 to 6.4 mutations compared to the variant used to initially inoculate the plant. Interestingly, it was possible to reconstitute and compare the sequence evolution of each of these clusters. In doing so, we identified several key mutations. This study provides a reliable pipeline for the treatment of viroid deep-sequencing. It also sheds new light on the extent of sequence variation that a viroid population can sustain, and which may give rise to a quasispecies.     

Conformational Properties of Polyglutamine Sequences in Guanidine Hydrochloride Solutions

Two sets of iso-1-cytochrome c variants have been prepared with N-terminal insertions of pure polyglutamine, i.e., PolyQ variants, or polyglutamine interrupted with lysine every sixth residue, i.e., Gln-rich variants. The polymer properties of these pure polyGln or Gln-rich sequences have been evaluated using equilibrium and kinetic His-heme loop formation methods for loop sizes ranging from 22 to 46 in 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). In 6.0 M GdnHCl, the scaling exponent, ν₃, for the pure polyGln sequences, is ∼1.7—significantly less than ν₃ ≈ 2.15 for the Gln-rich sequences. The stability of the His-heme loops becomes progressively greater for the pure polyGln sequences relative to the Gln-rich sequences as GdnHCl concentration decreases from 6.0 to 1.5 M. Thus, the context of the sequence effects the polymer properties of Gln repeats even in denaturing concentrations of GdnHCl. Comparison of data for the Gln-rich variants with previous results for Gly-rich and Ala-rich variants shows that ν₃ ∼ 2.2 for the Gln-rich, Gly-rich, and Ala-rich sequences in 6.0 M GdnHCl, whereas ν₃ remains unchanged at 3.0 M GdnHCl concentration for the Gln-rich and Ala-rich sequences but decreases to ∼1.7 for the Gly-rich sequences. Thus, the polymer properties of Gln-rich and Ala-rich sequences are less sensitive to solvent quality in denaturing solutions of GdnHCl than Gly-rich sequences. Evaluation of Flory’s characteristic ratio, C_n, for the Gln-rich and Ala-rich sequences relative to the Gly-rich sequences shows that Gln-rich sequences are stiffer than Ala-rich sequences at both 3.0 and 6.0 M GdnHCl.     

Dual roles for the Mss116 cofactor during splicing of the ai5γ group II intron

The autocatalytic group II intron ai5γ from Saccharomyces cerevisiae self-splices under high-salt conditions in vitro, but requires the assistance of the DEAD-box protein Mss116 in vivo and under near-physiological conditions in vitro. Here, we show that Mss116 influences the folding mechanism in several ways. By comparing intron precursor RNAs with long (∼300 nt) and short (∼20 nt) exons, we observe that long exon sequences are a major obstacle for self-splicing in vitro. Kinetic analysis indicates that Mss116 not only mitigates the inhibitory effects of long exons, but also assists folding of the intron core. Moreover, a mutation in conserved Motif III that impairs unwinding activity (SAT → AAA) only affects the construct with long exons, suggesting helicase unwinding during exon unfolding, but not in intron folding. Strong parallels between Mss116 and the related protein Cyt-19 from Neurospora crassa suggest that these proteins form a subclass of DEAD-box proteins that possess a versatile repertoire of diverse activities for resolving the folding problems of large RNAs.