High-Resolution Melting-Based Marker Development for Wheat Leaf Rust Resistance Gene Lr34.

Deploying adult plant resistance (APR) against rust diseases is an important breeding objective of most wheat-breeding programs. The gene Lr34 is an effective and widely deployed broad-spectrum APR gene in wheat against leaf rust fungus Puccinia triticina. Various molecular markers have been developed for Lr34, but they either require post-PCR handling processes or are not economical. Herein, we developed a high-resolution melting (HRM)-based diagnostic assay for Lr34 based on a 3-bp 'TTC' deletion in exon 11 of the resistant allele. The susceptible cultivar Thatcher (Tc) and the near-isogenic Thatcher line (RL6058) with Lr34 yielded distinct melting profiles and were differentiated with high reproducibility. For further validation, all three copies of Lr34 were cloned in plasmid vectors, and HRM analysis using individual and combination (equimolar mixture of three copies) homoeologs yielded distinct melting profiles. An additional layer of genotyping was provided by a LunaProbe assay. The allele-specific probes successfully distinguished the homoeologs but not Tc and RL6058. Furthermore, the practical deployment of the HRM assay was tested by running the marker on a set of breeding lines. When compared with a kompetitive allele-specific PCR (KASP) Lr34 assay, the HRM assay had similar genotyping results and was able to accurately differentiate the resistant and susceptible breeding lines. However, our HRM assay was unable to detect the heterozygote. To our knowledge, this is the first report of an HRM assay for genotyping a wheat rust resistance gene.

Pleiotropic Meta-Analysis of Cognition, Education, and Schizophrenia Differentiates Roles of Early Neurodevelopmental and Adult Synaptic Pathways.

Susceptibility to schizophrenia is inversely correlated with general cognitive ability at both the phenotypic and the genetic level. Paradoxically, a modest but consistent positive genetic correlation has been reported between schizophrenia and educational attainment, despite the strong positive genetic correlation between cognitive ability and educational attainment. Here we leverage published genome-wide association studies (GWASs) in cognitive ability, education, and schizophrenia to parse biological mechanisms underlying these results. Association analysis based on subsets (ASSET), a pleiotropic meta-analytic technique, allowed jointly associated loci to be identified and characterized. Specifically, we identified subsets of variants associated in the expected ("concordant") direction across all three phenotypes (i.e., greater risk for schizophrenia, lower cognitive ability, and lower educational attainment); these were contrasted with variants that demonstrated the counterintuitive ("discordant") relationship between education and schizophrenia (i.e., greater risk for schizophrenia and higher educational attainment). ASSET analysis revealed 235 independent loci associated with cognitive ability, education, and/or schizophrenia at p < 5 × 10-8. Pleiotropic analysis successfully identified more than 100 loci that were not significant in the input GWASs. Many of these have been validated by larger, more recent single-phenotype GWASs. Leveraging the joint genetic correlations of cognitive ability, education, and schizophrenia, we were able to dissociate two distinct biological mechanisms-early neurodevelopmental pathways that characterize concordant allelic variation and adulthood synaptic pruning pathways-that were linked to the paradoxical positive genetic association between education and schizophrenia. Furthermore, genetic correlation analyses revealed that these mechanisms contribute not only to the etiopathogenesis of schizophrenia but also to the broader biological dimensions implicated in both general health outcomes and psychiatric illness.

Age-induced and photoinduced changes in gene expression profiles in facial skin of Caucasian females across 6 decades of age.

