Analysis of East Asia genetic substructure using genome-wide SNP arrays

Forensic analysis and sequence variation of 133 STRs in the Hakka population

Introduction: Short Tandem Repeats (STRs) are highly valuable genetic markers in forensic science. However, the conventional PCR-CE technique has limitations, and the emergence of massively parallel sequencing (MPS) technology presents new opportunities for STR analysis. Yet, there is limited research on Chinese population diversity using MPS. Methods: In this study, we obtained genotype data for 52 A-STRs and 81 Y-STRs from the Hakka population in Meizhou, Guangdong, China, using the Forensic Analysis System Multiplecues SetB Kit on the MGISEQ-2000 platform. Results: Our findings demonstrate that these 133 STRs are highly efficient for forensic applications within the Meizhou Hakka population. Statistical analysis revealed Hobs values ranging from 0.61306 to 0.91083 and Hexp values ranging from 0.59156 to 0.91497 for A-STRs based on length polymorphism. For sequence polymorphism, Hobs values ranged from 0.61306 to 0.94586, and Hexp values fluctuated between 0.59156 and 0.94487. The CPE values were 1-5.0779620E-21 and 1-3.257436E-24 for length and sequence polymorphism, respectively, while the CPD values were 1-1.727007E-59 and 1-5.517015E-66, respectively. Among the 80 Y-STR loci, the HD values for length and sequence polymorphism were 0.99764282 and 0.99894195, respectively. The HMP values stood at 0.00418102 and 0.00288427, respectively, and the DC values were 0.75502742 and 0.83363803, respectively. For the 52 A-STR loci, we identified 554 and 989 distinct alleles based on length and sequence polymorphisms, respectively. For the 81 Y-STR loci, 464 and 652 unique alleles were detected at the length and sequence level, respectively. Population genetic analysis revealed that the Meizhou Hakka population has a close kinship relationship with the Asian populations THI and KOR based on length polymorphism data of A-STRs. Conversely, based on length polymorphism data of Y-STRs, the Meizhou Hakka population has the closest kinship relationship with the Henan Han population. Discussion: Overall, the variation information of repeat region sequences significantly enhances the forensic identification efficacy of STR genetic markers, providing an essential database for forensic individual and paternity testing in this region. Additionally, the data generated by our study will serve as a vital resource for research into the genetic structure and historical origins of the Meizhou Hakka population.

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

Genomic analysis reveals geography rather than culture as the predominant factor shaping genetic variation in northern Kenyan human populations

OBJECTIVES

The aim of this study was to characterize the genetic relationships within and among four neighboring ethnolinguistic groups in northern Kenya in light of cultural relationships to understand the extent to which geography and culture shape patterns of genetic variation.

MATERIALS AND METHODS

We collected DNA and demographic information pertaining to aspects of social identity and heritage from 572 individuals across the Turkana, Samburu, Waso Borana, and Rendille of northern Kenya. We sampled individuals across a total of nine clans from these four groups and, additionally, three territorial sections within the Turkana and successfully genotyped 376 individuals.

RESULTS

Here we report that geography predominately shapes genetic variation within and among human groups in northern Kenya. We observed a clinal pattern of genetic variation that mirrors the overall geographic distribution of the individuals we sampled. We also found relatively higher rates of intermarriage between the Rendille and Samburu and evidence of gene flow between them that reflect these higher rates of intermarriage. Among the Turkana, we observed strong recent genetic substructuring based on territorial section affiliation. Within ethnolinguistic groups, we found that Y chromosome haplotypes do not consistently cluster by natal clan affiliation. Finally, we found that sampled populations that are geographically closer have lower genetic differentiation, and that cultural similarity does not predict genetic similarity as a whole across these northern Kenyan populations.

DISCUSSION

Overall, the results from this study highlight the importance of geography, even on a local geographic scale, in shaping observed patterns of genetic variation in human populations.

An efficient ancestry informative SNPs panel for further discriminating East Asian populations

In forensic genetics, the use of ancestry informative single-nucleotide polymorphisms (AISNPs) panels can narrow the direction of the investigation by estimating an individual's biogeographic ancestry. However, distinguishing subgroups within continental regions requires more specific panels. In this study, we screened 19 AISNPs from the 1000 Genomes Project (1KG) based on their F_ST values to distinguish target populations in East Asia and obtained genotypes through SNaPshot. The 19 AISNPs could divide the global population of the 1KG into five clusters and could further divide the East Asian population into four clusters: Japanese, Han Chinese, Dai Chinese, and Kinh in Ho Chi Minh City of Vietnam. In summary, the 19-AISNP panel may serve as a useful and cost-effective tool for forensic ancestry inference in East Asian populations at a finer scale.

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.

Effects of interracial experience on the race preferences of infants

Over their first year of life, infants express visual preferences for own- versus other-race faces. This developmental transition has primarily been investigated in monoracial societies where infants have limited personal or societal contact with other races. We investigated whether previously reported visual preferences for race generalize to a multiracial society (i.e., Singapore). In addition, we investigated effects of caregiver race on visual preferences for race. In Experiment 1, race preferences were measured at 3, 6, and 9 months of age for own-race (Chinese) versus other-race (Indian) faces in infants with no regular interaction with Indian-race individuals. Singaporean infants displayed a significant visual preference for Indian-race faces at each age group. Furthermore, infants raised with other-race caregivers demonstrated an age-related increase in other-race visual preferences. The visual preferences of infants for other-race faces were predicted by the extent of other-race contact. In Experiment 2, we confirmed that an other-race visual preference was not exclusive to Indian faces in a sample of 6-month-old Singaporean Chinese infants who demonstrated a similar other-race visual preference for Caucasian faces over Chinese faces. Findings are discussed in terms of the influence of other-race contact on visual preferences for race in infants.

Population Genetic Polymorphism of Skeletal Muscle Strength Related Genes in Five Ethnic Minorities in North China

Musculoskeletal performance is a complex trait influenced by environmental and genetic factors, and it has different manifestations in different populations. Heilongjiang province, located in northern China, is a multi-ethnic region with human cultures dating back to the Paleolithic Age. The Daur, Hezhen, Ewenki, Mongolian and Manchu ethnic groups in Heilongjiang province may have strong physical fitness to a certain extent. Based on the genetic characteristics of significant correlation between some important genes and skeletal muscle function, this study selected 23 SNPs of skeletal muscle strength-related genes and analyzed the distribution of these loci and genetic diversity in the five ethnic groups. Use Haploview (version 4.1) software to calculate the chi-square and the Hardy-Weinberg equilibrium to assess the difference between the two ethnic groups. Use R (version 4.0.2) software to perform principal component analysis of different ethnic groups. Use MEGA (version 7.0) software to construct the phylogenetic tree of different ethnic groups. Use POPGENE (version 1.32) software to calculate the heterozygosity and the F_ST values of 23 SNPs. Use Arlequin (version 3.5.2.2) software to analyze molecular variance (AMOVA) among 31 populations. The results showed that there was haplotype diversity of VDR, angiotensin-converting enzyme, ACTN3, EPO and IGF1 genes in the five ethnic groups, and there were genetic differences in the distribution of these genes in the five ethnic groups. Among them, the average gene heterozygosity (AVE_HET) of the 23 SNPs in the five populations was 0.398. The F_ST values of the 23 SNPs among the five ethnic groups varied from 0.0011 to 0.0137. According to the principal component analysis, the genetic distance of Daur, Mongolian and Ewenki is relatively close. According to the phylogenetic tree, the five ethnic groups are clustered together with the Asian population. These data will enrich existing genetic information of ethnic minorities.

Genetic ancestry plays a central role in population pharmacogenomics

