varLD: a program for quantifying variation in linkage disequilibrium patterns between populations

Assessing the Risk Stratification of Breast Cancer Polygenic Risk Scores in a Brazilian Cohort

Polygenic risk scores (PRSs) for breast cancer have a clear clinical utility in risk prediction. PRS transferability across populations and ancestry groups is hampered by population-specific factors, ultimately leading to differences in variant effects, such as linkage disequilibrium and differences in variant frequency (allele frequency differences). Thus, locally sourced population-based phenotypic and genomic data sets are essential to assess the validity of PRSs derived from signals detected across populations. This study assesses the transferability of a breast cancer PRS composed of 313 risk variants (313-PRS) in a Brazilian trihybrid admixed ancestries (European, African, and Native American) whole-genome sequenced cohort, the Rare Genomes Project. 313-PRS was computed in the Rare Genomes Project (n = 853) using the UK Biobank (UKBB; n = 264,307) as reference. The Brazilian cohorts have a high European ancestry (EA) component, with allele frequency differences and to a lesser extent linkage disequilibrium patterns similar to those found in EA populations. The 313-PRS distribution was found to be inflated when compared with that of the UKBB, leading to potential overestimation of PRS-based risk if EA is taken as a standard. However, case controls lead to equivalent predictive power when compared with UKBB-EA samples with area under the receiver operating characteristic curve values of 0.66 to 0.62 compared with 0.63 for UKBB.

Uncovering genes underlying coat color variation in indigenous cattle breeds through genome-wide positive selection

The identification of candidate genes related to pigmentation and under selective sweep provides insights into the genetic basis of pigmentation and the evolutionary forces that have shaped this variation. The selective sweep events in the genes responsible for normal coat color in Indian cattle groups are still unknown. To find coat color genes displaying signs of selective sweeps in the indigenous cattle, we compiled a list of candidate genes previously investigated for their association with coat color and pigmentation. After that, we performed a genome-wide scan of positive selection signatures using the BovineSNP50K Bead Chip in 187 individuals of seven indigenous breeds. We applied a wide range of methods to find evidence of selection, such as Tajima's D, CLR, iHS, varLD, ROH, and FST. We found a total of sixteen genes under selective sweep, that were involved in coat color and pigmentation physiology. These genes are CRIM1 in Gir, MC1R in Sahiwal, MYO5A, PMEL and POMC in Tharparkar, TYRP1, ERBB2, and ASIP in Red Sindhi, MITF, LOC789175, PAX3 and TYR in Ongole, and IRF2, SDR165 and, KIT in Nelore, ADAMTS19 in Hariana. These genes are related to melanin synthesis, the biology of melanocytes and melanosomes, and the migration and survival of melanocytes during development.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

Fully exploiting SNP arrays: a systematic review on the tools to extract underlying genomic structure

Single nucleotide polymorphisms (SNPs) are the most abundant type of genomic variation and the most accessible to genotype in large cohorts. However, they individually explain a small proportion of phenotypic differences between individuals. Ancestry, collective SNP effects, structural variants, somatic mutations or even differences in historic recombination can potentially explain a high percentage of genomic divergence. These genetic differences can be infrequent or laborious to characterize; however, many of them leave distinctive marks on the SNPs across the genome allowing their study in large population samples. Consequently, several methods have been developed over the last decade to detect and analyze different genomic structures using SNP arrays, to complement genome-wide association studies and determine the contribution of these structures to explain the phenotypic differences between individuals. We present an up-to-date collection of available bioinformatics tools that can be used to extract relevant genomic information from SNP array data including population structure and ancestry; polygenic risk scores; identity-by-descent fragments; linkage disequilibrium; heritability and structural variants such as inversions, copy number variants, genetic mosaicisms and recombination histories. From a systematic review of recently published applications of the methods, we describe the main characteristics of R packages, command-line tools and desktop applications, both free and commercial, to help make the most of a large amount of publicly available SNP data.

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

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

Genetic Risk Prediction of COVID-19 Susceptibility and Severity in the Indian Population

Host genetic variants can determine their susceptibility to COVID-19 infection and severity as noted in a recent Genome-wide Association Study (GWAS). Given the prominent genetic differences in Indian sub-populations as well as differential prevalence of COVID-19, here, we compute genetic risk scores in diverse Indian sub-populations that may predict differences in the severity of COVID-19 outcomes. We utilized the top 100 most significantly associated single-nucleotide polymorphisms (SNPs) from a GWAS by Pairo-Castineira et al. determining the genetic susceptibility to severe COVID-19 infection, to compute population-wise polygenic risk scores (PRS) for populations represented in the Indian Genome Variation Consortium (IGVC) database. Using a generalized linear model accounting for confounding variables, we found that median PRS was significantly associated (p < 2 x 10⁻¹⁶) with COVID-19 mortality in each district corresponding to the population studied and had the largest effect on mortality (regression coefficient = 10.25). As a control we repeated our analysis on randomly selected 100 non-associated SNPs several times and did not find significant association. Therefore, we conclude that genetic susceptibility may play a major role in determining the differences in COVID-19 outcomes and mortality across the Indian sub-continent. We suggest that combining PRS with other observed risk-factors in a Bayesian framework may provide a better prediction model for ascertaining high COVID-19 risk groups and to design more effective public health resource allocation and vaccine distribution schemes.

Selection for environmental variance of litter size in rabbits involves genes in pathways controlling animal resilience

Background

Environmental variance (V_E) is partially under genetic control, which means that the V_E of individuals that share the same environment can differ because they have different genotypes. Previously, a divergent selection experiment for V_E of litter size (LS) during 13 generations in rabbit yielded a successful response and revealed differences in resilience between the divergent lines. The aim of the current study was to identify signatures of selection in these divergent lines to better understand the molecular mechanisms and pathways that control V_E of LS and animal resilience. Three methods (F_ST, ROH and varLD) were used to identify signatures of selection in a set of 473 genotypes from these rabbit lines (377) and a base population (96). A whole-genome sequencing (WGS) analysis was performed on 54 animals to detect genes with functional mutations.

Results

By combining signatures of selection and WGS data, we detected 373 genes with functional mutations in their transcription units, among which 111 had functions related to the immune system, stress response, reproduction and embryo development, and/or carbohydrate and lipid metabolism. The genes TTC23L, FBXL20, GHDC, ENSOCUG00000031631, SLC18A1, CD300LG, MC2R, and ENSOCUG00000006264 were particularly relevant, since each one carried a functional mutation that was fixed in one of the rabbit lines and absent in the other line. In the 3ʹUTR region of the MC2R and ENSOCUG00000006264 genes, we detected a novel insertion/deletion (INDEL) variant.

Conclusions

Our findings provide further evidence in favour of V_E as a measure of animal resilience. Signatures of selection were identified for V_E of LS in genes that have a functional mutation in their transcription units and are mostly implicated in the immune response and stress response pathways. However, the real implications of these genes for V_E and animal resilience will need to be assessed through functional analyses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00653-y.

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

Selection in Australian Thoroughbred horses acts on a locus associated with early two-year old speed

