Natural positive selection and north-south genetic diversity in East Asia

Noninvasive Prenatal Diagnosis of SEA-Thalassemia by Combining 1000 Genomes Database and Relative Haplotype Dosage

To explore a noninvasive method for diagnosis of SEA-thalassemia and to investigate whether the regional factors affect the accuracy of this method. The method involved using a public database and bioinformatics software to construct parental haplotypes for proband and predicting fetal genotypes using relative haplotype dosage. We screened and downloaded sequencing data of couples who were both SEA-thalassemia carriers from the China National Genebank public data platform, and matched the sequencing data format with that of the reference panel using Ubuntu system tools. We then used Beagle software to construct parental haplotypes, predicted fetal haplotypes by relative haplotype dosage. Finally, we used Hidden Markov Model and Viterbi algorithm to determine fetal pathogenic haplotypes. All noninvasive fetal genotype diagnosis results were compared with gold standard gap-PCR electrophoresis results. Our method was successful in diagnosing 13 families with SEA-thalassemia carriers. The best diagnostic results were obtained when Southern Chinese Han was used as the reference panel, and 10 families showed full agreement between our noninvasive diagnostic results and the gap-PCR electrophoresis results. The accuracy of our method was higher when using a Chinese Han as the reference panel for haplotype construction in the Southern Chinese Han region as opposed to Beijing Chinese region. The combined use of public databases and relative haplotype dosage for diagnosing SEA-thalassemia is a feasible approach. Our method produces the best noninvasive diagnostic results when the test samples and population reference panel are closely matched in both ethnicity and geography. When constructing parental haplotypes with our method, it is important to consider the effect of region in addition to population background alone.

High-resolution HLA genotyping identifies risk alleles in both class I and II for primary autoimmune neutropenia in early childhood in a Danish cohort

HLA studies in patients with autoimmune neutropenia (AIN) have shown very consistent results for the association with HLA class II alleles at low resolution. This study aimed to examine the association of both HLA class I and class II at high resolution to clarify the contribution of risk alleles to the disease. A total of 107 AIN patients were genotyped for six loci of HLA class I (HLA-A, -B and -C) and class II (HLA-DRB1, -DQB1, and -DPB1) genes by a high-resolution (3-field, 6-digit) analysis and compared with HLA typing of 1000 healthy controls. Compared with the controls, the allele frequencies were significantly higher in AIN patients for A*02:17:01G, C*01:02:01G, DRB1*10:01:01G, DRB1*14:01:01G, DRB1*16:01:01G, DQB1*05:02:01G, and DQB1*05:03:01G but lower significant for C*03:04:01G, DRB1*04:01:01G, DRB1*13:02:01G, DQB1*03:02:01G, and DQB1*06:04:01G. Frequently associated two-locus haplotypes were found to be DRB1*10:01:01G-DQB1*05:01:01G and DRB1*16:01:01G-DQB1*05:02:01G, while the S2 (Q- or D-KRAA) shared epitope (SE) was associated with lower risk. A unique association with HLA alleles was observed between patients with specific anti-HNA-1a antibodies and broad-reacting anti-FcγRIIIb. Anti-HNA-1a antibody-positive patients were associated with C*01:02:01G, DRB1*01:01:01G, DRB1*16:01:01G, DQB1*05:01:01G, DQB1*05:02:01G, DQB1*06:04:01G, and DPB1*10:01:01G; the two-locus haplotypes DRB1*01:01:01G-DQB1*05:01:01G and DRB1*16:01:01G-DQB1*05:02:01G; and the S3P (Q- or R-RRAA) SE. Anti-FcγRIIIb antibody-positive patients were associated with the alleles A*02:17:01G, DRB1*10:01:01G, and DQB1*05:02:01G; the haplotypes DRB1*10:01:01G-DQB1*05:01:01G and DRB1*11:01:02G-DQB1*05:02:01G; and the S3D (DRRAA) SE. The different associations regarding FcγRIIIb antibody specificities could indicate disease heterogeneity.

Identifying signatures of natural selection in Indian populations

In this study, we present the results of a genome-wide scan for signatures of positive selection using data from four tribal groups (Kokana, Warli, Bhil, and Pawara) and two caste groups (Deshastha Brahmin and Kunbi Maratha) from West of the Maharashtra State In India, as well as two samples of South Asian ancestry from the 1KG project (Gujarati Indian from Houston, Texas and Indian Telugu from UK). We used an outlier approach based on different statistics, including PBS, xpEHH, iHS, CLR, Tajima's D, as well as two recently developed methods: Graph-aware Retrieval of Selective Sweeps (GRoSS) and Ascertained Sequentially Markovian Coalescent (ASMC). In order to minimize the risk of false positives, we selected regions that are outliers in all the samples included in the study using more than one method. We identified putative selection signals in 107 regions encompassing 434 genes. Many of the regions overlap with only one gene. The signals observed using microarray-based data are very consistent with our analyses using high-coverage sequencing data, as well as those identified with a novel coalescence-based method (ASMC). Importantly, at least 24 of these genomic regions have been identified in previous selection scans in South Asian populations or in other population groups. Our study highlights genomic regions that may have played a role in the adaptation of anatomically modern humans to novel environmental conditions after the out of Africa migration.

HLA-Cw1 and Psoriasis

