Worldwide genetic variation in dopamine and serotonin pathway genes: implications for association studies

Genome-wide association meta-analysis highlights light-induced signaling as a driver for refractive error

Refractive errors, including myopia, are the most frequent eye disorders worldwide and an increasingly common cause of blindness. This genome-wide association meta-analysis in 160,420 participants and replication in 95,505 participants increased the number of established independent signals from 37 to 161 and showed high genetic correlation between Europeans and Asians (>0.78). Expression experiments and comprehensive in silico analyses identified retinal cell physiology and light processing as prominent mechanisms, and also identified functional contributions to refractive-error development in all cell types of the neurosensory retina, retinal pigment epithelium, vascular endothelium and extracellular matrix. Newly identified genes implicate novel mechanisms such as rod-and-cone bipolar synaptic neurotransmission, anterior-segment morphology and angiogenesis. Thirty-one loci resided in or near regions transcribing small RNAs, thus suggesting a role for post-transcriptional regulation. Our results support the notion that refractive errors are caused by a light-dependent retina-to-sclera signaling cascade and delineate potential pathobiological molecular drivers.

Candidate gene study reveals DRD1 and DRD2 as putative interacting risk factors for youth depression

Alterations in the monoaminergic neurotransmission systems are suspected to be involved in the etiology of neuropsychiatric disorders, including depression. The role of these pathways in the risk of developing depressive symptoms during childhood or adolescence is still not completely clear. This study sought to identify putative genetic factors in genes of serotonergic and dopaminergic systems modulating the level of manifestation of depressive symptoms in children and adolescents. We analyzed 170 single nucleotide polymorphisms (SNPs) in 21 candidate dopaminergic and serotonergic genes in a non-clinical sample of 410 Costa Rican participants of ages between 7 and 18 years, assessing the severity of depressive symptoms through the Child Depression Inventory (CDI). Genotypic and haplotypic associations, as well as epistatic effects, were examined. A significant interaction effect was detected between rs1039089 in conjunction with rs877138 located upstream of the dopamine D1 receptor (DRD1) and the dopamine D2 receptor (DRD2) genes respectively, although no evidence was found for any single variant or haplotype related to a differential liability. This newly described genetic interaction among putative regulatory regions of dopamine receptors could affect the level of manifestation of depressive symptoms through an imbalance of D1-D2 heteromers and modulation of cognitive processes.

MicroRNA Genetic Variation: From Population Analysis to Functional Implications of Three Allele Variants Associated with Cancer

Nucleotide variants in microRNA regions have been associated with disease; nevertheless, few studies still have addressed the allele-dependent effect of these changes. We studied microRNA genetic variation in human populations and found that while low-frequency variants accumulate indistinctly in microRNA regions, the mature and seed regions tend to be depleted of high-frequency variants, probably as a result of purifying selection. Comparison of pairwise population fixation indexes among regions showed that the seed had higher population fixation indexes than the other regions, suggesting the existence of local adaptation in the seed region. We further performed functional studies of three microRNA variants associated with cancer (rs2910164:C > G in MIR146A, rs11614913:C > T in MIR196A2, and rs3746444:A > G in both MIR499A and MIR499B). We found differences in the expression between alleles and in the regulation of several genes involved in cancer, such as TP53, KIT, CDH1, CLH, and TERT, which may result in changes in regulatory networks related to tumorigenesis. Furthermore, luciferase-based assays showed that MIR499A could be regulating the cadherin CDH1 and the cell adhesion molecule CLH1 in an allele-dependent fashion. A better understanding of the effect of microRNA variants associated with disease could be key in our way to a more personalized medicine.

Beyond main effects of gene-sets: harsh parenting moderates the association between a dopamine gene-set and child externalizing behavior

BACKGROUND

In a longitudinal cohort study, we investigated the interplay of harsh parenting and genetic variation across a set of functionally related dopamine genes, in association with children's externalizing behavior. This is one of the first studies to employ gene-based and gene-set approaches in tests of Gene by Environment (G × E) effects on complex behavior. This approach can offer an important alternative or complement to candidate gene and genome-wide environmental interaction (GWEI) studies in the search for genetic variation underlying individual differences in behavior.

METHODS

Genetic variants in 12 autosomal dopaminergic genes were available in an ethnically homogenous part of a population-based cohort. Harsh parenting was assessed with maternal (n = 1881) and paternal (n = 1710) reports at age 3. Externalizing behavior was assessed with the Child Behavior Checklist (CBCL) at age 5 (71 ± 3.7 months). We conducted gene-set analyses of the association between variation in dopaminergic genes and externalizing behavior, stratified for harsh parenting.

RESULTS

The association was statistically significant or approached significance for children without harsh parenting experiences, but was absent in the group with harsh parenting. Similarly, significant associations between single genes and externalizing behavior were only found in the group without harsh parenting. Effect sizes in the groups with and without harsh parenting did not differ significantly. Gene-environment interaction tests were conducted for individual genetic variants, resulting in two significant interaction effects (rs1497023 and rs4922132) after correction for multiple testing.

CONCLUSION

