Population-specific Effect Research Articles

Integrating genome-wide association studies and population genomics analysis reveals the genetic architecture of growth and backfat traits in pigs.

Growth and fat deposition are complex traits, which can affect economical income in the pig industry. Due to the intensive artificial selection, a significant genetic improvement has been observed for growth and fat deposition in pigs. Here, we first investigated genomic-wide association studies (GWAS) and population genomics (e.g., selection signature) to explore the genetic basis of such complex traits in two Large White pig lines (n = 3,727) with the GeneSeek GGP Porcine HD array (n = 50,915 SNPs). Ten genetic variants were identified to be associated with growth and fatness traits in two Large White pig lines from different genetic backgrounds by performing both within-population GWAS and cross-population GWAS analyses. These ten significant loci represented eight candidate genes, i.e., NRG4, BATF3, IRS2, ANO1, ANO9, RNF152, KCNQ5, and EYA2. One of them, ANO1 gene was simultaneously identified for both two lines in BF100 trait. Compared to single-population GWAS, cross-population GWAS was less effective for identifying SNPs with population-specific effect, but more powerful for detecting SNPs with population-shared effects. We further detected genomic regions specifically selected in each of two populations, but did not observe a significant enrichment for the heritability of growth and backfat traits in such regions. In summary, the candidate genes will provide an insight into the understanding of the genetic architecture of growth-related traits and backfat thickness, and may have a potential use in the genomic breeding programs in pigs.

Developmental temperature affects thermal dependence of locomotor activity in Drosophila.

In their natural environments, animals have to cope with fluctuations in numerous abiotic and biotic factors, and phenotypic plasticity can facilitate survival under such variable conditions. However, organisms may differ substantially in the ability to adjust their phenotypes in response to external factors. Here, we investigated how developmental temperature affects the thermal performance curve for locomotor activity in adult fruit flies (Drosophila melanogaster). We examined the thermal dependence of spontaneous activity in individuals originating from two natural populations (from tropical (India) and temperate climate zone (Slovakia)) that developed at three different temperatures (19°C, 25°C, and 29°C). Firstly, we found that developmental temperature has a significant impact on overall activity - flies that developed at high temperature (29°C) were, on average, less active than individuals that developed at lower temperatures. Secondly, developmental acclimation had a population-specific effect on the thermal optimum for activity. Whereas the optimal temperature was not affected by thermal conditions experienced during development in flies from India, developmental temperature shifted thermal optimum in flies from Slovakia. Thirdly, high developmental temperature broadened performance breadth in flies from the Indian population but narrowed it in individuals from the Slovak population. Finally, we did not detect a consistent effect of acclimation temperature on circadian rhythms of spontaneous activity. Altogether, our results demonstrate that developmental temperature can alter different parameters (maximum performance, thermal optimum, performance breadth) of the thermal performance curve for spontaneous activity. Since adult fruit flies are highly vagile, this sensitivity of locomotion to developmental conditions may be an important factor affecting fitness in changing environments.

Mendelian randomization study reveals a population-specific putative causal effect of type 2 diabetes in risk of cataract.

