Whole genome sequencing of 35 individuals provides insights into the genetic architecture of Korean population

Modern Japanese ancestry-derived variants reveal the formation process of the current Japanese regional gradations

Modern Japanese people have two major ancestral populations: indigenous Jomon hunter-gatherers and continental East Asian farmers. To determine the formation process of the current Japanese population, we developed a detection method for variants derived from ancestral populations using a summary statistic, the ancestry marker index (AMI). We applied AMI to modern Japanese population samples and identified 208,648 single nucleotide polymorphisms (SNPs) that were likely derived from the Jomon people (Jomon-derived variants). Analysis of Jomon-derived variants in 10,842 modern Japanese individuals recruited from all over Japan revealed that the admixture proportions of the Jomon people varied between prefectures, probably owing to the prehistoric population size difference. The estimated allele frequencies of genome-wide SNPs in the ancestral populations of the modern Japanese suggested their adaptive phenotypic characteristics to their respective livelihoods. Based on our findings, we propose a formation model for the genotypic and phenotypic gradations of the current Japanese archipelago populations.

Multi-ancestry genome-wide association analyses identify novel genetic mechanisms in rheumatoid arthritis

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10-8), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA.

A data set of variants derived from 1455 clinical and research exomes is efficient in variant prioritization for early-onset monogenic disorders in Indians

Given the genomic uniqueness, a local data set is most desired for Indians, who are underrepresented in existing public databases. We hypothesize patients with rare monogenic disorders and their family members can provide a reliable source of common variants in the population. Exome sequencing (ES) data from families with rare Mendelian disorders was aggregated from five centers in India. The dataset was refined by excluding related individuals and removing the disease-causing variants (refined cohort). The efficiency of these data sets was assessed in a new set of 50 exomes against gnomAD and GenomeAsia. Our original cohort comprised 1455 individuals from 1203 families. The refined cohort had 836 unrelated individuals that retained 1,251,064 variants with 181,125 population-specific and 489,618 common variants. The allele frequencies from our cohort helped to define 97,609 rare variants in gnomAD and 44,520 rare variants in GenomeAsia as common variants in our population. Our variant dataset provided an additional 1.7% and 0.1% efficiency for prioritizing heterozygous and homozygous variants respectively for rare monogenic disorders. We observed additional 19 genes/human knockouts. We list carrier frequency for 142 recessive disorders. This is a large and useful resource of exonic variants for Indians. Despite limitations, datasets from patients are efficient tools for variant prioritization in a resource-limited setting.

Gastric Microbiome Diversities in Gastric Cancer Patients from Europe and Asia Mimic the Human Population Structure and Are Partly Driven by Microbiome Quantitative Trait Loci

The human gastrointestinal tract harbors approximately 100 trillion microorganisms with different microbial compositions across geographic locations. In this work, we used RNASeq data from stomach samples of non-disease (164 individuals from European ancestry) and gastric cancer patients (137 from Europe and Asia) from public databases. Although these data were intended to characterize the human expression profiles, they allowed for a reliable inference of the microbiome composition, as confirmed from measures such as the genus coverage, richness and evenness. The microbiome diversity (weighted UniFrac distances) in gastric cancer mimics host diversity across the world, with European gastric microbiome profiles clustering together, distinct from Asian ones. Despite the confirmed loss of microbiome diversity from a healthy status to a cancer status, the structured profile was still recognized in the disease condition. In concordance with the parallel host-bacteria population structure, we found 16 human loci (non-synonymous variants) in the European-descendent cohorts that were significantly associated with specific genera abundance. These microbiome quantitative trait loci display heterogeneity between population groups, being mainly linked to the immune system or cellular features that may play a role in enabling microbe colonization and inflammation.

Association of DOCK8, IL17RA, and KLK12 Polymorphisms with Atopic Dermatitis in Koreans

Background

