A conifer genome spruces up plant phylogenomics

The Norway spruce genome provides key insights into the evolution of plant genomes, leading to testable new hypotheses about conifer, gymnosperm, and vascular plant evolution.In the past year a burst of plant genome sequences have been published, providing enhanced phylogenetic coverage of green plants (Figure ​(Figure1)1) and inclusion of new agricultural, ecological, and evolutionary models. Collectively, these sequences are revealing some extraordinary structural and evolutionary attributes in plant genomes. Perhaps most surprising is the exceptionally high frequency of whole-genome duplication (WGD): nearly every genome that has been analyzed has borne the signature of one or more WGDs, with particularly notable events having occurred in the common ancestors of seed plants, of angiosperms, and of core eudicots (the latter ''WGD'' represents two WGDs in close succession) [1,2]. Given this tendency for plant genomes to duplicate and then return to an essentially diploid genetic system (an example is the cotton genomes, which have accumulated the effects of perhaps 15 WGDs [3]), the conservation of genomes in terms of gene number, chromosomal organization, and gene content is astonishing. From the publication of the first plant genome, Arabidopsis thaliana [4], the number of inferred genes has been between 25,000 and 30,000, with many gene families shared across all land plants, although the number of members and patterns of expansion and contraction vary. Furthermore, conserved synteny has been detected across the genomes of diverse angiosperms, despite WGDs, diploidization, and millions of years of evolution.Open in a separate windowFigure 1Simplified phylogeny of land plants, showing major clades and their component lineages. Asterisks indicate species (or lineage) for which whole-genome sequence (or sequences) is (are) available. Increases and decreases in genome size are shown by arrows.Despite the proliferation of genome sequences available for angiosperms, genome-level data for both ferns (and their relatives, collectively termed monilophytes; Figure ​Figure1)1) and gymnosperms have been conspicuously lacking - until recently, with the publication of the genome sequence of the gymnosperm Norway spruce (Picea abies) [5]. The large genome sizes for both monilophytes and gymnosperms have discouraged attempts at genome sequencing and assembly, whereas the smaller genome size of angiosperms has resulted in more genome sequences being available (Table ​(Table1)1) [6]. Because of this limited phylogenetic sample, our understanding of the timing and phylogenetic positions of WGDs, the core number of plant genes, possible conserved syntenic regions, and patterns of expansion and contraction of gene families across both tracheophytes (vascular plants) and across all land plants is imperfect. This sampling problem is particularly acute in analyses of the genes and genomes of seed plants; many hundreds of genes are present in angiosperms that are not present in mosses or lycophytes, but whether these genes arose in the common ancestor of seed plants or of angiosperms cannot be determined without a gymnosperm genome sequence. The Norway spruce genome therefore offers tremendous power, not only for understanding the structure and evolution of conifer genomes, but also as a reference for interpreting gene and genome evolution in angiosperms.

Table 1

Genome sizes in land plants

Peptidoglycan Fine Structure of the Radiotolerant Bacterium Deinococcus radiodurans Sark

