Bao-Zhu Yang

Brain-Derived Neurotrophic Factor-5-HTTLPR Gene Interactions and Environmental Modifiers of Depression in Children

BACKGROUND Child abuse and genotype interact to contribute to risk for depression in children. This study examined gene-by-gene and gene-by-environment interactions. METHODS The study included 196 children: 109 maltreated and 87 nonmaltreated comparison subjects. Measures of psychiatric symptomatology and social supports were obtained using standard research instruments, and serotonin transporter (5-HTTLPR) (locus SLC6A4) and brain-derived neurotrophic factor (BDNF) (variant val66met) genotypes were obtained from saliva-derived DNA specimens. Population structure was controlled by means of ancestral proportion scores computed based on genotypes of ancestry informative markers in the entire sample. RESULTS There was a significant three-way interaction between BDNF genotype, 5-HTTLPR, and maltreatment history in predicting depression. Children with the met allele of the BDNF gene and two short alleles of 5-HTTLPR had the highest depression scores, but the vulnerability associated with these two genotypes was only evident in the maltreated children. A significant four-way interaction also emerged, with social supports found to further moderate risk for depression. CONCLUSIONS To the best of our knowledge, this is the first investigation to demonstrate a gene-by-gene interaction conveying vulnerability to depression. The current data also show a protective effect of social supports in ameliorating genetic and environmental risk for psychopathology.

Multiple Independent Loci at Chromosome 15q25.1 Affect Smoking Quantity: a Meta-Analysis and Comparison with Lung Cancer and COPD

Recently, genetic association findings for nicotine dependence, smoking behavior, and smoking-related diseases converged to implicate the chromosome 15q25.1 region, which includes the CHRNA5-CHRNA3-CHRNB4 cholinergic nicotinic receptor subunit genes. In particular, association with the nonsynonymous CHRNA5 SNP rs16969968 and correlates has been replicated in several independent studies. Extensive genotyping of this region has suggested additional statistically distinct signals for nicotine dependence, tagged by rs578776 and rs588765. One goal of the Consortium for the Genetic Analysis of Smoking Phenotypes (CGASP) is to elucidate the associations among these markers and dichotomous smoking quantity (heavy versus light smoking), lung cancer, and chronic obstructive pulmonary disease (COPD). We performed a meta-analysis across 34 datasets of European-ancestry subjects, including 38,617 smokers who were assessed for cigarettes-per-day, 7,700 lung cancer cases and 5,914 lung-cancer-free controls (all smokers), and 2,614 COPD cases and 3,568 COPD-free controls (all smokers). We demonstrate statistically independent associations of rs16969968 and rs588765 with smoking (mutually adjusted p-values<10−35 and <10−8 respectively). Because the risk alleles at these loci are negatively correlated, their association with smoking is stronger in the joint model than when each SNP is analyzed alone. Rs578776 also demonstrates association with smoking after adjustment for rs16969968 (p<10−6). In models adjusting for cigarettes-per-day, we confirm the association between rs16969968 and lung cancer (p<10−20) and observe a nominally significant association with COPD (p = 0.01); the other loci are not significantly associated with either lung cancer or COPD after adjusting for rs16969968. This study provides strong evidence that multiple statistically distinct loci in this region affect smoking behavior. This study is also the first report of association between rs588765 (and correlates) and smoking that achieves genome-wide significance; these SNPs have previously been associated with mRNA levels of CHRNA5 in brain and lung tissue.

Brain derived neurotrophic factor (BDNF) gene variants and Alzheimer's disease, affective disorders, posttraumatic stress disorder, schizophrenia, and substance dependence

