Linkage Disequilibrium Diversity Research Articles

Osteoporosis GWAS-implicated DNM3 locus contextually regulates osteoblastic and chondrogenic fate of mesenchymal stem/progenitor cells through oscillating miR-199a-5p levels.

Genome wide association study (GWAS)-implicated bone mineral density (BMD) signals have been shown to localize in cis-regulatory regions of distal effector genes using 3D genomic methods. Detailed characterization of such genes can reveal novel causal genes for BMD determination. Here, we elected to characterize the "DNM3" locus on chr1q24, where the long non-coding RNA DNM3OS and the embedded microRNA MIR199A2 (miR-199a-5p) are implicated as effector genes contacted by the region harboring variation in linkage disequilibrium with BMD-associated sentinel single nucleotide polymorphism, rs12041600. During osteoblast differentiation of human mesenchymal stem/progenitor cells (hMSC), miR-199a-5p expression was temporally decreased and correlated with the induction of osteoblastic transcription factors RUNX2 and Osterix. Functional relevance of miR-199a-5p downregulation in osteoblastogenesis was investigated by introducing miR-199a-5p mimic into hMSC. Cells overexpressing miR-199a-5p depicted a cobblestone-like morphological change and failed to produce BMP2-dependent extracellular matrix mineralization. Mechanistically, a miR-199a-5p mimic modified hMSC propagated normal SMAD1/5/9 signaling and expressed osteoblastic transcription factors RUNX2 and Osterix but depicted pronounced upregulation of SOX9 and enhanced expression of essential chondrogenic genes ACAN, COMP, and COL10A1. Mineralization defects, morphological changes, and enhanced chondrogenic gene expression associated with miR-199a-5p mimic over-expression were restored with miR-199a-5p inhibitor suggesting specificity of miR-199a-5p in chondrogenic fate specification. The expression of both the DNM3OS and miR-199a-5p temporally increased and correlated with hMSC chondrogenic differentiation. Although miR-199a-5p overexpression failed to further enhance chondrogenesis, blocking miR-199a-5p activity significantly reduced chondrogenic pellet size, extracellular matrix deposition, and chondrogenic gene expression. Taken together, our results indicate that oscillating miR-199a-5p levels dictate hMSC osteoblast or chondrocyte terminal fate. Our study highlights a functional role of miR-199a-5p as a BMD effector gene at the DNM3 BMD GWAS locus, where patients with cis-regulatory genetic variation which increases miR-199a-5p expression could lead to reduced osteoblast activity.

Establishment of NGS-based HLA 9-locus haplotypes in the Eastern Han Chinese population highlights the role of HLA-DP in donor selection for transplantation.

We collected HLA typing data from 653 families in the Eastern Han Chinese population. HLA-A, B, C, DRB1, DRB3, DRB4, DRB5, DQA1, DQB1, DPA1, and DPB1 (HLA-11 loci) typing of 1781 subjects was performed using a commercial next-generation sequencing (NGS) method in our laboratory. The phasing of haplotypes in each family was determined by Mendelian segregation. Haplotype analysis revealed 1634 different haplotypes among a total of 2230 haplotypes. The predominant haplotype was A*30:01-C*06:02-B*13:02-DRB1*07:01-DRB4*01:03-DQA1*02:01-DQB1*02:02-DPA1*02:01-DPB1*17:01 (HF=4.04%), followed by A*02:07-C*01:02-B*46:01-DRB1*09:01-DRB4*01:03-DQA1*03:02-DQB1*03:03-DPA1*02:02-DPB1*05:01 (HF=1.84%) and A*33:03-C*03:02-B*58:01-DRB1*03:01-DRB3*02:02-DQA1*05:01-DQB1*02:01-DPA1*01:03-DPB1*04:01 (HF = 1.48%), accounting for 7.35% of the total. Meanwhile 76.41% of all haplotypes were observed only once or twice (HF < 0.1%). Different from HLA-DRB3/4/5 and DQA1 loci, DP linkage markedly increased haplotype variation by 34.82% based on the 5-locus haplotype. The much weaker linkage disequilibrium (LD) of DQB1-DPB1 indicated the reason. We observed 10 analyzable recombination events, most of which occurred at DP loci. Even with the same common 5-locus haplotype, HLA-DP linkage alters the haplotype diversity and frequency. Analysis of related haplotype assignment and unrelated recipient-donor pairs matching at the 9-locus haplotype revealed that HLA-DP affects the donor selection strategy. Haplotype study of a large sample size using NGS identified linkage haplotypes beyond the 5 loci. LD, recombination events, and haplotype variation caused by DP loci emphasized that HLA 9-locus haplotype matching should be considered in donor selection, particularly the effect of DP loci. The finding lays the foundation for further studies on the effect of HLA-DP mismatch on transplantation.