Peptidoglycan from Deinococcus radiodurans was analyzed by high-performance liquid chromatography and mass spectrometry. The monomeric subunit was: N-acetylglucosamine–N-acetylmuramic acid–l-Ala–d-Glu-(γ)–l-Orn-[(δ)Gly-Gly]–d-Ala–d-Ala. Cross-linkage was mediated by (Gly)₂ bridges, and glycan strands were terminated in (1→6)anhydro-muramic acid residues. Structural relations with the phylogenetically close Thermus thermophilus are discussed.The gram-positive bacterium Deinococcus radiodurans is remarkable because of its extreme resistance to ionizing radiation (14). Phylogenetically the closest relatives of Deinococcus are the extreme thermophiles of the genus Thermus (4, 11). In 16S rRNA phylogenetic trees, the genera Thermus and Deinococcus group together as one of the older branches in bacterial evolution (11). Both microorganisms have complex cell envelopes with outer membranes, S-layers, and ornithine-Gly-containing mureins (7, 12, 19, 20, 22, 23). However, Deinococcus and Thermus differ in their response to the Gram reaction, having positive and negative reactions, respectively (4, 14). The murein structure for Thermus thermophilus HB8 has been recently elucidated (19). Here we report the murein structure of Deinococcus radiodurans with similar detail.D. radiodurans Sark (23) was used in the present study. Cultures were grown in Luria-Bertani medium (13) at 30°C with aeration. Murein was purified and subjected to amino acid and high-performance liquid chromatography (HPLC) analyses as previously described (6, 9, 10, 19). For further analysis muropeptides were purified, lyophilized, and desalted as reported elsewhere (6, 19). Purified muropeptides were subjected to plasma desorption linear time-of-flight mass spectrometry (PDMS) as described previously (1, 5, 16, 19). Positive and negative ion mass spectra were obtained on a short linear ²⁵²californium time-of-flight instrument (BioIon AB, Uppsala, Sweden). The acceleration voltage was between 17 and 19 kV, and spectra were accumulated for 1 to 10 million fission events. Calibration of the mass spectra was done in the positive ion mode with H⁺ and Na⁺ ions and in the negative ion mode with H⁻ and CN⁻ ions. Calculated m/z values are based on average masses.Amino acid analysis of muramidase (Cellosyl; Hoechst, Frankfurt am Main, Germany)-digested sacculi (50 μg) revealed Glu, Orn, Ala, and Gly as the only amino acids in the muramidase-solubilized material. Less than 3% of the total Orn remained in the muramidase-insoluble fraction, indicating an essentially complete solubilization of murein.Muramidase-digested murein samples (200 μg) were analyzed by HPLC as described in reference 19. The muropeptide pattern (Fig. ​(Fig.1)1) was relatively simple, with five dominating components (DR5 and DR10 to DR13 [Fig. 1]). The muropeptides resolved by HPLC were collected, desalted, and subjected to PDMS. The results are presented in Table ​Table11 compared with the m/z values calculated for best-matching muropeptides made up of N-acetylglucosamine (GlucNAc), N-acetylmuramic acid (MurNAc), and the amino acids detected in the murein. The more likely structures are shown in Fig. ​Fig.1.1. According to the m/z values, muropeptides DR1 to DR7 and DR9 were monomers; DR8, DR10, and DR11 were dimers; and DR12 and DR13 were trimers. The best-fitting structures for DR3 to DR8, DR11, and DR13 coincided with muropeptides previously characterized in T. thermophilus HB8 (19) and had identical retention times in comparative HPLC runs. The minor muropeptide DR7 (Fig. ​(Fig.1)1) was the only one detected with a d-Ala–d-Ala dipeptide and most likely represents the basic monomeric subunit. The composition of the major cross-linked species DR11 and DR13 confirmed that cross-linking is mediated by (Gly)₂ bridges, as proposed previously (20). Open in a separate windowFIG. 1HPLC muropeptide elution patterns of murein purified from D. radiodurans. Muramidase-digested murein samples were subjected to HPLC analysis, and the A₂₀₄ of the eluate was recorded. The most likely structures for each muroeptide as deduced by PDMS are shown. The position of residues in brackets is the most likely one as deduced from the structures of other muropeptides but could not be formally demonstrated. R = GlucNac–MurNac–l-Ala–d-Glu-(γ)→.

TABLE 1

Calculated and measured m/z values for the molecular ions of the major muropeptides from D. radiodurans

Cross-Subtype Neutralization Sensitivity despite Monoclonal Antibody Resistance among Early Subtype A,C, and D Envelope Variants of Human Immunodeficiency Virus Type 1

