KRGDB: the large-scale variant database of 1722 Koreans based on whole genome sequencing

The efficacy and safety of eculizumab in patients and the role of C5 polymorphisms

Eculizumab is an orphan drug with indications for extremely rare autoimmune disorders. It is primarily prescribed for use in patients with paroxysmal nocturnal hemoglobinuria and atypical hemolytic uremic syndrome; but is also highly effective in the treatment of myasthenia gravis, among others. By binding to the C5 protein in the complement system, eculizumab effectively inhibits cellular hemolysis and autoimmune reactions. Despite this effective treatment, some patients reported no improvement in symptoms. Genetic sequencing revealed three distinct C5 mutations in the non-responders and these polymorphisms appeared to be most prevalent among Japanese, Korean and African populations. Here, we present an overview of the current and potential future applications of eculizumab, as well as the disadvantages of eculizumab treatment in patients with C5 polymorphisms.

Evaluating the TaqMan Jra-Genotyping Method for Rapidly Predicting the Presence of Anti-Jra Antibodies

Background

The Jra antigen is a high-prevalence red blood cell (RBC) antigen. Reports on cases of fatal hemolytic disease of the fetus and newborn and acute hemolytic transfusion reactions suggest that antibodies against Jra (anti-Jra) have potential clinical significance. Identifying anti-Jra is challenging owing to a lack of commercially available antisera. We developed an alternative approach to rapidly predict the presence of anti-Jra using the TaqMan single-nucleotide polymorphism (SNP)-genotyping method.

Methods

Residual peripheral blood samples from 10 patients suspected of having the anti-Jra were collected. Two samples with confirmed Jr(a-) RBCs and anti-Jra were used to validate the TaqMan genotyping assay by comparing the genotyping results with direct sequencing. The accuracy of the assay in predicting the presence of anti-Jra was verified through crossmatching with in-house Jr(a-) O+ RBCs.

Results

The TaqMan-genotyping method was validated with two Jr(a-) RBC- and anti-Jra-confirmed samples that showed concordant Jra genotyping and direct sequencing results. Jra genotyping for the remaining samples and crossmatching the serum samples with inhouse Jr(a-) O+ RBCs showed consistent results.

Conclusions

We validated a rapid, simple, accurate, and cost-effective method for predicting the presence of anti-Jra using a TaqMan-based SNP-genotyping assay. Implementing this method in routine practice in clinical laboratories will assist in solving difficult problems regarding alloantibodies to high-prevalence RBC antigens and ultimately aid in providing safe and timely transfusions and proper patient care.

Subset of retinoblastoma tumours is associated with BRCA1/2 mutations

BACKGROUND

We investigated the potential association between pathogenic BRCA1/2 variants and retinoblastoma pathogenicity.

METHODS

In this single-centre, retrospective case series, we performed hereditary cancer panel tests using blood samples for patients with retinoblastoma diagnosed between March 2017 and October 2021. Bioinformatics prediction tools were then used to conduct in silico pathogenicity assessments for patients with BRCA1/2 family variants, in addition to the American College of Medical Genetics and Genomics (ACMG) variant classification. One patient with a germline BRCA1 variant was analysed with whole-genome sequencing (WGS), mutational signature analysis and methylation analysis for RB1 and BRCA using the patient's tumour and blood samples.

RESULTS

Of 30 retinoblastoma patients who underwent panel sequencing, six (20%) were found to carry germline variants in the BRCA1/2 or BRIP1 genes. Among these six patients, two had pathogenic or likely pathogenic variants as per the ACMG variant classification. Additionally, three patients showed potential pathogenic BRCA1/2 family variants through further analysis with alternative bioinformatics prediction tools. In the WGS analysis of a tumour from a patient with a germline likely pathogenic BRCA1 variant in one allele, we observed the loss of one RB1 allele due to a large deletion. No somatic non-synonymous mutations or frameshift indels were detected in the RB1 locus of the remaining allele. This sample also showed BRCA1 gene promoter hypermethylation in the tumour, indicating additional epigenetic silencing.

CONCLUSION

This study demonstrated that some retinoblastoma patients harboured germline BRCA1/2 family variants, which may be associated with the development of retinoblastoma along with RB1 mutations.

Discovery of Pathogenic Variants Associated with Idiopathic Recurrent Pregnancy Loss Using Whole-Exome Sequencing

Idiopathic recurrent pregnancy loss (RPL) is defined as at least two pregnancy losses before 20 weeks of gestation. Approximately 5% of pregnant couples experience idiopathic RPL, which is a heterogeneous disease with various causes including hormonal, chromosomal, and intrauterine abnormalities. Although how pregnancy loss occurs is still unknown, numerous biological factors are associated with the incidence of pregnancy loss, including genetic variants. Whole-exome sequencing (WES) was conducted on blood samples from 56 Korean patients with RPL and 40 healthy controls. The WES data were aligned by means of bioinformatic analysis, and the detected variants were annotated using machine learning tools to predict the pathogenicity of protein alterations. Each indicated variant was confirmed using Sanger sequencing. A replication study was also conducted in 112 patients and 114 controls. The Variant Effect Scoring Tool, Combined Annotation Dependent Depletion tool, Sorting Intolerant from Tolerant annotation tool, and various databases detected 10 potential variants previously associated with spontaneous abortion genes in patients by means of a bioinformatic analysis of WES data. Several variants were detected in more than one patient. Interestingly, several of the detected genes were functionally clustered, including some with a secretory function (mucin 4; MUC4; rs200737893 G>A and hyaluronan-binding protein 2; HABP2; rs542838125 G>T), in which growth arrest-specific 2 Like 2 (GAS2L2; rs140842796 C>T) and dynamin 2 (DNM2; rs763894364 G>A) are functionally associated with cell protrusion and the cytoskeleton. ATP Binding Cassette Subfamily C Member 6 (ABCC6) was the only gene with two variants. HABP2 (rs542838125 G>T), MUC4 (rs200737893 G>A), and GAS2L2 (rs140842796 C>T) were detected in only the patient group in the replication study. The combination of WES and machine learning tools is a useful method to detect potential variants associated with RPL. Using bioinformatic tools, we found 10 potential variants in 9 genes. WES data from patients are needed to better understand the causes of RPL.