Genome-wide estimates of genetic diversity, inbreeding and effective size of experimental and commercial rainbow trout lines undergoing selective breeding

BackgroundSelective breeding is a relatively recent practice in aquaculture species compared to terrestrial livestock. Nevertheless, the genetic variability of farmed salmonid lines, which have been selected for several generations, should be assessed. Indeed, a significant decrease in genetic variability due to high selection intensity could have occurred, potentially jeopardizing the long-term genetic progress as well as the adaptive capacities of populations facing change(s) in the environment. Thus, it is important to evaluate the impact of selection practices on genetic diversity to limit future inbreeding. The current study presents an analysis of genetic diversity within and between six French rainbow trout (Oncorhynchus mykiss) experimental or commercial lines based on a medium-density single nucleotide polymorphism (SNP) chip and various molecular genetic indicators: fixation index (FST), linkage disequilibrium (LD), effective population size (Ne) and inbreeding coefficient derived from runs of homozygosity (ROH).ResultsOur results showed a moderate level of genetic differentiation between selected lines (FST ranging from 0.08 to 0.15). LD declined rapidly over the first 100 kb, but then remained quite high at long distances, leading to low estimates of Ne in the last generation ranging from 24 to 68 depending on the line and methodology considered. These results were consistent with inbreeding estimates that varied from 10.0% in an unselected experimental line to 19.5% in a commercial line, and which are clearly higher than corresponding estimates in ruminants or pigs. In addition, strong variations in LD and inbreeding were observed along the genome that may be due to differences in local rates of recombination or due to key genes that tended to have fixed favorable alleles for domestication or production.ConclusionsThis is the first report on ROH for any aquaculture species. Inbreeding appeared to be moderate to high in the six French rainbow trout lines, due to founder effects at the start of the breeding programs, but also likely to sweepstakes reproductive success in addition to selection for the selected lines. Efficient management of inbreeding is a major goal in breeding programs to ensure that populations can adapt to future breeding objectives and SNP information can be used to manage the rate at which inbreeding builds up in the fish genome.

Extensive Ethnic Variation and Linkage Disequilibrium at the FCGR2/3 Locus: Different Genetic Associations Revealed in Kawasaki Disease.

The human Fc-gamma receptors (FcγRs) link adaptive and innate immunity by binding immunoglobulin G (IgG). All human low-affinity FcγRs are encoded by the FCGR2/3 locus containing functional single nucleotide polymorphisms (SNPs) and gene copy number variants. This locus is notoriously difficult to genotype and high-throughput methods commonly used focus on only a few SNPs. We performed multiplex ligation-dependent probe amplification for all relevant genetic variations at the FCGR2/3 locus in >4,000 individuals to define linkage disequilibrium (LD) and allele frequencies in different populations. Strong LD and extensive ethnic variation in allele frequencies was found across the locus. LD was strongest for the FCGR2C-ORF haplotype (rs759550223+rs76277413), which leads to expression of FcγRIIc. In Europeans, the FCGR2C-ORF haplotype showed strong LD with, among others, rs201218628 (FCGR2A-Q27W, r2 = 0.63). LD between these two variants was weaker (r2 = 0.17) in Africans, whereas the FCGR2C-ORF haplotype was nearly absent in Asians (minor allele frequency <0.005%). The FCGR2C-ORF haplotype and rs1801274 (FCGR2A-H131R) were in weak LD (r2 = 0.08) in Europeans. We evaluated the importance of ethnic variation and LD in Kawasaki Disease (KD), an acute vasculitis in children with increased incidence in Asians. An association of rs1801274 with KD was previously shown in ethnically diverse genome-wide association studies. Now, we show in 1,028 European KD patients that the FCGR2C-ORF haplotype, although nearly absent in Asians, was more strongly associated with susceptibility to KD than rs1801274 in Europeans. Our data illustrate the importance of interpreting findings of association studies concerning the FCGR2/3 locus with knowledge of LD and ethnic variation.