Recent studies have pointed out the essential role of genetic ancestry in population pharmacogenetics. In this study, we analyzed the whole-genome sequencing data from The 1000 Genomes Project (Phase 3) and the pharmacogenetic information from Drug Bank, PharmGKB, PharmaADME, and Biotransformation. Here we show that ancestry-informative markers are enriched in pharmacogenetic loci, suggesting that trans-ancestry differentiation must be carefully considered in population pharmacogenetics studies. Ancestry-informative pharmacogenetic loci are located in both protein-coding and non-protein-coding regions, illustrating that a whole-genome analysis is necessary for an unbiased examination over pharmacogenetic loci. Finally, those ancestry-informative pharmacogenetic loci that target multiple drugs are often a functional variant, which reflects their importance in biological functions and pathways. In summary, we develop an efficient algorithm for an ultrahigh-dimensional principal component analysis. We create genetic catalogs of ancestry-informative markers and genes. We explore pharmacogenetic patterns and establish a high-accuracy prediction panel of genetic ancestry. Moreover, we construct a genetic ancestry pharmacogenomic database Genetic Ancestry PhD (http://hcyang.stat.sinica.edu.tw/databases/genetic_ancestry_phd/).

Hsin-Chou Yang et al. examine population structure in several genomic databases and identify that pharmacogenetic loci are enriched for markers of genetic ancestry. Their results suggest that genetic ancestry must be carefully considered in population pharmacogenetics studies.

The number of SNPs required for distinguishing Japanese from other East Asians

In some cases, it is necessary to estimate the national origin of an unknown subject in forensic medicine. The use of single nucleotide polymorphism (SNP) markers appears to be very effective for this purpose, since genome-wide SNP genotype data of many human populations are publicly available. In this study, we examined the number of SNPs that could objectively and accurately distinguish Japanese subjects (1KG-JPT) from the other East Asians (1KG-CDX, -CHB, -CHS, and -KHV) using the combination of principal component analysis and hierarchical cluster analysis. A computer simulation showed that approximately 3000 randomly selected SNPs were enough for the discrimination. Our results suggest that at least a 0.024% coverage is needed in the next generation sequencing experiment to objectively determine whether an unknown person is Japanese or not if the amount of DNA sample from him/her is insufficient or the quality is low.

Influence of Paleolithic range contraction, admixture and long-distance dispersal on genetic gradients of modern humans in Asia

Cavalli-Sforza and coauthors originally explored the genetic variation of modern humans throughout the world and observed an overall east-west genetic gradient in Asia. However, the specific environmental and population genetics processes causing this gradient were not formally investigated and promoted discussion in recent studies. Here we studied the influence of diverse environmental and population genetics processes on Asian genetic gradients and identified which could have produced the observed gradient. To do so, we performed extensive spatially-explicit computer simulations of genetic data under the following scenarios: (a) variable levels of admixture between Paleolithic and Neolithic populations, (b) migration through long-distance dispersal (LDD), (c) Paleolithic range contraction induced by the last glacial maximum (LGM), and (d) Neolithic range expansions from one or two geographic origins (the Fertile Crescent and the Yangzi and Yellow River Basins). Next, we estimated genetic gradients from the simulated data and we found that they were sensible to the analysed processes, especially to the range contraction induced by LGM and to the number of Neolithic expansions. Some scenarios were compatible with the observed east-west genetic gradient, such as the Paleolithic expansion with a range contraction induced by the LGM or two Neolithic range expansions from both the east and the west. In general, LDD increased the variance of genetic gradients among simulations. We interpreted the obtained gradients as a consequence of both allele surfing caused by range expansions and isolation by distance along the vast east-west geographic axis of this continent.

Differentiation of Hispanic biogeographic ancestry with 80 ancestry informative markers

Ancestry informative single nucleotide polymorphisms (SNPs) can identify biogeographic ancestry (BGA); however, population substructure and relatively recent admixture can make differentiation difficult in heterogeneous Hispanic populations. Utilizing unrelated individuals from the Genomic Origins and Admixture in Latinos dataset (GOAL, n = 160), we designed an 80 SNP panel (Setser80) that accurately depicts BGA through STRUCTURE and PCA. We compared our Setser80 to the Seldin and Kidd panels via resampling simulations, which models data based on allele frequencies. We incorporated Admixed American 1000 Genomes populations (1000 G, n = 347), into a combined populations dataset to determine robustness. Using multinomial logistic regression (MLR), we compared the 3 panels on the combined dataset and found overall MLR classification accuracies: 93.2% Setser80, 87.9% Seldin panel, 71.4% Kidd panel. Naïve Bayesian classification had similar results on the combined dataset: 91.5% Setser80, 84.7% Seldin panel, 71.1% Kidd panel. Although Peru and Mexico were absent from panel design, we achieved high classification accuracy on the combined populations for Peru (MLR = 100%, naïve Bayes = 98%), and Mexico (MLR = 90%, naïve Bayes = 83.4%) as evidence of the portability of the Setser80. Our results indicate the Setser80 SNP panel can reliably classify BGA for individuals of presumed Hispanic origin.

Insights Into Breast Cancer in the East vs the West: A Review

Importance

During the past few decades, the incidence of breast cancer (BC) has been increasing rapidly in East Asia, and BC is currently the most common cancer in several countries. The rising incidence is likely related to changing lifestyle and environmental factors in addition to the increase in early diagnosis with BC awareness and screening. The understanding and management of BC are generally based on research and data from the West. However, emerging differences in BC epidemiology and tumor and host biology in Asian populations may be clinically relevant.

Observations

A higher proportion of premenopausal BCs occur in Asia, although this factor is possibly an age-cohort effect. Although the relative frequencies of different immunohistochemical subtypes of BC may be similar between the East and West, the higher prevalence of luminal B subtypes with more frequent mutations in TP53 may be confounded by disparities in early detection. In addition, Asian BCs appear to harbor a more immune-active microenvironment than BCs in the West. The spectra of germline mutations in BC predisposition genes and single-nucleotide polymorphisms contributing to BC risk vary with ethnicity as well. Differences in tolerability of certain cytotoxic and targeted agents used in BC treatment may be associated with pharmacogenomic factors, whereas the lower body mass of the average woman in East Asia may contribute to higher toxicities from drugs administered at fixed doses. Phenotypic characteristics, such as lower breast volume, may influence the type of surgery performed in East Asian women. On the other hand, increased breast density may affect the sensitivity of mammography in detecting BCs, limiting the benefits of screening mammography.

Conclusions and Relevance

Breast cancer has become a major health problem in Asia. The inclusion of more women from Asia in clinical trials and epidemiologic and translational studies may help unravel the interethnic heterogeneity of BCs and elucidate the complex interplay between environmental and intrinsic factors in its pathogenesis. These insights may help to refine prevention, diagnosis, and management strategies for BC in the setting of ethnic diversity.

Inconsistency in race and ethnic classification in pharmacogenetics studies and its potential clinical implications

Introduction

Racial and ethnic categories are frequently used in pharmacogenetics literature to stratify patients; however, these categories can be inconsistent across different studies. To address the ongoing debate on the applicability of traditional concepts of race and ethnicity in the context of precision medicine, we aimed to review the application of current racial and ethnic categories in pharmacogenetics and its potential impact on clinical care.

Methods

One hundred and three total pharmacogenetics papers involving the CYP2C9, CYP2C19, and CYP2D6 genes were analyzed for their country of origin, racial, and ethnic categories used, and allele frequency data. Correspondence between the major continental racial categories promulgated by National Institutes of Health (NIH) and those reported by the pharmacogenetics papers was evaluated.

Results

The racial and ethnic categories used in the papers we analyzed were highly heterogeneous. In total, we found 66 different racial and ethnic categories used which fall under the NIH race category “White”, 47 different racial and ethnic categories for “Asian”, and 62 different categories for “Black”. The number of categories used varied widely based on country of origin: Japan used the highest number of different categories for “White” with 17, Malaysia used the highest number for “Asian” with 24, and the US used the highest number for “Black” with 28. Significant variation in allele frequency between different ethnic subgroups was identified within 3 major continental racial categories.

Conclusion

Our analysis showed that racial and ethnic classification is highly inconsistent across different papers as well as between different countries. Evidence-based consensus is necessary for optimal use of self-identified race as well as geographical ancestry in pharmacogenetics. Common taxonomy of geographical ancestry which reflects specifics of particular countries and is accepted by the entire scientific community can facilitate reproducible pharmacogenetic research and clinical implementation of its results.

A single nucleotide polymorphism panel for individual identification and ancestry assignment in Caucasians and four East and Southeast Asian populations using a machine learning classifier

Single nucleotide polymorphism (SNP) profiling is an effective means of individual identification and ancestry inferences in forensic genetics. This study established a SNP panel for the simultaneous individual identification and ancestry assignment of Caucasian and four East and Southeast Asian populations. We analyzed 220 SNPs (125 autosomal, 17 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs) of the DNA samples from 563 unrelated individuals of five populations (89 Caucasian, 234 Taiwanese Han, 90 Filipino, 79 Indonesian and 71 Vietnamese) and 18 degraded DNA samples. Informativeness for assignment (In) was used to select ancestry informative SNPs (AISNPs). A machine learning classifier, support vector machine (SVM), was used for ancestry assignment. Of the 220 SNPs, 62 were individual identification SNPs (IISNPs) (51 autosomal and 11 X-chromosomal SNPs) and 191 were AISNPs (100 autosomal, 13 X-chromosomal, 30 Y-chromosomal, and 48 mitochondrial SNPs). The 51 autosomal IISNPs offered cumulative random match probabilities (cRMPs) ranging from 1.56 × 10^-21 to 3.16 × 10^-22 among these five populations. Using AISNPs with the SVM, the overall accuracy rate of ancestry inference achieved in the testing dataset between Caucasian, Taiwanese Han, and Filipino populations was 88.9%, whereas it was 70.0% between Caucasians and each of the four East and Southeast Asian populations. For the 18 degraded DNA samples with incomplete profiling, the accuracy rate of ancestry assignment was 94.4%. We have developed a 220-SNP panel for simultaneous individual identification and ethnic origin differentiation between Caucasian and the four East and Southeast Asian populations. This SNP panel may assist with DNA analysis of forensic casework.

Human ancestry indentification under resource constraints -- what can one chromosome tell us about human biogeographical ancestry?

BACKGROUND

While continental level ancestry is relatively simple using genomic information, distinguishing between individuals from closely associated sub-populations (e.g., from the same continent) is still a difficult challenge.

METHODS

We study the problem of predicting human biogeographical ancestry from genomic data under resource constraints. In particular, we focus on the case where the analysis is constrained to using single nucleotide polymorphisms (SNPs) from just one chromosome. We propose methods to construct such ancestry informative SNP panels using correlation-based and outlier-based methods.

RESULTS

We accessed the performance of the proposed SNP panels derived from just one chromosome, using data from the 1000 Genome Project, Phase 3. For continental-level ancestry classification, we achieved an overall classification rate of 96.75% using 206 single nucleotide polymorphisms (SNPs). For sub-population level ancestry prediction, we achieved an average pairwise binary classification rates as follows: subpopulations in Europe: 76.6% (58 SNPs); Africa: 87.02% (87 SNPs); East Asia: 73.30% (68 SNPs); South Asia: 81.14% (75 SNPs); America: 85.85% (68 SNPs).

CONCLUSION

Our results demonstrate that one single chromosome (in particular, Chromosome 1), if carefully analyzed, could hold enough information for accurate prediction of human biogeographical ancestry. This has significant implications in terms of the computational resources required for analysis of ancestry, and in the applications of such analyses, such as in studies of genetic diseases, forensics, and soft biometrics.

Genome-wide mapping of quantitative trait loci in admixed populations using mixed linear model and Bayesian multiple regression analysis

Background

Population stratification and cryptic relationships have been the main sources of excessive false-positives and false-negatives in population-based association studies. Many methods have been developed to model these confounding factors and minimize their impact on the results of genome-wide association studies. In most of these methods, a two-stage approach is applied where: (1) methods are used to determine if there is a population structure in the sample dataset and (2) the effects of population structure are corrected either by modeling it or by running a separate analysis within each sub-population. The objective of this study was to evaluate the impact of population structure on the accuracy and power of genome-wide association studies using a Bayesian multiple regression method.

Methods

We conducted a genome-wide association study in a stochastically simulated admixed population. The genome was composed of six chromosomes, each with 1000 markers. Fifteen segregating quantitative trait loci contributed to the genetic variation of a quantitative trait with heritability of 0.30. The impact of genetic relationships and breed composition (BC) on three analysis methods were evaluated: single marker simple regression (SMR), single marker mixed linear model (MLM) and Bayesian multiple-regression analysis (BMR). Each method was fitted with and without BC. Accuracy, power, false-positive rate and the positive predictive value of each method were calculated and used for comparison.

Results

SMR and BMR, both without BC, were ranked as the worst and the best performing approaches, respectively. Our results showed that, while explicit modeling of genetic relationships and BC is essential for models SMR and MLM, BMR can disregard them and yet result in a higher power without compromising its false-positive rate.

Conclusions

This study showed that the Bayesian multiple-regression analysis is robust to population structure and to relationships among study subjects and performs better than a single marker mixed linear model approach.

Japanese GWAS identifies variants for bust-size, dysmenorrhea, and menstrual fever that are eQTLs for relevant protein-coding or long non-coding RNAs

Traits related to primary and secondary sexual characteristics greatly impact females during puberty and day-to-day adult life. Therefore, we performed a GWAS analysis of 11,348 Japanese female volunteers and 22 gynecology-related phenotypic variables, and identified significant associations for bust-size, menstrual pain (dysmenorrhea) severity, and menstrual fever. Bust-size analysis identified significant association signals in CCDC170-ESR1 (rs6557160; P = 1.7 × 10-16) and KCNU1-ZNF703 (rs146992477; P = 6.2 × 10-9) and found that one-third of known European-ancestry associations were also present in Japanese. eQTL data points to CCDC170 and ZNF703 as those signals' functional targets. For menstrual fever, we identified a novel association in OPRM1 (rs17181171; P = 2.0 × 10-8), for which top variants were eQTLs in multiple tissues. A known dysmenorrhea signal near NGF replicated in our data (rs12030576; P = 1.1 × 10-19) and was associated with RP4-663N10.1 expression, a putative lncRNA enhancer of NGF, while a novel dysmenorrhea signal in the IL1 locus (rs80111889; P = 1.9 × 10-16) contained SNPs previously associated with endometriosis, and GWAS SNPs were most significantly associated with IL1A expression. By combining regional imputation with colocalization analysis of GWAS/eQTL signals along with integrated annotation with epigenomic data, this study further refines the sets of candidate causal variants and target genes for these known and novel gynecology-related trait loci.

Japaneseplex: A forensic SNP assay for identification of Japanese people using Japanese-specific alleles

It is sometimes necessary to determine whether a forensic biological sample came from a Japanese person. In this study, we developed a 60-locus SNP assay designed for the differentiation of Japanese people from other East Asians using entirely and nearly Japanese-specific alleles. This multiplex assay consisted of 6 independent PCR reactions followed by single nucleotide extension. The average number and standard deviation of Japanese-specific alleles possessed by an individual were 0.81 ± 0.93 in 108 Koreans from Seoul, 8.87 ± 2.89 in 103 Japanese from Tottori, 17.20 ± 3.80 in 88 Japanese from Okinawa, and 0 in 220 Han Chinese from Wuxi and Changsha. The Koreans had 0-4 Japanese-specific alleles per individual, whereas the Japanese had 4-26 Japanese-specific alleles. Almost all Japanese were distinguished from the Koreans and other people by the factorial correspondence and principal component analyses. The Snipper program was also useful to estimate the degree of Japaneseness. The method described here was successfully applied to the differentiation of Japanese from non-Japanese people in forensic cases. This Japanese-specific SNP assay was named Japaneseplex.

Robust inference of population structure from next-generation sequencing data with systematic differences in sequencing

Motivation

Inferring population structure is important for both population genetics and genetic epidemiology. Principal components analysis (PCA) has been effective in ascertaining population structure with array genotype data but can be difficult to use with sequencing data, especially when low depth leads to uncertainty in called genotypes. Because PCA is sensitive to differences in variability, PCA using sequencing data can result in components that correspond to differences in sequencing quality (read depth and error rate), rather than differences in population structure. We demonstrate that even existing methods for PCA specifically designed for sequencing data can still yield biased conclusions when used with data having sequencing properties that are systematically different across different groups of samples (i.e. sequencing groups). This situation can arise in population genetics when combining sequencing data from different studies, or in genetic epidemiology when using historical controls such as samples from the 1000 Genomes Project.

Results

To allow inference on population structure using PCA in these situations, we provide an approach that is based on using sequencing reads directly without calling genotypes. Our approach is to adjust the data from different sequencing groups to have the same read depth and error rate so that PCA does not generate spurious components representing sequencing quality. To accomplish this, we have developed a subsampling procedure to match the depth distributions in different sequencing groups, and a read-flipping procedure to match the error rates. We average over subsamples and read flips to minimize loss of information. We demonstrate the utility of our approach using two datasets from 1000 Genomes, and further evaluate it using simulation studies.

Availability and implementation

TASER-PC software is publicly available at http://web1.sph.emory.edu/users/yhu30/software.html.

Contact

yijuan.hu@emory.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome-wide association study of self-reported food reactions in Japanese identifies shrimp and peach specific loci in the HLA-DR/DQ gene region

Food allergy is an increasingly important health problem in the world. Several genome-wide association studies (GWAS) focused on European ancestry samples have identified food allergy-specific loci in the HLA class II region. We conducted GWAS of self-reported reactivity with common foods using the data from 11011 Japanese women and identified shrimp and peach allergy-specific loci in the HLA-DR/DQ gene region tagged by rs74995702 (P = 6.30 × 10⁻¹⁷, OR = 1.91) and rs28359884 (P = 2.3 × 10⁻¹², OR = 1.80), respectively. After HLA imputation using a Japanese population-specific reference, the most strongly associated haplotype was HLA-DRB1*04:05-HLA-DQB1*04:01 for shrimp allergy (P = 3.92 × 10⁻¹⁹, OR = 1.99) and HLA-DRB1*09:01-HLA-DQB1*03:03 for peach allergy (P = 1.15 × 10⁻⁷, OR = 1.68). Additionally, both allergies’ associated variants were eQTLs for several HLA genes, with HLA-DQA2 the single eQTL gene shared between the two traits. Our study suggests that allergy to certain foods may be related to genetic differences that tag both HLA alleles having particular epitope binding specificities as well as variants modulating expression of particular HLA genes. Investigating this further could increase our understanding of food allergy aetiology and potentially lead to better therapeutic strategies for allergen immunotherapies.

Korean Variant Archive (KOVA): a reference database of genetic variations in the Korean population

Despite efforts to interrogate human genome variation through large-scale databases, systematic preference toward populations of Caucasian descendants has resulted in unintended reduction of power in studying non-Caucasians. Here we report a compilation of coding variants from 1,055 healthy Korean individuals (KOVA; Korean Variant Archive). The samples were sequenced to a mean depth of 75x, yielding 101 singleton variants per individual. Population genetics analysis demonstrates that the Korean population is a distinct ethnic group comparable to other discrete ethnic groups in Africa and Europe, providing a rationale for such independent genomic datasets. Indeed, KOVA conferred 22.8% increased variant filtering power in addition to Exome Aggregation Consortium (ExAC) when used on Korean exomes. Functional assessment of nonsynonymous variant supported the presence of purifying selection in Koreans. Analysis of copy number variants detected 5.2 deletions and 10.3 amplifications per individual with an increased fraction of novel variants among smaller and rarer copy number variable segments. We also report a list of germline variants that are associated with increased tumor susceptibility. This catalog can function as a critical addition to the pre-existing variant databases in pursuing genetic studies of Korean individuals.

Population diversity of the genetically determined TTR expression in human tissues and its implications in TTR amyloidosis

BACKGROUND

Transthyretin (TTR) amyloidosis is a hereditary disease with a complex genotype-phenotype correlation. We conducted a literature survey to define the clinical landscape of TTR amyloidosis across populations worldwide. Then, we investigated whether the genetically determined TTR expression differs among human populations, contributing to the differences observed in patients. Polygenic scores for genetically determined TTR expression in 14 clinically relevant tissues were constructed using data from the GTEx (Genotype-Tissue Expression) project and tested in the samples from the 1,000 Genomes Project.

RESULTS

We observed differences among the ancestral groups and, to a lesser extent, among the investigated populations within the ancestry groups. Scandinavian populations differed in their genetically determined TTR expression of skeletal muscle tissue with respect to Southern Europeans (p = 6.79*10^-6). This is in line with epidemiological data related to Swedish and Portuguese TTR Val30Met endemic areas. Familial amyloidotic cardiomyopathy (TTR deposits occur primarily in heart tissues) presents clinical variability among human populations, a finding that agrees with the among-ancestry diversity of genetically determined TTR expression in heart tissues (i.e., Atrial Appendage p = 4.55*10^-28; Left Ventricle p = 6.54*10^-35).

CONCLUSIONS

Genetically determined TTR expression varied across human populations. This might contribute to the genotype-phenotype correlation of TTR amyloidosis.

Multiple genetic variants associated with primary biliary cirrhosis in a Han Chinese population

Multiple genome-wide association studies of primary biliary cirrhosis (PBC) in both European and Japanese ancestries have shown significant associations of many genetic loci contributing to the susceptibility to PBC. Major differences in susceptibility loci between these two population groups were observed. In this study, we examined whether the most significant loci observed in either European and/or Japanese cohorts are associated with PBC in a Han Chinese population. In 1070 PBC patients and 1198 controls, we observed highly significant associations at CD80 (rs2293370, P = 2.67 × 10(-8)) and TNFSF15 (rs4979462, P = 3.86 × 10(-8)) and significant associations at 17q12-21 (rs9303277), PDGFB (rs715505), NF-κB1 (rs7665090), IL12RB2 (rs11209050), and STAT4 (rs7574865; all corrected P values <0.01). However, no association was observed for POU2AF1 (rs4938534), IL12A (rs485499 and rs2366408), IL7R (rs6897932), CXCR5 (rs715412), SOCS1 (rs725613), and TNFRSF1A (rs1800693). STAT4 (rs7574865) was strongly associated after additional control samples were analyzed. Our study is the first large-scale genetic analysis in a Han Chinese PBC cohort. These results do not only reflect that Han Chinese PBC patients share common genetic susceptibility genes with both their Japanese and European counterparts but also suggest a distinctly different genetic susceptibility profile.

Characterizing Race/Ethnicity and Genetic Ancestry for 100,000 Subjects in the Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort

Using genome-wide genotypes, we characterized the genetic structure of 103,006 participants in the Kaiser Permanente Northern California multi-ethnic Genetic Epidemiology Research on Adult Health and Aging Cohort and analyzed the relationship to self-reported race/ethnicity. Participants endorsed any of 23 race/ethnicity/nationality categories, which were collapsed into seven major race/ethnicity groups. By self-report the cohort is 80.8% white and 19.2% minority; 93.8% endorsed a single race/ethnicity group, while 6.2% endorsed two or more. Principal component (PC) and admixture analyses were generally consistent with prior studies. Approximately 17% of subjects had genetic ancestry from more than one continent, and 12% were genetically admixed, considering only nonadjacent geographical origins. Self-reported whites were spread on a continuum along the first two PCs, indicating extensive mixing among European nationalities. Self-identified East Asian nationalities correlated with genetic clustering, consistent with extensive endogamy. Individuals of mixed East Asian-European genetic ancestry were easily identified; we also observed a modest amount of European genetic ancestry in individuals self-identified as Filipinos. Self-reported African Americans and Latinos showed extensive European and African genetic ancestry, and Native American genetic ancestry for the latter. Among 3741 genetically identified parent-child pairs, 93% were concordant for self-reported race/ethnicity; among 2018 genetically identified full-sib pairs, 96% were concordant; the lower rate for parent-child pairs was largely due to intermarriage. The parent-child pairs revealed a trend toward increasing exogamy over time; the presence in the cohort of individuals endorsing multiple race/ethnicity categories creates interesting challenges and future opportunities for genetic epidemiologic studies.

Ancestry informative markers for distinguishing between Thai populations based on genome-wide association datasets

The main purpose of this work was to identify a set of AIMs that stratify the genetic structure and diversity of the Thai population from a high-throughput autosomal genome-wide association study. In this study, more than one million SNPs from the international HapMap database and the Thai depression genome-wide association study have been examined to identify ancestry informative markers (AIMs) that distinguish between Thai populations. An efficient strategy is proposed to identify and characterize such SNPs and to test high-resolution SNP data from international HapMap populations. The best AIMs are identified to stratify the population and to infer genetic ancestry structure. A total of 124 AIMs were clearly clustered geographically across the continent, whereas only 89 AIMs stratified the Thai population from East Asian populations. Finally, a set of 273 AIMs was able to distinguish northern from southern Thai subpopulations. These markers will be of particular value in identifying the ethnic origins in regions where matching by self-reports is unavailable or unreliable, which usually occurs in real forensic cases.

Cytokine gene variations associated with trait and state anxiety in oncology patients and their family caregivers

PURPOSE

Anxiety is common among cancer patients and their family caregivers (FCs) and is associated with poorer outcomes. Recently, associations between inflammation and anxiety were identified. However, the relationship between variations in cytokine genes and anxiety warrants investigation. Therefore, phenotypic and genotypic characteristics associated with trait and state anxiety were evaluated in a sample of 167 oncology patients with breast, prostate, lung, or brain cancer and 85 of their FCs.

METHODS

Using multiple regression analyses, the associations between participants' demographic and clinical characteristics as well as variations in cytokine genes and trait and state anxiety were evaluated.

RESULTS

In the bivariate analyses, a number of phenotypic characteristics were associated with both trait and state anxiety (e.g., age, functional status). However, some associations were specific only to trait anxiety (e.g., number of comorbid conditions) or state anxiety (e.g., participation with a FC). Variations in three cytokine genes (i.e., interleukin (IL) 1 beta, IL1 receptor 2 (IL1R2), nuclear factor kappa beta 2 (NFKB2)) were associated with trait anxiety, and variations in two genes (i.e., IL1R2, tumor necrosis factor alpha (TNFA)) were associated with state anxiety.

CONCLUSIONS

These findings suggest that both trait and state anxiety need to be assessed in oncology patients and their FCs. Furthermore, variations in cytokine genes may contribute to higher levels of anxiety in oncology patients and their FCs.

Association between the NOTCH4 gene rs3131296 polymorphism with schizophrenia risk in the Chinese Zhuang population and Chinese Han population

BACKGROUND

Schizophrenia (SZ) is a common severe psychiatric disorder and a complex polygenic inherited disease that has not yet been fully interpreted. Heredity was proven to play an important role in the development of SZ. The association between the NOTCH4 gene rs3131296 polymorphism and SZ was reported to reach significance at the genome-wide level; therefore, it is necessary to replicate this association in other different populations.

METHODS

To evaluate the association of the NOTCH4 gene rs3131296 polymorphism with the risk for SZ, and to explore whether a significant association could be replicated in different ethnic groups of China, we conducted this case-control study on 282 SZ cases (188 Han and 94 Zhuang) and 282 controls (188 Han and 94 Zhuang) among the Chinese Zhuang and Han populations.

RESULTS

The results showed no statistically significant difference in the genotype or allele frequencies of the NOTCH4 gene variant rs3131296 between SZ patients and healthy controls in either the Zhuang or Han samples (p > 0.05). In addition, no significant difference was found in genotype or allele frequencies of the NOTCH4 gene variant rs3131296 between cases and controls in the combined samples including Zhuang and Han samples.

CONCLUSIONS

Our study failed to replicate the significant association between the NOTCH4 gene rs3131296 polymorphism and the risk for SZ.

Genetic, psychosocial, and demographic factors associated with social disinhibition in Mexican-origin youth

INTRODUCTION

The genetic heritability for sensation-seeking tendencies ranges from 40 to 60%. Sensation-seeking behaviors typically manifest during adolescence and are associated with alcohol and cigarette experimentation in adolescents. Social disinhibition is an aspect of sensation-seeking that is closely tied to cigarette and alcohol experimentation.

METHODS

We examined the contribution of candidate genes to social disinhibition among 1132 Mexican origin youth in Houston, Texas, adjusting for established demographic and psychosocial risk factors. Saliva samples were obtained at baseline in 2005-06, and social disinhibition and other psychosocial data were obtained in 2008-09. Participants were genotyped for 672 functional and tagging SNPs potentially related to sensation-seeking, risk-taking, smoking, and alcohol use.

RESULTS

Six SNPs were significantly associated with social disinhibition scores, after controlling for false discovery and adjusting for population stratification and relevant demographic/psychosocial characteristics. Minor alleles for three of the SNPs (rs1998220 on OPRM1; rs9534511 on HTR2A; and rs4938056 on HTR3B) were associated with increased risk of social disinhibition, while minor alleles for the other three SNPs (rs1003921 on KCNC1; rs16116 downstream of NPY; and rs16870286 on LINC00518) exhibited a protective effect. Age, linguistic acculturation, thrill and adventure-seeking, and drug and alcohol-seeking were all significantly positively associated with increased risk of social disinhibition in a multivariable model (P < 0.001).

CONCLUSIONS

These results add to our knowledge of genetic risk factors for social disinhibition. Additional research is needed to verify whether these SNPs are associated with social disinhibition among youth of different ethnicities and nationalities, and to elucidate whether and how these SNPs functionally contribute to social disinhibition.

Applying genome-wide gene-based expression quantitative trait locus mapping to study population ancestry and pharmacogenetics

BACKGROUND

Gene-based analysis has become popular in genomic research because of its appealing biological and statistical properties compared with those of a single-locus analysis. However, only a few, if any, studies have discussed a mapping of expression quantitative trait loci (eQTL) in a gene-based framework. Neither study has discussed ancestry-informative eQTL nor investigated their roles in pharmacogenetics by integrating single nucleotide polymorphism (SNP)-based eQTL (s-eQTL) and gene-based eQTL (g-eQTL).

RESULTS

In this g-eQTL mapping study, the transcript expression levels of genes (transcript-level genes; T-genes) were correlated with the SNPs of genes (sequence-level genes; S-genes) by using a method of gene-based partial least squares (PLS). Ancestry-informative transcripts were identified using a rank-score-based multivariate association test, and ancestry-informative eQTL were identified using Fisher's exact test. Furthermore, key ancestry-predictive eQTL were selected in a flexible discriminant analysis. We analyzed SNPs and gene expression of 210 independent people of African-, Asian- and European-descent. We identified numerous cis- and trans-acting g-eQTL and s-eQTL for each population by using PLS. We observed ancestry information enriched in eQTL. Furthermore, we identified 2 ancestry-informative eQTL associated with adverse drug reactions and/or drug response. Rs1045642, located on MDR1, is an ancestry-informative eQTL (P = 2.13E-13, using Fisher's exact test) associated with adverse drug reactions to amitriptyline and nortriptyline and drug responses to morphine. Rs20455, located in KIF6, is an ancestry-informative eQTL (P = 2.76E-23, using Fisher's exact test) associated with the response to statin drugs (e.g., pravastatin and atorvastatin). The ancestry-informative eQTL of drug biotransformation genes were also observed; cross-population cis-acting expression regulators included SPG7, TAP2, SLC7A7, and CYP4F2. Finally, we also identified key ancestry-predictive eQTL and established classification models with promising training and testing accuracies in separating samples from close populations.

CONCLUSIONS

In summary, we developed a gene-based PLS procedure and a SAS macro for identifying g-eQTL and s-eQTL. We established data archives of eQTL for global populations. The program and data archives are accessible at http://www.stat.sinica.edu.tw/hsinchou/genetics/eQTL/HapMapII.htm. Finally, the results from our investigations regarding the interrelationship between eQTL, ancestry information, and pharmacodynamics provide rich resources for future eQTL studies and practical applications in population genetics and medical genetics.

Nasal Morphology of the Chinese: Three-Dimensional Reference Values for Rhinoplasty

OBJECTIVE

To determine normative nasal measurements for Chinese young adults, conditioned on demographics.

STUDY DESIGN

A cross-sectional descriptive study.

SETTING

A university hospital.

SUBJECTS AND METHODS

Three-dimensional (3D) photographs were captured from 103 Chinese subjects between 18 and 35 years of age using a commercial stereophotographic system. Anthropometric landmarks were identified on these 3D surface images, and measurements suitable for nasal analysis were performed and contrasted against established Caucasian norms. Gender differences in anthropometric dimensions were also analyzed.

RESULTS

Normative data for these measurements are made available. Linear nasal measurements, except those for mid-columella length, were significantly larger in men than in women; further, the nasal tip angle and nasofrontal angle were significantly larger in Chinese women. Contrasts of these new data against published Caucasian norms revealed dimensions that differ for these 2 groups. The Chinese normative mean values for morphological nose width, nasal tip angle, nasofrontal angle, and alar slope angle exceeded those reported for North American Caucasians.

CONCLUSIONS

Gender-specific normative data for the Chinese nose were established in this study to provide a useful tool for surgeons in dealing with rhinoplasty. Moreover, the Chinese nasal anthropometric measurements in this study are broader and flatter than those reported for North American Caucasians.

Genetic associations of relaxin: preterm birth and premature rupture of fetal membranes

OBJECTIVE

Relaxin H2 (RLN2) is a systemic hormone (sRLN) that is produced by the corpus luteum, whereas decidual RLN (dRLN) acts only locally. Elevated sRLN is associated with spontaneous preterm birth (sPTB) and elevated dRLN with preterm premature rupture of membranes (PPROM). Associations were sought between single nucleotide polymorphisms (SNPs) in the RLN2 promoter with levels of dRLN and sRLN in Filipino patients with sPTB, PPROM, or normal term delivery.

STUDY DESIGN

Stringent selection of women with sPTB (n = 20) or PPROM (n = 20) and term control subjects (n = 20) was made from >8000 samples from Filipino patients who delivered at 34-36 weeks' gestation. Twelve SNPs were genotyped on maternal blood, with 9 excluded based on the high linkage disequilibrium or being the same as in the control population. Quantitative immunocytochemistry on parietal decidual tissue was performed (n = 60); sRLN was measured by enzyme-linked immunosorbent assay in a subset of patients (n = 21).

RESULTS

SNP rs4742076 was associated significantly with PPROM (P < .001) and increased expression of dRLN (P < .001). The genotype TT had increased dRLN in PPROM (P < .05). SNP rs3758239 was associated significantly with both PPROM and sPTB (P < .01), and genotype AA had increased dRLN expression (P < .05). The sRLN showed a trend of higher levels in PPROM and sPTB, but was not significant.

CONCLUSION

SNP rs4742076 in the RLN2 promoter was associated with increased dRLN expression and PPROM; SNP rs3758239 was associated with both PPROM and sPTB in these Filipino patients. Specific homozygous genotypes were identified for both SNPs and were shown to be associated with increased dRLN tissue expression.

Normative findings for periocular anthropometric measurements among Chinese young adults in Hong Kong

Measurement of periocular structures is of value in several clinical specialties including ophthalmology, optometry, medical and clinical genetics, oculoplastic surgery, and traumatology. Therefore we aimed to determine the periocular anthropometric norms for Chinese young adults using a noninvasive 3D stereophotography system. Craniofacial images using the 3dMDface system were acquired for 103 Chinese subjects (51 males and 52 females) between the ages of 18 and 35 years. Anthropometric landmarks were identified on these digital images according to standard definitions, and linear distances between these landmarks were calculated. It was found that ocular measurements were significantly larger in Chinese males than females for intercanthal width, biocular width, and eye fissure lengths. No gender differences were found in the eye fissure height and the canthal index which ranged between 43 and 44. Both right and left eye fissure height-length ratios were significantly larger in females. This is the first study to employ 3D stereophotogrammetry to create a database of anthropometric normative data for periocular measurements. These data would be useful for clinical interpretation of periocular pathology and serve as reference values when planning aesthetic and posttraumatic surgical interventions.

GeneOnEarth: fitting genetic PC plots on the globe

Principal component (PC) plots have become widely used to summarize genetic variation of individuals in a sample. The similarity between genetic distance in PC plots and geographical distance has shown to be quite impressive. However, in most situations, individual ancestral origins are not precisely known or they are heterogeneously distributed; hence, they are hardly linked to a geographical area. We have developed GeneOnEarth, a user-friendly web-based tool to help geneticists to understand whether a linear isolation-by-distance model may apply to a genetic data set; thus, genetic distances among a set of individuals resemble geographical distances among their origins. Its main goal is to allow users to first apply a by-view Procrustes method to visually learn whether this model holds. To do that, the user can choose the exact geographical area from an on line 2D or 3D world map by using, respectively, Google Maps or Google Earth, and rotate, flip, and resize the images. GeneOnEarth can also compute the optimal rotation angle using Procrustes analysis and assess statistical evidence of similarity when a different rotation angle has been chosen by the user. An online version of GeneOnEarth is available for testing and using purposes at http://bios.ugr.es/GeneOnEarth.

Distinct HLA allele and haplotype distributions in four ethnic groups of China

Distinct human leukocyte antigen (HLA) allele and haplotype distributions occur in the northern and southern Han populations of China. However, different ethnic groups in China show limited regional distributions for many HLA alleles and haplotypes. Therefore, it is necessary and meaningful to study the differences in HLA allele and haplotype distribution for northern and southern ethnic groups of China. A total of 428 unrelated individuals from the Lisu, Nu, Tu and Yugur ethnic populations were genotyped for HLA-A, -B, -C and -DRB1 alleles using the PCR-Luminex typing method. The frequencies of HLA alleles and statistically inferred haplotypes were calculated. A total of 29 HLA-A, 54 HLA-B, 27 HLA-C and 41 HLA-DRB1 alleles were spread throughout these four populations with distinct allele and deduced haplotype frequencies between populations. Some alleles and deduced haplotypes exhibited significantly different distributions between northern (Tu and Yugur) and southern groups (Lisu and Nu). A phylogenetic tree and principal component analysis were used to compare the HLA polymorphism between our dataset and 19 other eastern and southeastern Asian populations. This analysis showed that Lisu and Nu belong to a cluster of southern ethnic groups, while Tu and Yugur are most closely related to other northern groups. Thus, distinct ethnic population histories were revealed by analyzing HLA allelic polymorphisms with the HLA profiles of the Lisu and Nu southern Chinese ethnic groups clearly different from the Tu and Yugur northern ethnic groups. The results will be useful for future association studies of infectious disease and contribute toward a more efficient search of organ/tissue matches for transplantation.

ETHNOPRED: a novel machine learning method for accurate continental and sub-continental ancestry identification and population stratification correction

BACKGROUND

Population stratification is a systematic difference in allele frequencies between subpopulations. This can lead to spurious association findings in the case-control genome wide association studies (GWASs) used to identify single nucleotide polymorphisms (SNPs) associated with disease-linked phenotypes. Methods such as self-declared ancestry, ancestry informative markers, genomic control, structured association, and principal component analysis are used to assess and correct population stratification but each has limitations. We provide an alternative technique to address population stratification.

RESULTS

We propose a novel machine learning method, ETHNOPRED, which uses the genotype and ethnicity data from the HapMap project to learn ensembles of disjoint decision trees, capable of accurately predicting an individual's continental and sub-continental ancestry. To predict an individual's continental ancestry, ETHNOPRED produced an ensemble of 3 decision trees involving a total of 10 SNPs, with 10-fold cross validation accuracy of 100% using HapMap II dataset. We extended this model to involve 29 disjoint decision trees over 149 SNPs, and showed that this ensemble has an accuracy of ≥ 99.9%, even if some of those 149 SNP values were missing. On an independent dataset, predominantly of Caucasian origin, our continental classifier showed 96.8% accuracy and improved genomic control's λ from 1.22 to 1.11. We next used the HapMap III dataset to learn classifiers to distinguish European subpopulations (North-Western vs. Southern), East Asian subpopulations (Chinese vs. Japanese), African subpopulations (Eastern vs. Western), North American subpopulations (European vs. Chinese vs. African vs. Mexican vs. Indian), and Kenyan subpopulations (Luhya vs. Maasai). In these cases, ETHNOPRED produced ensembles of 3, 39, 21, 11, and 25 disjoint decision trees, respectively involving 31, 502, 526, 242 and 271 SNPs, with 10-fold cross validation accuracy of 86.5% ± 2.4%, 95.6% ± 3.9%, 95.6% ± 2.1%, 98.3% ± 2.0%, and 95.9% ± 1.5%. However, ETHNOPRED was unable to produce a classifier that can accurately distinguish Chinese in Beijing vs. Chinese in Denver.

CONCLUSIONS

ETHNOPRED is a novel technique for producing classifiers that can identify an individual's continental and sub-continental heritage, based on a small number of SNPs. We show that its learned classifiers are simple, cost-efficient, accurate, transparent, flexible, fast, applicable to large scale GWASs, and robust to missing values.

Association of a functional polymorphism in the 3'-untranslated region of SPI1 with systemic lupus erythematosus

OBJECTIVE

SPI1, also referred to as PU.1, is an Ets family transcription factor that interacts with IRF2, IRF4, and IRF8. In view of the significance of the type I interferon pathway in systemic lupus erythematosus (SLE), this study was undertaken to investigate a possible association between SPI1 polymorphisms and SLE.

METHODS

A case-control association study was performed using 6 tag single-nucleotide polymorphisms (SNPs), as well as a SNP located upstream of SPI1 previously found to be associated with acute myelogenous leukemia, in 400 Japanese patients with SLE and 450 healthy controls. Resequencing of all exons and known regulatory regions was performed to identify functional polymorphisms. Association of genotype and SPI1 expression was examined using the GENEVAR database and reporter assays.

RESULTS

A significant association was detected in 2 SNPs in intron 2 (rs10769258 and rs4752829) (P = 0.005 and P = 0.008, respectively, under the dominant model). The association was stronger in patients with nephropathy. Resequencing identified a potentially functional polymorphism in the 3'-untranslated region (3'-UTR), rs1057233, which was in strong linkage disequilibrium with the SNPs in intron 2. The number of risk alleles at rs1057233 was strongly correlated with SPI1 messenger RNA (mRNA) level in the database analysis (P = 0.0002), and was confirmed by a reporter assay. Interestingly, rs1057233 alters a target sequence for microRNA hsa-miR-569 (miR-569). Transfection experiments demonstrated that miR-569 inhibits expression of a reporter construct with the 3'-UTR sequence containing the nonrisk allele but not the risk allele.