Genome-wide association study implicates the role of TBXAS1 in the pathogenesis of depressive symptoms among the Korean population

Although depression is an emerging disorder affecting many people worldwide, most genetic studies have been performed in European descent populations. Herein, a genome-wide association study (GWAS) was conducted in Korean population to elucidate the genomic loci associated with depressive symptoms. Two independent cohorts were used as discovery datasets, which consisted of 6474 (1484 cases and 4990 controls) and 1654 (557 cases and 1097 controls) Korean participants, respectively. The participants were divided into case and control groups based on the Beck Depression Inventory (BDI). Meta-analysis using the two cohorts revealed that rs6945590 was significantly associated with the risk of depressive symptoms [P = 2.83 × 10-8; odds ratio (OR) = 1.23; 95% confidence interval (CI): 1.15-1.33]. This association was validated in other independent cohorts which were another Korean cohort (258 cases and 1757 controls) and the East Asian study of the Psychiatric Genomics Consortium (PGC) (12,455 cases and 85,548 controls). The predicted expression levels of thromboxane A synthase 1 gene (TBXAS1), which encodes the enzyme thromboxane A synthase 1 and participates in the arachidonic acid (AA) cascade, was significantly decreased in the whole blood tissues of the participants with depressive symptoms. Furthermore, Mendelian randomization (MR) analysis showed a causal association between TBXAS1 expression and the risk of depressive symptoms. In conclusion, as the number of risk alleles (A) of rs6945590 increased, TBXAS1 expression decreased, which subsequently caused an increase in the risk of depressive symptoms.

Carrier frequency and incidence of alpha-mannosidosis: population database-based study-focus on the East Asian and Korean population

Background: Alpha-mannosidosis caused by mutations in the MAN2B1 gene is a rare genetic disorder characterized by physical abnormalities and intellectual disabilities. The objective of this study was to analyze the carrier frequency and estimated incidence of alpha-mannosidosis in East Asian populations, as limited data exists on its incidence in this group. Methods: In this study, a total of 125,748 exomes from the gnomAD database was analyzed. Additionally, 5,305 data from the KOVA and 1,722 data from the KRGDB, both representing Korean populations, were included. Results: The global carrier frequency of alpha-mannosidosis in gnomAD was 0.23%; the highest carrier frequency was observed in the Finnish at 0.49%, and East Asians had the second highest carrier frequency at 0.30%. Globally, the approximate incidence of alpha-mannosidosis was calculated at 1 in 784,535, l in 166,801 Europeans (Finnish), and l in 431,689 East Asians. By integrating the data from the 8,936 Koreans in gnomAD Korean, KOVA and KRGDB, the carrier frequency of alpha-mannosidosis in the Korean population was 0.04% and estimated incidence was 1 in 19,963,024. Conclusion: This study is the first to investigate the carrier frequencies of alpha-mannosidosis in East Asians and Koreans, including specific subpopulations, utilizing gnomAD and the Korean genomic database. The variant spectrum of MAN2B1 genes in East Asians showed significant differences compared to other ethnic groups. Our data provide valuable reference information for future investigations into alpha-mannosidosis, aiding in understanding the genetic diversity and specific variants associated with the condition in East Asian populations.

A Common Functional Variant at the Enhancer of the Rheumatoid Arthritis Risk Gene ORMDL3 Regulates its Expression Through Allele-Specific JunD Binding

Genome-wide association studies (GWASs) have identified over 100 loci associated with rheumatoid arthritis (RA); however, the functionally affected genes and the underlying molecular mechanisms contributing to these associations are often unknown. In this study, we conducted an integrative genomic analysis incorporating multiple "omics" data and identified a functional regulatory DNA variant, rs56199421, and a plausible mechanism by which it regulates the expression of a putative RA risk gene, ORMDL Sphingolipid Biosynthesis Regulator 3 (ORMDL3). The T allele of rs56199421, located in the enhancer region of ORMDL3, exhibited stronger direct binding ability than the other C allele of rs56199421 did in vitro with the transcription factor JunD and demonstrated higher transcriptional activity. Moreover, the T allele of rs56199421 is associated with elevated RA risk, and ORMDL3 expression is increased in RA patients. Thus, these findings suggest that the T allele of rs56199421 enhances JunD transcription factor binding, increases enhancer activity, and elevates the expression of the RA risk gene ORMDL3.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00107-z.

Epigenetic Alteration of H3K27me3 as a Possible Oncogenic Mechanism of Central Neurocytoma