Psoriasis is a chronic inflammatory skin condition with regional and ethnic differences in its prevalence and clinical manifestations. Human leukocyte antigen (HLA)-Cw6 is the disease allele conferring the greatest risk to psoriasis, but its prevalence is lower in Asian individuals. Recent studies have found associations between HLA-Cw1 and some Asian populations with psoriasis, especially Southern Chinese. HLA-Cw6 was associated with type I early-onset psoriasis, guttate psoriasis, Koebner phenomenon, and better response to methotrexate, interleukin (IL)-12/23, IL-17, and IL-23 targeting drugs. In contrast, HLA-Cw1 positivity has been associated with erythrodermic psoriasis, pustular psoriasis, and the axial type of psoriatic arthritis. Furthermore, HLA-Cw1 was more frequently associated with high-need patients who did not respond to conventional therapies. No known trigger factor nor autoantigen has been identified for HLA-Cw1 positivity. However, HLA-Cw1 has been linked to some viral agents. For example, cytotoxic T lymphocytes recognize multiple cytomegalovirus pp65-derived epitopes presented by HLA alleles, including HLA-C*01:02. In addition, cytomegalovirus can lead to severe exacerbation of psoriatic skin disease. The proposed interaction between viral infection, HLA-Cw1, and psoriasis is through the killer cell immunoglobulin-like receptors of natural killer cells. Given the diverse nature of psoriasis pathogenesis and the difference in HLA-Cw prevalence in different racial groups, more studies are needed to confirm the role of HLA-Cw1 in psoriasis.

Tumor infiltrating lymphocyte signature is associated with single nucleotide polymorphisms and predicts survival in esophageal squamous cell carcinoma patients

Purpose: Esophageal cancer is the sixth leading cause of cancer-related death worldwide, and is associated with a poor prognosis. Stromal tumor infiltrating lymphocytes (sTIL) and certain single nucleotide polymorphisms (SNPs) have been found to be predictive of patient survival. In this study, we explored the association between SNPs and sTIL regarding the predictability of disease-free survival in patients with esophageal squamous cell carcinoma (ESCC).

Materials and methods: We collected 969 pathologically confirmed ESCC patients from 2010 to 2013 and genotyped 101 SNPs from 59 genes. The number of sTIL for each patient was determined using an automatic algorithm. A Kruskal-Wallis test was used to determine the association between genotype and sTIL. The genotypes and clinical factors related to survival were analyzed using a Kaplan-Meier curve, Cox proportional hazards model, and log-rank test.

Results: The median age of the patients was 67 (42-85 years), there was a median follow-up of 851.5 days and 586 patients died. The univariable analysis showed that 10 of the 101 SNPs were associated with sTIL. Six SNPs were also associated with disease-free survival. A multivariable analysis revealed that sTIL, rs1801131, rs25487, and rs8030672 were independent prognostic markers for ESCC patients. The model combining SNPs, clinical characteristics and sTIL outperformed the model with clinical characteristics alone for predicting outcomes in ESCC patients.

Conclusion: We discovered 10 SNPs associated with sTIL in ESCC and we built a model of sTIL, SNPs and clinical characteristics with improved prediction of survival in ESCC patients.

Limited Evidence for Selection at the FADS Locus in Native American Populations

The FADS locus contains the genes FADS1 and FADS2 that encode enzymes involved in the synthesis of long-chain polyunsaturated fatty acids. This locus appears to have been a repeated target of selection in human evolution, likely because dietary input of long-chain polyunsaturated fatty acids varied over time depending on environment and subsistence strategy. Several recent studies have identified selection at the FADS locus in Native American populations, interpreted as evidence for adaptation during or subsequent to the passage through Beringia. Here, we show that these signals are confounded by independent selection—postdating the split from Native Americans—in the European and, possibly, the East Asian populations used in the population branch statistic test. This is supported by direct evidence from ancient DNA that one of the putatively selected haplotypes was already common in Northern Eurasia at the time of the separation of Native American ancestors. An explanation for the present-day distribution of the haplotype that is more consistent with the data is that Native Americans retain the ancestral state of Paleolithic Eurasians. Another haplotype at the locus may reflect a secondary selection signal, although its functional impact is unknown.

Genome-wide analyses reveal footprints of divergent selection and popping-related traits in CIMMYT's maize inbred lines

Popcorn (Zea mays L. var. Everta) is the most ancient type of cultivated maize. However, there is little known about the genetics of popping-related traits based on genotyping-by-sequencing (GBS) technology. Here, we characterized the phenotypic variation for seven popping-related traits in maize kernels among 526 CIMMYT inbred lines (CMLs). In total, 155 083 high-quality single nucleotide polymorphism (SNP) markers were identified by a GBS approach. Several trait-associated loci were detected by genome-wide association study for color, popping expansion volume, shape, pericarp, flotation index, floury/vitreous, and protein content, explaining a majority of the observed phenotypic variance, and these were validated by a diverse panel comprising 764 tropical landrace accessions. Sixty two of the identified loci were recognized to have undergone selection. On average, there was a 55.27% frequency for alleles that promote popping in CMLs. Our work not only pinpoints previously unknown loci for popping-related traits, but also reveals that many of these loci have undergone selection. Beyond establishing a new benchmark for the genetics of popcorn, our study provides a foundation for gene discovery and breeding. It also presents evidence to investigate the role of a gradual loss of popping ability as a by-product of diversification of culinary uses throughout the evolution of teosinte-to-modern maize.

