A genome-wide search for quantitative trait loci affecting the cortical surface area and thickness of Heschl's gyrus

Localization Pattern of Dispatched Homolog 2 (DISP2) in the Central and Enteric Nervous System

Dispatched homolog (DISP) proteins have been implicated in the regulation of hedgehog signaling during embryologic development. Although DISP2 has recently been associated with neuronal development and control of cognitive functions, its localization pattern in the mammalian central and peripheral nervous system has not yet been investigated. In this study, the Disp2 expression profile was assessed in human tissues from publicly available transcriptomic datasets. The DISP2 localization pattern was further characterized in the human and rat central nervous system (CNS), as well as within the colonic enteric nervous system (ENS) using dual-label immunohistochemistry. Colocalization of DISP2 with neuronal and glial markers was additionally analyzed in murine primary ENS culture. At transcriptomic level, DISP2 expression was predominant in neuronal cell types of the CNS and ENS. DISP2 immunoreactivity was mainly located within PGP9.5-positive neurons rather than in S100-positive glial cells throughout the nervous system. Investigation of human and rat brain tissues, colonic specimens, and murine ENS primary cultures revealed that DISP2 was located in neuronal cell somata, as well as along neuronal processes both in the human and murine CNS and ENS. Our results indicate that DISP2 is prominently localized within neuronal cells of the CNS and ENS and support putative functions of DISP2 in these tissues.

Causal analysis of plasma IL-8 on carotid intima media thickness, a measure of subclinical atherosclerosis

BACKGROUND

We investigated the causality of IL-8 on carotid intima-media thickness (c-IMT), a measure of sub-clinical atherosclerosis.

METHODS

The IMPROVE is a multicenter European study (n = 3,711). The association of plasma IL-8 with c-IMT (mm) was estimated by quantile regression. Genotyping was performed using the Illumina CardioMetabo and Immuno chips. Replication was attempted in three independent studies and a meta-analysis was performed using a random model.

RESULTS

In IMPROVE, each unit increase in plasma IL-8 was associated with an increase in median c-IMT measures (all p<0·03) in multivariable analyses. Linear regression identified rs117518778 and rs8057084 as associated with IL-8 levels and with measures of c-IMT. The two SNPs were combined in an IL-8-increasing genetic risk that showed causality of IL-8 on c-IMT in IMPROVE and in the UK Biobank (n = 22,179). The effect of IL-8 on c-IMT measures was confirmed in PIVUS (n = 1,016) and MDCCC (n = 6,103). The association of rs8057084 with c-IMT was confirmed in PIVUS and UK Biobank with a pooled estimate effect (β) of -0·006 with 95%CI (-0·008- -0·003).

CONCLUSION

Our results indicate that genetic variants associated with plasma IL-8 also associate with c-IMT. However, we cannot infer causality of this association, as these variants lie outside of the IL8 locus.

In-vivo data-driven parcellation of Heschl's gyrus using structural connectivity

The human auditory cortex around Heschl’s gyrus (HG) exhibits diverging patterns across individuals owing to the heterogeneity of its substructures. In this study, we investigated the subregions of the human auditory cortex using data-driven machine-learning techniques at the individual level and assessed their structural and functional profiles. We studied an openly accessible large dataset of the Human Connectome Project and identified the subregions of the HG in humans using data-driven clustering techniques with individually calculated imaging features of cortical folding and structural connectivity information obtained via diffusion magnetic resonance imaging tractography. We characterized the structural and functional profiles of each HG subregion according to the cortical morphology, microstructure, and functional connectivity at rest. We found three subregions. The first subregion (HG₁) occupied the central portion of HG, the second subregion (HG₂) occupied the medial-posterior-superior part of HG, and the third subregion (HG₃) occupied the lateral-anterior-inferior part of HG. The HG₃ exhibited strong structural and functional connectivity to the association and paralimbic areas, and the HG₁ exhibited a higher myelin density and larger cortical thickness than other subregions. A functional gradient analysis revealed a gradual axis expanding from the HG₂ to the HG₃. Our findings clarify the individually varying structural and functional organization of human HG subregions and provide insights into the substructures of the human auditory cortex.

Kalirin as a Novel Treatment Target for Cognitive Dysfunction in Schizophrenia

The cognitive dysfunction experienced by patients with schizophrenia represents a major unmet clinical need. We believe that enhancing synaptic function and plasticity by targeting kalirin may provide a novel means to remediate these symptoms. Karilin (a protein encoded by the KALRN gene) has multiple functional domains, including two Dbl homology (DH) guanine exchange factor (GEF) domains, which act to enhance the activity of the Rho family guanosine triphosphate (GTP)-ases. Here, we provide an overview of kalirin's roles in brain function and its therapeutic potential in schizophrenia. We outline how it mediates diverse effects via a suite of distinct isoforms that couple to members of the Rho GTPase family to regulate synapse formation and stabilisation, and how genomic and post-mortem data implicate it in schizophrenia. We then review the current state of knowledge about the influence of kalirin on brain function at a systems level, based largely on evidence from transgenic mouse models, which support its proposed role in regulating dendritic spine function and plasticity. We demonstrate that, whilst the GTPases are classically considered to be 'undruggable', targeting kalirin and other Rho GEFs provides a means to indirectly modulate their activity. Finally, we integrate across the information presented to assess the therapeutic potential of kalirin for schizophrenia and highlight the key outstanding questions required to advance it in this capacity; namely, the need for more information about the diversity and function of its isoforms, how these change across neurodevelopment, and how they affect brain function in vivo.