Central neurocytoma (CN) is a low-grade neuronal tumor that mainly arises from the lateral ventricle (LV). This tumor remains poorly understood in the sense that no driver gene aberrations have been identified thus far. We investigated immunomarkers in fetal and adult brains and 45 supratentorial periventricular tumors to characterize the biomarkers, cell of origin, and tumorigenesis of CN. All CNs occurred in the LV. A minority involved the third ventricle, but none involved the fourth ventricle. As expected, next-generation sequencing performed using a brain-tumor-targeted gene panel in 7 CNs and whole exome sequencing in 5 CNs showed no driver mutations. Immunohistochemically, CNs were robustly positive for FGFR3 (100%), SSTR2 (92%), TTF-1 (Nkx2.1) (88%), GLUT-1 (84%), and L1CAM (76%), in addition to the well-known markers of CN, synaptophysin (100%) and NeuN (96%). TTF-1 was also positive in subependymal giant cell astrocytomas (100%, 5/5) and the pituicyte tumor family, including pituicytoma and spindle cell oncocytoma (100%, 5/5). Interestingly, 1 case of LV subependymoma (20%, 1/5) was positive for TTF-1, but all LV ependymomas were negative (0/5 positive). Because TTF-1-positive cells were detected in the medial ganglionic eminence around the foramen of Monro of the fetal brain and in the subventricular zone of the LV of the adult brain, CN may arise from subventricular TTF-1-positive cells undergoing neuronal differentiation. H3K27me3 loss was observed in all CNs and one case (20%) of LV subependymoma, suggesting that chromatin remodeling complexes or epigenetic alterations may be involved in the tumorigenesis of all CNs and some ST-subependymomas. Further studies are required to determine the exact tumorigenic mechanism of CN.

Comparison of Risk Allele Frequencies of Psoriasis-Associated Single-Nucleotide Polymorphisms in Different Population Groups

BACKGROUND

The prevalence of psoriasis differs by population, and it appears to be more common among Europeans than in East Asians. Recent genome-wide association studies (GWAS) have identified alleles that increase the risk of psoriasis, and these alleles may present different frequencies in different geographic regions.

OBJECTIVE

We aimed to gain insights into the causes of differences in disease frequencies according to populations and the factors affecting prevalence and pattern differences.

METHODS

We collected a total of 147 psoriasis-associated single-nucleotide polymorphisms (SNPs) from the GWAS catalog and compared the allele frequency differences in 27 populations using public population frequency in the 1000 Genomes Project phase 3 (n=2,504) and the Korean Reference Genome Database (n=1,722). Additionally, we calculated the composited genetic risk scores across the population groups.

RESULTS

There were distinct patterns of allele frequencies in different population groups. In many cases, East Asians exhibited allele frequencies opposite to that of Europeans. The genetic risk score was higher in Europeans (average: 0.487) and Americans (average: 0.492) than in East Asians (average: 0.471). The prevalence of psoriasis correlated with the average genetic risk score of the population.

CONCLUSION

We observed a difference in the allele frequencies of psoriasis-associated SNPs between the studied populations. This result suggests that the difference in the prevalence of psoriasis between population groups can be interpreted to some extent by the genotype.

Pathogenic Interleukin-10 Receptor Alpha Variants in Humans - Balancing Natural Selection and Clinical Implications

Balancing natural selection is a process by which genetic variants arise in populations that are beneficial to heterozygous carriers, but pathogenic when homozygous. We systematically investigated the prevalence, structural, and functional consequences of pathogenic IL10RA variants that are associated with monogenic inflammatory bowel disease. We identify 36 non-synonymous and non-sense variants in the IL10RA gene. Since the majority of these IL10RA variants have not been functionally characterized, we performed a systematic screening of their impact on STAT3 phosphorylation upon IL-10 stimulation. Based on the geographic accumulation of confirmed pathogenic IL10RA variants in East Asia and in Northeast China, the distribution of infectious disorders worldwide, and the functional evidence of IL-10 signaling in the pathogenesis, we identify Schistosoma japonicum infection as plausible selection pressure driving variation in IL10RA. Consistent with this is a partially augmented IL-10 response in peripheral blood mononuclear cells from heterozygous variant carriers. A parasite-driven heterozygote advantage through reduced IL-10 signaling has implications for health care utilization in regions with high allele frequencies and potentially indicates pathogen eradication strategies that target IL-10 signaling.

Comparison of histological and molecular features of pediatric-type follicular lymphoma and pediatric nodal marginal zone lymphoma

Pediatric-type follicular lymphoma and pediatric nodal marginal zone lymphoma are pediatric B cell lymphomas with similar clinical characteristics but distinct histological features. We investigated the differences between pediatric-type follicular lymphoma and pediatric nodal marginal zone lymphoma by comparing their histological and molecular characteristics. A total of 5 pediatric-type follicular lymphoma and 11 pediatric nodal marginal zone lymphoma patients were included in the study. In the histological review, 5 of the 16 cases showed overlapping morphological features of pediatric-type follicular lymphoma and pediatric nodal marginal zone lymphoma; hence, they were reclassified as "mixed type." In molecular analysis, using panel-based massively parallel sequencing, MAP2K1, TNFRSF14, and IRF8 mutations were found in 6, 3, and 2 of the 11 pediatric nodal marginal zone lymphoma patients, respectively, and IRF8 mutation was found in one of the five pediatric-type follicular lymphoma patients. There were no significant differences in genetic alterations established from the histologically reclassified diagnosis as well as the initial diagnosis. Pediatric-type follicular lymphoma and pediatric nodal marginal zone lymphoma showed morphological overlap in some cases, and no difference between the two was found upon molecular analysis. These findings suggest the possibility that pediatric-type follicular lymphoma and pediatric nodal marginal zone lymphoma are single entity pediatric B-cell lymphoma with broad morphological spectrum.