The epidemiological association between type 2 diabetes and cataract has been well established. However, it remains unclear whether the two diseases share a genetic basis, and if so, whether this reflects a putative causal relationship. We used East Asian population-based genome-wide association studies (GWAS) summary statistics of type 2 diabetes (Ncase = 36 614, Ncontrol = 155 150) and cataract (Ncase = 24 622, Ncontrol = 187 831) to comprehensively investigate the shared genetics between the two diseases. We performed: (i) linkage disequilibrium score regression (LDSC) and heritability estimation from summary statistics (ρ-HESS) to estimate the genetic correlation and local genetic correlation pattern between type 2 diabetes and cataract; (ii) multiple Mendelian randomization (MR) analyses to infer the putative causality between type 2 diabetes and cataract; and (iii) summary-data-based Mendelian randomization (SMR) to identify candidate risk genes underling the putative causality. Moreover, to investigate the extent of the population-specific genetic effect size underlying the shared genetics between type 2 diabetes and cataract, we applied the same analytical pipeline to perform a comparative analysis on European population-based GWAS of type 2 diabetes (Ncase = 62 892, Ncontrol = 596 424) and cataract (Ncase = 5045, Ncontrol = 356 096). Using East Asian population-based GWAS summary data, we observed a strong genetic correlation [rg = 0.58, 95% confidence interval (CI) = 0.33, 0.83), P-value = 5.60 × 10-6] between type 2 diabetes and cataract. Both ρ-HESS and multiple MR methods consistently showed a putative causal effect of type 2 diabetes on cataract, with estimated liability-scale MR odds ratios (ORs) at around 1.10 (95% CI = 1.06, 1.17). In contrast, no evidence supports a causal effect of cataract on type 2 diabetes. SMR analysis identified two novel genes MIR4453HG (βSMR = -0.34, 95% CI = -0.46, -0.22, P-value = 6.41 × 10-8) and KCNK17 (βSMR = -0.07, 95% CI = -0.09, -0.05, P-value = 2.49 × 10-10), whose expression levels were likely involved in the putative causality of type 2 diabetes on cataract. On the contrary, our comparative analysis on European population provided universally weak evidence on the genetic correlation and causal relationship between the two diseases. Our results provided robust evidence supporting a putative causal effect of type 2 diabetes on the risk of cataract in East Asians, and revealed potential genetic heterogeneity in the shared genetics underlying type 2 diabetes and cataract between East Asians and Europeans. These findings posed new paths on guiding the prevention and early-stage diagnosis of cataract in type 2 diabetes patients.

Population-specific causal disease effect sizes in functionally important regions impacted by selection

Many diseases exhibit population-specific causal effect sizes with trans-ethnic genetic correlations significantly less than 1, limiting trans-ethnic polygenic risk prediction. We develop a new method, S-LDXR, for stratifying squared trans-ethnic genetic correlation across genomic annotations, and apply S-LDXR to genome-wide summary statistics for 31 diseases and complex traits in East Asians (average N = 90K) and Europeans (average N = 267K) with an average trans-ethnic genetic correlation of 0.85. We determine that squared trans-ethnic genetic correlation is 0.82× (s.e. 0.01) depleted in the top quintile of background selection statistic, implying more population-specific causal effect sizes. Accordingly, causal effect sizes are more population-specific in functionally important regions, including conserved and regulatory regions. In regions surrounding specifically expressed genes, causal effect sizes are most population-specific for skin and immune genes, and least population-specific for brain genes. Our results could potentially be explained by stronger gene-environment interaction at loci impacted by selection, particularly positive selection.

Portable polygenic risk

Human Genomics Most genome-wide association studies (GWASs) have been performed on individuals of European descent. This imbalance has led to questions regarding how useful the findings of such studies, such as a polygenic risk score (PRS), are for individuals from other populations. Amariuta et al. examined 707 cell type–specific epigenetic regulatory annotations that had been identified in 111 GWASs investigating polygenic diseases and traits at transcription factor–binding sites in individuals of European and East Asian origin. Their method allowed them to prioritize those variants that affect gene regulation and improve transancestry PRS. Although population-specific causal effect sizes could not be determined, the results of this study indicate that information regarding functional genetic variation associated with disease and gene regulation is transferable across human populations. Nat. Genet. 52 , 1346 (2020).

Confirmation of GLRA3 as a susceptibility locus for albuminuria in Finnish patients with type 1 diabetes

Urinary albumin excretion is an early sign of diabetic kidney disease, affecting every third individual with diabetes. Despite substantial estimated heritability, only variants in the GLRA3 gene have been genome-wide significantly associated (p-value < 5 × 10−8) with diabetic albuminuria, in Finnish individuals with type 1 diabetes; However, replication attempt in non-Finnish Europeans with type 1 diabetes showed nominally significant association in the opposite direction, suggesting a population-specific effect, but simultaneously leaving the finding controversial. In this study, the association between the common rs10011025 variant in the GLRA3 locus, and albuminuria, was confirmed in 1259 independent Finnish individuals with type 1 diabetes (p = 0.0013), and meta-analysis of all Finnish individuals yielded a genome-wide significant association. The association was particularly pronounced in subjects not reaching the treatment target for blood glucose levels (HbA1c > 7%; N = 2560, p = 1.7 × 10−9). Even though further studies are needed to pinpoint the causal variants, dissecting the association at the GLRA3 locus may uncover novel molecular mechanisms for diabetic albuminuria irrespective of population background.