CONCLUSION

Our findings indicate that a SNP in the 3'-UTR of SPI1 is associated with elevated SPI1 mRNA level and with susceptibility to SLE.

Sensation-seeking genes and physical activity in youth

Many studies examining genetic influences on physical activity (PA) have evaluated the impact of single nucleotide polymorphisms (SNPs) related to the development of lifestyle-related chronic diseases, under the hypothesis that they would be associated with PA. However, PA is a multidetermined behavior and associated with a multitude of health consequences. Thus, examining a broader range of candidate genes associated with a broader range of PA correlates may provide new insights into the genetic underpinnings of PA. In this study, we focus on one such correlate - sensation-seeking behavior. Participants (N = 1130 Mexican origin youth) provided a saliva sample and data on PA and sensation-seeking tendencies in 2008-2009. Participants were genotyped for 630 functional and tagging variants in the dopamine, serotonin and cannabinoid pathways. Overall 30% of participants (males - 37.6% and females - 22.0%) reported ≥60 min of PA on 5 of 7 days. After adjusting for gender, age and population stratification, and applying the Bayesian False Discovery Probability approach for assessing noteworthiness, four gene variants were significantly associated with PA. In a multivariable model, being male, having higher sensation-seeking tendencies and at least one copy of the minor allele for SNPs in angiotensin I-converting enzyme gene [ACE; rs8066276 odds ratio (OR) = 1.44; P = 0.012] and tryptophan hydroxylase 2 gene (TPH2; rs11615016 OR = 1.73; P = 0.021) were associated with increased likelihood of meeting PA recommendations. Participants with at least one copy of the minor allele for SNPs in synaptosomal-associated protein 25 gene (SNAP25; rs363035 OR = 0.53; P = 0.005) and cannabinoid receptor 1 gene (CNR1; rs6454672 OR = 0.62; P = 0.022) have decreased likelihood of meeting PA recommendations. Our findings extend current knowledge of the complex relationship between PA and possible genetic underpinnings.