Genetic variation at the locus encoding the brain derived neurotrophic factor (BDNF) has been implicated in some neuropsychiatric disorders such as Alzheimers disease (AD), affective disorders (AFDs), schizophrenia, and substance dependence. We therefore performed a mutation scan of the BDNF gene to identify novel gene variants and examined the association between BDNF variants and several neuropsychiatric phenotypes in European American (EA) subjects and controls. Using denaturing high performance liquid chromatography (dHPLC), we identified a novel variant (G−712A) in the putative promoter region. This variant and two previously reported BDNF SNPs (C270T and Val66Met) were genotyped in 295 patients with AD, 108 with AFDs, 96 with posttraumatic stress disorder (PTSD), 84 with schizophrenia, 327 with alcohol and/or drug dependence, and 250 normal control subjects. No association was found between these three BDNF gene variants and AD, AFDs, PTSD, or schizophrenia. However, there was a nominally higher frequency of the G−712A G‐allele and the G/G genotype in subjects with substance dependence than in controls (Allele: χ2 = 4.080, df = 1, P = 0.043; Genotype: χ2 = 7.225, df = 2, P = 0.027). Although after correction for multiple testing, the findings are not considered significant (threshold P‐value was set at 0.020 by the program SNPSpD), logistic regression analyses confirmed the modest association between SNP G−712A and substance dependence, when the sex and age of subjects were taken into consideration. The negative results for AFDs, PTSD, and schizophrenia could be due to the low statistical power. Further study with larger samples is warranted.

Genetic and environmental predictors of early alcohol use.

BACKGROUND The goal of the current investigation was to examine genetic and environmental predictors of early alcohol use, a potent predictor of later alcohol dependence. METHODS This study represents an add-on project to an investigation examining the efficacy of an intervention for maltreated children entering out-of-home care. Predictors of early alcohol use include the following: maltreatment, family loading for alcohol or substance-use disorders, and serotonin transporter genotype (5-HTTLPR; locus SLC6A4). Participants included 127 subjects: 76 maltreated children and 51 demographically matched community controls. RESULTS At follow-up, 29% of the maltreated children reported alcohol use, a rate more than seven times the rate observed in controls. Maltreated children also drank alcohol, on average, more than 2 years earlier than controls (11.2 vs. 13.5 years). Early alcohol use was predicted by maltreatment, 5-HTTLPR, and a gene by environment interaction, with increased risk for early alcohol use associated with the s-allele. Psychopathology at baseline, severity of maltreatment, and poor mother-child relations also predicted early alcohol use. CONCLUSIONS Maltreated children are at high risk for psychiatric, alcohol, and substance abuse problems. Examination of genetic and environmental risk and protective factors can help identify those who are most vulnerable and help guide prevention and intervention efforts.

The OPRD1 and OPRK1 loci in alcohol or drug dependence: OPRD1 variation modulates substance dependence risk.

Eleven single-nucleotide polymorphisms (SNPs) spanning OPRD1 were examined in 1063 European Americans (EAs) (620 cases with substance dependence (SD), including 557 with alcohol dependence (AD), 225 with cocaine dependence (CD) and 111 with opioid dependence (OD), and 443 controls). Nominally significant associations (P<0.05) of five SNPs with SD were observed; only the association of the non-synonymous variant G80T with OD remained significant after correction for multiple testing using SNPSpD. Haplotype analyses with six tag SNPs indicated that a specific haplotype GCAACT, which harbors G80T G-allele and C921T C-allele, was significantly associated with AD (χ2=14.82, degrees of freedom (d.f.)=1, P<0.001), CD (χ2=9.19, d.f.=1, P=0.002) and OD (χ2=20.68, d.f.=1, P<0.001). Logistic regression analyses, with sex and age being considered, demonstrated that this haplotype had a risk effect on AD (P=0.03, β=1.86, odds ratio (OR)=6.43) and especially on OD (P<0.001, β=3.92, OR=50.57). Moreover, seven SNPs covering OPRK1 were examined in the majority of the above subjects (390 cases, including 327 AD, 177 CD and 97 OD subjects, and 358 controls). Although no significant differences in allele, genotype or haplotype frequency distributions were seen between cases and controls, a specific OPRK1 haplotype, GGCTTCT, was significantly associated with AD (χ2=8.12, d.f.=1, P=0.004). Logistic regression analyses also revealed its risk effect on AD (P=0.009, β=1.06, OR=2.90). Population stratification artifact was not observed in the sample. Taken together, our findings supported a positive association between OPRD1 variants and SD, and a positive haplotypic association between OPRK1 and AD in EAs.