Getting Started with the IDG KMC Datasets and Tools

The Illuminating the Druggable Genome (IDG) consortium is a National Institutes of Health (NIH) Common Fund program designed to enhance our knowledge of under-studied proteins, more specifically, proteins unannotated within the three most commonly drug-targeted protein families: G-protein coupled receptors, ion channels, and protein kinases. Since 2014, the IDG Knowledge Management Center (IDG-KMC) has generated several open-access datasets and resources that jointly serve as a highly translational machine-learning-ready knowledgebase focused on human protein-coding genes and their products. The goal of the IDG-KMC is to develop comprehensive integrated knowledge for the druggable genome to illuminate the uncharacterized or poorly annotated portion of the druggable genome. The tools derived from the IDG-KMC provide either user-friendly visualizations or ways to impute the knowledge about potential targets using machine learning strategies. In the following protocols, we describe how to use each web-based tool to accelerate illumination in under-studied proteins. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Interacting with the Pharos user interface Basic Protocol 2: Accessing the data in Harmonizome Basic Protocol 3: The ARCHS4 resource Basic Protocol 4: Making predictions about gene function with PrismExp Basic Protocol 5: Using Geneshot to illuminate knowledge about under-studied targets Basic Protocol 6: Exploring under-studied targets with TIN-X Basic Protocol 7: Interacting with the DrugCentral user interface Basic Protocol 8: Estimating Anti-SARS-CoV-2 activities with DrugCentral REDIAL-2020 Basic Protocol 9: Drug Set Enrichment Analysis using Drugmonizome Basic Protocol 10: The Drugmonizome-ML Appyter Basic Protocol 11: The Harmonizome-ML Appyter Basic Protocol 12: GWAS target illumination with TIGA Basic Protocol 13: Prioritizing kinases for lists of proteins and phosphoproteins with KEA3 Basic Protocol 14: Converting PubMed searches to drug sets with the DrugShot Appyter.

Altered Heschl's gyrus duplication pattern in first-episode schizophrenia

BACKGROUND

Reduced gray matter volumes in the superior temporal gyrus and its subregions, such as Heschl's gyrus (HG) and the planum temporale (PT), have been reported in schizophrenia (Sz). However, it remains unclear whether patients exhibit an altered sulcogyral pattern on the superior temporal plane.

METHODS

This magnetic resonance imaging study examined the distribution of HG duplication patterns [i.e., single HG, common stem duplication (CSD), or complete posterior duplication (CPD)] and their relationships with clinical variables and gray matter volumes in the HG and PT of 64 first-episode (FE) patients with Sz and 64 healthy controls.

RESULTS

The prevalence of duplicated HG patterns was significantly higher and gray matter volumes in the HG and PT of both hemispheres were smaller in FESz patients than in healthy controls. The right CPD pattern in the FESz group was associated with less severe positive symptoms. In the FESz and control groups, CSD and CPD patterns correlated with larger volumes in the HG and PT, respectively.

CONCLUSION

The present results revealed an altered HG duplication pattern at the earliest phase of Sz, which may reflect early neurodevelopmental anomalies. However, reduced HG and PT volumes in the FESz were not explained by this sulcogyral pattern only, supporting the complex superior temporal pathology of Sz.

Genome-wide association study reveals new insights into the heritability and genetic correlates of developmental dyslexia

Developmental dyslexia (DD) is a learning disorder affecting the ability to read, with a heritability of 40-60%. A notable part of this heritability remains unexplained, and large genetic studies are warranted to identify new susceptibility genes and clarify the genetic bases of dyslexia. We carried out a genome-wide association study (GWAS) on 2274 dyslexia cases and 6272 controls, testing associations at the single variant, gene, and pathway level, and estimating heritability using single-nucleotide polymorphism (SNP) data. We also calculated polygenic scores (PGSs) based on large-scale GWAS data for different neuropsychiatric disorders and cortical brain measures, educational attainment, and fluid intelligence, testing them for association with dyslexia status in our sample. We observed statistically significant (p  < 2.8 × 10^-6) enrichment of associations at the gene level, for LOC388780 (20p13; uncharacterized gene), and for VEPH1 (3q25), a gene implicated in brain development. We estimated an SNP-based heritability of 20-25% for DD, and observed significant associations of dyslexia risk with PGSs for attention deficit hyperactivity disorder (at p_T = 0.05 in the training GWAS: OR = 1.23[1.16; 1.30] per standard deviation increase; p  = 8 × 10^-13), bipolar disorder (1.53[1.44; 1.63]; p = 1 × 10^-43), schizophrenia (1.36[1.28; 1.45]; p = 4 × 10^-22), psychiatric cross-disorder susceptibility (1.23[1.16; 1.30]; p = 3 × 10^-12), cortical thickness of the transverse temporal gyrus (0.90[0.86; 0.96]; p = 5 × 10^-4), educational attainment (0.86[0.82; 0.91]; p = 2 × 10^-7), and intelligence (0.72[0.68; 0.76]; p = 9 × 10^-29). This study suggests an important contribution of common genetic variants to dyslexia risk, and novel genomic overlaps with psychiatric conditions like bipolar disorder, schizophrenia, and cross-disorder susceptibility. Moreover, it revealed the presence of shared genetic foundations with a neural correlate previously implicated in dyslexia by neuroimaging evidence.