BACKGROUND: Intrinsic and extrinsic factors, including ultraviolet irradiation, lead to visible signs of skin aging. OBJECTIVE: We evaluated molecular changes occurring in photoexposed and photoprotected skin of white women 20 to 74 years of age, some of whom appeared substantially younger than their chronologic age. METHODS: Histologic and transcriptomics profiling were conducted on skin biopsy samples of photoexposed (face and dorsal forearm) or photoprotected (buttocks) body sites from 158 women. 23andMe genotyping determined genetic ancestry. RESULTS: Gene expression and ontologic analysis revealed progressive changes from the 20s to the 70s in pathways related to oxidative stress, energy metabolism, senescence, and epidermal barrier; these changes were accelerated in the 60s and 70s. The gene expression patterns from the subset of women who were younger-appearing were similar to those in women who were actually younger. LIMITATIONS: Broader application of these findings (eg, across races and Fitzpatrick skin types) will require further studies. CONCLUSIONS: This study demonstrates a wide range of molecular processes in skin affected by aging, providing relevant targets for improving the condition of aging skin at different life stages and defining a molecular pattern of epidermal gene expression in women who appear younger than their chronologic age.

Multi-Trait Analysis of GWAS and Biological Insights Into Cognition: A Response to Hill (2018).

Hill (Twin Research and Human Genetics, Vol. 21, 2018, 84-88) presented a critique of our recently published paper in Cell Reports entitled 'Large-Scale Cognitive GWAS Meta-Analysis Reveals Tissue-Specific Neural Expression and Potential Nootropic Drug Targets' (Lam et al., Cell Reports, Vol. 21, 2017, 2597-2613). Specifically, Hill offered several interrelated comments suggesting potential problems with our use of a new analytic method called Multi-Trait Analysis of GWAS (MTAG) (Turley et al., Nature Genetics, Vol. 50, 2018, 229-237). In this brief article, we respond to each of these concerns. Using empirical data, we conclude that our MTAG results do not suffer from 'inflation in the FDR [false discovery rate]', as suggested by Hill (Twin Research and Human Genetics, Vol. 21, 2018, 84-88), and are not 'more relevant to the genetic contributions to education than they are to the genetic contributions to intelligence'.

A Common Polymorphism in SCN2A Predicts General Cognitive Ability through Effects on PFC Physiology.

Here we provide novel convergent evidence across three independent cohorts of healthy adults (n = 531), demonstrating that a common polymorphism in the gene encoding the α2 subunit of neuronal voltage-gated type II sodium channels (SCN2A) predicts human general cognitive ability or "g." Using meta-analysis, we demonstrate that the minor T allele of a common polymorphism (rs10174400) in SCN2A is associated with significantly higher "g" independent of gender and age. We further demonstrate using resting-state fMRI data from our discovery cohort (n = 236) that this genetic advantage may be mediated by increased capacity for information processing between the dorsolateral PFC and dorsal ACC, which support higher cognitive functions. Collectively, these findings fill a gap in our understanding of the genetics of general cognitive ability and highlight a specific neural mechanism through which a common polymorphism shapes interindividual variation in "g."

Identifying nootropic drug targets via large-scale cognitive GWAS and transcriptomics.

Broad-based cognitive deficits are an enduring and disabling symptom for many patients with severe mental illness, and these impairments are inadequately addressed by current medications. While novel drug targets for schizophrenia and depression have emerged from recent large-scale genome-wide association studies (GWAS) of these psychiatric disorders, GWAS of general cognitive ability can suggest potential targets for nootropic drug repurposing. Here, we (1) meta-analyze results from two recent cognitive GWAS to further enhance power for locus discovery; (2) employ several complementary transcriptomic methods to identify genes in these loci that are credibly associated with cognition; and (3) further annotate the resulting genes using multiple chemoinformatic databases to identify "druggable" targets. Using our meta-analytic data set (N = 373,617), we identified 241 independent cognition-associated loci (29 novel), and 76 genes were identified by 2 or more methods of gene identification. Actin and chromatin binding gene sets were identified as novel pathways that could be targeted via drug repurposing. Leveraging our transcriptomic and chemoinformatic databases, we identified 16 putative genes targeted by existing drugs potentially available for cognitive repurposing.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes.