Carrier frequency and incidence estimation of RPE65-associated inherited retinal diseases in East Asian population by population database-based analysis

BACKGROUND

Inherited retinal diseases (IRDs) are clinically and genetically heterogenous disorders leading to visual impairment and blindness. Because gene therapy for RPE65-associated IRDs was recently approved, it is necessary to predict the carrier frequency and prevalence for RPE65-associated IRDs. This study aimed to analyze the carrier frequency and expected incidence of RPE65-associated IRDs in East Asians and Koreans using exome data from the Genome Aggregation Database (gnomAD) and the Korean Reference Genome Database (KRGDB).

METHODS

We analyzed 9,197 exomes for East Asian populations from gnomAD comprising 1,909 Korean and 1,722 Korean genomes from KRGDB. All identified RPE65 variants were classified according to the 2015 American College of Medical Genetics and Genomics and the Association for Molecular Pathology guidelines.

RESULTS

The total carrier frequencies of East Asians and Koreans from both gnomAD and KRGDB were 0.10% (11/10,919) and 0.06% (2/3,631), respectively. The estimated incidence of RPE65-associated IRDs was 1/3,941,308 in East Asians and 1/13,184,161 in Koreans.

CONCLUSION

This study identified carrier frequencies of RPE65-associated IRDs in East Asians and Koreans using gnomAD and KRGDB. We confirmed that the carrier frequency of RPE65-associated IRDs patients was low in Koreans among all East Asian populations, and the incidence was also predicted to be lower than in other East Asian populations. The variant spectrum of RPE65 gene in East Asian and Korean populations differed greatly from those of other ethnic groups.

The first broad replication study of SNPs and a pilot genome-wide association study for androgenetic alopecia in Asian populations

BACKGROUND

Many candidate genes for androgenetic alopecia (AGA) have been identified in studies of the Caucasians and some Asian populations.

AIMS

This study aimed to confirm the known susceptibility genes reported in previous studies and find additional candidate genes for high-risk individuals for AGA in Korean population.

PATIENTS/METHODS

We recapitulated the previously reported SNPs and identified the novel Korean AGA risk genetic variants using a Korean hospital-based AGA case and control samples. The population was consisting of 494 individuals (275 AGA cases and 146 controls). Using the 800 K SNPs of precision medical research array (PMRA SNP microarray chip) and imputation-based SNPs, 12 previous GWAS reports for AGA and a total of 62 160 SNPs were examined in our study samples. Also, we conducted the genome-wide association study (GWAS) by the logistic regression analyses for AGA cases and controls with controlling the age as the covariates.

RESULTS

Among the 62 160 SNPs, a total of 1143 SNPs in 76 gene regions showed weak replication tendency with the p-values <0.05 and same direction of effects. Additionally, the GWAS results showed 110 SNPs in 13 independent regions with the suggestive p-values <1.00 × 10^-5 . The most significantly replicated SNP resided on chromosome 20, which were similar to other AGA replication studies including Chinese study. The GWAS identified two SNPs (rs11010734 and rs2420640) increasing the risk for AGA in our study population.

CONCLUSIONS

Our study would be a reference of the non-European studies to better understand AGA in different populations and ancestral contexts.

A single arm phase Ib/II trial of first-line pembrolizumab, trastuzumab and chemotherapy for advanced HER2-positive gastric cancer

In this multi-center phase II trial, we evaluated the efficacy and safety of a quadruplet regimen (pembrolizumab, trastuzumab, and doublet chemotherapy) as first-line therapy for unresectable or metastatic human epidermal growth factor receptor 2 (HER2)-positive advanced gastric cancer (AGC) (NCT02901301). The primary endpoints were recommended phase 2 dose (RP2D) for phase Ib and objective response rate (ORR) for phase II. The secondary endpoints included progression-free survival (PFS), overall survival (OS), duration of response, time to response and safety. Without dose-limiting or unexpected toxicities, the starting dose in the phase Ib trial was selected as RP2D. In 43 patients, the primary endpoint was achieved: the objective response rate was 76.7% (95% confidence interval [CI]: 61.4-88.2), with complete and partial responses in 14% and 62.8% of patients, respectively. The median progression-free survival, overall survival, and duration of response were 8.6 months, 19.3 months, and 10.8 months, respectively. No patients discontinued pembrolizumab because of immune-related adverse events. Programmed death ligand-1 status was not related to survival. Post hoc analyses of pretreatment tumor specimens via targeted sequencing indicated that ERBB2 amplification, RTK/RAS pathway alterations, and high neoantigen load corrected by HLA-B were positively related to survival. The current quadruplet regimen shows durable efficacy and safety for patients with HER2-positive AGC.

Genetic Screening in Korean Patients with Frontotemporal Dementia Syndrome

Background: Frontotemporal dementia (FTD) syndrome is a genetically heterogeneous group of diseases. However, pathogenic variants in the chromosome 9 open reading frame 72 (C9orf72), microtubule-associated protein tau (MAPT), and progranulin (GRN) genes are mainly associated with genetic FTD in Caucasian populations. Objective: To understand the genetic background of Korean patients with FTD syndrome. Methods: We searched for pathogenic variants of 52 genes related to FTD, amyotrophic lateral sclerosis, familial Alzheimer’s disease, and other dementias, and hexanucleotide repeats of the C9orf72 gene in 72 Korean patients with FTD using whole exome sequencing and the repeat-primed polymerase chain reaction, respectively. Results: One likely pathogenic variant, p.G706R of MAPT, in a patient with behavioral variant FTD (bvFTD) and 13 variants of uncertain significance (VUSs) in nine patients with FTD were identified. Of these VUSs, M232R of the PRNP gene, whose role in pathogenicity is controversial, was also found in two patients with bvFTD. Conclusions: These results indicate that known pathogenic variants of the three main FTD genes (MAPT, GRN, and C9orf72) in Western countries are rare in Korean FTD patients.