KALRN: A central regulator of synaptic function and synaptopathies

The synaptic regulator, kalirin, plays a key role in synaptic plasticity and formation of dendritic arbors and spines. Dysregulation of the KALRN gene has been linked to various neurological disorders, including autism spectrum disorder, Alzheimer's disease, schizophrenia, addiction and intellectual disabilities. Both genetic and molecular studies highlight the importance of normal KALRN expression for healthy neurodevelopment and function. This review aims to give an in-depth analysis of the structure and molecular mechanisms of kalirin function, particularly within the brain. These data are correlated to genetic evidence of patient mutations within KALRN and animal models of Kalrn that together give insight into the manner in which this gene may be involved in neurodevelopment and the etiology of disease. The emerging links to human disease from post-mortem, genome wide association (GWAS) and exome sequencing studies are examined to highlight the disease relevance of kalirin, particularly in neurodevelopmental diseases. Finally, we will discuss efforts to pharmacologically regulate kalirin protein activity and the implications of such endeavors for the treatment of human disease. As multiple disease states arise from deregulated synapse formation and altered KALRN expression and function, therapeutics may be developed to provide control over KALRN activity and thus synapse dysregulation. As such, a detailed understanding of how kalirin regulates neuronal development, and the manner in which kalirin dysfunction promotes neurological disease, may support KALRN as a valuable therapeutic avenue for future pharmacological intervention.

Identification of susceptibility variants to benign childhood epilepsy with centro-temporal spikes (BECTS) in Chinese Han population

BACKGROUND

Benign Childhood Epilepsy with Centro-temporal Spikes (BECTS) is the most common form of idiopathic epilepsy in children, accounting for up to 23% of pediatric epilepsy. The pathogenesis of BECTS is unknown, but it is thought that genetic factors play a role in susceptibility to the disease.

METHODS

To investigate the role of common genetic variants in BECTS pathogenesis, a 2-stage genome-wide association study (GWAS) was performed in 1,800 Chinese Han BECTS patients, and 7,090 healthy controls. Genetic findings were used in a Mendelian Randomization study in the UK Biobank dataset to investigate the potential role of smoking in BECTS.

FINDINGS

Definitive evidence of a role for common-variant heritability was demonstrated, with heritability of BECTS of >10% observed even with conservative disease prevalence assumptions. Although no individual locus achieved genome-wide significance, twelve loci achieved suggestive evidence of association (5 × 10^-8<P<10^-5). Using combined genetic and brain tissue gene expression data analyzed by Summary-data-based Mendelian Randomization (SMR), causative association of BECTS was demonstrated with SNP rs1948 and the CHRNA5 t3603436 transcript (P_eqtl = 2·10 × 10^-12, P_smr = 7·9 × 10^-5). This finding indicates rs1948 is significantly associated with BECTS through effects on expression of CHRNA5 in brain tissue. The identification of novel loci suggests involvements of KALRN and the CHRNA5-A3-B4 cluster in BECTS. Using a generalized SMR approach we demonstrate that maternal smoking around birth is significantly associated with increased risk of BECTS (odds ratio = 3·90, P = 0·0099).

INTERPRETATION

This study shows that BECTS risk is at least partially heritable and due to common genetic variants. Additionally, we demonstrate that BECTS risk is substantially increased by maternal smoking around birth.

Enhancer Locus in ch14q23.1 Modulates Brain Asymmetric Temporal Regions Involved in Language Processing

Identifying the genes that contribute to the variability in brain regions involved in language processing may shed light on the evolution of brain structures essential to the emergence of language in Homo sapiens. The superior temporal asymmetrical pit (STAP), which is not observed in chimpanzees, represents an ideal phenotype to investigate the genetic variations that support human communication. The left STAP depth was significantly associated with a predicted enhancer annotation located in the 14q23.1 locus, between DACT1 and KIAA0586, in the UK Biobank British discovery sample (N = 16 515). This association was replicated in the IMAGEN cohort (N = 1726) and the UK Biobank non-British validation sample (N = 2161). This genomic region was also associated to a lesser extent with the right STAP depth and the formation of sulcal interruptions, “plis de passage,” in the bilateral STAP but not with other structural brain MRI phenotypes, highlighting its notable association with the superior temporal regions. Diffusion MRI emphasized an association with the fractional anisotropy of the left auditory fibers of the corpus callosum and with networks involved in linguistic processing in resting-state functional MRI. Overall, this evidence demonstrates a specific relationship between this locus and the establishment of the superior temporal regions that support human communication.