Identification of positive selection signatures in pigs by comparing linkage disequilibrium variances.

Selection affects the patterns of linkage disequilibrium (LD) around the site of a beneficial allele with an increase in LD among the hitchhiking alleles. Comparing the differences in regional LD between pig populations could help to identify putative genomic regions with potential adaptations for economic traits. In this study, using Illumina Porcine SNP60K BeadChip genotyping data from 207 Chinese indigenous, 117 South American village and 408 Large White pigs, we estimated the variation of genome-wide LD between populations using the varld program. The top 0.1% standardized VarLD scores were used as a criterion for all comparisons, and compared with LD blocks, a total of four selection signatures on Sus scrofa chromosome (SSC) 7, 9, 13 and 14 were identified in all populations. These signatures overlapped with quantitative trait loci for linoleic acid content, age at puberty, number of muscle fibers per unit area, hip structure and body weight traits in pigs. Among them, one of the signatures (56.5-56.6Mb on SSC7) in Large White pigs harbored the ADAMTSL3 gene, which is known to affect body length. The findings of this study seem to point toward recent selection in different pig populations. Further investigations are encouraged to confirm the selection signatures detected by varld in the present study.

Evaluating signatures of selection through variation in linkage disequilibrium among different cattle breeds

The aim of presented investigation was to assess the overlap between selection signature discovered through analysis of variation in linkage disequilibrium and reported genomic regions associated with economic and traits of biological importance in cattle populations. The differences across Slovak Pinzgau, Austrian Pinzgau, Simmental and Holstein cattle and thus genome signatures of production and adaptation were found. The highest peak (top 0.01 percentile) was observed between Slovak and Austrian on chromosome 23, between Slovak Pinzgau and Simmental on chromosome 4 and between Slovak Pinzgau and Holstein on chromosomes 1, 7 and 20. Many candidate genes found have a known role in milk production (casein genes CSN1S1, CSN2, CSN1S2, CSN3; ABCG2, HBEGF, CAPN3, DGAT1, TG, GHR), reproduction (MGAT1, FGF1), feed efficiency (R3HDM1, ZRANB3), fertility (SPOCK1) and immune response (HSPA9, CD14, ARAP3, PCDH). Results of this study could be the basis for implementation of genomic selection programs in the Slovak Pinzgau cattle.

Temporal dynamics of linkage disequilibrium in two populations of bighorn sheep.

Linkage disequilibrium (LD) is the nonrandom association of alleles at two markers. Patterns of LD have biological implications as well as practical ones when designing association studies or conservation programs aimed at identifying the genetic basis of fitness differences within and among populations. However, the temporal dynamics of LD in wild populations has received little empirical attention. In this study, we examined the overall extent of LD, the effect of sample size on the accuracy and precision of LD estimates, and the temporal dynamics of LD in two populations of bighorn sheep (Ovis canadensis) with different demographic histories. Using over 200 microsatellite loci, we assessed two metrics of multi-allelic LD, D′, and χ′2. We found that both populations exhibited high levels of LD, although the extent was much shorter in a native population than one that was founded via translocation, experienced a prolonged bottleneck post founding, followed by recent admixture. In addition, we observed significant variation in LD in relation to the sample size used, with small sample sizes leading to depressed estimates of the extent of LD but inflated estimates of background levels of LD. In contrast, there was not much variation in LD among yearly cross-sections within either population once sample size was accounted for. Lack of pronounced interannual variability suggests that researchers may not have to worry about interannual variation when estimating LD in a population and can instead focus on obtaining the largest sample size possible.

