Multi-ethnic fine-mapping of 14 central adiposity loci

Human and rat skeletal muscle single-nuclei multi-omic integrative analyses nominate causal cell types, regulatory elements, and SNPs for complex traits

Skeletal muscle accounts for the largest proportion of human body mass, on average, and is a key tissue in complex diseases and mobility. It is composed of several different cell and muscle fiber types. Here, we optimize single-nucleus ATAC-seq (snATAC-seq) to map skeletal muscle cell-specific chromatin accessibility landscapes in frozen human and rat samples, and single-nucleus RNA-seq (snRNA-seq) to map cell-specific transcriptomes in human. We additionally perform multi-omics profiling (gene expression and chromatin accessibility) on human and rat muscle samples. We capture type I and type II muscle fiber signatures, which are generally missed by existing single-cell RNA-seq methods. We perform cross-modality and cross-species integrative analyses on 33,862 nuclei and identify seven cell types ranging in abundance from 59.6% to 1.0% of all nuclei. We introduce a regression-based approach to infer cell types by comparing transcription start site-distal ATAC-seq peaks to reference enhancer maps and show consistency with RNA-based marker gene cell type assignments. We find heterogeneity in enrichment of genetic variants linked to complex phenotypes from the UK Biobank and diabetes genome-wide association studies in cell-specific ATAC-seq peaks, with the most striking enrichment patterns in muscle mesenchymal stem cells (∼3.5% of nuclei). Finally, we overlay these chromatin accessibility maps on GWAS data to nominate causal cell types, SNPs, transcription factor motifs, and target genes for type 2 diabetes signals. These chromatin accessibility profiles for human and rat skeletal muscle cell types are a useful resource for nominating causal GWAS SNPs and cell types.

Trans-ethnic meta-analysis of rare variants in sequencing association studies

Trans-ethnic meta-analysis is a powerful tool for detecting novel loci in genetic association studies. However, in the presence of heterogeneity among different populations, existing gene-/region-based rare variants meta-analysis methods may be unsatisfactory because they do not consider genetic similarity or dissimilarity among different populations. In response, we propose a score test under the modified random effects model for gene-/region-based rare variants associations. We adapt the kernel regression framework to construct the model and incorporate genetic similarities across populations into modeling the heterogeneity structure of the genetic effect coefficients. We use a resampling-based copula method to approximate asymptotic distribution of the test statistic, enabling efficient estimation of p-values. Simulation studies show that our proposed method controls type I error rates and increases power over existing approaches in the presence of heterogeneity. We illustrate our method by analyzing T2D-GENES consortium exome sequence data to explore rare variant associations with several traits.

Genetics of Body Fat Distribution: Comparative Analyses in Populations with European, Asian and African Ancestries

Preferential fat accumulation in visceral vs. subcutaneous depots makes obese individuals more prone to metabolic complications. Body fat distribution (FD) is regulated by genetics. FD patterns vary across ethnic groups independent of obesity. Asians have more and Africans have less visceral fat compared with Europeans. Consequently, Asians tend to be more susceptible to type 2 diabetes even with lower BMIs when compared with Europeans. To date, genome-wide association studies (GWAS) have identified more than 460 loci related to FD traits. However, the majority of these data were generated in European populations. In this review, we aimed to summarize recent advances in FD genetics with a focus on comparisons between European and non-European populations (Asians and Africans). We therefore not only compared FD-related susceptibility loci identified in three ethnicities but also discussed whether known genetic variants might explain the FD pattern heterogeneity across different ancestries. Moreover, we describe several novel candidate genes potentially regulating FD, including NID2, HECTD4 and GNAS, identified in studies with Asian populations. It is of note that in agreement with current knowledge, most of the proposed FD candidate genes found in Asians belong to the group of developmental genes.

A lead candidate functional single nucleotide polymorphism within the WARS2 gene associated with waist-hip-ratio does not alter RNA stability

We have prioritised a single nucleotide polymorphism (SNP) rs2645294 as one candidate functional SNP in the TBX15-WARS2 waist-hip-ratio locus using posterior probability analysis. This SNP is located in the 3' untranslated region of the WARS2 (tryptophanyl tRNA synthetase 2, mitochondrial) gene with which it has an expression quantitative trait in subcutaneous white adipose tissue. We show that transcripts of the WARS2 gene in a human white adipose cell line, heterozygous for the rs2645294 SNP, showed allelic imbalance. We tested whether the rs2645294 SNP altered WARS2 RNA stability using three different methods: actinomycin-D inhibition and RNA decay, mature and nascent RNA analysis and luciferase reporter assays. We found no evidence of a difference in RNA stability between the rs2645294 alleles indicating that the allelic expression imbalance was likely due to transcriptional regulation.

Osteoporosis- and obesity-risk interrelationships: an epigenetic analysis of GWAS-derived SNPs at the developmental gene TBX15