Thoroughbred horse racing is a global sport with major hubs in Europe, North America, Australasia and Japan. Regional preferences for certain traits have resulted in phenotypic variation that may result from adaptation to the local racing ecosystem. Here, we test the hypothesis that genes selected for regional phenotypic variation may be identified by analysis of selection signatures in pan-genomic SNP genotype data. Comparing Australian to non-Australian Thoroughbred horses (n = 99), the most highly differentiated loci in a composite selection signals (CSS) analysis were on ECA6 (34.75–34.85 Mb), ECA14 (33.2–33.52 Mb and 35.52–36.94 Mb) and ECA16 (24.28–26.52 Mb) in regions containing candidate genes for exercise adaptations including cardiac function (ARHGAP26, HBEGF, SRA1), synapse development and locomotion (APBB3, ATXN7, CLSTN3), stress response (NR3C1) and the skeletal muscle response to exercise (ARHGAP26, NDUFA2). In a genome-wide association study for field-measured speed in two-year-olds (n = 179) SNPs contained within the single association peak (33.2–35.6 Mb) overlapped with the ECA14 CSS signals and spanned a protocadherin gene cluster. Association tests using higher density SNP genotypes across the ECA14 locus identified a SNP within the PCDHGC5 gene associated with elite racing performance (n = 922). These results indicate that there may be differential selection for racing performance under racing and management conditions that are specific to certain geographic racing regions. In Australia breeders have principally selected horses for favourable genetic variants at loci containing genes that modulate behaviour, locomotion and skeletal muscle physiology that together appear to be contributing to early two-year-old speed.

Identification of deleterious and regulatory genomic variations in known asthma loci

Background

Candidate gene and genome-wide association studies have identified hundreds of asthma risk loci. The majority of associated variants, however, are not known to have any biological function and are believed to represent markers rather than true causative mutations. We hypothesized that many of these associated markers are in linkage disequilibrium (LD) with the elusive causative variants.

Methods

We compiled a comprehensive list of 449 asthma-associated variants previously reported in candidate gene and genome-wide association studies. Next, we identified all sequence variants located within the 305 unique genes using whole-genome sequencing data from the 1000 Genomes Project. Then, we calculated the LD between known asthma variants and the sequence variants within each gene. LD variants identified were then annotated to determine those that are potentially deleterious and/or functional (i.e. coding or regulatory effects on the encoded transcript or protein).

Results

We identified 10,130 variants in LD (r² > 0.6) with known asthma variants. Annotations of these LD variants revealed that several have potentially deleterious effects including frameshift, alternate splice site, stop-lost, and missense. Moreover, 24 of the LD variants have been reported to regulate gene expression as expression quantitative trait loci (eQTLs).

Conclusions

This study is proof of concept that many of the genetic loci previously associated with complex diseases such as asthma are not causative but represent markers of disease, which are in LD with the elusive causative variants. We hereby report a number of potentially deleterious and regulatory variants that are in LD with the reported asthma loci. These reported LD variants could account for the original association signals with asthma and represent the true causative mutations at these loci.

Electronic supplementary material

The online version of this article (10.1186/s12931-018-0953-2) contains supplementary material, which is available to authorized users.

A genome-wide scan for signatures of selection in Azeri and Khuzestani buffalo breeds

Background

Identification of genomic regions that have been targets of selection may shed light on the genetic history of livestock populations and help to identify variation controlling commercially important phenotypes. The Azeri and Kuzestani buffalos are the most common indigenous Iranian breeds which have been subjected to divergent selection and are well adapted to completely different regions. Examining the genetic structure of these populations may identify genomic regions associated with adaptation to the different environments and production goals.

Results

A set of 385 water buffalo samples from Azeri (N = 262) and Khuzestani (N = 123) breeds were genotyped using the Axiom® Buffalo Genotyping 90 K Array. The unbiased fixation index method (F_ST) was used to detect signatures of selection. In total, 13 regions with outlier F_ST values (0.1%) were identified. Annotation of these regions using the UMD3.1 Bos taurus Genome Assembly was performed to find putative candidate genes and QTLs within the selected regions. Putative candidate genes identified include FBXO9, NDFIP1, ACTR3, ARHGAP26, SERPINF2, BOLA-DRB3, BOLA-DQB, CLN8, and MYOM2.

Conclusions

Candidate genes identified in regions potentially under selection were associated with physiological pathways including milk production, cytoskeleton organization, growth, metabolic function, apoptosis and domestication-related changes include immune and nervous system development. The QTL identified are involved in economically important traits in buffalo related to milk composition, udder structure, somatic cell count, meat quality, and carcass and body weight.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4759-x) contains supplementary material, which is available to authorized users.

Risk of nontyphoidal Salmonella bacteraemia in African children is modified by STAT4

Nontyphoidal Salmonella (NTS) is a major cause of bacteraemia in Africa. The disease typically affects HIV-infected individuals and young children, causing substantial morbidity and mortality. Here we present a genome-wide association study (180 cases, 2677 controls) and replication analysis of NTS bacteraemia in Kenyan and Malawian children. We identify a locus in STAT4, rs13390936, associated with NTS bacteraemia. rs13390936 is a context-specific expression quantitative trait locus for STAT4 RNA expression, and individuals carrying the NTS-risk genotype demonstrate decreased interferon-γ (IFNγ) production in stimulated natural killer cells, and decreased circulating IFNγ concentrations during acute NTS bacteraemia. The NTS-risk allele at rs13390936 is associated with protection against a range of autoimmune diseases. These data implicate interleukin-12-dependent IFNγ-mediated immunity as a determinant of invasive NTS disease in African children, and highlight the shared genetic architecture of infectious and autoimmune disease.

Regional evaluation of childhood acute lymphoblastic leukemia genetic susceptibility loci among Japanese

Genome-wide association studies (GWAS) performed mostly in populations of European and Hispanic ancestry have confirmed an inherited genetic basis for childhood acute lymphoblastic leukemia (ALL), but these associations are less clear in other races/ethnicities. DNA samples from ALL patients (aged 0–19 years) previously enrolled onto a Tokyo Children’s Cancer Study Group trial were collected during 2013–2015, and underwent single nucleotide polymorphism (SNP) microarray genotyping resulting in 527 B-cell ALL for analysis. Cases and control data for 3,882 samples from the Nagahama Study Group and Aichi Cancer Center Study were combined, and association analyses across 10 previous GWAS-identified regions were performed after targeted SNP imputation. Linkage disequilibrium (LD) patterns in Japanese and other populations were evaluated using the varLD score based on 1000 Genomes data. Risk associations for ARID5B (rs10821936, OR = 1.84, P = 6 × 10⁻¹⁷) and PIP4K2A (rs7088318, OR = 0.76, P = 2 × 10⁻⁴) directly transferred to Japanese, and the IKZF1 association was detected by an alternate SNP (rs1451367, OR = 1.52, P = 2 × 10⁻⁶). Marked regional LD differences between Japanese and Europeans was observed for most of the remaining loci for which associations did not transfer, including CEBPE, CDKN2A, CDKN2B, and ELK3. This study represents a first step towards characterizing the role of genetic susceptibility in childhood ALL risk in Japanese.

Testing Domestication Scenarios of Lima Bean (Phaseolus lunatus L.) in Mesoamerica: Insights from Genome-Wide Genetic Markers

Plant domestication can be seen as a long-term process that involves a complex interplay among demographic processes and evolutionary forces. Previous studies have suggested two domestication scenarios for Lima bean in Mesoamerica: two separate domestication events, one from gene pool MI in central-western Mexico and another one from gene pool MII in the area Guatemala-Costa Rica, or a single domestication from gene pool MI in central-western Mexico followed by post-domestication gene flow with wild populations. In this study we evaluated the genetic structure of the wild gene pool and tested these two competing domestication scenarios of Lima bean in Mesoamerica by applying an ABC approach to a set of genome-wide SNP markers. The results confirm the existence of three gene pools in wild Lima bean, two Mesoamerican gene pools (MI and MII) and the Andean gene pool (AI), and suggest the existence of another gene pool in central Colombia. The results indicate that although both domestication scenarios may be supported by genetic data, higher statistical support was given to the single domestication scenario in central-western Mexico followed by admixture with wild populations. Domestication would have involved strong founder effects reflected in loss of genetic diversity and increased LD levels in landraces. Genomic regions affected by selection were detected and these may harbor candidate genes related to domestication.