BACKGROUND: A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat. RESULTS: Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed. CONCLUSIONS: In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Improving Prediction Accuracy Using Multi-allelic Haplotype Prediction and Training Population Optimization in Wheat.

The use of haplotypes may improve the accuracy of genomic prediction over single SNPs because haplotypes can better capture linkage disequilibrium and genomic similarity in different lines and may capture local high-order allelic interactions. Additionally, prediction accuracy could be improved by portraying population structure in the calibration set. A set of 383 advanced lines and cultivars that represent the diversity of the University of Minnesota wheat breeding program was phenotyped for yield, test weight, and protein content and genotyped using the Illumina 90K SNP Assay. Population structure was confirmed using single SNPs. Haplotype blocks of 5, 10, 15, and 20 adjacent markers were constructed for all chromosomes. A multi-allelic haplotype prediction algorithm was implemented and compared with single SNPs using both k-fold cross validation and stratified sampling optimization. After confirming population structure, the stratified sampling improved the predictive ability compared with k-fold cross validation for yield and protein content, but reduced the predictive ability for test weight. In all cases, haplotype predictions outperformed single SNPs. Haplotypes of 15 adjacent markers showed the best improvement in accuracy for all traits; however, this was more pronounced in yield and protein content. The combined use of haplotypes of 15 adjacent markers and training population optimization significantly improved the predictive ability for yield and protein content by 14.3 (four percentage points) and 16.8% (seven percentage points), respectively, compared with using single SNPs and k-fold cross validation. These results emphasize the effectiveness of using haplotypes in genomic selection to increase genetic gain in self-fertilized crops.

Genetic architecture of agronomic and quality traits in a nested association mapping population of spring wheat.

Germplasm collections are rich sources of genetic variation to improve crops for many valuable traits. Nested association mapping (NAM) populations can overcome the limitations of genome-wide association studies (GWAS) in germplasm collections by reducing the effect of population structure. We exploited the genetic diversity of the USDA-ARS wheat (Triticum aestivum L.) core collection by developing the Spring Wheat Multiparent Introgression Population (SWMIP). To develop this population, twenty-five core parents were crossed and backcrossed to the Minnesota spring wheat cultivar RB07. The NAM population and 26 founder parents were genotyped using genotyping-by-sequencing and phenotyped for heading date, height, test weight, and grain protein content. After quality control, 20,312 markers with physical map positions were generated for 2,038 recombinant inbred lines (RILs). The number of RILs in each family varied between 58 and 96. Three GWAS models were utilized for quantitative trait loci (QTL) detection and accounted for known family stratification, genetic kinship, and both covariates. GWAS was performed on the whole population and also by bootstrap sampling of an equal number of RILs from each family. Greater power of QTL detection was achieved by treating families equally through bootstrapping. In total 16, 15, 12, and 13 marker-trait associations (MTAs) were identified for heading date, height, test weight, and grain protein content, respectively. Some of these MTAs were coincident with major genes known to control the traits, but others were novel and contributed by the wheat core parents. The SWMIP will be a valuable source of genetic variation for spring wheat breeding.

Bisphosphonates and bone fractures in long-term kidney transplant recipients.

BACKGROUND: There is little information on the role of bisphosphonates and bone mineral density (BMD) measurements for the follow-up and management of bone loss and fractures in long-term kidney transplant recipients. METHODS: To address this question, we retrospectively studied 554 patients who had two BMD measurements after the first year posttransplant and compared outcomes in patients treated, or not with bisphosphonates between the two BMD assessments. Kaplan-Meier survival and stepwise Cox regression analyses were performed to examine fracture-free survival rates and the risk-factors associated with fractures. RESULTS: The average time (+/-SE) between transplant and the first BMD was 1.2+/-0.05 years. The time interval between the two BMD measurements was 2.5+/-0.05 years. There were 239 and 315 patients in the no-bisphosphonate and bisphosphonate groups, respectively. Treatment was associated with significant preservation of bone loss at the femoral neck (HR 1.56, 95% CI 1.21-2.06, P=0.0007). However, there was no association between bone loss at the femoral neck and fractures regardless of bisphosphonate therapy. Stepwise Cox regression analyses showed that type-1 diabetes, baseline femoral neck T-score, interleukin-2 receptor blockade, and proteinuria (HR 2.02, 0.69, 0.4, 1.23 respectively, P<0.01), but not bisphosphonates, were associated with the risk of fracture. CONCLUSIONS: Bisphosphonates may prevent bone loss in long-term kidney transplant recipients. However, these data suggest a limited role for the initiation of therapy after the first posttransplant year to prevent fractures.