Bridging Integrator 1 (BIN1) Genotype Effects on Working Memory, Hippocampal Volume, and Functional Connectivity in Young Healthy Individuals.

Alzheimer's disease (AD) is the most common form of dementia and exhibits a considerable level of heritability. The bridging integrator 1 (BIN1) gene has recently been identified in several large genome-wide association studies (GWAS) as the second most important risk locus for AD following apolipoprotein E (APOE). However, how and when the established genetic risk locus BIN1 rs744373 confers risk to late-onset AD has yet to be determined. Here using an imaging genetic strategy in large-sample Chinese subjects, we show that healthy homozygous carriers of the rs744373 risk allele exhibit worse high-load working memory (WM) performance, larger hippocampal volume and lower functional connectivity between the bilateral hippocampus and the right dorsolateral prefrontal cortex (DLPFC), mirroring clinical evidence of disturbed memory and connectivity in patients. Our findings demonstrate that rs744373 itself or a variation in linkage disequilibrium may provide a neurogenetic mechanism for BIN1 while further validating the possibility of combining genetic and neuroimaging strategies to monitor individuals at risk for AD.

Genome-wide linkage disequilibrium in nine-spined stickleback populations.

Variation in the extent and magnitude of genome-wide linkage disequilibrium (LD) among populations residing in different habitats has seldom been studied in wild vertebrates. We used a total of 109 microsatellite markers to quantify the level and patterns of genome-wide LD in 13 Fennoscandian nine-spined stickleback (Pungitius pungitius) populations from four (viz. marine, lake, pond, and river) different habitat types. In general, high magnitude (D’ > 0.5) of LD was found both in freshwater and marine populations, and the magnitude of LD was significantly greater in inland freshwater than in marine populations. Interestingly, three coastal freshwater populations located in close geographic proximity to the marine populations exhibited similar LD patterns and genetic diversity as their marine neighbors. The greater levels of LD in inland freshwater compared with marine and costal freshwater populations can be explained in terms of their contrasting demographic histories: founder events, long-term isolation, small effective sizes, and population bottlenecks are factors likely to have contributed to the high levels of LD in the inland freshwater populations. In general, these findings shed new light on the patterns and extent of variation in genome-wide LD, as well as the ecological and evolutionary factors driving them.

Targeted next-generation sequencing and fine linkage disequilibrium mapping reveals association of PNPLA3 and PARVB with the severity of nonalcoholic fatty liver disease

The genomic regions containing PNPLA3, SAMM50 and PARVB are susceptibility loci for the development and progression of nonalcoholic fatty liver disease (NAFLD). In order to search for all common variations in this region, we amplified the genomic DNA of 28 NAFLD patients by long-range PCR, covering the entire susceptibility region and sequenced the DNA using indexed multiplex next-generation sequencing. We found 329 variations, including four novel variations. Fine mapping of variations including insertion/deletions was performed for 540 NAFLD patients (488 with nonalcoholic steatohepatitis (NASH) and 52 with simple steatosis) and 1012 control subjects. HaploView analysis showed that linkage disequilibrium (LD) block 1 and 2 occurred in PNPLA3, block 3 in SAMM50 and block 4 in PARVB. Variations in LD blocks 1-4 were significantly associated with NAFLD as compared with control subjects (P<1 × 10(-8)). Variations in LD block 2 were significantly associated with the NAFLD activity score (NAS), aspartate aminotransferase and alanine aminotransferase. Variations in LD block 1 were significantly associated with the fibrosis stage. The strongest associations were observed for variations in LD block 4, with NASH as compared with simple steatosis (P=7.1 × 10(-6)) and NAS (P=3.4 × 10(-6)). Our results suggested that variations, including insertion/deletions, in PARVB, as well as those in PNPLA3, are important in the progression of NAFLD.

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