Learning the properties of adaptive regions with functional data analysis

Identifying regions of positive selection in genomic data remains a challenge in population genetics. Most current approaches rely on comparing values of summary statistics calculated in windows. We present an approach termed SURFDAWave, which translates measures of genetic diversity calculated in genomic windows to functional data. By transforming our discrete data points to be outputs of continuous functions defined over genomic space, we are able to learn the features of these functions that signify selection. This enables us to confidently identify complex modes of natural selection, including adaptive introgression. We are also able to predict important selection parameters that are responsible for shaping the inferred selection events. By applying our model to human population-genomic data, we recapitulate previously identified regions of selective sweeps, such as OCA2 in Europeans, and predict that its beneficial mutation reached a frequency of 0.02 before it swept 1,802 generations ago, a time when humans were relatively new to Europe. In addition, we identify BNC2 in Europeans as a target of adaptive introgression, and predict that it harbors a beneficial mutation that arose in an archaic human population that split from modern humans within the hypothesized modern human-Neanderthal divergence range.

Whole Mitochondrial DNA Sequencing Analysis in 47 Han Populations in Southwest China

Background

Mitochondrial DNA (mtDNA) sequencing has been used in many areas, including forensic genetics. Due to the rapid development of sequencing technology, whole mtDNA sequencing is now possible and may be used in epidemiological and forensic studies. This study aimed to use whole mtDNA sequencing to investigate 47 Chongqing Han populations in southwest China and the diversity in the mtGenome reference data.

Material/Methods

The mtDNA of 47 Chongqing Han populations was generated using the Ion Torrent Personal Genome Machine (PGM) system. The extent of the effects of the mtDNA on the subpopulations was investigated and compared with six other populations from published studies. Pairwise fixation index (FST), a measure of population differentiation due to genetic structure, were calculated. Analysis of molecular variance (AMOVA) was performed, and 1257 hypervariable region data sets were added to the principal component analysis (PCA).

Results

The whole mtDNA sequencing data of 47 southwest Chinese Han populations were successfully recovered. Expanding the sequencing rage increased the discrimination power of mtDNA from three-times to 25-times based on different populations. The subpopulation effects showed 20 times the differences in match probability when compared with south China regions.

Conclusions

Whole mtDNA sequencing distinguished between individuals from 47 Chongqing Han populations in southwest China and has potential applications that include high-quality forensic identification.

A Comprehensive Map of Genetic Variation in the World's Largest Ethnic Group-Han Chinese

As are most non-European populations, the Han Chinese are relatively understudied in population and medical genetics studies. From low-coverage whole-genome sequencing of 11,670 Han Chinese women we present a catalog of 25,057,223 variants, including 548,401 novel variants that are seen at least 10 times in our data set. Individuals from this data set came from 24 out of 33 administrative divisions across China (including 19 provinces, 4 municipalities, and 1 autonomous region), thus allowing us to study population structure, genetic ancestry, and local adaptation in Han Chinese. We identified previously unrecognized population structure along the East-West axis of China, demonstrated a general pattern of isolation-by-distance among Han Chinese, and reported unique regional signals of admixture, such as European influences among the Northwestern provinces of China. Furthermore, we identified a number of highly differentiated, putatively adaptive, loci (e.g., MTHFR, ADH7, and FADS, among others) that may be driven by immune response, climate, and diet in the Han Chinese. Finally, we have made available allele frequency estimates stratified by administrative divisions across China in the Geography of Genetic Variant browser for the broader community. By leveraging the largest currently available genetic data set for Han Chinese, we have gained insights into the history and population structure of the world's largest ethnic group.

Architecture of polymorphisms in the human genome reveals functionally important and positively selected variants in immune response and drug transporter genes

Background

Genetic polymorphisms can contribute to phenotypic differences amongst individuals, including disease risk and drug response. Characterization of genetic polymorphisms that modulate gene expression and/or protein function may facilitate the identification of the causal variants. Here, we present the architecture of genetic polymorphisms in the human genome focusing on those predicted to be potentially functional/under natural selection and the pathways that they reside.

Results

In the human genome, polymorphisms that directly affect protein sequences and potentially affect function are the most constrained variants with the lowest single-nucleotide variant (SNV) density, least population differentiation and most significant enrichment of rare alleles. SNVs which potentially alter various regulatory sites, e.g. splicing regulatory elements, are also generally under negative selection.

Interestingly, genes that regulate the expression of transcription/splicing factors and histones are conserved as a higher proportion of these genes is non-polymorphic, contain ultra-conserved elements (UCEs) and/or has no non-synonymous SNVs (nsSNVs)/coding INDELs. On the other hand, major histocompatibility complex (MHC) genes are the most polymorphic with SNVs potentially affecting the binding of transcription/splicing factors and microRNAs (miRNA) exhibiting recent positive selection (RPS). The drug transporter genes carry the most number of potentially deleterious nsSNVs and exhibit signatures of RPS and/or population differentiation. These observations suggest that genes that interact with the environment are highly polymorphic and targeted by RPS.

Conclusions

In conclusion, selective constraints are observed in coding regions, master regulator genes, and potentially functional SNVs. In contrast, genes that modulate response to the environment are highly polymorphic and under positive selection.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0175-1) contains supplementary material, which is available to authorized users.