Network-based meta-analysis and the candidate gene association studies reveal novel ethnicity-specific variants in MFSD3 and MRPL43 associated with dementia with Lewy bodies

Dementia with Lewy bodies (DLB) is the second most common form of neurodegenerative dementia in elderly people, following Alzheimer's disease. Only three genes, SNCA (α-synuclein), APOE (apolipoprotein E), and GBA (glucosylceramidase), have been convincingly demonstrated to be associated with DLB. Here, we applied whole-genome sequencing to blood samples from 61 DLB patients and 45 cognitively normal controls. We used accumulation of candidate mutations to detect novel DLB-associated genes. Subsequent single nucleotide polymorphism (SNP) genotyping and association studies in a large number of samples from Japanese individuals revealed novel heterozygous variants in MFSD3 (rs143475431, c.888T>A:p.C296*; n = 5,421, p = 0.00063) and MRPL43 (chr10:102746730, c.241A>C:p.N81H; n = 4,782, p = 0.0029). We further found that the MFSD3 variant increased plasma levels of butyrylcholinesterase (n = 1,206, p = 0.029). We believe that our findings will contribute to the understanding of DLB and provide insight into its pathogenic mechanism for future studies.

CD47 overexpression is common in intestinal non-GCB type diffuse large B-cell lymphoma and associated with 18q21 gain

The CD47/signal regulatory protein α pathway is an emerging immune checkpoint that is a new therapeutic target. We investigated CD47 expression in diffuse large B-cell lymphoma (DLBCL) of various subtypes and organs. Moreover, the relationship between CD47 expression and genetic alterations was analyzed using panel-based massively parallel sequencing (next-generation sequencing [NGS]). CD8, CD68, and CD47 immunohistochemical staining were performed on 238 patients with DLBCL. CD47 was scored according to intensity on a 5-level scale, and CD8 and CD68 were quantitatively evaluated using QuPath software. Panel-based NGS was performed in 37 patients. In CD8 and CD68 quantitative analyses by organs, intestinal DLBCL showed significantly lower cytotoxic T-cell infiltration than that in others (P < .001). The CD47-high group comprised 24 of 58 (41.4%) patients in the group with DLBCL from intestine and 15 of 180 (8.3%) patients in the group with DLBCL from other organs (P < .001). The 18q21 gain/amplification was found in 10 of 37 patients, and all of them were CD47-high. Intestinal CD47-high DLBCL occurred in terminal ileum to ascending colon and was restricted to nongerminal center B-cell type. In the survival analyses, the prognosis of nonintestinal CD47-high DLBCL was poorer than that of intestinal CD47-high DLBCL (P = .025). CD47-high DLBCL was closely associated with 18q21 gain/amplification and showed a high prevalence in intestine. We propose to classify CD47-high DLBCL into intestinal and nonintestinal types. Further studies are necessary to assess whether the constellation of features seen here is reproducible and sufficient to consider primary intestinal DLBCL as a distinct biological entity.

The First Korean Siblings With Adult-Onset 4H Leukodystrophy Related to Nonsynonymous POLR3B Mutations

Objectives

4H leukodystrophy is a rare autosomal recessive hypomyelinating disorder characterized by several combinations of motor dysfunction, abnormal dentition, and ophthalmic and endocrine abnormalities. To date, only a single Korean case report of pediatric leukodystrophy caused by the POLR1C sequence variation has been published, while there are no reports on the POLR3B, POLR3A, or POLR3K variants.

Methods

Genetic tests of Korean sibling pairs with primary amenorrhea due to normosmic isolated hypogonadotropic hypogonadism and cognitive or behavioral symptoms were performed by whole-exome sequencing (WES). The WES results were validated by direct Sanger sequencing.

Results

We identified biallelic variations in the POLR3B gene of p.Tyr685* and p.Tyr746Cys, which have not been associated with 4H leukodystrophy. Both sequence variants lie in the hybrid-binding domain of the protein RPC2. The protein structure analysis predicted that cysteine substitution of the phylogenetically conserved amino acid tyrosine can cause destabilization.

Discussion

The siblings reported are the first POLR3B-related hypomyelinating leukodystrophy cases in Korea. Our report expands the mutational spectrum of 4H leukodystrophy and suggests that it is mandatory to consider its diagnostic possibility in adult patients presenting with primary amenorrhea and mild cognitive or behavioral symptoms.

Pathogenic NOTCH3 Variants Are Frequent Among the Korean General Population

OBJECTIVE

This study aimed to determine the frequency of pathogenic NOTCH3 variants among Koreans.

METHODS

In this cross-sectional study, we queried for pathogenic NOTCH3 variants in 2 Korean public genome databases: the Korean Reference Genome Database (KRGDB) and the Korean Genome Project (Korea1K). In addition, we screened the 3 most common pathogenic NOTCH3 variants (p.Arg75Pro, p.Arg544Cys, and p.Arg578Cys) for 1,000 individuals on Jeju Island, where the largest number of patients with cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) have been reported in Korea.

RESULTS

