Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2015 Oct 1:7:90.
doi: 10.1186/s13073-015-0211-x.

Concept and design of a genome-wide association genotyping array tailored for transplantation-specific studies

Yun R Li  1   2 Jessica van Setten  3 Shefali S Verma  4 Yontao Lu  5 Michael V Holmes  6 Hui Gao  2   6 Monkol Lek  7   8 Nikhil Nair  2   6 Hareesh Chandrupatla  2   6 Baoli Chang  2   6 Konrad J Karczewski  7   8 Chanel Wong  2   6 Maede Mohebnasab  2 Eyas Mukhtar  2   6 Randy Phillips  2   6 Vinicius Tragante  3 Cuiping Hou  2 Laura Steel  2   6 Takesha Lee  2   6 James Garifallou  2 Toumy Guettouche  2 Hongzhi Cao  9   10 Weihua Guan  11 Aubree Himes  2   6 Jacob van Houten  2 Andrew Pasquier  2 Reina Yu  2 Elena Carrigan  2 Michael B Miller  12 David Schladt  13 Abdullah Akdere  1 Ana Gonzalez  1 Kelsey M Llyod  1 Daniel McGinn  1 Abhinav Gangasani  1 Zach Michaud  1 Abigail Colasacco  1 James Snyder  2 Kelly Thomas  2 Tiancheng Wang  2 Baolin Wu  11 Alhusain J Alzahrani  14 Amein K Al-Ali  15 Fahad A Al-Muhanna  15 Abdullah M Al-Rubaish  15 Samir Al-Mueilo  15 Dimitri S Monos  2   16 Barbara Murphy  17 Kim M Olthoff  6 Cisca Wijmenga  18 Teresa Webster  5 Malek Kamoun  16 Suganthi Balasubramanian  19 Matthew B Lanktree  6 William S Oetting  20 Pablo Garcia-Pavia  21 Daniel G MacArthur  7   8 Paul I W de Bakker  22 Hakon Hakonarson  2 Kelly A Birdwell  23 Pamala A Jacobson  24 Marylyn D Ritchie  4 Folkert W Asselbergs  3   25   26 Ajay K Israni  27 Abraham Shaked  6 Brendan J Keating  28   29   30   31
Affiliations

Concept and design of a genome-wide association genotyping array tailored for transplantation-specific studies

Yun R Li et al. Genome Med. .

Abstract

Background: In addition to HLA genetic incompatibility, non-HLA difference between donor and recipients of transplantation leading to allograft rejection are now becoming evident. We aimed to create a unique genome-wide platform to facilitate genomic research studies in transplant-related studies. We designed a genome-wide genotyping tool based on the most recent human genomic reference datasets, and included customization for known and potentially relevant metabolic and pharmacological loci relevant to transplantation.

Methods: We describe here the design and implementation of a customized genome-wide genotyping array, the 'TxArray', comprising approximately 782,000 markers with tailored content for deeper capture of variants across HLA, KIR, pharmacogenomic, and metabolic loci important in transplantation. To test concordance and genotyping quality, we genotyped 85 HapMap samples on the array, including eight trios.

Results: We show low Mendelian error rates and high concordance rates for HapMap samples (average parent-parent-child heritability of 0.997, and concordance of 0.996). We performed genotype imputation across autosomal regions, masking directly genotyped SNPs to assess imputation accuracy and report an accuracy of >0.962 for directly genotyped SNPs. We demonstrate much higher capture of the natural killer cell immunoglobulin-like receptor (KIR) region versus comparable platforms. Overall, we show that the genotyping quality and coverage of the TxArray is very high when compared to reference samples and to other genome-wide genotyping platforms.

Conclusions: We have designed a comprehensive genome-wide genotyping tool which enables accurate association testing and imputation of ungenotyped SNPs, facilitating powerful and cost-effective large-scale genotyping of transplant-related studies.

PubMed Disclaimer

Figures