Genetic resolution of applied biosystems™ precision ID Ancestry panel for seven Asian populations

Massively parallel sequencing (MPS) offers additional information in cases that lack reference samples for comparison. The HID-Ion AmpliSeq Ancestry Panel is a forensic multiplex platform consisting of 165 autosomal markers designed to provide biogeographic ancestry information. We analyzed seven different population groups from Asia to assess the accuracy and reliability of analysis, using this panel. In this study, we have designated 750 unrelated Asians, from southern China (n = 99), Beijing (n = 100), Japan (n = 101), Korea (n = 100), Vietnam (n = 100), Nepal (n = 100), India (n = 51), and Pakistan (n = 99). The likelihood ratios of 750 Asians were calculated, using the Torrent Server and the HID SNP Genotyper Plugin Version 4.3.2. The results reveal that all Northeast Asians (China, Japan, and Korea), and Vietnamese, were predicted as East Asians with the highest ethnicity likelihood values from reference data. Although the samples from Nepal, India and Pakistan (Southwest Asians), were predominantly predicted as South Asians, there were seven cases of which results revealed as Europeans, with similar likelihood patterns. Population structure analysis indicated that Southwest Asians have a genetic profile that is distinguishable from those of other Asian populations. This panel was validated for potential usefulness in forensic analysis and the findings could be used as a basis for building additional Asian specific reference databases. Nevertheless, additional analyses comprising larger sample sizes will be necessary, especially Southeast Asians, to fully understand the Asian population structure, and to discriminate them with further details.

Recurrent deletions of the X chromosome linked CNV64, CNV67, and CNV69 shows geographic differences across China and no association with idiopathic infertility in men

A recent study found that three recurrent deletions of X chromosome linked copy number variations (CNVs), CNV64, CNV67 and CNV69 were associated with idiopathic male infertility in Spanish and Italian populations, especially CNV67 resembling the azoospermia factor deletions. That merits further investigations among different populations. This study was conducted to examine the prevalence of the three CNVs deletions and their associations with idiopathic male infertility in Chinese Han population. The present study included a large population of 1550 Chinese Han subjects recruited between 2014 and 2016. In total, 714 infertile participants were diagnosed as idiopathic infertility with different conditions (288 with non-obstructive azoospermia, 210 oligozoospermia and 216 asthenospermia) and 836 fertile participants (vasectomized men). The fertile participants were recruited from the representative areas: the north (Hebei and Shanxi), center (Hubei and Jiangsu), and south (Guangdong) of China. All patients were recruited from Hubei province. A multiplex PCR system was established to screen the deletion of the three CNVs, and deletion was confirmed by general PCR. Similar rates of these deletions were observed in infertile men and fertile participants (Hubei), and among the different conditions of infertility. Moreover, CNV64 and CNV67 map distribution geographically differed across China. The three CNVs in fertile groups of other regions were similar, except for Guangdong. No association between the three CNVs deletions and idiopathic male infertility was observed. CNV67 is rare in central China, albeit large sample size study for confirmation is warranted. It seems that the association between these CNVs deletions and idiopathic male infertility is ethnic dependent. There is still need to screen the CNVs deletions in other ethnicities. We suggested to consider the stratification patterns and geographic differences when prescribing CNVs deletions screening as a test in male infertility.

Establishing multiple omics baselines for three Southeast Asian populations in the Singapore Integrative Omics Study

The Singapore Integrative Omics Study provides valuable insights on establishing population reference measurement in 364 Chinese, Malay, and Indian individuals. These measurements include > 2.5 millions genetic variants, 21,649 transcripts expression, 282 lipid species quantification, and 284 clinical, lifestyle, and dietary variables. This concept paper introduces the depth of the data resource, and investigates the extent of ethnic variation at these omics and non-omics biomarkers. It is evident that there are specific biomarkers in each of these platforms to differentiate between the ethnicities, and intra-population analyses suggest that Chinese and Indians are the most biologically homogeneous and heterogeneous, respectively, of the three groups. Consistent patterns of correlations between lipid species also suggest the possibility of lipid tagging to simplify future lipidomics assays. The Singapore Integrative Omics Study is expected to allow the characterization of intra-omic and inter-omic correlations within and across all three ethnic groups through a systems biology approach.The Singapore Genome Variation projects characterized the genetics of Singapore's Chinese, Malay, and Indian populations. The Singapore Integrative Omics Study introduced here goes further in providing multi-omic measurements in individuals from these populations, including genetic, transcriptome, lipidome, and lifestyle data, and will facilitate the study of common diseases in Asian communities.

Chinese SLE Treatment and Research Group (CSTAR) registry:VIII. Influence of socioeconomic and geographical variables on disease phenotype and activity in Chinese patients with SLE

AIM

The aim of this study was to estimate the influence of socioeconomic and geographical variables on disease phenotype and activity of systemic lupus erythematosus (SLE) in a Chinese population.

METHODS

Data on 904 therapy-naïve SLE patients from the Chinese SLE Treatment and Research Group (CSTAR), which is a nation-wide database of SLE patients, were used to evaluate the influence of socioeconomic status (SES) and geographical regions on disease phenotypes. The influence of regional climate and economic variables (i.e., data from the database of the Chinese Statistical Yearly Annals) on disease activity and nephropathy was further analyzed in a multivariate analysis to identify some of the environmental factors contributing to SLE.