The pathogenic NOTCH3 variant (p.Arg544Cys) was found in 0.12% of sequences in the KRGDB, and 3 pathogenic variants (p.Arg75Pro, p.Arg182Cys, and p.Arg544Cys) were present in 0.44% of the Korea1K database. Of the 1,000 individuals on Jeju Island, we found 2 cysteine-altering NOTCH3 variants (p.Arg544Cys variant in 9 and p.Arg578Cys in 1 individual) in 1.00% of the participants (95% confidence interval: 0.48%-1.83%). The presence of cysteine-altering NOTCH3 variants was significantly associated with a history of stroke (p < 0.001).

DISCUSSION

Pathogenic NOTCH3 variants are frequently found in the general Korean population. Such a high prevalence of pathogenic variants could threaten the brain health of tens of thousands to hundreds of thousands of older adults in Korea.

Real-world data analysis of patients with cancer of unknown primary

Cancer of unknown primary (CUP) is a heterogeneous malignancy in which the primary site of the tumor cannot be identified through standard work-up. The survival outcome of CUP is generally poor, and there is no consensus for treatment. Here, we comprehensively analyzed the real-world data of 218 patients with CUP (median age, 62 years [range, 19-91]; male, 62.3%). Next-generation sequencing was conducted in 22 (10%) patients, one of whom showed level 1 genetic alteration. Most (60.3%) patients were treated with empirical cytotoxic chemotherapy, and two patients received targeted therapy based on the NGS results. The median OS was 8.3 months (95% confidence interval [CI] 6.2-11.4), and the median progression-free survival of patients treated with chemotherapy was 4.4 months (95% CI 3.4-5.3). In multivariate Cox regression analysis, Eastern Cooperative Oncology Group performance status (ECOG PS) of 0 or 1 and localized disease were significantly associated with favorable survival outcomes. Collectively, we found that CUP patients had a poor prognosis after standard treatment, and those with localized disease who received local treatment and those with better PS treated with multiple lines of chemotherapy had better survival outcomes. Targeted therapies based on NGS results are expected to improve survival outcomes.

Genetic profile of primary plasma cell leukemia in Korea: comparison with plasma cell myeloma

Plasma cell leukemia (PCL) is clinically and genetically distinct from multiple myeloma (MM), despite controversies regarding the disease definition. To determine the distinct features of PCL, the genetic property of primary PCL (pPCL) was compared with that of secondary PCL (sPCL) and MM. In patients with pPCL, Eighty-nine non-synonymous mutations were observed in 68 genes. The most frequently mutated genes were TP53, TSC2, and TYK2. In comparison with genetic abnormalities of sPCL and MM, 45 genes were present only in pPCL while 28 genes were only in sPCL and 22 genes only in MM. Among the common genes between pPCL and MM, a higher prevalence of TP53 was observed in pPCL, compared to MM (p < 0.05), while similar, compared to sPCL (p = 0.64). In summary, pPCL patients showed a higher level of genetic heterogeneity and distinctive genetic signature in their mutational profile compared to patients with MM and sPCL.

Comparison of tumor mutation burden of 300 various non-Hodgkin lymphomas using panel based massively parallel sequencing

BACKGROUND

Tumor mutation burden is an emerging biomarker for immunotherapy. Although several clinical trials for immunotherapy in lymphoma have been carried out, the mutation burden of various lymphomas is not well known yet. Thus, the objective of this study was to compare tumor mutation burden of various non-Hodgkin lymphomas using panel based massively parallel sequencing.

METHODS

We conducted 405 gene panel based massively parallel sequencing of 300 non-Hodgkin lymphomas and investigate the number of SNV/Indel in each lymphoma.

RESULTS

The number of SNV/Indel was higher in mature B-cell lymphoma than in mature T- and NK-cell lymphoma. (P < 0.001) The number of SNV/Indel in primary mediastinal large B-cell lymphoma and primary diffuse large B-cell lymphoma of the central nervous system was the highest, which was significantly higher than that in diffuse large B-cell lymphoma, not otherwise specified (DLBCL NOS).(P = 0.030 and P = 0.008, respectively) The SNV/Indel number in EBV-positive DLBCL NOS was significantly lower than that in DLBCL NOS. (P = 0.048) Peripheral T-cell lymphoma, NOS showed no significant difference in the number of SNV/Indel from extranodal NK/T-cell lymphoma, nasal type (P = 0.942) or angioimmunoblastic T-cell lymphoma (P = 0.739). The number of SNV/Indel in anaplastic large cell lymphoma, ALK-positive was significantly lower than that in anaplastic large cell lymphoma, ALK-negative (P = 0.049). It was the lowest among all the lymphomas considered.

CONCLUSION

Various lymphomas have different mutation burdens. Thus, tumor mutation burden can be used as a promising biomarker for immunotherapy in lymphomas.

Genetic Polymorphisms of Pesticide-Metabolizing Enzymes and Transporters in Agricultural Workers and Thyroid Hormone Levels

PURPOSE

Chronic exposure to pesticides has been associated with thyroid dysfunction owing to their endocrine disruption ability. Genetic variations in genes encoding phase I and II enzymes and phase III transporters are partly responsible for individual responses to chemical pesticides. This study investigated the association between variations in genes involved in pesticide metabolism and altered thyroid hormone concentrations.

METHODS

We assessed thyroid-stimulating hormone (TSH), free triiodothyronine (FT3), and free thyroxine (FT4) in organic agriculture workers (n = 216) and workers who used chemical pesticides (n = 229). A questionnaire was used to collect sociodemographic, pesticide exposure, and health status data. Blood samples were analyzed for TSH, FT3, and FT4. Genomic DNA was extracted and genotyped using the TaqMan real-time PCR genotyping assay and restriction fragment length polymorphism method for 15 metabolically related genes.