Early-onset and severe atopic dermatitis (AD) in patients increase the probability of the development of allergic rhinitis or asthma. Treatment and prevention strategies in infants and young children with AD are targeted toward treating the symptoms, restoring skin barrier functions, and reducing the absorption of environmental allergens in an attempt to attenuate or block the onset of asthma and food allergy.

Objective

Given that the initiating events in AD remain poorly understood, identifying those at risk and implementing strategies to prevent AD is necessary.

Methods

Whole-exome sequencing (WES) was performed in a 43 control group and a disease group with 20 AD patients without atopic march (AM) and 20 with AM. Sanger sequencing was carried out to validate found variants in cohorts.

Results

DOCK8, IL17RA, and KLK12 single-nucleotide polymorphisms were identified by WES as missense mutations: c.1289C>A, p.P97T (rs529208); c.1685C>A, p.P562G (rs12484684); and c.457+27>C, rs3745540, respectively. A case-control study show that total immunoglobulin E (IgE) level was significantly increased in the AA genotype of DOCK8 compared to the CA genotype in allergic patients. The rs12484684 of IL17RA increased risk of adult-onset AD (odds ratio: 1.63) compared to the control for (A) allele frequency. AD and AM Patients with the IL17RA CA genotype also had elevated IgE levels. rs3745540 of KLK12 was associated with AD in dominant model (odds ratio: 2.86).

Conclusion

DOCK8 (rs529208), IL17RA (rs12484684), and KLK12 (rs3745540), were identified using a new WES filtering method. the result suggests that polymorphism of DOCK8 and IL17RA might be related to increase the total IgE level.

A Vietnamese human genetic variation database