RESULTS

Phenotype differences in SLE were found between patients with different SES. In the multivariate regression analysis that considered regional economic and climate factors, after adjusting for age, gender and disease duration, lower educational levels (β: -0.122, P = 0.001), lower availability of medical technical personnel (MTP) per 1000 people (β: -0.192, P < 0.001), higher annual duration of sunshine (β: 0.178, P = 0.010) and residence in southern China (β: 0.165, P = 0.001) were significantly associated with higher disease activity scores. Living in southern China (odds ratio [OR]: 1.907, P = 0.002) seemed to be a significant risk factor for nephropathy. A high number of MTP per 1000 people (OR: 0.951, P = 0.014) and relatively high temperatures were determined to be protective factors for nephropathy (OR: 0.946, P = 0.011).

CONCLUSIONS

The phenotype pattern and disease activity varied between SLE patients by SES and other regional variables, such as educational level, availability of medical technical personnel, annual sunshine duration and regional temperature.

Characterising private and shared signatures of positive selection in 37 Asian populations

The Asian Diversity Project (ADP) assembled 37 cosmopolitan and ethnic minority populations in Asia that have been densely genotyped across over half a million markers to study patterns of genetic diversity and positive natural selection. We performed population structure analyses of the ADP populations and divided these populations into four major groups based on their genographic information. By applying a highly sensitive algorithm haploPS to locate genomic signatures of positive selection, 140 distinct genomic regions exhibiting evidence of positive selection in at least one population were identified. We examined the extent of signal sharing for regions that were selected in multiple populations and observed that populations clustered in a similar fashion to that of how the ancestry clades were phylogenetically defined. In particular, populations predominantly located in South Asia underwent considerably different adaptation as compared with populations from the other geographical regions. Signatures of positive selection present in multiple geographical regions were predicted to be older and have emerged prior to the separation of the populations in the different regions. In contrast, selection signals present in a single population group tended to be of lower frequencies and thus can be attributed to recent evolutionary events.

Fast Principal-Component Analysis Reveals Convergent Evolution of ADH1B in Europe and East Asia

Searching for genetic variants with unusual differentiation between subpopulations is an established approach for identifying signals of natural selection. However, existing methods generally require discrete subpopulations. We introduce a method that infers selection using principal components (PCs) by identifying variants whose differentiation along top PCs is significantly greater than the null distribution of genetic drift. To enable the application of this method to large datasets, we developed the FastPCA software, which employs recent advances in random matrix theory to accurately approximate top PCs while reducing time and memory cost from quadratic to linear in the number of individuals, a computational improvement of many orders of magnitude. We apply FastPCA to a cohort of 54,734 European Americans, identifying 5 distinct subpopulations spanning the top 4 PCs. Using the PC-based test for natural selection, we replicate previously known selected loci and identify three new genome-wide significant signals of selection, including selection in Europeans at ADH1B. The coding variant rs1229984(∗)T has previously been associated to a decreased risk of alcoholism and shown to be under selection in East Asians; we show that it is a rare example of independent evolution on two continents. We also detect selection signals at IGFBP3 and IGH, which have also previously been associated to human disease.

Mapping the genetic diversity of HLA haplotypes in the Japanese populations

Japan has often been viewed as an Asian country that possesses a genetically homogenous community. The basis for partitioning the country into prefectures has largely been geographical, although cultural and linguistic differences still exist between some of the districts/prefectures, especially between Okinawa and the mainland prefectures. The Major Histocompatibility Complex (MHC) region has consistently emerged as the most polymorphic region in the human genome, harbouring numerous biologically important variants; nevertheless the presence of population-specific long haplotypes hinders the imputation of SNPs and classical HLA alleles. Here, we examined the extent of genetic variation at the MHC between eight Japanese populations sampled from Okinawa, and six other prefectures located in or close to the mainland of Japan, specifically focusing at the haplotypes observed within each population, and what the impact of any variation has on imputation. Our results indicated that Okinawa was genetically farther to the mainland Japanese than were Gujarati Indians from Tamil Indians, while the mainland Japanese from six prefectures were more homogeneous than between northern and southern Han Chinese. The distribution of haplotypes across Japan was similar, although imputation was most accurate for Okinawa and several mainland prefectures when population-specific panels were used as reference.

Fifteen years of genomewide scans for selection: trends, lessons and unaddressed genetic sources of complication

Genomewide scans for natural selection (GWSS) have become increasingly common over the last 15 years due to increased availability of genome-scale genetic data. Here, we report a representative survey of GWSS from 1999 to present and find that (i) between 1999 and 2009, 35 of 49 (71%) GWSS focused on human, while from 2010 to present, only 38 of 83 (46%) of GWSS focused on human, indicating increased focus on nonmodel organisms; (ii) the large majority of GWSS incorporate interpopulation or interspecific comparisons using, for example F(ST), cross-population extended haplotype homozygosity or the ratio of nonsynonymous to synonymous substitutions; (iii) most GWSS focus on detection of directional selection rather than other modes such as balancing selection; and (iv) in human GWSS, there is a clear shift after 2004 from microsatellite markers to dense SNP data. A survey of GWSS meant to identify loci positively selected in response to severe hypoxic conditions support an approach to GWSS in which a list of a priori candidate genes based on potential selective pressures are used to filter the list of significant hits a posteriori. We also discuss four frequently ignored determinants of genomic heterogeneity that complicate GWSS: mutation, recombination, selection and the genetic architecture of adaptive traits. We recommend that GWSS methodology should better incorporate aspects of genomewide heterogeneity using empirical estimates of relevant parameters and/or realistic, whole-chromosome simulations to improve interpretation of GWSS results. Finally, we argue that knowledge of potential selective agents improves interpretation of GWSS results and that new methods focused on correlations between environmental variables and genetic variation can help automate this approach.

