Whole Exome Sequencing in Patients with Phenotypically Associated Familial Intracranial Aneurysm

Objective@#Familial intracranial aneurysms (FIAs) are found in approximately 6%–20% of patients with intracranial aneurysms (IAs), suggesting that genetic predisposition likely plays a role in its pathogenesis. The aim of this study was to identify possible IA-associated variants using whole exome sequencing (WES) in selected Korean families with FIA. @*Materials and Methods@#Among the 26 families in our institutional database with two or more IA-affected first-degree relatives, three families that were genetically enriched (multiple, early onset, or common site involvement within the families) for IA were selected for WES. Filtering strategies, including a family-based approach and knowledge-based prioritization, were applied to derive possible IA-associated variants from the families. A chromosomal microarray was performed to detect relatively large chromosomal abnormalities. @*Results@#Thirteen individuals from the three families were sequenced, of whom seven had IAs. We noted three rare, potentially deleterious variants (PLOD3 c.1315G>A, NTM c.968C>T, and CHST14 c.58C>T), which are the most promising candidates among the 11 potential IA-associated variants considering gene-phenotype relationships, gene function, co-segregation, and variant pathogenicity. Microarray analysis did not reveal any significant copy number variants in the families. @*Conclusion@#Using WES, we found that rare, potentially deleterious variants in PLOD3, NTM, and CHST14 genes are likely responsible for the subsets of FIAs in a cohort of Korean families.

Identification of rare coding variants associated with Kawasaki disease by whole exome sequencing

Kawasaki disease (KD) is an acute pediatric vasculitis that affects genetically susceptible infants and children. To identify coding variants that influence susceptibility to KD, we conducted whole exome sequencing of 159 patients with KD and 902 controls, and performed a replication study in an independent 586 cases and 732 controls. We identified five rare coding variants in five genes (FCRLA, PTGER4, IL17F, CARD11, and SIGLEC10) associated with KD (odds ratio [OR], 1.18–4.41; p = 0.0027–0.031). We also performed association analysis in 26 KD patients with coronary artery aneurysms (CAAs; diameter > 5 mm) and 124 patients without CAAs (diameter < 3 mm), and identified another five rare coding variants in five genes (FGFR4, IL31RA, FNDC1, MMP8, and FOXN1), which may be associated with CAA (OR, 3.89–37.3; p = 0.0058–0.0261). These results provide insights into new candidate genes and genetic variants potentially involved in the development of KD and CAA.

Rare Neurovascular Diseases in Korea: Classification and Related Genetic Variants

Rare neurovascular diseases (RNVDs) have not been well-recognized in Korea. They involve the central nervous system and greatly affect the patients’ lives. However, these diseases are difficult to diagnose and treat due to their rarity and incurability. We established a list of RNVDs by referring to the previous literature and databases worldwide to better understand the diseases and their current management status. We categorized 68 RNVDs based on their pathophysiology and clinical manifestations and estimated the prevalence of each disease in Korea. Recent advances in genetic, molecular, and developmental research have enabled further understanding of these RNVDs. Herein, we review each disease, while considering its classification based on updated pathologic mechanisms, and discuss the management status of RNVD in Korea.

IgA Levels Are Associated with Coronary Artery Lesions in Kawasaki Disease

Background and Objectives@#Kawasaki disease (KD) is an acute systemic vasculitis that affects the coronary arteries. Abnormal immune reactions are thought to contribute to disease pathogenesis. The effect of immunoglobulin (Ig) isotype (IgG, IgA, IgM, and IgE) on inflammatory data and clinical outcomes of patients with KD was examined. @*Methods@#Ig levels in 241 patients with KD were measured during the acute, subacute, convalescent, and normal phases of the disease. @*Results@#Compared with reference Ig values, IgG, IgA, and IgM levels were significantly higher in the subacute phase, while IgE levels were elevated in 73.9% (178/241) of patients with KD in all clinical phases. However, high IgE levels were not associated with clinical outcomes, including intravenous immunoglobulin unresponsiveness and coronary artery lesions (CALs).Significantly more CALs were observed in the high IgA group than in the normal IgA group (44.7% vs. 20.8%, respectively; p<0.01). In addition, IgA levels in the acute phase (p=0.038) were 2.2-fold higher, and those in the subacute phase were 1.7-fold higher (p <0.001), in the CAL group than in the non-CAL group. IgA concentrations increased along with the size of the coronary artery aneurysm (p <0.001). Furthermore, there was a strong correlation between IgA levels and CAL size (r=0.435, p<0.001), with a high odds ratio of 2.58 (p=0.022). @*Conclusions@#High IgA levels in patients with KD are prognostic for the risk of CALs.