The history of human populations in the Japanese Archipelago inferred from genome-wide SNP data with a special reference to the Ainu and the Ryukyuan populations

The Japanese Archipelago stretches over 4000 km from north to south, and is the homeland of the three human populations; the Ainu, the Mainland Japanese and the Ryukyuan. The archeological evidence of human residence on this Archipelago goes back to >30 000 years, and various migration routes and root populations have been proposed. Here, we determined close to one million single-nucleotide polymorphisms (SNPs) for the Ainu and the Ryukyuan, and compared these with existing data sets. This is the first report of these genome-wide SNP data. Major findings are: (1) Recent admixture with the Mainland Japanese was observed for more than one third of the Ainu individuals from principal component analysis and frappe analyses; (2) The Ainu population seems to have experienced admixture with another population, and a combination of two types of admixtures is the unique characteristics of this population; (3) The Ainu and the Ryukyuan are tightly clustered with 100% bootstrap probability followed by the Mainland Japanese in the phylogenetic trees of East Eurasian populations. These results clearly support the dual structure model on the Japanese Archipelago populations, though the origins of the Jomon and the Yayoi people still remain to be solved.

Population pharmacokinetics of modafinil acid and estimation of the metabolic conversion of modafinil into modafinil acid in 5 major ethnic groups of China