Dissecting ancestry genomic background in substance dependence genome-wide association studies

AIMS

To understand the role of ancestral genomic background in substance dependence (SD) genome-wide association studies (GWAS), we analyzed population diversity at genetic loci associated with SD traits and evaluated its effect on GWAS outcomes.

MATERIALS & METHODS

We investigated 24 genes with variants associated with SD by GWAS; and 82 loci with putative subordinate roles with respect to SD-associated genes.

RESULTS

We observed high ancestry-related frequency differences in common functional alleles in GWAS relevant genes and their interactive partners. Common functional alleles with high frequency differences demonstrated significant effects on the GWAS outcomes.

CONCLUSION

Population differences in SD GWAS outcomes seem not to be influenced by general variation across the genome, but by ancestry-related local haplotype structures at SD-associated loci.

Revisiting classic clines in Drosophila melanogaster in the age of genomics

Adaptation to spatially varying environments has been studied for decades, but advances in sequencing technology are now enabling researchers to investigate the landscape of genetic variation underlying this adaptation genome wide. In this review we highlight some of the decades-long research on local adaptation in Drosophila melanogaster from well-studied clines in North America and Australia. We explore the evidence for parallel adaptation and identify commonalities in the genes responding to clinal selection across continents as well as discussing instances where patterns differ among clines. We also investigate recent studies utilizing whole-genome data to identify clines in D. melanogaster and several other systems. Although connecting segregating genomic variation to variation in phenotypes and fitness remains challenging, clinal genomics is poised to increase our understanding of local adaptation and the selective pressures that drive the extensive phenotypic diversity observed in nature.

Fine-scale population structure of Malays in Peninsular Malaysia and Singapore and implications for association studies

Fine scale population structure of Malays - the major population in Malaysia, has not been well studied. This may have important implications for both evolutionary and medical studies. Here, we investigated the population sub-structure of Malay involving 431 samples collected from all states from peninsular Malaysia and Singapore. We identified two major clusters of individuals corresponding to the north and south peninsular Malaysia. On an even finer scale, the genetic coordinates of the geographical Malay populations are in correlation with the latitudes (R² = 0.3925; P = 0.029). This finding is further supported by the pairwise F_ST of Malay sub-populations, of which the north and south regions showed the highest differentiation (F_{ST [North–south]} = 0.0011). The collective findings therefore suggest that population sub-structure of Malays are more heterogenous than previously expected even within a small geographical region, possibly due to factors like different genetic origins, geographical isolation, could result in spurious association as demonstrated in our analysis. We suggest that cautions should be taken during the stage of study design or interpreting the association signals in disease mapping studies which are expected to be conducted in Malay population in the near future.

Inherited metabolic diseases in the Southern Chinese population: spectrum of diseases and estimated incidence from recurrent mutations

Inherited metabolic diseases (IMDs) are a large group of rare genetic diseases. The spectrum and incidences of IMDs differ among populations, which has been well characterised in Caucasians but much less so in Chinese. In a setting of a University Hospital Metabolic Clinic in Hong Kong, over 100 patients with IMDs have been seen during a period of 13 years (from 1997 to 2010). The data were used to define the spectrum of diseases in the Southern Chinese population. Comparison with other populations revealed a unique spectrum of common IMDs. Furthermore, the incidence of the common IMDs was estimated by using population carrier frequencies of known recurrent mutations. Locally common diseases (their estimated incidence) include (1) glutaric aciduria type 1 (∼1/60,000), (2) multiple carboxylase deficiency (∼1/60,000), (3) primary carnitine deficiency (∼1/60,000), (4) carnitine-acylcarnitine translocase deficiency (∼1/60,000), (5) glutaric aciduria type 2 (∼1/22,500), (6) citrin deficiency (∼1/17,000), (7) tetrahydrobiopterin-deficient hyperphenylalaninaemia due to 6-pyruvoyl-tetrahydropterin synthase deficiency (∼1/60,000), (8) glycogen storage disease type 1 (∼1/150,000). In addition, ornithine carbamoyltransferase deficiency and X-linked adrenoleukodystrophy are common X-linked diseases. Findings of the disease spectrum and treatment outcome are summarised here which may be useful for clinical practice. In addition, data will also be useful for policy makers in planning of newborn screening programs and resource allocation.

Characterizing the genetic differences between two distinct migrant groups from Indo-European and Dravidian speaking populations in India

Background

India is home to many ethnically and linguistically diverse populations. It is hypothesized that history of invasions by people from Persia and Central Asia, who are referred as Aryans in Hindu Holy Scriptures, had a defining role in shaping the Indian population canvas. A shift in spoken languages from Dravidian languages to Indo-European languages around 1500 B.C. is central to the Aryan Invasion Theory. Here we investigate the genetic differences between two sub-populations of India consisting of: (1) The Indo-European language speaking Gujarati Indians with genome-wide data from the International HapMap Project; and (2) the Dravidian language speaking Tamil Indians with genome-wide data from the Singapore Genome Variation Project.