TASH: Toolbox for the Automated Segmentation of Heschl's gyrus

Auditory cortex volume and shape differences have been observed in the context of phonetic learning, musicianship and dyslexia. Heschl’s gyrus, which includes primary auditory cortex, displays large anatomical variability across individuals and hemispheres. Given this variability, manual labelling is the gold standard for segmenting HG, but is time consuming and error prone. Our novel toolbox, called ‘Toolbox for the Automated Segmentation of HG’ or TASH, automatically segments HG in brain structural MRI data, and extracts measures including its volume, surface area and cortical thickness. TASH builds upon FreeSurfer, which provides an initial segmentation of auditory regions, and implements further steps to perform finer auditory cortex delineation. We validate TASH by showing significant relationships between HG volumes obtained using manual labelling and using TASH, in three independent datasets acquired on different scanners and field strengths, and by showing good qualitative segmentation. We also present two applications of TASH, demonstrating replication and extension of previously published findings of relationships between HG volumes and (a) phonetic learning, and (b) musicianship. In sum, TASH effectively segments HG in a fully automated and reproducible manner, opening up a wide range of applications in the domains of expertise, disease, genetics and brain plasticity.

Variants in MED12L, encoding a subunit of the mediator kinase module, are responsible for intellectual disability associated with transcriptional defect

Purpose

Mediator is a multiprotein complex that allows the transfer of genetic information from DNA binding proteins to the RNA polymerase II during transcription initiation. MED12L is a subunit of the kinase module, which is one of the four sub-complexes of the mediator complex. Other subunits of the kinase module have been already implicated in intellectual disability, namely MED12, MED13L, MED13 and CDK19.

Methods

We describe an international cohort of seven affected individuals harboring variants involving MED12L identified by array CGH, exome or genome sequencing.

Results

All affected individuals presented with intellectual disability and/or developmental delay, including speech impairment. Other features included autism spectrum disorder, aggressive behavior, corpus callosum abnormality and mild facial morphological features. Three individuals had a MED12L deletion or duplication. The other four individuals harbored single nucleotide variants (one nonsense, one frameshift and two splicing variants). Functional analysis confirmed a moderate and significant alteration of RNA synthesis in two individuals.

Conclusion

Overall data suggest that MED12L haploinsufficiency is responsible for intellectual disability and transcriptional defect. Our findings confirm that the integrity of this kinase module is a critical factor for neurological development.

The DCBLD receptor family: emerging signaling roles in development, homeostasis and disease

The discoidin, CUB, and LCCL domain-containing (DCBLD) receptor family are composed of the type-I transmembrane proteins DCBLD1 and DCBLD2 (also ESDN and CLCP1). These proteins are highly conserved across vertebrates and possess similar domain structure to that of neuropilins, which act as critical co-receptors in developmental processes. Although DCBLD1 remains largely uncharacterized, the functional and mechanistic roles of DCBLD2 are emerging. This review provides a comprehensive discussion of this presumed receptor family, ranging from structural and signaling aspects to their associations with cancer, physiology, and development.

Significantly altered peripheral blood cell DNA methylation profile as a result of immediate effect of metformin use in healthy individuals

BACKGROUND

Metformin is a widely prescribed antihyperglycemic agent that has been also associated with multiple therapeutic effects in various diseases, including several types of malignancies. There is growing evidence regarding the contribution of the epigenetic mechanisms in reaching metformin's therapeutic goals; however, the effect of metformin on human cells in vivo is not comprehensively studied. The aim of our study was to examine metformin-induced alterations of DNA methylation profiles in white blood cells of healthy volunteers, employing a longitudinal study design.

RESULTS

Twelve healthy metformin-naïve individuals where enrolled in the study. Genome-wide DNA methylation pattern was estimated at baseline, 10 h and 7 days after the start of metformin administration. The whole-genome DNA methylation analysis in total revealed 125 differentially methylated CpGs, of which 11 CpGs and their associated genes with the most consistent changes in the DNA methylation profile were selected: POFUT2, CAMKK1, EML3, KIAA1614, UPF1, MUC4, LOC727982, SIX3, ADAM8, SNORD12B, VPS8, and several differentially methylated regions as novel potential epigenetic targets of metformin. The main functions of the majority of top-ranked differentially methylated loci and their representative cell signaling pathways were linked to the well-known metformin therapy targets: regulatory processes of energy homeostasis, inflammatory responses, tumorigenesis, and neurodegenerative diseases.