AIM

To describe the population pharmacokinetic profile of modafinil acid and to compare the extent of metabolism of modafinil into modafinil acid in 5 major ethnic groups (Han, Mongolian, Korean, Uygur, and Hui) of China.

METHODS

In a multi-center, open-label, single dose clinical trial, 49 healthy volunteers from the 5 ethnic groups received 200 mg of modafinil orally. Blood samples for pharmacokinetic evaluation of modafinil and modafinil acid were drawn before and at different time after the administration. Systematic population pharmacokinetic (PopPK) modeling for modafinil acid was conducted, integrating with our previous PopPK model for modafinil. The influence of ethnicity, gender, height, body weight and body mass index (BMI) was estimated. The extent of metabolism of modafinil into modafinil acid, expressed as the relative conversion fraction, was estimated and compared among the 5 ethnic groups.

RESULTS

When combined with the PopPK model of modafinil, the concentration of modafinil acid versus time profile was best described with a one-compartment model. The typical clearance and volume of distribution for modafinil acid were 4.94 (l/h) and 2.73 (l), respectively. The Korean group had 25% higher clearance, and the Uygur and Hui groups had 12% higher clearance than the Han group. The median for the relative conversion fraction was 0.53 for Koreans, and 0.24 for the other 4 ethnicities.