Results

We implemented three population genetics measures to identify genomic regions that are significantly differentiated between the two Indian populations originating from the north and south of India. These measures singled out genomic regions with: (i) SNPs exhibiting significant variation in allele frequencies in the two Indian populations; and (ii) differential signals of positive natural selection as quantified by the integrated haplotype score (iHS) and cross-population extended haplotype homozygosity (XP-EHH). One of the regions that emerged spans the SLC24A5 gene that has been functionally shown to affect skin pigmentation, with a higher degree of genetic sharing between Gujarati Indians and Europeans.

Conclusions

Our finding points to a gene-flow from Europe to north India that provides an explanation for the lighter skin tones present in North Indians in comparison to South Indians.

Evaluating the possibility of detecting evidence of positive selection across Asia with sparse genotype data from the HUGO Pan-Asian SNP Consortium

Background

The HUGO Pan-Asian SNP Consortium (PASNP) has generated a genetic resource of almost 55,000 autosomal single nucleotide polymorphisms (SNPs) across more than 1,800 individuals from 73 urban and indigenous populations in Asia. This has offered valuable insights into the correlation between the genetic ancestry of these populations with major linguistic systems and geography. Here, we attempt to understand whether adaptation to local climate, diet and environment partly explains the genetic variation present in these populations by investigating the genomic signatures of positive selection.

Results

To evaluate the impact to the selection analyses due to the considerably lower SNP density as compared to other population genetics resources such as the International HapMap Project (HapMap) or the Singapore Genome Variation Project, we evaluated the extent of haplotype phasing switch errors and the consistency of selection signals from three haplotype-based approaches (iHS, XP-EHH, haploPS) when the HapMap data is thinned to a similar density as PASNP. We subsequently applied haploPS to detect and characterize positive selection in the PASNP populations, identifying 59 genomics regions that were selected in at least one PASNP populations. A cluster analysis on the basis of these 59 signals showed that indigenous populations such as the Negrito from Malaysia and Philippines, the China Hmong, and the Taiwan Ami and Atayal shared more of these signals. We also reported evidence of a positive selection signal encompassing the beta globin gene in the Taiwan Ami and Atayal that was distinct from the signal in the HapMap Africans, suggesting the possibility of convergent evolution at this locus due to malarial selection.

Conclusions

We established that the lower SNP content of the PASNP data conferred weaker ability to detect signatures of positive selection, but the availability of the new approach haploPS retained modest power. Out of all the populations in PASNP, we identified only 59 signals, suggesting a strong need for high-density population-level genotyping data or sequencing data in order to achieve a comprehensive survey of positive selection in Asian populations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-332) contains supplementary material, which is available to authorized users.

HLA variation reveals genetic continuity rather than population group structure in East Asia

Genetic differences between Northeast Asian (NEA) and Southeast Asian (SEA) populations have been observed in numerous studies. At the among-population level, despite a clear north-south differentiation observed for many genetic markers, debates were led between abrupt differences and a continuous pattern. At the within-population level, whether NEA or SEA populations have higher genetic diversity is also highly controversial. In this study, we analyzed a large set of HLA data from East Asia in order to map the genetic variation among and within populations in this continent and to clarify the distribution pattern of HLA lineages and alleles. We observed a genetic differentiation between NEA and SEA populations following a continuous pattern from north to south, and we show a significant and continuous decrease of HLA diversity by the same direction. This continuity is shaped by clinal distributions of many HLA lineages and alleles with increasing or decreasing frequencies along the latitude. These results bring new evidence in favor of the "overlapping model" proposed previously for East Asian peopling history, whereby modern humans migrated eastward from western Eurasia via two independent routes along each side of the Himalayas and, later, overlapped in East Asia across open land areas. Our study strongly suggests that intensive gene flow between NEA and SEA populations occurred and shaped the latitude-related continuous pattern of genetic variation and the peculiar HLA lineage and allele distributions observed in this continent. Probably for a very long period, the exact duration of these events remains to be estimated.

A genome-wide association study of atopic dermatitis identifies loci with overlapping effects on asthma and psoriasis

Atopic dermatitis (AD) is the most common dermatological disease of childhood. Many children with AD have asthma and AD shares regions of genetic linkage with psoriasis, another chronic inflammatory skin disease. We present here a genome-wide association study (GWAS) of childhood-onset AD in 1563 European cases with known asthma status and 4054 European controls. Using Illumina genotyping followed by imputation, we generated 268 034 consensus genotypes and in excess of 2 million single nucleotide polymorphisms (SNPs) for analysis. Association signals were assessed for replication in a second panel of 2286 European cases and 3160 European controls. Four loci achieved genome-wide significance for AD and replicated consistently across all cohorts. These included the epidermal differentiation complex (EDC) on chromosome 1, the genomic region proximal to LRRC32 on chromosome 11, the RAD50/IL13 locus on chromosome 5 and the major histocompatibility complex (MHC) on chromosome 6; reflecting action of classical HLA alleles. We observed variation in the contribution towards co-morbid asthma for these regions of association. We further explored the genetic relationship between AD, asthma and psoriasis by examining previously identified susceptibility SNPs for these diseases. We found considerable overlap between AD and psoriasis together with variable coincidence between allergic rhinitis (AR) and asthma. Our results indicate that the pathogenesis of AD incorporates immune and epidermal barrier defects with combinations of specific and overlapping effects at individual loci.