BackgroundSignatures 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.MethodsBy 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.ResultsFor 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.ConclusionsThe 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.

Peach genetic resources: diversity, population structure and linkage disequilibrium

BackgroundPeach (Prunus persica (L.) Batsch) is one of the most important model fruits in the Rosaceae family. Native to the west of China, where peach has been domesticated for more than 4,000 years, its cultivation spread from China to Persia, Mediterranean countries and to America. Chinese peach has had a major impact on international peach breeding programs due to its high genetic diversity. In this research, we used 48 highly polymorphic SSRs, distributed over the peach genome, to investigate the difference in genetic diversity, and linkage disequilibrium (LD) among Chinese cultivars, and North American and European cultivars, and the evolution of current peach cultivars.ResultsIn total, 588 alleles were obtained with 48 SSRs on 653 peach accessions, giving an average of 12.25 alleles per locus. In general, the average value of observed heterozygosity (0.47) was lower than the expected heterozygosity (0.60). The separate analysis of groups of accessions according to their origin or reproductive strategies showed greater variability in Oriental cultivars, mainly due to the high level of heterozygosity in Chinese landraces. Genetic distance analysis clustered the cultivars into two main groups: one included four wild related Prunus, and the other included most of the Oriental and Occidental landraces and breeding cultivars. STRUCTURE analysis assigned 469 accessions to three subpopulations: Oriental (234), Occidental (174), and Landraces (61). Nested STRUCTURE analysis divided the Oriental subpopulation into two different subpopulations: ‘Yu Lu’ and ‘Hakuho’. The Occidental breeding subpopulation was also subdivided into nectarine and peach subpopulations. Linkage disequilibrium (LD) analysis in each of these subpopulations showed that the percentage of linked (r2 > 0.1) intra-chromosome comparisons ranged between 14% and 47%. LD decayed faster in Oriental (1,196 Kbp) than in Occidental (2,687 Kbp) samples. In the ‘Yu Lu’ subpopulation there was considerable LD extension while no variation of LD with physical distance was observed in the landraces. From the first STRUCTURE result, LG1 had the greatest proportion of alleles in LD within all three subpopulations.ConclusionsOur study demonstrates a high level of genetic diversity and relatively fast decay of LD in the Oriental peach breeding program. Inclusion of Chinese landraces will have a greater effect on increasing genetic diversity in Occidental breeding programs. Fingerprinting with genotype data for all 658 cultivars will be used for accession management in different germplasms. A higher density of markers are needed for association mapping in Oriental germplasm due to the low extension of LD. Population structure and evaluation of LD provides valuable information for GWAS experiment design in peach.

CEBPG Exhibits Allele-Specific Expression in Human Bronchial Epithelial Cells

Inter-individual variation in CCAAT/enhancer binding protein gamma (CEBPG) transcript expression in normal human bronchial epithelial cells (NBEC) is associated with predisposition to lung cancer. We hypothesize that this inter-individual variation is in part explained by cis-acting genetic variation in CEBPG. To test this hypothesis we measured transcript expression derived from each parental copy of CEBPG (ie, allele-specific expression; ASE). There was a significant 2.9-fold higher cell cycle-specific variation in ASE of CEBPG rs2772 A compared to C allele (P < 0.001). In 20% of NBEC samples, CEBPG rs2772 A allele was expressed on average 2.10 fold greater than rs2772 C allele. These data support the hypothesis that genetic variation in linkage disequilibrium with rs2772 influences regulation of CEBPG transcript expression through a trans-effect downstream of RNA polymerase II transcription and confirm that cis-acting genetic variation contributes to inter-individual variation in CEBPG transcript expression in NBEC, which is associated with variation in lung cancer risk.

Similarity in recombination rate and linkage disequilibrium at CYP2C and CYP2D cytochrome P450 gene regions among Europeans indicates signs of selection and no advantage of using tagSNPs in population isolates