The human immunodeficiency virus type 1 (HIV-1) variants that are transmitted to newly infected individuals are the primary targets of interventions, such as vaccines and microbicides, aimed at preventing new infections. Newly acquired subtype A, B, and C variants have been the focus of neutralization studies, although many of these viruses, particularly of subtypes A and B, represent viruses circulating more than a decade ago. In order to better represent the global diversity of transmitted HIV-1 variants, an additional 31 sexually transmitted Kenyan HIV-1 env genes, representing several recent infections with subtype A, as well as subtypes A/D, C, and D, were cloned, and their neutralization profiles were characterized. Most env variants were resistant to neutralization by the monoclonal antibodies (MAbs) b12, 4E10, 2F5, and 2G12, suggesting that targeting the epitopes of these MAbs may not be effective against variants that are spreading in areas of endemicity. However, significant cross-subtype neutralization by plasma was observed, indicating that there may be other epitopes, not yet defined by the limited available MAbs, which could be recognized more broadly.Most effective viral vaccines are thought to provide protection primarily by stimulating neutralizing antibodies (NAbs) to clear cell-free virus (25, 27). Because protection by NAbs requires recognition of common viral epitopes, the extreme genetic diversity of human immunodeficiency virus type 1 (HIV-1) presents a particular challenge to NAb-based vaccine approaches. Therefore, a critical starting point for studies of immune-mediated protection against HIV-1 is a collection of newly transmitted HIV-1 variants, particularly from areas of endemicity, such as sub-Saharan Africa, in order to determine whether vaccines are appropriately targeted to common epitopes from these relevant transmitted strains.During HIV-1 transmission, a bottleneck allows only one or a few variants to be transmitted to a newly infected individual (6, 9, 16, 29, 34, 37, 39), and the sensitivity of these early transmitted strains to antibody-mediated neutralization is therefore of particular interest. Newly transmitted HIV-1 variants have demonstrated significant heterogeneity in their neutralization phenotypes both within and between subtypes (2, 3, 6-8, 11, 13-15, 22, 30, 32, 36). Panels of sexually transmitted HIV-1 envelope variants (based on the envelope gene, env) have been characterized, including subtype B variants from North America, Trinidad, and Europe, subtype C variants from South Africa and Zambia, and subtype A variants from Kenya collected between 1994 and 1996 (2, 14, 15). Here, we characterize an additional 31 envelope variants from 14 subjects with sexually transmitted HIV-1 who were infected in Kenya, where subtypes A, C, and D circulate, between 1993 and 2005 (24, 31).The env genes were cloned from samples drawn 14 to 391 (median, 65) days postinfection from individuals enrolled in a prospective cohort of high-risk women in Mombasa, Kenya (19-21). Demographic characteristics of the subjects are summarized in Table ​Table1;1; the timing of first infection was determined by both HIV-1 serology and HIV RNA testing as described previously (12). All of the subjects were presumably infected by male-to-female transmission and displayed a range of plasma viral loads at the time of env gene cloning (Table ​(Table1).1). For most individuals, full-length env genes were cloned from uncultured peripheral blood mononuclear cell (PBMC) DNA, though for two individuals, clones were obtained from DNA following short-term coculture with donor PBMCs (Table ​(Table1).1). env genes were cloned by single-copy nested PCR with primers and PCR conditions as described previously (4, 17). We tested env genes for their ability to mediate infection by transfecting env plasmid DNA into 293T cells along with an env-deficient HIV-1 subtype A proviral plasmid, Q23Δenv, to make pseudoviral particles (17). More than 80 env clones were obtained from 16 subjects; less than one-half were functional on the basis of the infectivity of pseudoviral particles in a single-round infection of TZM-bl cells (AIDS Research and Reference Reagent Program, National Institutes of Health), as observed previously for env genes cloned from proviral sequences (17); a lower fraction of functional env genes have been reported from plasma (18). We focused on the proviral sequences here because they presumably best represent the sequence closest to that of the transmitted strains. The 31 functional env variants are described in Table ​Table11.

TABLE 1.

Demographic characteristics, diversities, gp120 variable-region lengths, numbers of PNGS, and accession numbers of cloned env variants

Niemann-Pick disease,type B with TRAP-positive storage cells and secondary sea blue histiocytosis

We present 2 cases of Niemann Pick disease, type B with secondary sea-blue histiocytosis. Strikingly, in both cases the Pick cells were positive for tartrate resistant acid phosphatase, a finding hitherto described only in Gaucher cells. This report highlights the importance of this finding as a potential cytochemical diagnostic pitfall in the diagnosis of Niemann Pick disease.Key words: Niemann pick disease, Gaucher disease, tartrate resistant acid phosphatase, sea blue histiocytosis.We present two unrelated patients who were referred to the Hematology OPD from Gastroenterology during work-up of long-standing splenomegaly 2 years apart and whose details are presented in