Functional diversity of the glutathione peroxidase gene family among human populations: implications for genetic predisposition to disease and drug response

Aim: To analyze the human genetic variation of glutathione peroxidases (GPX), estimating the functional differences among human populations and suggesting interethnic differences in predisposition to disease and drug response. Materials & methods: Using 1000 Genomes Project data, we analyzed 723 GPX variants in 1092 individuals belonging to 14 populations. Combining functional prediction analyses of coding and noncoding variants, we developed a method to estimate haplotype functionality. Results: GPX rare variants have a higher functional impact than common variants. The frequency among Asian patients of haplotypes associated with normal functionality is significantly higher for GPX1 and lower for GPX3 than for non-Asian patients; no adaptation signals in GPX1 and GPX3 were found in Asian patients. Conclusion: GPX1 and GPX3 differences may be associated with alterations in antioxidant capacity and redox regulation, which suggests diverse susceptibility to complex disease and diverse response to relevant drugs in Asians compared with individuals with other ethnic origins.

Original submitted 7 February 2013; Revision submitted 16 May 2013

Detecting and characterizing genomic signatures of positive selection in global populations

Natural selection is a significant force that shapes the architecture of the human genome and introduces diversity across global populations. The question of whether advantageous mutations have arisen in the human genome as a result of single or multiple mutation events remains unanswered except for the fact that there exist a handful of genes such as those that confer lactase persistence, affect skin pigmentation, or cause sickle cell anemia. We have developed a long-range-haplotype method for identifying genomic signatures of positive selection to complement existing methods, such as the integrated haplotype score (iHS) or cross-population extended haplotype homozygosity (XP-EHH), for locating signals across the entire allele frequency spectrum. Our method also locates the founder haplotypes that carry the advantageous variants and infers their corresponding population frequencies. This presents an opportunity to systematically interrogate the whole human genome whether a selection signal shared across different populations is the consequence of a single mutation process followed subsequently by gene flow between populations or of convergent evolution due to the occurrence of multiple independent mutation events either at the same variant or within the same gene. The application of our method to data from 14 populations across the world revealed that positive-selection events tend to cluster in populations of the same ancestry. Comparing the founder haplotypes for events that are present across different populations revealed that convergent evolution is a rare occurrence and that the majority of shared signals stem from the same evolutionary event.

A panel of ancestry informative markers to estimate and correct potential effects of population stratification in Han Chinese

Population stratification acts as a confounding factor in genetic association studies and may lead to false-positive or false-negative results. Previous studies have analyzed the genetic substructures in Han Chinese population, the largest ethnic group in the world comprising ∼20% of the global human population. In this study, we examined 5540 Han Chinese individuals with about 1 million single-nucleotide polymorphisms (SNPs) and screened a panel of ancestry informative markers (AIMs) to facilitate the discerning and controlling of population structure in future association studies on Han Chinese. Based on genome-wide data, we first confirmed our previous observation of the north-south differentiation in Han Chinese population. Second, we developed a panel of 150 validated SNP AIMs to determine the northern or southern origin of each Han Chinese individual. We further evaluated the performance of our AIMs panel in association studies in simulation analysis. Our results showed that this AIMs panel had sufficient power to discern and control population stratification in Han Chinese, which could significantly reduce false-positive rates in both genome-wide association studies (GWAS) and candidate gene association studies (CGAS). We suggest this AIMs panel be genotyped and used to control and correct population stratification in the study design or data analysis of future association studies, especially in CGAS which is the most popular approach to validate previous reports on genetic associations of diseases in post-GWAS era.

Efficiency of trans-ethnic genome-wide meta-analysis and fine-mapping

Genome-wide association studies have seen unprecedented success in identifying genetic loci that correlate with disease susceptibility and severity. Early phases of these studies have predominantly been performed in the Caucasian populations. The next phase in medical genetics is to extend the exploration across genetically diverse populations to leverage on larger sample sizes for locating smaller effects that may be present in most human populations. However, discoveries from these studies do not actually reveal the underlying functional changes to the human genome, but only point to broad regions stipulated by the extent of linkage disequilibrium (LD). Fine-mapping the functional variants can, however, be hampered by extensive LD, which can yield multiple perfect surrogates that are not distinguishable from the underlying causal variants, although several studies have illustrated the value of relying on multiple genetically diverse populations to narrow the candidate regions where the functional variants can be found in. Here, we explore the efficiency of trans-ethnic meta-analysis in discovering genetic association and in fine-mapping the causal variants by asking: are there any population diversity metrics that will be useful for: (i) identifying the populations or genomic regions where meta-analysis are likely to be more successful for discovering associations?; (ii) identifying the populations or loci to perform deep targeted sequencing for the purpose of fine-mapping causal variants? Our results indicate that simple metrics like the F(ST) or the population specificity of haplotypes are useful in trans-ethnic meta-analyses, while the degree of haplotype sharing and LD variation are informative of the efficiency in trans-ethnic fine-mapping.