CONCLUSIONS

Here we demonstrate for the first time the immediate effect of short-term metformin administration at therapeutic doses on epigenetic regulation in human white blood cells. These findings suggest the DNA methylation process as one of the mechanisms involved in the action of metformin, thereby revealing novel targets and directions of the molecular mechanisms underlying the various beneficial effects of metformin.

TRIAL REGISTRATION

EU Clinical Trials Register, 2016-001092-74. Registered 23 March 2017, https://www.clinicaltrialsregister.eu/ctr-search/trial/2016-001092-74/LV .

2q24 deletion in a 9-month old girl with anal atresia, hearing impairment, and hypotonia

Deletion of 2q24.2 is a rare cytogenetic aberration in patients, exhibiting heterogeneous clinical features, and common phenotypes included developmental delay, intellectual disability, hypotonia, and mild dysmorphic features. Hearing impairment and anal atresia are rarely described. Here we described a 9-month-old female patient with hypotonia in all four limbs, developmental delay, and intellectual disability. In addition, congenital anal atresia was diagnosed and treated after birth, and hearing impairment was found in right ear. Single nucleotide polymorphisms (SNP) array detected a 5.2 Mb deletion on 2q24.2q24.3, including 19 genes (ITGB6; TBR1; SLC4A10; KCNH7 SCN3A; SCN2A et al.). Among these genes, it is affirmative that TBR1 is a causative gene for intellectual disability; however, the pathogenic genes of other phenotypes remain unclear. We briefly review the knowledge of genes likely involved in these clinical features, including hearing impairment, anal atresia, and developmental delay.

COBLL1 modulates cell morphology and facilitates androgen receptor genomic binding in advanced prostate cancer

Androgen receptor (AR) signaling is essential for prostate cancer progression and acquiring resistance to hormone therapy. However, the molecular pathogenesis through AR activation has not been fully understood. We performed integrative transcriptomic analysis to compare the AR program in a castration-resistant prostate cancer (CRPC) model with that in their parental hormone-sensitive cells. We found that the gene cordon-bleu-like 1 (COBLL1) is highly induced by AR in CRPC model cells. The expression of COBLL1 that possesses an actin-binding domain is up-regulated in clinical prostate cancer tissues and is associated with a poor prognosis for prostate cancer patients. COBLL1 is involved in the cancer cell morphogenesis to a neuron-like cell shape observed in the CRPC model cells, promoting cell growth and migration. Moreover, nuclear COBLL1 interacts with AR to enhance complex formation with CDK1 and facilitates AR phosphorylation for genomic binding in CRPC model cells. Thus, our findings showed the mechanistic relevance of cordon-bleu proteins during the AR-mediated progression to CRPC.

A Schizophrenia-Linked KALRN Coding Variant Alters Neuron Morphology, Protein Function, and Transcript Stability

BACKGROUND

Large-scale genetic studies have revealed that rare sequence variants, including single nucleotide variants (SNVs), in glutamatergic synaptic genes are enriched in schizophrenia patients. However, the majority are too rare to show any association with disease and have not been examined functionally. One such SNV, KALRN-P2255T, displays a penetrance that greatly exceeds that of previously identified schizophrenia-associated SNVs. Therefore, we sought to characterize its effects on the function of kalirin (Kal)-9, a dual Ras-related C3 botulinum toxin substrate 1 and Ras homologue gene family, member A (RhoA) guanine nucleotide exchange factor, upregulated in human schizophrenia brain tissue.

METHODS

Kal9 was overexpressed in primary rat cortical neurons or human embryonic kidney 293 (HEK293) cells. The effects of the P2255T variant on dendritic branching, dendritic spine morphology, protein and messenger RNA stability, and catalytic activity were examined.

RESULTS

Kal9-P2255T leads to diminished basal dendritic branching and dendritic spine size, compared with wild-type Kal9. The P2255T SNV directly affected Kal9 protein function, causing increased RhoA activation in HEK293 cells, but had no effect on Ras-related C3 botulinum toxin substrate 1 activation. Consistent with human postmortem findings, we found that Kal9-P2255T protein levels were higher than those of wild-type Kal9 in neurons. Increased messenger RNA stability was detected in HEK293 cells, indicating that this was the cause of the higher protein levels. When analyzed together, increased intrinsic RhoA guanine nucleotide exchange factor catalytic activity combined with increased messenger RNA expression led to net enhancement of RhoA activation, known to negatively impact neuronal morphology.

CONCLUSIONS

Taken together, our data reveal a novel mechanism for disease-associated SNVs and provide a platform for modeling morphological changes in mental disorders.

Neuroanatomical profiles of bilingual children

