eQTL of KCNK2 regionally influences the brain sulcal widening: evidence from 15,597 UK Biobank participants with neuroimaging data

Validation of TREK1 ion channel activators as an immunomodulatory and neuroprotective strategy in neuroinflammation

Modulation of two-pore domain potassium (K_2P) channels has emerged as a novel field of therapeutic strategies as they may regulate immune cell activation and metabolism, inflammatory signals, or barrier integrity. One of these ion channels is the TWIK-related potassium channel 1 (TREK1). In the current study, we report the identification and validation of new TREK1 activators. Firstly, we used a modified potassium ion channel assay to perform high-throughput-screening of new TREK1 activators. Dose-response studies helped to identify compounds with a high separation between effectiveness and toxicity. Inside-out patch-clamp measurements of Xenopus laevis oocytes expressing TREK1 were used for further validation of these activators regarding specificity and activity. These approaches yielded three substances, E1, B3 and A2 that robustly activate TREK1. Functionally, we demonstrated that these compounds reduce levels of adhesion molecules on primary human brain and muscle endothelial cells without affecting cell viability. Finally, we studied compound A2 via voltage-clamp recordings as this activator displayed the strongest effect on adhesion molecules. Interestingly, A2 lacked TREK1 activation in the tested neuronal cell type. Taken together, this study provides data on novel TREK1 activators that might be employed to pharmacologically modulate TREK1 activity.

Genetic variations within human gained enhancer elements affect human brain sulcal morphology

The expansion of the cerebral cortex is one of the most distinctive changes in the evolution of the human brain. Cortical expansion and related increases in cortical folding may have contributed to emergence of our capacities for high-order cognitive abilities. Molecular analysis of humans, archaic hominins, and non-human primates has allowed identification of chromosomal regions showing evolutionary changes at different points of our phylogenetic history. In this study, we assessed the contributions of genomic annotations spanning 30 million years to human sulcal morphology measured via MRI in more than 18,000 participants from the UK Biobank. We found that variation within brain-expressed human gained enhancers, regulatory genetic elements that emerged since our last common ancestor with Old World monkeys, explained more trait heritability than expected for the left and right calloso-marginal posterior fissures and the right central sulcus. Intriguingly, these are sulci that have been previously linked to the evolution of locomotion in primates and later on bipedalism in our hominin ancestors.

Genetic map of regional sulcal morphology in the human brain from UK biobank data

Genetic associations with macroscopic brain structure can provide insights into brain function and disease. However, specific associations with measures of local brain folding are largely under-explored. Here, we conducted large-scale genome- and exome-wide associations of regional cortical sulcal measures derived from magnetic resonance imaging scans of 40,169 individuals in UK Biobank. We discovered 388 regional brain folding associations across 77 genetic loci, with genes in associated loci enriched for expression in the cerebral cortex, neuronal development processes, and differential regulation during early brain development. We integrated brain eQTLs to refine genes for various loci, implicated several genes involved in neurodevelopmental disorders, and highlighted global genetic correlations with neuropsychiatric phenotypes. We provide an interactive 3D visualisation of our summary associations, emphasising added resolution of regional analyses. Our results offer new insights into the genetic architecture of brain folding and provide a resource for future studies of sulcal morphology in health and disease.

Browsing Multiple Subjects When the Atlas Adaptation Cannot Be Achieved via a Warping Strategy

Brain mapping studies often need to identify brain structures or functional circuits into a set of individual brains. To this end, multiple atlases have been published to represent such structures based on different modalities, subject sets, and techniques. The mainstream approach to exploit these atlases consists in spatially deforming each individual data onto a given atlas using dense deformation fields, which supposes the existence of a continuous mapping between atlases and individuals. However, this continuity is not always verified, and this "iconic" approach has limits. We present in this study an alternative, complementary, "structural" approach, which consists in extracting structures from the individual data, and comparing them without deformation. A "structural atlas" is thus a collection of annotated individual data with a common structure nomenclature. It may be used to characterize structure shape variability across individuals or species, or to train machine learning systems. This study exhibits Anatomist, a powerful structural 3D visualization software dedicated to building, exploring, and editing structural atlases involving a large number of subjects. It has been developed primarily to decipher the cortical folding variability; cortical sulci vary enormously in both size and shape, and some may be missing or have various topologies, which makes iconic approaches inefficient to study them. We, therefore, had to build structural atlases for cortical sulci, and use them to train sulci identification algorithms. Anatomist can display multiple subject data in multiple views, supports all kinds of neuroimaging data, including compound structural object graphs, handles arbitrary coordinate transformation chains between data, and has multiple display features. It is designed as a programming library in both C++ and Python languages, and may be extended or used to build dedicated custom applications. Its generic design makes all the display and structural aspects used to explore the variability of the cortical folding pattern work in other applications, for instance, to browse axonal fiber bundles, deep nuclei, functional activations, or other kinds of cortical parcellations. Multimodal, multi-individual, or inter-species display is supported, and adaptations to large scale screen walls have been developed. These very original features make it a unique viewer for structural atlas browsing.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.

