Korean Genome Project: 1094 Korean personal genomes with clinical information

Germline functional variants contribute to somatic mutation and outcomes in neuroblastoma

Germline genetic context may play a significant role in the development and evolution of cancer, particularly in childhood cancers such as neuroblastoma. This study investigates the role of putatively functional germline variants in neuroblastoma, even if they do not directly increase disease risk. Our whole-exome sequencing analysis of 125 patients with neuroblastoma reveals a positive correlation between germline variant burden and somatic mutations. Moreover, patients with higher germline variant burden exhibit worse outcomes. Similar findings are observed in the independent neuroblastoma cohort where a higher germline variant burden correlates with a higher somatic mutational burden and a worse overall survival outcome. However, contrasting results emerge in adult-onset cancer, emphasizing the importance of germline genetics in neuroblastoma. The enrichment of putatively functional germline variants in cancer predisposition genes is borderline significant when compared to healthy populations (P = 0.077; Odds Ratio, 1.45; 95% confidence intervals, 0.94-2.21) and significantly more pronounced against adult-onset cancers (P = 0.016; Odds Ratio, 2.13; 95% confidence intervals, 1.10-3.91). Additionally, the presence of these variants proves to have prognostic significance in neuroblastoma (log-rank P < 0.001), and combining germline with clinical risk factors notably improves survival predictions.

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.

Quantification method of ctDNA using cell-free DNA methylation profile for noninvasive screening and monitoring of colon cancer

BACKGROUND

Colon cancer ranks as the second most lethal form of cancer globally. In recent years, there has been active investigation into using the methylation profile of circulating tumor DNA (ctDNA), derived from blood, as a promising indicator for diagnosing and monitoring colon cancer.

RESULTS

We propose a liquid biopsy-based epigenetic method developed by utilizing 49 patients and 260 healthy controls methylation profile data to screen and monitor colon cancer. Our method initially identified 901 colon cancer-specific hypermethylated (CaSH) regions in the tissues of the 49 cancer patients. We then used these CaSH regions to accurately quantify the amount of circulating tumor DNA (ctDNA) in the blood samples of these same patients, utilizing cell-free DNA methylation profiles. Notably, the methylation profiles of ctDNA in the blood exhibited high sensitivity (82%) and specificity (93%) in distinguishing patients with colon cancer from the control group, with an area under the curve of 0.903. Furthermore, we confirm that our method for ctDNA quantification is effective for monitoring cancer patients and can serve as a valuable tool for postoperative prognosis.

CONCLUSIONS

This study demonstrated a successful application of the quantification of ctDNA among cfDNA using the original cancer tissue-derived CaSH region for screening and monitoring colon cancer.

Brief Report: Risk Variants Could Inform Early Neurodevelopmental Outcome in Children with Developmental Disabilities

The aim of this study was to examine genetic variations underlying the early neurodevelopmental outcome of developmental disabilities (DDs). The study cohort consisted of 191 children with DDs. Diagnosis was assessed at baseline and at the index age (72-84 months). Exome sequencing was conducted and putative risk variants were identified. According to the diagnostic improvement, children were categorized into the improvement group (n = 19) and the non-improvement group (n = 172). Compared to the non-improvement group, the improvement group had lower number of risk variants in known DD genes and haploinsufficient genes, and lower number of overall putative risk variants. Our results may serve as a preliminary basis for developing a model that informs clinical outcome by sequencing analysis.

Comprehensive multiomics analysis reveals distinct differences between pediatric choroid plexus papilloma and carcinoma

Choroid plexus tumors (CPTs) are intraventricular tumors derived from the choroid plexus epithelium and occur frequently in children. The aim of this study was to investigate the genomic and epigenomic characteristics of CPT and identify the differences between choroid plexus papilloma (CPP) and choroid plexus carcinoma (CPC). We conducted multiomics analyses of 20 CPT patients including CPP and CPC. Multiomics analysis included whole-genome sequencing, whole-transcriptome sequencing, and methylation sequencing. Mutually exclusive TP53 and EPHA7 point mutations, coupled with the amplification of chromosome 1, were exclusively identified in CPC. In contrast, amplification of chromosome 9 was specific to CPP. Differential gene expression analysis uncovered a significant overexpression of genes related to cell cycle regulation and epithelial-mesenchymal transition pathways in CPC compared to CPP. Overexpression of genes associated with tumor metastasis and progression was observed in the CPC subgroup with leptomeningeal dissemination. Furthermore, methylation profiling unveiled hypomethylation in major repeat regions, including long interspersed nuclear elements, short interspersed nuclear elements, long terminal repeats, and retrotransposons in CPC compared to CPP, implying that the loss of epigenetic silencing of transposable elements may play a role in tumorigenesis of CPC. Finally, the differential expression of AK1, regulated by both genomic and epigenomic factors, emerged as a potential contributing factor to the histological difference of CPP against CPC. Our results suggest pronounced genomic and epigenomic disparities between CPP and CPC, providing insights into the pathogenesis of CPT at the molecular level.

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.

Human pangenome analysis of sequences missing from the reference genome reveals their widespread evolutionary, phenotypic, and functional roles

Nonreference sequences (NRSs) are DNA sequences present in global populations but absent in the current human reference genome. However, the extent and functional significance of NRSs in the human genomes and populations remains unclear. Here, we de novo assembled 539 genomes from five genetically divergent human populations using long-read sequencing technology, resulting in the identification of 5.1 million NRSs. These were merged into 45284 unique NRSs, with 29.7% being novel discoveries. Among these NRSs, 38.7% were common across the five populations, and 35.6% were population specific. The use of a graph-based pangenome approach allowed for the detection of 565 transcript expression quantitative trait loci on NRSs, with 426 of these being novel findings. Moreover, 26 NRS candidates displayed evidence of adaptive selection within human populations. Genes situated in close proximity to or intersecting with these candidates may be associated with metabolism and type 2 diabetes. Genome-wide association studies revealed 14 NRSs to be significantly associated with eight phenotypes. Additionally, 154 NRSs were found to be in strong linkage disequilibrium with 258 phenotype-associated SNPs in the GWAS catalogue. Our work expands the understanding of human NRSs and provides novel insights into their functions, facilitating evolutionary and biomedical researches.

Synergistic toxicity with copper contributes to NAT2-associated isoniazid toxicity

Anti-tuberculosis (AT) medications, including isoniazid (INH), can cause drug-induced liver injury (DILI), but the underlying mechanism remains unclear. In this study, we aimed to identify genetic factors that may increase the susceptibility of individuals to AT-DILI and to examine genetic interactions that may lead to isoniazid (INH)-induced hepatotoxicity. We performed a targeted sequencing analysis of 380 pharmacogenes in a discovery cohort of 112 patients (35 AT-DILI patients and 77 controls) receiving AT treatment for active tuberculosis. Pharmacogenome-wide association analysis was also conducted using 1048 population controls (Korea1K). NAT2 and ATP7B genotypes were analyzed in a replication cohort of 165 patients (37 AT-DILI patients and 128 controls) to validate the effects of both risk genotypes. NAT2 ultraslow acetylators (UAs) were found to have a greater risk of AT-DILI than other genotypes (odds ratio [OR] 5.6 [95% confidence interval; 2.5-13.2], P = 7.2 × 10-6). The presence of ATP7B gene 832R/R homozygosity (rs1061472) was found to co-occur with NAT2 UA in AT-DILI patients (P = 0.017) and to amplify the risk in NAT2 UA (OR 32.5 [4.5-1423], P = 7.5 × 10-6). In vitro experiments using human liver-derived cell lines (HepG2 and SNU387 cells) revealed toxic synergism between INH and Cu, which were strongly augmented in cells with defective NAT2 and ATP7B activity, leading to increased mitochondrial reactive oxygen species generation, mitochondrial dysfunction, DNA damage, and apoptosis. These findings link the co-occurrence of ATP7B and NAT2 genotypes to the risk of INH-induced hepatotoxicity, providing novel mechanistic insight into individual AT-DILI susceptibility. Yoon et al. showed that individuals who carry NAT2 UAs and ATP7B 832R/R genotypes are at increased risk of developing isoniazid hepatotoxicity, primarily due to the increased synergistic toxicity between isoniazid and copper, which exacerbates mitochondrial dysfunction-related apoptosis.