Rare Neurovascular Diseases in Korea: Classification and Related Genetic Variants

Rare neurovascular diseases (RNVDs) have not been well-recognized in Korea. They involve the central nervous system and greatly affect the patients’ lives. However, these diseases are difficult to diagnose and treat due to their rarity and incurability. We established a list of RNVDs by referring to the previous literature and databases worldwide to better understand the diseases and their current management status. We categorized 68 RNVDs based on their pathophysiology and clinical manifestations and estimated the prevalence of each disease in Korea. Recent advances in genetic, molecular, and developmental research have enabled further understanding of these RNVDs. Herein, we review each disease, while considering its classification based on updated pathologic mechanisms, and discuss the management status of RNVD in Korea.

IgA Levels Are Associated with Coronary Artery Lesions in Kawasaki Disease

Background and Objectives@#Kawasaki disease (KD) is an acute systemic vasculitis that affects the coronary arteries. Abnormal immune reactions are thought to contribute to disease pathogenesis. The effect of immunoglobulin (Ig) isotype (IgG, IgA, IgM, and IgE) on inflammatory data and clinical outcomes of patients with KD was examined. @*Methods@#Ig levels in 241 patients with KD were measured during the acute, subacute, convalescent, and normal phases of the disease. @*Results@#Compared with reference Ig values, IgG, IgA, and IgM levels were significantly higher in the subacute phase, while IgE levels were elevated in 73.9% (178/241) of patients with KD in all clinical phases. However, high IgE levels were not associated with clinical outcomes, including intravenous immunoglobulin unresponsiveness and coronary artery lesions (CALs).Significantly more CALs were observed in the high IgA group than in the normal IgA group (44.7% vs. 20.8%, respectively; p<0.01). In addition, IgA levels in the acute phase (p=0.038) were 2.2-fold higher, and those in the subacute phase were 1.7-fold higher (p <0.001), in the CAL group than in the non-CAL group. IgA concentrations increased along with the size of the coronary artery aneurysm (p <0.001). Furthermore, there was a strong correlation between IgA levels and CAL size (r=0.435, p<0.001), with a high odds ratio of 2.58 (p=0.022). @*Conclusions@#High IgA levels in patients with KD are prognostic for the risk of CALs.

Assessment of the Clinical Heterogeneity of Kawasaki Disease Using Genetic Variants of BLK and FCGR2A

BACKGROUND AND OBJECTIVES@#Patients with Kawasaki disease (KD) are clinically heterogeneous because its diagnosis is based solely on clinical observation and there are no definitive biomarkers. We dissected the clinical heterogeneity of KD patients using the KD-associated genetic variants.@*METHODS@#We performed a genetic association analysis in several KD subgroups categorized by clinical characteristics using the KD-associated variants of the B lymphoid tyrosine kinase (BLK; rs6993775) and Fc gamma receptor II a (FCGR2A; rs1801274) in a large number of case (n=1,011) and control (n=4,533) samples.@*RESULTS@#BLK and FCGR2A were very significantly associated with KD in Korean KD patients (odds ratio [OR],1.48; p=4.63×10⁻¹¹ for BLK, and OR, 1.26; p=1.42×10⁻⁴ for FCGR2A). However, in KD subgroup analysis, we found that neither BLK nor FCGR2A were associated with either incomplete Kawasaki disease (iKD) type patients or those older than 5 years of age (p>0.2), suggesting that patients with iKD or those older than 5 years of age are a unique subgroup of KD. In genetic association analysis after excluding iKD patients and those older than 5 years old, we found that BLK was associated with all KD subgroups, whereas FCGR2A was specifically associated with male KD patients younger than 1 year of age (OR, 2.22; p=2.35×10⁻⁵).@*CONCLUSIONS@#KD is a clinically and genetically heterogeneous disease. These findings will provide new insights into the clinical and genetic heterogeneity of KD.

Assessment of the Clinical Heterogeneity of Kawasaki Disease Using Genetic Variants of BLK and FCGR2A

