CN104673892B - 开发恒河猴t细胞免疫组库的引物组、高通量测序方法及应用 - Google Patents
开发恒河猴t细胞免疫组库的引物组、高通量测序方法及应用 Download PDFInfo
- Publication number
- CN104673892B CN104673892B CN201410843543.XA CN201410843543A CN104673892B CN 104673892 B CN104673892 B CN 104673892B CN 201410843543 A CN201410843543 A CN 201410843543A CN 104673892 B CN104673892 B CN 104673892B
- Authority
- CN
- China
- Prior art keywords
- rhesus macaque
- tcr
- seq
- cell immune
- primer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 210000001744 T-lymphocyte Anatomy 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000012165 high-throughput sequencing Methods 0.000 title abstract description 11
- 241000282560 Macaca mulatta Species 0.000 claims abstract description 67
- 239000012634 fragment Substances 0.000 claims description 34
- 239000002773 nucleotide Substances 0.000 claims description 27
- 125000003729 nucleotide group Chemical group 0.000 claims description 27
- 230000000692 anti-sense effect Effects 0.000 claims description 18
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 210000004027 cell Anatomy 0.000 claims description 11
- 238000002156 mixing Methods 0.000 claims description 9
- 238000007622 bioinformatic analysis Methods 0.000 claims description 6
- 210000002966 serum Anatomy 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 5
- 238000001962 electrophoresis Methods 0.000 claims description 4
- 238000010079 rubber tapping Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 238000004445 quantitative analysis Methods 0.000 claims description 2
- 238000013467 fragmentation Methods 0.000 claims 1
- 238000006062 fragmentation reaction Methods 0.000 claims 1
- 238000012163 sequencing technique Methods 0.000 abstract description 39
- 108090000623 proteins and genes Proteins 0.000 abstract description 38
- 238000004458 analytical method Methods 0.000 abstract description 28
- 241000282414 Homo sapiens Species 0.000 abstract description 24
- 238000007403 mPCR Methods 0.000 abstract description 16
- 201000010099 disease Diseases 0.000 abstract description 9
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 9
- 208000035473 Communicable disease Diseases 0.000 abstract description 6
- 238000011144 upstream manufacturing Methods 0.000 abstract description 6
- 210000005259 peripheral blood Anatomy 0.000 abstract description 5
- 239000011886 peripheral blood Substances 0.000 abstract description 5
- 230000001404 mediated effect Effects 0.000 abstract description 3
- 230000001363 autoimmune Effects 0.000 abstract description 2
- 230000002349 favourable effect Effects 0.000 abstract 2
- 230000001737 promoting effect Effects 0.000 abstract 1
- 150000001413 amino acids Chemical group 0.000 description 18
- 238000011160 research Methods 0.000 description 16
- 108020004414 DNA Proteins 0.000 description 14
- 238000009826 distribution Methods 0.000 description 13
- 239000013612 plasmid Substances 0.000 description 11
- 235000001014 amino acid Nutrition 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 9
- 108091008874 T cell receptors Proteins 0.000 description 8
- 102000016266 T-Cell Antigen Receptors Human genes 0.000 description 8
- 239000003814 drug Substances 0.000 description 8
- 238000003780 insertion Methods 0.000 description 8
- 230000037431 insertion Effects 0.000 description 8
- 241000894007 species Species 0.000 description 8
- 230000004087 circulation Effects 0.000 description 7
- 230000002068 genetic effect Effects 0.000 description 7
- 108020003175 receptors Proteins 0.000 description 7
- 102000005962 receptors Human genes 0.000 description 7
- 230000006798 recombination Effects 0.000 description 7
- 238000005215 recombination Methods 0.000 description 7
- 241000282553 Macaca Species 0.000 description 6
- 239000000203 mixture Substances 0.000 description 6
- 101150008942 J gene Proteins 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 230000000295 complement effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000002998 immunogenetic effect Effects 0.000 description 5
- 241000894006 Bacteria Species 0.000 description 4
- 101100510671 Rattus norvegicus Lnpep gene Proteins 0.000 description 4
- 108010087408 alpha-beta T-Cell Antigen Receptors Proteins 0.000 description 4
- 230000003321 amplification Effects 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 239000008280 blood Substances 0.000 description 4
- 239000003292 glue Substances 0.000 description 4
- 208000015181 infectious disease Diseases 0.000 description 4
- 210000004698 lymphocyte Anatomy 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 241000282693 Cercopithecidae Species 0.000 description 3
- 101150097493 D gene Proteins 0.000 description 3
- 241000252212 Danio rerio Species 0.000 description 3
- 108060003951 Immunoglobulin Proteins 0.000 description 3
- 241000288906 Primates Species 0.000 description 3
- 208000007660 Residual Neoplasm Diseases 0.000 description 3
- 102100026967 T cell receptor beta chain MC.7.G5 Human genes 0.000 description 3
- 108700042077 T-Cell Receptor beta Genes Proteins 0.000 description 3
- 101150117115 V gene Proteins 0.000 description 3
- 239000000427 antigen Substances 0.000 description 3
- 108091007433 antigens Proteins 0.000 description 3
- 102000036639 antigens Human genes 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 102000018358 immunoglobulin Human genes 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000000746 purification Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 201000008827 tuberculosis Diseases 0.000 description 3
- 238000002255 vaccination Methods 0.000 description 3
- 208000023275 Autoimmune disease Diseases 0.000 description 2
- 208000001490 Dengue Diseases 0.000 description 2
- 206010012310 Dengue fever Diseases 0.000 description 2
- 208000006664 Precursor Cell Lymphoblastic Leukemia-Lymphoma Diseases 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 2
- 241000703392 Tribec virus Species 0.000 description 2
- 102100022203 Tumor necrosis factor receptor superfamily member 25 Human genes 0.000 description 2
- 241000287433 Turdus Species 0.000 description 2
- 241000700605 Viruses Species 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 2
- 238000000246 agarose gel electrophoresis Methods 0.000 description 2
- 210000002203 alpha-beta t lymphocyte Anatomy 0.000 description 2
- 238000010171 animal model Methods 0.000 description 2
- 210000003719 b-lymphocyte Anatomy 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012350 deep sequencing Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 208000025729 dengue disease Diseases 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000005304 joining Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000035800 maturation Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 150000007523 nucleic acids Chemical class 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 235000015097 nutrients Nutrition 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000003752 polymerase chain reaction Methods 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000012552 review Methods 0.000 description 2
- 230000003612 virological effect Effects 0.000 description 2
- BSDCIRGNJKZPFV-GWOFURMSSA-N (2r,3s,4r,5r)-2-(hydroxymethyl)-5-(2,5,6-trichlorobenzimidazol-1-yl)oxolane-3,4-diol Chemical group O[C@@H]1[C@H](O)[C@@H](CO)O[C@H]1N1C2=CC(Cl)=C(Cl)C=C2N=C1Cl BSDCIRGNJKZPFV-GWOFURMSSA-N 0.000 description 1
- 108091032973 (ribonucleotides)n+m Proteins 0.000 description 1
- OPIFSICVWOWJMJ-AEOCFKNESA-N 5-bromo-4-chloro-3-indolyl beta-D-galactoside Chemical compound O[C@@H]1[C@@H](O)[C@@H](O)[C@@H](CO)O[C@H]1OC1=CNC2=CC=C(Br)C(Cl)=C12 OPIFSICVWOWJMJ-AEOCFKNESA-N 0.000 description 1
- 208000030507 AIDS Diseases 0.000 description 1
- 241000251468 Actinopterygii Species 0.000 description 1
- 208000024893 Acute lymphoblastic leukemia Diseases 0.000 description 1
- 208000014697 Acute lymphocytic leukaemia Diseases 0.000 description 1
- 229920001817 Agar Polymers 0.000 description 1
- 241000282709 Aotus trivirgatus Species 0.000 description 1
- 208000010839 B-cell chronic lymphocytic leukemia Diseases 0.000 description 1
- 101100263837 Bovine ephemeral fever virus (strain BB7721) beta gene Proteins 0.000 description 1
- 101150107439 CDR3 gene Proteins 0.000 description 1
- 241000288950 Callithrix jacchus Species 0.000 description 1
- 241000282552 Chlorocebus aethiops Species 0.000 description 1
- 108091026890 Coding region Proteins 0.000 description 1
- 108020004705 Codon Proteins 0.000 description 1
- 108090000695 Cytokines Proteins 0.000 description 1
- 102000004127 Cytokines Human genes 0.000 description 1
- 108010008286 DNA nucleotidylexotransferase Proteins 0.000 description 1
- 102100033215 DNA nucleotidylexotransferase Human genes 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 101100316840 Enterobacteria phage P4 Beta gene Proteins 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 241000997865 Histrio histrio Species 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 101000911390 Homo sapiens Coagulation factor VIII Proteins 0.000 description 1
- 206010061598 Immunodeficiency Diseases 0.000 description 1
- 208000029462 Immunodeficiency disease Diseases 0.000 description 1
- COLNVLDHVKWLRT-QMMMGPOBSA-N L-phenylalanine Chemical compound OC(=O)[C@@H](N)CC1=CC=CC=C1 COLNVLDHVKWLRT-QMMMGPOBSA-N 0.000 description 1
- 208000031422 Lymphocytic Chronic B-Cell Leukemia Diseases 0.000 description 1
- 241000282521 Macaca assamensis Species 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 241000699666 Mus <mouse, genus> Species 0.000 description 1
- 241000699670 Mus sp. Species 0.000 description 1
- 206010029260 Neuroblastoma Diseases 0.000 description 1
- 108091028043 Nucleic acid sequence Proteins 0.000 description 1
- 206010033128 Ovarian cancer Diseases 0.000 description 1
- 206010061535 Ovarian neoplasm Diseases 0.000 description 1
- 238000012408 PCR amplification Methods 0.000 description 1
- 241000282579 Pan Species 0.000 description 1
- 241000282577 Pan troglodytes Species 0.000 description 1
- 241001127637 Plantago Species 0.000 description 1
- 241000223801 Plasmodium knowlesi Species 0.000 description 1
- 108091008109 Pseudogenes Proteins 0.000 description 1
- 102000057361 Pseudogenes Human genes 0.000 description 1
- 208000025747 Rheumatic disease Diseases 0.000 description 1
- 241000288960 Saguinus oedipus Species 0.000 description 1
- 238000012300 Sequence Analysis Methods 0.000 description 1
- 241000713311 Simian immunodeficiency virus Species 0.000 description 1
- 230000005867 T cell response Effects 0.000 description 1
- 108700042075 T-Cell Receptor Genes Proteins 0.000 description 1
- 210000000662 T-lymphocyte subset Anatomy 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000008272 agar Substances 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000007720 allelic exclusion Effects 0.000 description 1
- 102000006707 alpha-beta T-Cell Antigen Receptors Human genes 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005784 autoimmunity Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000002619 cancer immunotherapy Methods 0.000 description 1
- 210000002421 cell wall Anatomy 0.000 description 1
- 230000007969 cellular immunity Effects 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 208000032852 chronic lymphocytic leukemia Diseases 0.000 description 1
- 238000010367 cloning Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 239000002299 complementary DNA Substances 0.000 description 1
- 230000009514 concussion Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 235000018417 cysteine Nutrition 0.000 description 1
- XUJNEKJLAYXESH-UHFFFAOYSA-N cysteine Natural products SCC(N)C(O)=O XUJNEKJLAYXESH-UHFFFAOYSA-N 0.000 description 1
- 230000003013 cytotoxicity Effects 0.000 description 1
- 231100000135 cytotoxicity Toxicity 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 239000000539 dimer Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000017188 evasion or tolerance of host immune response Effects 0.000 description 1
- 230000010429 evolutionary process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000684 flow cytometry Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 210000004475 gamma-delta t lymphocyte Anatomy 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 239000008236 heating water Substances 0.000 description 1
- 210000002443 helper t lymphocyte Anatomy 0.000 description 1
- 238000011134 hematopoietic stem cell transplantation Methods 0.000 description 1
- 239000000833 heterodimer Substances 0.000 description 1
- 102000057593 human F8 Human genes 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 210000000987 immune system Anatomy 0.000 description 1
- 201000007294 immune system cancer Diseases 0.000 description 1
- 208000026278 immune system disease Diseases 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000007813 immunodeficiency Effects 0.000 description 1
- 230000016784 immunoglobulin production Effects 0.000 description 1
- 230000036046 immunoreaction Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- 206010022000 influenza Diseases 0.000 description 1
- 230000003834 intracellular effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- BPHPUYQFMNQIOC-NXRLNHOXSA-N isopropyl beta-D-thiogalactopyranoside Chemical compound CC(C)S[C@@H]1O[C@H](CO)[C@H](O)[C@H](O)[C@H]1O BPHPUYQFMNQIOC-NXRLNHOXSA-N 0.000 description 1
- 238000012177 large-scale sequencing Methods 0.000 description 1
- 208000032839 leukemia Diseases 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 210000000207 lymphocyte subset Anatomy 0.000 description 1
- 208000003747 lymphoid leukemia Diseases 0.000 description 1
- 201000004792 malaria Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 108020004999 messenger RNA Proteins 0.000 description 1
- 238000010369 molecular cloning Methods 0.000 description 1
- 210000004940 nucleus Anatomy 0.000 description 1
- 230000008506 pathogenesis Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 210000003819 peripheral blood mononuclear cell Anatomy 0.000 description 1
- COLNVLDHVKWLRT-UHFFFAOYSA-N phenylalanine Natural products OC(=O)C(N)CC1=CC=CC=C1 COLNVLDHVKWLRT-UHFFFAOYSA-N 0.000 description 1
- 239000013600 plasmid vector Substances 0.000 description 1
- 210000004180 plasmocyte Anatomy 0.000 description 1
- 238000007747 plating Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000000770 proinflammatory effect Effects 0.000 description 1
- 238000012175 pyrosequencing Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 229940047431 recombinate Drugs 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 210000003289 regulatory T cell Anatomy 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010839 reverse transcription Methods 0.000 description 1
- 238000003757 reverse transcription PCR Methods 0.000 description 1
- 206010039073 rheumatoid arthritis Diseases 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000028327 secretion Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000002864 sequence alignment Methods 0.000 description 1
- 208000021160 simian immunodeficiency virus infection Diseases 0.000 description 1
- YEENEYXBHNNNGV-XEHWZWQGSA-M sodium;3-acetamido-5-[acetyl(methyl)amino]-2,4,6-triiodobenzoate;(2r,3r,4s,5s,6r)-2-[(2r,3s,4s,5r)-3,4-dihydroxy-2,5-bis(hydroxymethyl)oxolan-2-yl]oxy-6-(hydroxymethyl)oxane-3,4,5-triol Chemical compound [Na+].CC(=O)N(C)C1=C(I)C(NC(C)=O)=C(I)C(C([O-])=O)=C1I.O[C@H]1[C@H](O)[C@@H](CO)O[C@]1(CO)O[C@@H]1[C@H](O)[C@@H](O)[C@H](O)[C@@H](CO)O1 YEENEYXBHNNNGV-XEHWZWQGSA-M 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000002992 thymic effect Effects 0.000 description 1
- 210000001541 thymus gland Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 210000003171 tumor-infiltrating lymphocyte Anatomy 0.000 description 1
- 229960005486 vaccine Drugs 0.000 description 1
- 230000029812 viral genome replication Effects 0.000 description 1
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明提供了一种开发恒河猴T细胞免疫组库的引物组、高通量测序方法及应用。本发明通过针对恒河猴TCR的可变区V区(variable)设置了26条上游引物序列，针对TCR的连接区(joining)J区设置了13条下游引物序列，通过多重PCR扩增出外周血或者其他免疫组织中高达107左右数量级的TCR序列，制备高通量测序文库并测序，通过生物信息学对恒河猴TCR基因免疫组库进行了全面分析。本发明提供的技术方案对构建恒河猴分析模型意义重大，有利于促进恒河猴模型在研究自体免疫系统疾病和传染病中的应用，也有利于人类对T细胞介导的疾病有更深入和广泛的认识。
Description
技术领域
本发明属于分子生物学领域，特别是涉及一种开发恒河猴T细胞免疫组库的引物组、高通量测序方法及应用。
背景技术
T细胞受体(TCR)在细胞免疫中发挥着非常重要的作用，TCR属于异源二聚体分子，在循环系统中，超过95％的TCR的属于αβ类型。在人类外周血中，有超过2×107特异的TCRα和TCRβ配对。所述的免疫组库的复杂性是由淋巴细胞的成熟过程中，染色体上三类基因片段，包括V、D、J基因重排以及三种基因片段连接处的插入缺失生成的。该关键区域也被成为CDR3，是主要的抗原识别位点，也称为互补决定区3(CDR3)，它是TCR基因可变区中最多样和复杂的区域。
恒河猴(猕猴)这种动物模型被广泛用于研究许多重要的人类疾病，尤其是用于研究T淋巴细胞在发病机制中的作用。因此，该非人灵长类的TCR序列的表征是非常重要的。早在1992年，Gene Levinson et al.等，通过传统RT-PCR和克隆的方法，在隆恒河猴的外周血测到了23种TCRβ重排序列。后来，Emma E.M.Jaeg又在这类猕猴中发现了若干其他的TCRβ序列。通过序列比对，我们发现这些TCRβ在猕猴，黑猩猩和人类中的多样性和结构非常相似，具有很高的同源性。这些前期的工作，一共在恒河猴发现了超过23种Vβ基因，为后续建立更完整的TCRβ序列提供了至关重要的一个公共序列库。然而，由于传统克隆技术的低通量性质，在过去的二十年中恒河猴中发现的有关TCR序列的数目大约也就几百条，远远少于预估的107，这也成为了限制我们了解T细胞在免疫疾病中的演化过程的重要因素。
近年来，基于下一代测序这个强大的新技术在免疫系统中的应用，一个新的领域，称为免疫基因组库测序的技术诞生了。我们可以将循环血液中来自单个样品的数以百万计的TCR序列，在单次测序中全部一次性读取。这项技术，最早是在2009年，Weinstein等人，应用高通量测序，估计斑马鱼抗体库的大小，用这种方法，他们在每条鱼中检测到了1200至3700条抗体序列。之后Freeman等人，将该方法用于检测人T细胞检测，并且检测到了33664种不同的TCRβ克隆。随后由于测序方法的改进和通量的增加，Wang等在单个健康人中检测到了113290条独特的TCRβCDR3序列。运用统计模型，估计出的T细胞受体库的总多样性在106左右。运用这种免疫组库测序方法，科研人员检测到了抗体的不同过程的动态变化，并用于追踪淋巴细胞成熟过程以及在不同年龄组之间的差异。最近，科研人员用这项技术检测各种人类疾病。例如，检测各种白血病微小残留病变，他的灵敏度较传统的流式细胞分析法更灵敏，更可靠。这种方法也用于研究免疫系统癌症，传染性疾病如登革热，和自身免疫疾病特别是类风湿性关节炎。
尽管免疫组库测序技术在人类基础研究领域中已经有广泛应用，但这种技术在其他模式生物的使用非常有限，仅限于小鼠和斑马鱼。猕猴是一种重要的动物模型，尤其是对人类疾病的研究。超过70种传染病的研究都在恒河猴模型中进行，特别是由人类免疫缺陷病毒(HIV)引起的获得性免疫缺陷综合症 (AIDS)，结核病(TB)和疟疾。因此，建立对恒河猴TCRβ进行大规模测序的方案对构建该物种分析模型十分必要。
发明内容
鉴于此，本发明提供了一种开发恒河猴T细胞免疫组库的引物组、高通量测序方法及应用。
第一方面，本发明提供了一种开发恒河猴T细胞免疫组库的引物组，包括26条上游引物和13条下游引物，其中，所述26条上游引物为如SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列，所述13条上游引物为如SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列。
本发明通过针对恒河猴TCR的可变区V区(variable)设置至少26条上游引物序列，针对TCR的连接区(joining)J区设置至少13条下游引物序列，通过多重PCR扩增出目的链的PCR产物，获得高通量测序文库。
本发明采用至少13条下游引物序列与所述至少26条上游引物序列随机组合成配对引物，然后进行多重PCR反应，扩增TCRβ区。
本发明针对TCR所有的64类V和13类J基因家族进行了比对分析，26条上游引物序列针对TCR的CDR3区上游(即FR3区)，13条下游引物针对J基因，扩增TCR的CDR3区。
优选地，所述的上游引物为3’端比SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列多或少1～3个碱基的核苷酸序列，所述下游引物的为3’端比SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列多或少1～3个碱基的核苷酸序列；
所述上游引物或下游引物多出的1～3个碱基为与目的TCR互补的碱基。
本发明所述的与设计目的TCR互补的碱基是指：针对恒河猴TCR的可变区V区(variable)或针对TCR的连接区(joining)J区设计引物。
第二方面，本发明提供了一种检测恒河猴T细胞免疫组库的高通量测序方法，包括如下步骤：
提取恒河猴基因组DNA或总RNA；
采用开发恒河猴T细胞免疫组库的引物组，进行多重PCR反应；
将所得多重PCR产物进行高通量测序；
其中，所述开发恒河猴T细胞免疫组库的引物组包括26条上游引物和13条下游引物，其中，所述26条上游引物为如SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列，所述13条上游引物为如SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列。
优选地，所述多重PCR反应的体系中，26条上游引物等摩尔混合；13条下游引物等摩尔混合。
优选地，所述多重PCR反应的体系中，模板量为1～3ug/50ul体系。
优选地，所述多重PCR反应程序为：
优选地，所述多重PCR反应结束后，电泳，割胶回收片段长度为100-150bp的DNA片段。
优选地，所述的上游引物为3’端比SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列多或少1～3个碱基的核苷酸序列，所述下游引物的为3’端比SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列多或少1～3个碱基的核苷酸序列；
所述上游引物或下游引物多出的1～3个碱基为与目的TCR互补的碱基。
本发明提供的TCR高通量测序文库能获得的长度和数量丰富的TCR序列，有利于TCR序列的多态性程度分析及TCRCDR3区长度多态性分布分析。
优选地，所述T淋巴细胞受体(TCR)待测序列为采用RNA试剂盒提取人外周血单个核细胞的总RNA或者获得的基因组DNA。
总RNA样品需要先进行逆转录，再进行多重PCR。
第三方面，本发明提供了一种构建恒河猴T细胞免疫组库的方法，包括：
提取恒河猴基因组DNA或总RNA；
采用开发恒河猴T细胞免疫组库的引物组，进行多重PCR反应；
将所得多重PCR产物进行高通量测序；
生物信息学分析，对高通量测序结果进行定量分析；
其中，所述开发恒河猴T细胞免疫组库的引物组包括26条上游引物和13条下游引物，其中，所述26条上游引物为如SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列，所述13条上游引物为如SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列。
在高通量测序平台的基础上，通过对恒河猴TCR基因免疫组库进行全面的生物信息学分析，获得了TCR在VDJ重组时的基因偏好性(usage patterns)，基因组合、连接多样性信息，以及CDR3序列中氨基酸的偏好性(usage patterns)、CDR3氨基酸序列的长度多样性以及连接处N端碱基的特性。正是这些因素形成数量庞大且种类多样的TCR受体库。
优选地，所述生物信息学分析采用了发明人自编的在线生物学分析软件“iRAP”，网址为： http://sustc-genome.org.cn/irap。
优选地，所述多重PCR反应的体系中，26条上游引物等摩尔混合；13条下游引物等摩尔混合。
优选地，所述的上游引物为3’端比SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列多或少1～3个碱基的核苷酸序列，所述下游引物的为3’端比SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列多或少1～3个碱基的核苷酸序列；
所述上游引物或下游引物多出的1～3个碱基为与目的TCR互补的碱基。
优选地，所述多重PCR反应的体系中，模板量为1～3ug/50ul体系。
优选地，所述多重PCR反应程序为：
优选地，所述多重PCR反应结束后，电泳，割胶回收片段长度为100-150bp的DNA片段。
第四方面，本发明提供了一种如第一方面所述的开发恒河猴T细胞免疫组库的引物组在检测恒河猴T细胞受体、对恒河猴T细胞受体进行高通量测序、构建恒河猴T细胞免疫组库、以及制备恒河猴T细胞免疫组库试剂盒中的应用。
本发明提供的技术方案对构建恒河猴分析模型意义重大，有利于促进恒河猴模型在研究自体免疫系统疾病和传染病中的应用，也有利于人类对T细胞介导的疾病有更深入和广泛的认识。
采用本发明的技术方案，达到如下有益效果：
1)获得了恒河猴1.26(百万)条TCR序列；
2)获得了恒河猴643,570条特异性TCR克隆；
2)获得了恒河猴270,557条特异的TCR CDR3序列，这大大扩展了对该物种免疫TCR的认识。而目前IMGT和NCBI数据库所有恒河猴TCR相关的序列只有几百条。
附图说明
图1-a为TCRβVDJ重组、TCR文库构建以及生物信息流程示意图；
图1-b为高通量测序序列的统计分析结果；
图1-c为高通量测序TCRβ序列的频率直方图；
图2-a TCRβV基因使用频率分布
图2-b TCRβJ基因使用频率分布
图2-c TCRβD基因使用频率分布
图2-d为整个VDJ基因组的3D分布图；
图2-e为不同VDJ重组频率和的数量对应分布情况；
图2-f为测序深度和VDJ重组片段饱和度分析结果；
图3-a TCRβ的CDR3氨基酸长度分布和不同J基因使用分布；
图3-b为12aa长度的TCRβ的CDR3氨基酸组成成份分析；
图3-c为TCRβ的CDR3序列中，不同氨基酸采用的密码子情况；
图3-d为不同TCRβ的CDR3序列的出现频率图；
图3-e为CDR3序列组大小和多样性统计分析；
图4-a为VDJ重组片段区域的13个分区；
图4-b为VDJ重组片段区域中，插入或缺失的片段长度分布情况；
图5-a为PCR偏差的分析；
图5-b为测序误差的评估。
具体实施方式
图1-a为TCRβVDJ重组、TCR文库构建以及生物信息流程示意图；图1-a显示TCRβ的多样性来源 于VDJ重组以及连接多样性，我们通过设计多重PCR引物组扩增TCR基因，并进行测序和分析；根据图1，本发明提供了如下实施例：
材料及试剂说明：
恒河猴：10年大的母猴，来源于广州生物医药与健康研究院(GIBH)，饲养按行业标准。非特殊说明，本发明实施例采用的试剂均为市售商品，本发明实施例采用的数据库均为公开的在线数据库。
实施例1
T淋巴细胞受体(TCR)基因组DNA样品的制备，包括如下步骤：
(1)收集的新鲜外周血样本各10毫升(ml)，按LymphoPrep试剂盒(Axis-shield,Cat.No.AS1114544UK)说明书操作，获得相对较纯的PBMC；
(2)采用PureLink Genomic DNA Mini Kit(Life Technology,Cat.No:K1820-00)试剂盒提取步骤(1)所得细胞的基因组DNA，并用Nanodrop2000(Thermo)测定DNA的浓度及纯度，然后保存基因组DNA。
实施例2
采用开发恒河猴T细胞免疫组库的引物组进行多重PCR，构建TCRβ高通量测序文库，进行高通量测序，包括如下步骤：
(1)以实施例1所得基因组DNA为扩增模板，取TCR引物，再采用QIAGEN公司Multiplex PCR试剂盒(货号：206143)，按试剂盒说明书配置多重PCR体系，其中，每个前向引物等摩尔混合，引物总浓度是10微摩尔，每个反向引物等摩尔混合，引物总浓度是10微摩尔，模板量可以调整，本实施例中采用3ug，1ug也可以做出来，但是多样性会下降。
该体系中的上游引物为如SEQ ID NO:1～SEQ ID NO:26所示引物的等摩尔混合物；下游引物为如SEQ ID NO:27～SEQ ID NO:39所示引物的等摩尔混合物。
(2)多重PCR
为方便送测，若无特别说明，本发明实施例进行多重PCR时，分别在上游引物和下游引物加上测序接头，具体为：在如SEQ ID NO:1～SEQ ID NO:26所示的引物的5’端分别接上illumina测序公司的上游引物接头序列，在下游引物的5’端分别接上illumina测序公司的下游引物接头序列，具体步骤参照illumina高通量测序文库构建说明书；
再按下述多重PCR的条件设置PCR仪器程序，进行多重PCR：
72℃7min，PCR结束后，4℃保存PCR产物并电泳检测，挑选片段长度约为120bp的TCRβ片段，割胶回收，得到纯化后的TCRβ片段，胶回收步骤采用QIAGEN公司QIAquick胶纯化试剂盒，按常规实验室操作进行；Nanodrop 2000测试DNA浓度，并送公司进行高通量测序(采用Illumina hiseq2000测序，2*100pair-end)。
本发明实施例设计26条上游引物，13条下游引物对TCRβ免疫组重链进行分析：
1)设计引物：针对TCR所有的64类V和13类J基因进行了比对分析，采用Oligo 7.0和MFEprimer-2.0对引物二聚体以及茎环错配进行分析，在TCR的CDR3区上游(即FR3区)设置了上游引物，针对J基因下游设计反向引物，扩增TCRβ基因。
本实施例提供的引物组覆盖了大部分TCRβVDJ重组片段。由于很小的序列变化将导致引物扩增效果显著降低，本发明选取了扩增效果最佳的引物组。
2)多重PCR产物大小在120bp左右。
3)采用Illumina Hiseq 2*100pair-end测序。
4)对高通量测序结果的生物信息学分析采用本发明人自研在线软件ImmuneRepertoire Analysis Pipeline(iRAP，http://www.sustc-genome.org.cn/irap/)
生物信息学分析结果如下：
一、TCR基因序列的分析
采用本发明的引物组以及多重PCR构库后，高通量测序得到大约1.69(百万)条序列。经过和TCRβ比对分析，得到图1-b所示结果，图1-b为高通量测序序列的统计分析结果，由图1-b可知，1.69(百万)条序列中，1.26(百万)条为TCRβ序列，其中，1.07(百万)条为可正常编码氨基酸序列的CDR3序列，另有0.19(百万)条在CDR3不能正常编码序列，T细胞中含有的out-of-frame序列，如果不被nonsense-mediated decay(NMD)途径沉默或破坏，将会引起严重疾病。
CDR3区域是TCRβ基因中多样性最丰富的区域，是抗原结合位点的关键，本次分析，我们发现了0.23(百万)条CDR3氨基酸序列，这可以和人类CDR3多样性相匹敌。图1-c为高通量测序TCRβ序列的频率直方图，结果显示，25.2％的TCR序列只有1个拷贝，所有TCR CDR3的寡克隆指数(oligoclonality index，OI)为0.63，表明CDR3序列拷贝数的分布非常不一致，变化幅度很大。
二、VDJ基因出现频率的规律
根据IMGT数据库可知，恒河猴具有64种V,2种D和13种J基因片段。本实施例参考现有数据库，通过比对分析每一条高通量测序所得TCRβ序列，获得57种V,2种D和13种J基因片段，由此可知，本发明提供的多重PCR引物组覆盖了大部分VDJ基因片段。
根据分析，V、D、J基因片段的丰度非常不一致，有些片段丰度显著高于其他片段，结果如图2-a～c所示。图2-a TCRβV基因使用频率分布，图2-b TCRβJ基因使用频率分布，图2-c TCRβD基因使用频率分布。横坐标为不同的V,D,J基因家族，不同的V基因家族，纵坐标为每类V基因家族的序列占 总序列数(Total reads)的百分比，可以看出：丰度最高的V基因片段为V6-1,V6-3和V6-4，丰度排前8的V基因片段共占所有V基因的50％。同样，J基因中，J2-7占所有J基因的27％。此外，我们发现，D1的丰度高于D2。图2-a～c的现象可能是受重组过程中染色体空间位置影响，即空间位置越近的越容易重组在一起。
由上可知，恒河猴具有1664中可能的VDJ重组方式(64V×2D×13J＝1,664VDJ)，本发明实施例获得了1094中VDJ重组，覆盖率达66％，结果如图2-d所示。
图2-d为整个VDJ基因组的3D分布图，3个坐标轴代表了所有可能的VDJ值，每一个(V,D,J)坐标代表一个特异的VDJ重组方式，图中小球的体积代表了该点对应的VDJ重组方式的数量。红、绿、蓝(竖轴从上至下)分别代表D2、D1和未确定的D，由于TCR D是较复杂的的区域，可能存在不稳定的缺失或转化，导致在VDJ重组过程中的D基因缺失，与此相符的，本发明分析发现，在对D基因的测序结果进行分析后，大概有16.8％的D基因测序结果比对不到确定的、对应的D基因，这些基因归类为未确定的D。
由于VDJ3种基因片段丰度的不同，我们预测VDJ重组种类也有所差异，与此相符的是，如图2-d中，在所有1094种重复方式中，不同VDJ重组方式的丰度在0.01％-1.85％之间变化，丰度最高的VDJ重组(V10-1/D2/J2-7)具有12057条测序所得序列，然而，最小丰度的VDJ重组至对应1条测序所得序列。
进一步地，如图2-e所示，图2-e为不同VDJ重组频率和数量对应分布情况。发明人随机在测序所得基因文库中选择一定比例的基因片段，统计了不同VDJ重组种类的数量。由图2-e可知，VDJ重组种类丰度变化曲线具有长尾巴分布规律。
图2-f为为测序深度和VDJ重组片段饱和度分析结果。
通过多重采样分析结果显示目前的测序深度对于TCR的多样性已经接近饱和(Fig.2f).
三、CDR3的特点
CDR3氨基酸序列中，1个氨基酸(aa)的变化可能导致受体构象的改变，因此，CDR3氨基酸序列长度的变化可以反映出CDR3基因连接区的多样性。本发明实施例中，恒河猴CDR3的定义采用人类中一样的定义，即，从TRBV最后一个半胱氨酸到TRBJ片段的FGXG结构域中的苯丙氨酸之间的序列。可以发现，CDR3的aa长度为10-16(核苷酸序列30-48bp)，其中，84％的CDR3序列长度为11-13aa之间，如图3-a所示。图3-a为TCRβ的CDR3氨基酸长度分布和不同J基因使用分布；横坐标为氨基酸序列的长度，纵坐标为某种长度CDR3序列的数量，其中，柱子中的不同颜色代表了不同的J基因，由图3-a可知，所有的TCRβ的CDR3氨基酸序列均采用了J2-1基因；此外，12个aa长度的TCRβ的CDR3序列数量最庞大；图3-a还显示了部分具有13个氨基酸(aa)的序列的数量。
我们以12aa长度的TCRβ的CDR3序列为目标，采用Crooks et al.中的WebLogo法，研究CDR3氨基酸组成组成成份，如图3-b所示。图3-b为12aa长度的TCRβ的CDR3序列的组成情况。
图3-c反应了TCRβ的CDR3区域中氨基酸的种类，以及不同氨基酸采用的密码子情况；
图3-d为不同TCRβCDR3序列的丰度图；由图3-d可知，多数CDR3序列少于10个拷贝。
图3-e为CDR3序列组大小和多样性统计分析。
该统计是用于分析我们的测序深度是否足够用于评估目前的TCR多样性，结果显示我们的测序深度已经接近饱和，足够用于评估。
四、恒河猴的总TCRβCDR3多样性分析
CDR3多样性是正常免疫反应的基础。现有研究表明人类TCR□中特异的CDR3有340,000种。为了确定我们的测序是否覆盖恒河猴主要的特异性CDR3种类，我们参照FisherRA et al.中提到的采样-重采样统计分析方法对样本的大小和多样性进行了分析。结果如图3-e所示，相比由该方法预测的262,412种特异性CDR3序列，我们实际获得的特异性CDR3序列为235,573中，覆盖了总特异性CDR3序列的90％。
五、恒河猴连接区的多样性分析
恒河猴TCR多样性不仅来源于VDJ重组，很多TCR还来源于CDR3区域V-D和D-J之间连接区核苷酸的插入或缺失，此类插入或缺失不依赖于模板，而是通过末端脱氧核苷酸转移酶实现，在本发明中，我们参照IMGT数据库的定义，将恒河猴的TCRβ的连接区划分为13个区域，如图4-a所示。图4-a为VDJ重组片段区域的13个分区；我们对所有270,557种CDR3的核苷酸序列的每个区域均进行研究，分析了不同插入或缺失序列的长度分布情况，如图4-b所示。图4-b为VDJ重组片段区域的13个分区中，插入或缺失的片段长度分布情况，x轴表示13种片段，y轴表示VDJ重组过程中，插入或缺失的片段长度，z轴表示某种特定片段长度所占比例。z轴正方形为插入，负方向为缺失。N1和N2插入片段的平均长度分别为4bp和3bp。
注：
本发明测序结果分析方法：Illumina Hiseq 2000(paired-end 2*100)测序得到2,645,432条序列，我们用FLASH软件拼接有重叠的paired-end序列。拼接后的结果采用IMGT-High V-QUEST比对V、D和J基因。同源性低的从测序所得文库中去除。CDR3区域的起点和重点、编码框以及丰度的定义来源于IMGT。
我们用采样-重采样的方式研究整个CDR3免疫组库的大小，第一次取样中，占总CDR3免疫组库1％ 的序列中，获得953条特异性的CDR3序列，第二次取样中，占总CDR3免疫组库1％的序列中，获得885条特异性的CDR3序列，第三次取样中，占总CDR3免疫组库1％的序列中，获得9831条特异性的CDR3序列，我们总结出对数衰减值(R2＝0.97)。所有样品计算后，将所有特异性CDR3序列求和，可预测出整个CDR3免疫组库的大小。
我们用寡克隆指数(oligoclonality index，OI,见式1)来评估CDR3序列丰度的分布，如果是单克隆则OI值为1，如果不止1个克隆，OI值记为0，式1的定义基于Gini相关因素：
其中，S代表所有的TCRβCDR3克隆
Si表示特定TCRβCDR3克隆的丰度
测序所得的TCRβCDR3克隆
表示TCRβCDR3克隆的相对丰度
表示递减顺序下，TCRβCDR3克隆的累加相对丰度。
实施例3
矫正多重PCR引物扩增误差的方法，包括：
本发明的目的之一是分析TCR序列的丰度情况，然而，PCR扩增过程中，由于PCR引物对不同基因片段的扩增偏好，会影响TCR序列丰度的真实情况。
为了评估这类系统误差，我们采用了Robins et al中的方法，我们按照实施例1的方法提取1只恒河猴的基因组，均分为6管，3管采用25循环的多重PCR(参照实施例2)，3管采用30个循环，每管扩增产物均由Illumina Miseq测序，在25循环的实验组中，共发现10,345条特异的TCR CDR3序列，15.3％were also found in the 30-cycle PCR experiment.
我们将25循环组中发现的19658条测序结果和30循环组中发现的32647条测序结果进行分析，发现了如图5-a所示的线性关系。图5-a为PCR偏差的分析。对比25循环组和30循环组，假设误差全部来源于PCR误差，则平均每个PCR循环引入1.033，因此，整个30循环组引入的误差为1.03330＝2.7。
矫正高通量测序误差的方法，包括：
(1)制备多重PCR模板
(1a)构建质粒文库
参照实施例1和实施例2步骤(1)-(2)，进行多重PCR；
采用Invitrogen公司TOPO克隆试剂盒将实施例2步骤(2)胶回收纯化后的PCR片段与TA克隆的质粒载体连接；然后取2-3μl连接产物，加入至100μl JM109感受态细胞中，轻轻混匀，冰中静置30分钟后，再在42℃水浴加热60秒钟，然后迅速置于冰上1-2分钟；往转化产物中加入890μlLB培养基(amp-)，37℃、175rpm振荡培养60分钟，然后将培养液以4000rmp离心1min，弃上清若干，剩下的菌液混匀后取3/4的量在含有X-Gal(10ul)、IPTG(100ng20ul)、Amp的L-琼脂平板培养基上涂布，待菌液完全被培养基吸收后，37度倒置培养12-16h，形成单菌落，计数白色、蓝色菌落；
挑取每个平板上单个白色菌落，加到600ul LB培养液(Amp+)，恒温孵育箱225rpm37度震荡培养2-5h，然后取2.5ul菌液做PCR(25ul体系)；
将得到的PCR产物进行琼脂糖凝胶电泳，确认载体中插入片段的长度大小，测序后选取阳性克隆组成TCR质粒文库；
(1b)挑选26种质粒混合制备多重PCR模板
在(1a)所得TCR质粒文库中挑选出26种质粒，挑选标准为：根据BCR的C区、V、D和J区的特征，所得克隆与IMGT数据库中的序列进行生物信息学分析并分类，26种质粒为26种V区beta基因质粒，序列已知。
将该26种TCR质粒等摩尔混合后，作为多重PCR模板，每种质粒的添加数量是已知的，即每条DNA的数量是已知的；
(2)设置多重PCR引物
多重PCR引物为实施例2步骤1中多重PCR引物；
(3)多重PCR
取步骤(1)制备的多重PCR模板与多重PCR引物，参照实施例2步骤2的方法进行多重PCR，不同点在于设置3个实验组，每组平行3管，PCR循环次数分别设置为20、25、30，分别得到的PCR产物，将得到的PCR产物进行2％琼脂糖凝胶电泳，确认载体中插入片段的长度大小100-150bp，胶回收(Qiagen,Cat.No.28704,Valencia,CA)；
(4)测序分析
将所有步骤3的PCR产物送测(Illumina Miseq 2*150)，测序结果和原始质粒比对，计算误差(方法参照Robins et el.)。结果如图5-b所示。图5-b为测序误差的评价结果。
由图5-b可知，95.1％的测序结果和原始质粒完全一致，4.9％的测序结果在CDR3区域有至少1个bp的错误，经过计算，得到CDR3的测序误差为0.0186％/bp。
注：
本发明人采用的引物组：
本发明参考文献：
1Freeman,J.D.,Warren,R.L.,Webb,J.R.,Nelson,B.H.&Holt,R.A.Profilingthe T-cell receptor beta-chain repertoire by massively parallelsequencing.Genome research 19,1817-1824,doi:10.1101/gr.092924.109(2009).
2Arstila,T.P.et al.A direct estimate of the human alphabeta T cellreceptor diversity.Science 286,958-961(1999).
3Levinson,G.,Hughes,A.L.&Letvin,N.L.Sequence and diversity of rhesusmonkey T-cell receptor beta chain genes.Immunogenetics 35,75-88(1992).
4Jaeger,E.E.,Bontrop,R.E.&Lanchbury,J.S.Structure,diversity,andevolution of the T-cell receptor VB gene repertoire inprimates.Immunogenetics 40,184-191(1994).
5Hali,M.A.,Jaeger,E.E.,Bontrop,R.E.&Lanchbury,J.S.Characterization ofa novel human T cell receptor beta chain variable region family bytransspecies DNA hybridization.European journal of immunology 24,641-645,doi:10.1002/eji.1830240323(1994).
6Allen,T.M.,Lanchbury,J.S.,Hughes,A.L.&Watkins,D.I.The T-cellreceptor beta chain-encoding gene repertoire of a New World primate species,the cotton-top tamarin.Immunogenetics 45,151-160(1996).
7Uccelli,A.et al.Characterization of the TCRB chain repertoire in theNew World monkey Callithrix jacchus.Journal of immunology 158,1201-1207(1997).
8Currier,J.R.,Stevenson,K.S.&Robinson,M.A.Molecular cloning andsequencing of cDNA encoding eight new rhesus macaque TCRBV genefamilies.Immunogenetics 51,387-391(2000).
9Meyer-Olson,D.et al.Analysis of the TCR beta variable generepertoire in chimpanzees:identification of functional homologs to humanpseudogenes.Journal of immunology 170,4161-4169(2003).
10Moncada,C.A.et al.The T-cell receptor in primates:identifying andsequencing new owl monkey TRBV gene sub-groups.Immunogenetics 57,42-52,doi:10.1007/s00251-004-0758-y(2005).
11Weinstein,J.A.,Jiang,N.,White,R.A.,3rd,Fisher,D.S.&Quake,S.R.High-throughput sequencing of the zebrafish antibody repertoire.Science 324,807-810,doi:10.1126/science.1170020(2009).
12Wang,C.et al.High throughput sequencing reveals a complex patternof dynamic interrelationships among human T cell subsets.Proceedings of theNational Academy of Sciences of the United States of America 107,1518-1523,doi:10.1073/pnas.0913939107(2010).
13Robins,H.S.et al.Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells.Blood 114,4099-4107,doi:10.1182/blood-2009-04-217604(2009).
14Jiang,N.et al.Lineage structure of the human antibody repertoire inresponse to influenza vaccination. Science translational medicine 5,171ra119,doi:10.1126/scitranslmed.3004794(2013).
15Wu,Y.C.,Kipling,D.&Dunn-Walters,D.K.Age-Related Changes in HumanPeripheral Blood IGH Repertoire Following Vaccination.Frontiers in immunology3,193,doi:10.3389/fimmu.2012.00193(2012).
16Ademokun,A.et al.Vaccination-induced changes in human B-cellrepertoire and pneumococcal IgM and IgA antibody at different ages.Aging cell10,922-930, doi:10.1111/j.1474-9726.2011.00732.x(2011).
17Sherwood，A.M.et al.Deep sequencing of the human TCRgamma andTCRbeta repertoires suggests that TCRbeta rearranges after alphabeta andgammadelta T cell commitment.Science translational medicine 3,90ra61,doi:10.1126/scitranslmed.3002536(2011).
18Robins,H.Immunosequencing: applications of immune repertoire deepsequencing.Current opinion in immunology 25,646-652(2013).
19Wu,D.et al.High-throughput sequencing detects minimal residualdisease in acute T lymphoblastic leukemia.Science translational medicine 4,134ra163,doi:10.1126/scitranslmed.3003656(2012).
20Logan,A.C.et al.High-throughput VDJ sequencing for quantificationof minimal residual disease in chronic lymphocytic leukemia and immunereconstitution assessment.Proceedings of the National Academy of Sciences ofthe United States of America 108,21194-21199,doi:10.1073/pnas.1118357109(2011).
21Boyd,S.D.et al.Measurement and clinical monitoring of humanlymphocyte clonality by massively parallel VDJ pyrosequencing.Sciencetranslational medicine 1,12ra23(2009).
22van Heijst,J.W.et al.Quantitative assessment of T cell repertoirerecovery after hematopoietic stem cell transplantation.Nature medicine 19,372-377,doi:10.1038/nm.3100(2013).
23Grupp,S.A.et al.Adoptive transfer of autologous T cells improves T-cell repertoire diversity and long-term B-cell function in pediatric patientswith neuroblastoma.Clinical cancer research : an official journal of theAmerican Association for Cancer Research 18,6732-6741,doi:10.1158/1078-0432.CCR-12-1432(2012).
24Brentjens,R.J.et al.CD19-targeted T cells rapidly induce molecularremissions in adults with chemotherapy-refractory acute lymphoblasticleukemia.Science translational medicine 5,177ra138,doi:10.1126/scitranslmed.3005930(2013).
25Turtle,C.J.&Riddell,S.R.Genetically retargeting CD8+ lymphocytesubsets for cancer immunotherapy.Curr Opin Immunol 23,299-305,doi:10.1016/j.coi.2010.12.012(2011).
26Emerson,R.O.et al.High-throughput sequencing of T cell receptorsreveals a homogeneous repertoire of tumor-infiltrating lymphocytes in ovariancancer.The Journal of pathology,doi:10.1002/path.4260(2013).
27Maecker,H.T.et al.New tools for classification and monitoring ofautoimmune diseases.Nature reviews.Rheumatology 8,317-328,doi:10.1038/nrrheum.2012.66(2012).
28Parameswaran,P.et al.Convergent antibody signatures in humandengue.Cell hostµbe 13,691-700,doi:10.1016/j.chom.2013.05.008(2013).
29Klarenbeek,P.L.et al.Inflamed target tissue provides a specificniche for highly expanded T-cell clones in early human autoimmunedisease.Annals of the rheumatic diseases 71,1088-1093,doi:10.1136/annrheumdis-2011-200612(2012).
30DeKosky,B.J.et al.High-throughput sequencing of the paired humanimmunoglobulin heavy and light chain repertoire.Nature biotechnology 31,166-169,doi:10.1038/nbt.2492(2013).
31Reddy,S.T.et al.Monoclonal antibodies isolated without screening byanalyzing the variable-gene repertoire of plasma cells.Nature biotechnology28,965-969,doi:10.1038/nbt.1673(2010).
32Genolet,R.，Stevenson,B.J.,Farinelli,L.,Osteras,M.&Luescher,I.F.Highly diverse TCRalpha chain repertoire of pre-immune CD8(+) T cellsreveals new insights in gene recombination.The EMBO journal31,4247-4248，doi:10.1038/emboj.2012.277(2012).
33Gardner,M.B.&Luciw,P.A.Macaque models of human infectiousdisease.ILAR journal/National Research Council,Institute of Laboratory AnimalResources 49,220-255(2008).
34Price,D.A.et al.T cell receptor recognition motifs govern immuneescape patterns in acute SIV infection.Immunity 21,793-803,doi:10.1016/j.immuni.2004.10.010(2004).
35Kwofie,T.&Miura,T.Increased Virus Replication and Cytotoxicity ofNon-pathogenic Simian Human Immuno Deficiency Viruses-NM-3rN After SerialPassage in a Monkey-Derived Cell Line.Annals of medical and health sciencesresearch 3,55-61,doi:10.4103/2141-9248.109490(2013).
36Marcondes,M.C.et al.Simian immunodeficiency virus-induced CD4+ Tcell deficits in cytokine secretion profile are dependent on monkeyorigin.Viral immunology 19,679-689,doi:10.1089/vim.2006.19.679(2006).
37Murayama,Y.,Mukai,R.,Inoue-Murayama,M.&Yoshikawa,Y.An African greenmonkey lacking peripheral CD4 lymphocytes that retains helper T cell activityand coexists with SIVagm.Clinical and experimental immunology 117,504-512(1999).
38Sircar,P.,Furr,K.L.,Dorosh,L.A.&Letvin,N.L.Clonal repertoires ofvirus-specific CD8+T lymphocytes are shared in mucosal and systemiccompartments during chronic simian immunodeficiency virus infection in rhesusmonkeys.Journal of immunology 185,2191-2199,doi:10.4049/jimmunol.1001340(2010).
39Ribi,E.et al.Efficacy of mycobacterial cell walls as a vaccineagainst airborne tuberculosis in the Rheusus monkey.The Journal of infectiousdiseases 123,527-538(1971).
40Schmidt-Ullrich,R.,Wallach,D.F.&Lightholder,J.Two Plasmodiumknowlesi-specific antigens on the surface of schizont-infected Rhesus monkeyerythrocytes induce antibody production in immune hosts.The Journal ofexperimental medicine 150,86-99(1979).
41Mahowald,G.K.,Mahowald,M.A.,Moon,C.,Khor，B.&Sleckman,B.P.Out-of-frame T cell receptor beta transcripts are eliminated by multiple pathways invivo.PloS one 6,e21627,doi:10.1371/journal.pone.0021627(2011).
42Brogna,S.&Wen,J.Nonsense-mediated mRNA decay(NMD) mechanisms.Naturestructural&molecular biology 16,107-113,doi:10.1038/nsmb.1550(2009).
43Gillet,N.A.et al.The host genomic environment of the provirusdetermines the abundance of HTLV-1-infected T-cell clones.Blood 117,3113-3122,doi:10.1182/blood-2010-10-312926(2011).
44Giudicelli,V.,Chaume,D.&Lefranc,M.P.IMGT/GENE-DB: a comprehensivedatabase for human and mouse immunoglobulin and T cell receptor genes.NucleicAcids Res 33,D256-261,doi:10.1093/nar/gki010(2005).
45Brochet,X.,Lefranc,M.P.&Giudicelli,V.IMGT/V-QUEST: the highlycustomized and integrated system for IG and TR standardized V-J and V-D-Jsequence analysis.Nucleic acids research 36,W503-508,doi:10.1093/nar/gkn316(2008).
46Giudicelli,V.,Brochet,X.&Lefranc,M.P.IMGT/V-QUEST: IMGTstandardized analysis of the immunoglobulin(IG) and T cell receptor(TR)nucleotide sequences.Cold Spring Harbor protocols 2011,695-715,doi:10.1101/pdb.prot5633(2011).
47Genolet,R.,Stevenson,B.J.,Farinelli,L.,Osteras,M.&Luescher,I.F.Highly diverse TCRalpha chain repertoire of pre-immune CD8(+)T cellsreveals new insights in gene recombination.The EMBO journal 31,1666-1678,doi:10.1038/emboj.2012.48(2012).
48Kondilis-Mangum,H.D.,Shih,H.Y.,Mahowald,G.,Sleckman,B.P.&Krangel,M.S.Regulation of TCRbeta allelic exclusion by gene segment proximity andaccessibility.Journal of immunology 187,6374-6381,doi:10.4049/jimmunol.1102611(2011).
49Hsieh,C.S.,Lee,H.M.&Lio,C.W.Selection of regulatory T cells in thethymus.Nature reviews.Immunology 12,157-167,doi:10.1038/nri3155(2012).
50Naramura,M.,Kole,H.K.,Hu,R.J.&Gu,H.Altered thymic positiveselection and intracellular signals in Cbl-deficient mice.Proceedings of theNational Academy of Sciences of the United States of America 95,15547-15552(1998).
51Clute,S.C.et al.Broad cross-reactive TCR repertoires recognizingdissimilar Epstein-Barr and influenzaA virus epitopes.Journal of immunology185,6753-6764,doi:10.4049/jimmunol.1000812(2010).
52Degauque,N.et al.The blood of healthy individuals exhibits CD8 Tcells with a highly altered TCR Vb repertoire but with an unmodifiedphenotype.PloS one 6,e21240,doi:10.1371/journal.pone.0021240(2011).
53Cabaniols,J.P.,Fazilleau,N.,Casrouge,A.,Kourilsky,P.&Kanellopoulos,J.M.Most alpha/beta T cell receptor diversity is due to terminaldeoxynucleotidyl transferase.The Journal of experimental medicine194,1385-1390(2001).
54Crooks,G.E.,Hon，G.,Chandonia,J.M.&Brenner,S.E.WebLogo: a sequencelogo generator.Genome research 14,1188-1190,doi:10.1101/gr.849004(2004).
55Efron B,T.R.Estimating the number of un- seen species: How manywords did Shakespeare know？Biometrika 63,435-447(1976).
56Arif,S.et al.Autoreactive T cell responses show proinflammatorypolarization in diabetes but a regulatory phenotype in health.The Journal ofclinical investigation 113,451-463,doi:10.1172/JCI19585(2004).
57Arnaout，R.et al.High-resolution description of antibody heavy-chainrepertoires in humans.PloS one 6,e22365(2011).
58Gillet,N.A.et al.The host genomic environment of the provirusdetermines the abundance of HTLV-1–infected T-cell clones.Blood 117,3113-3122(2011).
Claims (8)
1.一种开发恒河猴T细胞免疫组库的引物组，其特征在于，包括26条上游引物和13条下游引物，其中，所述26条上游引物为如SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列，所述13条下游引物为如SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列。
2.一种检测恒河猴T细胞免疫组库的高通量测序方法，其特征在于，包括如下步骤：
提取恒河猴基因组DNA或总RNA；
采用开发恒河猴T细胞免疫组库的引物组，进行多重PCR反应；
将所得多重PCR产物进行高通量测序；
其中，所述开发恒河猴T细胞免疫组库的引物组包括26条上游引物和13条下游引物，其中，所述26条上游引物为如SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列，所述13条下游引物为如SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列。
3.如权利要求2所述的检测恒河猴T细胞免疫组库的高通量测序方法，其特征在于，所述多重PCR反应的体系中，26条上游引物等摩尔混合；13条下游引物等摩尔混合。
4.如权利要求2所述的检测恒河猴T细胞免疫组库的高通量测序方法，其特征在于，所述多重PCR反应的体系中，模板量为1～3ug/50ul体系。
5.如权利要求2所述的检测恒河猴T细胞免疫组库的高通量测序方法，其特征在于，所述多重PCR反应程序为：
6.如权利要求2所述的检测恒河猴T细胞免疫组库的高通量测序方法，其特征在于，所述多重PCR反应结束后，电泳，割胶回收片段长度为100-150bp的DNA片段。
7.一种构建恒河猴T细胞免疫组库的方法，其特征在于，包括：
提取恒河猴基因组DNA或总RNA；
采用开发恒河猴T细胞免疫组库的引物组，进行多重PCR反应；
将所得多重PCR产物进行高通量测序；
生物信息学分析，对高通量测序结果进行定量分析；
其中，所述开发恒河猴T细胞免疫组库的引物组包括26条上游引物和13条下游引物，其中，所述26条上游引物为如SEQ ID NO:1～SEQ ID NO:26所示的核苷酸序列，所述13条下游引物为如SEQ ID NO:27～SEQ ID NO:39所示的核苷酸序列。
8.一种如权利要求1所述的开发恒河猴T细胞免疫组库的引物组在检测恒河猴T细胞受体、对恒河猴T细胞受体进行高通量测序、构建恒河猴T细胞免疫组库、以及制备构建恒河猴T细胞免疫组库的试剂盒中的应用。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410843543.XA CN104673892B (zh) | 2014-12-30 | 2014-12-30 | 开发恒河猴t细胞免疫组库的引物组、高通量测序方法及应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410843543.XA CN104673892B (zh) | 2014-12-30 | 2014-12-30 | 开发恒河猴t细胞免疫组库的引物组、高通量测序方法及应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104673892A CN104673892A (zh) | 2015-06-03 |
| CN104673892B true CN104673892B (zh) | 2017-07-25 |
Family
ID=53309478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410843543.XA Expired - Fee Related CN104673892B (zh) | 2014-12-30 | 2014-12-30 | 开发恒河猴t细胞免疫组库的引物组、高通量测序方法及应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104673892B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105095687B (zh) * | 2015-06-26 | 2018-10-26 | 南方科技大学 | 一种免疫组库数据分析方法及终端 |
| CN105087789B (zh) * | 2015-08-10 | 2018-08-03 | 北京吉因加科技有限公司 | 一种检测血浆cfDNA中BCR和TCR免疫组库的方法 |
| CN105154440B (zh) * | 2015-08-14 | 2016-11-30 | 深圳市瀚海基因生物科技有限公司 | 一种基于高通量测序构建白血病微小残留病灶tcr文库的多重pcr引物和方法 |
| CN106834472A (zh) * | 2015-08-14 | 2017-06-13 | 深圳市瀚海基因生物科技有限公司 | Bcr多样性检测试剂盒及应用 |
| CN107955831A (zh) * | 2016-10-13 | 2018-04-24 | 深圳华大基因研究院 | 用于淋巴细胞定量检测的标记物和淋巴细胞定量检测方法 |
| CN109554440A (zh) * | 2018-12-26 | 2019-04-02 | 山东艾克韦生物技术有限公司 | 多重引物组及利用该引物组基于高通量测序构建人类t细胞免疫组库的方法 |
| CN109972209B (zh) * | 2019-01-18 | 2022-05-03 | 南开大学 | 基于乳液pcr的一步法将单链抗体库重排为全长抗体库的方法 |
| CN111261226B (zh) * | 2020-03-12 | 2020-10-27 | 江苏先声医学诊断有限公司 | 基于ngs的微小残留病变自动化测序分析方法及装置 |
| CN116640855A (zh) * | 2023-04-13 | 2023-08-25 | 华南农业大学 | 一组用于单细胞v（d）j测序的猪源tcr序列引物及其应用 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2062982A1 (fr) * | 2007-11-26 | 2009-05-27 | ImmunID | Procédé d'étude de la diversité combinatoire V(D)J |
| CN103710454B (zh) * | 2013-12-31 | 2016-02-17 | 南方科技大学 | Tcr或bcr高通量测序的方法及利用标签序列矫正多重pcr引物偏差的方法 |
-
2014
- 2014-12-30 CN CN201410843543.XA patent/CN104673892B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN104673892A (zh) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104673892B (zh) | 开发恒河猴t细胞免疫组库的引物组、高通量测序方法及应用 | |
| Robins | Immunosequencing: applications of immune repertoire deep sequencing | |
| US11053543B2 (en) | System and methods for massively parallel analysis of nucleic acids in single cells | |
| Afik et al. | Targeted reconstruction of T cell receptor sequence from single cell RNA-seq links CDR3 length to T cell differentiation state | |
| US11041202B2 (en) | Method of identifying human compatible T cell receptors specific for an antigenic target | |
| CN103710454B (zh) | Tcr或bcr高通量测序的方法及利用标签序列矫正多重pcr引物偏差的方法 | |
| CN106086013B (zh) | 一种用于核酸富集捕获的探针及设计方法 | |
| CN105087789B (zh) | 一种检测血浆cfDNA中BCR和TCR免疫组库的方法 | |
| US20120058902A1 (en) | Method of measuring adaptive immunity | |
| Boyd et al. | High‐throughput DNA sequencing analysis of antibody repertoires | |
| CN105506746A (zh) | 构建可变区测序文库的方法及确定可变区核酸序列的方法 | |
| US20210317522A1 (en) | Phenotypic and molecular characterisation of single cells | |
| CN107779495B (zh) | T细胞抗原受体多样性测序文库的构建方法以及试剂盒 | |
| CN106755410A (zh) | 一种基于高通量测序同时检测t细胞和b细胞免疫组库的方法 | |
| WO2019180271A1 (en) | Method of diagnosing celiac disease | |
| Wang et al. | A comprehensive analysis of the T and B lymphocytes repertoire shaped by HIV vaccines | |
| JP2019521642A (ja) | 免疫レパートリー発掘 | |
| CN107058484B (zh) | 一种应用于高通量测序同时检测t细胞和b细胞免疫组库的引物组合及试剂盒 | |
| Tu | Recovery of T cell receptor variable sequences from 3'barcoded single-cell RNA sequencing libraries | |
| CN110139952B (zh) | 骆驼科抗体可变区免疫组库构建的引物组合及应用 | |
| Morgan | Integrated, single-cell analysis of transcriptional phenotype and clonotypic identity | |
| Hu | Identifying, Investigating, and Modulating Anti-CD19 CAR T Cells for Large B-Cell Lymphoma | |
| Heather | T-cell receptor repertoire sequencing in health and disease | |
| CA3181835A1 (en) | Lymphocyte clonality determination |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170725 Termination date: 20181230 |