Effects of Schizophrenia Polygenic Risk Scores on Brain Activity and Performance During Working Memory Subprocesses in Healthy Young Adults.

Recent work has begun to shed light on the neural correlates and possible mechanisms of polygenic risk for schizophrenia. Here, we map a schizophrenia polygenic risk profile score (PRS) based on genome-wide association study significant loci onto variability in the activity and functional connectivity of a frontoparietal network supporting the manipulation versus maintenance of information during a numerical working memory (WM) task in healthy young adults (n = 99, mean age = 19.8). Our analyses revealed that higher PRS was associated with hypoactivity of the dorsolateral prefrontal cortex (dlPFC) during the manipulation but not maintenance of information in WM (r2 = .0576, P = .018). Post hoc analyses revealed that PRS-modulated dlPFC hypoactivity correlated with faster reaction times during WM manipulation (r2 = .0967, P = .002), and faster processing speed (r2 = .0967, P = .003) on a separate behavioral task. These PRS-associated patterns recapitulate dlPFC hypoactivity observed in patients with schizophrenia during central executive manipulation of information in WM on this task.

Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence.

Intelligence is highly heritable1 and a major determinant of human health and well-being2. Recent genome-wide meta-analyses have identified 24 genomic loci linked to variation in intelligence3-7, but much about its genetic underpinnings remains to be discovered. Here, we present a large-scale genetic association study of intelligence (n = 269,867), identifying 205 associated genomic loci (190 new) and 1,016 genes (939 new) via positional mapping, expression quantitative trait locus (eQTL) mapping, chromatin interaction mapping, and gene-based association analysis. We find enrichment of genetic effects in conserved and coding regions and associations with 146 nonsynonymous exonic variants. Associated genes are strongly expressed in the brain, specifically in striatal medium spiny neurons and hippocampal pyramidal neurons. Gene set analyses implicate pathways related to nervous system development and synaptic structure. We confirm previous strong genetic correlations with multiple health-related outcomes, and Mendelian randomization analysis results suggest protective effects of intelligence for Alzheimer's disease and ADHD and bidirectional causation with pleiotropic effects for schizophrenia. These results are a major step forward in understanding the neurobiology of cognitive function as well as genetically related neurological and psychiatric disorders.

Reward-related ventral striatum activity links polygenic risk for attention-deficit/hyperactivity disorder to problematic alcohol use in young adulthood.

BACKGROUND: Problematic alcohol use in adolescence and adulthood is a common and often debilitating correlate of childhood attention-deficit/hyperactivity disorder (ADHD). Converging evidence suggests that ADHD and problematic alcohol use share a common additive genetic basis, which may be mechanistically related to reward-related brain function. In the current study, we examined whether polygenic risk for childhood ADHD is linked to problematic alcohol use in young adulthood through alterations in reward-related activity of the ventral striatum, a neural hub supporting appetitive behaviors and reinforcement learning. METHODS: Genomic, neuroimaging, and self-report data were available for 404 non-Hispanic European-American participants who completed the ongoing Duke Neurogenetics Study. Polygenic risk scores for childhood ADHD were calculated based on a genome-wide association study meta-analysis conducted by the Psychiatric Genomics Consortium and tested for association with reward-related ventral striatum activity, measured using a number-guessing functional magnetic resonance imaging paradigm, and self-reported problematic alcohol use. A mediational model tested whether ventral striatum activity indirectly links polygenic risk for ADHD to problematic alcohol use. RESULTS: Despite having no main effect on problematic alcohol use, polygenic risk for childhood ADHD was indirectly associated with problematic alcohol use through increased reward-related ventral striatum activity. CONCLUSIONS: Individual differences in reward-related brain function may, at least in part, mechanistically link polygenic risk for childhood ADHD to problematic alcohol use.