CONCLUSION

Ethnicity has significant influence on the clearance of modafinil acid. When patients in the 5 ethnic groups are administered drugs or prodrugs catalyzed by esterases and/or amidases, the variability in the extent of drug metabolism should be considered.

Association of common genetic variants with breast cancer risk and clinicopathological characteristics in a Chinese population

Genome-wide association studies (GWAS) have identified various genetic susceptibility loci for breast cancer based mainly on European-ancestry populations. Differing linkage disequilibrium patterns exist between European and Asian populations, and thus GWAS-identified single nucleotide polymorphisms (SNPs) in one population may not be of significance in another population. In order to explore the role of breast cancer susceptibility variants in a Chinese population of Southern Chinese descent, we analyzed 22 SNPs for 1,191 breast cancer cases and 1,534 female controls. Associations between the SNPs and clinicopathological features were also investigated. In addition, we evaluated the combined effects of associated SNPs by constructing risk models. Eight SNPs were associated with an elevated breast cancer risk. Rs2046210/6q25.1 increased breast cancer risk via an additive model [per-allele odds ratio (OR) = 1.43, 95 % confidence interval (CI) = 1.26-1.62], and was associated with estrogen receptor (ER)-positive (per-allele OR = 1.39, 95 % CI = 1.20-1.61) and ER-negative (per-allele OR = 1.55, 95 % CI = 1.28-1.89) disease. Rs2046210 was also associated with stage 1, stage 2, and stage 3 disease, with per-allele ORs of 1.38 (1.14-1.68), 1.48 (1.25-1.74), and 1.58 (1.28-1.94), respectively. Four SNPs mapped to 10q26.13/FGFR2 were associated with increased breast cancer risk via an additive model with per-allelic risks (95 % CI) of 1.26 (1.12-1.43) at rs1219648, 1.22 (1.07-1.38) at rs2981582, 1.21 (1.07-1.36) at rs2981579, and 1.18 (1.04-1.35) at rs11200014. Variants of rs7696175/TLR1, TLR6, rs13281615/8q24, and rs16886165/MAP3K1 were also associated with increased breast cancer risk, with per-allele ORs (95 % CI) of 1.16 (1.00-1.34), 1.15 (1.02-1.29), and 1.15 (1.01-1.29), respectively. Five SNPs associated with breast cancer risk predominantly among ER-positive tumors (rs2981582/FGFR2, rs4415084/MRPS30, rs1219648/FGFR2, rs2981579/FGFR2, and rs11200014/FGFR2). Among our Chinese population, the risk of developing breast cancer increased by 90 % for those with a combination of 6 or more risk alleles, compared to patients with ≤3 risk alleles.