BACKGROUND AND OBJECTIVES: Patients with Kawasaki disease (KD) are clinically heterogeneous because its diagnosis is based solely on clinical observation and there are no definitive biomarkers. We dissected the clinical heterogeneity of KD patients using the KD-associated genetic variants. METHODS: We performed a genetic association analysis in several KD subgroups categorized by clinical characteristics using the KD-associated variants of the B lymphoid tyrosine kinase (BLK; rs6993775) and Fc gamma receptor II a (FCGR2A; rs1801274) in a large number of case (n=1,011) and control (n=4,533) samples. RESULTS: BLK and FCGR2A were very significantly associated with KD in Korean KD patients (odds ratio [OR],1.48; p=4.63×10⁻¹¹ for BLK, and OR, 1.26; p=1.42×10⁻⁴ for FCGR2A). However, in KD subgroup analysis, we found that neither BLK nor FCGR2A were associated with either incomplete Kawasaki disease (iKD) type patients or those older than 5 years of age (p>0.2), suggesting that patients with iKD or those older than 5 years of age are a unique subgroup of KD. In genetic association analysis after excluding iKD patients and those older than 5 years old, we found that BLK was associated with all KD subgroups, whereas FCGR2A was specifically associated with male KD patients younger than 1 year of age (OR, 2.22; p=2.35×10⁻⁵). CONCLUSIONS: KD is a clinically and genetically heterogeneous disease. These findings will provide new insights into the clinical and genetic heterogeneity of KD.

Genome-Wide Analysis Reveals Four Novel Loci for Attention-Deficit Hyperactivity Disorder in Korean Youths

OBJECTIVES: The molecular mechanisms underlying attention-deficit hyperactivity disorder (ADHD) remain unclear. Therefore, this study aimed to identify the genetic susceptibility loci for ADHD in Korean children with ADHD. We performed a case-control and a family-based genome-wide association study (GWAS), as well as genome-wide quantitative trait locus (QTL) analyses, for two symptom traits. METHODS: A total of 135 subjects (71 cases and 64 controls), for the case-control analysis, and 54 subjects (27 probands and 27 unaffected siblings), for the family-based analysis, were included. RESULTS: The genome-wide QTL analysis identified four single nucleotide polymorphisms (SNPs) (rs7684645 near APELA, rs12538843 near YAE1D1 and POU6F2, rs11074258 near MCTP2, and rs34396552 near CIDEA) that were significantly associated with the number of inattention symptoms in ADHD. These SNPs showed possible association with ADHD in the family-based GWAS, and with hyperactivity-impulsivity in genome-wide QTL analyses. Moreover, association signals in the family-based QTL analysis for the number of inattention symptoms were clustered near genes IL10, IL19, SCL5A9, and SKINTL. CONCLUSION: We have identified four QTLs with genome-wide significance and several promising candidates that could potentially be associated with ADHD (CXCR4, UPF1, SETD5, NALCN-AS1, ERC1, SOX2-OT, FGFR2, ANO4, and TBL1XR1). Further replication studies with larger sample sizes are needed.

Identification of LEF1 as a Susceptibility Locus for Kawasaki Disease in Patients Younger than 6 Months of Age

Kawasaki disease (KD) is an acute febrile vasculitis predominately affecting infants and children. The dominant incidence age of KD is from 6 months to 5 years of age, and the incidence is unusual in those younger than 6 months and older than 5 years of age. We tried to identify genetic variants specifically associated with KD in patients younger than 6 months or older than 5 years of age. We performed an age-stratified genome-wide association study using the Illumina HumanOmni1-Quad BeadChip data (296 cases vs. 1,000 controls) and a replication study (1,360 cases vs. 3,553 controls) in the Korean population. Among 26 candidate single nucleotide polymorphisms (SNPs) tested in replication study, only a rare nonsynonymous SNP (rs4365796: c.1106C>T, p.Thr369Met) in the lymphoid enhancer binding factor 1 (LEF1) gene was very significantly associated with KD in patients younger than 6 months of age (odds ratio [OR], 3.07; p(combined) = 1.10 × 10⁻⁵), whereas no association of the same SNP was observed in any other age group of KD patients. The same SNP (rs4365796) in the LEF1 gene showed the same direction of risk effect in Japanese KD patients younger than 6 months of age, although the effect was not statistically significant (OR, 1.42; p = 0.397). This result indicates that the LEF1 gene may play an important role as a susceptibility gene specifically affecting KD patients younger than 6 months of age.

Evaluation of Digital PCR as a Technique for Monitoring Acute Rejection in Kidney Transplantation