RESULTS

Significant differences in the TSH (1.58 vs 1.12 µIU/mL) and FT3 (0.34 vs 0.31 ng/dL) concentrations between the chemical and organic worker groups were observed. The frequencies of all single nucleotide polymorphisms were in Hardy-Weinberg equilibrium and were mostly consistent with Asian populations. The findings showed the association between SNPs of enzymes and transporters and TSH, FT3, and FT4. The odd ratio and adjusted odd ratio (with sex, age, smoking status, alcohol consumption and exposure parameters) for subclinical thyroid disease by the variant alleles CYP1A1 rs1048943, CYP2B6 rs2279343, CYP2C19 rs4244285, NAT2 rs1799931, and PON1 rs662 in the chemical workers compared with the organic workers were found (P values < 0.05).

CONCLUSION

This is the first study to assess gene-environment interactions in Thai agricultural workers by investigating disruptions of the hypothalamic-pituitary-thyroid axis. The investigated SNP profiles revealed several gene-thyroid hormone associations in which even low levels of pesticide exposure could disturb thyroid homeostasis. These findings provide a foundation for planning future studies investigating associations between complex diseases and occupational pesticide exposure.

New insights into hallux valgus by whole exome sequencing study

The traditional view is that the occurrence and development of hallux valgus (HV) are mainly due to environmental factors. Recent studies have suggested the large contribution of genetic heritability to HV, but it remains elusive about the genetic variants underlying the development of HV. To gain knowledge about the molecular mechanisms of HV pathogenesis by genetic approach, whole exome sequencing studies were performed in 10 individuals (7 affected by HV and 3 unaffected) from three independent families. Specific mutations were found to be related to the pathogenesis of HV and conform to the laws of inheritance. A total of 36 genes with functional candidate single nucleotide variants were identified. Genetic predisposition plays an important role in the development of HV. Interestingly, some of these genes are related to chronic arthritis, such as the complement encoding gene C7, or are related to long toe or long fingers, such as TTN, COL6A3, LARS, FIG4, and CBS. This study identified rare potentially pathogenic mutations represented by genes related to digital anomalies and chronic arthritis underlying the familial types of HV, which acquired new insights into the genetic and physiological foundations of HV, thereby might improve accurate prevention and drug development for HV.

KCNV2-Associated Retinopathy: Genetics, Electrophysiology, and Clinical Course-KCNV2 Study Group Report 1

PURPOSE

To investigate genetics, electrophysiology, and clinical course of KCNV2-associated retinopathy in a cohort of children and adults.

STUDY DESIGN

This was a multicenter international clinical cohort study.

METHODS

Review of clinical notes and molecular genetic testing. Full-field electroretinography (ERG) recordings, incorporating the international standards, were reviewed and quantified and compared with age and recordings from control subjects.

RESULTS

In total, 230 disease-associated alleles were identified from 117 patients, corresponding to 75 different KCNV2 variants, with 28 being novel. The mean age of onset was 3.9 years old. All patients were symptomatic before 12 years of age (range, 0-11 years). Decreased visual acuity was present in all patients, and 4 other symptoms were common: reduced color vision (78.6%), photophobia (53.5%), nyctalopia (43.6%), and nystagmus (38.6%). After a mean follow-up of 8.4 years, the mean best-corrected visual acuity (BCVA ± SD) decreased from 0.81 ± 0.27 to 0.90 ± 0.31 logarithm of minimal angle of resolution. Full-field ERGs showed pathognomonic waveform features. Quantitative assessment revealed a wide range of ERG amplitudes and peak times, with a mean rate of age-associated reduction indistinguishable from the control group. Mean amplitude reductions for the dark-adapted 0.01 ERG, dark-adapted 10 ERG a-wave, and LA 3.0 30 Hz and LA3 ERG b-waves were 55%, 21%, 48%, and 74%, respectively compared with control values. Peak times showed stability across 6 decades.

CONCLUSION

In KCNV2-associated retinopathy, full-field ERGs are diagnostic and consistent with largely stable peripheral retinal dysfunction. Report 1 highlights the severity of the clinical phenotype and established a large cohort of patients, emphasizing the unmet need for trials of novel therapeutics.

Analysis of BRCA1/2 variants of unknown significance in the prospective Korean Hereditary Breast Cancer study

Genetic testing for BRCA1 and BRCA2 is crucial in diagnosing hereditary breast and ovarian cancer syndromes and has increased with the development of multigene panel tests. However, results classified as variants of uncertain significance (VUS) present challenges to clinicians in attempting to choose an appropriate management plans. We reviewed a total of 676 breast cancer patients included in the Korean Hereditary Breast Cancer (KOHBRA) study with a VUS on BRCA mutation tests between November 2007 and April 2013. These results were compared to the ClinVar database. We calculated the incidence and odds ratios for these variants using the Korean Reference Genome Database. A total of 58 and 91 distinct VUS in BRCA1 and BRCA2 were identified in the KOHBRA study (comprising 278 and 453 patients, respectively). A total of 27 variants in the KOHBRA study were not registered in the Single Nucleotide Polymorphism database. Among BRCA1 VUSs, 20 were reclassified as benign or likely benign, four were reclassified as pathogenic or likely pathogenic, and eight remained as VUSs according to the ClinVar database. Of the BRCA2 VUSs, 25 were reclassified as benign or likely benign, two were reclassified as pathogenic or likely pathogenic, and 33 remained as VUS according to the ClinVar database. There were 12 variants with conflicting interpretations of pathogenicity for BRCA1 and 18 for BRCA2. Among them, p.Leu1780Pro showed a particularly high odds ratio. Six pathogenic variants and one conflicting variant identified using ClinVar could be reclassified as pathogenic variants in this study. Using updated ClinVar information and calculating odds ratios can be helpful when reclassifying VUSs in BRCA1/2.