Orchestration of Ion Channels and Transporters in Neocortical Development and Neurological Disorders

Electrical activity plays crucial roles in neural circuit formation and remodeling. During neocortical development, neurons are generated in the ventricular zone, migrate to their correct position, elongate dendrites and axons, and form synapses. In this review, we summarize the functions of ion channels and transporters in neocortical development. Next, we discuss links between neurological disorders caused by dysfunction of ion channels (channelopathies) and neocortical development. Finally, we introduce emerging optical techniques with potential applications in physiological studies of neocortical development and the pathophysiology of channelopathies.

The Two-Pore Domain Potassium Channel TREK-1 Promotes Blood-Brain Barrier Breakdown and Exacerbates Neuronal Death After Focal Cerebral Ischemia in Mice

Earlier studies have shown the neuroprotective role of TWIK-related K⁺ channel 1 (TREK-1) in global cerebral and spinal cord ischemia, while its function in focal cerebral ischemia has long been debated. This study used TREK-1-deficient mice to directly investigate the role of TREK-1 after focal cerebral ischemia. First, immunofluorescence assays in the mouse cerebral cortex indicated that TREK-1 expression was mostly abundant in astrocytes, neurons, and oligodendrocyte precursor cells but was low in myelinating oligodendrocytes, microglia, or endothelial cells. TREK-1 deficiency did not affect brain weight and morphology or the number of neurons, astrocytes, or microglia but did increase glial fibrillary acidic protein (GFAP) expression in astrocytes of the cerebral cortex. The anatomy of the major cerebral vasculature, number and structure of brain micro blood vessels, and blood-brain barrier integrity were unaltered. Next, mice underwent 60 min of focal cerebral ischemia and 72 h of reperfusion induced by the intraluminal suture method. TREK-1-deficient mice showed less neuronal death, smaller infarction size, milder blood-brain barrier (BBB) breakdown, reduced immune cell invasion, and better neurological function. Finally, the specific pharmacological inhibition of TREK-1 also decreased infarction size and improved neurological function. These results demonstrated that TREK-1 might play a detrimental rather than beneficial role in focal cerebral ischemia, and inhibition of TREK-1 would be a strategy to treat ischemic stroke in the clinic.

Advanced high-throughput plant phenotyping techniques for genome-wide association studies: A review

Linking phenotypes and genotypes to identify genetic architectures that regulate important traits is crucial for plant breeding and the development of plant genomics. In recent years, genome-wide association studies (GWASs) have been applied extensively to interpret relationships between genes and traits. Successful GWAS application requires comprehensive genomic and phenotypic data from large populations. Although multiple high-throughput DNA sequencing approaches are available for the generation of genomics data, the capacity to generate high-quality phenotypic data is lagging far behind. Traditional methods for plant phenotyping mostly rely on manual measurements, which are laborious, inaccurate, and time-consuming, greatly impairing the acquisition of phenotypic data from large populations. In contrast, high-throughput phenotyping has unique advantages, facilitating rapid, non-destructive, and high-throughput detection, and, in turn, addressing the shortcomings of traditional methods. Aim of Review: This review summarizes the current status with regard to the integration of high-throughput phenotyping and GWAS in plants, in addition to discussing the inherent challenges and future prospects. Key Scientific Concepts of Review: High-throughput phenotyping, which facilitates non-contact and dynamic measurements, has the potential to offer high-quality trait data for GWAS and, in turn, to enhance the unraveling of genetic structures of complex plant traits. In conclusion, high-throughput phenotyping integration with GWAS could facilitate the revealing of coding information in plant genomes.

Automatic recognition of specific local cortical folding patterns