Variation in Adult Plant Phenotypes and Partitioning among Seed and Stem-Borne Roots across Brachypodium distachyon Accessions to Exploit in Breeding Cereals for Well-Watered and Drought Environments

Seedling roots enable plant establishment. Their small phenotypes are measured routinely. Adult root systems are relevant to yield and efficiency, but phenotyping is challenging. Root length exceeds the volume of most pots. Field studies measure partial adult root systems through coring or use seedling roots as adult surrogates. Here, we phenotyped 79 diverse lines of the small grass model Brachypodium distachyon to adults in 50-cm-long tubes of soil with irrigation; a subset of 16 lines was droughted. Variation was large (total biomass, ×8; total root length [TRL], ×10; and root mass ratio, ×6), repeatable, and attributable to genetic factors (heritabilities ranged from approximately 50% for root growth to 82% for partitioning phenotypes). Lines were dissected into seed-borne tissues (stem and primary seminal axile roots) and stem-borne tissues (tillers and coleoptile and leaf node axile roots) plus branch roots. All lines developed one seminal root that varied, with branch roots, from 31% to 90% of TRL in the well-watered condition. With drought, 100% of TRL was seminal, regardless of line because nodal roots were almost always inhibited in drying topsoil. Irrigation stimulated nodal roots depending on genotype. Shoot size and tillers correlated positively with roots with irrigation, but partitioning depended on genotype and was plastic with drought. Adult root systems of B. distachyon have genetic variation to exploit to increase cereal yields through genes associated with partitioning among roots and their responsiveness to irrigation. Whole-plant phenotypes could enhance gain for droughted environments because root and shoot traits are coselected.Adult plant root systems are relevant to the size and efficiency of seed yield. They supply water and nutrients for the plant to acquire biomass, which is positively correlated to the harvest index (allocation to seed grain), and the stages of flowering and grain development. Modeling in wheat (Triticum aestivum) suggested that an extra 10 mm of water absorbed by such adult root systems during grain filling resulted in an increase of approximately 500 kg grain ha⁻¹ (Manschadi et al., 2006). This was 25% above the average annual yield of wheat in rain-fed environments of Australia. This number was remarkably close to experimental data obtained in the field in Australia (Kirkegaard et al., 2007). Together, these modeling and field experiments have shown that adult root systems are critical for water absorption and grain yield in cereals, such as wheat, emphasizing the importance of characterizing adult root systems to identify phenotypes for productivity improvements.Most root phenotypes, however, have been described for seedling roots. Seedling roots are essential for plant establishment, and hence, the plant’s potential to set seed. For technical reasons, seedlings are more often screened than adult plants because of the ease of handling smaller plants and the high throughput. Seedling-stage phenotyping may also improve overall reproducibility of results because often, growth media are soil free. Seedling soil-free root phenotyping conditions are well suited to dissecting fine and sensitive mechanisms, such as lateral root initiation (Casimiro et al., 2003; Péret et al., 2009a, 2009b). A number of genes underlying root processes have been identified or characterized using seedlings, notably with the dicotyledonous models Arabidopsis (Arabidopsis thaliana; Mouchel et al., 2004; Fitz Gerald et al., 2006; Yokawa et al., 2013) and Medicago truncatula (Laffont et al., 2010) and the cereals maize (Zea mays; Hochholdinger et al., 2001) and rice (Oryza sativa; Inukai et al., 2005; Kitomi et al., 2008).Extrapolation from seedling to adult root systems presents major questions (Hochholdinger and Zimmermann, 2008; Chochois et al., 2012; Rich and Watt, 2013). Are phenotypes in seedling roots present in adult roots given developmental events associated with aging? Is expression of phenotypes correlated in seedling and adult roots if time compounds effects of growth rates and growth conditions on roots? Watt et al. (2013) showed in wheat seedlings that root traits in the laboratory and field correlated positively but that neither correlated with adult root traits in the field. Factors between seedling and adult roots seemed to be differences in developmental stage and the time that growing roots experience the environment.Seedling and adult root differences may be larger in grasses than dicotyledons. Grass root systems have two developmental components: seed-borne (seminal) roots, of which a number emerge at germination and continue to grow and branch throughout the plant life, and stem-borne (nodal or adventitious) roots, which emerge from around the three-leaf stage and continue to emerge, grow, and branch throughout the plant life. Phenotypes and traits of adult root systems of grasses, which include the major cereal crops wheat, rice, and maize, are difficult to predict in seedling screens and ideally identified from adult root systems first (Gamuyao et al., 2012).Phenotyping of adult roots is possible in the field using trenches (Maeght et al., 2013) or coring (Wasson et al., 2014). A portion of the root system is captured with these methods. Alternatively, entire adult root systems can be contained within pots dug into the ground before sowing. These need to be large; field wheat roots, for example, can reach depths greater than 1.5 m depending on genotype and environment. This method prevents root-root interactions that occur under normal field sowing of a plant canopy and is also a compromise.A solution to the problem of phenotyping adult cereal root systems is a model for monocotyledon grasses: Brachypodium distachyon. B. distachyon is a small-stature grass with a small genome that is fully sequenced (Vogel et al., 2010). It has molecular tools equivalent to those available in Arabidopsis (Draper et al., 2001; Brkljacic et al., 2011; Mur et al., 2011). The root system of B. distachyon reference line Bd21 is more similar to wheat than other model and crop grasses (Watt et al., 2009). It has a seed-borne primary seminal root (PSR) that emerges from the embryo at seed germination and multiple stem-borne coleoptile node axile roots (CNRs) and leaf node axile roots (LNRs), also known as crown roots or adventitious roots, that emerge at about three leaves through to grain development. Branch roots emerge from all root types. There are no known anatomical differences between root types of wheat and B. distachyon (Watt et al., 2009). In a recent study, we report postflowering root growth in B. distachyon line Bd21-3, showing that this model can be used to answer questions relevant to the adult root systems of grasses (Chochois et al., 2012).In this study, we used B. distachyon to identify adult plant phenotypes related to the partitioning among seed-borne and stem-borne shoots and roots for the genetic improvement of well-watered and droughted cereals (Fig. 1; Krassovsky, 1926; Navara et al., 1994), nitrogen, phosphorus (Tennant, 1976; Brady et al., 1995), oxygen (Wiengweera and Greenway, 2004), soil hardness (Acuna et al., 2007), and microorganisms (Sivasithamparam et al., 1978). Of note is the study by Krassovsky (1926), which was the first, to our knowledge, to show differences in function related to water. Krassovsky (1926) showed that seminal roots of wheat absorbed almost 2 times the water as nodal roots per unit dry weight but that nodal roots absorbed a more diluted nutrient solution than seminal roots. Krassovsky (1926) also showed by removing seminal or nodal roots as they emerged that “seminal roots serve the main stem, while nodal roots serve the tillers” (Krassovsky, 1926). Volkmar (1997) showed, more recently, in wheat that nodal and seminal roots may sense and respond to drought differently. In millet (Pennisetum glaucum) and sorghum (Sorghum bicolor), Rostamza et al. (2013) found that millet was able to grow nodal roots in a dryer soil than sorghum, possibly because of shoot and root vigor.Open in a separate windowFigure 1.B. distachyon plant scanned at the fourth leaf stage, with the root and shoot phenotypes studied indicated. Supplemental Table S1.

