Evidence for the role of EPHX2 gene variants in anorexia nervosa

Fast individual ancestry inference from DNA sequence data leveraging allele frequencies for multiple populations

Background

Estimation of individual ancestry from genetic data is useful for the analysis of disease association studies, understanding human population history and interpreting personal genomic variation. New, computationally efficient methods are needed for ancestry inference that can effectively utilize existing information about allele frequencies associated with different human populations and can work directly with DNA sequence reads.

Results

We describe a fast method for estimating the relative contribution of known reference populations to an individual’s genetic ancestry. Our method utilizes allele frequencies from the reference populations and individual genotype or sequence data to obtain a maximum likelihood estimate of the global admixture proportions using the BFGS optimization algorithm. It accounts for the uncertainty in genotypes present in sequence data by using genotype likelihoods and does not require individual genotype data from external reference panels. Simulation studies and application of the method to real datasets demonstrate that our method is significantly times faster than previous methods and has comparable accuracy. Using data from the 1000 Genomes project, we show that estimates of the genome-wide average ancestry for admixed individuals are consistent between exome sequence data and whole-genome low-coverage sequence data. Finally, we demonstrate that our method can be used to estimate admixture proportions using pooled sequence data making it a valuable tool for controlling for population stratification in sequencing based association studies that utilize DNA pooling.

Conclusions

Our method is an efficient and versatile tool for estimating ancestry from DNA sequence data and is available from https://sites.google.com/site/vibansal/software/iAdmix.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0418-7) contains supplementary material, which is available to authorized users.

Allele frequencies of variants in ultra conserved elements identify selective pressure on transcription factor binding

Ultra-conserved genes or elements (UCGs/UCEs) in the human genome are extreme examples of conservation. We characterized natural variations in 2884 UCEs and UCGs in two distinct populations; Singaporean Chinese (n = 280) and Italian (n = 501) by using a pooled sample, targeted capture, sequencing approach. We identify, with high confidence, in these regions the abundance of rare SNVs (MAF<0.5%) of which 75% is not present in dbSNP137. UCEs association studies for complex human traits can use this information to model expected background variation and thus necessary power for association studies. By combining our data with 1000 Genome Project data, we show in three independent datasets that prevalent UCE variants (MAF>5%) are more often found in relatively less-conserved nucleotides within UCEs, compared to rare variants. Moreover, prevalent variants are less likely to overlap transcription factor binding site. Using SNPfold we found no significant influence of RNA secondary structure on UCE conservation. All together, these results suggest UCEs are not under selective pressure as a stretch of DNA but are under differential evolutionary pressure on the single nucleotide level.

A genome-wide association study of anorexia nervosa

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

The role of leptin, melanocortin, and neurotrophin system genes on body weight in anorexia nervosa and bulimia nervosa

OBJECTIVE

Although low weight is a key factor contributing to the high mortality in anorexia nervosa (AN), it is unclear how AN patients sustain low weight compared with bulimia nervosa (BN) patients with similar psychopathology. Studies of genes involved in appetite and weight regulation in eating disorders have yielded variable findings, in part due to small sample size and clinical heterogeneity. This study: (1) assessed the role of leptin, melanocortin, and neurotrophin genetic variants in conferring risk for AN and BN; and (2) explored the involvement of these genes in body mass index (BMI) variations within AN and BN.

METHOD

Our sample consisted of 745 individuals with AN without a history of BN, 245 individuals with BN without a history of AN, and 321 controls. We genotyped 20 markers with known or putative function among genes selected from leptin, melanocortin, and neurotrophin systems.

RESULTS

There were no significant differences in allele frequencies among individuals with AN, BN, and controls. AGRP rs13338499 polymorphism was associated with lowest illness-related BMI in those with AN (p = 0.0013), and NTRK2 rs1042571 was associated with highest BMI in those with BN (p = 0.0018).

DISCUSSION

To our knowledge, this is the first study to address the issue of clinical heterogeneity in eating disorder genetic research and to explore the role of known or putatively functional markers in genes regulating appetite and weight in individuals with AN and BN. If replicated, our results may serve as an important first step toward gaining a better understanding of weight regulation in eating disorders.

Effect of soluble epoxide hydrolase polymorphism on substrate and inhibitor selectivity and dimer formation

Epoxy FAs (EpFAs) are important lipid mediators that are mainly metabolized by soluble epoxide hydrolase (sEH). Thus, sEH inhibition is a promising therapeutic target to treat numerous ailments. Several sEH polymorphisms result in amino acid substitutions and alter enzyme activity. K55R and R287Q are associated with inflammatory, cardiovascular, and metabolic diseases. R287Q seems to affect sEH activity through reducing formation of a catalytically active dimer. Thus, understanding how these SNPs affect the selectivity of sEH for substrates and inhibitors is of potential clinical importance. We investigated the selectivity of four sEH SNPs toward a series of EpFAs and inhibitors. We found that the SNPs alter the catalytic activity of the enzyme but do not alter the relative substrate and inhibitor selectivity. We also determined their dimer/monomer constants (KD/M). The WT sEH formed a very tight dimer, with a KD/M in the low picomolar range. Only R287Q resulted in a large change of the KD/M However, human tissue concentrations of sEH suggest that it is always in its dimer form independently of the SNP. These results suggest that the different biologies associated with K55R and R287Q are not explained by alteration in dimer formation or substrate selectivity.

Association claims in the sequencing era

Since the completion of the Human Genome Project, the field of human genetics has been in great flux, largely due to technological advances in studying DNA sequence variation. Although community-wide adoption of statistical standards was key to the success of genome-wide association studies, similar standards have not yet been globally applied to the processing and interpretation of sequencing data. It has proven particularly challenging to pinpoint unequivocally disease variants in sequencing studies of polygenic traits. Here, we comment on a number of factors that may contribute to irreproducible claims of association in scientific literature and discuss possible steps that we can take towards cultural change.

A genome-wide association study on common SNPs and rare CNVs in anorexia nervosa

Anorexia nervosa (AN) is a mental illness with high mortality that most commonly afflicts adolescent female individuals. Clinical symptoms include chronic food refusal, weight loss and body image distortions. We carried out a genome-wide association study on 1033 AN cases and 3733 pediatric control subjects, all of whom were of European ancestry and were genotyped on the Illumina HumanHap610 platform (Illumina, San Diego, CA, USA). We confirmed that common single-nucleotide polymorphisms (SNPs) within OPRD1 (rs533123, P=0.0015) confer risk for AN, and obtained suggestive evidence that common SNPs near HTR1D (rs7532266, P=0.04) confer risk for restricting-type AN specifically. However, no SNPs reached genome-wide significance in our data, whereas top association signals were detected near ZNF804B, CSRP2BP, NTNG1, AKAP6 and CDH9. In parallel, we performed genome-wide analysis on copy number variations (CNVs) using the signal intensity data from the SNP arrays. We did not find evidence that AN cases have more CNVs than control subjects, nor do they have over-representation of rare or large CNVs. However, we identified several regions with rare CNVs that were only observed in AN cases, including a recurrent 13q12 deletion (1.5 Mb) disrupting SCAS in two cases, and CNVs disrupting the CNTN6/CNTN4 region in several AN cases. In conclusion, our study suggests that both common SNPs and rare CNVs may confer genetic risk to AN. These results point to intriguing genes that await further validation in independent cohorts for confirmatory roles in AN.