A meta-analysis of genome-wide association studies of asthma in Puerto Ricans

Puerto Ricans are disproportionately affected with asthma in the USA. In this study, we aim to identify genetic variants that confer susceptibility to asthma in Puerto Ricans.We conducted a meta-analysis of genome-wide association studies (GWAS) of asthma in Puerto Ricans, including participants from: the Genetics of Asthma in Latino Americans (GALA) I-II, the Hartford-Puerto Rico Study and the Hispanic Community Health Study. Moreover, we examined whether susceptibility loci identified in previous meta-analyses of GWAS are associated with asthma in Puerto Ricans.The only locus to achieve genome-wide significance was chromosome 17q21, as evidenced by our top single nucleotide polymorphism (SNP), rs907092 (OR 0.71, p=1.2×10^-12) at IKZF3 Similar to results in non-Puerto Ricans, SNPs in genes in the same linkage disequilibrium block as IKZF3 (e.g. ZPBP2, ORMDL3 and GSDMB) were significantly associated with asthma in Puerto Ricans. With regard to results from a meta-analysis in Europeans, we replicated findings for rs2305480 at GSDMB, but not for SNPs in any other genes. On the other hand, we replicated results from a meta-analysis of North American populations for SNPs at IL1RL1, TSLP and GSDMB but not for IL33Our findings suggest that common variants on chromosome 17q21 have the greatest effects on asthma in Puerto Ricans.

A survey of sub-Saharan gene flow into the Mediterranean at risk loci for coronary artery disease

This study tries to find detectable signals of gene flow of Sub-Saharan origin into the Mediterranean in four genomic regions previously associated with coronary artery disease. A total of 366 single-nucleotide polymorphisms were genotyped in 772 individuals from 10 Mediterranean countries. Population structure analyses were performed, in which a noticeable Sub-Saharan component was found in the studied samples. The overall percentage of this Sub-Saharan component presents differences between the two Mediterranean coasts. D-statistics suggest possible Sub-Saharan introgression into one of the studied genomic regions (10q11). We also found differences in linkage disequilibrium patterns between the two Mediterranean coasts, possibly attributable to differential Sub-Saharan admixture. Our results confirm the potentially important role of human demographic history when performing epidemiological studies.

On the performance of tests for the detection of signatures of selection: a case study with the Spanish autochthonous beef cattle populations

Background

Procedures for the detection of signatures of selection can be classified according to the source of information they use to reject the null hypothesis of absence of selection. Three main groups of tests can be identified that are based on: (1) the analysis of the site frequency spectrum, (2) the study of the extension of the linkage disequilibrium across the length of the haplotypes that surround the polymorphism, and (3) the differentiation among populations. The aim of this study was to compare the performance of a subset of these procedures by using a dataset on seven Spanish autochthonous beef cattle populations.

Results

Analysis of the correlations between the logarithms of the statistics that were obtained by 11 tests for detecting signatures of selection at each single nucleotide polymorphism confirmed that they can be clustered into the three main groups mentioned above. A factor analysis summarized the results of the 11 tests into three canonical axes that were each associated with one of the three groups. Moreover, the signatures of selection identified with the first and second groups of tests were shared across populations, whereas those with the third group were more breed-specific. Nevertheless, an enrichment analysis identified the metabolic pathways that were associated with each group; they coincided with canonical axes and were related to immune response, muscle development, protein biosynthesis, skin and pigmentation, glucose metabolism, fat metabolism, embryogenesis and morphology, heart and uterine metabolism, regulation of the hypothalamic–pituitary–thyroid axis, hormonal, cellular cycle, cell signaling and extracellular receptors.

Conclusions

We show that the results of the procedures used to identify signals of selection differed substantially between the three groups of tests. However, they can be classified using a factor analysis. Moreover, each canonical factor that coincided with a group of tests identified different signals of selection, which could be attributed to processes of selection that occurred at different evolutionary times. Nevertheless, the metabolic pathways that were associated with each group of tests were similar, which suggests that the selection events that occurred during the evolutionary history of the populations probably affected the same group of traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0258-1) contains supplementary material, which is available to authorized users.

Trans-ethnic Meta-analysis and Functional Annotation Illuminates the Genetic Architecture of Fasting Glucose and Insulin

Knowledge of the genetic basis of the type 2 diabetes (T2D)-related quantitative traits fasting glucose (FG) and insulin (FI) in African ancestry (AA) individuals has been limited. In non-diabetic subjects of AA (n = 20,209) and European ancestry (EA; n = 57,292), we performed trans-ethnic (AA+EA) fine-mapping of 54 established EA FG or FI loci with detailed functional annotation, assessed their relevance in AA individuals, and sought previously undescribed loci through trans-ethnic (AA+EA) meta-analysis. We narrowed credible sets of variants driving association signals for 22/54 EA-associated loci; 18/22 credible sets overlapped with active islet-specific enhancers or transcription factor (TF) binding sites, and 21/22 contained at least one TF motif. Of the 54 EA-associated loci, 23 were shared between EA and AA. Replication with an additional 10,096 AA individuals identified two previously undescribed FI loci, chrX FAM133A (rs213676) and chr5 PELO (rs6450057). Trans-ethnic analyses with regulatory annotation illuminate the genetic architecture of glycemic traits and suggest gene regulation as a target to advance precision medicine for T2D. Our approach to utilize state-of-the-art functional annotation and implement trans-ethnic association analysis for discovery and fine-mapping offers a framework for further follow-up and characterization of GWAS signals of complex trait loci.

Evaluation of type 2 diabetes genetic risk variants in Chinese adults: findings from 93,000 individuals from the China Kadoorie Biobank

AIMS/HYPOTHESIS

Genome-wide association studies (GWAS) have discovered many risk variants for type 2 diabetes. However, estimates of the contributions of risk variants to type 2 diabetes predisposition are often based on highly selected case-control samples, and reliable estimates of population-level effect sizes are missing, especially in non-European populations.

METHODS

The individual and cumulative effects of 59 established type 2 diabetes risk loci were measured in a population-based China Kadoorie Biobank (CKB) study of 93,000 Chinese adults, including >7,100 diabetes cases.

RESULTS

Association signals were directionally consistent between CKB and the original discovery GWAS: of 56 variants passing quality control, 48 showed the same direction of effect (binomial test, p = 2.3 × 10(-8)). We observed a consistent overall trend towards lower risk variant effect sizes in CKB than in case-control samples of GWAS meta-analyses (mean 19-22% decrease in log odds, p ≤ 0.0048), likely to reflect correction of both 'winner's curse' and spectrum bias effects. The association with risk of diabetes of a genetic risk score, based on lead variants at 25 loci considered to act through beta cell function, demonstrated significant interactions with several measures of adiposity (BMI, waist circumference [WC], WHR and percentage body fat [PBF]; all p interaction < 1 × 10(-4)), with a greater effect being observed in leaner adults.

CONCLUSIONS/INTERPRETATION

Our study provides further evidence of shared genetic architecture for type 2 diabetes between Europeans and East Asians. It also indicates that even very large GWAS meta-analyses may be vulnerable to substantial inflation of effect size estimates, compared with those observed in large-scale population-based cohort studies.

ACCESS TO RESEARCH MATERIALS