The goal of the present study was to examine differences in cortical thickness, cortical surface area, and subcortical volume between bilingual children who are highly proficient in two languages (i.e., English and Spanish) and bilingual children who are mainly proficient in one of the languages (i.e., Spanish). All children (N = 49) learned Spanish as a native language (L1) at home and English as a second language (L2) at school. Proficiency of both languages was assessed using the standardized Woodcock Language Proficiency Battery. Five-minute high-resolution anatomical scans were acquired with a 3-Tesla scanner. The degree of discrepancy between L1 and L2 proficiency was used to classify the children into two groups: children with balanced proficiency and children with unbalanced proficiency. The groups were comparable on language history, parental education, and other variables except English proficiency. Values of cortical thickness and surface area of the transverse STG, IFG-pars opercularis, and MFG, as well as subcortical volume of the caudate and putamen, were extracted from FreeSurfer. Results showed that children with balanced bilingualism had thinner cortices of the left STG, left IFG, left MFG and a larger bilateral putamen, whereas unbalanced bilinguals showed thicker cortices of the same regions and a smaller putamen. Additionally, unbalanced bilinguals with stronger foreign accents in the L2 showed reduced surface areas of the MFG and STS bilaterally. The results suggest that balanced/unbalanced bilingualism is reflected in different neuroanatomical characteristics that arise from biological and/or environmental factors.

Genome-wide association study of pigmentary traits (skin and iris color) in individuals of East Asian ancestry

Background

Currently, there is limited knowledge about the genetics underlying pigmentary traits in East Asian populations. Here, we report the results of the first genome-wide association study of pigmentary traits (skin and iris color) in individuals of East Asian ancestry.

Methods

We obtained quantitative skin pigmentation measures (M-index) in the inner upper arm of the participants using a portable reflectometer (N = 305). Quantitative measures of iris color (expressed as L*, a* and b* CIELab coordinates) were extracted from high-resolution iris pictures (N = 342). We also measured the color differences between the pupillary and ciliary regions of the iris (e.g., iris heterochromia). DNA samples were genotyped with Illumina’s Infinium Multi-Ethnic Global Array (MEGA) and imputed using the 1000 Genomes Phase 3 samples as reference haplotypes.

Results

For skin pigmentation, we did not observe any genome-wide significant signal. We followed-up in three independent Chinese samples the lead SNPs of five regions showing multiple common markers (minor allele frequency ≥ 5%) with good imputation scores and suggestive evidence of association (p-values < 10⁻⁵). One of these markers, rs2373391, which is located in an intron of the ZNF804B gene on chromosome 7, was replicated in one of the Chinese samples (p = 0.003). For iris color, we observed genome-wide signals in the OCA2 region on chromosome 15. This signal is driven by the non-synonymous rs1800414 variant, which explains 11.9%, 10.4% and 6% of the variation observed in the b*, a* and L* coordinates in our sample, respectively. However, the OCA2 region was not associated with iris heterochromia.

Discussion

Additional genome-wide association studies in East Asian samples will be necessary to further disentangle the genetic architecture of pigmentary traits in East Asian populations.

Genome-wide association study of language performance in Alzheimer's disease

Language impairment is common in prodromal stages of Alzheimer's disease (AD) and progresses over time. However, the genetic architecture underlying language performance is poorly understood. To identify novel genetic variants associated with language performance, we analyzed brain MRI and performed a genome-wide association study (GWAS) using a composite measure of language performance from the Alzheimer's Disease Neuroimaging Initiative (ADNI; n=1560). The language composite score was associated with brain atrophy on MRI in language and semantic areas. GWAS identified GLI3 (GLI family zinc finger 3) as significantly associated with language performance (p<5×10^-8). Enrichment of GWAS association was identified in pathways related to nervous system development and glutamate receptor function and trafficking. Our results, which warrant further investigation in independent and larger cohorts, implicate GLI3, a developmental transcription factor involved in patterning brain structures, as a putative gene associated with language dysfunction in AD.

P2Y Receptors in Immune Response and Inflammation

Metabotropic pyrimidine and purine nucleotide receptors (P2Y receptors) are expressed in virtually all cells with implications in very diverse biological functions, including the well-established platelet aggregation (P2Y12), but also immune regulation and inflammation. The classical P2Y receptors bind nucleotides and are encoded by eight genes with limited sequence homology, while phylogenetically related receptors (e.g., P2Y12-like) recognize lipids and peptides, but also nucleotide derivatives. Growing lines of evidence suggest an important function of P2Y receptors in immune cell differentiation and maturation, migration, and cell apoptosis. Here, we give a perspective on the P2Y receptors' molecular structure and physiological importance in immune cells, as well as the related diseases and P2Y-targeting therapies. Extensive research is being undertaken to find modulators of P2Y receptors and uncover their physiological roles. We anticipate the medical applications of P2Y modulators and their immune relevance.

A genome wide association study of fast beta EEG in families of European ancestry

BACKGROUND

Differences in fast beta (20-28Hz) electroencephalogram (EEG) oscillatory activity distinguish some individuals with psychiatric and substance use disorders, suggesting that it may be a useful endophenotype for studying the genetics of disorders characterized by neural hyper-excitability. Despite the high heritability estimates provided by twin and family studies, there have been relatively few genetic studies of beta EEG, and to date only one genetic association finding has replicated (i.e., GABRA2).