Child abuse and epigenetic mechanisms of disease risk

BACKGROUND Child abuse is highly prevalent and associated with increased risk for a range of health problems, including cancer, cardiovascular disease, diabetes, psychiatric disorders, and other health problems. Little is currently known about the mechanism by which early adversity confers risk for health problems later in life. PURPOSE To determine if there are epigenetic differences associated with child maltreatment that may help explain association between adverse childhood experiences and later health problems. METHODS As part of a study examining genetic and environmental factors associated with depression, saliva DNA specimens were collected on 96 maltreated children removed from their parents due to abuse or neglect and 96 demographically matched control children between 2003 and 2010. In 2011, the Illumina 450K BeadChip was used on stored DNA specimens and analyzed to examine whole-genome methylation differences between maltreated and control children. RESULTS After controlling for multiple comparisons, maltreated and control children had significantly different methylation values at 2868 CpG sites (p<5.0 × 10(-7), all sites; average methylation difference per site=17%; range=1%-62%). The gene set contained numerous markers of diseases and biological processes related to the health problems associated with early childhood adversity. CONCLUSIONS Although replication is required, this study suggests that epigenetic mechanisms may be associated with risk for health problems later in life in maltreated children. This study lays the groundwork for future studies examining health and methylation measures to further characterize the role of epigenetic mechanisms in conferring risk for medical problems in individuals with histories of early adversity.

ADH4 Gene Variation is Associated with Alcohol Dependence and Drug Dependence in European Americans: Results from HWD Tests and Case–Control Association Studies

The alcohol dehydrogenase (ADH) family constitutes one of the key sets of enzymes responsible for the oxidation of alcohol. The ADH4 gene, an important member of this family, is a functional and positional candidate for alcohol dependence. The present study aimed to investigate the relationship between ADH4 gene variation and alcohol dependence and drug dependence in European-Americans (EAs) and African-Americans (AAs). Seven single nucleotide polymorphisms (SNPs) spanning the ADH4 gene were genotyped in 365 healthy controls (317 EAs and 48 AAs) and 561 subjects (400 EAs and 161 AAs) affected with alcohol dependence and/or drug dependence (436 with alcohol dependence; 356 with drug dependence). Hardy–Weinberg equilibrium (HWE) for the genotype frequency distributions of these markers was tested in all phenotype groups to evaluate association between ADH4 gene variation and phenotypes and to fine-map the disease risk locus. The allele, genotype, and haplotype frequency distributions of these markers were compared between cases and controls to confirm the associations. The genotype frequency distributions of ADH4 markers were in HWE in EA controls, but were in Hardy–Weinberg disequilibrium (HWD) (ie, deviation from HWE) in EA cases. Among all markers, SNP2 (rs1042363) at exon 9 or SNP6 (rs1800759) at the promoter showed the greatest degree of HWD, among patients with either alcohol dependence or drug dependence. Significant differences between EA cases and controls were seen for genotype (10−6<global p<0.044), but not any allele or haplotype, frequency distributions for all seven ADH4 markers. These findings suggest that ADH4 genotypes predispose to alcohol dependence and drug dependence in a recessive manner, a predisposition that is population specific. SNP2 or SNP6 was the marker genetically closest to the functional risk loci for both alcohol dependence and drug dependence.

Variant Callers for Next-Generation Sequencing Data: A Comparison Study