Details of how to access China Kadoorie Biobank data and details of the data release schedule are available from www.ckbiobank.org/site/Data+Access .

A Meta-Assembly of Selection Signatures in Cattle

Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.

Pharmacogenetics of drug-metabolizing enzymes in US Hispanics

Although the Hispanic population is continuously growing in the United States, they are underrepresented in pharmacogenetic studies. This review addresses the need for compiling available pharmacogenetic data in US Hispanics, discussing the prevalence of clinically relevant polymorphisms in pharmacogenes encoding for drug-metabolizing enzymes. CYP3A5*3 (0.245-0.867) showed the largest frequency in a US Hispanic population. A higher prevalence of CYP2C9*3, CYP2C19*4, and UGT2B7 IVS1+985 A>G was observed in US Hispanic vs. non-Hispanic populations. We found interethnic and intraethnic variability in frequencies of genetic polymorphisms for metabolizing enzymes, which highlights the need to define the ancestries of participants in pharmacogenetic studies. New approaches should be integrated in experimental designs to gain knowledge about the clinical relevance of the unique combination of genetic variants occurring in this admixed population. Ethnic subgroups in the US Hispanic population may harbor variants that might be part of multiple causative loci or in linkage-disequilibrium with functional variants. Pharmacogenetic studies in Hispanics should not be limited to ascertain commonly studied polymorphisms that were originally identified in their parental populations. The success of the Personalized Medicine paradigm will depend on recognizing genetic diversity between and within US Hispanics and the uniqueness of their genetic backgrounds.

Analysis of Genomic Regions Associated With Coronary Artery Disease Reveals Continent-Specific Single Nucleotide Polymorphisms in North African Populations

Background

In recent years, several genomic regions have been robustly associated with coronary artery disease (CAD) in different genome-wide association studies (GWASs) conducted mainly in people of European descent. These kinds of data are lacking in African populations, even though heart diseases are a major cause of premature death and disability.

Methods

Here, 384 single nucleotide polymorphisms (SNPs) in the top four CAD risk regions (1p13, 1q41, 9p21, and 10q11) were genotyped in 274 case-control samples from Morocco and Tunisia, with the aim of analyzing for the first time if the associations found in European populations were transferable to North Africans.

Results

The results indicate that, as in Europe, these four genetic regions are also important for CAD risk in North Africa. However, the individual SNPs associated with CAD in Africa are different from those identified in Europe in most cases (1p13, 1q41, and 9p21). Moreover, the seven risk variants identified in North Africans are efficient in discriminating between cases and controls in North African populations, but not in European populations.

Conclusions

This study indicates a disparity in markers associated to CAD susceptibility between North Africans and Europeans that may be related to population differences in the chromosomal architecture of these risk regions.

Trans-ancestry genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation

We carried out a trans-ancestry genome-wide association and replication study of blood pressure phenotypes among up to 320,251 individuals of East Asian, European and South Asian ancestry. We find genetic variants at 12 new loci to be associated with blood pressure (P = 3.9 × 10(-11) to 5.0 × 10(-21)). The sentinel blood pressure SNPs are enriched for association with DNA methylation at multiple nearby CpG sites, suggesting that, at some of the loci identified, DNA methylation may lie on the regulatory pathway linking sequence variation to blood pressure. The sentinel SNPs at the 12 new loci point to genes involved in vascular smooth muscle (IGFBP3, KCNK3, PDE3A and PRDM6) and renal (ARHGAP24, OSR1, SLC22A7 and TBX2) function. The new and known genetic variants predict increased left ventricular mass, circulating levels of NT-proBNP, and cardiovascular and all-cause mortality (P = 0.04 to 8.6 × 10(-6)). Our results provide new evidence for the role of DNA methylation in blood pressure regulation.

Potential Signals of Natural Selection in the Top Risk Loci for Coronary Artery Disease: 9p21 and 10q11

BACKGROUND

Coronary artery disease (CAD) is a complex disease and the leading cause of death in the world. Populations of different ancestry do not always share the same risk markers. Natural selective processes may be the cause of some of the population differences detected for specific risk mutations.

OBJECTIVE

In this study, 384 single nucleotide polymorphisms (SNPs) located in four genomic regions associated with CAD (1p13, 1q41, 9p21 and 10q11) are analysed in a set of 19 populations from Europe, Middle East and North Africa and also in Asian and African samples from the 1000 Genomes Project. The aim of this survey is to explore for the first time whether the genetic variability in these genomic regions is better explained by demography or by natural selection.

RESULTS

The results indicate significant differences in the structure of genetic variation and in the LD patterns among populations that probably explain the population disparities found in markers of susceptibility to CAD.

CONCLUSIONS

The results are consistent with potential signature of positive selection in the 9p21 region and of balancing selection in the 9p21 and 10q11. Specifically, in Europe three CAD risk markers in the 9p21 region (rs9632884, rs1537371 and rs1333042) show consistent signals of positive selection. The results of this study are consistent with a potential selective role of CAD in the configuration of genetic diversity in current human populations.

Detection of selection signatures in Piemontese and Marchigiana cattle, two breeds with similar production aptitudes but different selection histories

BACKGROUND

Domestication and selection are processes that alter the pattern of within- and between-population genetic variability. They can be investigated at the genomic level by tracing the so-called selection signatures. Recently, sequence polymorphisms at the genome-wide level have been investigated in a wide range of animals. A common approach to detect selection signatures is to compare breeds that have been selected for different breeding goals (i.e. dairy and beef cattle). However, genetic variations in different breeds with similar production aptitudes and similar phenotypes can be related to differences in their selection history.

METHODS

In this study, we investigated selection signatures between two Italian beef cattle breeds, Piemontese and Marchigiana, using genotyping data that was obtained with the Illumina BovineSNP50 BeadChip. The comparison was based on the fixation index (Fst), combined with a locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach. In addition, analyses of Fst were carried out to confirm candidate genes. In particular, data were processed using the varLD method, which compares the regional variation of linkage disequilibrium between populations.

RESULTS

Genome scans confirmed the presence of selective sweeps in the genomic regions that harbour candidate genes that are known to affect productive traits in cattle such as DGAT1, ABCG2, CAPN3, MSTN and FTO. In addition, several new putative candidate genes (for example ALAS1, ABCB8, ACADS and SOD1) were detected.

CONCLUSIONS

This study provided evidence on the different selection histories of two cattle breeds and the usefulness of genomic scans to detect selective sweeps even in cattle breeds that are bred for similar production aptitudes.

F-box/LRR-repeat protein 7 is genetically associated with Alzheimer's disease

Objective

In the context of late-onset Alzheimer’s disease (LOAD) over 20 genes have been identified but, aside APOE, all show small effect sizes, leaving a large part of the genetic component unexplained. Admixed populations, such as Caribbean Hispanics, can provide a valuable contribution because of their unique genetic profile and higher incidence of the disease. We aimed to identify novel loci associated with LOAD.

Methods

About 4514 unrelated Caribbean Hispanics (2451 cases and 2063 controls) were selected for genome-wide association analysis. Significant loci were further tested in the expanded cohort that also included related family members (n = 5300). Two AD-like transgenic mice models (J20 and rTg4510) were used to study gene expression. Independent data sets of non-Hispanic Whites and African Americans were used to further validate findings, along with publicly available brain expression data sets.

Results