Linkage disequilibrium (LD) and recombination rate variations are known to vary considerably between human genome regions and populations mostly because of the combined effects of mutation, recombination, and demographic history. Thus, the pattern of LD is a key issue to disentangle variants associated with complex traits. Here, we aim to describe the haplotype structure and LD variation at the pharmacogenetically relevant cytochrome P450 CYP2C and CYP2D gene regions among European populations. To assess the haplotype structure, LD pattern, and recombination rate variations in the clinically significant CYP2C and CYP2D regions, we genotyped 143 single-nucleotide polymorphisms (SNPs) across these two genome regions in a diverse set of 11 European population samples and one sub-Saharan African sample. Our results showed extended patterns of LD and in general a low rate of recombination at these loci, and a low degree of allele differentiation for the two cytochrome P450 regions among Europeans, with the exception of the Sami and the Finns as European outliers. The Sami sample showed reduced haplotype diversity and higher LD for the two cytochrome P450 regions than the other Europeans, a feature that is proposed to enhance the LD mapping of underlying common complex traits. However, recombination hotspots and LD blocks at these two regions showed highly consistent structures across Europeans including Finns and Sami. Moreover, we showed that the CEPH sample has significantly higher tag transferability among Europeans and a more efficient tagging of both the rare CYP2C9 and the common CYP2C19 functional variants than the Sami. Our data set included CYP2C9*3 (rs1057910) and CYP2C19*2 (rs4244285) enzyme activity-altering variants associated in a recent genome-wide study with acenocoumarol-induced and warfarin-induced anticoagulation or to the antiplatelet effect of clopidogrel, respectively. Including these known activity-altering variants, we showed the haplotype variation and high derived allele frequencies of novel recently identified acenocoumarol genome-wide associated SNPs at CYP2C9 (rs4086116) and CYP2C18 (rs12772169, rs1998591, rs2104543, rs1042194) loci in a comprehensive set of 11 European populations. Furthermore, a significant frequency difference of a CYP2C19*2 gene mutation causing variable drug reactions was observed among Europeans. The CEPH sample representing the general European population as such in the HapMap project seems to be the optimal population sample for the LD mapping of common complex traits among Europeans. Nevertheless, it is still argued that the unique pattern of LD in the Sami may offer an advantage for further association mapping, especially if multiple rare variants play a role in disease etiology. However, besides the activity-altering CYP2C9*3 (rs1057910) and CYP2C19*2 (rs4244285) variants, the high derived allele frequencies of novel recently identified acenocoumarol genome-wide associated SNPs at CYP2C9 (rs4086116) and CYP2C18 (rs12772169, rs1998591, rs2104543, rs1042194) loci variants indicated that the CYP2C region may have been influenced by selection. Thus, this fine-scale haplotype map of the CYP2C and CYP2D regions may help to choose markers for further association mapping of complex pharmacogenetic traits at these loci.

Genetic variability of glutathione S-transferase enzymes in human populations: Functional inter-ethnic differences in detoxification systems

Glutathione S-Transferase enzymes (GSTs) constitute the principal Phase II superfamily which plays a key role in cellular detoxification and in other biological processes. Studies of GSTs have revealed that genetic polymorphisms are present in these enzymes and that some of these are Loss-of-Function (LoF) variants, which affect enzymatic functions and are related to different aspects of human health.The aim of this study was to analyze functional genetic differences in GST enzymes among human populations. Attention was focused on LoF polymorphisms of GSTA1, GSTM1, GSTO1, GSTO2, GSTP1 and GSTT1 genes. These LoF variants were analyzed in 668 individuals belonging to six human groups with different ethnic backgrounds: Amhara and Oromo from Ethiopia; Colorado and Cayapa Amerindians and African Ecuadorians from Ecuador; and one sample from central Italy. The HapMap database was used to compare our data with reference populations and to analyze the haplotype and Linkage Disequilibrium diversity in different ethnic groups.Our results highlighted that ethnicity strongly affects the genetic variability of GST enzymes. In particular, GST haplotypes/variants with functional impact showed significant differences in human populations, according to their ethnic background. These data underline that human populations have different structures in detoxification genes, suggesting that these ethnic differences influence disease risk or response to drugs and therefore have implications for genetic association studies involving GST enzymes.In conclusion, our investigation provides data about the distribution of important LoF variants in GST genes in human populations. This information may be useful for designing and interpreting genetic association studies.

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

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. 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. 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. 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.

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.

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.