RNF213 in moyamoya disease: Genotype-phenotype association and the underlying mechanism

ABSTRACT

Moyamoya disease (MMD) is a cerebrovascular disorder characterized by a steno-occlusive internal carotid artery and compensatory vascular network formation. Although the precise pathogenic mechanism remains elusive, genetic association studies have identified RNF213 as the principal susceptibility gene for MMD, with the single nucleotide polymorphism p.R4810K recognized as the founder variant predominantly in the Asian populations. Distinct genotype-phenotype correlations are observable in RNF213-related MMD. The clinical manifestations linked to p.R4810K bear commonalities within Asian cohort, including familial predisposition, earlier age of onset, ischemic episodes, and involvement of the posterior cerebral artery (PCA). However, despite these shared phenotypic characteristics, there is significant heterogeneity in RNF213-related MMD presentations. This diversity manifests as variations across ethnic groups, inconsistent clinical symptoms and prognosis, and occurrence of other vasculopathies involving RNF213. This heterogeneity, in conjunction with the observed low disease penetrance of RNF213 mutations, suggests that the presence of these mutations may not be sufficient to cause MMD, underscoring the potential influence of other genetic or environmental factors. Although the current research might not have fully identified these additional contributors, experimental evidence points toward the involvement of RNF213 in angiogenesis, lipid metabolism, and the immune response. Future research is required to unveil the molecular mechanisms and identify the factors that synergize with RNF213 in the pathogenesis of MMD.

Korea4K: whole genome sequences of 4,157 Koreans with 107 phenotypes derived from extensive health check-ups

BACKGROUND

Phenome-wide association studies (PheWASs) have been conducted on Asian populations, including Koreans, but many were based on chip or exome genotyping data. Such studies have limitations regarding whole genome-wide association analysis, making it crucial to have genome-to-phenome association information with the largest possible whole genome and matched phenome data to conduct further population-genome studies and develop health care services based on population genomics.

RESULTS

Here, we present 4,157 whole genome sequences (Korea4K) coupled with 107 health check-up parameters as the largest genomic resource of the Korean Genome Project. It encompasses most of the variants with allele frequency >0.001 in Koreans, indicating that it sufficiently covered most of the common and rare genetic variants with commonly measured phenotypes for Koreans. Korea4K provides 45,537,252 variants, and half of them were not present in Korea1K (1,094 samples). We also identified 1,356 new genotype-phenotype associations that were not found by the Korea1K dataset. Phenomics analyses further revealed 24 significant genetic correlations, 14 pleiotropic associations, and 127 causal relationships based on Mendelian randomization among 37 traits. In addition, the Korea4K imputation reference panel, the largest Korean variants reference to date, showed a superior imputation performance to Korea1K across all allele frequency categories.

CONCLUSIONS

Collectively, Korea4K provides not only the largest Korean genome data but also corresponding health check-up parameters and novel genome-phenome associations. The large-scale pathological whole genome-wide omics data will become a powerful set for genome-phenome level association studies to discover causal markers for the prediction and diagnosis of health conditions in future studies.

Microbiome analysis of circulating bacterial extracellular vesicles in obsessive-compulsive disorder

AIM

The present study examined the microbiome abundance and composition of drug-naive or drug-free patients with obsessive-compulsive disorder (OCD) compared with healthy controls. In addition, in the OCD group, the microbiome composition was compared between early-onset and late-onset OCD.

METHODS

Serum samples were collected from 89 patients with OCD and 107 age- and sex-matched healthy controls. Bacterial DNA was isolated from bacteria-derived extracellular vesicles in serum and then amplified and quantified using primers specific to the V3-V4 hypervariable region of the 16S ribosomal RNA gene. The 16S ribosomal DNA gene amplicon sequencing was performed.

RESULTS

The pooled estimate showed that α-diversity was significantly reduced in patients with OCD compared with that in healthy controls (PShannon  = 0.00015). In addition, a statistically significant difference was observed in β-diversity between patients with OCD and healthy controls at the order (P = 0.012), family (P = 0.003), genus (P < 0.001), and species (P = 0.005) levels. In the microbiome composition, Pseudomonas, Caulobacteraceae (f), Streptococcus, Novosphingobium, and Enhydrobacter at the genus level were significantly less prevalent in patients with OCD than in controls. In addition, among patients with OCD, the microbial composition in the early-onset versus late-onset types was significantly different with respect to the genera Corynebacterium and Pelomonas.

CONCLUSION

The present study showed an aberrant microbiome in patients with OCD, suggesting a role of the microbiota-brain interaction in the pathophysiology of OCD. Further longitudinal studies with larger sample sizes adjusting for various confounders are warranted.

Genotype imputation accuracy and the quality metrics of the minor ancestry in multi-ancestry reference panels

Large-scale imputation reference panels are currently available and have contributed to efficient genome-wide association studies through genotype imputation. However, whether large-size multi-ancestry or small-size population-specific reference panels are the optimal choices for under-represented populations continues to be debated. We imputed genotypes of East Asian (180k Japanese) subjects using the Trans-Omics for Precision Medicine reference panel and found that the standard imputation quality metric (Rsq) overestimated dosage r2 (squared correlation between imputed dosage and true genotype) particularly in marginal-quality bins. Variance component analysis of Rsq revealed that the increased imputed-genotype certainty (dosages closer to 0, 1 or 2) caused upward bias, indicating some systemic bias in the imputation. Through systematic simulations using different template switching rates (θ value) in the hidden Markov model, we revealed that the lower θ value increased the imputed-genotype certainty and Rsq; however, dosage r2 was insensitive to the θ value, thereby causing a deviation. In simulated reference panels with different sizes and ancestral diversities, the θ value estimates from Minimac decreased with the size of a single ancestry and increased with the ancestral diversity. Thus, Rsq could be deviated from dosage r2 for a subpopulation in the multi-ancestry panel, and the deviation represents different imputed-dosage distributions. Finally, despite the impact of the θ value, distant ancestries in the reference panel contributed only a few additional variants passing a predefined Rsq threshold. We conclude that the θ value substantially impacts the imputed dosage and the imputation quality metric value.

Genetic influences on alcohol flushing in East Asian populations

BACKGROUND

Although it is known that variation in the aldehyde dehydrogenase 2 (ALDH2) gene family influences the East Asian alcohol flushing response, knowledge about other genetic variants that affect flushing symptoms is limited.

METHODS

We performed a genome-wide association study meta-analysis and heritability analysis of alcohol flushing in 15,105 males of East Asian ancestry (Koreans and Chinese) to identify genetic associations with alcohol flushing. We also evaluated whether self-reported flushing can be used as an instrumental variable for alcohol intake.