A novel locus, rs75002042 in FBXL7 (5p15.1), was found genome-wide significant in the case–control cohort (odd ratio [OR] = 0.61, P = 6.19E-09) and confirmed in the related members cohorts (OR = 0.63, P = 4.7E-08). Fbxl7 protein was overexpressed in both AD-like transgenic mice compared to wild-type littermates. Publicly available microarray studies also showed significant overexpression of Fbxl7 in LOAD brains compared to nondemented controls. single-nucleotide polymorphism (SNP) rs75002042 was in complete linkage disequilibrium with other variants in two independent non-Hispanic White and African American data sets (0.0005 < P < 0.02) used for replication.

Interpretation

FBXL7, encodes a subcellular protein involved in phosphorylation-dependent ubiquitination processes and displays proapoptotic activity. F-box proteins also modulate inflammation and innate immunity, which may be important in LOAD pathogenesis. Further investigations are needed to validate and understand its role in this and other populations.

Trans-ethnic genome-wide association studies: advantages and challenges of mapping in diverse populations

Genome-wide association studies (GWASs) are the method most often used by geneticists to interrogate the human genome, and they provide a cost-effective way to identify the genetic variants underpinning complex traits and diseases. Most initial GWASs have focused on genetically homogeneous cohorts from European populations given the limited availability of ethnic minority samples and so as to limit population stratification effects. Transethnic studies have been invaluable in explaining the heritability of common quantitative traits, such as height, and in examining the genetic architecture of complex diseases, such as type 2 diabetes. They provide an opportunity for large-scale signal replication in independent populations and for cross-population meta-analyses to boost statistical power. In addition, transethnic GWASs enable prioritization of candidate genes, fine-mapping of functional variants, and potentially identification of SNPs associated with disease risk in admixed populations, by taking advantage of natural differences in genomic linkage disequilibrium across ethnically diverse populations. Recent efforts to assess the biological function of variants identified by GWAS have highlighted the need for large-scale replication, meta-analyses and fine-mapping across worldwide populations of ethnically diverse genetic ancestries. Here, we review recent advances and new approaches that are important to consider when performing, designing or interpreting transethnic GWASs, and we highlight existing challenges, such as the limited ability to handle heterogeneity in linkage disequilibrium across populations and limitations in dissecting complex architectures, such as those found in recently admixed populations.

A gene-based analysis of variants in the serum/glucocorticoid regulated kinase (SGK) genes with blood pressure responses to sodium intake: the GenSalt Study

BACKGROUND

Serum and glucocorticoid regulated kinase (SGK) plays a critical role in the regulation of renal sodium transport. We examined the association between SGK genes and salt sensitivity of blood pressure (BP) using single-marker and gene-based association analysis.

METHODS

A 7-day low-sodium (51.3 mmol sodium/day) followed by a 7-day high-sodium intervention (307.8 mmol sodium/day) was conducted among 1,906 Chinese participants. BP measurements were obtained at baseline and each intervention using a random-zero sphygmomanometer. Additive associations between each SNP and salt-sensitivity phenotypes were assessed using a mixed linear regression model to account for family dependencies. Gene-based analyses were conducted using the truncated p-value method. The Bonferroni-method was used to adjust for multiple testing in all analyses.

RESULTS

In single-marker association analyses, SGK1 marker rs2758151 was significantly associated with diastolic BP (DBP) response to high-sodium intervention (P = 0.0010). DBP responses (95% confidence interval) to high-sodium intervention for genotypes C/C, C/T, and T/T were 2.04 (1.57 to 2.52), 1.79 (1.42 to 2.16), and 0.85 (0.30 to 1.41) mmHg, respectively. Similar trends were observed for SBP and MAP responses although not significant (P = 0.15 and 0.0026, respectively). In addition, gene-based analyses demonstrated significant associations between SGK1 and SBP, DBP and MAP responses to high sodium intervention (P = 0.0002, 0.0076, and 0.00001, respectively). Neither SGK2 nor SGK3 were associated with the salt-sensitivity phenotypes in single-maker or gene-based analyses.

CONCLUSIONS

The current study identified association of the SGK1 gene and BP salt-sensitivity in the Han Chinese population. Further studies are warranted to identify causal SGK1 gene variants.

Fast principal component analysis of large-scale genome-wide data

Principal component analysis (PCA) is routinely used to analyze genome-wide single-nucleotide polymorphism (SNP) data, for detecting population structure and potential outliers. However, the size of SNP datasets has increased immensely in recent years and PCA of large datasets has become a time consuming task. We have developed flashpca, a highly efficient PCA implementation based on randomized algorithms, which delivers identical accuracy in extracting the top principal components compared with existing tools, in substantially less time. We demonstrate the utility of flashpca on both HapMap3 and on a large Immunochip dataset. For the latter, flashpca performed PCA of 15,000 individuals up to 125 times faster than existing tools, with identical results, and PCA of 150,000 individuals using flashpca completed in 4 hours. The increasing size of SNP datasets will make tools such as flashpca essential as traditional approaches will not adequately scale. This approach will also help to scale other applications that leverage PCA or eigen-decomposition to substantially larger datasets.

Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle

BACKGROUND

Signatures of selection are regions in the genome that have been preferentially increased in frequency and fixed in a population because of their functional importance in specific processes. These regions can be detected because of their lower genetic variability and specific regional linkage disequilibrium (LD) patterns.

METHODS

By comparing the differences in regional LD variation between dairy and beef cattle types, and between indicine and taurine subspecies, we aim at finding signatures of selection for production and adaptation in cattle breeds. The VarLD method was applied to compare the LD variation in the autosomal genome between breeds, including Angus and Brown Swiss, representing taurine breeds, and Nelore and Gir, representing indicine breeds. Genomic regions containing the top 0.01 and 0.1 percentile of signals were characterized using the UMD3.1 Bos taurus genome assembly to identify genes in those regions and compared with previously reported selection signatures and regions with copy number variation.

RESULTS

For all comparisons, the top 0.01 and 0.1 percentile included 26 and 165 signals and 17 and 125 genes, respectively, including TECRL, BT.23182 or FPPS, CAST, MYOM1, UVRAG and DNAJA1.

CONCLUSIONS

The VarLD method is a powerful tool to identify differences in linkage disequilibrium between cattle populations and putative signatures of selection with potential adaptive and productive importance.

Accounting for linkage disequilibrium in association analysis of diverse populations

The National Human Genome Research Institute's catalog of published genome-wide association studies (GWAS) lists over 10,000 genetic variants collectively associated with over 800 human diseases or traits. Most of these GWAS have been conducted in European-ancestry populations. Findings gleaned from these studies have led to identification of disease-associated loci and biologic pathways involved in disease etiology. In multiple instances, these genomic findings have led to the development of novel medical therapies or evidence for prescribing a given drug as the appropriate treatment for a given individual beyond phenotypic appearances or socially defined constructs of race or ethnicity. Such findings have implications for populations throughout the globe and GWAS are increasingly being conducted in more diverse populations. A major challenge for investigators seeking to follow up genomic findings between diverse populations is discordant patterns of linkage disequilibrium (LD). We provide an overview of common measures of LD and opportunities for their use in novel methods designed to address challenges associated with following up GWAS conducted in European-ancestry populations in African-ancestry populations or, more generally, between populations with discordant LD patterns. We detail the strengths and weaknesses associated with different approaches. We also describe application of these strategies in follow-up studies of populations with concordant LD patterns (replication) or discordant LD patterns (transferability) as well as fine-mapping studies. We review application of these methods to a variety of traits and diseases.

Genome-wide association study meta-analysis reveals transethnic replication of mean arterial and pulse pressure loci