The study of local cortical folding patterns showed links with psychiatric illnesses as well as cognitive functions. Despite the tools now available to visualize cortical folds in 3D, manually classifying local sulcal patterns is a time-consuming and tedious task. In fact, 3D visualization of folds helps experts to identify different sulcal patterns but fold variability is so high that the distinction between these patterns sometimes requires the definition of complex criteria, making manual classification difficult and not reliable. However, the assessment of the impact of these patterns on the functional organization of the cortex could benefit from the study of large databases, especially when studying rare patterns. In this paper, several algorithms for the automatic classification of fold patterns are proposed to allow morphological studies to be extended and confirmed on such large databases. Three methods are proposed, the first based on a Support Vector Machine (SVM) classifier, the second on the Scoring by Non-local Image Patch Estimator (SNIPE) approach and the third based on a 3D Convolution Neural Network (CNN). These methods are generic enough to be applicable to a wide range of folding patterns. They are tested on two types of patterns for which there is currently no method to automatically identify them: the Anterior Cingulate Cortex (ACC) patterns and the Power Button Sign (PBS). The two ACC patterns are almost equally present whereas PBS is a particularly rare pattern in the general population. The three models proposed achieve balanced accuracies of approximately 80% for ACC patterns classification and 60% for PBS classification. The CNN-based model is more interesting for the classification of ACC patterns thanks to its rapid execution. However, SVM and SNIPE-based models are more effective in managing unbalanced problems such as PBS recognition.

Genome-wide haplotype association study in imaging genetics using whole-brain sulcal openings of 16,304 UK Biobank subjects

Neuroimaging-genetics cohorts gather two types of data: brain imaging and genetic data. They allow the discovery of associations between genetic variants and brain imaging features. They are invaluable resources to study the influence of genetics and environment in the brain features variance observed in normal and pathological populations. This study presents a genome-wide haplotype analysis for 123 brain sulcus opening value (a measure of sulcal width) across the whole brain that include 16,304 subjects from UK Biobank. Using genetic maps, we defined 119,548 blocks of low recombination rate distributed along the 22 autosomal chromosomes and analyzed 1,051,316 haplotypes. To test associations between haplotypes and complex traits, we designed three statistical approaches. Two of them use a model that includes all the haplotypes for a single block, while the last approach considers each haplotype independently. All the statistics produced were assessed as rigorously as possible. Thanks to the rich imaging dataset at hand, we used resampling techniques to assess False Positive Rate for each statistical approach in a genome-wide and brain-wide context. The results on real data show that genome-wide haplotype analyses are more sensitive than single-SNP approach and account for local complex Linkage Disequilibrium (LD) structure, which makes genome-wide haplotype analysis an interesting and statistically sound alternative to the single-SNP counterpart.

Intergenerational Transmission of Cortical Sulcal Patterns from Mothers to their Children

Intergenerational effects are described as the genetic, epigenetic, as well as pre- and postnatal environmental influence parents have on their offspring's behavior, cognition, and brain. During fetal brain development, the primary cortical sulci emerge with a distinctive folding pattern that are under strong genetic influence and show little change of this pattern throughout postnatal brain development. We examined intergenerational transmission of cortical sulcal patterns by comparing primary sulcal patterns between children (N = 16, age 5.5 ± 0.81 years, 8 males) and their biological mothers (N = 15, age 39.72 ± 4.68 years) as well as between children and unrelated adult females. Our graph-based sulcal pattern comparison method detected stronger sulcal pattern similarity for child-mother pairs than child-unrelated pairs, where higher similarity between child-mother pairs was observed mostly for the right lobar regions. Our results also show that child-mother versus child-unrelated pairs differ for daughters and sons with a trend toward significance, particularly for the left hemisphere lobar regions. This is the first study to reveal significant intergenerational transmission of cortical sulcal patterns, and our results have important implications for the study of the heritability of complex behaviors, brain-based disorders, the identification of biomarkers, and targets for interventions.

A slower rate of sulcal widening in the brains of the nondemented oldest old