METHOD

In a sample of 1564 individuals from 117 families of European Ancestry (EA) drawn from the Collaborative Study on the Genetics of Alcoholism (COGA), we performed a Genome-Wide Association Study (GWAS) on resting-state fronto-central fast beta EEG power, adjusting regression models for family relatedness, age, sex, and ancestry. To further characterize genetic findings, we examined the functional and behavioral significance of GWAS findings.

RESULTS

Three intronic variants located within DSE (dermatan sulfate epimerase) on 6q22 were associated with fast beta EEG at a genome wide significant level (p<5×10-8). The most significant SNP was rs2252790 (p<2.6×10-8; MAF=0.36; β=0.135). rs2252790 is an eQTL for ROS1 expressed most robustly in the temporal cortex (p=1.2×10-6) and for DSE/TSPYL4 expressed most robustly in the hippocampus (p=7.3×10-4; β=0.29). Previous studies have indicated that DSE is involved in a network of genes integral to membrane organization; gene-based tests indicated that several variants within this network (i.e., DSE, ZEB2, RND3, MCTP1, and CTBP2) were also associated with beta EEG (empirical p<0.05), and of these genes, ZEB2 and CTBP2 were associated with DSM-V Alcohol Use Disorder (AUD; empirical p<0.05).'

DISCUSSION

In this sample of EA families enriched for AUDs, fast beta EEG is associated with variants within DSE on 6q22; the most significant SNP influences the mRNA expression of DSE and ROS1 in hippocampus and temporal cortex, brain regions important for beta EEG activity. Gene-based tests suggest evidence of association with related genes, ZEB2, RND3, MCTP1, CTBP2, and beta EEG. Converging data from GWAS, gene expression, and gene-networks presented in this study provide support for the role of genetic variants within DSE and related genes in neural hyperexcitability, and has highlighted two potential candidate genes for AUD and/or related neurological conditions: ZEB2 and CTBP2. However, results must be replicated in large, independent samples.

Regulatory Variants Modulate Protein Kinase C α (PRKCA) Gene Expression in Human Heart

PURPOSE

Protein kinase C α (PRKCA) is involved in multiple functions and has been implicated in heart failure risks and treatment outcomes. This study aims to identify regulatory variants affecting PRKCA expression in human heart, and evaluate attributable risk of heart disease.

METHODS

mRNA expression quantitative trait loci (eQTLs) were extracted from the Genotype and Tissue Expression Project (GTEx). Allelic mRNA ratios were measured in 51 human heart tissues to identify cis-acting regulatory variants. Potential regulatory regions were tested with luciferase reporter gene assays and further evaluated in GTEx and genome-wide association studies.

RESULTS

Located in a region with robust enhancer activity in luciferase reporter assays, rs9909004 (T > C, minor allele frequency =0.47) resides in a haplotype displaying strong eQTLs for PRKCA in heart (p = 1.2 × 10^-23). The minor C allele is associated with both decreased PRKCA mRNA expression and decreased risk of phenotypes characteristic of heart failure in GWAS analyses (QT interval p = 3.0 × 10^-14). While rs9909004 is the likely regulatory variant, other variants in high linkage disequilibrium cannot be excluded. Distinct regulatory variants appear to affect expression in other tissues.

CONCLUSIONS

The haplotype carrying rs9909004 influences PRKCA expression in the heart and is associated with traits linked to heart failure, potentially affecting therapy of heart failure.

Common Genetic Variant in VIT Is Associated with Human Brain Asymmetry

Brain asymmetry varies across individuals. However, genetic factors contributing to this normal variation are largely unknown. Here we studied variation of cortical surface area asymmetry in a large sample of subjects. We performed principal component analysis (PCA) to capture correlated asymmetry variation across cortical regions. We found that caudal and rostral anterior cingulate together account for a substantial part of asymmetry variation among individuals. To find SNPs associated with this subset of brain asymmetry variation we performed a genome-wide association study followed by replication in an independent cohort. We identified one SNP (rs11691187) that had genome-wide significant association (P_Combined = 2.40e-08). The rs11691187 is in the first intron of VIT. In a follow-up analysis, we found that VIT gene expression is associated with brain asymmetry in six donors of the Allen Human Brain Atlas. Based on these findings we suggest that VIT contributes to normal brain asymmetry variation. Our results can shed light on disorders associated with altered brain asymmetry.

Regional correlations between cortical thickness and surface area asymmetries: A surface-based morphometry study of 250 adults