We conducted a genome-wide association study meta-analysis of mean arterial pressure and pulse pressure among 26,600 East Asian participants (stage 1) followed by replication study of up to 28,783 participants (stage 2). For novel loci, statistical significance was determined by a P<5.0×10(-8) in joint analysis of stage 1 and stage 2 data. For loci reported by the previous mean arterial and pulse pressure genome-wide association study meta-analysis in Europeans, evidence of transethnic replication was determined by consistency in effect direction and a Bonferroni-corrected P<1.4×10(-3). No novel loci were identified by the current study. Five independent mean arterial pressure variants demonstrated robust evidence for transethnic replication including rs17249754 at ATP2B1 (P=7.5×10(-15)), rs2681492 at ATP2B1 (P=3.4×10(-7)), rs11191593 at NT5C2 (1.1×10(-6)), rs3824755 at CYP17A1 (P=1.2×10(-6)), and rs13149993 at FGF5 (P=2.4×10(-4)). Two additional variants showed suggestive evidence of transethnic replication (consistency in effect direction and P<0.05), including rs319690 at MAP4 (P=0.014) and rs1173771 at NPR3 (P=0.018). For pulse pressure, robust evidence of replication was identified for 2 independent variants, including rs17249754 at ATP2B1 (P=1.2×10(-5)) and rs11191593 at NT5C2 (P=1.1×10(-3)), with suggestive evidence of replication among an additional 2 variants including rs3824755 at CYP17A1 (P=6.1×10(-3)) and rs2681492 at ATP2B1 (P=9.0×10(-3)). Replicated variants demonstrated consistency in effect sizes between East Asian and European samples, with effect size differences ranging from 0.03 to 0.24 mm Hg for mean arterial pressure and from 0.03 to 0.21 mm Hg for pulse pressure. In conclusion, we present the first evidence of transethnic replication of several mean arterial and pulse pressure loci in an East Asian population.

Are C-reactive protein associated genetic variants associated with serum levels and retinal markers of microvascular pathology in Asian populations from Singapore?

Introduction

C-reactive protein (CRP) levels are associated with cardiovascular disease and systemic inflammation. We assessed whether CRP-associated loci were associated with serum CRP and retinal markers of microvascular disease, in Asian populations.

Methods

Genome-wide association analysis (GWAS) for serum CRP was performed in East-Asian Chinese (N = 2,434) and Malays (N = 2,542) and South-Asian Indians (N = 2,538) from Singapore. Leveraging on GWAS data, we assessed, in silico, association levels among the Singaporean datasets for 22 recently identified CRP-associated loci. At loci where directional inconsistencies were observed, quantification of inter-ethnic linkage disequilibrium (LD) difference was determined. Next, we assessed association for a variant at CRP and retinal vessel traits [central retinal artery equivalent (CRAE) and central retinal vein equivalent (CRVE)] in a total of 24,132 subjects of East-Asian, South-Asian and European ancestry.

Results

Serum CRP was associated with SNPs in/near APOE, CRP, HNF1A and LEPR (p-values ≤4.7×10⁻⁸) after meta-analysis of Singaporean populations. Using a candidate-SNP approach, we further replicated SNPs at 4 additional loci that had been recently identified to be associated with serum CRP (IL6R, GCKR, IL6 and IL1F10) (p-values ≤0.009), in the Singaporean datasets. SNPs from these 8 loci explained 4.05% of variance in serum CRP. Two SNPs (rs2847281 and rs6901250) were detected to be significant (p-value ≤0.036) but with opposite effect directions in the Singaporean populations as compared to original European studies. At these loci we did not detect significant inter-population LD differences. We further did not observe a significant association between CRP variant and CRVE or CRAE levels after meta-analysis of all Singaporean and European datasets (p-value >0.058).

Conclusions

Common variants associated with serum CRP, first detected in primarily European studies, are also associated with CRP levels in East-Asian and South-Asian populations. We did not find a causal link between CRP and retinal measures of microvascular disease.

High trans-ethnic replicability of GWAS results implies common causal variants

Genome-wide association studies (GWAS) have detected many disease associations. However, the reported variants tend to explain small fractions of risk, and there are doubts about issues such as the portability of findings over different ethnic groups or the relative roles of rare versus common variants in the genetic architecture of complex disease. Studying the degree of sharing of disease-associated variants across populations can help in solving these issues. We present a comprehensive survey of GWAS replicability across 28 diseases. Most loci and SNPs discovered in Europeans for these conditions have been extensively replicated using peoples of European and East Asian ancestry, while the replication with individuals of African ancestry is much less common. We found a strong and significant correlation of Odds Ratios across Europeans and East Asians, indicating that underlying causal variants are common and shared between the two ancestries. Moreover, SNPs that failed to replicate in East Asians map into genomic regions where Linkage Disequilibrium patterns differ significantly between populations. Finally, we observed that GWAS with larger sample sizes have detected variants with weaker effects rather than with lower frequencies. Our results indicate that most GWAS results are due to common variants. In addition, the sharing of disease alleles and the high correlation in their effect sizes suggest that most of the underlying causal variants are shared between Europeans and East Asians and that they tend to map close to the associated marker SNPs.

Describing and identifying the genetic variants that increase risk for complex diseases remains a central focus of human genetics and is fundamental for the emergent field of personalized medicine. Over the last six years, GWAS have revolutionized the field, discovering hundreds of disease loci. However, with only a handful of exceptions, the causal variants that generate the associations unveiled by GWAS have not been identified, and their frequency and degree of sharing across populations remains unknown. Here, we present a comprehensive comparison of GWAS results designed to try to understand the nature of causal variants. By examining the results of GWAS for 28 diseases that have been performed with peoples of European, East Asian, and African ancestries, we conclude that a large fraction of associations are caused by common causal variants that should map relatively close to the associated markers. Our results indicate that many of the disease risk variants discovered by GWAS are shared across Eurasians.

Genome-wide association study in a Chinese population identifies a susceptibility locus for type 2 diabetes at 7q32 near PAX4

AIMS/HYPOTHESIS

Most genetic variants identified for type 2 diabetes have been discovered in European populations. We performed genome-wide association studies (GWAS) in a Chinese population with the aim of identifying novel variants for type 2 diabetes in Asians.

METHODS

We performed a meta-analysis of three GWAS comprising 684 patients with type 2 diabetes and 955 controls of Southern Han Chinese descent. We followed up the top signals in two independent Southern Han Chinese cohorts (totalling 10,383 cases and 6,974 controls), and performed in silico replication in multiple populations.

RESULTS

We identified CDKN2A/B and four novel type 2 diabetes association signals with p < 1 × 10(-5) from the meta-analysis. Thirteen variants within these four loci were followed up in two independent Chinese cohorts, and rs10229583 at 7q32 was found to be associated with type 2 diabetes in a combined analysis of 11,067 cases and 7,929 controls (p meta = 2.6 × 10(-8); OR [95% CI] 1.18 [1.11, 1.25]). In silico replication revealed consistent associations across multiethnic groups, including five East Asian populations (p meta = 2.3 × 10(-10)) and a population of European descent (p = 8.6 × 10(-3)). The rs10229583 risk variant was associated with elevated fasting plasma glucose, impaired beta cell function in controls, and an earlier age at diagnosis for the cases. The novel variant lies within an islet-selective cluster of open regulatory elements. There was significant heterogeneity of effect between Han Chinese and individuals of European descent, Malaysians and Indians.

CONCLUSIONS/INTERPRETATION

Our study identifies rs10229583 near PAX4 as a novel locus for type 2 diabetes in Chinese and other populations and provides new insights into the pathogenesis of type 2 diabetes.

Detecting loci under recent positive selection in dairy and beef cattle by combining different genome-wide scan methods