Resminostat and sorafenib combination therapy for advanced hepatocellular carcinoma in patients previously untreated with systemic chemotherapy.

252 Background: Resminostat is an oral hydroxamate-type inhibitor of class I, IIB, and IV histone deacetylases. A European Phase II study of second-line combination therapy with resminostat and sorafenib for hepatocellular carcinoma (HCC) in patients (pts) revealed a promising improvement in overall survival (OS). Here we report the findings on safety and efficacy of an Asian Phase I/II study on first-line combination therapy with sorafenib and resminostat in HCC pts. Methods: Pts with advanced or metastatic HCC considered Child-Pugh A and ECOG 0/1 were enrolled in Japan and Korea. Sorafenib was administered at 400 mg (bid) in both Phase I and II. Resminostat was administered on days 1 to 5 every 14 days. In Phase I, the dose of resminostat was escalated from 400 mg/day (DL1) to 600 mg/day (DL2). In Phase II, pts were randomly assigned to sorafenib monotherapy or sorafenib/resminostat combination therapy at a ratio of 1:1. The primary endpoint was time to progression (TTP). Tumor response was assessed according to RECIST version 1.1 every 6 weeks. Results: A total of 9 pts were enrolled in Phase I (DL1, 3 pts; DL2, 6 pts). Higher incidences of G3-4 toxicities, including one DLT (G4 thrombocytopenia), were observed at DL2. Therefore, DL1 was determined as the recommended dose for Phase II. A total of 170 pts were enrolled in Phase II. The median TTP was 2.8 months in the combination and control arm, respectively (HR: 0.984). No significant difference was observed in the median OS. Retrospective analysis revealed favorable results for the combination option in certain subgroups: for example, HBV+ (TTP: HR, 0.630; OS: HR, 0.846); no prior therapy (TTP: HR, 0.629; OS: HR, 0.590); and platelet count &gt; = 151.000 (TTP: HR, 0.646; OS: HR, 0.509). Conclusions: Although the primary endpoint was not reached in this Phase II all-comer HCC study, the results of the subgroup analysis suggest a population-specific effect for the combination therapy, especially in one which is HBV+. This warrants the further development of this combination as first-line therapy in a well-defined subset of pts with advanced HCC. Clinical trial information: NCT02400788.

The population determines whether and how life-history traits vary between reproductive events in an insect with maternal care.

The last reproductive event of a female is often associated with major changes in terms of both maternal and offspring life-history traits. However, the nature of these changes and the importance of population-specific environmental constraints in shaping their expression are difficult to predict and, as a consequence, poorly understood. Here, we investigated whether and how life-history traits vary between reproductive events and whether this variation is population-dependent in the European earwig Forficula auricularia. In this insect species, females produce up to two clutches during their lifetime and express extensive forms of maternal care. We conducted a common garden experiment, in which we measured 11 life-history traits of the first and second clutches of 132 females originating from three populations. Our results showed that clutch size was higherand the level of care expressed towards juveniles lowerin second as compared to the first clutches in all three populations. In contrast, we found a population-specific effect on whether and how the reproductive event shaped juvenile quality and a trade-off between egg developmental time and female weight at hatching. Overall, these findings emphasise that the last reproductive event of a female entails both positive and negative effects on various life-history traits of the female herself and her clutch of juveniles. Moreover, our study stresses the importance of population idiosyncrasies on the expression and nature of such cohort-specific effects.

Genome-wide meta-analysis identifies multiple novel associations and ethnic heterogeneity of psoriasis susceptibility.