Genetic variation analyses indicate conserved SARS-CoV-2-host interaction and varied genetic adaptation in immune response factors in modern human evolution

Coronavirus disease 2019 (COVID‐19), caused by severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2), is a pandemic as of early 2020. Upon infection, SARS‐CoV‐2 attaches to its receptor, that is, angiotensin‐converting enzyme 2 (ACE2), on the surface of host cells and is then internalized into host cells via enzymatic machineries. This subsequently stimulates immune response factors. Since the host immune response and severity of COVID‐19 vary among individuals, genetic risk factors for severe COVID‐19 cases have been investigated. Our research group recently conducted a survey of genetic variants among SARS‐CoV‐2‐interacting molecules across populations, noting near absence of difference in allele frequency spectrum between populations in these genes. Recent genome‐wide association studies have identified genetic risk factors for severe COVID‐19 cases in a segment of chromosome 3 that involves six genes encoding three immune‐regulatory chemokine receptors and another three molecules. The risk haplotype seemed to be inherited from Neanderthals, suggesting genetic adaptation against pathogens in modern human evolution. Therefore, SARS‐CoV‐2 uses highly conserved molecules as its virion interaction, whereas its immune response appears to be genetically biased in individuals to some extent. We herein review the molecular process of SARS‐CoV‐2 infection as well as our further survey of genetic variants of its related immune effectors. We also discuss aspects of modern human evolution.

Exome Chip Analysis of 14,026 Koreans Reveals Known and Newly Discovered Genetic Loci Associated with Type 2 Diabetes Mellitus

BACKGROUND

Most loci associated with type 2 diabetes mellitus (T2DM) discovered to date are within noncoding regions of unknown functional significance. By contrast, exonic regions have advantages for biological interpretation.

METHODS

We analyzed the association of exome array data from 14,026 Koreans to identify susceptible exonic loci for T2DM. We used genotype information of 50,543 variants using the Illumina exome array platform.

RESULTS

In total, 7 loci were significant with a Bonferroni adjusted P=1.03×10-6. rs2233580 in paired box gene 4 (PAX4) showed the highest odds ratio of 1.48 (P=1.60×10-10). rs11960799 in membrane associated ring-CH-type finger 3 (MARCH3) and rs75680863 in transcobalamin 2 (TCN2) were newly identified loci. When we built a model to predict the incidence of diabetes with the 7 loci and clinical variables, area under the curve (AUC) of the model improved significantly (AUC=0.72, P<0.05), but marginally in its magnitude, compared with the model using clinical variables (AUC=0.71, P<0.05). When we divided the entire population into three groups-normal body mass index (BMI; <25 kg/m2), overweight (25≤ BMI <30 kg/m2), and obese (BMI ≥30 kg/m2) individuals-the predictive performance of the 7 loci was greatest in the group of obese individuals, where the net reclassification improvement was highly significant (0.51; P=8.00×10-5).

CONCLUSION

We found exonic loci having a susceptibility for T2DM. We found that such genetic information is advantageous for predicting T2DM in a subgroup of obese individuals.

A survey of genetic variants in SARS-CoV-2 interacting domains of ACE2, TMPRSS2 and TLR3/7/8 across populations

The COVID-19 pandemic highlighted healthcare disparities in multiple countries. As such morbidity and mortality vary significantly around the globe between populations and ethnic groups. Underlying medical conditions and environmental factors contribute higher incidence in some populations and a genetic predisposition may play a role for severe cases with respiratory failure. Here we investigated whether genetic variation in the key genes for viral entry to host cells-ACE2 and TMPRSS2-and sensing of viral genomic RNAs (i.e., TLR3/7/8) could explain the variation in incidence across diverse ethnic groups. Overall, these genes are under strong selection pressure and have very few nonsynonymous variants in all populations. Genetic determinant for the binding affinity between SARS-CoV-2 and ACE2 does not show significant difference between populations. Non-genetic factors are likely to contribute differential population characteristics affected by COVID-19. Nonetheless, a systematic mutagenesis study on the receptor binding domain of ACE2 is required to understand the difference in host-viral interaction across populations.

No association between SCN9A and monogenic human epilepsy disorders

Many studies have demonstrated the clinical utility and importance of epilepsy gene panel testing to confirm the specific aetiology of disease, enable appropriate therapeutic interventions, and inform accurate family counselling. Previously, SCN9A gene variants, in particular a c.1921A>T p.(Asn641Tyr) substitution, have been identified as a likely autosomal dominant cause of febrile seizures/febrile seizures plus and other monogenic seizure phenotypes indistinguishable from those associated with SCN1A, leading to inclusion of SCN9A on epilepsy gene testing panels. Here we present serendipitous findings of genetic studies that identify the SCN9A c.1921A>T p.(Asn641Tyr) variant at high frequency in the Amish community in the absence of such seizure phenotypes. Together with findings in UK Biobank these data refute an association of SCN9A with epilepsy, which has important clinical diagnostic implications.