Early detection and proper management of kidney rejection are crucial for the long-term health of a transplant recipient. Recipients are normally monitored by serum creatinine measurement and sometimes with graft biopsies. Donor-derived cell-free deoxyribonucleic acid (cfDNA) in the recipient's plasma and/or urine may be a better indicator of acute rejection. We evaluated digital PCR (dPCR) as a system for monitoring graft status using single nucleotide polymorphism (SNP)-based detection of donor DNA in plasma or urine. We compared the detection abilities of the QX200, RainDrop, and QuantStudio 3D dPCR systems. The QX200 was the most accurate and sensitive. Plasma and/or urine samples were isolated from 34 kidney recipients at multiple time points after transplantation, and analyzed by dPCR using the QX200. We found that donor DNA was almost undetectable in plasma DNA samples, whereas a high percentage of donor DNA was measured in urine DNA samples, indicating that urine is a good source of cfDNA for patient monitoring. We found that at least 24% of the highly polymorphic SNPs used to identify individuals could also identify donor cfDNA in transplant patient samples. Our results further showed that autosomal, sex-specific, and mitochondrial SNPs were suitable markers for identifying donor cfDNA. Finally, we found that donor-derived cfDNA measurement by dPCR was not sufficient to predict a patient's clinical condition. Our results indicate that donor-derived cfDNA is not an accurate predictor of kidney status in kidney transplant patients.

Is Whole Exome Sequencing Clinically Practical in the Management of Pediatric Crohn's Disease?

BACKGROUND/AIMS: The aim of this study was to identify the profile of rare variants associated with Crohn's disease (CD) using whole exome sequencing (WES) analysis of Korean children with CD and to evaluate whether genetic profiles could provide information during medical decision making. METHODS: DNA samples from 18 control individuals and 22 patients with infantile, very-early and early onset CD of severe phenotype were used for WES. Genes were filtered using panels of inflammatory bowel disease (IBD)-associated genes and genes of primary immunodeficiency (PID) and monogenic IBD. RESULTS: Eighty-one IBD-associated variants and 35 variants in PID genes were revealed by WES. The most frequently occurring variants were carried by nine (41%) and four (18.2%) CD probands and were ATG16L2 (rs11235604) and IL17REL (rs142430606), respectively. Twenty-four IBD-associated variants and 10 PID variants were predicted to be deleterious and were identified in the heterozygous state. However, their functions were unknown with the exception of a novel p.Q111X variant in XIAP (X chromosome) of a male proband. CONCLUSIONS: The presence of many rare variants of unknown significance limits the clinical applicability of WES for individual CD patients. However, WES in children may be beneficial for distinguishing CD secondary to PID.

JAK2 V617F, MPL, and CALR Mutations in Korean Patients with Essential Thrombocythemia and Primary Myelofibrosis

Mutations in the calreticulin gene, CALR, have recently been discovered in subsets of patients with essential thrombocythemia (ET) or primary myelofibrosis (PMF). We investigated Korean patients with ET and PMF to determine the prevalence, and clinical and laboratory correlations of CALR/JAK2/MPL mutations. Among 84 ET patients, CALR mutations were detected in 23 (27.4%) and were associated with higher platelet counts (P=0.006) and lower leukocyte counts (P=0.035) than the JAK2 V617F mutation. Among 50 PMF patients, CALR mutations were detected in 11 (22.0%) and were also associated with higher platelet counts (P=0.035) and trended to a lower rate of cytogenetic abnormalities (P=0.059) than the JAK2 V617F mutation. By multivariate analysis, triple-negative status was associated with shorter overall survival (HR, 7.0; 95% CI, 1.6-31.1, P=0.01) and leukemia-free survival (HR, 6.3; 95% CI, 1.8-22.0, P=0.004) in patients with PMF. The type 1 mutation was the most common (61.1%) type among all patients with CALR mutations, and tended toward statistical predominance in PMF patients. All 3 mutations were mutually exclusive and were never detected in patients with other myeloid neoplasms showing thrombocytosis. CALR mutations characterize a distinct group of Korean ET and PMF patients. Triple-negative PMF patients in particular have an unfavorable prognosis, which supports the idea that triple-negative PMF is a molecularly high-risk disease.

Consortium-Based Genetic Studies of Kawasaki Disease in Korea: Korean Kawasaki Disease Genetics Consortium