Psoriasis is a common inflammatory skin disease with complex genetics and different degrees of prevalence across ethnic populations. Here we present the largest trans-ethnic genome-wide meta-analysis (GWMA) of psoriasis in 15,369 cases and 19,517 controls of Caucasian and Chinese ancestries. We identify four novel associations at LOC144817, COG6, RUNX1 and TP63, as well as three novel secondary associations within IFIH1 and IL12B. Fine-mapping analysis of MHC region demonstrates an important role for all three HLA class I genes and a complex and heterogeneous pattern of HLA associations between Caucasian and Chinese populations. Further, trans-ethnic comparison suggests population-specific effect or allelic heterogeneity for 11 loci. These population-specific effects contribute significantly to the ethnic diversity of psoriasis prevalence. This study not only provides novel biological insights into the involvement of immune and keratinocyte development mechanism, but also demonstrates a complex and heterogeneous genetic architecture of psoriasis susceptibility across ethnic populations.

Single nucleotide polymorphisms in candidate genes and dengue severity in children: A case–control, functional and meta-analysis study

Dengue is an arthropod-borne emerging viral disease with high morbidity and mortality risk in tropical countries like Brazil. Clinical manifestations are vast, ranging from asymptomatic to most severe forms of dengue such as shock. Previous data have shown that host genetics play a role in disease susceptibility and severity. Herein, we have tested the association of single nucleotide polymorphisms (SNPs) at TNF, IL10, MIF, DCSIGN, CLEC5A, NOD2, CCR5 and MRC1 as candidate genes using a matched case–control study design including 88 severe children cases of dengue patients and 335 healthy unrelated subjects that was also separated in IgG+ and IgG− controls. We demonstrated that the TT genotype of CLEC5A SNP (rs1285933 C>T) is associated with dengue severity (OR=2.25; p=0.03) and that GG genotype of −336G>A DCSIGN (CD209) SNP is associated with protection to severe dengue (OR=0.12; p=0.04). Both comparisons were borderline significant when cases were compared with IgG+ controls subgroup. Nevertheless, genotype–phenotype correlation was also assessed using serum levels of TNF from infected patients at the onset of dengue fever, and CT/TT carriers in CLEC5A secreted higher levels of TNF than CC individuals in 5–7days of infection. No significant difference was observed in TNF levels between genotypes GG versus AG/AA at DCSIGN promoter. Next, we performed a meta-analysis retrieving results from the literature for −336G>A DCSIGN and −308G>A TNF SNPs demonstrating that the consensus estimates of these SNPs indicated no association with dengue severity (when compared to Dengue fever) in the overall analysis. But, a subgroup analysis in the −336G>A DCSIGN, the G allele was associated with severe dengue susceptibility in Asians (ORallele=2.77; p=0.0001; ORcarriers=2.99; p=0.0001) and protection in Brazilians (ORallele=0.66; p=0.013). In summary, our results suggest that genetic variations at CLEC5A increase the risk and regulate TNF secretion in dengue severity among Brazilians. Also, combined data of the literature suggest population-specific effect of the −336 DCSIGN SNP more prominent in Asians and in a different direction than Brazilians.

Role of sepiapterin reductase gene at the PARK3 locus in Parkinson's disease

Sepiapterin reductase (SPR) gene is an enzyme which catalyses the final step of tetrahydrobiopterin synthesis (BH4) and was implicated in Parkinson's disease (PD) pathogenesis as a candidate gene for PARK3 locus. A number of studies yielded association of the PARK3 locus with PD, and SPR knockout mice were shown to display parkinsonian features. To evaluate the role of SPR gene polymorphisms in diverse populations in PD, we performed collaborative analyses in the Genetic Epidemiology of Parkinson Disease (GEO-PD) Consortium. A total of 5 single nucleotide polymorphisms (3 in the promoter region and 2 in the 3' untranslated region [UTR]) were genotyped. Fixed as well as random effect models were used to provide summary risk estimates of SPR variants. A total of 19 sites provided data for 6547 cases and 9321 controls. Overall odds ratio estimates varied from 0.92 to 1.01. No overall association with the SPR gene using either fixed effect or random effect model was observed in the studied population. I(2) Metric varied from 0% to 36.2%. There was some evidence for an association for participants of North European/Scandinavian descent with the strongest signal for rs1876487 (odds ratio = 0.82; p value = 0.003). Interestingly, families which were used to map the PARK3 locus, have Scandinavian ancestry suggesting a founder effect. In conclusion, this large association study for the SPR gene revealed no association for PD worldwide. However, taking the initial mapping of the PARK3 into account, the role of a population-specific effect warrants consideration in future studies.

