Whole genome sequencing analysis of multidrug-resistant tuberculosis in Singapore, 2006-2018.

There were 290 multidrug-resistant (MDR)-TB cases diagnosed in Singapore from 2006 to 2018. Eighty-one percent were foreign-born. Spoligotyping and MIRU-VNTR methods identified 108 patients in 24 clusters. The Beijing spoligotype accounted for 22 clusters. Whole genome sequencing (WGS) analysis reduced the number of clustered patients and clusters to 43 and nine respectively. One MIRU cluster was redefined into three WGS clusters. All the clusters had foreign-born source cases. Forty percent of local-born, versus 9% of foreign-born, MDR-TB cases belonged to WGS clusters. WGS more accurately elucidated potential MDR-TB transmission which was overestimated by conventional genotyping methods in Singapore.

Risk for ACPA-positive rheumatoid arthritis is driven by shared HLA amino acid polymorphisms in Asian and European populations.

Previous studies have emphasized ethnically heterogeneous human leukocyte antigen (HLA) classical allele associations to rheumatoid arthritis (RA) risk. We fine-mapped RA risk alleles within the major histocompatibility complex (MHC) in 2782 seropositive RA cases and 4315 controls of Asian descent. We applied imputation to determine genotypes for eight class I and II HLA genes to Asian populations for the first time using a newly constructed pan-Asian reference panel. First, we empirically measured high imputation accuracy in Asian samples. Then we observed the most significant association in HLA-DRß1 at amino acid position 13, located outside the classical shared epitope (Pomnibus = 6.9 × 10(-135)). The individual residues at position 13 have relative effects that are consistent with published effects in European populations (His > Phe > Arg > Tyr ≅ Gly > Ser)--but the observed effects in Asians are generally smaller. Applying stepwise conditional analysis, we identified additional independent associations at positions 57 (conditional Pomnibus = 2.2 × 10(-33)) and 74 (conditional Pomnibus = 1.1 × 10(-8)). Outside of HLA-DRß1, we observed independent effects for amino acid polymorphisms within HLA-B (Asp9, conditional P = 3.8 × 10(-6)) and HLA-DPß1 (Phe9, conditional P = 3.0 × 10(-5)) concordant with European populations. Our trans-ethnic HLA fine-mapping study reveals that (i) a common set of amino acid residues confer shared effects in European and Asian populations and (ii) these same effects can explain ethnically heterogeneous classical allelic associations (e.g. HLA-DRB1*09:01) due to allele frequency differences between populations. Our study illustrates the value of high-resolution imputation for fine-mapping causal variants in the MHC.

A genome-wide association study of hepatitis B vaccine response in an Indonesian population reveals multiple independent risk variants in the HLA region.

We performed a two-stage genome-wide association study (GWAS) of antibody titer in 3614 hepatitis B vaccine recipients from Indonesia's Riau Archipelago, leading to the identification of at least three independent signals within the human leukocyte antigen (HLA) complex. These appear to implicate HLA-DR [rs3135363; P= 6.53 × 10(-22); odds ratio (OR) = 1.53, 95% confidence interval (CI) = 1.35-1.74]; HLA-DP, previously associated with the risk of chronic hepatitis B infection (rs9277535; P= 2.91 × 10(-12); OR = 0.72, 95% CI = 0.63-0.81); and a gene rich HLA Class III interval (rs9267665; P = 1.24 × 10(-17); OR = 2.05, CI = 1.64-2.57). The substantial overlap of these variants and those identified by GWAS of chronic hepatitis B infection confirms vaccine response as a model for infection, while suggesting that the vaccine is least effective in those most at risk of lifelong infection, following exposure to the virus.

Singapore Genome Variation Project: a haplotype map of three Southeast Asian populations.

The Singapore Genome Variation Project (SGVP) provides a publicly available resource of 1.6 million single nucleotide polymorphisms (SNPs) genotyped in 268 individuals from the Chinese, Malay, and Indian population groups in Southeast Asia. This online database catalogs information and summaries on genotype and phased haplotype data, including allele frequencies, assessment of linkage disequilibrium (LD), and recombination rates in a format similar to the International HapMap Project. Here, we introduce this resource and describe the analysis of human genomic variation upon agglomerating data from the HapMap and the Human Genome Diversity Project, providing useful insights into the population structure of the three major population groups in Asia. In addition, this resource also surveyed across the genome for variation in regional patterns of LD between the HapMap and SGVP populations, and for signatures of positive natural selection using two well-established metrics: iHS and XP-EHH. The raw and processed genetic data, together with all population genetic summaries, are publicly available for download and browsing through a web browser modeled with the Generic Genome Browser.

Identifying candidate causal variants via trans-population fine-mapping.