Population- and genome-specific patterns of linkage disequilibrium and SNP variation in spring and winter wheat (Triticum aestivum L.).

BackgroundSingle nucleotide polymorphisms (SNPs) are ideally suited for the construction of high-resolution genetic maps, studying population evolutionary history and performing genome-wide association mapping experiments. Here, we used a genome-wide set of 1536 SNPs to study linkage disequilibrium (LD) and population structure in a panel of 478 spring and winter wheat cultivars (Triticum aestivum) from 17 populations across the United States and Mexico.ResultsMost of the wheat oligo pool assay (OPA) SNPs that were polymorphic within the complete set of 478 cultivars were also polymorphic in all subpopulations. Higher levels of genetic differentiation were observed among wheat lines within populations than among populations. A total of nine genetically distinct clusters were identified, suggesting that some of the pre-defined populations shared significant proportion of genetic ancestry. Estimates of population structure (FST) at individual loci showed a high level of heterogeneity across the genome. In addition, seven genomic regions with elevated FST were detected between the spring and winter wheat populations. Some of these regions overlapped with previously mapped flowering time QTL. Across all populations, the highest extent of significant LD was observed in the wheat D-genome, followed by lower LD in the A- and B-genomes. The differences in the extent of LD among populations and genomes were mostly driven by differences in long-range LD ( > 10 cM).ConclusionsGenome- and population-specific patterns of genetic differentiation and LD were discovered in the populations of wheat cultivars from different geographic regions. Our study demonstrated that the estimates of population structure between spring and winter wheat lines can identify genomic regions harboring candidate genes involved in the regulation of growth habit. Variation in LD suggests that breeding and selection had a different impact on each wheat genome both within and among populations. The higher extent of LD in the wheat D-genome versus the A- and B-genomes likely reflects the episodes of recent introgression and population bottleneck accompanying the origin of hexaploid wheat. The assessment of LD and population structure in this assembled panel of diverse lines provides critical information for the development of genetic resources for genome-wide association mapping of agronomically important traits in wheat.

Identifying candidate causal variants via trans‐population fine‐mapping

Genome-wide association studies have discovered and confirmed a large number of loci that are implicated with disease susceptibility and severity. Polymorphisms that emerged from these studies are mostly indirectly associated to the phenotype, and the natural progression is to identify the causal variants that are functionally responsible for these association signals. Long stretches of high linkage disequilibrium (LD) benefitted the initial discovery phase in a genome-wide scan, allowing commercial genotyping products with imperfect coverage to detect genomic regions genuinely associated with the phenotype. However, regions of high LD confound the fine-mapping phase, as markers that are perfectly correlated to the causal variants display similar evidence of phenotypic association, hampering the process of differentiating the functional polymorphisms from neighboring surrogates. Here, we explore the potential of integrating information across different populations for narrowing the candidate region that a causal variant resides in, and compare the efficacy of this process of trans-population fine-mapping with the extent of variation in patterns of LD between the populations. In addition, we explore two different strategies for pooling data across multiple populations for the purpose of prioritizing the rankings of the causal variants. Our results clearly establish the benefits of trans-population analysis in reducing the number of possible candidates for the causal variants, particularly in genomic regions displaying strong evidence of inter-population LD variation. Directly integrating the statistical evidence by summing the test statistics outperforms the standard meta-analytic procedure. These findings have direct relevance to the design and analysis of ongoing fine-mapping studies.