Is FCGR2A a susceptibility gene to systemic lupus erythematosus in Chinese?

Recent genome-wide association scans and replication studies reinforce that FCGR2A is a susceptibility gene in systemic lupus erythematosus (SLE) in Caucasians. However, previous case control studies denied such conclusions in Chinese people. Besides genetic heterogeneity among different ethnicities, copy number variation (CNV), non-homogenous phenotypes and insufficient power may be confounders. We performed a case control study with 1066 Chinese (589 SLE patients and 477 healthy controls) and a meta-analysis based on 2328 SLE patients and 2313 healthy controls. FCGR2A CNV and FCGR2A131H/R [rs1801274] were detected by TaqMan assays. No variation of copy numbers of FCGR2A gene was found in Chinese. A further case control study suggested a dose-response character for FCGR2A131H/R and it affected disease activity, severity and prognosis. Finally, meta-analysis indicated FCGR2A that was a susceptibility gene to SLE in Chinese with an odds ratio of 1.094 and population attributable risk proportion of (PARP) 0.031. By an integrative strategy, we validate that FCGR2A bears no population-specific CNV. FCGR2A131H/R contributes to SLE susceptibility in Chinese, and affects disease activity, severity and prognosis. The undetected association in Chinese derives from under-power rather than any methodological obstacle due to CNV or population-specific genetic effect.

The SLC2A9 nonsynonymous Arg265His variant and gout: evidence for a population-specific effect on severity

IntroductionThe C allele of the nonsynonymous Arg265His (rs3733591) variant of SLC2A9 confers risk for gout in Han Chinese, Solomon Island and Japanese samples, with a stronger role in tophaceous gout. There is no evidence for an association with gout in Caucasian populations. In the present study, we tested rs3733591 for association with gout in New Zealand (NZ) Māori, Pacific Island and Caucasian samples.MethodsRs3733591 was genotyped across gout patients (n = 229, 232 and 327 NZ Māori, Pacific Island and Caucasian samples, respectively) and non-gout controls (n = 343, 174 and 638 Māori, Pacific Island and Caucasian samples, respectively). Further Caucasian sample sets consisting of 67 cases and 4,712 controls as well as 153 cases and 6,969 controls were obtained from the Framingham Heart Study and the Atherosclerosis Risk in Communities study, respectively. The Polynesian samples were analyzed according to Eastern and Western Polynesian ancestry.ResultsNo evidence for risk conferred by the C allele of rs3733591 with gout was found in the sample sets of NZ Māori (odd ratio (OR) = 0.98, P = 0.86), Eastern Polynesians (OR = 0.99, P = 0.92), Western Polynesians (OR = 1.16, P = 0.36) or combined Caucasians (OR = 1.15, P = 0.13). The C allele was significantly overrepresented in Māori tophaceous cases compared to cases without tophi (OR = 2.21, P = 0.008), but not in the other ancestral groupings.ConclusionsNoting that our study's power was limited for detecting weak genetic effects, we were unable to replicate associations of rs3733591 with gout in Eastern Polynesian, Western Polynesian and Caucasian samples. However, consistent with a previous study of Han Chinese and Solomon Island populations, our data suggest that rs3733591 could be a marker of severe gout in some populations. Our results also suggest that the effect of this variant is population-specific, further confirming population heterogeneity regarding the association of SLC2A9 with gout.

Detecting interacting genetic loci with effects on quantitative traits where the nature and order of the interaction are unknown