In order to perform large-scale genetic studies of Kawasaki disease (KD) in Korea, the Korean Kawasaki Disease Genetics Consortium (KKDGC) was formed in 2008 with 10 hospitals. Since the establishment of KKDGC, there has been a collection of clinical data from a total of 1198 patients, and approximately 5 mL of blood samples per patient (for genomic deoxyribonucleic acid and plasma isolation), using a standard clinical data collection form and a nation-wide networking system for blood sample pick-up. In the clinical risk factor analysis using the collected clinical data of 478 KD patients, it was found that incomplete KD type, intravenous immunoglobulin (IVIG) non-responsiveness, and long febrile days are major risk factors for coronary artery lesions development, whereas low serum albumin concentration is an independent risk factor for IVIG non-responsiveness. In addition, we identified a KD susceptibility locus at 1p31, a coronary artery aneurysm locus (KCNN2 gene), and the causal variant in the C-reactive protein (CRP) promoter region, as determining the increased CRP levels in KD patients, by means of genome-wide association studies. Currently, this consortium is continually collecting more clinical data and genomic samples to identify the clinical and genetic risk factors via a single nucleotide polymorphism chip and exome sequencing, as well as collaborating with several international KD genetics teams. The consortium-based approach for genetic studies of KD in Korea will be a very effective way to understand the unknown etiology and causal mechanism of KD, which may be affected by multiple genes and environmental factors.

VKORC1 and CYP2C9 Genotype Variations in Relation to Warfarin Dosing in Korean Stroke Patients / 대한뇌졸중학회지

BACKGROUND AND PURPOSE: Variant alleles of CYP2C9 and VKORC1 account for differences in anticoagulation response. We sought to establish a warfarin dosing formula for individualized target International Normalization Ratio of Prothrombin Times (INRs) using data from single nucleotide polymorphisms (SNPs) in VKORC1 and CYP2C9 in Korean patients. METHODS: Ischemic stroke patients displaying stable target INR for at least 3 months before enrollment were analyzed. Warfarin and vitamin K levels were measured to adjust for confounders. Phenotypes were defined using the 'warfarin response index' (WRI) defined as INR divided by the daily maintenance warfarin dose. We tested SNPs in CYP2C9 (3 sites: 430C>T (rs1799853), 1075A>C (rs1057910), 1076T>C) and VKORC1 (14 sites: 381C>T, 861C>A (rs17880887), 2653G>C, 3673A>G, 5496G>T, 5808T>G (r17882154), 6009C>T, 6484T>C (rs9934438), 6853C>T (rs17886369), 7566T>C, 8767G>C, 8814T>C, 9041G>A (rs17880624), and 9071G>T) using a standard sequencing method. Multivariate linear regression analysis was applied to establish the formula for warfarin dosage. RESULTS: All 204 patients had excellent drug compliance. The mean INR was 2.22 (+0.56) and mean daily maintenance dose of warfarin was 3.92 mg (+1.54). Patients with low WRI were younger (PC, PC, 3673G>A, 6484T>C, 6853C>G. 7566C>T, 9041G>A) (r2=1). Based on these data, a warfarin dosing formula was established. CONCLUSIONS: WRI is influenced by age, body mass index and SNPs in VKORC1 and CYP2C9 in Korean stroke patients. The obtained warfarin dosing formula may be clinically applicable.

Lack of Replication of Genetic Association with Body Mass Index Detected by Genome-wide Association Study

Obesity provokes many serious human diseases, including various cardiovascular diseases and diabetes. Body mass index (BMI) is a highly heritable trait that is broadly used to diagnose obesity. To identify genetic loci associated with obesity in Asians, we conducted a genome-wide association study (GWAS) of a population of Korean adults (n=6,742, age 40~60 years) and detected six BMI risk loci (TNR, FAM124B, RGS12, NFE2L3, MC4R and FTO) having p<1x10(-5). However, in the replication study, only melanocortin 4 receptor gene (MC4R) (rs9946888, p=4.58x10(-7)) was replicated with marginal significance (p<0.05) in the second cohort (n=5,102, age 40~60 years). This study indicates that each locus associated with BMI has very weak genetic effect.

Association of the X-linked Androgen Receptor Leu57Gln Polymorphism with Monomelic Amyotrophy