Genome-wide association studies have discovered and confirmed a large number of loci that are implicated with disease susceptibility and severity. Polymorphisms that emerged from these studies are mostly indirectly associated to the phenotype, and the natural progression is to identify the causal variants that are functionally responsible for these association signals. Long stretches of high linkage disequilibrium (LD) benefitted the initial discovery phase in a genome-wide scan, allowing commercial genotyping products with imperfect coverage to detect genomic regions genuinely associated with the phenotype. However, regions of high LD confound the fine-mapping phase, as markers that are perfectly correlated to the causal variants display similar evidence of phenotypic association, hampering the process of differentiating the functional polymorphisms from neighboring surrogates. Here, we explore the potential of integrating information across different populations for narrowing the candidate region that a causal variant resides in, and compare the efficacy of this process of trans-population fine-mapping with the extent of variation in patterns of LD between the populations. In addition, we explore two different strategies for pooling data across multiple populations for the purpose of prioritizing the rankings of the causal variants. Our results clearly establish the benefits of trans-population analysis in reducing the number of possible candidates for the causal variants, particularly in genomic regions displaying strong evidence of inter-population LD variation. Directly integrating the statistical evidence by summing the test statistics outperforms the standard meta-analytic procedure. These findings have direct relevance to the design and analysis of ongoing fine-mapping studies.

A genome-wide association study suggests contrasting associations in ACPA-positive versus ACPA-negative rheumatoid arthritis.

BACKGROUND: Rheumatoid arthritis (RA) can be divided into two major subsets based on the presence or absence of antibodies to citrullinated peptide antigens (ACPA). Until now, data from genome-wide association studies (GWAS) have only been published from ACPA-positive subsets of RA or from studies that have not separated the two subsets. The aim of the current study is to provide and compare GWAS data for both subsets. METHODS AND RESULTS: GWAS using the Illumina 300K chip was performed for 774 ACPA-negative patients with RA, 1147 ACPA-positive patients with RA and 1079 controls from the Swedish population-based case-control study EIRA. Imputation was performed which allowed comparisons using 1,723,056 single nucleotide polymorphisms (SNPs). No SNP achieved genome-wide significance (2.9 × 10⁻8) in the comparison between ACPA-negative RA and controls. A case-case association study was then performed between ACPA-negative and ACPA-positive RA groups. The major difference in this analysis was in the HLA region where 768 HLA SNPs passed the threshold for genome-wide significance whereas additional contrasting SNPs did not reach genome-wide significance. However, one SNP close to the RPS12P4 locus in chromosome 2 reached a p value of 2 × 106 and this locus can thus be considered as a tentative candidate locus for ACPA-negative RA. CONCLUSIONS: ACPA-positive and ACPA-negative RA display significant risk allele frequency differences which are mainly confined to the HLA region. The data provide further support for distinct genetic aetiologies of RA subsets and emphasise the need to consider them separately in genetic as well as functional studies of this disease.

Genomic copy number variations in three Southeast Asian populations.

Research on the role of copy number variations (CNVs) in the genetic risk of diseases in Asian populations has been hampered by a relative lack of reference CNV maps for Asian populations outside the East Asians. In this article, we report the population characteristics of CNVs in Chinese, Malay, and Asian Indian populations in Singapore. Using the Illumina Human 1M Beadchip array, we identify 1,174 CNV loci in these populations that corroborated with findings when the same samples were typed on the Affymetrix 6.0 platform. We identify 441 novel loci not previously reported in the Database of Genomic Variations (DGV). We observe a considerable number of loci that span all three populations and were previously unreported, as well as population-specific loci that are quite common in the respective populations. From this we observe the distribution of CNVs in the Asian Indian population to be considerably different from the Chinese and Malay populations. About half of the deletion loci and three-quarters of duplication loci overlap UCSC genes. Tens of loci show population differentiation and overlap with genes previously known to be associated with genetic risk of diseases. One of these loci is the CYP2A6 deletion, previously linked to reduced susceptibility to lung cancer.

TRAF1-C5 as a risk locus for rheumatoid arthritis--a genomewide study.

BACKGROUND: Rheumatoid arthritis has a complex mode of inheritance. Although HLA-DRB1 and PTPN22 are well-established susceptibility loci, other genes that confer a modest level of risk have been identified recently. We carried out a genomewide association analysis to identify additional genetic loci associated with an increased risk of rheumatoid arthritis. METHODS: We genotyped 317,503 single-nucleotide polymorphisms (SNPs) in a combined case-control study of 1522 case subjects with rheumatoid arthritis and 1850 matched control subjects. The patients were seropositive for autoantibodies against cyclic citrullinated peptide (CCP). We obtained samples from two data sets, the North American Rheumatoid Arthritis Consortium (NARAC) and the Swedish Epidemiological Investigation of Rheumatoid Arthritis (EIRA). Results from NARAC and EIRA for 297,086 SNPs that passed quality-control filters were combined with the use of Cochran-Mantel-Haenszel stratified analysis. SNPs showing a significant association with disease (P<1x10(-8)) were genotyped in an independent set of case subjects with anti-CCP-positive rheumatoid arthritis (485 from NARAC and 512 from EIRA) and in control subjects (1282 from NARAC and 495 from EIRA). RESULTS: We observed associations between disease and variants in the major-histocompatibility-complex locus, in PTPN22, and in a SNP (rs3761847) on chromosome 9 for all samples tested, the latter with an odds ratio of 1.32 (95% confidence interval, 1.23 to 1.42; P=4x10(-14)). The SNP is in linkage disequilibrium with two genes relevant to chronic inflammation: TRAF1 (encoding tumor necrosis factor receptor-associated factor 1) and C5 (encoding complement component 5). CONCLUSIONS: A common genetic variant at the TRAF1-C5 locus on chromosome 9 is associated with an increased risk of anti-CCP-positive rheumatoid arthritis.

