METHODS: We designed a 32-SNP panel for PGx testing in clinical laboratories. The variants were selected using the clinical annotations of the Pharmacogenomics Knowledgebase (PharmGKB) and include polymorphisms of CYP2C9, CYP2C19, CYP2D6, CYP3A5 and VKORC1 genes. The CYP2D6 gene allele quantification was determined simultaneously with TaqMan copy number assays targeting intron 2 and exon 9 regions. The genotyping results showed high call rate accuracy according to concordance with genotypes identified by independent analyses on Sequenome massarray and droplet digital PCR. Furthermore, 506 genomic samples across three major ethnic groups of Singapore (Malay, Indian and Chinese) were analysed on our workflow.
RESULTS: We found that 98% of our study subjects carry one or more CPIC actionable variants. The major alleles detected include CYP2C9*3, CYP2C19*2, CYP2D6*10, CYP2D6*36, CYP2D6*41, CYP3A5*3 and VKORC1*2. These translate into a high percentage of intermediate (IM) and poor metabolizer (PM) phenotypes for these genes in our population.
CONCLUSION: Genotyping may be useful to identify patients who are prone to drug toxicity with standard doses of drug therapy in our population. The simplicity and robustness of this PGx panel is highly suitable for use in a clinical laboratory.
RESULTS: We analyzed the whole-genome deep sequencing data (~ 30×) of five native trios from Peninsular Malaysia and North Borneo, and characterized the genomic variants, including single nucleotide variants (SNVs), small insertions and deletions (indels) and copy number variants (CNVs). We discovered approximately 6.9 million SNVs, 1.2 million indels, and 9000 CNVs in the 15 samples, of which 2.7% SNVs, 2.3% indels and 22% CNVs were novel, implying the insufficient coverage of population diversity in existing databases. We identified a higher proportion of novel variants in the Orang Asli (OA) samples, i.e., the indigenous people from Peninsular Malaysia, than that of the North Bornean (NB) samples, likely due to more complex demographic history and long-time isolation of the OA groups. We used the pedigree information to identify de novo variants and estimated the autosomal mutation rates to be 0.81 × 10- 8 - 1.33 × 10- 8, 1.0 × 10- 9 - 2.9 × 10- 9, and ~ 0.001 per site per generation for SNVs, indels, and CNVs, respectively. The trio-genomes also allowed for haplotype phasing with high accuracy, which serves as references to the future genomic studies of OA and NB populations. In addition, high-frequency inherited CNVs specific to OA or NB were identified. One example is a 50-kb duplication in DEFA1B detected only in the Negrito trios, implying plausible effects on host defense against the exposure of diverse microbial in tropical rainforest environment of these hunter-gatherers. The CNVs shared between OA and NB groups were much fewer than those specific to each group. Nevertheless, we identified a 142-kb duplication in AMY1A in all the 15 samples, and this gene is associated with the high-starch diet. Moreover, novel insertions shared with archaic hominids were identified in our samples.
CONCLUSION: Our study presents a full catalogue of the genome variants of the native Malaysian populations, which is a complement of the genome diversity in Southeast Asians. It implies specific population history of the native inhabitants, and demonstrated the necessity of more genome sequencing efforts on the multi-ethnic native groups of Malaysia and Southeast Asia.