As the methodologies available for the detection of positive selection from genomic data vary in terms of assumptions and execution, weak correlations are expected among them. However, if there is any given signal that is consistently supported across different methodologies, it is strong evidence that the locus has been under past selection. In this paper, a straightforward frequentist approach based on the Stouffer Method to combine P-values across different tests for evidence of recent positive selection in common variations, as well as strategies for extracting biological information from the detected signals, were described and applied to high density single nucleotide polymorphism (SNP) data generated from dairy and beef cattle (taurine and indicine). The ancestral Bovinae allele state of over 440,000 SNP is also reported. Using this combination of methods, highly significant (P<3.17×10⁻⁷) population-specific sweeps pointing out to candidate genes and pathways that may be involved in beef and dairy production were identified. The most significant signal was found in the Cornichon homolog 3 gene (CNIH3) in Brown Swiss (P = 3.82×10⁻¹²), and may be involved in the regulation of pre-ovulatory luteinizing hormone surge. Other putative pathways under selection are the glucolysis/gluconeogenesis, transcription machinery and chemokine/cytokine activity in Angus; calpain-calpastatin system and ribosome biogenesis in Brown Swiss; and gangliosides deposition in milk fat globules in Gyr. The composite method, combined with the strategies applied to retrieve functional information, may be a useful tool for surveying genome-wide selective sweeps and providing insights in to the source of selection.

Genomic study in Mexicans identifies a new locus for triglycerides and refines European lipid loci

BACKGROUND

The Mexican population and others with Amerindian heritage exhibit a substantial predisposition to dyslipidemias and coronary heart disease. Yet, these populations remain underinvestigated by genomic studies, and to date, no genome-wide association (GWA) studies have been reported for lipids in these rapidly expanding populations.

METHODS AND FINDINGS

We performed a two-stage GWA study for hypertriglyceridemia and low high-density lipoprotein cholesterol (HDL-C) in Mexicans (n=4361), and identified a novel Mexican-specific genome-wide significant locus for serum triglycerides (TGs) near the Niemann-Pick type C1 protein gene (p=2.43×10(-08)). Furthermore, three European loci for TGs (APOA5, GCKR and LPL), and four loci for HDL-C (ABCA1, CETP, LIPC and LOC55908) reached genome-wide significance in Mexicans. We used cross-ethnic mapping to narrow three European TG GWA loci, APOA5, MLXIPL, and CILP2 that were wide and contained multiple candidate variants in the European scan. At the APOA5 locus, this reduced the most likely susceptibility variants to one, rs964184. Importantly, our functional analysis demonstrated a direct link between rs964184 and postprandial serum apoAV protein levels, supporting rs964184 as the causative variant underlying the European and Mexican GWA signal. Overall, 52 of the 100 reported associations from European lipid GWA meta-analysis generalised to Mexicans. However, in 82 of the 100 European GWA loci, a different variant other than the European lead/best-proxy variant had the strongest regional evidence of association in Mexicans.

CONCLUSIONS

This first Mexican GWA study of lipids identified a novel GWA locus for high TG levels; used the interpopulation heterogeneity to significantly restrict three previously known European GWA signals, and surveyed whether the European lipid GWA SNPs extend to the Mexican population.

Transferability and fine mapping of type 2 diabetes loci in African Americans: the Candidate Gene Association Resource Plus Study

Type 2 diabetes (T2D) disproportionally affects African Americans (AfA) but, to date, genetic variants identified from genome-wide association studies (GWAS) are primarily from European and Asian populations. We examined the single nucleotide polymorphism (SNP) and locus transferability of 40 reported T2D loci in six AfA GWAS consisting of 2,806 T2D case subjects with or without end-stage renal disease and 4,265 control subjects from the Candidate Gene Association Resource Plus Study. Our results revealed that seven index SNPs at the TCF7L2, KLF14, KCNQ1, ADCY5, CDKAL1, JAZF1, and GCKR loci were significantly associated with T2D (P < 0.05). The strongest association was observed at TCF7L2 rs7903146 (odds ratio [OR] 1.30; P = 6.86 × 10⁻⁸). Locus-wide analysis demonstrated significant associations (P(emp) < 0.05) at regional best SNPs in the TCF7L2, KLF14, and HMGA2 loci as well as suggestive signals in KCNQ1 after correction for the effective number of SNPs at each locus. Of these loci, the regional best SNPs were in differential linkage disequilibrium (LD) with the index and adjacent SNPs. Our findings suggest that some loci discovered in prior reports affect T2D susceptibility in AfA with similar effect sizes. The reduced and differential LD pattern in AfA compared with European and Asian populations may facilitate fine mapping of causal variants at loci shared across populations.

A genome-wide association study identifies GRK5 and RASGRP1 as type 2 diabetes loci in Chinese Hans

Substantial progress has been made in identification of type 2 diabetes (T2D) risk loci in the past few years, but our understanding of the genetic basis of T2D in ethnically diverse populations remains limited. We performed a genome-wide association study and a replication study in Chinese Hans comprising 8,569 T2D case subjects and 8,923 control subjects in total, from which 10 single nucleotide polymorphisms were selected for further follow-up in a de novo replication sample of 3,410 T2D case and 3,412 control subjects and an in silico replication sample of 6,952 T2D case and 11,865 control subjects. Besides confirming seven established T2D loci (CDKAL1, CDKN2A/B, KCNQ1, CDC123, GLIS3, HNF1B, and DUSP9) at genome-wide significance, we identified two novel T2D loci, including G-protein-coupled receptor kinase 5 (GRK5) (rs10886471: P = 7.1 × 10(-9)) and RASGRP1 (rs7403531: P = 3.9 × 10(-9)), of which the association signal at GRK5 seems to be specific to East Asians. In nondiabetic individuals, the T2D risk-increasing allele of RASGRP1-rs7403531 was also associated with higher HbA(1c) and lower homeostasis model assessment of β-cell function (P = 0.03 and 0.0209, respectively), whereas the T2D risk-increasing allele of GRK5-rs10886471 was also associated with higher fasting insulin (P = 0.0169) but not with fasting glucose. Our findings not only provide new insights into the pathophysiology of T2D, but may also shed light on the ethnic differences in T2D susceptibility.

Transferability and fine-mapping of glucose and insulin quantitative trait loci across populations: CARe, the Candidate Gene Association Resource

AIMS/HYPOTHESIS

Hyperglycaemia disproportionately affects African-Americans (AfAs). We tested the transferability of 18 single-nucleotide polymorphisms (SNPs) associated with glycaemic traits identified in European ancestry (EuA) populations in 5,984 non-diabetic AfAs.

METHODS

We meta-analysed SNP associations with fasting glucose (FG) or insulin (FI) in AfAs from five cohorts in the Candidate Gene Association Resource. We: (1) calculated allele frequency differences, variations in linkage disequilibrium (LD), fixation indices (F(st)s) and integrated haplotype scores (iHSs); (2) tested EuA SNPs in AfAs; and (3) interrogated within ± 250 kb around each EuA SNP in AfAs.

RESULTS

Allele frequency differences ranged from 0.6% to 54%. F(st) exceeded 0.15 at 6/16 loci, indicating modest population differentiation. All iHSs were <2, suggesting no recent positive selection. For 18 SNPs, all directions of effect were the same and 95% CIs of association overlapped when comparing EuA with AfA. For 17 of 18 loci, at least one SNP was nominally associated with FG in AfAs. Four loci were significantly associated with FG (GCK, p = 5.8 × 10(-8); MTNR1B, p = 8.5 × 10(-9); and FADS1, p = 2.2 × 10(-4)) or FI (GCKR, p = 5.9 × 10(-4)). At GCK and MTNR1B the EuA and AfA SNPs represented the same signal, while at FADS1, and GCKR, the EuA and best AfA SNPs were weakly correlated (r(2) <0.2), suggesting allelic heterogeneity for association with FG at these loci.