Our findings are suggestive of G × E interplay, with associations between dopamine genes and externalizing behavior present in children without harsh parenting, but not in children with harsh parenting experiences. Harsh parenting may overrule the role of genetic factors in externalizing behavior. Gene-based and gene-set analyses offer promising new alternatives to analyses focusing on single candidate polymorphisms when examining the interplay between genetic and environmental factors.

The genetics of East African populations: a Nilo-Saharan component in the African genetic landscape

East Africa is a strategic region to study human genetic diversity due to the presence of ethnically, linguistically, and geographically diverse populations. Here, we provide new insight into the genetic history of populations living in the Sudanese region of East Africa by analysing nine ethnic groups belonging to three African linguistic families: Niger-Kordofanian, Nilo-Saharan and Afro-Asiatic. A total of 500 individuals were genotyped for 200,000 single-nucleotide polymorphisms. Principal component analysis, clustering analysis using ADMIXTURE, F_ST statistics, and the three-population test were used to investigate the underlying genetic structure and ancestry of the different ethno-linguistic groups. Our analyses revealed a genetic component for Sudanese Nilo-Saharan speaking groups (Darfurians and part of Nuba populations) related to Nilotes of South Sudan, but not to other Sudanese populations or other sub-Saharan populations. Populations inhabiting the North of the region showed close genetic affinities with North Africa, with a component that could be remnant of North Africans before the migrations of Arabs from Arabia. In addition, we found very low genetic distances between populations in genes important for anti-malarial and anti-bacterial host defence, suggesting similar selective pressures on these genes and stressing the importance of considering functional pathways to understand the evolutionary history of populations.

The dopamine receptor D4 gene and familial loading interact with perceived parenting in predicting externalizing behavior problems in early adolescence: the TRacking Adolescents' Individual Lives Survey (TRAILS)

Although externalizing behavior problems show in general a high stability over time, the course of externalizing behavior problems may vary from individual to individual. Our main goal was to investigate the predictive role of parenting on externalizing behavior problems. In addition, we investigated the potential moderating role of gender and genetic risk (operationalized as familial loading of externalizing behavior problems (FLE), and presence or absence of the dopamine receptor D4 (DRD4) 7-repeat and 4-repeat allele, respectively). Perceived parenting (rejection, emotional warmth, and overprotection) and FLE were assessed in a population-based sample of 1768 10- to 12-year-old adolescents. Externalizing behavior problems were assessed at the same age and 212 years later by parent report (CBCL) and self-report (YSR). DNA was extracted from blood samples. Parental emotional warmth predicted lower, and parental overprotection and rejection predicted higher levels of externalizing behavior problems. Whereas none of the parenting factors interacted with gender and the DRD4 7-repeat allele, we did find interaction effects with FLE and the DRD4 4-repeat allele. That is, the predictive effect of parental rejection was only observed in adolescents from low FLE families and the predictive effect of parental overprotection was stronger in adolescents not carrying the DRD4 4-repeat allele.

Positive selection in the chromosome 16 VKORC1 genomic region has contributed to the variability of anticoagulant response in humans

VKORC1 (vitamin K epoxide reductase complex subunit 1, 16p11.2) is the main genetic determinant of human response to oral anticoagulants of antivitamin K type (AVK). This gene was recently suggested to be a putative target of positive selection in East Asian populations. In this study, we genotyped the HGDP-CEPH Panel for six VKORC1 SNPs and downloaded chromosome 16 genotypes from the HGDP-CEPH database in order to characterize the geographic distribution of footprints of positive selection within and around this locus. A unique VKORC1 haplotype carrying the promoter mutation associated with AVK sensitivity showed especially high frequencies in all the 17 HGDP-CEPH East Asian population samples. VKORC1 and 24 neighboring genes were found to lie in a 505 kb region of strong linkage disequilibrium in these populations. Patterns of allele frequency differentiation and haplotype structure suggest that this genomic region has been submitted to a near complete selective sweep in all East Asian populations and only in this geographic area. The most extreme scores of the different selection tests are found within a smaller 45 kb region that contains VKORC1 and three other genes (BCKDK, MYST1 (KAT8), and PRSS8) with different functions. Because of the strong linkage disequilibrium, it is not possible to determine if VKORC1 or one of the three other genes is the target of this strong positive selection that could explain present-day differences among human populations in AVK dose requirement. Our results show that the extended region surrounding a presumable single target of positive selection should be analyzed for genetic variation in a wide range of genetically diverse populations in order to account for other neighboring and confounding selective events and the hitchhiking effect.

Mapping the epistatic network underlying murine reproductive fatpad variation