We report on the patterns of asymmetries of various MRI-derived cortical phenotypes, namely cortical thickness (CT), cortical surface area (CSA), cortical volume (CV) and sulcal depth (SULC), as well as on their cross-relationships. A surface-based approach was designed to obtain cortical asymmetry maps unbiased for inter-hemispheric structural positional misalignment. Accurate vertex-wise asymmetries of CT, CSA, CV, and SULC were so obtained in 250 individuals including 120 left-handers that had been selected from a larger population as having a typical leftward language lateralization. We found no significant effect of handedness on CT, CSA or CV asymmetries, although a trend for a significant difference in the SULC asymmetry of the Rolandic genu was present (leftward in right-handers, rightward in left-handers). At the hemispheric level, we found rightward CT and CSA asymmetries that were not correlated. At the regional level, asymmetry patterns of CT and CSA were marked by a spatial overlap of both concordant and opposite CT and CSA asymmetries. Half of these regions of overlap presented a significant association (positive or negative) between CT and CSA asymmetries. Strikingly, the 3 regions showing positive correlations between CT and CSA asymmetries were those known to exhibit robust asymmetries across methodologies, such as the leftward asymmetrical planum temporale and the rightward asymmetrical superior temporal sulcus and cingulate cortex. This study demonstrates that regional correlations between CT and CSA asymmetries are a characteristic of brain structural organization that could be of significance in the choice of structural markers in studies dealing with the genetic basis of brain lateralization.

Early developmental gene enhancers affect subcortical volumes in the adult human brain

Genome-wide association screens aim to identify common genetic variants contributing to the phenotypic variability of complex traits, such as human height or brain morphology. The identified genetic variants are mostly within noncoding genomic regions and the biology of the genotype-phenotype association typically remains unclear. In this article, we propose a complementary targeted strategy to reveal the genetic underpinnings of variability in subcortical brain volumes, by specifically selecting genomic loci that are experimentally validated forebrain enhancers, active in early embryonic development. We hypothesized that genetic variation within these enhancers may affect the development and ultimately the structure of subcortical brain regions in adults. We tested whether variants in forebrain enhancer regions showed an overall enrichment of association with volumetric variation in subcortical structures of >13,000 healthy adults. We observed significant enrichment of genomic loci that affect the volume of the hippocampus within forebrain enhancers (empirical P = 0.0015), a finding which robustly passed the adjusted threshold for testing of multiple brain phenotypes (cutoff of P < 0.0083 at an alpha of 0.05). In analyses of individual single nucleotide polymorphisms (SNPs), we identified an association upstream of the ID2 gene with rs7588305 and variation in hippocampal volume. This SNP-based association survived multiple-testing correction for the number of SNPs analyzed but not for the number of subcortical structures. Targeting known regulatory regions offers a way to understand the underlying biology that connects genotypes to phenotypes, particularly in the context of neuroimaging genetics. This biology-driven approach generates testable hypotheses regarding the functional biology of identified associations. Hum Brain Mapp 37:1788-1800, 2016. © 2016 Wiley Periodicals, Inc.

Early Childhood Depression and Alterations in the Trajectory of Gray Matter Maturation in Middle Childhood and Early Adolescence

IMPORTANCE

The trajectory of cortical gray matter development in childhood has been characterized by early neurogenesis and volume increase, peaking at puberty followed by selective elimination and myelination, resulting in volume loss and thinning. This inverted U-shaped trajectory, as well as cortical thickness, has been associated with cognitive and emotional function. Synaptic pruning-based volume decline has been related to experience-dependent plasticity in animals. To date, there have been no data to inform whether and how childhood depression might be associated with this trajectory.

OBJECTIVE

To examine the effects of early childhood depression, from the preschool age to the school age period, on cortical gray matter development measured across 3 waves of neuroimaging from late school age to early adolescence.

DESIGN, SETTING, AND PARTICIPANTS

Data were collected in an academic research setting from September 22, 2003, to December 13, 2014, on 193 children aged 3 to 6 years from the St Louis, Missouri, metropolitan area who were observed for up to 11 years in a longitudinal behavioral and neuroimaging study of childhood depression. Multilevel modeling was applied to explore the association between the number of childhood depression symptoms and prior diagnosis of major depressive disorder and the trajectory of gray matter change across 3 scan waves. Data analysis was conducted from October 29, 2014, to September 28, 2015.

MAIN OUTCOMES AND MEASURES

Volume, thickness, and surface area of cortical gray matter measured using structural magnetic resonance imaging at 3 scan waves.

RESULTS

Of the 193 children, 90 had a diagnosis of major depressive disorder; 116 children had 3 full waves of neuroimaging scans. Findings demonstrated marked alterations in cortical gray matter volume loss (slope estimate, -0.93 cm³; 95% CI, -1.75 to -0.10 cm³ per scan wave) and thinning (slope estimate, -0.0044 mm; 95% CI, -0.0077 to -0.0012 mm per scan wave) associated with experiencing an episode of major depressive disorder before the first magnetic resonance imaging scan. In contrast, no significant associations were found between development of gray matter and family history of depression or experiences of traumatic or stressful life events during this period.

CONCLUSIONS AND RELEVANCE

This study demonstrates an association between early childhood depression and the trajectory of cortical gray matter development in late school age and early adolescence. These findings underscore the significance of early childhood depression on alterations in neural development.