A quantitative comparison of the similarity between genes and geography in worldwide human populations

Multivariate statistical techniques such as principal components analysis (PCA) and multidimensional scaling (MDS) have been widely used to summarize the structure of human genetic variation, often in easily visualized two-dimensional maps. Many recent studies have reported similarity between geographic maps of population locations and MDS or PCA maps of genetic variation inferred from single-nucleotide polymorphisms (SNPs). However, this similarity has been evident primarily in a qualitative sense; and, because different multivariate techniques and marker sets have been used in different studies, it has not been possible to formally compare genetic variation datasets in terms of their levels of similarity with geography. In this study, using genome-wide SNP data from 128 populations worldwide, we perform a systematic analysis to quantitatively evaluate the similarity of genes and geography in different geographic regions. For each of a series of regions, we apply a Procrustes analysis approach to find an optimal transformation that maximizes the similarity between PCA maps of genetic variation and geographic maps of population locations. We consider examples in Europe, Sub-Saharan Africa, Asia, East Asia, and Central/South Asia, as well as in a worldwide sample, finding that significant similarity between genes and geography exists in general at different geographic levels. The similarity is highest in our examples for Asia and, once highly distinctive populations have been removed, Sub-Saharan Africa. Our results provide a quantitative assessment of the geographic structure of human genetic variation worldwide, supporting the view that geography plays a strong role in giving rise to human population structure.