A major challenge in translating findings from genome-wide association studies (GWAS) to biological mechanisms is pinpointing functional variants because only a very small percentage of variants associated with a given trait actually impact the trait. We used an extensive epigenetics, transcriptomics, and genetics analysis of the TBX15/WARS2 neighbourhood to prioritize this region's best-candidate causal variants for the genetic risk of osteoporosis (estimated bone density, eBMD) and obesity (waist-hip ratio or waist circumference adjusted for body mass index). TBX15 encodes a transcription factor that is important in bone development and adipose biology. Manual curation of 692 GWAS-derived variants gave eight strong candidates for causal SNPs that modulate TBX15 transcription in subcutaneous adipose tissue (SAT) or osteoblasts, which highly and specifically express this gene. None of these SNPs were prioritized by Bayesian fine-mapping. The eight regulatory causal SNPs were in enhancer or promoter chromatin seen preferentially in SAT or osteoblasts at TBX15 intron-1 or upstream. They overlap strongly predicted, allele-specific transcription factor binding sites. Our analysis suggests that these SNPs act independently of two missense SNPs in TBX15. Remarkably, five of the regulatory SNPs were associated with eBMD and obesity and had the same trait-increasing allele for both. We found that WARS2 obesity-related SNPs can be ascribed to high linkage disequilibrium with TBX15 intron-1 SNPs. Our findings from GWAS index, proxy, and imputed SNPs suggest that a few SNPs, including three in a 0.7-kb cluster, act as causal regulatory variants to fine-tune TBX15 expression and, thereby, affect both obesity and osteoporosis risk.

EDAR, LYPLAL1, PRDM16, PAX3, DKK1, TNFSF12, CACNA2D3, and SUPT3H gene variants influence facial morphology in a Eurasian population

In human society, the facial surface is visible and recognizable based on the facial shape variation which represents a set of highly polygenic and correlated complex traits. Understanding the genetic basis underlying facial shape traits has important implications in population genetics, developmental biology, and forensic science. A number of single nucleotide polymorphisms (SNPs) are associated with human facial shape variation, mostly in European populations. To bridge the gap between European and Asian populations in term of the genetic basis of facial shape variation, we examined the effect of these SNPs in a European-Asian admixed Eurasian population which included a total of 612 individuals. The coordinates of 17 facial landmarks were derived from high resolution 3dMD facial images, and 136 Euclidean distances between all pairs of landmarks were quantitatively derived. DNA samples were genotyped using the Illumina Infinium Global Screening Array and imputed using the 1000 Genomes reference panel. Genetic association between 125 previously reported facial shape-associated SNPs and 136 facial shape phenotypes was tested using linear regression. As a result, a total of eight SNPs from different loci demonstrated significant association with one or more facial shape traits after adjusting for multiple testing (significance threshold p < 1.28 × 10-3), together explaining up to 6.47% of sex-, age-, and BMI-adjusted facial phenotype variance. These included EDAR rs3827760, LYPLAL1 rs5781117, PRDM16 rs4648379, PAX3 rs7559271, DKK1 rs1194708, TNFSF12 rs80067372, CACNA2D3 rs56063440, and SUPT3H rs227833. Notably, the EDAR rs3827760 and LYPLAL1 rs5781117 SNPs displayed significant association with eight and seven facial phenotypes, respectively (2.39 × 10-5 < p < 1.28 × 10-3). The majority of these SNPs showed a distinct allele frequency between European and East Asian reference panels from the 1000 Genomes Project. These results showed the details of above eight genes influence facial shape variation in a Eurasian population.

Mutant Wars2 gene in spontaneously hypertensive rats impairs brown adipose tissue function and predisposes to visceral obesity

Brown adipose tissue (BAT) plays an important role in lipid and glucose metabolism in rodents and possibly also in humans. Identification of genes responsible for BAT function would shed light on underlying pathophysiological mechanisms of metabolic disturbances. Recent linkage analysis in the BXH/HXB recombinant inbred (RI) strains, derived from Brown Norway (BN) and spontaneously hypertensive rats (SHR), identified two closely linked quantitative trait loci (QTL) associated with glucose oxidation and glucose incorporation into BAT lipids in the vicinity of Wars2 (tryptophanyl tRNA synthetase 2 (mitochondrial)) gene on chromosome 2. The SHR harbors L53F WARS2 protein variant that was associated with reduced angiogenesis and Wars2 thus represents a prominent positional candidate gene. In the current study, we validated this candidate as a quantitative trait gene (QTG) using transgenic rescue experiment. SHR-Wars2 transgenic rats with wild type Wars2 gene when compared to SHR, showed more efficient mitochondrial proteosynthesis and increased mitochondrial respiration, which was associated with increased glucose oxidation and incorporation into BAT lipids, and with reduced weight of visceral fat. Correlation analyses in RI strains showed that increased activity of BAT was associated with amelioration of insulin resistance in muscle and white adipose tissue. In summary, these results demonstrate important role of Wars2 gene in regulating BAT function and consequently lipid and glucose metabolism.

Dopamine gene methylation patterns are associated with obesity markers and carbohydrate intake