Hypothalamic-pituitary-adrenal axis genetic variation and early stress moderates amygdala function.

Early life stress may precipitate psychopathology, at least in part, by influencing amygdala function. Converging evidence across species suggests that links between childhood stress and amygdala function may be dependent upon hypothalamic-pituitary-adrenal (HPA) axis function. Using data from college-attending non-Hispanic European-Americans (n=308) who completed the Duke Neurogenetics Study, we examined whether early life stress (ELS) and HPA axis genetic variation interact to predict threat-related amygdala function as well as psychopathology symptoms. A biologically-informed multilocus profile score (BIMPS) captured HPA axis genetic variation (FKBP5 rs1360780, CRHR1 rs110402; NR3C2 rs5522/rs4635799) previously associated with its function (higher BIMPS are reflective of higher HPA axis activity). BOLD fMRI data were acquired while participants completed an emotional face matching task. ELS and depression and anxiety symptoms were measured using the childhood trauma questionnaire and the mood and anxiety symptom questionnaire, respectively. The interaction between HPA axis BIMPS and ELS was associated with right amygdala reactivity to threat-related stimuli, after accounting for multiple testing (empirical-p=0.016). Among individuals with higher BIMPS (i.e., the upper 21.4%), ELS was positively coupled with threat-related amygdala reactivity, which was absent among those with average or low BIMPS. Further, higher BIMPS were associated with greater self-reported anxious arousal, though there was no evidence that amygdala function mediated this relationship. Polygenic variation linked to HPA axis function may moderate the effects of early life stress on threat-related amygdala function and confer risk for anxiety symptomatology. However, what, if any, neural mechanisms may mediate the relationship between HPA axis BIMPS and anxiety symptomatology remains unclear.

Large-Scale Cognitive GWAS Meta-Analysis Reveals Tissue-Specific Neural Expression and Potential Nootropic Drug Targets.

Here, we present a large (n = 107,207) genome-wide association study (GWAS) of general cognitive ability ("g"), further enhanced by combining results with a large-scale GWAS of educational attainment. We identified 70 independent genomic loci associated with general cognitive ability. Results showed significant enrichment for genes causing Mendelian disorders with an intellectual disability phenotype. Competitive pathway analysis implicated the biological processes of neurogenesis and synaptic regulation, as well as the gene targets of two pharmacologic agents: cinnarizine, a T-type calcium channel blocker, and LY97241, a potassium channel inhibitor. Transcriptome-wide and epigenome-wide analysis revealed that the implicated loci were enriched for genes expressed across all brain regions (most strongly in the cerebellum). Enrichment was exclusive to genes expressed in neurons but not oligodendrocytes or astrocytes. Finally, we report genetic correlations between cognitive ability and disparate phenotypes including psychiatric disorders, several autoimmune disorders, longevity, and maternal age at first birth.

Interactions Between Anandamide and Corticotropin-Releasing Factor Signaling Modulate Human Amygdala Function and Risk for Anxiety Disorders: An Imaging Genetics Strategy for Modeling Molecular Interactions.