Large scale human genome projects have created tremendous human genome databases for some well-studied populations. Vietnam has about 95 million people (the 14th largest country by population in the world) of which more than 86% are Kinh people. To date, genetic studies for Vietnamese people mostly rely on genetic information from other populations. Building a Vietnamese human genetic variation database is a must for properly interpreting Vietnamese genetic variants. To this end, we sequenced 105 whole genomes and 200 whole exomes of 305 unrelated Kinh Vietnamese (KHV) people. We also included 101 other previously published KHV genomes to build a Vietnamese human genetic variation database of 406 KHV people. The KHV database contains 24.81 million variants (22.47 million single nucleotide polymorphisms (SNPs) and 2.34 million indels) of which 0.71 million variants are novel. It includes more than 99.3% of variants with a frequency of >1% in the KHV population. Noticeably, the KHV database revealed 107 variants reported in the human genome mutation database as pathological mutations with a frequency above 1% in the KHV population. The KHV database (available at https://genomes.vn) would be beneficial for genetic studies and medical applications not only for the Vietnamese population but also for other closely related populations.

Analysis of whole Y-chromosome sequences reveals the Japanese population history in the Jomon period

The Jomon and the Yayoi are considered to be the two major ancestral populations of the modern mainland Japanese. The Jomon people, who inhabited mainland Japan, admixed with Yayoi immigrants from the Asian continent. To investigate the population history in the Jomon period (14,500–2,300 years before present [YBP]), we analyzed whole Y-chromosome sequences of 345 Japanese males living in mainland Japan. A phylogenetic analysis of East Asian Y chromosomes identified a major clade (35.4% of mainland Japanese) consisting of only Japanese Y chromosomes, which seem to have originated from indigenous Jomon people. A Monte Carlo simulation indicated that ~70% of Jomon males had Y chromosomes in this clade. The Bayesian skyline plots of 122 Japanese Y chromosomes in the clade detected a marked decrease followed by a subsequent increase in the male population size from around the end of the Jomon period to the beginning of the Yayoi period (2,300 YBP). The colder climate in the Late to Final Jomon period may have resulted in critical shortages of food for the Jomon people, who were hunter-gatherers, and the rice farming introduced by Yayoi immigrants may have helped the population size of the Jomon people to recover.

Population-scale genomics-Enabling precision public health

The current excitement for affordable genomics technologies and national precision medicine initiatives marks a turning point in worldwide healthcare practices. The last decade of global population sequencing efforts has defined the enormous extent of genetic variation in the human population resulting in insights into differential disease burden and response to therapy within and between populations. Population-scale pharmacogenomics helps to provide insights into the choice of optimal therapies and an opportunity to estimate, predict and minimize adverse events. Such an approach can potentially empower countries to formulate national selection and dosing policies for therapeutic agents thereby promoting public health with precision. We review the breadth and depth of worldwide population-scale sequencing efforts and its implications for the implementation of clinical pharmacogenetics toward making precision medicine a reality.

Detection of novel germline mutations in six breast cancer predisposition genes by targeted next-generation sequencing

In this study, a customized amplicon-based target sequencing panel was designed to enrich the whole exon regions of six genes associated with the risk of breast cancer. Targeted next-generation sequencing (NGS) was performed for 146 breast cancer patients (BC), 71 healthy women with a family history of breast cancer (high risk), and 55 healthy women without a family history of cancer (control). Sixteen possible disease-causing mutations on four genes were identified in 20 samples. The percentages of possible disease-causing mutation carriers in the BC group (8.9%) and in the high-risk group (8.5%) were higher than that in the control group (1.8%). The BRCA1 possible disease-causing mutation group had a higher prevalence in family history and triple-negative breast cancer, while the BRCA2 possible disease-causing mutation group was younger and more likely to develop axillary lymph node metastasis (P < 0.05). Among the 146 patients, 47 with a family history of breast cancer were also sequenced with another 14 moderate-risk genes. Three additional possible disease-causing mutations were found on PALB2, CHEK2, and PMS2 genes, respectively. The results demonstrate that the six-gene targeted NGS panel may provide an approach to assess the genetic risk of breast cancer and predict the clinical prognosis of breast cancer patients.

Population genetics-informed meta-analysis in seven genes associated with risk to dengue fever disease

Population genetics theory predicted that rare frequent markers would be the main contributors for heritability of complex diseases, but meta-analyses of genome-wide association studies are revealing otherwise common markers, present in all population groups, as the identified candidate genes. In this work, we applied a population-genetics informed meta-analysis to 10 markers located in seven genes said to be associated with dengue fever disease. Seven markers (in PLCE1, CD32, CD209, OAS1 and OAS3 genes) have high-frequency and the other three (in MICB and TNFA genes) have intermediate frequency. Most of these markers have high discriminatory power between population groups, but their frequencies follow the rules of genetic drift, and seem to have not been under strong selective pressure. There was a good agreement in directional consistency across trans-ethnic association signals, in East Asian and Latin American cohorts, with heterogeneity generated by randomness between studies and especially by low sample sizes. This led to confirm the following significant associations: with DF, odds ratio of 0.67 for TNFA-rs1800629-A; with DHF, 0.82 for CD32-rs1801274-G; with DSS, 0.55 for OAS3-rs2285933-G, 0.80 for PLCE1-rs2274223-G and 1.32 for MICB-rs3132468-C. The overall genetic risks confirmed sub-Saharan African populations and descendants as the best protected against the severer forms of the disease, while Southeast and Northeast Asians are the least protected ones. European and close neighbours are the best protected against dengue fever, while, again, Southeast and Northeast Asians are the least protected ones. These risk scores provide important predictive information for the largely naïve European and North American regions, as well as for Africa where misdiagnosis with other hemorrhagic diseases is of concern.

KoVariome: Korean National Standard Reference Variome database of whole genomes with comprehensive SNV, indel, CNV, and SV analyses

High-coverage whole-genome sequencing data of a single ethnicity can provide a useful catalogue of population-specific genetic variations, and provides a critical resource that can be used to more accurately identify pathogenic genetic variants. We report a comprehensive analysis of the Korean population, and present the Korean National Standard Reference Variome (KoVariome). As a part of the Korean Personal Genome Project (KPGP), we constructed the KoVariome database using 5.5 terabases of whole genome sequence data from 50 healthy Korean individuals in order to characterize the benign ethnicity-relevant genetic variation present in the Korean population. In total, KoVariome includes 12.7M single-nucleotide variants (SNVs), 1.7M short insertions and deletions (indels), 4K structural variations (SVs), and 3.6K copy number variations (CNVs). Among them, 2.4M (19%) SNVs and 0.4M (24%) indels were identified as novel. We also discovered selective enrichment of 3.8M SNVs and 0.5M indels in Korean individuals, which were used to filter out 1,271 coding-SNVs not originally removed from the 1,000 Genomes Project when prioritizing disease-causing variants. KoVariome health records were used to identify novel disease-causing variants in the Korean population, demonstrating the value of high-quality ethnic variation databases for the accurate interpretation of individual genomes and the precise characterization of genetic variations.

Association of PD-L1 Expression with Tumor-Infiltrating Immune Cells and Mutation Burden in High-Grade Neuroendocrine Carcinoma of the Lung

INTRODUCTION

The immune microenvironment of high-grade neuroendocrine carcinoma of the lung, including programmed death ligand 1 (PD-L1) expression, has not been well characterized.

METHODS

On the basis of immunohistochemistry (IHC) results, PD-L1 expression on tumor cells (TCs) and tumor-infiltrating immune cells (ICs) was scored as follows: TC0 and IC0 were defined as PD-L1 expression less than 1%, TC1 and IC1 as at least 1% but less than 10%, TC2 and IC2 as 10% or more but less than 50%, and TC3 and IC3 as 50% or more. Phosphatase and tensin homolog (PTEN) IHC was scored as either lost or retained expression. The Ion AmpliSeq Comprehensive Cancer Panel (ThermoFisher Scientific, Waltham, MA) was used to identify mutations in all coding exons of 409 cancer-related genes.

RESULTS

A total of 192 patients with large cell neuroendocrine carcinoma (LCNEC) (n = 72) and SCLC (n = 120) were studied. The prevalence of PD-L1 expression on TCs was 15.1% (29 of 192). IC infiltration and PD-L1 expression on ICs were observed in 34.4% of patients (66 of 192) and 31.3% of patients (60 of 192), respectively. The prevalence of IC infiltration and PD-L1 expression on IC were more strongly correlated with LCNEC than with SCLC (57.6% versus 23.3%, p < 0.01; 45.8% versus 22.5%, p < 0.01) and high nonsynonymous mutations (p = 0.05 and .04). PTEN loss was found in 9.5% of patients (18 of 189) and showed no correlation with PD-L1 expression. Progression-free survival was better in patients with IC infiltration than in those without IC infiltration (median 11.3 versus 6.8 months [p < 0.01]) and in patients with PD-L1 expression of IC1/2/3 than in those with expression of IC0 (median 11.3 versus 7.0 months [p = 0.03]).

CONCLUSION

These findings suggest that the PD-1/PD-L1 pathway is activated in the microenvironment of pulmonary high-grade neuroendocrine carcinoma and correlated with a higher mutation burden.

Demographic history and biologically relevant genetic variation of Native Mexicans inferred from whole-genome sequencing

Understanding the genetic structure of Native American populations is important to clarify their diversity, demographic history, and to identify genetic factors relevant for biomedical traits. Here, we show a demographic history reconstruction from 12 Native American whole genomes belonging to six distinct ethnic groups representing the three main described genetic clusters of Mexico (Northern, Southern, and Maya). Effective population size estimates of all Native American groups remained below 2,000 individuals for up to 10,000 years ago. The proportion of missense variants predicted as damaging is higher for undescribed (~ 30%) than for previously reported variants (~ 15%). Several variants previously associated with biological traits are highly frequent in the Native American genomes. These findings suggest that the demographic and adaptive processes that occurred in these groups shaped their genetic architecture and could have implications in biological processes of the Native Americans and Mestizos of today.

People of Mexico have diverse historical and genetic background. Here, Romero-Hidalgo and colleagues sequence whole genomes of Native Americans of Mexico, and show demographic history and genetic variation shared among subgroups of Native Americans.

Highly Variable Genomic Landscape of Endogenous Retroviruses in the C57BL/6J Inbred Strain, Depending on Individual Mouse, Gender, Organ Type, and Organ Location

Transposable repetitive elements, named the "TREome," represent ~40% of the mouse genome. We postulate that the germ line genome undergoes temporal and spatial diversification into somatic genomes in conjunction with the TREome activity. C57BL/6J inbred mice were subjected to genomic landscape analyses using a TREome probe from murine leukemia virus-type endogenous retroviruses (MLV-ERVs). None shared the same MLV-ERV landscape within each comparison group: (1) sperm and 18 tissues from one mouse, (2) six brain compartments from two females, (3) spleen and thymus samples from four age groups, (4) three spatial tissue sets from two females, and (5) kidney and liver samples from three females and three males. Interestingly, males had more genomic MLV-ERV copies than females; moreover, only in the males, the kidneys had higher MLV-ERV copies than the livers. Perhaps, the mouse-, gender-, and tissue/cell-dependent MLV-ERV landscapes are linked to the individual-specific and dynamic phenotypes of the C57BL/6J inbred population.

An ethnically relevant consensus Korean reference genome is a step towards personal reference genomes

Human genomes are routinely compared against a universal reference. However, this strategy could miss population-specific and personal genomic variations, which may be detected more efficiently using an ethnically relevant or personal reference. Here we report a hybrid assembly of a Korean reference genome (KOREF) for constructing personal and ethnic references by combining sequencing and mapping methods. We also build its consensus variome reference, providing information on millions of variants from 40 additional ethnically homogeneous genomes from the Korean Personal Genome Project. We find that the ethnically relevant consensus reference can be beneficial for efficient variant detection. Systematic comparison of human assemblies shows the importance of assembly quality, suggesting the necessity of new technologies to comprehensively map ethnic and personal genomic structure variations. In the era of large-scale population genome projects, the leveraging of ethnicity-specific genome assemblies as well as the human reference genome will accelerate mapping all human genome diversity.

Comparing genetic variants detected in the 1000 genomes project with SNPs determined by the International HapMap Consortium

Single-nucleotide polymorphisms (SNPs) determined based on SNP arrays from the international HapMap consortium (HapMap) and the genetic variants detected in the 1000 genomes project (1KGP) can serve as two references for genomewide association studies (GWAS). We conducted comparative analyses to provide a means for assessing concerns regarding SNP array-based GWAS findings as well as for realistically bounding expectations for next generation sequencing (NGS)-based GWAS. We calculated and compared base composition, transitions to transversions ratio, minor allele frequency and heterozygous rate for SNPs from HapMap and 1KGP for the 622 common individuals. We analysed the genotype discordance between HapMap and 1KGP to assess consistency in the SNPs from the two references. In 1KGP, 90.58% of 36,817,799 SNPs detected were not measured in HapMap. More SNPs with minor allele frequencies less than 0.01 were found in 1KGP than HapMap. The two references have low disc ordance (generally smaller than 0.02) in genotypes of common SNPs, with most discordance from heterozygous SNPs. Our study demonstrated that SNP array-based GWAS findings were reliable and useful, although only a small portion of genetic variances were explained. NGS can detect not only common but also rare variants, supporting the expectation that NGS-based GWAS will be able to incorporate a much larger portion of genetic variance than SNP arrays-based GWAS.

Mutational Analysis of Extranodal NK/T-Cell Lymphoma Using Targeted Sequencing with a Comprehensive Cancer Panel

Extranodal natural killer (NK)/T-cell lymphoma, nasal type (NKTCL), is a malignant disorder of cytotoxic lymphocytes of NK or T cells. It is an aggressive neoplasm with a very poor prognosis. Although extranodal NKTCL reportedly has a strong association with Epstein-Barr virus, the molecular pathogenesis of NKTCL has been unexplored. The recent technological advancements in next-generation sequencing (NGS) have made DNA sequencing cost- and time-effective, with more reliable results. Using the Ion Proton Comprehensive Cancer Panel, we sequenced 409 cancer-related genes to identify somatic mutations in five NKTCL tissue samples. The sequencing analysis detected 25 mutations in 21 genes. Among them, KMT2D, a histone modification-related gene, was the most frequently mutated gene (four of the five cases). This result was consistent with recent NGS studies that have suggested KMT2D as a novel driver gene in NKTCL. Mutations were also found in ARID1A, a chromatin remodeling gene, and TP53, which also recurred in recent NGS studies. We also found mutations in 18 novel candidate genes, with molecular functions that were potentially implicated in cancer development. We suggest that these genes may result in multiple oncogenic events and may be used as potential bio-markers of NKTCL in the future.

Genome-wide association study of antidepressant response: involvement of the inorganic cation transmembrane transporter activity pathway

BACKGROUND

Genome-wide association studies (GWAS) represent the current frontier in pharmacogenomics. Thousands of subjects of Caucasian ancestry have been included in previous GWAS investigating antidepressant response. GWAS focused on this phenotype are lacking in Asian populations.

METHODS

A sample of 109 major depressive disorder (MDD) patients of Korean origin in antidepressant treatment was collected. Phenotypes were response and remission according to the Hamilton Rating Scale for Depression (HRSD). Genome-wide genotyping was performed using the Illumina Human Omni2.5-8 platform. The same phenotypes were used in the STAR*D level 1 (n = 1677) for independent replication. In order to corroborate findings and increase the comparability between the two datasets, three levels of analysis (SNPs, genes and pathways) were carried out. Bonferroni correction, permutations, and replication across samples were used to reduce the risk of false positives.

RESULTS

Among the genes replicated across the two samples (permutated p < 0.05 in both of them), CTNNA3 appeared promising. The inorganic cation transmembrane transporter activity pathway (GO:0022890) was associated with antidepressant response in both samples (p = 2.9e-5 and p = 0.001 in the Korean and STAR*D samples, respectively) and this pathway included CACNA1A, CACNA1C, and CACNB2 genes.

CONCLUSIONS

The present study supported the involvement of genes coding for subunits of L-type voltage-gated calcium channel in antidepressant efficacy across different ethnicities but replication of findings is required before any definitive statement.

Identification of genomic mutations associated with clinical outcomes of induction chemotherapy in patients with head and neck squamous cell carcinoma

PURPOSE

We performed deep sequencing of target genes in head and neck squamous cell carcinoma (HNSCC) tumors to identify somatic mutations that are associated with induction chemotherapy (IC) response.

METHODS

Patients who were diagnosed with HNSCC were retrospectively identified. Patients who were treated with IC were divided into two groups: good responders and poor responders by tumor response and progression-free survival. Targeted gene sequencing for 2404 somatic mutations of 44 genes was performed on HNSCC tissues. Mutations with total coverage of <500 were excluded, and the cutoff for altered allele frequency was >10 %.

RESULTS

Of the 71 patients, 45 were treated upfront with IC. Mean total coverage was 1941 per locus, and 42.2 % of tumors had TP53 mutations. Thirty-three mutations in TP53, NOTCH3, FGFR2, FGFR3, ATM, EGFR, MET, PTEN, FBXW7, SYNE1, and SUFU were frequently altered in poor responders. Among the patients who were treated with IC, those with unfavorable genomic profiles had significantly poorer overall survival than those without unfavorable genomic profiles (hazard ratio 6.45, 95 % confidence interval 2.07-20.10, P < 0.001).

CONCLUSIONS

Comprehensive analysis of mutation frequencies identified unfavorable genomic profiles, and the patients without unfavorable genomic profiles can obtain clinical benefits from IC in patients with HNSCC.

Human genetic variation database, a reference database of genetic variations in the Japanese population

Whole-genome and -exome resequencing using next-generation sequencers is a powerful approach for identifying genomic variations that are associated with diseases. However, systematic strategies for prioritizing causative variants from many candidates to explain the disease phenotype are still far from being established, because the population-specific frequency spectrum of genetic variation has not been characterized. Here, we have collected exomic genetic variation from 1208 Japanese individuals through a collaborative effort, and aggregated the data into a prevailing catalog. In total, we identified 156 622 previously unreported variants. The allele frequencies for the majority (88.8%) were lower than 0.5% in allele frequency and predicted to be functionally deleterious. In addition, we have constructed a Japanese-specific major allele reference genome by which the number of unique mapping of the short reads in our data has increased 0.045% on average. Our results illustrate the importance of constructing an ethnicity-specific reference genome for identifying rare variants. All the collected data were centralized to a newly developed database to serve as useful resources for exploring pathogenic variations. Public access to the database is available at http://www.genome.med.kyoto-u.ac.jp/SnpDB/.

Genetics of Aldosterone-Producing Adenoma in Korean Patients

Objectives

Recently, somatic mutations in KCNJ5, ATP1A1, ATP2B3, and CACNA1D genes were found to be associated with the pathogenesis of aldosterone-producing adenoma (APA). This study aimed to investigate the prevalence of somatic mutations in KCNJ5, ATP1A1, ATP2B3, and CACNA1D and examine the correlations between these mutations and the clinical and biochemical characteristics in Korean patients with APA.

Methods

We performed targeted gene sequencing in 66 patients with APA to detect somatic mutations in these genes.

Results

Somatic KCNJ5 mutations were found in 47 (71.2%) of the 66 patients with APA (31 cases of p.G151R and 16 cases of p.L168R); these two mutations were mutually exclusive. Somatic mutations in the ATP1A1, ATP2B3, and CACNA1D genes were not observed. Somatic KCNJ5 mutations were more prevalent in female patients (66% versus 36.8%, respectively; P = 0.030). Moreover, patients with KCNJ5 mutations comprised a significantly higher proportion of patients younger than 35 years of age (19.1% versus 0%, respectively; P = 0.040). There were no significant differences in pre-operative blood pressure, plasma aldosterone, serum potassium, lateralization index, and adenoma size according to mutational status. Patients with KCNJ5 mutations were less likely to need antihypertensive medications after adrenalectomy compared with those without mutation (36.2% versus 63.2%; P = 0.045).

Conclusions

The present study demonstrated the high prevalence of somatic KCNJ5 mutations in Korean patients with APA. Carriers of somatic KCNJ5 mutations were more likely to be female. Early diagnosis and better therapeutic outcomes were associated with somatic KCNJ5 mutations in APA.

Alignment of Short Reads: A Crucial Step for Application of Next-Generation Sequencing Data in Precision Medicine

Precision medicine or personalized medicine has been proposed as a modernized and promising medical strategy. Genetic variants of patients are the key information for implementation of precision medicine. Next-generation sequencing (NGS) is an emerging technology for deciphering genetic variants. Alignment of raw reads to a reference genome is one of the key steps in NGS data analysis. Many algorithms have been developed for alignment of short read sequences since 2008. Users have to make a decision on which alignment algorithm to use in their studies. Selection of the right alignment algorithm determines not only the alignment algorithm but also the set of suitable parameters to be used by the algorithm. Understanding these algorithms helps in selecting the appropriate alignment algorithm for different applications in precision medicine. Here, we review current available algorithms and their major strategies such as seed-and-extend and q-gram filter. We also discuss the challenges in current alignment algorithms, including alignment in multiple repeated regions, long reads alignment and alignment facilitated with known genetic variants.