INTRODUCTION

Dopamine (DA) is a neurotransmitter that regulates the rewarding and motivational processes underlying food intake and eating behaviors. This study hypothesized associations of DNA methylation signatures at genes modulating DA signaling with obesity features, metabolic profiles, and dietary intake.

METHODS

An adult population within the Methyl Epigenome Network Association project was included (n = 473). DNA methylation levels in white blood cells were measured by microarray (450K). Differentially methylated genes were mapped within the dopaminergic synapse pathway using the KEGG reference database (map04728). Subsequently, network enrichment analyses were run in the pathDIP portal. Associations of methylation patterns with anthropometric markers of general (BMI) and abdominal obesity (waist circumference), the blood metabolic profile, and daily dietary intakes were screened.

RESULTS

After applying a correction for multiple comparisons, 12 CpG sites were strongly associated (p < 0.0001) with BMI: cg03489495 (ITPR3), cg22851378 (PPP2R2D), cg04021127 (PPP2R2D), cg22441882 (SLC18A1), cg03045635 (DRD5), cg23341970 (ITPR2), cg13051970 (DDC), cg08943004 (SLC6A3), cg20557710 (CACNA1C), cg24085522 (GNAL), cg16846691 (ITPR2), and cg09691393 (SLC6A3). Moreover, average methylation levels of these genes differed according to the presence or absence of abdominal obesity. Pathway analyses revealed a statistically significant contribution of the aforementioned genes to dopaminergic synapse transmission (p = 4.78E-08). Furthermore, SLC18A1 and SLC6A3 gene methylation signatures correlated with total energy (p < 0.001) and carbohydrate (p < 0.001) intakes.

CONCLUSIONS

The results of this investigation reveal that methylation status on DA signaling genes may underlie epigenetic mechanisms contributing to carbohydrate and calorie consumption and fat deposition.

Genetic variants in SEC16B are associated with body composition in black South Africans

Objective

The latest genome-wide association studies of obesity-related traits have identified several genetic loci contributing to body composition (BC). These findings have not been robustly replicated in African populations, therefore, this study aimed to assess whether European BC-associated gene loci played a similar role in a South African black population.

Methods

A replication and fine-mapping study was performed in participants from the Birth to Twenty cohort (N = 1,926) using the Metabochip. Measurements included body mass index (BMI), waist and hip circumference, waist-to-hip ratio (WHR), total fat mass, total lean mass and percentage fat mass (PFM).

Results