RESULTS

We identified variants in the region of ALDH2 strongly associated with alcohol flushing, replicating previous studies conducted in East Asian populations. Additionally, we identified variants in the alcohol dehydrogenase 1B (ADH1B) gene region associated with alcohol flushing. Several novel variants were identified after adjustment for the lead variants (ALDH2-rs671 and ADH1B-rs1229984), which need to be confirmed in larger studies. The estimated SNP-heritability on the liability scale was 13% (S.E. = 4%) for flushing, but the heritability estimate decreased to 6% (S.E. = 4%) when the effects of the lead variants were controlled for. Genetic instrumentation of higher alcohol intake using these variants recapitulated known associations of alcohol intake with hypertension. Using self-reported alcohol flushing as an instrument gave a similar association pattern of higher alcohol intake and cardiovascular disease-related traits (e.g. stroke).

CONCLUSION

This study confirms that ALDH2-rs671 and ADH1B-rs1229984 are associated with alcohol flushing in East Asian populations. Our findings also suggest that self-reported alcohol flushing can be used as an instrumental variable in future studies of alcohol consumption.

Genetic assessment of pathogenic germline alterations in lysosomal genes among Asian patients with pancreatic ductal adenocarcinoma

BACKGROUND

Lysosomes are closely linked to autophagic activity, which plays a vital role in pancreatic ductal adenocarcinoma (PDAC) biology. The survival of PDAC patients is still poor, and the identification of novel genetic factors for prognosis and treatment is highly required to prevent PDAC-related deaths. This study investigated the germline variants related to lysosomal dysfunction in patients with PDAC and to analyze whether they contribute to the development of PDAC.

METHODS

The germline putative pathogenic variants (PPV) in genes involved in lysosomal storage disease (LSD) was compared between patients with PDAC (n = 418) and healthy controls (n = 845) using targeted panel and whole-exome sequencing. Furthermore, pancreatic organoids from wild-type and KrasG12D mice were used to evaluate the effect of lysosomal dysfunction on PDAC development. RNA sequencing (RNA-seq) analysis was performed with established PDAC patient-derived organoids (PDOs) according to the PPV status.

RESULTS

The PPV in LSD-related genes was higher in patients with PDAC than in healthy controls (8.13 vs. 4.26%, Log2 OR = 1.65, P = 3.08 × 10-3). The PPV carriers of LSD-related genes with PDAC were significantly younger than the non-carriers (mean age 61.5 vs. 65.3 years, P = 0.031). We further studied a variant of the lysosomal enzyme, galactosylceramidase (GALC), which was the most frequently detected LSD variant in our cohort. Autophagolysosomal activity was hampered when GALC was downregulated, which was accompanied by paradoxically elevated autophagic flux. Furthermore, the number of proliferating Ki-67+ cells increased significantly in pancreatic organoids derived from Galc knockout KrasG12D mice. Moreover, GALC PPV carriers tended to show drug resistance in both PDAC cell line and PDAC PDO, and RNA-seq analysis revealed that various metabolism and gene repair pathways were upregulated in PDAC PDOs harboring a GALC variant.

CONCLUSIONS

Genetically defined lysosomal dysfunction is frequently observed in patients with young-onset PDAC. This might contribute to PDAC development by altering metabolism and impairing autophagolysosomal activity, which could be potentially implicated in therapeutic applications for PDAC.

The complete and fully-phased diploid genome of a male Han Chinese

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Personalised circulating tumour DNA assay with large-scale mutation coverage for sensitive minimal residual disease detection in colorectal cancer

BACKGROUND

Postoperative minimal residual disease (MRD) detection using circulating-tumour DNA (ctDNA) requires a highly sensitive analysis platform. We have developed a tumour-informed, hybrid-capture ctDNA sequencing MRD assay.

METHODS

Personalised target-capture panels for ctDNA detection were designed using individual variants identified in tumour whole-exome sequencing of each patient. MRD status was determined using ultra-high-depth sequencing data of plasma cell-free DNA. The MRD positivity and its association with clinical outcome were analysed in Stage II or III colorectal cancer (CRC).

RESULTS

In 98 CRC patients, personalised panels for ctDNA sequencing were built from tumour data, including a median of 185 variants per patient. In silico simulation showed that increasing the number of target variants increases MRD detection sensitivity in low fractions (<0.01%). At postoperative 3-week, 21.4% of patients were positive for MRD by ctDNA. Postoperative positive MRD was strongly associated with poor disease-free survival (DFS) (adjusted hazard ratio 8.40, 95% confidence interval 3.49-20.2). Patients with a negative conversion of MRD after adjuvant therapy showed significantly better DFS (P < 0.001).

CONCLUSION

Tumour-informed, hybrid-capture-based ctDNA assay monitoring a large number of patient-specific mutations is a sensitive strategy for MRD detection to predict recurrence in CRC.

Identification and validation of six acute myocardial infarction-associated variants, including a novel prognostic marker for cardiac mortality

Background

Acute myocardial infarction (AMI) is one of the leading causes of death worldwide, and approximately half of AMI-related deaths occur before the affected individual reaches the hospital. The present study aimed to identify and validate genetic variants associated with AMI and their role as prognostic markers.

Materials and methods

We conducted a replication study of 29 previously identified novel loci containing 85 genetic variants associated with early-onset AMI using a new independent set of 2,920 Koreans [88 patients with early- and 1,085 patients with late-onset AMI, who underwent percutaneous coronary intervention (PCI), and 1,747 healthy controls].

Results

Of the 85 previously reported early-onset variants, six were confirmed in our genome-wide association study with a false discovery rate of less than 0.05. Notably, rs12639023, a cis-eQTL located in the intergenic region between LINC02005 and CNTN3, significantly increased longitudinal cardiac mortality and recurrent AMI. CNTN3 is known to play a role in altering vascular permeability. Another variant, rs78631167, located upstream of PLAUR and known to function in fibrinolysis, was moderately replicated in this study. By surveying the nearby genomic region around rs78631167, we identified a significant novel locus (rs8109584) located 13 bp downstream of rs78631167. The present study showed that six of the early-onset variants of AMI are applicable to both early- and late-onset cases.

Conclusion

Our results confirm markers that can potentially be utilized to predict, screen, prevent, and treat candidate patients with AMI and highlight the potential of rs12639023 as a prognostic marker for cardiac mortality in AMI.

Genomic characteristics of triple negative apocrine carcinoma: a comparison to triple negative breast cancer

