Multi-Sample Pooling and Illumina Genome Analyzer Sequencing Methods to Determine Gene Sequence Variation for Database Development

Determination of sequence variation within a genetic locus to develop clinically relevant databases is critical for molecular assay design and clinical test interpretation, so multisample pooling for Illumina genome analyzer (GA) sequencing was investigated using the RET proto-oncogene as a model. Samples were Sanger-sequenced for RET exons 10, 11, and 13–16. Ten samples with 13 known unique variants (“singleton variants” within the pool) and seven common changes were amplified and then equimolar-pooled before sequencing on a single flow cell lane, generating 36 base reads. For comparison, a single “control” sample was run in a different lane. After alignment, a 24-base quality score-screening threshold and 3` read end trimming of three bases yielded low background error rates with a 27% decrease in aligned read coverage. Sequencing data were evaluated using an established variant detection method (percent variant reads), by the presented subtractive correction method, and with SNPSeeker software. In total, 41 variants (of which 23 were singleton variants) were detected in the 10 pool data, which included all Sanger-identified variants. The 23 singleton variants were detected near the expected 5% allele frequency (average 5.17%±0.90% variant reads), well above the highest background error (1.25%). Based on background error rates, read coverage, simulated 30, 40, and 50 sample pool data, expected singleton allele frequencies within pools, and variant detection methods; ≥30 samples (which demonstrated a minimum 1% variant reads for singletons) could be pooled to reliably detect singleton variants by GA sequencing.     

Diversity of Basidiomycetes in Michigan Agricultural Soils

We analyzed the communities of soil basidiomycetes in agroecosystems that differ in tillage history at the Kellogg Biological Station Long-Term Ecological Research site near Battle Creek, Michigan. The approach combined soil DNA extraction through a bead-beating method modified to increase recovery of fungal DNA, PCR amplification with basidiomycete-specific primers, cloning and restriction fragment length polymorphism screening of mixed PCR products, and sequencing of unique clones. Much greater diversity was detected than was anticipated in this habitat on the basis of culture-based methods or surveys of fruiting bodies. With “species” defined as organisms yielding PCR products with ≥99% identity in the 5′ 650 bases of the nuclear large-subunit ribosomal DNA, 241 “species” were detected among 409 unique basidiomycete sequences recovered. Almost all major clades of basidiomycetes from basidiomycetous yeasts and other heterobasidiomycetes through polypores and euagarics (gilled mushrooms and relatives) were represented, with a majority from the latter clade. Only 24 of 241 “species” had 99% or greater sequence similarity to named reference sequences in GenBank, and several clades with multiple “species” could not be identified at the genus level by phylogenetic comparisons with named sequences. The total estimated “species” richness for this 11.2-ha site was 367 “species” of basidiomycetes. Since >99% of the study area has not been sampled, the accuracy of our diversity estimate is uncertain. Replication in time and space is required to detect additional diversity and the underlying community structure.     

Rheostats and Toggle Switches for Modulating Protein Function

The millions of protein sequences generated by genomics are expected to transform protein engineering and personalized medicine. To achieve these goals, tools for predicting outcomes of amino acid changes must be improved. Currently, advances are hampered by insufficient experimental data about nonconserved amino acid positions. Since the property “nonconserved” is identified using a sequence alignment, we designed experiments to recapitulate that context: Mutagenesis and functional characterization was carried out in 15 LacI/GalR homologs (rows) at 12 nonconserved positions (columns). Multiple substitutions were made at each position, to reveal how various amino acids of a nonconserved column were tolerated in each protein row. Results showed that amino acid preferences of nonconserved positions were highly context-dependent, had few correlations with physico-chemical similarities, and were not predictable from their occurrence in natural LacI/GalR sequences. Further, unlike the “toggle switch” behaviors of conserved positions, substitutions at nonconserved positions could be rank-ordered to show a “rheostatic”, progressive effect on function that spanned several orders of magnitude. Comparisons to various sequence analyses suggested that conserved and strongly co-evolving positions act as functional toggles, whereas other important, nonconserved positions serve as rheostats for modifying protein function. Both the presence of rheostat positions and the sequence analysis strategy appear to be generalizable to other protein families and should be considered when engineering protein modifications or predicting the impact of protein polymorphisms.     

Comprehensive Characterization of Human Genome Variation by High Coverage Whole-Genome Sequencing of Forty Four Caucasians

Whole genome sequencing studies are essential to obtain a comprehensive understanding of the vast pattern of human genomic variations. Here we report the results of a high-coverage whole genome sequencing study for 44 unrelated healthy Caucasian adults, each sequenced to over 50-fold coverage (averaging 65.8×). We identified approximately 11 million single nucleotide polymorphisms (SNPs), 2.8 million short insertions and deletions, and over 500,000 block substitutions. We showed that, although previous studies, including the 1000 Genomes Project Phase 1 study, have catalogued the vast majority of common SNPs, many of the low-frequency and rare variants remain undiscovered. For instance, approximately 1.4 million SNPs and 1.3 million short indels that we found were novel to both the dbSNP and the 1000 Genomes Project Phase 1 data sets, and the majority of which (∼96%) have a minor allele frequency less than 5%. On average, each individual genome carried ∼3.3 million SNPs and ∼492,000 indels/block substitutions, including approximately 179 variants that were predicted to cause loss of function of the gene products. Moreover, each individual genome carried an average of 44 such loss-of-function variants in a homozygous state, which would completely “knock out” the corresponding genes. Across all the 44 genomes, a total of 182 genes were “knocked-out” in at least one individual genome, among which 46 genes were “knocked out” in over 30% of our samples, suggesting that a number of genes are commonly “knocked-out” in general populations. Gene ontology analysis suggested that these commonly “knocked-out” genes are enriched in biological process related to antigen processing and immune response. Our results contribute towards a comprehensive characterization of human genomic variation, especially for less-common and rare variants, and provide an invaluable resource for future genetic studies of human variation and diseases.     

Rheostat functional outcomes occur when substitutions are introduced at nonconserved positions that diverge with speciation

When amino acids vary during evolution, the outcome can be functionally neutral or biologically‐important. We previously found that substituting a subset of nonconserved positions, “rheostat” positions, can have surprising effects on protein function. Since changes at rheostat positions can facilitate functional evolution or cause disease, more examples are needed to understand their unique biophysical characteristics. Here, we explored whether “phylogenetic” patterns of change in multiple sequence alignments (such as positions with subfamily specific conservation) predict the locations of functional rheostat positions. To that end, we experimentally tested eight phylogenetic positions in human liver pyruvate kinase (hLPYK), using 10–15 substitutions per position and biochemical assays that yielded five functional parameters. Five positions were strongly rheostatic and three were non‐neutral. To test the corollary that positions with low phylogenetic scores were not rheostat positions, we combined these phylogenetic positions with previously‐identified hLPYK rheostat, “toggle” (most substitution abolished function), and “neutral” (all substitutions were like wild‐type) positions. Despite representing 428 variants, this set of 33 positions was poorly statistically powered. Thus, we turned to the in vivo phenotypic dataset for E. coli lactose repressor protein (LacI), which comprised 12–13 substitutions at 329 positions and could be used to identify rheostat, toggle, and neutral positions. Combined hLPYK and LacI results show that positions with strong phylogenetic patterns of change are more likely to exhibit rheostat substitution outcomes than neutral or toggle outcomes. Furthermore, phylogenetic patterns were more successful at identifying rheostat positions than were co‐evolutionary or eigenvector centrality measures of evolutionary change.     

Longitudinal Analysis of CCR5 and CXCR4 Usage in a Cohort of Antiretroviral Therapy-Na?ve Subjects with Progressive HIV-1 Subtype C Infection

HIV-1 subtype C (C-HIV) is responsible for most HIV-1 cases worldwide. Although the pathogenesis of C-HIV is thought to predominantly involve CCR5-restricted (R5) strains, we do not have a firm understanding of how frequently CXCR4-using (X4 and R5X4) variants emerge in subjects with progressive C-HIV infection. Nor do we completely understand the molecular determinants of coreceptor switching by C-HIV variants. Here, we characterized a panel of HIV-1 envelope glycoproteins (Envs) (n = 300) cloned sequentially from plasma of 21 antiretroviral therapy (ART)-naïve subjects who experienced progression from chronic to advanced stages of C-HIV infection, and show that CXCR4-using C-HIV variants emerged in only one individual. Mutagenesis studies and structural models suggest that the evolution of R5 to X4 variants in this subject principally involved acquisition of an “Ile-Gly” insertion in the gp120 V3 loop and replacement of the V3 “Gly-Pro-Gly” crown with a “Gly-Arg-Gly” motif, but that the accumulation of additional gp120 “scaffold” mutations was required for these V3 loop changes to confer functional effects. In this context, either of the V3 loop changes could confer possible transitional R5X4 phenotypes, but when present together they completely abolished CCR5 usage and conferred the X4 phenotype. Our results show that the emergence of CXCR4-using strains is rare in this cohort of untreated individuals with advanced C-HIV infection. In the subject where X4 variants did emerge, alterations in the gp120 V3 loop were necessary but not sufficient to confer CXCR4 usage.     

Evolutionary interactions between N-linked glycosylation sites in the HIV-1 envelope

The addition of asparagine (N)-linked polysaccharide chains (i.e., glycans) to the gp120 and gp41 glycoproteins of human immunodeficiency virus type 1 (HIV-1) envelope is not only required for correct protein folding, but also may provide protection against neutralizing antibodies as a “glycan shield.” As a result, strong host-specific selection is frequently associated with codon positions where nonsynonymous substitutions can create or disrupt potential N-linked glycosylation sites (PNGSs). Moreover, empirical data suggest that the individual contribution of PNGSs to the neutralization sensitivity or infectivity of HIV-1 may be critically dependent on the presence or absence of other PNGSs in the envelope sequence. Here we evaluate how glycan–glycan interactions have shaped the evolution of HIV-1 envelope sequences by analyzing the distribution of PNGSs in a large-sequence alignment. Using a “covarion”-type phylogenetic model, we find that the rates at which individual PNGSs are gained or lost vary significantly over time, suggesting that the selective advantage of having a PNGS may depend on the presence or absence of other PNGSs in the sequence. Consequently, we identify specific interactions between PNGSs in the alignment using a new paired-character phylogenetic model of evolution, and a Bayesian graphical model. Despite the fundamental differences between these two methods, several interactions are jointly identified by both. Mapping these interactions onto a structural model of HIV-1 gp120 reveals that negative (exclusive) interactions occur significantly more often between colocalized glycans, while positive (inclusive) interactions are restricted to more distant glycans. Our results imply that the adaptive repertoire of alternative configurations in the HIV-1 glycan shield is limited by functional interactions between the N-linked glycans. This represents a potential vulnerability of rapidly evolving HIV-1 populations that may provide useful glycan-based targets for neutralizing antibodies.     

Infectious Molecular Clones with the Nonhomologous Dimer Initiation Sequences Found in Different Subtypes of Human Immunodeficiency Virus Type 1 Can Recombine and Initiate a Spreading Infection In Vitro

Recombinant forms of human immunodeficiency virus type 1 (HIV-1) have been shown to be of major importance in the global AIDS pandemic. Viral RNA dimer formation mediated by the dimerization initiation sequence (DIS) is believed to be essential for viral genomic RNA packaging and therefore for RNA recombination. Here, we demonstrate that HIV-1 recombination and replication are not restricted by variant DIS loop sequences. Three DIS loop forms found among HIV-1 isolates, DIS (CG), DIS (TA), and DIS (TG), when introduced into deletion mutants of HIV-1 recombined efficiently, and the progeny virions replicated with comparable kinetics. A fourth DIS loop form, containing an artificial AAAAAA sequence disrupting the putative DIS loop-loop interactions [DIS (A6)], supported efficient recombination with DIS loop variants; however, DIS (A6) progeny virions exhibited a modest replication disadvantage in mixed cultures. Our studies indicate that the nonhomologous DIS sequences found in different HIV-1 subtypes are not a primary obstacle to intersubtype recombination.     

Prion Variants and Species Barriers Among Saccharomyces Ure2 Proteins

As hamster scrapie cannot infect mice, due to sequence differences in their PrP proteins, we find “species barriers” to transmission of the [URE3] prion in Saccharomyces cerevisiae among Ure2 proteins of S. cerevisiae, paradoxus, bayanus, cariocanus, and mikatae on the basis of differences among their Ure2p prion domain sequences. The rapid variation of the N-terminal Ure2p prion domains results in protection against the detrimental effects of infection by a prion, just as the PrP residue 129 Met/Val polymorphism may have arisen to protect humans from the effects of cannibalism. Just as spread of bovine spongiform encephalopathy prion variant is less impaired by species barriers than is sheep scrapie, we find that some [URE3] prion variants are infectious to another yeast species while other variants (with the identical amino acid sequence) are not. The species barrier is thus prion variant dependent as in mammals. [URE3] prion variant characteristics are maintained even on passage through the Ure2p of another species. Ure2p of Saccharomyces castelli has an N-terminal Q/N-rich “prion domain” but does not form prions (in S. cerevisiae) and is not infected with [URE3] from Ure2p of other Saccharomyces. This implies that conservation of its prion domain is not for the purpose of forming prions. Indeed the Ure2p prion domain has been shown to be important, though not essential, for the nitrogen catabolism regulatory role of the protein.     

Genomic Evolution of Porcine Reproductive and Respiratory Syndrome Virus (PRRSV) Isolates Revealed by Deep Sequencing

Most studies on PRRSV evolution have been limited to a particular region of the viral genome. A thorough genome-wide understanding of the impact of different mechanisms on shaping PRRSV genetic diversity is still lacking. To this end, deep sequencing was used to obtain genomic sequences of a diverse set of 16 isolates from a region of Hong Kong with a complex PRRSV epidemiological record. Genome assemblies and phylogenetic typing indicated the co-circulation of strains of both genotypes (type 1and type 2) with varying Nsp2 deletion patterns and distinct evolutionary lineages (“High Fever”-like and local endemic type). Recombination analyses revealed genomic breakpoints in structural and non-structural regions of genomes of both genotypes with evidence of many recombination events originating from common ancestors. Additionally, the high fold of coverage per nucleotide allowed the characterization of minor variants arising from the quasispecies of each strain. Overall, 0.56–2.83% of sites were found to be polymorphic with respect to cognate consensus genomes. The distribution of minor variants across each genome was not uniform indicating the influence of selective forces. Proportion of variants capable of causing an amino acid change in their respective codons ranged between 25–67% with many predicted to be non-deleterious. Low frequency deletion variants were also detected providing one possible mechanism for their sudden emergence as cited in previous reports.     

Analysis of sequence conservation at nucleotide resolution

One of the major goals of comparative genomics is to understand the evolutionary history of each nucleotide in the human genome sequence, and the degree to which it is under selective pressure. Ascertainment of selective constraint at nucleotide resolution is particularly important for predicting the functional significance of human genetic variation and for analyzing the sequence substructure of cis-regulatory sequences and other functional elements. Current methods for analysis of sequence conservation are focused on delineation of conserved regions comprising tens or even hundreds of consecutive nucleotides. We therefore developed a novel computational approach designed specifically for scoring evolutionary conservation at individual base-pair resolution. Our approach estimates the rate at which each nucleotide position is evolving, computes the probability of neutrality given this rate estimate, and summarizes the result in a Sequence CONservation Evaluation (SCONE) score. We computed SCONE scores in a continuous fashion across 1% of the human genome for which high-quality sequence information from up to 23 genomes are available. We show that SCONE scores are clearly correlated with the allele frequency of human polymorphisms in both coding and noncoding regions. We find that the majority of noncoding conserved nucleotides lie outside of longer conserved elements predicted by other conservation analyses, and are experiencing ongoing selection in modern humans as evident from the allele frequency spectrum of human polymorphism. We also applied SCONE to analyze the distribution of conserved nucleotides within functional regions. These regions are markedly enriched in individually conserved positions and short (<15 bp) conserved “chunks.” Our results collectively suggest that the majority of functionally important noncoding conserved positions are highly fragmented and reside outside of canonically defined long conserved noncoding sequences. A small subset of these fragmented positions may be identified with high confidence.     

Isolation and identification of lymphocytic and myelogenous leukemia-specific sequences in genomes of gibbon oncornaviruses

Five gibbon ape leukemia virus substrains (two from gibbons with lymphocytic leukemia and three from gibbons with myelogenous leukemia) were examined for unique genomic sequences specific for each form of leukemia. By using sequential adsorption procedures, the genome from each gibbon ape leukemia virus was fractionated into four sets of distinct nucleotide sequences. Based on their hybridization specificities toward DNAs of leukemic tissues, these sequences were designated as follows: (i) “COM,” (ii) “LYM” or “MYE,” (iii) “UNI,” and (iv) “UND.” The COM fraction represented sequences common to all of the viral genomes. The LYM fraction, which was isolated only from gibbon ape leukemia viruses associated with lymphocytic leukemia, represented genomic sequences associated with lymphocytic leukemia since the RNA hybridized at a 4- to 15-fold-higher rate to infected tissue DNA from lymphocytic leukemic gibbons than to infected tissue DNA from myelogenous leukemic gibbons. The MYE fraction, which was isolated only from gibbon ape leukemia viruses associated with myelogenous leukemia, represented genomic sequences associated with myelogenous leukemia since the RNA hybridized at a 5- to 15-fold-higher rate to infected tissue DNA from myelogenous leukemic gibbons than to infected tissue DNA from lymphocytic leukemic gibbons. The UNI fraction contained sequences unique to one virus substrain. The UND fraction contained sequences which varied depending upon the substrains involved in the adsorption procedures. These findings suggest that each gibbon ape leukemia virus examined in this study contains subgenomic sequences that are specifically identifiable only with the form of leukemia from which the virus was isolated.     

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.     

Mitochondrial DNA Variant Discovery and Evaluation in Human Cardiomyopathies through Next-Generation Sequencing

Mutations in mitochondrial DNA (mtDNA) may cause maternally-inherited cardiomyopathy and heart failure. In homoplasmy all mtDNA copies contain the mutation. In heteroplasmy there is a mixture of normal and mutant copies of mtDNA. The clinical phenotype of an affected individual depends on the type of genetic defect and the ratios of mutant and normal mtDNA in affected tissues. We aimed at determining the sensitivity of next-generation sequencing compared to Sanger sequencing for mutation detection in patients with mitochondrial cardiomyopathy. We studied 18 patients with mitochondrial cardiomyopathy and two with suspected mitochondrial disease. We “shotgun” sequenced PCR-amplified mtDNA and multiplexed using a single run on Roche''s 454 Genome Sequencer. By mapping to the reference sequence, we obtained 1,300× average coverage per case and identified high-confidence variants. By comparing these to >400 mtDNA substitution variants detected by Sanger, we found 98% concordance in variant detection. Simulation studies showed that >95% of the homoplasmic variants were detected at a minimum sequence coverage of 20× while heteroplasmic variants required >200× coverage. Several Sanger “misses” were detected by 454 sequencing. These included the novel heteroplasmic 7501T>C in tRNA serine 1 in a patient with sudden cardiac death. These results support a potential role of next-generation sequencing in the discovery of novel mtDNA variants with heteroplasmy below the level reliably detected with Sanger sequencing. We hope that this will assist in the identification of mtDNA mutations and key genetic determinants for cardiomyopathy and mitochondrial disease.     

The distribution of different classes of nuclear localization signals (NLSs) in diverse organisms and the utilization of the minor NLS-binding site inplantnuclear import factor importin-α

The specific recognition between the import receptor importin-α and the nuclear localization signals (NLSs) is crucial to ensure the selective transport of cargoes into the nucleus. NLSs contain 1 or 2 clusters of positively charged amino acids, which usually bind to the major (monopartite NLSs) or both minor and major NLS-binding sites (bipartite NLSs). In our recent study, we determined the structure of importin-α1a from rice (Oryza sativa), and made 2 observations that suggest an increased utilization of the minor NLS-binding site in this protein. First, unlike the mammalian protein, both the major and minor NLS-binding sites are auto-inhibited in the unliganded rice protein. Second, we showed that NLSs of the “plant-specific” class preferentially bind to the minor NLS-binding site of rice importin-α. Here, we show that a distinct group of “minor site-specific” NLSs also bind to the minor site of the rice protein. We further show a greater enrichment of proteins containing these “plant-specific” and “minor site-specific” NLSs in the rice proteome. However, the analysis of the distribution of different classes of NLSs in diverse eukaryotes shows that in all organisms, the minor site-specific NLSs are much less prevalent than the classical monopartite and bipartite NLSs.     

Chaperone Proteins Select and Maintain [PIN+] Prion Conformations in Saccharomyces cerevisiae

Prions are proteins that can adopt different infectious conformations known as “strains” or “variants,” each with a distinct, epigenetically inheritable phenotype. Mechanisms by which prion variants are determined remain unclear. Here we use the Saccharomyces cerevisiae prion Rnq1p/[PIN⁺] as a model to investigate the effects of chaperone proteins upon prion variant determination. We show that deletion of specific chaperone genes alters [PIN⁺] variant phenotypes, including [PSI⁺] induction efficiency, Rnq1p aggregate morphology/size and variant dominance. Mating assays demonstrate that gene deletion-induced phenotypic changes are stably inherited in a non-Mendelian manner even after restoration of the deleted gene, confirming that they are due to a bona fide change in the [PIN⁺] variant. Together, our results demonstrate a role for chaperones in regulating the prion variant complement of a cell.     

Similar Replicative Fitness Is Shared by the Subtype B and Unique BF Recombinant HIV-1 Isolates that Dominate the Epidemic in Argentina

The HIV-1 epidemic in South America is dominated by pure subtypes (mostly B and C) and more than 7 BF and BC recombinant forms. In Argentina, circulating recombinant forms (CRFs) comprised of subtypes B and F make up more than 50% of HIV infections. For this study, 28 HIV-1 primary isolates were obtained from patients in Buenos Aires, Argentina and initially classified into subtype B (n = 9, 32.1%), C (n = 1, 3.6%), and CRFs (n = 18, 64.3%) using partial pol and vpu-env sequences, which proved to be inconsistent and inaccurate for these phylogenetic analyses. Near full length genome sequences of these primary HIV-1 isolates revealed that nearly all intersubtype BF recombination sites were unique and countered previous “CRF” B/F classifications. The majority of these Argentinean HIV-1 isolates were CCR5-using but 4 had a dual/mixed tropism as predicted by both phenotypic and genotypic assays. Comparison of the replicative fitness of these BF primary HIV-1 isolates to circulating B, F, and C HIV-1 using pairwise competitions in peripheral blood mononuclear cells (PBMCs) indicated a similarity in fitness of these BF recombinants to subtypes B and F HIV-1 (of the same co-receptor usage) whereas subtype C HIV-1 was significantly less fit than all as previously reported. These results suggest that the multitude of BF HIV-1 strains present within the Argentinean population do not appear to have gained replicative fitness following recent B and F recombination events.     

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.     

Recognition of escape variants in ELISPOT does not always predict CD8+ T-cell recognition of simian immunodeficiency virus-infected cells expressing the same variant sequences

Human immunodeficiency virus (HIV)'s tremendous sequence variability is a major obstacle for the development of cytotoxic-T-lymphocyte-based vaccines, especially since much of this variability is selected for by CD8⁺ T cells. We investigated to what extent reactivity to escape variant peptides in standard enzyme-linked immunospot (ELISPOT) assays predicts the recognition of cells infected with corresponding escape variant viruses. Most of the variant peptides tested were recognized in standard ELISPOT and intracellular cytokine stain (ICS) assays. Functional avidity of epitope-specific T cells for some of the variants was, however, markedly reduced. These mutations which reduced avidity also abrogated recognition by epitope-specific CD8⁺ T cells in a viral suppression assay. Our results indicate that “cross-reactive” CD8⁺ T-cell responses identified in ELISPOT and ICS assays using a single high concentration of variant peptide often fail to predict the recognition of cells infected with variant viruses.