Clade-Specific Evolution Mediated by HLA-B*57/5801 in Human Immunodeficiency Virus Type 1 Clade A1 p24

HLA-B*57-mediated selection pressure leads to a typical escape pathway in human immunodeficiency virus type 1 (HIV-1) CD8 epitopes such as TW10. Whether this T242N pathway is shared by all clades remains unknown. We therefore assessed the nature of HLA-B*57 selection in a large, observational Kenyan cohort where clades A1 and D predominate. While T242N was ubiquitous in clade D HLA-B*57⁺ subjects, this mutation was rare (15%) in clade A1. Instead, P243T and I247L were selected by clade A1-infected HLA-B*57 subjects but not by HLA-B*5801⁺ subjects. Our data suggest that clade A1 consensus proline at Gag residue 243 might represent an inherent block to T242N escape in clade A1. We confirmed immunologically that P243T and I247L likely represent escape mutations. HLA-B*57 evolution also differed between clades in the KF11 and IW9 epitopes. A better understanding of clade-specific evolution is important for the development of HIV vaccines in regions with multiple clades.Human immunodeficiency virus type 1 (HIV-1) displays extreme genetic diversity, with nine clades (subtypes) described in group M, and frequent genomic recombination among and within the clades (7, 44). HIV is also capable of rapid evolution, which can lead to mutational escape from immune control (43). Escape from CD8⁺ T-cell responses occurs frequently in HIV-1 infection through mutations that affect epitope processing, HLA class I binding, and/or T-cell receptor recognition (23). In early HIV-1 infection, the majority of amino acid substitutions are associated HLA class I alleles (1). The timing and consequences of mutational escape from CD8⁺ T-cell responses vary considerably (8, 22).HLA-B*57, and to a lesser extent HLA-B*5801, has been associated with slower progression to AIDS in several studies (18, 27, 39), and HLA-B*5701 was associated with a lower viremia set point in a genome-wide association study (16). Several attributes of HLA-B*57-restricted CD8⁺ T-cell responses may contribute to their protectiveness, including dominant responses in acute infection (2), recognition of protective epitopes in HIV-1 p24 (33), better recognition of epitope variation (45), and retention of proliferative capability in chronic infection (24).HLA-B*57/5801 also exert powerful selection pressure on HIV to avoid CD8⁺ T-cell recognition. This was first demonstrated in the HLA-B*57-restricted TW10 epitope (TSTLQEQIGW [Gag_240-249]), which accounts for >30% of overall HIV-specific CD8⁺ T-cell responses in acutely infected HLA-B*57⁺ subjects (3). Escape in this epitope usually occurs early in infection, which coincidently is when HLA-B*57 is most protective (18). In clade B and C infections, >75 to 100% of HLA-B*57/5801⁺ subjects develop the T242N escape mutation, while HLA-B*57/5701-negative subjects rarely display polymorphism at this residue (5, 9, 10, 15, 32, 35, 38, 41). When T242N is transmitted to HLA-B*57/5801-negative subjects, it rapidly reverts to the consensus, suggesting that T242N is associated with a fitness defect (32, 35).While CD8⁺ T-cell cross-clade recognition has been tested extensively (6, 11, 19, 36, 48), few studies have addressed the possibility of clade-specific escape from CD8⁺ T-cell responses. This may be especially relevant where clade consensus sequences differ in immunologically relevant epitopes. Here we demonstrate in a large Kenyan cohort substantial differences in HLA-B*57/B*5801-mediated selection among HIV clades.Participants were enrolled from a Nairobi, Kenya-based cohort, and the relevant ethical review boards approved the study. HLA typing was performed as described previously (34). CD4 counts were measured longitudinally at biannual visits. Multiple and other clade infections were excluded. The HIV-1 p24 gene was amplified from proviral HIV DNA or RNA using a nested PCR approach and sequenced, and viral subtyping was carried out as described previously (42). Previously described HLA-B*57 epitopes IW9 (ISPRTLNAW), KF11 (KAFSPEVIPMF), and TW10 (TSTLQEQIGW) and selected variants were tested in immunological assays and described where relevant. Gamma interferon enzyme-linked immunospot (ELISPOT) assays were performed as described previously (37) using blood samples from HLA-B*57⁺ and -B*5801⁺ subjects. All peptides were tested at concentrations of 10 μg, 1 μg, 0.1 μg, and 0.01 μg/ml. Responses were considered positive if they were more than two times higher than that of the negative control and were measured at ≥100 spot-forming units ml⁻¹. Fisher''s exact test and chi-square analyses were used to determine differences among groups in categorical analyses. Mann-Whitney U tests were used to compare response magnitudes and disease progression between groups.We confirmed the protective effects of HLA-B*57 in clade A1 infection (mean of 9.9 years versus 7.8 years until CD4 counts were <200, P = 0.041) (Fig. ​(Fig.1).1). Slow progressors were overrepresented in HLA-B*57⁺ clade A1⁺ subjects (52.2%) compared to both HLA-B*5801⁺ clade A1⁺ (13.3%, P = 0.02) and HLA-B*57/5801-negative clade A1⁺ (27.8%, P = 0.028) subjects (Fig. ​(Fig.1b).1b). In contrast to what has been shown for other clades (2, 27), protection was not observed for clade A1-infected HLA-B*5801⁺ subjects (mean of 6.5 years versus 7.8 years until CD4 counts were <200, P > 0.3) (Fig. ​(Fig.11).Open in a separate windowFIG. 1.HLA-B*57, but not HLA-B*5801, is associated with a lower rate of disease progression in clade A1-infected subjects than that of the overall cohort. (a) Number of years from cohort entry until sequential CD4 counts fell below 200/μl. (b) Slow progressors (>10 years with CD4 counts of >200) were also more common in HLA-B*57⁺ clade A1-infected subjects than in those expressing HLA-B*5801 or neither.Stratification of TW10 (Gag_240-249) proviral sequences on the basis of HLA allele and clade revealed several differences in selection between clades A1 and D (Fig. ​(Fig.2a).2a). We observed the expected T242N substitution in 100% of HLA-B*57⁺ clade D-infected subjects (7/7), compared to only 14.7% variability at Gag residue 242 in HLA-B*57/5801-negative subjects (13/88, P = 3.26 × 10⁻⁹) (Fig. ​(Fig.2b).2b). In contrast, T242N was found infrequently in clade A1-infected HLA-B*57⁺ subjects (15%, 5/33, P = 0.0004). Instead, variants containing the mutations P243T and I247L were more frequently observed (both observed in 11/33 subjects). Overall, variation at residues 243 and 247 was more common in HLA-B*57⁺ subjects (51% and 15%, respectively; P = 2.92 × 10⁻⁶) than in HLA-B*57/5801-negative clade A1⁺ subjects (33% and 9%, respectively; P = 0.0008). Selection at both residues 243 and 247 was observed only in 2/33 HLA-B*57⁺ subjects, suggesting that these mutations are independent. Selection at residue 248, observed in clade B infection (32), was not evident in either clade A1 or D. While I247X selection has been described in other clades at low frequencies and in elite controllers (21, 40), HLA-B*57-mediated selection at Gag residue 243 has not yet been described. In summary, the T242N mutation, which is typical of other clades, does not appear to be the primary escape mutant in clade A1.Open in a separate windowFIG. 2.HLA-B*57-mediated selection in TW10 differs between clade A1 and clade D. (a) TW10 sequences were stratified by HLA-B*57, HLA-B*5801, or other alleles (HLA-B*57/5801⁻) and compared between clades A1 and D, based on the clade B consensus TW10 sequence. Each subject is represented by one sequence, and the numbers of subjects with a given sequence are shown in parentheses. A summary of variation from the TW10 consensus at Gag residues 242, 243, 247, and 248 is shown for HLA-B*57⁺ (b) and -B*5801⁺ (c) subjects. (b and c) Clade D is shown at the top, and clade A1 is shown at the bottom. Variation is shown in dark gray, and consensus is shown in light gray. (d) The proportions of clade A1-infected subjects with selection at Gag residue 242 only, and those with selection at residues 242 and 243 in combination, are shown. *, P < 0.05; **, P < 0.005; ***, P < 0.0005.Previous studies have suggested that HLA-B*5801 places selection pressure on TW10, similar to that of HLA-B*57 (35). Similar to clades B and C, selection of T242N was evident in HLA-B*5801⁺ clade D-infected subjects (TW10 variation in 8/11 HLA-B*5801⁺ subjects versus 13/88 HLA-B*57/5801-negative subjects; P = 0.0069) (Fig. ​(Fig.2c).2c). Limited T242N selection was observed in clade A1-infected HLA-B*5801⁺ subjects, and in contrast to HLA-B*57, there were no HLA-B*5801-associated substitutions at residues 243 and 247 in clade A1 (P values of 0.75 and 0.29, respectively) (Fig. ​(Fig.2c).2c). In summary, these data suggest that in addition to HLA-B*5801 not being associated with protection in clade A1 (Fig. ​(Fig.1),1), HLA-B*5801 does not select the HLA-B*57-associated clade A1 TW10 escape mutations.Inclusion of all clade A1 sequences with T242X substitutions (regardless of the HLA allele) reveals that in every case (10/10), there is an accompanying residue 243 mutation. Polymorphisms at these sites correlate very strongly (P = 3.71 × 10⁻⁸) (Fig. ​(Fig.2d).2d). Together, these data suggest that residue 242 polymorphism in clade A1 is incompatible with proline at residue 243, which is the clade A1 consensus.We next assessed the immunological implications of novel clade A1 variants in HLA-B*57⁺ (n = 12) and -B*5801⁺ (n = 6) subjects infected primarily by clade A1. Clade A1-infected subjects commonly made anamnestic, low-avidity responses to TW10. The majority of HLA-B*57⁺ subjects who recognized clade A1 TW10 did not respond to P243T or I247L in ELISPOT assays (Fig. ​(Fig.3),3), supporting the hypothesis that these represent escape mutations. Those who did recognize P243T and I247L had lower magnitude responses than those who recognized clade A1 TW10 at the 10-μg/ml peptide concentration (P of 0.0005 for both) (Fig. ​(Fig.3a).3a). Similarly, these variants were not well recognized by CD8⁺ T cells from HLA-B*5801⁺ subjects (Fig. ​(Fig.3b).3b). For two HLA-B*57⁺ subjects, P243T and I247L responses had lower avidity than clade A1 TW10 responses (Fig. ​(Fig.3c).3c). These data support the hypothesis that P243T and I247L likely represent escape mutations.Open in a separate windowFIG. 3.Peptides with novel TW10 clade A1-selected mutations are poorly recognized in ex vivo gamma interferon ELISPOT avidity assays, suggestive of escape mutations. ELISPOT responses to TW10 and variants at 10 μg/ml peptide by HLA-B*57⁺ (a) and HLA-B*5801⁺ (b) subjects are shown. (c) The functional avidity of TW10 and variants for two HLA-B*57⁺ subjects is shown, suggesting that P243T and I247L are less recognized than TW10, particularly at lower peptide concentrations. Sequence names are described in the text. (d) Elispot responses to A1 TW10 and T242N correlated at 10 μg/ml. SFU/m, spot-forming units/million PBMCs.Recognition of the clade B/D consensus (TSTLQEQIGW) was diminished compared to that of clade A1 TW10. However, despite the presumed absence of this variant in these subjects'' autologous sequences, the clade B/D escape variant (TSNLQEQIGW [T242N]) was recognized at a magnitude similar to that of the consensus clade A1 TW10 (r = 0.71, P = 0.0099) (Fig. ​(Fig.3d).3d). No responses to T242N/G248A were observed (not shown), as described previously (32). These data suggest that clade A1 and B TW10, and their escape variants, are immunologically distinct from one another.We next assessed whether clade-specific selection was evident in other immunodominant HLA-B*57 p24 epitopes that are commonly targeted in chronic clade B infection (2). In clade D IW9 (ISPRTLNAW [Gag_147-155]), variants containing the escape variant I147L (14) were more common in HLA-B*57⁺ subjects than in HLA-B*57/5801-negative subjects (86% and 30%, respectively; P = 0.0055) (Table ​(Table1).1). However, this variant was not selected in clade A1 (variation in 30% versus 21% subjects; P value was not significant), where leucine is the consensus. Interestingly, ELISPOT data indicated substantial cross-reactivity between 147L and 147I in clade A1-infected subjects (10 μg/ml, r = 0.987, P < 0.0001) (data not shown), suggesting that infection with an escape variant from one clade (clade D) does not necessarily preclude recognition of this epitope in another one (clade A1). Other amino acids (F, M, and P) were common in HLA-B*57⁺ subjects at residue 147 (>30% versus 3% in HLA-B*57/5801-negative subjects, P = 3.85 × 10⁻⁶). Although the immunological consequences are unknown, these HLA-associated substitutions could represent novel escape variants.

TABLE 1.

p24 sequences in HLA-B*57⁺ subjects infected by clades A1 and D