Brugada syndrome in Japan and Europe: a genome-wide association study reveals shared genetic architecture and new risk loci.

BACKGROUND AND AIMS: Brugada syndrome (BrS) is an inherited arrhythmia with a higher disease prevalence and more lethal arrhythmic events in Asians than in Europeans. Genome-wide association studies (GWAS) have revealed its polygenic architecture mainly in European populations. The aim of this study was to identify novel BrS-associated loci and to compare allelic effects across ancestries. METHODS: A GWAS was conducted in Japanese participants, involving 940 cases and 1634 controls, followed by a cross-ancestry meta-analysis of Japanese and European GWAS (total of 3760 cases and 11 635 controls). The novel loci were characterized by fine-mapping, gene expression, and splicing quantitative trait associations in the human heart. RESULTS: The Japanese-specific GWAS identified one novel locus near ZSCAN20 (P = 1.0 × 10-8), and the cross-ancestry meta-analysis identified 17 association signals, including six novel loci. The effect directions of the 17 lead variants were consistent (94.1%; P for sign test = 2.7 × 10-4), and their allelic effects were highly correlated across ancestries (Pearson's R = .91; P = 2.9 × 10-7). The genetic risk score derived from the BrS GWAS of European ancestry was significantly associated with the risk of BrS in the Japanese population [odds ratio 2.12 (95% confidence interval 1.94-2.31); P = 1.2 × 10-61], suggesting a shared genetic architecture across ancestries. Functional characterization revealed that a lead variant in CAMK2D promotes alternative splicing, resulting in an isoform switch of calmodulin kinase II-δ, favouring a pro-inflammatory/pro-death pathway. CONCLUSIONS: This study demonstrates novel susceptibility loci implicating potentially novel pathogenesis underlying BrS. Despite differences in clinical expressivity and epidemiology, the polygenic architecture of BrS was substantially shared across ancestries.

A diverse ancestrally-matched reference panel increases genotype imputation accuracy in a underrepresented population.

Variant imputation, a common practice in genome-wide association studies, relies on reference panels to infer unobserved genotypes. Multiple public reference panels are currently available with variations in size, sequencing depth, and represented populations. Currently, limited data exist regarding the performance of public reference panels when used in an imputation of populations underrepresented in the reference panel. Here, we compare the performance of various public reference panels: 1000 Genomes Project, Haplotype Reference Consortium, GenomeAsia 100 K, and the recent Trans-Omics for Precision Medicine (TOPMed) program, when used in an imputation of samples from the Thai population. Genotype yields were assessed, and imputation accuracies were examined by comparison with high-depth whole genome sequencing data of the same sample. We found that imputation using the TOPMed panel yielded the largest number of variants (~ 271 million). Despite being the smallest in size, GenomeAsia 100 K achieved the best imputation accuracy with a median genotype concordance rate of 0.97. For rare variants, GenomeAsia 100 K also offered the best accuracy, although rare variants were less accurately imputable than common variants (30.3% reduction in concordance rates). The high accuracy observed when using GenomeAsia 100 K is likely attributable to the diverse representation of populations genetically similar to the study cohort emphasizing the benefits of sequencing populations classically underrepresented in human genomics.

Metagenomic analysis of viral genes integrated in whole genome sequencing data of Thai patients with Brugada syndrome.

Brugada syndrome (BS) is an autosomal dominant inheritance cardiac arrhythmia disorder associated with sudden death in young adults. Thailand has the highest prevalence of BS worldwide, and over 60% of patients with BS still have unclear disease etiology. Here, we performeda new viral metagenome analysis pipeline called VIRIN and validated it with whole genome sequencing (WGS) data of HeLa cell lines and hepatocellular carcinoma. Then the VIRIN pipelinewas applied to identify viral integration positions from unmapped WGS data of Thai males, including 100 BS patients (case) and 100 controls. Even though the sample preparation had noviral enrichment step, we can identify several virus genes from our analysis pipeline. The predominance of human endogenous retrovirus K (HERV-K) viruses was found in both cases andcontrols by blastn and blastx analysis. This study is the first report on the full-length HERV-Kassembled genomes in the Thai population. Furthermore, the HERV-K integration breakpointpositions were validated and compared between the case and control datasets. Interestingly,Brugada cases contained HERV-K integration breakpoints at promoters five times more oftenthan controls. Overall, the highlight of this study is the BS-specific HERV-K breakpoint positionsthat were found at the gene coding region "NBPF11" (n = 9), "NBPF12" (n = 8) and longnon-coding RNA (lncRNA) "PCAT14" (n = 4) region. The genes and the lncRNA have been reported to be associated with congenital heart and arterial diseases. These findings provide another aspect of the BS etiology associated with viral genome integrations within the humangenome.

High prevalence of circulating DS-1-like human rotavirus A and genotype diversity in children with acute gastroenteritis in Thailand from 2016 to 2019.

BACKGROUND: Human rotavirus A (RVA) infection is the primary cause of acute gastroenteritis (AGE) in infants and young children worldwide, especially in children under 5 years of age and is a major public health problem causing severe diarrhea in children in Thailand. This study aimed to investigate the prevalence, genotype diversity, and molecular characterization of rotavirus infection circulating in children under 15 years of age diagnosed with AGE in Thailand from January 2016 to December 2019. METHODS: A total of 2,001 stool samples were collected from children with gastroenteritis (neonates to children <15 years of age) and tested for RVA by real-time polymerase chain reaction (RT-PCR). Amplified products were sequenced and submitted to an online genotyping tool for analysis. RESULTS: Overall, 301 (15.0%) stool samples were positive for RVA. RVA occurred most frequently among children aged 0-24 months. The seasonal incidence of rotavirus infection occurred typically in Thailand during the winter months (December-March). The G3P[8] genotype was identified as the most prevalent genotype (33.2%, 100/301), followed by G8P[8] (10.6%, 32/301), G9P[8] (6.3%, 19/301), G2P[4] (6.0%, 18/301), and G1P[6] (5.3%, 16/301). Uncommon G and P combinations such as G9P[4], G2P[8], G3P[4] and G3P[9] were also detected at low frequencies. In terms of genetic backbone, the unusual DS-1-like G3P[8] was the most frequently detected (28.2%, 85/301), and the phylogenetic analysis demonstrated high nucleotide identity with unusual DS-1-like G3P[8] detected in Thailand and several countries. CONCLUSIONS: A genetic association between RVA isolates from Thailand and other countries ought to be investigated given the local and global dissemination of rotavirus as it is crucial for controlling viral gastroenteritis, and implications for the national vaccination programs.

Phenotype prediction and characterization of 25 pharmacogenes in Thais from whole genome sequencing for clinical implementation.

Publicly available pharmacogenomics (PGx) databases enable translation of genotype data into clinically actionable information. As variation within pharmacogenes is population-specific, this study investigated the spectrum of 25 clinically relevant pharmacogenes in the Thai population (n = 291) from whole genome sequencing. The bioinformatics tool Stargazer was used for phenotype prediction, through assignment of alleles and detection of structural variation. Known and unreported potentially deleterious PGx variants were identified. Over 25% of Thais carried a high-risk diplotype in CYP3A5, CYP2C19, CYP2D6, NAT2, SLCO1B1, and UGT1A1. CYP2D6 structural variants accounted for 83.8% of all high-risk diplotypes. Of 39 known PGx variants identified, six variants associated with adverse drug reactions were common. Allele frequencies of CYP3A5*3 (rs776746), CYP2B6*6 (rs2279343), and NAT2 (rs1041983) were significantly higher in Thais than East-Asian and global populations. 121 unreported variants had potential to exert clinical impact, majority were rare and population-specific, with 60.3% of variants absent from gnomAD database. This study demonstrates the population-specific variation in clinically relevant pharmacogenes, the importance of CYP2D6 structural variation detection in the Thai population, and potential of unreported variants in explaining drug response. These findings are essential in development of dosing guidelines, PGx testing, clinical trials, and drugs.

Common and rare susceptibility genetic variants predisposing to Brugada syndrome in Thailand.

BACKGROUND: Mutations in SCN5A are rarely found in Thai patients with Brugada syndrome (BrS). Recent evidence suggested that common genetic variations may underlie BrS in a complex inheritance model. OBJECTIVE: The purpose of this study was to find common and rare/low-frequency genetic variants predisposing to BrS in persons in Thailand. METHODS: We conducted a genome-wide association study (GWAS) to explore the association of common variants in 154 Thai BrS cases and 432 controls. We sequenced SCN5A in 131 cases and 205 controls. Variants were classified according to current guidelines, and case-control association testing was performed for rare and low-frequency variants. RESULTS: Two loci were significantly associated with BrS. The first was near SCN5A/SCN10A (lead marker rs10428132; odds ratio [OR] 2.4; P = 3 × 10-10). Conditional analysis identified a novel independent signal in the same locus (rs6767797; OR 2.3; P = 2.7 × 10-10). The second locus was near HEY2 (lead marker rs3734634; OR 2.5; P = 7 × 10-9). Rare (minor allele frequency [MAF] <0.0001) coding variants in SCN5A were found in 8 of the 131 cases (6.1% in cases vs 2.0% in controls; P = .046; OR 3.3; 95% confident interval [CI] 1.0-11.1), but an enrichment of low-frequency (MAF<0.001 and >0.0001) variants also was observed in cases, with 1 variant (SCN5A: p.Arg965Cys) detected in 4.6% of Thai BrS patients vs 0.5% in controls (P = 0.015; OR 9.8; 95% CI 1.2-82.3). CONCLUSION: The genetic basis of BrS in Thailand includes a wide spectrum of variant frequencies and effect sizes. As previously shown in European and Japanese populations, common variants near SCN5A and HEY2 are associated with BrS in the Thai population, confirming the transethnic transferability of these 2 major BrS loci.

Genotype and clinical characteristics of congenital long QT syndrome in Thailand.

BACKGROUND: Congenital long QT syndrome (LQTS) is an inheritable arrhythmic disorder which is linked to at least 17 genes. The clinical characteristics and genetic mutations may be variable among different population groups and they have not yet been studied in Thai population. METHODS: Clinical characteristics were retrospectively reviewed from children and young adults with congenital long QT syndrome whose blood samples were sent for genotyping during 1998-2017. Sangers sequencing was used to sequentially identify KCNQ1 or KCNH2 genetic variants. Whole exome sequencing (WES) was used to identify variants in all other known LQTS genes. RESULTS: Of the 20 subjects (17 families), 45% were male, mean QTc was 550.3 ±â€¯68.8 msec (range 470-731 msec) and total Schwartz's score was 5.6 ±â€¯1.2 points (range 3-8 points). Fifty percent of patients had events at rest, 30% had symptoms after adrenergic mediated events, and 20% were asymptomatic. We discovered pathogenic and likely pathogenic genetic variants in KCNQ1, KCNH2, and SCN5A in 6 (35%), 4 (24%), and 2 (12%) families, respectively. One additional patient had variance of unknown significance (VUS) in KCNH2 and another one in ANK2. No pathogenic genetic variant was found in 3 patients (18%). Most patients received beta-blocker and 9 (45%) had ICD implanted. LQT1 patients were either asymptomatic or had stress-induced arrhythmia. Most of the LQT2 and LQT3 patients developed symptoms at rest or during sleep. CONCLUSIONS: Our patients with LQTS were mostly symptomatic at presentation. The genetic mutations were predominantly in LQT1, LQT2, and LQT3 genes.

Respiratory syncytial virus genotypes NA1, ON1, and BA9 are prevalent in Thailand, 2012-2015.

Respiratory syncytial virus (RSV) causes acute lower respiratory tract infection in infants and young children worldwide. To investigate the RSV burden in Thailand over four consecutive years (January 2012 to December 2015), we screened 3,306 samples obtained from children ≤5 years old with acute respiratory tract infection using semi-nested reverse-transcription polymerase chain reaction (RT-PCR). In all, 8.4% (277/3,306) of the specimens tested positive for RSV, most of which appeared in the rainy months of July to November. We then genotyped RSV by sequencing the G glycoprotein gene and performed phylogenetic analysis to determine the RSV antigenic subgroup. The majority (57.4%, 159/277) of the RSV belonged to subgroup A (RSV-A), of which NA1 genotype was the most common in 2012 while ON1 genotype became prevalent the following year. Among samples tested positive for RSV-B subgroup B (RSV-B) (42.6%, 118/277), most were genotype BA9 (92.6%, 87/94) with some BA10 and BA-C. Predicted amino acid sequence from the partial G region showed highly conserved N-linked glycosylation site at residue N237 among all RSV-A ON1 strains (68/68), and at residues N296 (86/87) and N310 (87/87) among RSV-B BA9 strains. Positive selection of key residues combined with notable sequence variations on the G gene contributed to the continued circulation of this rapidly evolving virus.

High prevalence of coxsackievirus A2 in children with herpangina in Thailand in 2015.

Coxsackievirus (CV) is a member of the genus Enterovirus and the family Picornaviridae. CV infection can cause herpangina, a disease characterized by multiple ulcers on the tonsils and soft palate affecting mostly young children. CV strains are categorized by serotypes. Unfortunately, serotypes responsible for infections in patients are often undetermined. This knowledge gap partly contributes to the ineffective prevention and control of CV-associated herpangina in Southeast Asia. To characterize the viral etiology of children presented with herpangina, 295 throat swabs were tested for human enterovirus infection. Using RT-PCR specific for the viral 5'UTR/VP2 and the VP1 regions, two most frequent CV types found in these samples were CV-A2 (33.33%, 40/120) and CV-A4 (15.8%, 19/120). Phylogenetic analysis of the VP1 gene demonstrated that the CV-A2 strains in this study not only were closely related to those previously identified in Asia and Europe, but the majority clustered into a distinct group. Thus, infection predominantly by CV-A2 and CV-A4 caused herpangina in 2015 in Thailand.

Complete genome sequence analysis of enterovirus 71 isolated from children with hand, foot, and mouth disease in Thailand, 2012-2014.

The complete genomic sequences of 14 enterovirus 71 (EV71) strains isolated from children with hand, foot, and mouth disease in Thailand from 2012 to 2014 were determined and compared to enterovirus group A prototypes. Phylogenetic analysis revealed that 13 strains resembled the B5 subgroup, while one strain from a fatal case designated THA_1219 belonged to the C4 subgroup. Similarity plot and bootscan analyses suggested that THA_1219 underwent recombination in the P2 and P3 regions. Full-genome data from this work will contribute to the study of evolution dynamics of EV71.

PREVALENCE OF HUMAN ENTEROVIRUS AMONG PATIENTS WITH HAND, FOOT, AND MOUTH DISEASE AND HERPANGINA IN THAILAND, 2013.

Human enterovirus (EV) infection causes hand, foot, and mouth disease (HFMD) and herpangina (HA). We studied the prevalence of enterovirus (EV) among patients with HFMD and HA in Thailand during 2013. We conducted a study in archived specimens of patients sent for screening for enterovirus. A total of 203 clinical specimens from 184 individuals with painful blister in the oropharynx and on the palms, soles, knees, elbows or buttock were examined by semi-nested polymerase chain reaction (PCR) for the 5'UTR and VP1 genes of EV. Eighty-six samples were positive: EV71 was detected in 14 (30%), CV-A8 in 12 (26%) and CV-A16 in 10 (21%). Classification of EV species detected revealed that 46 specimens were EV-A, 14 specimens were EV-B, 1 specimen was EV-D, and 16 specimens were positive for unclassified enterovirus. The majority of individuals with EV infection were aged 2-6 years. Multiple EV-A serotypes were detected among HFMD and HA patients in our study.

Prevalence and characterization of enterovirus infections among pediatric patients with hand foot mouth disease, herpangina and influenza like illness in Thailand, 2012.

Hand, foot, and mouth disease (HFMD) and herpangina are common infectious diseases caused by several genotypes of human enterovirus species A and frequently occurring in young children. This study was aimed at analyzing enteroviruses from patients with these diseases in Thailand in 2012. Detection and genotype determination of enteroviruses were accomplished by reverse transcription-polymerase chain reaction and sequencing of the VP1 region. Enterovirus-positive samples were differentiated into 17 genotypes (coxsackievirus A4 (CAV4), A5, A6, A8, A9, A10, A12, A16, A21, B1, B2, B4, B5, echovirus 7, 16, 25 and Enterovirus 71). The result showed CAV6 (33.5%), followed by CAV16 (9.4%) and EV71 (8.8%) as the most frequent genotypes in HFMD, CAV8 (19.3%) in herpangina and CAV6 (1.5%) in influenza like illness. Enterovirus infections were most prevalent during July with 34.4% in HFMD, 39.8% in herpangina and 1.6% in ILI. The higher enterovirus infection associated with HFMD and herpangina occurred in infants over one year-old. This represents the first report describing the circulation of multiple enteroviruses in Thailand.

Epidemiology and seroepidemiology of human enterovirus 71 among Thai populations.

BACKGROUND: Human enterovirus 71 (EV71) is an important pathogen caused large outbreaks in Asian-Pacific region with severe neurological complications and may lead to death in young children. Understanding of the etiological spectrum and epidemic changes of enterovirus and population's immunity against EV71 are crucial for the implementation of future therapeutic and prophylactic intervention. RESULTS: A total of 1,182 patients who presented with the symptoms of hand foot and mouth disease (67.3%) or herpangina (HA) (16.7%) and admitted to the hospitals during 2008-2013 were tested for enterovirus using pan-enterovirus PCR targeting 5'-untranslated region and specific PCR for viral capsid protein 1 gene. Overall, 59.7% were pan-enterovirus positive comprising 9.1% EV71 and 31.2% coxsackievirus species A (CV-A) including 70.5% CV-A6, 27.6% CV-A16, 1.1% CV-A10, and 0.8% CV-A5. HFMD and HA occurred endemically during 2008-2011. The number of cases increased dramatically in June 2012 with the percentage of the recently emerged CV-A6 significantly rose to 28.4%. Co-circulation between different EV71 genotypes was observed during the outbreak. Total of 161 sera obtained from healthy individuals were tested for neutralizing antibodies (NAb) against EV71 subgenotype B5 (EV71-B5) using microneutralization assay. The seropositive rate of EV71-B5 was 65.8%. The age-adjusted seroprevalence for individuals was found to be lowest in children aged >6 months to 2 years (42.5%). The seropositive rate remained relatively low in preschool children aged > 2 years to 6 years (48.3%) and thereafter increased sharply to more than 80% in individuals aged > 6 years. CONCLUSIONS: This study describes longitudinal data reflecting changing patterns of enterovirus prevalence over 6 years and demonstrates high seroprevalences of EV71-B5 NAb among Thai individuals. The rate of EV71 seropositive increased with age but without gender-specific significant difference. We identified that relative lower EV71 seropositive rate in early 2012 may demonstrate widely presented of EV71-B5 in the population before account for a large outbreak scale epidemic occurred in 2012 with due to a relatively high susceptibility of the younger population.