SNPs in several gene loci, including SEC16B (P_adj <  9.48 × 10⁻⁷), NEGR1 (P_adj < 1.64 × 10⁻⁶), FTO (P_adj < 2.91 × 10⁻⁵), TMEM18 (P_adj < 2.27 × 10⁻⁵), and WARS2(P_adj < 3.25 × 10⁻⁵) were similarly associated (albeit not at array-wide signficance (P ≤ 6.7 × 10⁻⁷) with various phenotypes including fat mass, PFM, WHR linked to BC in this African cohort, however the associations were driven by different sentinel SNPs. More importantly, DXA-derived BC measures revealed stronger genetic associations than simple anthropometric measures. Association signals generated in this study were shared by European and African populations, as well as unique to this African cohort. Moreover, sophisticated estimates like DXA measures enabled an enhanced characterisation of genetic associations for BC traits.

Conclusion

Results from this study suggest that in-depth genomic studies in larger African cohorts may reveal novel SNPs for body composition and adiposity, which will provide greater insight into the aetiology of obesity.

A fine-mapping study of central obesity loci incorporating functional annotation and imputation

A recent genome-wide association study (GWAS) of central obesity identified 27 loci, from sex-combined analysis, associated with waist-to-hip ratio adjusted for body-mass index (WHRadjBMI) in European-ancestry individuals. Nevertheless, the identified variants may not be the biological causal ones due to the presence of linkage disequilibrium (LD). To better understand the mechanisms underlying the identified loci from the GWAS meta-analysis, we first imputed summary statistics at GWAS loci to increase genetic resolution, and then we applied a Bayesian statistical fine-mapping method through PAINTOR, incorporating LD structure and functional annotations to select and prioritize the most plausible causal variants across WHRadjBMI-associated regions. Using adipose tissue- and cell-specific annotations that showed significant associations with WHRadjBMI, we identified 33 single-nucleotide polymorphisms (SNPs) from 27 sex-combined fine-mapping loci with posterior probability of causality greater than 0.9. Six of the selected 33 SNPs belong to at least one of the top five identified annotations. SNPs rs1440372 (SMAD6) and rs12608504 (JUND) are particularly important since they not only have associated functional annotations but are also GWA hits in the original study. Incorporation of functional annotations helps identify additional plausible causal variants, such as rs2213731 (DNM3-PIGC) and rs4531856 (JUND), that did not reach genome-wide significance in GWAS. Our results provide promising candidates for future functional validation experiments.

Genome-wide association study in Japanese females identifies fifteen novel skin-related trait associations

Skin trait variation impacts quality-of-life, especially for females from the viewpoint of beauty. To investigate genetic variation related to these traits, we conducted a GWAS of various skin phenotypes in 11,311 Japanese women and identified associations for age-spots, freckles, double eyelids, straight/curly hair, eyebrow thickness, hairiness, and sweating. In silico annotation with RoadMap Epigenomics epigenetic state maps and colocalization analysis of GWAS and GTEx Project eQTL signals provided information about tissue specificity, candidate causal variants, and functional target genes. Novel signals for skin-spot traits neighboured AKAP1/MSI2 (rs17833789; P = 2.2 × 10-9), BNC2 (rs10810635; P = 2.1 × 10-22), HSPA12A (rs12259842; P = 7.1 × 10-11), PPARGC1B (rs251468; P = 1.3 × 10-21), and RAB11FIP2 (rs10444039; P = 5.6 × 10-21). HSPA12A SNPs were the only protein-coding gene eQTLs identified across skin-spot loci. Double edged eyelid analysis identified that a signal around EMX2 (rs12570134; P = 8.2 × 10-15) was also associated with expression of EMX2 and the antisense-RNA gene EMX2OS in brain putamen basal ganglia tissue. A known hair morphology signal in EDAR was associated with both eyebrow thickness (rs3827760; P = 1.7 × 10-9) and straight/curly hair (rs260643; P = 1.6 × 10-103). Excessive hairiness signals' top SNPs were also eQTLs for TBX15 (rs984225; P = 1.6 × 10-8), BCL2 (rs7226979; P = 7.3 × 10-11), and GCC2 and LIMS1 (rs6542772; P = 2.2 × 10-9). For excessive sweating, top variants in two signals in chr2:28.82-29.05 Mb (rs56089836; P = 1.7 × 10-11) were eQTLs for either PPP1CB or PLB1, while a top chr16:48.26-48.45 Mb locus SNP was a known ABCC11 missense variant (rs6500380; P = 6.8 × 10-10). In total, we identified twelve loci containing sixteen association signals, of which fifteen were novel. These findings will help dermatologic researchers better understand the genetic underpinnings of skin-related phenotypic variation in human populations.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Current Challenges and New Opportunities for Gene-Environment Interaction Studies of Complex Diseases

Recently, many new approaches, study designs, and statistical and analytical methods have emerged for studying gene-environment interactions (G×Es) in large-scale studies of human populations. There are opportunities in this field, particularly with respect to the incorporation of -omics and next-generation sequencing data and continual improvement in measures of environmental exposures implicated in complex disease outcomes. In a workshop called "Current Challenges and New Opportunities for Gene-Environment Interaction Studies of Complex Diseases," held October 17-18, 2014, by the National Institute of Environmental Health Sciences and the National Cancer Institute in conjunction with the annual American Society of Human Genetics meeting, participants explored new approaches and tools that have been developed in recent years for G×E discovery. This paper highlights current and critical issues and themes in G×E research that need additional consideration, including the improved data analytical methods, environmental exposure assessment, and incorporation of functional data and annotations.

Evidence for Association between SH2B1 Gene Variants and Glycated Hemoglobin in Nondiabetic European American Young Adults: The Add Health Study

Glycated hemoglobin (HbA1c) is used to classify glycaemia and type 2 diabetes (T2D). Body mass index (BMI) is a predictor of HbA1c levels and T2D. We tested 43 established BMI and obesity loci for association with HbA1c in a nationally representative multiethnic sample of young adults from the National Longitudinal Study of Adolescent to Adult Health [Add Health: age 24-34 years; n = 5641 European Americans (EA); 1740 African Americans (AA); 1444 Hispanic Americans (HA)] without T2D, using two levels of covariate adjustment (Model 1: age, sex, smoking, and geographic region; Model 2: Model 1 covariates plus BMI). Bonferroni adjustment was made for 43 SNPs and we considered P < 0.0011 statistically significant. Means (SD) for HbA1c were 5.4% (0.3) in EA, 5.7% (0.4) in AA, and 5.5% (0.3) in HA. We observed significant evidence for association with HbA1c for two variants near SH2B1 in EA (rs4788102, P = 2.2 × 10(-4) ; rs7359397, P = 9.8 × 10(-4) ) for Model 1. Both results were attenuated after adjustment for BMI (rs4788102, P = 1.7 × 10(-3) ; rs7359397, P = 4.6 × 10(-3) ). No variant reached Bonferroni-corrected significance in AA or HA. These results suggest that SH2B1 polymorphisms are associated with HbA1c, largely independent of BMI, in EA young adults.

Genetic effects influencing risk for major depressive disorder in China and Europe

Major depressive disorder (MDD) is a common, complex psychiatric disorder and a leading cause of disability worldwide. Despite twin studies indicating its modest heritability (~30-40%), extensive heterogeneity and a complex genetic architecture have complicated efforts to detect associated genetic risk variants. We combined single-nucleotide polymorphism (SNP) summary statistics from the CONVERGE and PGC studies of MDD, representing 10 502 Chinese (5282 cases and 5220 controls) and 18 663 European (9447 cases and 9215 controls) subjects. We determined the fraction of SNPs displaying consistent directions of effect, assessed the significance of polygenic risk scores and estimated the genetic correlation of MDD across ancestries. Subsequent trans-ancestry meta-analyses combined SNP-level evidence of association. Sign tests and polygenic score profiling weakly support an overlap of SNP effects between East Asian and European populations. We estimated the trans-ancestry genetic correlation of lifetime MDD as 0.33; female-only and recurrent MDD yielded estimates of 0.40 and 0.41, respectively. Common variants downstream of GPHN achieved genome-wide significance by Bayesian trans-ancestry meta-analysis (rs9323497; log10 Bayes Factor=8.08) but failed to replicate in an independent European sample (P=0.911). Gene-set enrichment analyses indicate enrichment of genes involved in neuronal development and axonal trafficking. We successfully demonstrate a partially shared polygenic basis of MDD in East Asian and European populations. Taken together, these findings support a complex etiology for MDD and possible population differences in predisposing genetic factors, with important implications for future genetic studies.

Archaic Adaptive Introgression in TBX15/WARS2

A recent study conducted the first genome-wide scan for selection in Inuit from Greenland using single nucleotide polymorphism chip data. Here, we report that selection in the region with the second most extreme signal of positive selection in Greenlandic Inuit favored a deeply divergent haplotype that is closely related to the sequence in the Denisovan genome, and was likely introgressed from an archaic population. The region contains two genes, WARS2 and TBX15, and has previously been associated with adipose tissue differentiation and body-fat distribution in humans. We show that the adaptively introgressed allele has been under selection in a much larger geographic region than just Greenland. Furthermore, it is associated with changes in expression of WARS2 and TBX15 in multiple tissues including the adrenal gland and subcutaneous adipose tissue, and with regional DNA methylation changes in TBX15.

Signatures of Archaic Adaptive Introgression in Present-Day Human Populations

Comparisons of DNA from archaic and modern humans show that these groups interbred, and in some cases received an evolutionary advantage from doing so. This process-adaptive introgression-may lead to a faster rate of adaptation than is predicted from models with mutation and selection alone. Within the last couple of years, a series of studies have identified regions of the genome that are likely examples of adaptive introgression. In many cases, once a region was ascertained as being introgressed, commonly used statistics based on both haplotype as well as allele frequency information were employed to test for positive selection. Introgression by itself, however, changes both the haplotype structure and the distribution of allele frequencies, thus confounding traditional tests for detecting positive selection. Therefore, patterns generated by introgression alone may lead to false inferences of positive selection. Here we explore models involving both introgression and positive selection to investigate the behavior of various statistics under adaptive introgression. In particular, we find that the number and allelic frequencies of sites that are uniquely shared between archaic humans and specific present-day populations are particularly useful for detecting adaptive introgression. We then examine the 1000 Genomes dataset to characterize the landscape of uniquely shared archaic alleles in human populations. Finally, we identify regions that were likely subject to adaptive introgression and discuss some of the most promising candidate genes located in these regions.

Dissecting the genetics of complex traits using summary association statistics

During the past decade, genome-wide association studies (GWAS) have been used to successfully identify tens of thousands of genetic variants associated with complex traits and diseases. These studies have produced extensive repositories of genetic variation and trait measurements across large numbers of individuals, providing tremendous opportunities for further analyses. However, privacy concerns and other logistical considerations often limit access to individual-level genetic data, motivating the development of methods that analyse summary association statistics. Here, we review recent progress on statistical methods that leverage summary association data to gain insights into the genetic basis of complex traits and diseases.

Generalization and fine mapping of European ancestry-based central adiposity variants in African ancestry populations

BACKGROUND/OBJECTIVES

Central adiposity measures such as waist circumference (WC) and waist-to-hip ratio (WHR) are associated with cardiometabolic disorders independently of body mass index (BMI) and are gaining clinically utility. Several studies report genetic variants associated with central adiposity, but most utilize only European ancestry populations. Understanding whether the genetic associations discovered among mainly European descendants are shared with African ancestry populations will help elucidate the biological underpinnings of abdominal fat deposition.

SUBJECTS/METHODS

To identify the underlying functional genetic determinants of body fat distribution, we conducted an array-wide association meta-analysis among persons of African ancestry across seven studies/consortia participating in the Population Architecture using Genomics and Epidemiology (PAGE) consortium. We used the Metabochip array, designed for fine-mapping cardiovascular-associated loci, to explore novel array-wide associations with WC and WHR among 15 945 African descendants using all and sex-stratified groups. We further interrogated 17 known WHR regions for African ancestry-specific variants.

RESULTS

Of the 17 WHR loci, eight single-nucleotide polymorphisms (SNPs) located in four loci were replicated in the sex-combined or sex-stratified meta-analyses. Two of these eight independently associated with WHR after conditioning on the known variant in European descendants (rs12096179 in TBX15-WARS2 and rs2059092 in ADAMTS9). In the fine-mapping assessment, the putative functional region was reduced across all four loci but to varying degrees (average 40% drop in number of putative SNPs and 20% drop in genomic region). Similar to previous studies, the significant SNPs in the female-stratified analysis were stronger than the significant SNPs from the sex-combined analysis. No novel associations were detected in the array-wide analyses.

CONCLUSIONS

Of 17 previously identified loci, four loci replicated in the African ancestry populations of this study. Utilizing different linkage disequilibrium patterns observed between European and African ancestries, we narrowed the suggestive region containing causative variants for all four loci.

Genome-wide Trans-ethnic Meta-analysis Identifies Seven Genetic Loci Influencing Erythrocyte Traits and a Role for RBPMS in Erythropoiesis

Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.

PPARG-LYPLAL1 Multi-Allelic Combination Associated with Obesity and Overweight in Mexican Adolescent Females

OBJECTIVE

We studied multi-loci variants to identify the contribution of six candidate genes (ADIPOQ, CDH13, LYPLAL1, MC4R, PPARG and PGC1A) in the development of obesity and overweight.

DESIGN

We genotyped 404 chromosomes with eleven SNPs in Mexican female adolescents, who were subdivided into two groups (obesity-overweight and normal-weight) using the World Health Organization parameters. Genomic (800 chromosomes) and ancestral (208 chromosomes) controls were included to reduce the population bias. Anthropometric measurements, biochemical parameters, and caloric intake were obtained only in the groups of Mexican female adolescents.

RESULTS

A positive genotype-phenotype association was found that involves the multi-allelic combination of three risk alleles (one in PPARG and two in LYPLAL1) with obesity and overweight (OR=3.1, P=.010). This combination also exhibited a significant association with waist circumference (P=.030) and triglycerides levels (P=.030). These associations were supported by a logistic regression analysis adjusted for several confounding variables.

CONCLUSIONS

Our data suggest the joint participation of PPARG-LYPLAL1 genes in metabolic disorders development. Hence, these genes could act as potential biomarkers in obesity and overweight. Our findings underscore the complexity of metabolic disorders and provide evidence about the importance of multi-loci analysis to study complex diseases.

Resolving the etiology of atopic disorders by using genetic analysis of racial ancestry

Atopic dermatitis (AD), food allergy, allergic rhinitis, and asthma are common atopic disorders of complex etiology. The frequently observed atopic march from early AD to asthma, allergic rhinitis, or both later in life and the extensive comorbidity of atopic disorders suggest common causal mechanisms in addition to distinct ones. Indeed, both disease-specific and shared genomic regions exist for atopic disorders. Their prevalence also varies among races; for example, AD and asthma have a higher prevalence in African Americans when compared with European Americans. Whether this disparity stems from true genetic or race-specific environmental risk factors or both is unknown. Thus far, the majority of the genetic studies on atopic diseases have used populations of European ancestry, limiting their generalizability. Large-cohort initiatives and new analytic methods, such as admixture mapping, are currently being used to address this knowledge gap. Here we discuss the unique and shared genetic risk factors for atopic disorders in the context of ancestry variations and the promise of high-throughput "-omics"-based systems biology approach in providing greater insight to deconstruct their genetic and nongenetic etiologies. Future research will also focus on deep phenotyping and genotyping of diverse racial ancestry, gene-environment, and gene-gene interactions.

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.

Genome-wide association study of colorectal cancer in Hispanics

This manuscript describes the first large-scale genome-wide association study of colorectal cancer in Hispanics and Latinos. Our results demonstrate the broad replication of known susceptibility regions and the importance of fine-mapping in ethnic minority populations.

Genome-wide association studies (GWAS) have identified 58 susceptibility alleles across 37 regions associated with the risk of colorectal cancer (CRC) with P < 5×10⁻⁸. Most studies have been conducted in non-Hispanic whites and East Asians; however, the generalizability of these findings and the potential for ethnic-specific risk variation in Hispanic and Latino (HL) individuals have been largely understudied. We describe the first GWAS of common genetic variation contributing to CRC risk in HL (1611 CRC cases and 4330 controls). We also examine known susceptibility alleles and implement imputation-based fine-mapping to identify potential ethnicity-specific association signals in known risk regions. We discovered 17 variants across 4 independent regions that merit further investigation due to suggestive CRC associations (P < 1×10⁻⁶) at 1p34.3 (rs7528276; Odds Ratio (OR) = 1.86 [95% confidence interval (CI): 1.47–2.36); P = 2.5×10⁻⁷], 2q23.3 (rs1367374; OR = 1.37 (95% CI: 1.21–1.55); P = 4.0×10⁻⁷), 14q24.2 (rs143046984; OR = 1.65 (95% CI: 1.36–2.01); P = 4.1×10⁻⁷) and 16q12.2 [rs142319636; OR = 1.69 (95% CI: 1.37–2.08); P=7.8×10⁻⁷]. Among the 57 previously published CRC susceptibility alleles with minor allele frequency ≥1%, 76.5% of SNPs had a consistent direction of effect and 19 (33.3%) were nominally statistically significant (P < 0.05). Further, rs185423955 and rs60892987 were identified as novel secondary susceptibility variants at 3q26.2 (P = 5.3×10^–5) and 11q12.2 (P = 6.8×10⁻⁵), respectively. Our findings demonstrate the importance of fine mapping in HL. These results are informative for variant prioritization in functional studies and future risk prediction modeling in minority populations.

Evaluation of a Two-Stage Approach in Trans-Ethnic Meta-Analysis in Genome-Wide Association Studies

Meta-analysis of genome-wide association studies (GWAS) has achieved great success in detecting loci underlying human diseases. Incorporating GWAS results from diverse ethnic populations for meta-analysis, however, remains challenging because of the possible heterogeneity across studies. Conventional fixed-effects (FE) or random-effects (RE) methods may not be most suitable to aggregate multiethnic GWAS results because of violation of the homogeneous effect assumption across studies (FE) or low power to detect signals (RE). Three recently proposed methods, modified RE (RE-HE) model, binary-effects (BE) model and a Bayesian approach (Meta-analysis of Transethnic Association [MANTRA]), show increased power over FE and RE methods while incorporating heterogeneity of effects when meta-analyzing trans-ethnic GWAS results. We propose a two-stage approach to account for heterogeneity in trans-ethnic meta-analysis in which we clustered studies with cohort-specific ancestry information prior to meta-analysis. We compare this to a no-prior-clustering (crude) approach, evaluating type I error and power of these two strategies, in an extensive simulation study to investigate whether the two-stage approach offers any improvements over the crude approach. We find that the two-stage approach and the crude approach for all five methods (FE, RE, RE-HE, BE, MANTRA) provide well-controlled type I error. However, the two-stage approach shows increased power for BE and RE-HE, and similar power for MANTRA and FE compared to their corresponding crude approach, especially when there is heterogeneity across the multiethnic GWAS results. These results suggest that prior clustering in the two-stage approach can be an effective and efficient intermediate step in meta-analysis to account for the multiethnic heterogeneity.

Trans-ethnic study design approaches for fine-mapping

Studies that traverse ancestrally diverse populations may increase power to detect novel loci and improve fine-mapping resolution of causal variants by leveraging linkage disequilibrium differences between ethnic groups. The inclusion of African ancestry samples may yield further improvements because of low linkage disequilibrium and high genetic heterogeneity. We investigate the fine-mapping resolution of trans-ethnic fixed-effects meta-analysis for five type II diabetes loci, under various settings of ancestral composition (European, East Asian, African), allelic heterogeneity, and causal variant minor allele frequency. In particular, three settings of ancestral composition were compared: (1) single ancestry (European), (2) moderate ancestral diversity (European and East Asian), and (3) high ancestral diversity (European, East Asian, and African). Our simulations suggest that the European/Asian and European ancestry-only meta-analyses consistently attain similar fine-mapping resolution. The inclusion of African ancestry samples in the meta-analysis leads to a marked improvement in fine-mapping resolution.

Frailty and geography: should these two factors be added to the ABCDE contemporary guide to diabetes therapy?

On the road towards personalized treatments for type 2 diabetes, we suggest here that two parameters could be added to the ABCDE algorithm, 'F' for frailty and 'G' for geography. Indeed, the progressive ageing of population is causing a simultaneous increase of frailty worldwide. The identification of the optimal therapeutic approach is often difficult in frail subjects because of the complexity of 'frailty syndrome'. Nevertheless, given the relevance of diabetes in the development and progression of frailty, a safe and effective cure of diabetes is extremely important to guarantee a good medical outcome. There are few data about diabetes treatment in this delicate category of patients, and the choice of the appropriate therapy mostly remains a challenge. Moreover, type 2 diabetes affects more than 382 million people of different countries, races and ethnicities. To face the lack of solid evidence-based medicine for the treatment of diabetes in different ethnic groups, it is extremely important to increase knowledge about the different pathophysiology of diabetes according to ethnicity. In this way, a tailored approach to treatment of various ethnic groups living in the same or different regions can eventually be developed. Copyright © 2016 John Wiley & Sons, Ltd.

Integration of multiethnic fine-mapping and genomic annotation to prioritize candidate functional SNPs at prostate cancer susceptibility regions

Interpretation of biological mechanisms underlying genetic risk associations for prostate cancer is complicated by the relatively large number of risk variants (n = 100) and the thousands of surrogate SNPs in linkage disequilibrium. Here, we combined three distinct approaches: multiethnic fine-mapping, putative functional annotation (based upon epigenetic data and genome-encoded features), and expression quantitative trait loci (eQTL) analyses, in an attempt to reduce this complexity. We examined 67 risk regions using genotyping and imputation-based fine-mapping in populations of European (cases/controls: 8600/6946), African (cases/controls: 5327/5136), Japanese (cases/controls: 2563/4391) and Latino (cases/controls: 1034/1046) ancestry. Markers at 55 regions passed a region-specific significance threshold (P-value cutoff range: 3.9 × 10(-4)-5.6 × 10(-3)) and in 30 regions we identified markers that were more significantly associated with risk than the previously reported variants in the multiethnic sample. Novel secondary signals (P < 5.0 × 10(-6)) were also detected in two regions (rs13062436/3q21 and rs17181170/3p12). Among 666 variants in the 55 regions with P-values within one order of magnitude of the most-associated marker, 193 variants (29%) in 48 regions overlapped with epigenetic or other putative functional marks. In 11 of the 55 regions, cis-eQTLs were detected with nearby genes. For 12 of the 55 regions (22%), the most significant region-specific, prostate-cancer associated variant represented the strongest candidate functional variant based on our annotations; the number of regions increased to 20 (36%) and 27 (49%) when examining the 2 and 3 most significantly associated variants in each region, respectively. These results have prioritized subsets of candidate variants for downstream functional evaluation.

Leveraging Functional-Annotation Data in Trans-ethnic Fine-Mapping Studies

Localization of causal variants underlying known risk loci is one of the main research challenges following genome-wide association studies. Risk loci are typically dissected through fine-mapping experiments in trans-ethnic cohorts for leveraging the variability in the local genetic structure across populations. More recent works have shown that genomic functional annotations (i.e., localization of tissue-specific regulatory marks) can be integrated for increasing fine-mapping performance within single-population studies. Here, we introduce methods that integrate the strength of association between genotype and phenotype, the variability in the genetic backgrounds across populations, and the genomic map of tissue-specific functional elements to increase trans-ethnic fine-mapping accuracy. Through extensive simulations and empirical data, we have demonstrated that our approach increases fine-mapping resolution over existing methods. We analyzed empirical data from a large-scale trans-ethnic rheumatoid arthritis (RA) study and showed that the functional genetic architecture of RA is consistent across European and Asian ancestries. In these data, we used our proposed methods to reduce the average size of the 90% credible set from 29 variants per locus for standard non-integrative approaches to 22 variants.

LIN28A Modulates Splicing and Gene Expression Programs in Breast Cancer Cells

LIN28 is an evolutionarily conserved RNA-binding protein with critical functions in developmental timing and cancer. However, the molecular mechanisms underlying LIN28's oncogenic properties are yet to be described. RNA-protein immunoprecipitation coupled with genome-wide sequencing (RIP-Seq) analysis revealed significant LIN28 binding within 843 mRNAs in breast cancer cells. Many of the LIN28-bound mRNAs are implicated in the regulation of RNA and cell metabolism. We identify heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1), a protein with multiple roles in mRNA metabolism, as a LIN28-interacting partner. Subsequently, we used a custom computational method to identify differentially spliced gene isoforms in LIN28 and hnRNP A1 small interfering RNA (siRNA)-treated cells. The results reveal that these proteins regulate alternative splicing and steady-state mRNA expression of genes implicated in aspects of breast cancer biology. Notably, cells lacking LIN28 undergo significant isoform switching of the ENAH gene, resulting in a decrease in the expression of the ENAH exon 11a isoform. The expression of ENAH isoform 11a has been shown to be elevated in breast cancers that express HER2. Intriguingly, analysis of publicly available array data from the Cancer Genome Atlas (TCGA) reveals that LIN28 expression in the HER2 subtype is significantly different from that in other breast cancer subtypes. Collectively, our data suggest that LIN28 may regulate splicing and gene expression programs that drive breast cancer subtype phenotypes.

Identifying functional noncoding variants from genome-wide association studies for cardiovascular disease and related traits

PURPOSE OF REVIEW

Genome-wide association studies have identified many novel loci for cardiovascular disease and related traits. Attention is now shifting towards the analysis of these loci for causal variants, with a view to identify the novel mechanisms leading to disease.

RECENT FINDINGS

This review focuses on the approaches to identify causal, noncoding variants for coronary artery disease, lipid traits and other cardiovascular risk factors. Fine-mapping studies are discussed, along with the novel statistical approaches to produce 'credible sets'. The use of combining genome-wide association study datasets with experimental methods such as expression quantitative trait loci and allele-specific chromatin accessibility are explored, with recent examples discussed. Mapping long-range chromatin interactions and evolving genome-editing technologies such as clustered regularly interspaced short palindromic repeats combined with clustered regularly interspaced short palindromic repeats-associated (Cas9) nuclease promise to aid considerably the search for causal variants.

SUMMARY

Identification of causal variants for cardiovascular disease and related traits is still in the early stages, but with technologies evolving and increasingly relevant tissue samples undergoing analysis, there are favourable prospects that many new mechanisms for disease will be uncovered by the end of this decade.

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.