The spatial pattern of human genetic variation provides a basis for investigating the history of human migrations. Statistical techniques such as principal components analysis (PCA) and multidimensional scaling (MDS) have been used to summarize spatial patterns of genetic variation, typically by placing individuals on a two-dimensional map in such a way that pairwise Euclidean distances between individuals on the map approximately reflect corresponding genetic relationships. Although similarity between these statistical maps of genetic variation and the geographic maps of sampling locations is often observed, it has not been assessed systematically across different parts of the world. In this study, we combine genome-wide SNP data from more than 100 populations worldwide to perform a formal comparison between genes and geography in different regions. By examining a worldwide sample and samples from Europe, Sub-Saharan Africa, Asia, East Asia, and Central/South Asia, we find that significant similarity between genes and geography exists in general in different geographic regions and at different geographic levels. Surprisingly, the highest similarity is found in Asia, even though the geographic barrier of the Himalaya Mountains has created a discontinuity on the PCA map of genetic variation.

Integrative analysis of single nucleotide polymorphisms and gene expression efficiently distinguishes samples from closely related ethnic populations

Background

Ancestry informative markers (AIMs) are a type of genetic marker that is informative for tracing the ancestral ethnicity of individuals. Application of AIMs has gained substantial attention in population genetics, forensic sciences, and medical genetics. Single nucleotide polymorphisms (SNPs), the materials of AIMs, are useful for classifying individuals from distinct continental origins but cannot discriminate individuals with subtle genetic differences from closely related ancestral lineages. Proof-of-principle studies have shown that gene expression (GE) also is a heritable human variation that exhibits differential intensity distributions among ethnic groups. GE supplies ethnic information supplemental to SNPs; this motivated us to integrate SNP and GE markers to construct AIM panels with a reduced number of required markers and provide high accuracy in ancestry inference. Few studies in the literature have considered GE in this aspect, and none have integrated SNP and GE markers to aid classification of samples from closely related ethnic populations.

Results

We integrated a forward variable selection procedure into flexible discriminant analysis to identify key SNP and/or GE markers with the highest cross-validation prediction accuracy. By analyzing genome-wide SNP and/or GE markers in 210 independent samples from four ethnic groups in the HapMap II Project, we found that average testing accuracies for a majority of classification analyses were quite high, except for SNP-only analyses that were performed to discern study samples containing individuals from two close Asian populations. The average testing accuracies ranged from 0.53 to 0.79 for SNP-only analyses and increased to around 0.90 when GE markers were integrated together with SNP markers for the classification of samples from closely related Asian populations. Compared to GE-only analyses, integrative analyses of SNP and GE markers showed comparable testing accuracies and a reduced number of selected markers in AIM panels.

Conclusions

Integrative analysis of SNP and GE markers provides high-accuracy and/or cost-effective classification results for assigning samples from closely related or distantly related ancestral lineages to their original ancestral populations. User-friendly BIASLESS (Biomarkers Identification and Samples Subdivision) software was developed as an efficient tool for selecting key SNP and/or GE markers and then building models for sample subdivision. BIASLESS was programmed in R and R-GUI and is available online at http://www.stat.sinica.edu.tw/hsinchou/genetics/prediction/BIASLESS.htm.

The genetic association of DUSP6 with bipolar disorder and its effect on ERK activity

The dual-specificity phosphatase 6 (DUSP6) gene resides at chromosome location 12q22-23, which is one of the candidate loci for susceptibility to bipolar disorder and which encodes a phosphatase selective for extracellular signal-regulated kinase (ERK). Previously, we reported a positive association between the functional Leu114Val polymorphism (rs2279574) in DUSP6 and bipolar disorder. Given that the association between DUSP6 and the reported down-regulation of DUSP6 transcript in bipolar postmortem brains were sex-dimorphic, showing significance in women but not men, we performed two independent analyses in homogenous samples of male and female Korean patients with bipolar disorder or schizophrenia using samples enlarged from our previous report. Among the examined DUSP6 SNPs, five (rs769700, rs704076, rs770087, rs808820, and rs2279574) showed positive allelic associations, with the frequency of minor alleles (C, T, G, G, and G) in each SNP significantly increased in women with BD. Consequently, the "C-T-G-G-G" haplotype was significantly over-represented (P=0.016; OR=3.242), whereas the "T-G-T-A-T" haplotype was significantly under-represented (P=0.014; OR=0.697). We found no significant associations with DUSP6 SNPs in men with bipolar disorder or schizophrenia. We also investigated the functions of the functional SNPs' positive associations and found that Leu114Val (rs2279574; T/G) and Ser144Ala (rs770087; T/G) mutations in DUSP6 proteins reduced lithium-induced ERK1/2 phosphorylation in vitro, implicating the dominant active functions. Thus, DUSP6 may not only play important roles in the pathogenesis of bipolar disorder, particularly in women, but also affect the therapeutic response to lithium through modulating lithium's effects on intracellular signaling.