Monomelic amyotrophy (MA), also known as Hirayama disease, occurs mainly in young men and manifests as weakness and wasting of the muscles of the distal upper limbs. Here, we sought to identify a genetic basis for MA. Given the predominance of MA in males, we focused on candidate neurological disease genes located on the X chromosome, selecting two X-linked candidate genes, androgen receptor (AR) and ubiquitin-like modifier activating enzyme 1 (UBA1). Screening for genetic variants using patients' genomic DNA revealed three known genetic variants in the coding region of the AR gene: one nonsynonymous single-nucleotide polymorphism (SNP; rs78686797) encoding Leu57Gln, and two variants of polymorphic trinucleotide repeat segments that encode polyglutamine (CAG repeat; rs5902610) and polyglycine (GGC repeat; rs3138869) tracts. Notably, the Leu57Gln polymorphism was found in two patients with MA from 24 MA patients, whereas no variants were found in 142 healthy male controls. However, the numbers of CAG and GGC repeats in the AR gene were within the normal range. These data suggest that the Leu57Gln polymorphism encoded by the X-linked AR gene may contribute to the development of MA.

Association of a c.1084A>G (p.Thr362Ala)Variant in the DCTN4 Gene with Wilson Disease

PURPOSE: Wilson disease is an autosomal recessive disorder which causes excessive copper accumulation in the hepatic region. So far, ATP7B gene is the only disease-causing gene of Wilson disease known to date. However, ATP7B mutations have not been found in ~15% of the patients. This study was performed to identify any causative gene in Wilson disease patients without an ATP7B mutation in either allele. MATERIALS AND METHODS: The sequence of the coding regions and exon-intron boundaries of the five ATP7B-interacting genes, ATOX1, COMMD1, GLRX, DCTN4, and ZBTB16, were analyzed in the 12 patients with Wilson disease. RESULTS: Three nonsynonymous variants including c.1084A>G (p.Thr362Ala) in the exon 12 of the DCTN4 gene were identified in the patients examined. Among these, only p.Thr362Ala was predicted as possibly damaging protein function by in silico analysis. Examination of allele frequency of c.1084A>G (p.Thr362Ala) variant in the 176 patients with Wilson disease and in the 414 normal subjects revealed that the variant was more prevalent in the Wilson disease patients (odds ratio [OR]=3.14, 95% confidence interval=1.36-7.22, P=0.0094). CONCLUSION: Our result suggests that c.1084A>G (p.Thr362Ala) in the ATP7B-interacting DCTN4 gene may be associated with the pathogenesis of Wilson disease.

Exome Sequencing in Mendelian Disorders

More than 7,000 rare Mendelian diseases have been reported, but less than half of all rare monogenic disorders has been discovered. In addition, the majority of mutations that are known to cause Mendelian disorders are located in protein-coding regions. Therefore, exome sequencing is an efficient strategy to selectively sequence the coding regions of the human genome to identify novel genes associated with rare genetic disorders. The "exome" represents all of the exons in the human genome, constituting about 1.5% of the human genome. Exome sequencing is carried out by targeted capture and intense parallel sequencing. After the first report of successful exome sequencing for the identification of causal genes and mutations in Freeman Sheldon syndrome, exome sequencing has become a standard approach to identify genes in rare Mendelian disorders. Exome sequencing is also used to search the causal genes and variants in complex diseases. The successful use of exome sequencing in Mendelian disorders and complex diseases will facilitate the development of personalized genomic medicine.

Identification and Characterization of Human Genes Targeted by Natural Selection

The human genome has evolved as a consequence of evolutionary forces, such as natural selection. In this study, we investigated natural selection on the human genes by comparing the numbers of nonsynonymous(NS) and synonymous (S) mutations in individual genes. We initially collected all coding SNP data of all human genes from the public dbSNP. Among the human genes, we selected 3 different selection groups of genes: positively selected genes (NS/S > or = 3), negatively selected genes (NS/S < or = 1/3) and neutral selection genes (0.9 < NS/S < 1.1). We characterized human genes targeted by natural selection. Negatively selected human genes were markedly associated with disease occurrence, but not positively selected genes. Interestingly, positively selected genes displayed an increase in potentially deleterious nonsynonymous SNPs with an increased frequency of tryptophan and tyrosine residues, suggesting a correlation with protective effects against human disease. Furthermore, our nonsynonymous/synonymous ratio data imply that specific human genes, such as ALMS1 and SPTBN5 genes, are differentially selected among distinct populations. We confirmed that inferences of natural selection using the NS/S ratio can be used extensively to identify functional genes selected during the evolutionary adaptation process.