CONCLUSIONS/INTERPRETATION

Few glycaemic SNPs showed strict evidence of transferability from EuA to AfAs. Four loci were significantly associated in both AfAs and those with EuA after accounting for varying LD across ancestral groups, with new signals emerging to aid fine-mapping.

Regions of homozygosity in three Southeast Asian populations

The genomes of outbred populations were first shown in 2006 to contain regions of homozygosity (ROHs) of several megabases. Further studies have also investigated the characteristics of ROHs in healthy individuals in various populations but there are no studies on Singapore populations to date. This study aims to identify and investigate the characteristics of ROHs in three Singapore populations. A total of 268 samples (96 Chinese, 89 Malays and 83 Indians) are genotyped on Illumina Human 1 M Beadchip and Affymetrix Genome-Wide Human SNP Array 6.0. We use the PennCNV algorithm to detect ROHs. We report an abundance of ROHs (≥500 kb), with an average of more than one hundred regions per individual. On average, the Indian population has the lowest number of ROHs and smallest total length of ROHs per individual compared with the Chinese and Malay populations. We further investigate the relationship between the occurrence of ROHs and haplotype frequency, regional linkage disequilibrium (LD) and positive selection. Based on the results of this data set, we find that the frequency of occurrence of ROHs is positively associated with haplotype frequency and regional LD. The majority of regions detected for recent positive selection and regions with differential LD between populations overlap with the ROH loci. When we consider both the location of the ROHs and the allelic form of the ROHs, we are able to separate the populations by principal component analysis, demonstrating that ROHs contain information on population structure and the demographic history of a population.

Genome-wide association study in individuals of South Asian ancestry identifies six new type 2 diabetes susceptibility loci

We carried out a genome-wide association study of type-2 diabetes (T2D) in individuals of South Asian ancestry. Our discovery set included 5,561 individuals with T2D (cases) and 14,458 controls drawn from studies in London, Pakistan and Singapore. We identified 20 independent SNPs associated with T2D at P < 10(-4) for testing in a replication sample of 13,170 cases and 25,398 controls, also all of South Asian ancestry. In the combined analysis, we identified common genetic variants at six loci (GRB14, ST6GAL1, VPS26A, HMG20A, AP3S2 and HNF4A) newly associated with T2D (P = 4.1 × 10(-8) to P = 1.9 × 10(-11)). SNPs at GRB14 were also associated with insulin sensitivity (P = 5.0 × 10(-4)), and SNPs at ST6GAL1 and HNF4A were also associated with pancreatic beta-cell function (P = 0.02 and P = 0.001, respectively). Our findings provide additional insight into mechanisms underlying T2D and show the potential for new discovery from genetic association studies in South Asians, a population with increased susceptibility to T2D.

High-throughput single-nucleotide polymorphism discovery and the search for candidate genes for long-term SIVmac nonprogression in Chinese rhesus macaques (Macaca mulatta)

BACKGROUND

Genetic differences between Indian and Chinese rhesus macaques contribute to the phenotypic variance of clinical trials, including infection with SIVmac. The completion of the rhesus genome has facilitated the discovery of several thousand markers.

METHODS

We developed a genome-wide SNP map for rhesus macaques containing 3869 validated markers with an average distance of 0.88 Mb and used the program VarLD to identify genomic areas with significant differences in linkage disequilibrium (LD) between Indian-derived and Chinese rhesus macaques.

RESULTS

Forty-one statistically significant differences in LD between Chinese and Indian-origin rhesus were detected on chromosomes 1, 4, 5 and 11. The region of greatest LD difference was located on the proximal end of chromosome one, which also contained the genes ELAVL4, MAST2 and HIVEP3.

CONCLUSION

These genomic areas provide entry to more detailed studies of gene function. This method is also applicable to the study of differences in biomarkers between regional populations of other species.

Meta-analysis of genome-wide association studies identifies common variants associated with blood pressure variation in east Asians

We conducted a meta-analysis of genome-wide association studies of systolic (SBP) and diastolic (DBP) blood pressure in 19,608 subjects of east Asian ancestry from the AGEN-BP consortium followed up with de novo genotyping (n = 10,518) and further replication (n = 20,247) in east Asian samples. We identified genome-wide significant (P < 5 × 10(-8)) associations with SBP or DBP, which included variants at four new loci (ST7L-CAPZA1, FIGN-GRB14, ENPEP and NPR3) and a newly discovered variant near TBX3. Among the five newly discovered variants, we obtained significant replication in the independent samples for all of these loci except NPR3. We also confirmed seven loci previously identified in populations of European descent. Moreover, at 12q24.13 near ALDH2, we observed strong association signals (P = 7.9 × 10(-31) and P = 1.3 × 10(-35) for SBP and DBP, respectively) with ethnic specificity. These findings provide new insights into blood pressure regulation and potential targets for intervention.

Transferability of type 2 diabetes implicated loci in multi-ethnic cohorts from Southeast Asia

Recent large genome-wide association studies (GWAS) have identified multiple loci which harbor genetic variants associated with type 2 diabetes mellitus (T2D), many of which encode proteins not previously suspected to be involved in the pathogenesis of T2D. Most GWAS for T2D have focused on populations of European descent, and GWAS conducted in other populations with different ancestry offer a unique opportunity to study the genetic architecture of T2D. We performed genome-wide association scans for T2D in 3,955 Chinese (2,010 cases, 1,945 controls), 2,034 Malays (794 cases, 1,240 controls), and 2,146 Asian Indians (977 cases, 1,169 controls). In addition to the search for novel variants implicated in T2D, these multi-ethnic cohorts serve to assess the transferability and relevance of the previous findings from European descent populations in the three major ethnic populations of Asia, comprising half of the world's population. Of the SNPs associated with T2D in previous GWAS, only variants at CDKAL1 and HHEX/IDE/KIF11 showed the strongest association with T2D in the meta-analysis including all three ethnic groups. However, consistent direction of effect was observed for many of the other SNPs in our study and in those carried out in European populations. Close examination of the associations at both the CDKAL1 and HHEX/IDE/KIF11 loci provided some evidence of locus and allelic heterogeneity in relation to the associations with T2D. We also detected variation in linkage disequilibrium between populations for most of these loci that have been previously identified. These factors, combined with limited statistical power, may contribute to the failure to detect associations across populations of diverse ethnicity. These findings highlight the value of surveying across diverse racial/ethnic groups towards the fine-mapping efforts for the casual variants and also of the search for variants, which may be population-specific.

Type 2 diabetes mellitus (T2D) is a chronic disease which can lead to complications such as heart disease, stroke, hypertension, blindness due to diabetic retinopathy, amputations from peripheral vascular diseases, and kidney disease from diabetic nephropathy. The increasing prevalence and complications of T2D are likely to increase the health and economic burden of individuals, families, health systems, and countries. Our study carried out in three major Asian ethnic groups (Chinese, Malays, and Indians) in Singapore suggests that the findings of studies carried out in populations of European ancestry (which represents most studies to date) may be relevant to populations in Asia. However, our study also raises the possibility that different genes, and within the genes different variants, may confer susceptibility to T2D in these populations. These findings are particularly relevant in Asia, where the greatest growth of T2D is expected in the coming years, and emphasize the importance of studying diverse populations when trying to localize the regions of the genome associated with T2D. In addition, we may need to consider novel methods for combining data across populations.