Apocrine carcinoma is a rare breast cancer subtype. As such, the genomic characteristics of apocrine carcinoma with triple negative immunohistochemical results (TNAC), which has been treated as triple negative breast cancer (TNBC), have not been revealed. In this study, we evaluated the genomic characteristics of TNAC compared to TNBC with low Ki-67 (LK-TNBC). In the genetic analysis of 73 TNACs and 32 LK-TNBCs, the most frequently mutated driver gene in TNAC was TP53 (16/56, 28.6%), followed by PIK3CA (9/56, 16.1%), ZNF717 (8/56, 14.3%), and PIK3R1 (6/56, 10.71%). Mutational signature analysis showed enrichment of defective DNA mismatch repair (MMR)-related signatures (SBS6 and SBS21) and the SBS5 signature in TNAC, whereas an APOBEC activity-associated mutational signature (SBS13) was more prominent in LK-TNBC (Student's t test, p < 0.05). In intrinsic subtyping, 38.4% of TNACs were classified as luminal A, 27.4% as luminal B, 26.0% as HER2-enriched (HER2-E), 2.7% as basal, and 5.5% as normal-like. The basal subtype was the most dominant subtype (43.8%) in LK-TNBC (p < 0.001), followed by luminal B (21.9%), HER2-E (21.9%), and luminal A (12.5%). In the survival analysis, TNAC had a five-year disease-free survival (DFS) rate of 92.2% compared to 59.1% for LK-TNBC (P = 0.001) and a five-year overall survival (OS) rate of 95.3% compared to 74.6% for LK-TNBC (P = 0.0099). TNAC has different genetic characteristics and better survival outcomes than LK-TNBC. In particular, normal-like and luminal A subtypes in TNAC have much better DFS and OS than other intrinsic subtypes. Our findings are expected to impact medical practice for patients diagnosed with TNAC.

Integrative analysis of risk factors for immune-related adverse events of checkpoint blockade therapy in cancer

Immune-related adverse events (irAEs) induced by checkpoint inhibitors involve a multitude of different risk factors. Here, to interrogate the multifaceted underlying mechanisms, we compiled germline exomes and blood transcriptomes with clinical data, before and after checkpoint inhibitor treatment, from 672 patients with cancer. Overall, irAE samples showed a substantially lower contribution of neutrophils in terms of baseline and on-therapy cell counts and gene expression markers related to neutrophil function. Allelic variation of HLA-B correlated with overall irAE risk. Analysis of germline coding variants identified a nonsense mutation in an immunoglobulin superfamily protein, TMEM162. In our cohort and the Cancer Genome Atlas (TCGA) data, TMEM162 alteration was associated with higher peripheral and tumor-infiltrating B cell counts and suppression of regulatory T cells in response to therapy. We developed machine learning models for irAE prediction, validated using additional data from 169 patients. Our results provide valuable insights into risk factors of irAE and their clinical utility.

A heterozygous mutation in UBE2H in a patient with developmental delay leads to an aberrant brain development in zebrafish

BACKGROUND

Ubiquitin-related rare diseases are generally characterized by developmental delays and mental retardation, but the exact incidence or prevalence is not yet fully understood. The clinical application of next-generation sequencing for pediatric seizures and developmental delay of unknown causes has become common in studies aimed at identification of a causal gene in patients with ubiquitin-related rare diseases that cannot be diagnosed using conventional fluorescence in situ hybridization or chromosome microarray tests. Our study aimed to investigate the effects of ubiquitin-proteasome system on ultra-rare neurodevelopmental diseases, through functional identification of candidate genes and variants.

METHODS

In our present work, we carried out genome analysis of a patient with clinical phenotypes of developmental delay and intractable convulsion, to identify causal mutations. Further characterization of the candidate gene was performed using zebrafish, through gene knockdown approaches. Transcriptomic analysis using whole embryos of zebrafish knockdown morphants and additional functional studies identified downstream pathways of the candidate gene affecting neurogenesis.

RESULTS

Through trio-based whole-genome sequencing analysis, we identified a de novo missense variant of the ubiquitin system-related gene UBE2H (c.449C>T; p.Thr150Met) in the proband. Using zebrafish, we found that Ube2h is required for normal brain development. Differential gene expression analysis revealed activation of the ATM-p53 signaling pathway in the absence of Ube2h. Moreover, depletion of ube2h led to induction of apoptosis, specifically in the differentiated neural cells. Finally, we found that a missense mutation in zebrafish, ube2h (c.449C>T; p.Thr150Met), which mimics a variant identified in a patient with neurodevelopmental defects, causes aberrant Ube2h function in zebrafish embryos.

CONCLUSION

A de novo heterozygous variant in the UBE2H c.449C>T (p.Thr150Met) has been identified in a pediatric patient with global developmental delay and UBE2H is essential for normal neurogenesis in the brain.

Genetic investigation of 124 SNPs in a Myanmar population using the Precision ID Identity Panel and the Illumina MiSeq

Single nucleotide polymorphisms (SNPs) have become popular in forensic genetics as an alternative to short tandem repeats (STRs). The Precision ID Identity Panel (Thermo Fisher Scientific), consisting of 90 autosomal SNPs and 34 Y-chromosomal SNPs, enabled human identification studies on global populations through next-generation sequencing (NGS). However, most previous studies on the panel have used the Ion Torrent platform, and there are few reports on the Southeast Asian population. Here, a total of 96 unrelated males from Myanmar (Yangon) were analyzed with the Precision ID Identity Panel on a MiSeq (Illumina) using an in-house TruSeq compatible universal adapter and a custom variant caller, Visual SNP. The sequencing performance evaluated by locus balance and heterozygote balance was comparable to that of the Ion Torrent platform. For 90 autosomal SNPs, the combined match probability (CMP) was 6.994 × 10^-34, lower than that of 22 PowerPlex Fusion autosomal STRs (3.130 × 10^-26). For 34 Y-SNPs, 14 Y-haplogroups (mostly O2 and O1b) were observed. We found 51 cryptic variations (42 haplotypes) around target SNPs, of which haplotypes corresponding to 33 autosomal SNPs decreased CMP. Interpopulation analysis revealed that the Myanmar population is genetically closer to the East and Southeast Asian populations. In conclusion, the Precision ID Identity Panel can be successfully analyzed on the Illumina MiSeq and provides high discrimination power for human identification in the Myanmar population. This study broadened the accessibility of the NGS-based SNP panel by expanding the available NGS platforms and adopting a robust NGS data analysis tool.

Exome sequencing identified five novel USH2A variants in Korean patients with retinitis pigmentosa

BACKGROUND

Retinitis pigmentosa (RP) is an inherited disorder that causes progressive loss of vision. This study aimed to describe the possible causative variants of the USH2A gene in Korean RP families and their associated phenotypes.

MATERIALS AND METHODS

We recruited 94 RP families (220 subjects, including 94 probands and 126 family members) in a Korean cohort, and analyzed USH2A gene variants through whole-exome sequencing. The pathogenicity of the variants was classified according to American College of Medical Genetics and Genomics and Association for Molecular Pathology guidelines.

RESULTS

We found 14 USH2A disease-causing variants, including 5 novel variants. Disease causing variants were identified in 10 probands with RP, accounting for 10.6% (10/94) of the Korean RPs in the cohort. To visually represent the structural changes induced by novel variants, we modeled the three-dimensional structures of the wild-type and mutant proteins.

CONCLUSIONS

This study expands the spectrum of USH2A variants and provides information for future therapeutic strategies for RP.

Whole-genome sequencing reveals an association between small genomic deletions and an increased risk of developing Parkinson's disease

Single-nucleotide variants (SNVs) associated with Parkinson's disease (PD) have been investigated mainly through genome-wide association studies. However, other genomic alterations, including copy number variations, remain less explored. In this study, we conducted whole-genome sequencing of primary (310 PD patients and 100 healthy individuals) and independent (100 PD patients and 100 healthy individuals) cohorts from the Korean population to identify high-resolution small genomic deletions, gains, and SNVs. Global small genomic deletions and gains were found to be associated with an increased and decreased risk of PD development, respectively. Thirty significant locus deletions were identified in PD, with most being associated with an increased PD risk in both cohorts. Small genomic deletions in clustered loci located in the GPR27 region had high enhancer signals and showed the closest association with PD. GPR27 was found to be expressed specifically in brain tissue, and GPR27 copy number loss was associated with upregulated SNCA expression and downregulated dopamine neurotransmitter pathways. Clustering of small genomic deletions on chr20 in exon 1 of the GNAS isoform was detected. In addition, we found several PD-associated SNVs, including one in the enhancer region of the TCF7L2 intron, which exhibited a cis-acting regulatory mode and an association with the beta-catenin signaling pathway. These findings provide a global, whole-genome view of PD and suggest that small genomic deletions in regulatory domains contribute to the risk of PD development.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Genome-wide analyses of early-onset acute myocardial infarction identify 29 novel loci by whole genome sequencing

Early-onset acute myocardial infarction (AMI) may have a higher genetic predisposition than late-onset AMI. The present study aimed to identify and characterize germline variants that affect early-onset AMI using whole-genome sequencing (WGS). We performed a genome-wide association study based on the WGS of 1239 Koreans, including 596 early-onset AMI patients and 643 healthy individuals. Patients with AMI who underwent percutaneous coronary intervention (PCI) caused by atherothrombotic occlusive lesions were included in the study. A total of 29 novel loci were found to be associated with early-onset AMI. These loci are involved in thrombosis, fibrinolysis, inflammation, and lipid metabolism. One of the associated single nucleotide variants (SNVs), rs1614576, located upstream of PRKCB, is known to be associated with thrombus formation. Additionally, the results revealed a novel locus, rs78631167, located upstream of PLAUR which plays a critical role in regulating plasminogen activation and is related to fibrinolysis. The association between early-onset AMI and rs9357455, which is located upstream of PHACTR1 and regulates inflammation in AMI, was found. Moreover, we identified a lipid metabolism related genetic risk locus, rs5072, in the APOA1-AS gene. This study provides new evidence supporting the genetic association between early-onset AMI and thrombosis and fibrinolysis, as well as inflammation and lipid metabolism, by analyzing the whole-genome of 596 patients with early-onset AMI who have been treated with PCI. Our findings highlight potential genetic markers for the prediction and management of AMI, as well as for understanding the etiology of AMI.

A Low-Frequency APOB p.(Pro955Ser) Variant Contributes to the Severity of/Variability in Familial Hypercholesterolemia

CONTEXT

Heterozygous familial hypercholesterolemia (HeFH) is caused by a rare pathogenic variant in the LDLR, APOB, and PCSK9 genes. However, the causative variants in these genes have not been identified in approximately 40% of HeFH patients.

OBJECTIVE

Our aim was to identify novel (or additional) genes/variants that contribute to HeFH.

METHODS

Whole-exome sequencing was performed for 215 family members from 122 families with HeFH without pathogenic variants in the LDLR or PCSK9 genes.

RESULTS

We could not find novel causative familial hypercholesterolemia (FH) genes/variants by family analysis. Next, we examined all APOB variants. Twenty-four nonsynonymous APOB variants were identified. The allele frequencies of the c.2863C > T:p.(Pro955Ser) variant in the HeFH probands and the general Japanese population were 0.15 and 0.034, respectively [odds ratio 4.9 (95% CI 3.4-7.1); P = 6.9 × 10-13]. The patients harboring the c.2863C > T:p.(Pro955Ser) variant accounted for 9.8% (n = 63) of unrelated patients with HeFH (n = 645). The penetrance of the c.2863C > T:p.(Pro955Ser) variant was low in the pedigree-based genetic analysis. In an in vitro assay, low-density lipoprotein (LDL) uptake from patients with the homozygous c.2863C > T:p.(Pro955Ser) variant was 44% of the LDL uptake from control subjects, and it was similar to that of the LDL uptake from patients with the known pathogenic heterozygous p.(Arg3527Gln) variant.

CONCLUSIONS

The low-frequency APOB c.2863C > T:p.(Pro955Ser) variant is not an FH-causative variant, but it has a moderate effect size in HeFH. These findings suggest that the combination of the APOB c.2863C > T:p.(Pro955Ser) variant and age, environmental factors, or other genetic factors contributes to the severity of or variability in the HeFH phenotype.

A database of 5305 healthy Korean individuals reveals genetic and clinical implications for an East Asian population

Despite substantial advances in disease genetics, studies to date have largely focused on individuals of European descent. This limits further discoveries of novel functional genetic variants in other ethnic groups. To alleviate the paucity of East Asian population genome resources, we established the Korean Variant Archive 2 (KOVA 2), which is composed of 1896 whole-genome sequences and 3409 whole-exome sequences from healthy individuals of Korean ethnicity. This is the largest genome database from the ethnic Korean population to date, surpassing the 1909 Korean individuals deposited in gnomAD. The variants in KOVA 2 displayed all the known genetic features of those from previous genome databases, and we compiled data from Korean-specific runs of homozygosity, positively selected intervals, and structural variants. In doing so, we found loci, such as the loci of ADH1A/1B and UHRF1BP1, that are strongly selected in the Korean population relative to other East Asian populations. Our analysis of allele ages revealed a correlation between variant functionality and evolutionary age. The data can be browsed and downloaded from a public website ( https://www.kobic.re.kr/kova/ ). We anticipate that KOVA 2 will serve as a valuable resource for genetic studies involving East Asian populations.

PharmaKoVariome database for supporting genetic testing

Pharmacogenomics (PGx) provides information about routine precision medicine, based on the patient's genotype. However, many of the available information about human allele frequencies, and about clinical drug-gene interactions, is based on American and European populations. PharmaKoVariome database was constructed to support genetic testing for safe prescription and drug development. It consolidated and stored 2507 diseases, 11 459 drugs and 61 627 drug-target or druggable genes from public databases. PharmaKoVariome precomputed ethnic-specific abundant variants for approximately 120 M single-nucleotide variants of drug-target or druggable genes. A user can search by gene symbol, drug name, disease and reference SNP ID number (rsID) to statistically analyse the frequency of ethnical variations, such as odds ratio and P-values for related genes. In an example study, we observed five Korean-enriched variants in the CYP2B6 and CYP2D6 genes, one of which (rs1065852) is known to be incapable of metabolizing drug. It is also shown that 4-6% of North and East Asians have risk factors for drugs metabolized by the CYP2D6 gene. Therefore, PharmaKoVariome is a useful database for pharmaceutical or diagnostic companies for developing diagnostic technologies that can be applied in the Asian PGx industry. Database URL: http://www.pharmakovariome.com/.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Precision Oncology in Canada: Converting Vision to Reality with Lessons from International Programs

Canada's healthcare system, like others worldwide, is immersed in a process of evolution, attempting to adapt conventional frameworks of health technology assessment (HTA) and funding models to a new landscape of precision medicine in oncology. In particular, the need for real-world evidence in Canada is not matched by the necessary infrastructure and technologies required to integrate genomic and clinical data. Since healthcare systems in many developed nations face similar challenges, we adopted a solutions-based approach and conducted a search of worldwide programs in personalized medicine, with an emphasis on precision oncology. This search strategy included review articles published between 1 January 2016 and 1 March 2021 and hand-searches of their reference lists for relevant publications back to 1 December 2005. Thirty-nine initiatives across 37 countries in Europe, Australasia, Africa, and the Americas had the potential to lead to real-world data (RWD) on the clinical utility of oncology biomarkers. We highlight four initiatives with helpful lessons for Canada: Genomic Medicine France 2025, UNICANCER, the German Medical Informatics Initiative, and CANCER-ID. Among the 35 other programs evaluated, the main themes included the need for collaboration and systems to support data harmonization across multiple jurisdictions. In order to generate RWD in precision oncology that will prove acceptable to HTA bodies, Canada must take a national approach to biomarker strategy and unite all stakeholders at the highest level to overcome jurisdictional and technological barriers.

Human genetics: The dual origin of Three Kingdoms period Koreans

The genetic history of Koreans remains poorly understood due to a lack of ancient DNA. A new paleo-genomic study shows that population stratification in 4^th-5^th century South Korean populations was linked to a varied proportion of indigenous Jomon-related ancestry, which does not survive in present-day Koreans.

Infantile esotropia in a family with TUBB3 mutation associated congenital fibrosis of extraocular muscles

BACKGROUND

The TUBB3 gene has been reported to be associated with type 3 congenital fibrosis of the extraocular muscles (CFEOM). The clinical features of CFEOM3 that are linked to TUBB3 mutations are diverse, ranging from mild ptosis and limitation of extraocular movement to severe ocular motility problems and central nervous system abnormalities.

MATERIALS AND METHODS

This was a single retrospective case report.

RESULT

This case report describes a patient with infantile esotropia, who had a heterozygous variant in TUBB3 c.904 G > A (p.A302T) known to cause CFEOM3 and her family members, who presented with manifestations associated with CFEOM3.

CONCLUSION

Given the diverse clinical features of CFEOM3, the possibility of the occurrence of CFEOM3 should be considered when there is a congenital abnormality of extraocular muscle movement and a positive family history.

Northeastern Asian and Jomon-related genetic structure in the Three Kingdoms period of Gimhae, Korea

The genetic history of prehistoric and protohistoric Korean populations is not well understood because only a small number of ancient genomes are available. Here, we report the first paleogenomic data from the Korean Three Kingdoms period, a crucial point in the cultural and historic formation of Korea. These data comprise eight shotgun-sequenced genomes from ancient Korea (0.7×-6.1× coverage). They were derived from two archeological sites in Gimhae: the Yuha-ri shell mound and the Daesung-dong tumuli, the latter being the most important funerary complex of the Gaya confederacy. All individuals are from between the 4th and 5th century CE and are best modeled as an admixture between a northern China Bronze Age genetic source and a source of Jomon-related ancestry that shares similarities with the present-day genomes from Japan. The observed substructure and proportion of Jomon-related ancestry suggest the presence of two genetic groups within the population and diversity among the Gaya population. We could not correlate the genomic differences between these two groups with either social status or sex. All the ancient individuals' genomic profiles, including phenotypically relevant SNPs associated with hair and eye color, facial morphology, and myopia, imply strong genetic and phenotypic continuity with modern Koreans for the last 1,700 years.

EASplex: A panel of 308 AISNPs for East Asian ancestry inference using next generation sequencing

Ancestry inference is useful in many scientific fields, such as forensic genetics, medical genetics, and molecular archaeology. Various ancestry inferring methods have been released for major continental populations. However, few reports refer to sub-populations within the East Asian population using hundreds of ancestry informative SNPs (AISNPs). In this study, we developed a 308-AISNP panel (EASplex NGS DNA panel) using multiplex PCR and next generation sequencing (NGS). This panel included 56 SNPs relevant for the continent-level ancestry inference and 252 Japanese-, Korean-, and/or Han Chinese-specific AISNPs to address the ancestry inference of global populations and regional populations among Japanese (JPT), Korean minority (CHK), and Han Chinese (CHH). A total of 87 CHK and 59 CHH samples were used to check the performance of the EASplex NGS DNA panel. By analyzing 146 profiles of samples with JPT and CHH data from Beijing and South China in 1000 genomes project, the following results were obtained: (1) the 146 tested samples were correctly assigned to the East Asian group; (2) the paired population assignment rate was 99.73% for JPT and CHH, 95% for JPT and CHK, and 90.11% for CHK and CHH; and (3) the whole population assignment was 92.14% for the JPT, CHK, and CHH data. Overall, the EASplex NGS DNA panel displayed informativeness for continental ancestry inference and regional ancestry inference among JPT, CHH, and CHK and has the potential for use in forensic and genetic population studies.

Genomic analyses of 10,376 individuals in the Westlake BioBank for Chinese (WBBC) pilot project

We initiate the Westlake BioBank for Chinese (WBBC) pilot project with 4,535 whole-genome sequencing (WGS) individuals and 5,841 high-density genotyping individuals, and identify 81.5 million SNPs and INDELs, of which 38.5% are absent in dbSNP Build 151. We provide a population-specific reference panel and an online imputation server ( https://wbbc.westlake.edu.cn/ ) which could yield substantial improvement of imputation performance in Chinese population, especially for low-frequency and rare variants. By analyzing the singleton density of the WGS data, we find selection signatures in SNX29, DNAH1 and WDR1 genes, and the derived alleles of the alcohol metabolism genes (ADH1A and ADH1B) emerge around 7,000 years ago and tend to be more common from 4,000 years ago in East Asia. Genetic evidence supports the corresponding geographical boundaries of the Qinling-Huaihe Line and Nanling Mountains, which separate the Han Chinese into subgroups, and we reveal that North Han was more homogeneous than South Han.

Shared genetic architectures of subjective well-being in East Asian and European ancestry populations

Subjective well-being (SWB) has been explored in European ancestral populations; however, whether the SWB genetic architecture is shared across populations remains unclear. We conducted a cross-population genome-wide association study for SWB using samples from Korean (n = 110,919) and European (n = 563,176) ancestries. Five ancestry-specific loci and twelve cross-ancestry significant genomic loci were identified. One novel locus (rs12298541 near HMGA2) associated with SWB was also identified through the European meta-analysis. Significant cross-ancestry genetic correlation for SWB between samples was observed. Polygenic risk analysis in an independent Korean cohort (n = 22,455) demonstrated transferability between populations. Significant correlations between SWB and major depressive disorder, and significant enrichment of central nervous system-related polymorphisms heritability in both ancestry populations were found. Hence, large-scale cross-ancestry genome-wide association studies can advance our understanding of SWB genetic architecture and mental health.

The Thousand Polish Genomes-A Database of Polish Variant Allele Frequencies

Although Slavic populations account for over 4.5% of world inhabitants, no centralised, open-source reference database of genetic variation of any Slavic population exists to date. Such data are crucial for clinical genetics, biomedical research, as well as archeological and historical studies. The Polish population, which is homogenous and sedentary in its nature but influenced by many migrations of the past, is unique and could serve as a genetic reference for the Slavic nations. In this study, we analysed whole genomes of 1222 Poles to identify and genotype a wide spectrum of genomic variation, such as small and structural variants, runs of homozygosity, mitochondrial haplogroups, and de novo variants. Common variant analyses showed that the Polish cohort is highly homogenous and shares ancestry with other European populations. In rare variant analyses, we identified 32 autosomal-recessive genes with significantly different frequencies of pathogenic alleles in the Polish population as compared to the non-Finish Europeans, including C2, TGM5, NUP93, C19orf12, and PROP1. The allele frequencies for small and structural variants, calculated for 1076 unrelated individuals, are released publicly as The Thousand Polish Genomes database, and will contribute to the worldwide genomic resources available to researchers and clinicians.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

BACKGROUND

Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population.

METHOD

Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation.

RESULTS

A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively.

CONCLUSIONS

Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

Mutational Analysis of Triple-Negative Breast Cancer Using Targeted Kinome Sequencing

Purpose

Triple-negative breast cancer (TNBC) does not have defined therapeutic targets and is currently treated with chemotherapy only. Kinase dysregulation triggers cancer cell proliferation and metastasis and is a crucial therapeutic target for cancer. In this study, targeted kinome sequencing of TNBC tumors was performed to assess the association between kinome gene alterations and disease outcomes in TNBC.

Methods

A kinome gene panel consisting of 612 genes was used for the targeted sequencing of 166 TNBC samples and matched normal tissues. Analyses of the significantly mutated genes were performed. Genomic differences between Asian and non-Asian patients with TNBC were evaluated using two Asian TNBC datasets (from Seoul National University Hospital [SNUH] and Fudan University Shanghai Cancer Center [FUSCC]) and three non-Asian TNBC datasets (The Cancer Genome Atlas [TCGA], METABRIC, and Gustave Roussy). The prognostic value of kinome gene mutations was evaluated using tumor mutational burden (TMB) and oncogenic pathway analyses. Mutational profiles from the TCGA were used for validation.

Results

The significantly mutated genes included TP53 (60% of patients), PIK3CA (21%), BRCA2 (8%), and ATM (8%). Compared with data from non-Asian public databases, the mutation rates of PIK3CA p.H1047R/Q were significantly higher in the SNUH cohort (p = 0.003, 0.048, and 0.032, respectively). This was verified using the FUSCC dataset (p = 0.003, 0.078, and 0.05, respectively). The TMB-high group showed a trend toward longer progression-free survival in our cohort and the TCGA TNBC cohort (p = 0.041 and 0.195, respectively). Kinome gene alterations in the Wnt pathway in patients with TNBC were associated with poor survival in both datasets (p = 0.002 and 0.003, respectively).

Conclusion

Comprehensive analyses of kinome gene alterations in TNBC revealed genomic alterations that offer therapeutic targets and should help identify high-risk patients more precisely in future studies.

Analyzing the Korean reference genome with meta-imputation increased the imputation accuracy and spectrum of rare variants in the Korean population

Genotype imputation is essential for enhancing the power of association-mapping and discovering rare and indels that are missed by most genotyping arrays. Imputation analysis can be more accurate with a population-specific reference panel or a multi-ethnic reference panel with numerous samples. The National Institute of Health, Republic of Korea, initiated the Korean Reference Genome (KRG) project to identify variants in whole-genome sequences of ∼20,000 Korean participants. In the pilot phase, we analyzed the data from 1,490 participants. The genetic characteristics and imputation performance of the KRG were compared with those of the 1,000 Genomes Project Phase 3, GenomeAsia 100K Project, ChinaMAP, NARD, and TOPMed reference panels. For comparison analysis, genotype panels were artificially generated using whole-genome sequencing data from combinations of four different ancestries (Korean, Japanese, Chinese, and European) and two population-specific optimized microarrays (Korea Biobank Array and UK Biobank Array). The KRG reference panel performed best for the Korean population (R ² = 0.78-0.84, percentage of well-imputed is 91.9% for allele frequency >5%), although the other reference panels comprised a larger number of samples with genetically different background. By comparing multiple reference panels and multi-ethnic genotype panels, optimal imputation was obtained using reference panels from genetically related populations and a population-optimized microarray. Indeed, the reference panels of KRG and TOPMed showed the best performance when applied to the genotype panels of KBA (R ² = 0.84) and UKB (R ² = 0.87), respectively. Using a meta-imputation approach to merge imputation results from different reference panels increased the imputation accuracy for rare variants (∼7%) and provided additional well-imputed variants (∼20%) with comparable imputation accuracy to that of the KRG. Our results demonstrate the importance of using a population-specific reference panel and meta-imputation to assess a substantial number of accurately imputed rare variants.

A retrospective single-centered, comprehensive targeted genetic sequencing analysis of prognostic survival using tissues from Korean patients with metastatic renal cell carcinoma after targeted therapy

Purpose

To identify candidate gene mutations to significantly predict the risk of survival prognosis after treatment with systemic first-line targeted therapy (TT) in metastatic renal cell carcinoma (mRCC) patients.

Materials and Methods

Between 2005 and 2017, 168 triplet-tissue block samples from 56 mRCC patients were selected for targeted gene sequencing (TGS). Fifty-six patients’ medical records including overall survival (OS) and progression-free survival (PFS) at the time of mRCC diagnosis were evaluated. The patients were grouped into favorable (>12 months/>3 years), intermediate (3–12/12–36 months), and poor groups according to their PFS/OS (<3 months/<12 months). We identified any significant therapeutic targeted genes relating to the survival with a significance at p<0.050.

Results

The first line therapeutic response showed 1.8% complete remission, 14.2% partial response, 42.9% stable disease, and 41.1% progressive disease. Among the overall TGS results, the cumulative effect of CDH1, and/or PTK2 genes significantly reflected the therapeutic responses in terms of PFS/OS; CDH1 and PTK2 mutations were associated with poor prognostic outcomes (p<0.050). Among only triplet-quality check passed tissues, the SGO2, BRAF, URB1, and NEDD1 mutated genes significantly correlated with OS. Regarding metastasis, patients with liver metastasis had the worst OS (p=0.050). The combinational mutation number from these two candidate genes in the liver metastatic samples with mutated EGFR2 and FABP7 also showed a significantly worse OS than those with other metastatic lesions (p<0.050).

Conclusions

This study reports several significant mutated genes related to the survival prognosis in mRCC patients treated with first-line TT.

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.

DNA bridges: A novel platform for single-molecule sequencing and other DNA-protein interaction applications

DNA molecular combing is a technique that stretches thousands of long individual DNA molecules (up to 10 Mbp) into a parallel configuration on surface. It has previously been proposed to sequence these molecules by synthesis. However, this approach poses two critical challenges: 1-Combed DNA molecules are overstretched and therefore a nonoptimal substrate for polymerase extension. 2-The combing surface sterically impedes full enzymatic access to the DNA backbone. Here, we introduce a novel approach that attaches thousands of molecules to a removable surface, with a tunable stretching factor. Next, we dissolve portions of the surface, leaving the DNA molecules suspended as 'bridges'. We demonstrate that the suspended molecules are enzymatically accessible, and we have used an enzyme to incorporate labeled nucleotides, as predicted by the specific molecular sequence. Our results suggest that this novel platform is a promising candidate to achieve high-throughput sequencing of Mbp-long molecules, which could have additional genomic applications, such as the study of other protein-DNA interactions.

Structural variants in the Chinese population and their impact on phenotypes, diseases and population adaptation

A complete characterization of genetic variation is a fundamental goal of human genome research. Long-read sequencing has improved the sensitivity of structural variant discovery. Here, we conduct the long-read sequencing-based structural variant analysis for 405 unrelated Chinese individuals, with 68 phenotypic and clinical measurements. We discover a landscape of 132,312 nonredundant structural variants, of which 45.2% are novel. The identified structural variants are of high-quality, with an estimated false discovery rate of 3.2%. The concatenated length of all the structural variants is approximately 13.2% of the human reference genome. We annotate 1,929 loss-of-function structural variants affecting the coding sequence of 1,681 genes. We discover rare deletions in HBA1/HBA2/HBB associated with anemia. Furthermore, we identify structural variants related to immunity which differentiate the northern and southern Chinese populations. Our study describes the landscape of structural variants in the Chinese population and their contribution to phenotypes and disease.

Development and Validation of Targeted Gene Sequencing Panel Based Companion Diagnostic for Korean Patients with Solid Tumors

Simple Summary

We have developed and analytically validated the Korean Pan-cancer Companion Diagnostic (CDX) Panel to apply targeted anticancer drugs to Korean patients based on the molecular characteristics of tumors using tumor samples without matched patient normal samples. The panel included 31 genes with reported single nucleotide variants, 9 genes with reported copy number variations, and 15 genes with predictive responses to targeted drugs under clinical testing, enabling the panel to be analyzed for the targets of 30 targeted anticancer drugs. It is cost-effective and optimized for cancer type-specific therapy in Korean cancer patients across solid cancer types while minimizing the limitations of existing approaches. This gene screening method is expected to reduce test turnaround time and cost, making it a balanced approach to investigate solid cancer-related gene regions.

Abstract

Recently, several panels using two representative targeting methods have been developed but they do not reflect racial specificity, especially for Asians. We have developed and analytically validated the Korean Pan-cancer Companion Diagnostic (CDX) Panel to apply targeted anticancer drugs to Korean patients based on the molecular characteristics of tumors using tumor samples without matched patient normal samples. The panel included 31 genes with reported single nucleotide variants, 9 genes with reported copy number variations, and 15 genes with predictive responses to targeted drugs under clinical testing, enabling the panel to be analyzed for the targets of 30 targeted anticancer drugs. It is cost-effective and optimized for cancer type-specific therapy in Korean cancer patients across solid cancer types while minimizing the limitations of existing approaches. In addition, the optimized filtering protocol for somatic variants from tumor-only samples enables researchers to use this panel without matched normal samples. To verify the panel, 241 frozen tumor tissues and 71 formalin-fixed paraffin-embedded (FFPE) samples from several institutes were registered. This gene screening method is expected to reduce test turnaround time and cost, making it a balanced approach to investigate solid cancer-related gene regions.

Genomic sequencing and functional analyses identify MAP4K3/GLK germline and somatic variants associated with systemic lupus erythematosus

OBJECTIVES

MAP4K3 (GLK) overexpression in T cells induces interleukin (IL)-17A production and autoimmune responses. GLK overexpressing T-cell population is correlated with severity of human systemic lupus erythematosus (SLE); however, it is unclear how GLK is upregulated in patients with SLE.

METHODS

We enrolled 181 patients with SLE and 250 individuals without SLE (93 healthy controls and 157 family members of patients with SLE) in two independent cohorts from different hospitals/cities. Genomic DNAs of peripheral blood mononuclear cells were subjected to next-generation sequencing to identify GLK gene variants. The functional consequences of the identified GLK germline or somatic variants were investigated using site-directed mutagenesis and cell transfection, followed by reporter assays, mass spectrometry, immunoblotting, coimmunoprecipitation, and in situ proximity ligation assays.

RESULTS

We identified 58 patients with SLE from Cohort #1 and #2 with higher frequencies of a somatic variant (chr2:39 477 124 A>G) in GLK 3'-untranslated region (UTR); these patients with SLE showed increased serum anti-double-stranded DNA levels and decreased serum C3/C4 levels. This somatic variant in 3'-UTR enhanced GLK mRNA levels in T cells. In addition, we identified five patients with SLE with GLK (A410T) germline variant in Cohort #1 and #2, as well as two other patients with SLE with GLK (K650R) germline variant in Cohort #1. Another GLK germline variant, A579T, was also detected in one patient with SLE from Cohort #2. Both GLK (A410T) and GLK (K650R) mutants inhibited GLK ubiquitination induced by the novel E3 ligase makorin ring-finger protein 4 (MKRN4), leading to GLK protein stabilisation.

CONCLUSIONS

Multiple GLK germline and somatic variants cause GLK induction by increasing mRNA or protein stability in patients with SLE.

Regional TMPRSS2 V197M Allele Frequencies Are Correlated with COVID-19 Case Fatality Rates

Coronavirus disease, COVID-19 (coronavirus disease 2019), caused by SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), has a higher case fatality rate in European countries than in others, especially East Asian ones. One potential explanation for this regional difference is the diversity of the viral infection efficiency. Here, we analyzed the allele frequencies of a nonsynonymous variant rs12329760 (V197M) in the TMPRSS2 gene, a key enzyme essential for viral infection and found a significant association between the COVID-19 case fatality rate and the V197M allele frequencies, using over 200,000 present-day and ancient genomic samples. East Asian countries have higher V197M allele frequencies than other regions, including European countries which correlates to their lower case fatality rates. Structural and energy calculation analysis of the V197M amino acid change showed that it destabilizes the TMPRSS2 protein, possibly negatively affecting its ACE2 and viral spike protein processing.

Paired comparisons of mutational profiles before and after brachytherapy in asian uveal melanoma patients

Uveal melanoma(UM) is the most common primary intraocular malignancy in adults. However, the incidence of UM in Asia is 10 to 20 times less than in Western populations. Therefore, for the first time, we report our whole exome sequencing (WES) data analysis to discover differences in the molecular features of Asian and Western UM, and to determine the disparities between the primary tumor before brachytherapy and enucleated samples after brachytherapy. WES of 19 samples (13 primary tumors, 5 enucleation samples after brachytherapy, and 1 liver metastasis) from 13 patients diagnosed with UM and treated between 2007 and 2019 at the Yonsei University Health System (YUHS) were analyzed using bioinformatics pipelines. We identified significantly altered genes in Asian UM and changes in mutational profiles before and after brachytherapy using various algorithms. GNAQ, BAP1, GNA11, SF3B1 and CYSLTR2 were significantly mutated in Asian UM, which is similar that reported frequently in previous Western-based UM studies. There were also similar copy number alterations (M3, 1p loss, 6p gain, 8q gain) in both groups. In paired comparisons of the same patients, DICER1 and LRP1B were distinctly mutated only in tumor samples obtained after brachytherapy using rare-variant association tests (P = 0.01, 0.01, respectively). The mutational profiles of Asian UM were generally similar to the data from previous Western-based studies. DICER1 and LRP1B were newly mutated genes with statistical significance in the regrowth samples after brachytherapy compared to the primary tumors, which may be related to resistance to brachytherapy.

Discovery of genomic variation across a generation

Over the past 30 years (the timespan of a generation), advances in genomics technologies have revealed tremendous and unexpected variation in the human genome and have provided increasingly accurate answers to long-standing questions of how much genetic variation exists in human populations and to what degree the DNA complement changes between parents and offspring. Tracking the characteristics of these inherited and spontaneous (or de novo) variations has been the basis of the study of human genetic disease. From genome-wide microarray and next-generation sequencing scans, we now know that each human genome contains over 3 million single nucleotide variants when compared with the ~ 3 billion base pairs in the human reference genome, along with roughly an order of magnitude more DNA—approximately 30 megabase pairs (Mb)—being ‘structurally variable’, mostly in the form of indels and copy number changes. Additional large-scale variations include balanced inversions (average of 18 Mb) and complex, difficult-to-resolve alterations. Collectively, ~1% of an individual’s genome will differ from the human reference sequence. When comparing across a generation, fewer than 100 new genetic variants are typically detected in the euchromatic portion of a child’s genome. Driven by increasingly higher-resolution and higher-throughput sequencing technologies, newer and more accurate databases of genetic variation (for instance, more comprehensive structural variation data and phasing of combinations of variants along chromosomes) of worldwide populations will emerge to underpin the next era of discovery in human molecular genetics.