7
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: not found
      • Article: not found

      Cross-reactive public TCR sequences undergo positive selection in the human thymic repertoire

      Read this article at

      ScienceOpenPublisherPMC
      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          <p class="first" id="d5279657e314">We studied human T cell repertoire formation using high-throughput T cell receptor β (TCRβ) complementarity-determining region 3 (CDR3) sequencing in immunodeficient mice receiving human hematopoietic stem cells (HSCs) and human thymus grafts. Replicate humanized mice generated diverse and highly divergent repertoires. We observed repertoire narrowing and increased CDR3β sharing during thymocyte selection. Whereas hydrophobicity analysis implicated self-peptides in positive selection of the overall repertoire, positive selection favored shorter shared sequences that had reduced hydrophobicity at positions 6 and 7 of CDR3βs, suggesting weaker interactions with self-peptides than were observed with unshared sequences, possibly allowing escape from negative selection. Sharing was similar between autologous and allogeneic thymi and occurred between different cell subsets. Shared sequences were enriched for allo–cross-reactive CDR3βs and for type 1 diabetes–associated autoreactive CDR3βs. Single-cell TCR sequencing showed increased sharing of CDR3αs compared with CDR3βs between mice. Our data collectively implicate preferential positive selection for shared human CDR3βs that are highly cross-reactive. Although previous studies suggested a role for recombination bias in producing “public” sequences in mice, our study is the first to our knowledge to demonstrate a role for thymic selection. Our results implicate positive selection for promiscuous TCRβ sequences that probably evade negative selection, given their low affinity for self-ligands, in the abundance of “public” human TCRβ sequences. </p>

          Related collections

          Most cited references51

          • Record: found
          • Abstract: found
          • Article: not found

          Membrane protein folding and stability: physical principles.

          Stably folded membrane proteins reside in a free energy minimum determined by the interactions of the peptide chains with each other, the lipid bilayer hydrocarbon core, the bilayer interface, and with water. The prediction of three-dimensional structure from sequence requires a detailed understanding of these interactions. Progress toward this objective is summarized in this review by means of a thermodynamic framework for describing membrane protein folding and stability. The framework includes a coherent thermodynamic formalism for determining and describing the energetics of peptide-bilayer interactions and a review of the properties of the environment of membrane proteins--the bilayer milieu. Using a four-step thermodynamic cycle as a guide, advances in three main aspects of membrane protein folding energetics are discussed: protein binding and folding in bilayer interfaces, transmembrane helix insertion, and helix-helix interactions. The concepts of membrane protein stability that emerge provide insights to fundamental issues of protein folding.
            Bookmark
            • Record: found
            • Abstract: found
            • Article: found
            Is Open Access

            A Single Autoimmune T Cell Receptor Recognizes More Than a Million Different Peptides*

            Introduction The mammalian T cell receptor (TCR) 6 orchestrates immune responses by binding to foreign peptides presented at the cell surface in the context of major histocompatibility complex (MHC) molecules. Recognition is mediated by the highly variable complementarity-determining regions of the αβ TCR (1, 2). A priori, the TCR repertoire must be broad enough to respond to all foreign peptides that can bind to self-MHC molecules (3). If this were not the case, then pathogens could rapidly evolve to exploit such deficiencies in immune coverage. Current estimates of human αβ TCR diversity suggest that there are 1015 distinct peptide-MHCs). This represents a particular biochemical challenge to the immune system because the TCR, unlike the B cell receptor, cannot undergo affinity maturation in the form of somatic hypermutation. It is unclear how the limited naïve T cell pool responds to a multitude of ligands that it has never encountered before and cannot adapt to at the protein sequence level. One possibility is that T cell immunity inherently features an extremely high level of receptor degeneracy, enabling each TCR to recognize multiple peptides. However, clonal selection theory suggests that individual T cells are specific for a single peptide-MHC molecule with recognition of alternative ligands unlikely. In contrast, studies published in the 1990s demonstrated that T cells can recognize several peptides (5–11). Since then, observations of TCR degeneracy have continued to accumulate in the literature (12–19). In addition, other studies have shown that TCRs can recognize distinct peptides in the context of non-self-MHC, a phenomenon known as alloreactivity (20–25). The majority of previous studies of TCR degeneracy have examined sets of between 2 and 200 peptides, with one recent study examining ∼4,000 peptides (19). Given that the entire universe of decamer peptides alone comprises >1013 distinct amino acid sequences, the proportion of the peptide universe at this length that has been examined in the most comprehensive study to date (19) remains extremely small ( one million distinct peptides. As such, the TCR represents one of the most remarkable biological receptors and by far the most promiscuous known. EXPERIMENTAL PROCEDURES Generation and Maintenance of an Autoimmune CD8+ T Cell Clone The 1E6 CD8+ T cell clone specific for the human leukocyte antigen (HLA) A*0201-restricted autoantigen preproinsulin peptide ALWGPDPAAA (PPI15–24) was generated as described previously (26). Decamer Combinatorial Peptide Library (CPL) Scan The decamer combinatorial peptide library (Pepscan) contains a total of 9.36 × 1012 (= (10 + 19) ×199) different decamer peptides and is divided into 200 different peptide mixtures (see Fig. 1). In every peptide mixture, one position has a fixed l-amino acid residue and all other positions are degenerate, with the possibility of any one of 19 natural l-amino acids being incorporated in each individual position (cysteine is excluded). Each library mixture consists of 3.2 × 1011 (199) different decamer peptides in approximately equimolar concentrations. For CPL screening, 1E6 CD8+ T cells were washed and rested overnight in RPMI 1640 medium containing 100 units/ml penicillin, 100 μg/ml streptomycin, 2 mm l-glutamine, and 2% heat-inactivated fetal calf serum (all Invitrogen). In 96-well U-bottom plates, 6 × 104 C1R-A2 cells were incubated with various peptide library mixtures (at 100 μg/ml) in duplicate for 2 h at 37 °C. Following peptide pulsing, 3 × 104 1E6 CD8+ T cells were added, and the assay was incubated overnight at 37 °C. Subsequently, the supernatant was harvested and assayed for MIP1β by ELISA according to the manufacturer's instructions (R&D Systems). CD8+ T Cell Effector Function Assays: MIP1β ELISA Individual peptides were assayed for agonist activity by MIP1β ELISA as described above. Functional sensitivity is expressed by the pEC50 of each peptide with respect to the TCR. This is defined as −1 × the base 10 logarithm (p) of the 50% efficacy concentration (EC50); a greater functional sensitivity is indicated by a larger pEC50 value, which was estimated as described in supplemental Equations. Overview of Sampling Approaches Used to Quantify TCR Degeneracy Although the TCR has an appreciable degeneracy, it is still specific enough that sampling peptides at random will most likely result in less than ∼10 strong agonists for every ∼10,000 peptides sampled. For this reason, we employed conditioned sampling. In one approach, we sampled from motif-restricted peptide sets. This results in lower bound estimates of the actual number of agonist ligands, in that any agonist not fitting the prescribed motif is excluded from the sample. In general, more stringent motifs exclude more agonists but provide better resolution at the high pEC50 end of the curve, and vice versa; for this reason, a range of motifs of varying stringency was used. In a second approach (CPL-based importance sampling), we sampled from the entire peptide universe with bias toward peptides that were likely to elicit a response, then estimated a true distribution by applying a correction weighting to the observations (i.e. dividing back by the bias). Sampling Equations Motif-restricted Sampling A sampling motif specifies, for each of the m positions in an m-mer peptide, one or more amino acid residues that may occur at that position. If np alternative residues have been specified at position p, the probability that a sampled peptide has a given residue at position p equals 1/np if the given residue is one of the np given alternatives, and zero otherwise. Consider a sample of n peptides, consisting of peptides P [1], P [2], …, P [n], with measured functional sensitivities pEC50 [1], pEC50 [2], …, pEC50 [n]. The cross-reactivity of the TCR is expressed by the number of peptides whose pEC50 exceeds a given value ω. The degeneracy of the TCR is then represented by determining this number for a range of ω values, estimated as follows. Here, N[pEC50 > ω] is the number of peptides that have a pEC50 for the TCR that exceeds ω, and if sample peptide j has a pEC50 that exceeds ω; otherwise. The sampling motif typically excludes a number of agonists (unless np = 20 at every position, i.e. the “universal motif”, which was not used). This exclusion means that the quantity on the right in Equation 1 always underestimates the true N[pEC50 > ω]. This means that the motif-based method provides a lower boundary to TCR degeneracy, i.e. a conservative estimate. CPL-based Importance Sampling The idea of sampling motifs can be generalized to the well known strategy of importance sampling by specifying, for all positions and amino acids, the probability that a given amino acid occurs at a given position. Then, the probability of drawing a peptide P that has amino acid residue P(p) at position p is given by where ℙ[P(p)] is the probability of drawing amino acid P at position p. Standard rules of probability stipulate that for each fixed p, the ℙ[P(p)] sum to unity over the 20 amino acids. The CPL was used to generate distributions that could be expected to bias the sample toward good agonists. However, to correct for the bias in the estimate, observations must be weighted by the reciprocal of the bias. Accordingly, degeneracy was estimated as follows: where H is the entropy of the sampling distribution, defined by The factor n exp{H} in Equation 5 is an estimate of the effective sample size. This estimate is based on a standard result concerning the probability of sampling a fixed set of n items. A well known theorem (27) states that this probability is close to exp{−n H}. On the other hand, in unbiased sampling from a set of size N, the probability of obtaining any given n-element sample is exactly equal to N − n . Together, these observations indicate that exp{H} can be interpreted as the “diameter” of the population from which the peptide is drawn. Consequently, an n-element sample of m-mer peptides probes an effective set size of the approximation is accurate when n is tiny compared with exp{−H}, as is the case in the experiments reported here. RESULTS CPL Screening Reveals the Potential for TCR Degeneracy To determine the extent of T cell cross-reactivity, we probed the peptide recognition degeneracy of the autoimmune CD8+ T cell clone IE6 using a CPL comprising 9.36 × 1012 different decamer peptides (Fig. 1). Using this approach, we were able to scan every amino acid at every position of the peptide within a random residue “backbone” and build a detailed picture of the molecular landscape preferred by the 1E6 TCR. The 1E6 clone was generated from a patient with type 1 diabetes and is the only documented example of an autoreactive CD8+ T cell that can kill human pancreatic islet β-cells (26). Killing is mediated via recognition of residues 15–24 (ALWGPDPAAA) of the autoantigen preproinsulin bound to HLA A*0201 on the β-cell surface. The majority of HLA A*0201+ patients with type I diabetes recognize the preproinsulin 15–24 epitope (26), and HLA A*0201 is known to confer an increased risk of disease (28, 29). The number of amino acids that were recognized by the 1E6 clone was restricted in the central region of the peptide (residues 4–6), suggesting that this TCR makes the majority of its peptide contacts with these residues. In contrast, recognition was highly degenerate at the remaining positions. This degeneracy was confirmed by the ability of 1E6 T cells to recognize a panel of peptides robustly with any of the 20 natural proteogenic l-amino acids at peptide position 8 (supplemental Fig. S1), with half of the substitutions leading to increased levels of functional sensitivity. Similar results were obtained with corresponding scans at other degenerate positions (data not shown). The CPL scan results also revealed that the index peptide is suboptimal in all positions outside the central region (residues 4–6). Thus, positional peptide degeneracy is extreme at 7 of 10 positions, hinting at the potential for a single TCR to recognize a multitude of different amino acid combinations. FIGURE 1. Decamer CPL scan of the 1E6 CD8+ T cell clone. A, 6 × 104 C1R-A2 cells were pulsed in duplicate with each mixture from a decamer CPL (100 μg/ml) at 37 °C. After 2 h, 3 × 104 1E6 CD8+ T cells were added and incubated overnight. Supernatant was harvested and assayed for MIP1β. B, data from A are displayed as a box plot summary. The index insulin peptide sequence is shown below the boxes in black. Quantifying the Number of Decamer Peptides That Can Be Recognized by a Single TCR We next sought to quantify the number of decamer peptides that can be recognized by the 1E6 clone. The total number of decamers that can be made by combinations of the 20 natural proteogenic l-amino acids is 2010 (1.02 × 1013). Only a small proportion (∼1–3%) of all peptides are predicted to bind any given MHC (3, 19), although our own experiments predict that this percentage could be far greater for HLA A*0201 (data not shown). Even with the most conservative estimates of MHC binding (19), the number of potential antigenic HLA A*0201-restricted decamer peptides is still extremely large (1.02 × 1011) and precludes screening all possibilities in T cell recognition assays. To overcome this problem, we screened sets of 30 peptides sampled (Mathematica®; supplemental Fig. S5) from larger motif-restricted or CPL-based importance-sampled sets differing in total size from 225 to 1.66 × 108 individual peptides as described below. Quantifying T Cell Cross-reactivity Using Motif-restricted Sampling Motif-restricted peptide sets were designed based on CPL evidence for amino acid preference at each position of the peptide. First, we screened 30 peptides sampled at random from a total set size of 225 (Motif I, RQWGPDP{A/C/D/F/H/I/K/L/M/N/P/R/S/V/Y}{A/C/G/H/I/K/L/M/N/P/Q/R/S/T/V}A; Fig. 2 and supplemental Fig. S2). Values of pEC50 (−1 × the base 10 logarithm of the EC50) as a measure of functional sensitivity were estimated for all peptides using simultaneous curve fitting (as described in supplemental Equation S1, Fig. S6, and Table S3). Accordingly, increases in functional sensitivity translate into increases in the pEC50 value. The 1E6 clone recognized all peptides within this subset efficiently, with 24 of 30 peptides eliciting greater levels of functional sensitivity than the index peptide. A further 30 peptides were sampled at random from a total set size of 5,776 (Motif II, RQWGP{D/F}{P/F}XX{A/I/L/V}, where X denotes any one of the amino acids excluding cysteine (Fig. 3 and supplemental Fig. S2 and Table S3). One peptide from this subset was recognized poorly; 16 were recognized with pEC50 >7, with a total of 8 peptides recognized more efficiently than the index peptide (Fig. 3 and supplemental Fig. S2). A further two motif-restricted sets of increasing degeneracy were screened (Motif III, RQXGPDXXXA, total set size 194; and Motif IV, XQXGPDXXXV, total set size 195; X denotes any one of the amino acids excluding cysteine; Figs. 4 and 5 and supplemental Fig. S3 and Table S3). The extent of peptide recognition when peptides are sampled at random from large subsets demonstrates the considerable degree of degeneracy exhibited by the clonal 1E6 TCR. FIGURE 2. Recognition of 30 peptides sampled at random from a large peptide set (motif, RQWGPDP{A/C/D/F/H/I/K/L/M/N/P/R/S/V/Y}{A/C/G/H/I/K/L/M/N/P/Q/R/S/T/V}A; total set size = 225). 6 × 104 C1R-A2 cells were pulsed with peptides at various concentrations. After 2 h, 3 × 104 1E6 CD8+ T cells were added and incubated overnight. Supernatant was harvested and assayed for MIP1β. Each panel displays titrations of five different peptides relative to index. A, titrations of peptides with the highest functional sensitivities. F, titrations of peptides with the lowest functional sensitivities. Error bars, S.D. from the mean of two replicates. pEC50 values for each peptide are displayed in supplemental Table S3. FIGURE 3. Recognition of 30 peptides sampled at random from a large peptide set (motif,RQWGP{D/F}{P/F}XX{A/I/L/V}; total set size = 5,776). Details are as described for Fig. 2. FIGURE 4. Recognition of 30 peptides sampled at random from a large peptide set (motif, RQXGPDXXXA; total set size = 194). Details are as described for Fig. 2. FIGURE 5. Recognition of 30 peptides sampled at random from a large peptide set (motif, XQXGPDXXXV; total set size = 195). Details are as described for Fig. 2. Quantifying T Cell Cross-reactivity Using CPL-based Importance Sampling A limitation of sampling from motif-restricted sets is that strong ligands will invariably be excluded, resulting in an underestimate of the true number of agonists (i.e. a lower bound estimate of the number of different peptides that a single TCR can recognize). To obtain a more accurate estimate of true T cell cross-reactivity, we employed CPL-based importance sampling, which makes no assumptions about TCR contact or MHC class I-binding residues. Importance sampling ensures that every peptide has a chance of being sampled (although cysteine was excluded to avoid the potential for oxidation), but incorporates bias toward strong agonists predicted using CPL scan data. This bias is adjusted for to yield unbiased estimates of agonist numbers (see under “Sampling Equations”). Raw data from the primary CPL scan were modified as described in supplemental Table S1 and subsequently normalized (supplemental Table S2) to provide a peptide sampling distribution biasing the sample toward good agonists. The chance of picking a peptide is the product of the normalized weights assigned to each amino acid residue at each given peptide position. Two sets of 30 peptides were drawn from an effective set size of 1.66 × 108 (calculated from the sampling entropy; see under “Sampling Equations,” Equation 7). Of a total of 60 peptides, 34 were recognized efficiently with a pEC50 > 7, and 22 were better agonists than the index peptide (Figs. 6 and 7 and supplemental Fig. S4 and Table S3). Interestingly, just this set of 60 peptides identified four peptides with functional sensitivities 100-fold better than the index peptide, despite differing from the index peptide sequence (ALWGPDPAAA) at 6 (YQFGPDFPIA, KQFGPDFPTA) or 7 (RQFGPDFPKL, RQFGPDFPTI) positions (Figs. 6 and 7 and supplemental Fig. S4 and Table S3). Thus, CPL-based importance sampling demonstrates that a large proportion of peptides from a biased set of 1.66 × 108 peptides would be recognized and that these recognized peptides can differ considerably from the index peptide. FIGURE 6. Recognition of 30 peptides drawn from a CPL-based importance sampling set with effective size = 1. 66 × 108 (calculated from the sampling entropy) (first set). Details are as described for Fig. 2. FIGURE 7. Recognition of 30 peptides drawn from a CPL-based importance sampling set with effective size = 1. 66 × 108 (calculated from the sampling entropy) (second set). Details are as described for Fig. 2. A Single TCR Can Recognize More Than One Million Different Peptides The pEC50 expresses the potency of a ligand, often referred to as “functional avidity” (supplemental Fig. S6). Cross-reactivity can be quantified precisely by specifying the number of ligands that the TCR recognizes with a pEC50 of at least a given value ω. This number decreases as ω increases; an insight into the nature of TCR degeneracy is afforded by plotting this agonist number as a function of ω. Estimation of this number was performed using Equations 1 and 5, resulting in Fig. 8. The motif-restricted estimate (solid lines in Fig. 8) is a lower boundary, which becomes tighter as the degeneracy of the motif increases; however, this advantage is offset by the decreasing chances of finding good agonists in a sample of 30 peptides. The estimates based on the biased sampling (dashed lines in Fig. 8) indicate that in the order of one million agonists exist for 1E6 that are as least as good as the index peptide. For comparison, the curve derived from TCR activation theory (30, 31) is also shown (gray dashed line in Fig. 8; supplemental Equation S2). FIGURE 8. 1E6 CD8+ T cells can recognize more than one million different decamer peptides. I, RQWGPDP{A/C/D/F/H/I/K/L/M/N/P/R/S/V/Y}{A/C/G/H/I/K/L/M/N/P/Q/R/S/T/V}A (set size 225; 30 peptides sampled at random). II, RQWGP{D/F}{P/F}XX{A/I/L/V} (set size 5,776; 30 peptides sampled at random). III, RQXGPDXXXA (set size 194; 30 peptides sampled at random). IV, XQXGPDXXXV (set size 195; 30 peptides sampled at random). In the motifs, X denotes any one of the 19 amino acids excluding cysteine. Va and Vb, two replicates of a biased sampling set (effective set size 1.66 × 108, calculated from the sampling entropy); each set of 30 peptides was sampled with bias toward strong agonists, where the bias weights were based on the primary CPL scan. Relative functional sensitivities (pEC50 − pEC50 index) are plotted as survivor curves. Gray dashed line, theoretical curve (supplemental Equation S2, with α = 4.5, β = 10, γ = 2, N 0 = 208 anchorable peptides). The biased samples (black dashed and dotted lines) estimate the TCR degeneracy spectrum, whereas the motif-based samples (black solid lines) provide a lower boundary to the TCR degeneracy spectrum. DISCUSSION Despite the huge potential importance of TCR degeneracy to human health, there has never been a comprehensive attempt to quantify the number of peptides that can be recognized by a single TCR. To address this issue, we examined the extent of cross-reactivity exhibited by a single autoimmune T cell clone with pathogenic relevance in human type I diabetes (1E6). Our analysis demonstrates that the 1E6 TCR can recognize ∼500 peptides within a factor 2 of the optimal agonist (i.e. peptides that have a functional sensitivity that is at least 50% of the functional sensitivity of the optimal agonist). An estimated 60,000 peptides have a functional sensitivity within a factor 10 of the optimal agonist, and ∼1.3 × 106 peptides are within a factor 100 of the optimal agonist. These considerations are especially significant given that the functional sensitivity of 1E6 CD8+ T cells for the index peptide is at least 100-fold lower than the optimal agonist; this is illustrated by RQFGPDFPTI, the functional sensitivity of which is 100-fold better than the index. Almost 10 million peptides are within a factor 1,000 of the optimal agonist, but such weak agonists will not generally be physiologically significant unless presented at very high copy numbers. Taken together, these results indicate that the 1E6 TCR has >one million significant peptide agonists at concentrations with the potential to be physiologically relevant. When putting TCR degeneracy into perspective, it is important to realize that individual TCRs capable of recognizing 106 decamer peptides still only have a less than 1 in 107 chance of cross-reacting with any peptide chosen at random from the entire decamer peptide universe (∼1013). A high level of cross-reactivity is therefore amply compatible with the degree of specificity required for self/non-self determination. Furthermore, the number of decamer peptides that it is possible to make from the entire human proteome (excluding post-translational amino acid modifications) is only one millionth of the possible peptide universe at this length. Functional recognition of ∼106 decamer peptides by a single TCR translates into a frequency of cross-reactivity of 1:100,000 (assuming that 1% of peptides bind to MHC), which is likely to be the most accurate estimate of this parameter to date due to the comprehensive nature of our approach. The probability of cross-reactivity with any individual peptide is an important consideration in terms of viral escape, bystander activation and autoimmune side effects, and the results presented here fit well with theoretical considerations of T cell immunity (3). It should be noted, however, that the degeneracy curves are estimates based on samples that constitute a small fraction of the number of possible peptide ligands and, hence, should be considered as depicting an order-of-magnitude estimate that is supported by the lower bounds inferred from motif-based samples. The 1E6 CD8+ T cell clone was chosen for these studies to highlight the huge potential for T cell cross-reactivity as a possible cause of autoimmunity. In support of the generality of our findings, we have also observed high levels of degeneracy at some positions in a 9-amino acid residue nonautoimmune epitope (32). Furthermore, the recognition of longer, MHC class II-restricted peptides by CD4+ T cells with a “TCR footprint” of similar size could ensure that these cells incorporate the capacity to recognize tens, or possibly even hundreds, of millions of peptides at physiologically relevant surface densities (33). The reality of T cell cross-reactivity, as quantified here, has far reaching implications. It provides an explanation for how a limited pool of TCRs can provide the broad antigenic coverage that is required for effective immunity. In addition, the extent of TCR degeneracy suggests that almost any peptide antigen can be improved for any given cognate TCR, in the sense of there being at least one stronger agonist than the original index peptide, thereby providing scope for rational therapeutic interventions based on the directed manipulation of T cell immunity. Supplementary Material Supplemental Data
              Bookmark
              • Record: found
              • Abstract: found
              • Article: not found

              Overlap and effective size of the human CD8+ T cell receptor repertoire.

              Diversity in T lymphocyte antigen receptors is generated by somatic rearrangement of T cell receptor (TCR) genes and is concentrated within the third complementarity-determining region 3 (CDR3) of each chain of the TCR heterodimer. We sequenced the CDR3 regions from millions of rearranged TCR beta chain genes in naïve and memory CD8(+) T cells of seven adults. The CDR3 sequence repertoire realized in each individual is strongly biased toward specific V(beta)-J(beta) pair utilization, dominated by sequences containing few inserted nucleotides, and drawn from a defined subset comprising less than 0.1% of the estimated 5 x 10(11) possible sequences. Surprisingly, the overlap in the naïve CD8(+) CDR3 sequence repertoires of any two of the individuals is approximately 7000-fold larger than predicted and appears to be independent of the degree of human leukocyte antigen matching.
                Bookmark

                Author and article information

                Journal
                Journal of Clinical Investigation
                American Society for Clinical Investigation
                0021-9738
                1558-8238
                April 29 2019
                April 29 2019
                June 3 2019
                May 13 2019
                May 13 2019
                June 3 2019
                : 129
                : 6
                : 2446-2462
                Article
                10.1172/JCI124358
                6546456
                30920391
                c7f3f8e8-a643-4672-a83a-0a7fedc9addb
                © 2019
                History

                Comments

                Comment on this article