Next generation sequencing (NGS) has been leading the genetic study of human disease into an era of unprecedented productivity. Many bioinformatics pipelines have been developed to call variants from NGS data. The performance of these pipelines depends crucially on the variant caller used and on the calling strategies implemented. We studied the performance of four prevailing callers, SAMtools, GATK, glftools and Atlas2, using single-sample and multiple-sample variant-calling strategies. Using the same aligner, BWA, we built four single-sample and three multiple-sample calling pipelines and applied the pipelines to whole exome sequencing data taken from 20 individuals. We obtained genotypes generated by Illumina Infinium HumanExome v1.1 Beadchip for validation analysis and then used Sanger sequencing as a “gold-standard” method to resolve discrepancies for selected regions of high discordance. Finally, we compared the sensitivity of three of the single-sample calling pipelines using known simulated whole genome sequence data as a gold standard. Overall, for single-sample calling, the called variants were highly consistent across callers and the pairwise overlapping rate was about 0.9. Compared with other callers, GATK had the highest rediscovery rate (0.9969) and specificity (0.99996), and the Ti/Tv ratio out of GATK was closest to the expected value of 3.02. Multiple-sample calling increased the sensitivity. Results from the simulated data suggested that GATK outperformed SAMtools and glfSingle in sensitivity, especially for low coverage data. Further, for the selected discrepant regions evaluated by Sanger sequencing, variant genotypes called by exome sequencing versus the exome array were more accurate, although the average variant sensitivity and overall genotype consistency rate were as high as 95.87% and 99.82%, respectively. In conclusion, GATK showed several advantages over other variant callers for general purpose NGS analyses. The GATK pipelines we developed perform very well.

Genetic variants of Nogo-66 Receptor with possible association to schizophrenia block myelin inhibition of axon growth

In schizophrenia, genetic predisposition has been linked to chromosome 22q11 and myelin-specific genes are misexpressed in schizophrenia. Nogo-66 receptor 1 (NGR or RTN4R) has been considered to be a 22q11 candidate gene for schizophrenia susceptibility because it encodes an axonal protein that mediates myelin inhibition of axonal sprouting. Confirming previous studies, we found that variation at the NGR locus is associated with schizophrenia in a Caucasian case-control analysis, and this association is not attributed to population stratification. Within a limited set of schizophrenia-derived DNA samples, we identified several rare NGR nonconservative coding sequence variants. Neuronal cultures demonstrate that four different schizophrenia-derived NgR1 variants fail to transduce myelin signals into axon inhibition, and function as dominant negatives to disrupt endogenous NgR1. This provides the first evidence that certain disease-derived human NgR1 variants are dysfunctional proteins in vitro. Mice lacking NgR1 protein exhibit reduced working memory function, consistent with a potential endophenotype of schizophrenia. For a restricted subset of individuals diagnosed with schizophrenia, the expression of dysfunctional NGR variants may contribute to increased disease risk.

The Role and Challenges of Exome Sequencing in Studies of Human Diseases

Recent advances in next-generation sequencing technologies have transformed the genetics study of human diseases; this is an era of unprecedented productivity. Exome sequencing, the targeted sequencing of the protein-coding portion of the human genome, has been shown to be a powerful and cost-effective method for detection of disease variants underlying Mendelian disorders. Increasing effort has been made in the interest of the identification of rare variants associated with complex traits in sequencing studies. Here we provided an overview of the application fields for exome sequencing in human diseases. We describe a general framework of computation and bioinformatics for handling sequencing data. We then demonstrate data quality and agreement between exome sequencing and exome microarray (chip) genotypes using data collected on the same set of subjects in a genetic study of panic disorder. Our results show that, in sequencing data, the data quality was generally higher for variants within the exonic target regions, compared to that outside the target regions, due to the target enrichment. We also compared genotype concordance for variant calls obtained by exome sequencing vs. exome genotyping microarrays. The overall consistency rate was >99.83% and the heterozygous consistency rate was >97.55%. The two platforms share a large amount of agreement over low frequency variants in the exonic regions, while exome sequencing provides much more information on variants not included on exome genotyping microarrays. The results demonstrate that exome sequencing data are of high quality and can be used to investigate the role of rare coding variants in human diseases.