BACKGROUND: Preclinical models reveal that stress-induced amygdala activity and impairment in fear extinction reflect reductions in anandamide driven by corticotropin-releasing factor receptor type 1 (CRF1) potentiation of the anandamide catabolic enzyme fatty acid amide hydrolase. METHODS: Here, we provide clinical translation for the importance of these molecular interactions using an imaging genetics strategy to examine whether interactions between genetic polymorphisms associated with differential anandamide (FAAH rs324420) and CRF1 (CRHR1 rs110402) signaling modulate amygdala function and anxiety disorder diagnosis. RESULTS: Analyses revealed that individuals with a genetic background predicting relatively high anandamide and CRF1 signaling exhibited blunted basolateral amygdala habituation, which further mediated increased risk for anxiety disorders among these same individuals. CONCLUSIONS: The convergence of preclinical and clinical data suggests that interactions between anandamide and CRF1 represent a fundamental molecular mechanism regulating amygdala function and anxiety. Our results further highlight the potential of imaging genetics to powerfully translate complex preclinical findings to clinically meaningful human phenotypes.

PER1 rs3027172 Genotype Interacts with Early Life Stress to Predict Problematic Alcohol Use, but Not Reward-Related Ventral Striatum Activity.

Increasing evidence suggests that the circadian and stress regulatory systems contribute to alcohol use disorder (AUD) risk, which may partially arise through effects on reward-related neural function. The C allele of the PER1 rs3027172 single nucleotide polymorphism (SNP) reduces PER1 expression in cells incubated with cortisol and has been associated with increased risk for adult AUD and problematic drinking among adolescents exposed to high levels of familial psychosocial adversity. Using data from undergraduate students who completed the ongoing Duke Neurogenetics Study (DNS) (n = 665), we tested whether exposure to early life stress (ELS; Childhood Trauma Questionnaire) moderates the association between rs3027172 genotype and later problematic alcohol use (Alcohol Use Disorders Identification Test) as well as ventral striatum (VS) reactivity to reward (card-guessing task while functional magnetic resonance imaging data were acquired). Initial analyses found that PER1 rs3027172 genotype interacted with ELS to predict both problematic drinking and VS reactivity; minor C allele carriers, who were also exposed to elevated ELS reported greater problematic drinking and exhibited greater ventral striatum reactivity to reward-related stimuli. When gene × covariate and environment × covariate interactions were controlled for, the interaction predicting problematic alcohol use remained significant (p < 0.05, corrected) while the interaction predicting VS reactivity was no longer significant. These results extend our understanding of relationships between PER1 genotype, ELS, and problematic alcohol use, and serve as a cautionary tale on the importance of controlling for potential confounders in studies of moderation including gene × environment interactions.

Virtual research visits and direct-to-consumer genetic testing in Parkinson's disease.