Meta-analysis identifies common variants associated with body mass index in east Asians.

Multiple genetic loci associated with obesity or body mass index (BMI) have been identified through genome-wide association studies conducted predominantly in populations of European ancestry. We performed a meta-analysis of associations between BMI and approximately 2.4 million SNPs in 27,715 east Asians, which was followed by in silico and de novo replication studies in 37,691 and 17,642 additional east Asians, respectively. We identified ten BMI-associated loci at genome-wide significance (P < 5.0 × 10(-8)), including seven previously identified loci (FTO, SEC16B, MC4R, GIPR-QPCTL, ADCY3-DNAJC27, BDNF and MAP2K5) and three novel loci in or near the CDKAL1, PCSK1 and GP2 genes. Three additional loci nearly reached the genome-wide significance threshold, including two previously identified loci in the GNPDA2 and TFAP2B genes and a newly identified signal near PAX6, all of which were associated with BMI with P < 5.0 × 10(-7). Findings from this study may shed light on new pathways involved in obesity and demonstrate the value of conducting genetic studies in non-European populations.

Genome-wide association study identifies susceptibility loci for dengue shock syndrome at MICB and PLCE1.

Hypovolemic shock (dengue shock syndrome (DSS)) is the most common life-threatening complication of dengue. We conducted a genome-wide association study of 2,008 pediatric cases treated for DSS and 2,018 controls from Vietnam. Replication of the most significantly associated markers was carried out in an independent Vietnamese sample of 1,737 cases and 2,934 controls. SNPs at two loci showed genome-wide significant association with DSS. We identified a susceptibility locus at MICB (major histocompatibility complex (MHC) class I polypeptide-related sequence B), which was within the broad MHC region on chromosome 6 but outside the class I and class II HLA loci (rs3132468, P(meta) = 4.41 × 10(-11), per-allele odds ratio (OR) = 1.34 (95% confidence interval: 1.23-1.46)). We identified associated variants within PLCE1 (phospholipase C, epsilon 1) on chromosome 10 (rs3765524, P(meta) = 3.08 × 10(-10), per-allele OR = 0.80 (95% confidence interval: 0.75-0.86)). We identify two loci associated with susceptibility to DSS in people with dengue, suggesting possible mechanisms for this severe complication of dengue.

Genome-wide association identifies multiple ulcerative colitis susceptibility loci.

Ulcerative colitis is a chronic, relapsing inflammatory condition of the gastrointestinal tract with a complex genetic and environmental etiology. In an effort to identify genetic variation underlying ulcerative colitis risk, we present two distinct genome-wide association studies of ulcerative colitis and their joint analysis with a previously published scan, comprising, in aggregate, 2,693 individuals with ulcerative colitis and 6,791 control subjects. Fifty-nine SNPs from 14 independent loci attained an association significance of P < 10(-5). Seven of these loci exceeded genome-wide significance (P < 5 x 10(-8)). After testing an independent cohort of 2,009 cases of ulcerative colitis and 1,580 controls, we identified 13 loci that were significantly associated with ulcerative colitis (P < 5 x 10(-8)), including the immunoglobulin receptor gene FCGR2A, 5p15, 2p16 and ORMDL3 (orosomucoid1-like 3). We confirmed association with 14 previously identified ulcerative colitis susceptibility loci, and an analysis of acknowledged Crohn's disease loci showed that roughly half of the known Crohn's disease associations are shared with ulcerative colitis. These data implicate approximately 30 loci in ulcerative colitis, thereby providing insight into disease pathogenesis.

Association of skin barrier genes within the PSORS4 locus is enriched in Singaporean Chinese with early-onset psoriasis.

Psoriasis (OMIM#177900) is a common polygenic skin disorder affecting approximately 2% of the northern European population and 0.1% of the Han Chinese. Psoriasis patients suffer from chronic skin inflammation, manifested by erythematous scaly lesions. PSORS1-PSORS9 have been confirmed as psoriasis susceptibility loci in independent genetic studies on predominantly Caucasian populations, with psoriasis susceptibility loci (PSORS1, PSORS9) and additional loci at 9q33-34 and 2p22.3-11.2 reported in Han Chinese patients. In this study, we show the association of PSORS4 with psoriasis in Singaporean Chinese. Dense genotyping of single-nucleotide polymorphism-tagging candidate genes within the epidermal differentiation complex revealed significant association in the proximity of the involucrin gene (IVL); the strongest association was seen in early-onset psoriasis patients (P=0.0014). A follow-up genome-wide association screen localized the psoriasis susceptibility region to approximately 360 kb along chromosome 1 in the vicinity of IVL, small proline-rich region (SPRR) and proline-rich region 9 (PRR9) genes. The study of interactions between the causative variant(s) in this locus will provide insights into a possible role for epidermal barrier formation in the pathogenesis of psoriasis.