The relationships between aging and brain morphology have been reported in many previous structural brain studies. However, the trajectories of successful brain aging in the extremely old remain underexplored. In the limited research on the oldest old, covering individuals aged 85 years and older, there are very few studies that have focused on the cortical morphology, especially cortical sulcal features. In this paper, we measured sulcal width and depth as well as cortical thickness from T1-weighted scans of 290 nondemented community-dwelling participants aged between 76 and 103 years. We divided the participants into young old (between 76 and 84; mean = 80.35±2.44; male/female = 76/88) and oldest old (between 85 and 103; mean = 91.74±5.11; male/female = 60/66) groups. The results showed that most of the examined sulci significantly widened with increased age and that the rates of sulcal widening were lower in the oldest old. The spatial pattern of the cortical thinning partly corresponded with that of sulcal widening. Compared to females, males had significantly wider sulci, especially in the oldest old. This study builds a foundation for future investigations of neurocognitive disorders and neurodegenerative diseases in the oldest old, including centenarians.

The reliability and heritability of cortical folds and their genetic correlations across hemispheres

Cortical folds help drive the parcellation of the human cortex into functionally specific regions. Variations in the length, depth, width, and surface area of these sulcal landmarks have been associated with disease, and may be genetically mediated. Before estimating the heritability of sulcal variation, the extent to which these metrics can be reliably extracted from in-vivo MRI must be established. Using four independent test-retest datasets, we found high reliability across the brain (intraclass correlation interquartile range: 0.65-0.85). Heritability estimates were derived for three family-based cohorts using variance components analysis and pooled (total N > 3000); the overall sulcal heritability pattern was correlated to that derived for a large population cohort (N > 9000) calculated using genomic complex trait analysis. Overall, sulcal width was the most heritable metric, and earlier forming sulci showed higher heritability. The inter-hemispheric genetic correlations were high, yet select sulci showed incomplete pleiotropy, suggesting hemisphere-specific genetic influences.

Automatic labeling of cortical sulci using patch- or CNN-based segmentation techniques combined with bottom-up geometric constraints

The extreme variability of the folding pattern of the human cortex makes the recognition of cortical sulci, both automatic and manual, particularly challenging. Reliable identification of the human cortical sulci in its entirety, is extremely difficult and is practiced by only a few experts. Moreover, these sulci correspond to more than a hundred different structures, which makes manual labeling long and fastidious and therefore limits access to large labeled databases to train machine learning. Here, we seek to improve the current model proposed in the Morphologist toolbox, a widely used sulcus recognition toolbox included in the BrainVISA package. Two novel approaches are proposed: patch-based multi-atlas segmentation (MAS) techniques and convolutional neural network (CNN)-based approaches. Both are currently applied for anatomical segmentations because they embed much better representations of inter-subject variability than approaches based on a single template atlas. However, these methods typically focus on voxel-wise labeling, disregarding certain geometrical and topological properties of interest for sulcus morphometry. Therefore, we propose to refine these approaches with domain specific bottom-up geometric constraints provided by the Morphologist toolbox. These constraints are utilized to provide a single sulcus label to each topologically elementary fold, the building blocks of the pattern recognition problem. To eliminate the shortcomings associated with the Morphologist's pre-segmentation into elementary folds, we complement this regularization scheme using a top-down perspective which triggers an additional cleavage of the elementary folds when required. All the newly proposed models outperform the current Morphologist model, the most efficient being a CNN U-Net-based approach which carries out sulcus recognition within a few seconds.

Brain age prediction using deep learning uncovers associated sequence variants

Machine learning algorithms can be trained to estimate age from brain structural MRI. The difference between an individual’s predicted and chronological age, predicted age difference (PAD), is a phenotype of relevance to aging and brain disease. Here, we present a new deep learning approach to predict brain age from a T1-weighted MRI. The method was trained on a dataset of healthy Icelanders and tested on two datasets, IXI and UK Biobank, utilizing transfer learning to improve accuracy on new sites. A genome-wide association study (GWAS) of PAD in the UK Biobank data (discovery set: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=12378$$\end{document}N=12378, replication set: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N=4456$$\end{document}N=4456) yielded two sequence variants, rs1452628-T (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta =-0.08$$\end{document}β=−0.08, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P=1.15\times{10}^{-9}$$\end{document}P=1.15×10−9) and rs2435204-G (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta =0.102$$\end{document}β=0.102, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P=9.73\times 1{0}^{-12}$$\end{document}P=9.73×10−12). The former is near KCNK2 and correlates with reduced sulcal width, whereas the latter correlates with reduced white matter surface area and tags a well-known inversion at 17q21.31 (H2).

Machine learning algorithms can be trained to estimate age from brain structural MRI. Here, the authors introduce a new deep-learning-based age prediction approach, and then carry out a GWAS of the difference between predicted and chronological age, revealing two associated variants.