OBJECTIVE: The purpose of this study was to conduct a proof-of-concept study to evaluate remote recruitment and assessment of individuals ("virtual research visits") with Parkinson's disease who have pursued direct-to-consumer genetic testing. METHODS: Participants in 23andMe's "Parkinson's Research Community" were contacted by 23andMe. Fifty willing participants living in 23 states underwent a remote, standardized assessment including cognitive and motor tests by a neurologist via video conferencing and then completed a survey. Primary outcomes assessed were (a) proportion of participants who completed the remote assessments; (b) level of agreement (using Cohen's kappa coefficient) of patient-reported data with that of a neurologist; and (c) interest in future virtual research visits. RESULTS: The self-reported diagnosis of Parkinson's disease was confirmed in all cases (k = 1.00). The level of agreement for age of symptom onset (k = 0.97) and family history (k = 0.85) was very good but worse for falling (k = 0.59), tremor (k = 0.56), light-headedness (k = 0.31), and urine control (k = 0.15). Thirty-eight (76%) of the 50 participants completed a post-assessment survey, and 87% of respondents said they would be more or much more willing to participate in future clinical trials if they could do research visits remotely. CONCLUSION: Remote clinical assessments of individuals with known genotypes were conducted nationally and rapidly from a single site, confirmed self-reported diagnosis, and were received favorably. Direct-to-consumer genetic testing and virtual research visits together may enable characterization of genotype and phenotype for geographically diverse populations.

FRAS1-related extracellular matrix 3 (FREM3) single-nucleotide polymorphism effects on gene expression, amygdala reactivity and perceptual processing speed: An accelerated aging pathway of depression risk.

The A allele of the FRAS1-related extracellular matrix protein 3 (FREM3) rs7676614 single nucleotide polymorphism (SNP) was linked to major depressive disorder (MDD) in an early genome-wide association study (GWAS), and to symptoms of psychomotor retardation in a follow-up investigation. In line with significant overlap between age- and depression-related molecular pathways, parallel work has shown that FREM3 expression in postmortem human brain decreases with age. Here, we probe the effect of rs7676614 on amygdala reactivity and perceptual processing speed, both of which are altered in depression and aging. Amygdala reactivity was assessed using a face-matching BOLD fMRI paradigm in 365 Caucasian participants in the Duke Neurogenetics Study (DNS) (192 women, mean age 19.7 ± 1.2). Perceptual processing speed was indexed by reaction times in the same task and the Trail Making Test (TMT). The effect of rs7676614 on FREM3 mRNA brain expression levels was probed in a postmortem cohort of 169 Caucasian individuals (44 women, mean age 50.8 ± 14.9). The A allele of rs7676614 was associated with blunted amygdala reactivity to faces, slower reaction times in the face-matching condition (p < 0.04), as well as marginally slower performance on TMT Part B (p = 0.056). In the postmortem cohort, the T allele of rs6537170 (proxy for the rs7676614 A allele), was associated with trend-level reductions in gene expression in Brodmann areas 11 and 47 (p = 0.066), reminiscent of patterns characteristic of older age. The low-expressing allele of another FREM3 SNP (rs1391187) was similarly associated with reduced amygdala reactivity and slower TMT Part B speed, in addition to reduced BA47 activity and extraversion (p < 0.05). Together, these results suggest common genetic variation associated with reduced FREM3 expression may confer risk for a subtype of depression characterized by reduced reactivity to environmental stimuli and slower perceptual processing speed, possibly suggestive of accelerated aging.

Monoacylglycerol lipase (MGLL) polymorphism rs604300 interacts with childhood adversity to predict cannabis dependence symptoms and amygdala habituation: Evidence from an endocannabinoid system-level analysis.

Despite evidence for heritable variation in cannabis involvement and the discovery of cannabinoid receptors and their endogenous ligands, no consistent patterns have emerged from candidate endocannabinoid (eCB) genetic association studies of cannabis involvement. Given interactions between eCB and stress systems and associations between childhood stress and cannabis involvement, it may be important to consider childhood adversity in the context of eCB-related genetic variation. We employed a system-level gene-based analysis of data from the Comorbidity and Trauma Study (N = 1,558) to examine whether genetic variation in six eCB genes (anabolism: DAGLA, DAGLB, NAPEPLD; catabolism: MGLL, FAAH; binding: CNR1; SNPs N = 65) and childhood sexual abuse (CSA) predict cannabis dependence symptoms. Significant interactions with CSA emerged for MGLL at the gene level (p = .009), and for rs604300 within MGLL (ΔR2 = .007, p < .001), the latter of which survived SNP-level Bonferroni correction and was significant in an additional sample with similar directional effects (N = 859; ΔR2 = .005, p = .026). Furthermore, in a third sample (N = 312), there was evidence that rs604300 genotype interacts with early life adversity to predict threat-related basolateral amygdala habituation, a neural phenotype linked to the eCB system and addiction (ΔR2 = .013, p = .047). Rs604300 may be related to epigenetic modulation of MGLL expression. These results are consistent with rodent models implicating 2-arachidonoylglycerol (2-AG), an endogenous cannabinoid metabolized by the enzyme encoded by MGLL, in the etiology of stress adaptation related to cannabis dependence, but require further replication.