Fig. 1
Fig. 1
Key content modules included in the TxArray. Marker modules are divided into those in common with the UK Biobank Axiom Array, Axiom Biobank Genotyping Array, and those unique to the TxArray along with module descriptions and counts of unique markers in each category. The goal of preserving the specific markers already covered on the two existing biobank arrays is maximize shared coverage for future cross-study joint and meta-analyses that may utilize these popular state-of-the art genotyping platforms
Fig. 2
Fig. 2
Modular content in the TxArray targeting the MHC and custom-designed for transplant-focused association studies. a MHC modules included in the TxArray and the total count of markers either directly tiled or tagged. Note that the marker counts are not unique as some modules overlap. A total of 10,820 and 13,428 unique markers across the MHC and extended MHC are included in the TxArray. b Major categories of transplant-specific markers included in the array. See methods and supplemental tables for details of design and content considerations
Fig. 3
Fig. 3
Comparison of coverage of 1000 genomes project reference panel between TxArray and other genome-wide genotyping platforms for variants with MAF >0.05 (a) and >0.01 (b). Coverage (ordinate) for the 1000 Genomes Project Phase I integrated reference panel was assessed using maximum r 2 (abscissa), at an MAF cutoff of 0.05 (a) and 0.01 (b). Populations included: (1) African ancestry (AAM): Yoruba in Ibadan, Nigeria (YRI) and Americans of African Ancestry in SouthWest, USA (ASW); (2) Admixed American (AMR): Colombians from Medellin, Colombia (CLM), Mexican Ancestry from Los Angeles USA (MXL), and Puerto Ricans from Puerto Rico (PUR); (3) Asian (ASN): Han Chinese in Beijing (CHB), Southern Han Chinese (CHS), Japanese in Tokyo, Japan (JPT); and (4) European ancestry (CTI): Utah residents with ancestry from Northern Western Europe (CTU), Central and Eastern European (CEU), and Toscani in Italia (TSI), as described in the HapMap and 1KGP. The platforms compared include the TxArray using 767,203 SNPs passing QC. ILMN_1M and ILMN660 refer to Illumina’s Infinium one million and the Illumina 660K genotyping platforms. Affy_6.0 refers to the Affymetrix 6.0 SNP chip containing approximately 906,600 SNPs
Fig. 4
Fig. 4
Comparison of coverage between TxArray and ILMN_1M genotyping platforms across exonic regions, the extended MHC and the KIR-encoding locus. a Coverage (ordinate) for all exonic markers and UTR region markers in the 1000 genomes reference panel was assessed using max r2 (abscissa), at an MAF cutoff of 0.05 (a) and 0.01 (b), in (1) European ancestry ((CEU) and Tuscany in Italia (TSI)); (2) African ancestry (AAM) (Yoruba in Ibadan, Nigeria (YRI)) and Americans of African Ancestry in SouthWest, USA (ASW); (3) Admixed American (AMR) (Colombians from Medellin, Colombia (CLM), Mexican Ancestry from Los Angeles, USA (MXL), and Puerto Ricans from Puerto Rico (PUR)); and (4) Asian (ASN) (Han Chinese in Beijing (CHB), Southern Han Chinese (CHS), Japanese in Tokyo, Japan (JPT)) HapMap and 1KGP individuals. The compared platforms include the TxArray using 767,203 SNPs that passed manufacturing and standard genotyping QC. ILMN_1M refer to Illumina’s Infinium one million SNP GWAS array. b Comparison of coverage across variants within KIR-encoding regions using the TxArray (TX) or the Illumina 1M (ILMN_1M) genotyping platforms across the four major HapMap populations (European (CTU): CEU+TSI; AAM: ASW+YRI; AMR: CLM+MXL+PUR; ASN: CHB+CHS+JPT). Coverage is based on mean r2 of variants included in the 1000 genomes phase I reference panel with a MAF of >0.01 (top) or >0.05 (bottom). KIR genes included: (KIR2DP1, KIR2DL1, KIR2DL2, KIR2DL3, KIR2DL4, KIR2DL5A, KIR2DL5B, KIR2DS1, KIR2DS2, KIR2DS3, KIR2DS4, KIR2DS5, KIR3DL1, KIR3DL2, KIR3DL3, KIR3DP1, KIR3DS1, KIR3DX1). Coverage was compared for either all KIR region markers (left) or only those in exonic regions (right). c Comparison of coverage across the extended MHC (25,500,000–34,000,000) using either the TxArray (TX) or the Illumina 1M (ILMN_1M) genotyping platforms across the four major HAPMAP populations (CTU: CEU/TSI; AAM: ASW/YRI; AMR: CLM/MXL/PUR; ASN: CHB/CHS/JPT). Coverage rate is calculated based on the mean achieved r2 for variants included in the 1000 Genomes Project (1KGP) Phase I reference panel with a MAF of >0.01 (left) or >0.05 (right)

Similar articles

Cited by

References

    1. Organ Procurement and Transplantation Network. Available at: http://optn.transplant.hrsa.gov.
    1. Stehlik J, Edwards LB, Kucheryavaya AY, Aurora P, Christie JD, Kirk R, et al. The Registry of the International Society for Heart and Lung Transplantation: twenty-seventh official adult heart transplant report--2010. J Heart Lung Transplant. 2010;29:1089–103. doi: 10.1016/j.healun.2010.08.007. - DOI - PubMed
    1. Burton CM, Iversen M, Carlsen J, Mortensen J, Andersen CB, Steinbrüchel D, et al. Acute cellular rejection is a risk factor for bronchiolitis obliterans syndrome independent of post-transplant baseline FEV1. J Heart Lung Transplant. 2009;28:888–93. doi: 10.1016/j.healun.2009.04.022. - DOI - PubMed
    1. Terasaki PI. Deduction of the fraction of immunologic and non-immunologic failure in cadaver donor transplants. Clin Transpl. 2003;449–52. - PubMed
    1. Kaplan B, Qazi Y, Wellen JR. Strategies for the management of adverse events associated with mTOR inhibitors. Transplant Rev (Orlando) 2014;28:126–33. doi: 10.1016/j.trre.2014.03.002. - DOI - PubMed

Publication types