Genome-wide mapping analyses are now commonplace in many species and several networks of interacting loci have been reported. However, relatively few details regarding epistatic interactions and their contribution to complex trait variation in multicellular organisms are available and the identification of positional candidate loci for epistatic QTL (epiQTL) is hampered, especially in mammals, by the limited genetic resolution inherent in most study designs. Here we further investigate the genetic architecture of reproductive fatpad weight in mice using the F(10) generation of the LG,SM advanced intercross (AI) line. We apply multiple mapping techniques including a single-locus model, locus-specific composite interval mapping (CIM), and tests for multiple QTL per chromosome to the 12 chromosomes known to harbor single-locus QTL (slQTL) affecting obesity in this cross. We also perform a genome-wide scan for pairwise epistasis. Using this combination of approaches we detect 199 peaks spread over all 19 autosomes, which potentially contribute to trait variation including all eight original F(2) loci (Adip1-8), novel slQTL peaks on chromosomes 7 and 9, and several novel epistatic loci. Extensive epistasis is confirmed involving both slQTL confidence intervals (C.I.) as well as regions that show no significant additive or dominance effects. These results provide important new insights into mapping complex genetic architectures and the role of epistasis in complex trait variation.

AGT M235T genotype/anxiety interaction and gender in the HyperGEN study

Background

Both anxiety and elevated heart rate (HR) have been implicated in the development of hypertension. The HyperGen cohort, consisting of siblings with severe and mild hypertension, an age-matched random sample of persons from the same base populations, and unmedicated adult offspring of the hypertensive siblings (N = 1,002 men and 987 women), was analyzed for an association of the angiotenisinogen AGTM235T genotype (TT, MT, MM) with an endophenotype, heart rate (HR) in high and low anxious groups.

Methodology

The interaction of AGTM genotype with anxiety, which has been independently associated with hypertension, was investigated adjusting for age, hypertension status, smoking, alcohol consumption, beta blocker medication, body mass index, physical activity and hours of television viewing (sedentary life style).

Principal Findings

Although there was no main effect of genotype on HR in men or women, high anxious men with the TT genotype had high HR, whereas high anxious men with the MM genotype had low HR. In women, HR was inversely associated with anxiety but there was no interaction with genotype.

Conclusion/Significance

The results suggest that high anxiety in men with the TT genotype may increase risk for hypertension whereas the MM genotype may be protective in high anxious men. This type of gene x environment interaction may be one reason why genome wide association studies sometimes fail to replicate. The locus may be important only in combination with certain environmental factors.

From 'omics' to complex disease: a systems biology approach to gene-environment interactions in cancer

Background

Cancer is a complex disease that involves a sequence of gene-environment interactions in a progressive process that cannot occur without dysfunction in multiple systems, including DNA repair, apoptotic and immune functions. Epigenetic mechanisms, responding to numerous internal and external cues in a dynamic ongoing exchange, play a key role in mediating environmental influences on gene expression and tumor development.

Hypothesis

The hypothesis put forth in this paper addresses the limited success of treatment outcomes in clinical oncology. It states that improvement in treatment efficacy requires a new paradigm that focuses on reversing systemic dysfunction and tailoring treatments to specific stages in the process. It requires moving from a reductionist framework of seeking to destroy aberrant cells and pathways to a transdisciplinary systems biology approach aimed at reversing multiple levels of dysfunction.

Conclusion

Because there are many biological pathways and multiple epigenetic influences working simultaneously in the expression of cancer phenotypes, studying individual components in isolation does not allow an adequate understanding of phenotypic expression. A systems biology approach using new modeling techniques and nonlinear mathematics is needed to investigate gene-environment interactions and improve treatment efficacy. A broader array of study designs will also be required, including prospective molecular epidemiology, immune competent animal models and in vitro/in vivo translational research that more accurately reflects the complex process of tumor initiation and progression.

Variants of the serotonin transporter gene and NEO-PI-R Neuroticism: No association in the BLSA and SardiNIA samples

The polymorphism in the serotonin transporter gene promoter region (5-HTTLPR) is by far the most studied variant hypothesized to influence Neuroticism-related personality traits. The results of previous studies have been mixed and appear moderated by the personality questionnaire used. Studies that used the TCI to assess Harm Avoidance or the EPQ to assess Neuroticism have found no association with the 5-HTTLPR. However, studies that used the NEO-PI-R or related instruments (NEO-PI, NEO-FFI) to measure Neuroticism have found some evidence of association. This study examines the association of variants in the serotonin transporter gene in a sample from a genetically isolated population within Sardinia (Italy) that is several times larger than previous samples that used the NEO-PI-R (N = 3,913). The association was also tested in a sample (N = 548) from the Baltimore Longitudinal Study of Aging (BLSA), in which repeated NEO-PI-R assessments were obtained. In the SardiNIA sample, we found no significant association of the 5-HTTLPR genotypes with Neuroticism or its facets (Anxiety, Angry-Hostility, Depression, Self-Consciousness, Impulsiveness, and Vulnerability). In the BLSA sample, we found lower scores on Neuroticism traits for the heterozygous group, which is inconsistent with previous studies. We also examined eight SNPs in the SardiNIA (N = 3,972) and nine SNPs in the BLSA (N = 1,182) that map within or near the serotonin transporter gene (SLC6A4), and found no association. Along with other large studies that used different phenotypic measures and found no association, this study substantially increases the evidence against a link between 5-HTT variants and Neuroticism-related traits.

Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

BACKGROUND

It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8.

RESULTS

Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates.

CONCLUSION

The "isolated population" label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.