Standard techniques for single marker quantitative trait mapping perform poorly in detecting complex interacting genetic influences. When a genetic marker interacts with other genetic markers and/or environmental factors to influence a quantitative trait, a sample of individuals will show different effects according to their exposure to other interacting factors. This paper presents a Bayesian mixture model, which effectively models heterogeneous genetic effects apparent at a single marker. We compute approximate Bayes factors which provide an efficient strategy for screening genetic markers (genome-wide) for evidence of a heterogeneous effect on a quantitative trait. We present a simulation study which demonstrates that the approximation is good and provide a real data example which identifies a population-specific genetic effect on gene expression in the HapMap CEU and YRI populations. We advocate the use of the model as a strategy for identifying candidate interacting markers without any knowledge of the nature or order of the interaction. The source of heterogeneity can be modeled as an extension.

No evidence of association of FLJ10986 and ITPR2 with ALS in a large German cohort

A recent genome-wide association study (GWAS) found significant association of six single nucleotide polymorphisms (SNPs) in the gene FLJ10986 with sporadic amyotrophic lateral sclerosis (SALS). Another independent GWAS reported significant association of one SNP in the gene inositol 1,4,5-triphosphate receptor 2 (ITPR2) with SALS. These studies provided conflicting results. We examined the six most significant SNPs in FLJ10986 and one SNP in ITPR2 in a large cohort consisting of 595 SALS cases and 681 controls ascertained from Germany. Our results did not provide evidence for the association of these SNPs with SALS, suggesting a possible population-specific effect for FLJ10986 and ITPR2 that do not modulate the risk for SALS in the German population.

Mutation screening and association study of the candidate prostate cancer susceptibility genesMSR1,PTEN, andKLF6

MSR1, PTEN, and KLF6 have been implicated as candidate susceptibility genes for prostate tumorigenesis. Three hundred Jewish prostate cancer patients were screened for alterations in these genes. MSR1 was conserved in all patients. PTEN screening revealed a novel missense mutation and a silent change. Five KLF6 alterations were detected in 17 patients, including Q160X, the only nonsense KLF6 germline mutation described to date in a cancer patient. The KLF6 IVS1-27G>A polymorphism, recently associated with prostate cancer risk, was detected in 11.9% of the patients and 17.3% of the controls (P = 0.043). IVS1-27A allele frequency was significantly lower in prostate cancer patients (P = 0.030), specifically in Ashkenazi patients (P = 0.047) compared to controls. We found no evidence that MSR1 and PTEN germline mutations are associated with prostate cancer risk in Jews. The negative association between KLF6 IVS1-27A and prostate cancer risk supports a population-specific effect of susceptibility alleles in prostate tumorigenesis.

A Population-Based Validation of the American Joint Committee on Cancer Melanoma Staging System

A major revision of the American Joint Committee on Cancer (AJCC) stages for melanoma was implemented in 2002 after its validation in multinational cohorts including patients from cancer centers and cooperative groups. This staging system has not been validated in a US population-based cohort. We used 41,417 patients with primary invasive cutaneous melanoma diagnosed between 1988 and 2001 from the National Cancer Institute's Surveillance, Epidemiology, and End Results (SEER) cancer registry to validate the revised AJCC staging system. Survival rates computed from stage-specific Kaplan-Meier curves (time to melanoma-specific death) were compared with the survival rates from 17,600 patients in the original AJCC validation study. In the SEER cohort, 65% of reported melanomas were < or = 1.00 mm in thickness and 8.7% were more than 4.00 mm compared with 39% and 10% in the AJCC cohort (P < .001), respectively. AJCC stages were able to discriminate among SEER patient groups with different prognosis. However, SEER survival rates were significantly higher than those in the AJCC study and notably so in patients with T1a lesions (< or = 1 mm without ulceration). This population-specific effect remained significant after controlling for lesion thickness in all substages except stage IIA. Although this national population-based study validates the most recent revision of AJCC stages for melanoma, it emphasizes that survival rates are population specific and found them to be generally higher for SEER compared with AJCC patients. Population-specific survival rates should be used in study designs and decisions about patient-specific interventions.