A latent clinical-anatomical dimension relating metabolic syndrome to brain structure and cognition

The link between metabolic syndrome (MetS) and neurodegenerative as well as cerebrovascular conditions holds substantial implications for brain health in at-risk populations. This study elucidates the complex relationship between MetS and brain health by conducting a comprehensive examination of cardiometabolic risk factors, brain morphology, and cognitive function in 40,087 individuals. Multivariate, data-driven statistics identified a latent dimension linking more severe MetS to widespread brain morphological abnormalities, accounting for up to 71% of shared variance in the data. This dimension was replicable across sub-samples. In a mediation analysis, we could demonstrate that MetS-related brain morphological abnormalities mediated the link between MetS severity and cognitive performance in multiple domains. Employing imaging transcriptomics and connectomics, our results also suggest that MetS-related morphological abnormalities are linked to the regional cellular composition and macroscopic brain network organization. By leveraging extensive, multi-domain data combined with a dimensional stratification approach, our analysis provides profound insights into the association of MetS and brain health. These findings can inform effective therapeutic and risk mitigation strategies aimed at maintaining brain integrity.

Inhibition of CERS1 in skeletal muscle exacerbates age-related muscle dysfunction

Age-related muscle wasting and dysfunction render the elderly population vulnerable and incapacitated, while underlying mechanisms are poorly understood. Here, we implicate the CERS1 enzyme of the de novo sphingolipid synthesis pathway in the pathogenesis of age-related skeletal muscle impairment. In humans, CERS1 abundance declines with aging in skeletal muscle cells and, correlates with biological pathways involved in muscle function and myogenesis. Furthermore, CERS1 is upregulated during myogenic differentiation. Pharmacological or genetic inhibition of CERS1 in aged mice blunts myogenesis and deteriorates aged skeletal muscle mass and function, which is associated with the occurrence of morphological features typical of inflammation and fibrosis. Ablation of the CERS1 orthologue lagr-1 in Caenorhabditis elegans similarly exacerbates the age-associated decline in muscle function and integrity. We discover genetic variants reducing CERS1 expression in human skeletal muscle and Mendelian randomization analysis in the UK biobank cohort shows that these variants reduce muscle grip strength and overall health. In summary, our findings link age-related impairments in muscle function to a reduction in CERS1, thereby underlining the importance of the sphingolipid biosynthesis pathway in age-related muscle homeostasis.

Resting-state alterations in behavioral variant frontotemporal dementia are related to the distribution of monoamine and GABA neurotransmitter systems

Background

Aside to clinical changes, behavioral variant frontotemporal dementia (bvFTD) is characterized by progressive structural and functional alterations in frontal and temporal regions. We examined if there is a selective vulnerability of specific neurotransmitter systems in bvFTD by evaluating the link between disease-related functional alterations and the spatial distribution of specific neurotransmitter systems and their underlying gene expression levels.

Methods

Maps of fractional amplitude of low-frequency fluctuations (fALFF) were derived as a measure of local activity from resting-state functional magnetic resonance imaging for 52 bvFTD patients (mean age = 61.5 ± 10.0 years; 14 females) and 22 healthy controls (HC) (mean age = 63.6 ± 11.9 years; 13 females). We tested if alterations of fALFF in patients co-localize with the non-pathological distribution of specific neurotransmitter systems and their coding mRNA gene expression. Furthermore, we evaluated if the strength of co-localization is associated with the observed clinical symptoms.

Results

Patients displayed significantly reduced fALFF in frontotemporal and frontoparietal regions. These alterations co-localized with the distribution of serotonin (5-HT1b and 5-HT2a) and γ-aminobutyric acid type A (GABAa) receptors, the norepinephrine transporter (NET), and their encoding mRNA gene expression. The strength of co-localization with NET was associated with cognitive symptoms and disease severity of bvFTD.

Conclusions

Local brain functional activity reductions in bvFTD followed the distribution of specific neurotransmitter systems indicating a selective vulnerability. These findings provide novel insight into the disease mechanisms underlying functional alterations. Our data-driven method opens the road to generate new hypotheses for pharmacological interventions in neurodegenerative diseases even beyond bvFTD.

Funding

This study has been supported by the German Consortium for Frontotemporal Lobar Degeneration, funded by the German Federal Ministry of Education and Research (BMBF; grant no. FKZ01GI1007A).

Huntingtin recruits KIF1A to transport synaptic vesicle precursors along the mouse axon to support synaptic transmission and motor skill learning

Neurotransmitters are released at synapses by synaptic vesicles (SVs), which originate from SV precursors (SVPs) that have traveled along the axon. Because each synapse maintains a pool of SVs, only a small fraction of which are released, it has been thought that axonal transport of SVPs does not affect synaptic function. Here, studying the corticostriatal network both in microfluidic devices and in mice, we find that phosphorylation of the Huntingtin protein (HTT) increases axonal transport of SVPs and synaptic glutamate release by recruiting the kinesin motor KIF1A. In mice, constitutive HTT phosphorylation causes SV over-accumulation at synapses, increases the probability of SV release, and impairs motor skill learning on the rotating rod. Silencing KIF1A in these mice restored SV transport and motor skill learning to wild-type levels. Axonal SVP transport within the corticostriatal network thus influences synaptic plasticity and motor skill learning.

Biological brain age prediction using machine learning on structural neuroimaging data: Multi-cohort validation against biomarkers of Alzheimer's disease and neurodegeneration stratified by sex

Brain-age can be inferred from structural neuroimaging and compared to chronological age (brain-age delta) as a marker of biological brain aging. Accelerated aging has been found in neurodegenerative disorders like Alzheimer's disease (AD), but its validation against markers of neurodegeneration and AD is lacking. Here, imaging-derived measures from the UK Biobank dataset (N=22,661) were used to predict brain-age in 2,314 cognitively unimpaired (CU) individuals at higher risk of AD and mild cognitive impaired (MCI) patients from four independent cohorts with available biomarker data: ALFA+, ADNI, EPAD, and OASIS. Brain-age delta was associated with abnormal amyloid-β, more advanced stages (AT) of AD pathology and APOE-ε4 status. Brain-age delta was positively associated with plasma neurofilament light, a marker of neurodegeneration, and sex differences in the brain effects of this marker were found. These results validate brain-age delta as a non-invasive marker of biological brain aging in non-demented individuals with abnormal levels of biomarkers of AD and axonal injury.

Promoting Fc-Fc interactions between anti-capsular antibodies provides strong immune protection against Streptococcus pneumoniae

Streptococcus pneumoniae is the leading cause of community-acquired pneumonia and an important cause of childhood mortality. Despite the introduction of successful vaccines, the global spread of both non-vaccine serotypes and antibiotic-resistant strains reinforces the development of alternative therapies against this pathogen. One possible route is the development of monoclonal antibodies (mAbs) that induce killing of bacteria via the immune system. Here, we investigate whether mAbs can be used to induce killing of pneumococcal serotypes for which the current vaccines show unsuccessful protection. Our study demonstrates that when human mAbs against pneumococcal capsule polysaccharides (CPS) have a poor capacity to induce complement activation, a critical process for immune protection against pneumococci, their activity can be strongly improved by hexamerization-enhancing mutations. Our data indicate that anti-capsular antibodies may have a low capacity to form higher-order oligomers (IgG hexamers) that are needed to recruit complement component C1. Indeed, specific point mutations in the IgG-Fc domain that strengthen hexamerization strongly enhance C1 recruitment and downstream complement activation on encapsulated pneumococci. Specifically, hexamerization-enhancing mutations E430G or E345K in CPS6-IgG strongly potentiate complement activation on S. pneumoniae strains that express capsular serotype 6 (CPS6), and the highly invasive serotype 19A strain. Furthermore, these mutations improve complement activation via mAbs recognizing CPS3 and CPS8 strains. Importantly, hexamer-enhancing mutations enable mAbs to induce strong opsonophagocytic killing by human neutrophils. Finally, passive immunization with CPS6-IgG1-E345K protected mice from developing severe pneumonia. Altogether, this work provides an important proof of concept for future optimization of antibody therapies against encapsulated bacteria.

High-resolution quantitative and functional MRI indicate lower myelination of thin and thick stripes in human secondary visual cortex

The characterization of cortical myelination is essential for the study of structure-function relationships in the human brain. However, knowledge about cortical myelination is largely based on post-mortem histology, which generally renders direct comparison to function impossible. The repeating pattern of pale-thin-pale-thick stripes of cytochrome oxidase (CO) activity in the primate secondary visual cortex (V2) is a prominent columnar system, in which histology also indicates different myelination of thin/thick versus pale stripes. We used quantitative magnetic resonance imaging (qMRI) in conjunction with functional magnetic resonance imaging (fMRI) at ultra-high field strength (7 T) to localize and study myelination of stripes in four human participants at sub-millimeter resolution in vivo. Thin and thick stripes were functionally localized by exploiting their sensitivity to color and binocular disparity, respectively. Resulting functional activation maps showed robust stripe patterns in V2 which enabled further comparison of quantitative relaxation parameters between stripe types. Thereby, we found lower longitudinal relaxation rates (R1) of thin and thick stripes compared to surrounding gray matter in the order of 1-2%, indicating higher myelination of pale stripes. No consistent differences were found for effective transverse relaxation rates (R2*). The study demonstrates the feasibility to investigate structure-function relationships in living humans within one cortical area at the level of columnar systems using qMRI.

SNTA1 gene rescues ion channel function and is antiarrhythmic in cardiomyocytes derived from induced pluripotent stem cells from muscular dystrophy patients

Background

Patients with cardiomyopathy of Duchenne Muscular Dystrophy (DMD) are at risk of developing life-threatening arrhythmias, but the mechanisms are unknown. We aimed to determine the role of ion channels controlling cardiac excitability in the mechanisms of arrhythmias in DMD patients.

Methods

To test whether dystrophin mutations lead to defective cardiac NaV1.5-Kir2.1 channelosomes and arrhythmias, we generated iPSC-CMs from two hemizygous DMD males, a heterozygous female, and two unrelated control males. We conducted studies including confocal microscopy, protein expression analysis, patch-clamping, non-viral piggy-bac gene expression, optical mapping and contractility assays.

Results

Two patients had abnormal ECGs with frequent runs of ventricular tachycardia. iPSC-CMs from all DMD patients showed abnormal action potential profiles, slowed conduction velocities, and reduced sodium (INa) and inward rectifier potassium (IK1) currents. Membrane NaV1.5 and Kir2.1 protein levels were reduced in hemizygous DMD iPSC-CMs but not in heterozygous iPSC-CMs. Remarkably, transfecting just one component of the dystrophin protein complex (α1-syntrophin) in hemizygous iPSC-CMs from one patient restored channelosome function, INa and IK1 densities, and action potential profile in single cells. In addition, α1-syntrophin expression restored impulse conduction and contractility and prevented reentrant arrhythmias in hiPSC-CM monolayers.

Conclusions

We provide the first demonstration that iPSC-CMs reprogrammed from skin fibroblasts of DMD patients with cardiomyopathy have a dysfunction of the NaV1.5-Kir2.1 channelosome, with consequent reduction of cardiac excitability and conduction. Altogether, iPSC-CMs from patients with DMD cardiomyopathy have a NaV1.5-Kir2.1 channelosome dysfunction, which can be rescued by the scaffolding protein α1-syntrophin to restore excitability and prevent arrhythmias.

Funding

Supported by National Institutes of Health R01 HL122352 grant; 'la Caixa' Banking Foundation (HR18-00304); Fundación La Marató TV3: Ayudas a la investigación en enfermedades raras 2020 (LA MARATO-2020); Instituto de Salud Carlos III/FEDER/FSE; Horizon 2020 - Research and Innovation Framework Programme GA-965286 to JJ; the CNIC is supported by the Instituto de Salud Carlos III (ISCIII), the Ministerio de Ciencia e Innovación (MCIN) and the Pro CNIC Foundation), and is a Severo Ochoa Center of Excellence (grant CEX2020-001041-S funded by MICIN/AEI/10.13039/501100011033). American Heart Association postdoctoral fellowship 19POST34380706s to JVEN. Israel Science Foundation to OB and MA [824/19]. Rappaport grant [01012020RI]; and Niedersachsen Foundation [ZN3452] to OB; US-Israel Binational Science Foundation (BSF) to OB and TH [2019039]; Dr. Bernard Lublin Donation to OB; and The Duchenne Parent Project Netherlands (DPPNL 2029771) to OB. National Institutes of Health R01 AR068428 to DM and US-Israel Binational Science Foundation Grant [2013032] to DM and OB.

Genome-scale modeling of yeast metabolism: retrospectives and perspectives

Yeasts have been widely used for production of bread, beer and wine, as well as for production of bioethanol, but they have also been designed as cell factories to produce various chemicals, advanced biofuels and recombinant proteins. To systematically understand and rationally engineer yeast metabolism, genome-scale metabolic models (GEMs) have been reconstructed for the model yeast Saccharomyces cerevisiae and nonconventional yeasts. Here, we review the historical development of yeast GEMs together with their recent applications, including metabolic flux prediction, cell factory design, culture condition optimization and multi-yeast comparative analysis. Furthermore, we present an emerging effort, namely the integration of proteome constraints into yeast GEMs, resulting in models with improved performance. At last, we discuss challenges and perspectives on the development of yeast GEMs and the integration of proteome constraints.

Rotational dynamics in motor cortex are consistent with a feedback controller

Recent studies have identified rotational dynamics in motor cortex (MC), which many assume arise from intrinsic connections in MC. However, behavioral and neurophysiological studies suggest that MC behaves like a feedback controller where continuous sensory feedback and interactions with other brain areas contribute substantially to MC processing. We investigated these apparently conflicting theories by building recurrent neural networks that controlled a model arm and received sensory feedback from the limb. Networks were trained to counteract perturbations to the limb and to reach toward spatial targets. Network activities and sensory feedback signals to the network exhibited rotational structure even when the recurrent connections were removed. Furthermore, neural recordings in monkeys performing similar tasks also exhibited rotational structure not only in MC but also in somatosensory cortex. Our results argue that rotational structure may also reflect dynamics throughout the voluntary motor system involved in online control of motor actions.

Coupling between motor cortex and striatum increases during sleep over long-term skill learning

The strength of cortical connectivity to the striatum influences the balance between behavioral variability and stability. Learning to consistently produce a skilled action requires plasticity in corticostriatal connectivity associated with repeated training of the action. However, it remains unknown whether such corticostriatal plasticity occurs during training itself or 'offline' during time away from training, such as sleep. Here, we monitor the corticostriatal network throughout long-term skill learning in rats and find that non-rapid-eye-movement (NREM) sleep is a relevant period for corticostriatal plasticity. We first show that the offline activation of striatal NMDA receptors is required for skill learning. We then show that corticostriatal functional connectivity increases offline, coupled to emerging consistent skilled movements, and coupled cross-area neural dynamics. We then identify NREM sleep spindles as uniquely poised to mediate corticostriatal plasticity, through interactions with slow oscillations. Our results provide evidence that sleep shapes cross-area coupling required for skill learning.

Highly contiguous assemblies of 101 drosophilid genomes

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.

Spatiotemporal control of cell cycle acceleration during axolotl spinal cord regeneration

Axolotls are uniquely able to resolve spinal cord injuries, but little is known about the mechanisms underlying spinal cord regeneration. We previously found that tail amputation leads to reactivation of a developmental-like program in spinal cord ependymal cells (Rodrigo Albors et al., 2015), characterized by a high-proliferation zone emerging 4 days post-amputation (Rost et al., 2016). What underlies this spatiotemporal pattern of cell proliferation, however, remained unknown. Here, we use modeling, tightly linked to experimental data, to demonstrate that this regenerative response is consistent with a signal that recruits ependymal cells during ~85 hours after amputation within ~830 μm of the injury. We adapted Fluorescent Ubiquitination-based Cell Cycle Indicator (FUCCI) technology to axolotls (AxFUCCI) to visualize cell cycles in vivo. AxFUCCI axolotls confirmed the predicted appearance time and size of the injury-induced recruitment zone and revealed cell cycle synchrony between ependymal cells. Our modeling and imaging move us closer to understanding bona fide spinal cord regeneration.

A statistical framework for assessing pharmacological responses and biomarkers using uncertainty estimates

High-throughput testing of drugs across molecular-characterised cell lines can identify candidate treatments and discover biomarkers. However, the cells' response to a drug is typically quantified by a summary statistic from a best-fit dose-response curve, whilst neglecting the uncertainty of the curve fit and the potential variability in the raw readouts. Here, we model the experimental variance using Gaussian Processes, and subsequently, leverage uncertainty estimates to identify associated biomarkers with a new Bayesian framework. Applied to in vitro screening data on 265 compounds across 1074 cancer cell lines, our models identified 24 clinically established drug-response biomarkers, and provided evidence for six novel biomarkers by accounting for association with low uncertainty. We validated our uncertainty estimates with an additional drug screen of 26 drugs, 10 cell lines with 8 to 9 replicates. Our method is applicable to any dose-response data without replicates, and improves biomarker discovery for precision medicine.

Influenza virus infection induces widespread alterations of host cell splicing

Influenza A viruses (IAVs) use diverse mechanisms to interfere with cellular gene expression. Although many RNA-seq studies have documented IAV-induced changes in host mRNA abundance, few were designed to allow an accurate quantification of changes in host mRNA splicing. Here, we show that IAV infection of human lung cells induces widespread alterations of cellular splicing, with an overall increase in exon inclusion and decrease in intron retention. Over half of the mRNAs that show differential splicing undergo no significant changes in abundance or in their 3' end termination site, suggesting that IAVs can specifically manipulate cellular splicing. Among a randomly selected subset of 21 IAV-sensitive alternative splicing events, most are specific to IAV infection as they are not observed upon infection with VSV, induction of interferon expression or induction of an osmotic stress. Finally, the analysis of splicing changes in RED-depleted cells reveals a limited but significant overlap with the splicing changes in IAV-infected cells. This observation suggests that hijacking of RED by IAVs to promote splicing of the abundant viral NS1 mRNAs could partially divert RED from its target mRNAs. All our RNA-seq datasets and analyses are made accessible for browsing through a user-friendly Shiny interface (http://virhostnet.prabi.fr:3838/shinyapps/flu-splicing or https://github.com/cbenoitp/flu-splicing).

Vascular control of the CO2/H+-dependent drive to breathe

Respiratory chemoreceptors regulate breathing in response to changes in tissue CO2/H+. Blood flow is a fundamental determinant of tissue CO2/H+, yet little is known regarding how regulation of vascular tone in chemoreceptor regions contributes to respiratory behavior. Previously, we showed in rat that CO2/H+-vasoconstriction in the retrotrapezoid nucleus (RTN) supports chemoreception by a purinergic-dependent mechanism (Hawkins et al., 2017). Here, we show in mice that CO2/H+ dilates arterioles in other chemoreceptor regions, thus demonstrating CO2/H+ vascular reactivity in the RTN is unique. We also identify P2Y2 receptors in RTN smooth muscle cells as the substrate responsible for this response. Specifically, pharmacological blockade or genetic deletion of P2Y2 from smooth muscle cells blunted the ventilatory response to CO2, and re-expression of P2Y2 receptors only in RTN smooth muscle cells fully rescued the CO2/H+ chemoreflex. These results identify P2Y2 receptors in RTN smooth muscle cells as requisite determinants of respiratory chemoreception.

Anatomy of nerve fiber bundles at micrometer-resolution in the vervet monkey visual system

Although the primate visual system has been extensively studied, detailed spatial organization of white matter fiber tracts carrying visual information between areas has not been fully established. This is mainly due to the large gap between tracer studies and diffusion-weighted MRI studies, which focus on specific axonal connections and macroscale organization of fiber tracts, respectively. Here we used 3D polarization light imaging (3D-PLI), which enables direct visualization of fiber tracts at micrometer resolution, to identify and visualize fiber tracts of the visual system, such as stratum sagittale, inferior longitudinal fascicle, vertical occipital fascicle, tapetum and dorsal occipital bundle in vervet monkey brains. Moreover, 3D-PLI data provide detailed information on cortical projections of these tracts, distinction between neighboring tracts, and novel short-range pathways. This work provides essential information for interpretation of functional and diffusion-weighted MRI data, as well as revision of wiring diagrams based upon observations in the vervet visual system.

FUS ALS-causative mutations impair FUS autoregulation and splicing factor networks through intron retention

Mutations in the RNA-binding protein FUS cause amyotrophic lateral sclerosis (ALS), a devastating neurodegenerative disease. FUS plays a role in numerous aspects of RNA metabolism, including mRNA splicing. However, the impact of ALS-causative mutations on splicing has not been fully characterized, as most disease models have been based on overexpressing mutant FUS, which will alter RNA processing due to FUS autoregulation. We and others have recently created knockin models that overcome the overexpression problem, and have generated high depth RNA-sequencing on FUS mutants in parallel to FUS knockout, allowing us to compare mutation-induced changes to genuine loss of function. We find that FUS-ALS mutations induce a widespread loss of function on expression and splicing. Specifically, we find that mutant FUS directly alters intron retention levels in RNA-binding proteins. Moreover, we identify an intron retention event in FUS itself that is associated with its autoregulation. Altered FUS levels have been linked to disease, and we show here that this novel autoregulation mechanism is altered by FUS mutations. Crucially, we also observe this phenomenon in other genetic forms of ALS, including those caused by TDP-43, VCP and SOD1 mutations, supporting the concept that multiple ALS genes interact in a regulatory network.

Numerical magnitude, rather than individual bias, explains spatial numerical association in newborn chicks

We associate small numbers with the left and large numbers with the right side of space. Recent evidence from human newborns and non-human animals has challenged the primary role assigned to culture, in determining this spatial numerical association (SNA). Nevertheless, the effect of individual spatial biases has not been considered in previous research. Here, we tested the effect of numerical magnitude in SNA and we controlled for itablendividual biases. We trained 3-day-old chicks (Gallus gallus) on a given numerical magnitude (5). Then chicks could choose between two identical, left or right, stimuli both representing either 2, 8, or 5 elements. We computed the percentage of Left-sided Choice (LC). Numerical magnitude, but not individual lateral bias, explained LC: LC2 vs. 2>LC5 vs. 5>LC8 vs. 8. These findings suggest that SNA originates from pre-linguistic precursors, and pave the way to the investigation of the neural correlates of the number space association.

Tex19.1 inhibits the N-end rule pathway and maintains acetylated SMC3 cohesin and sister chromatid cohesion in oocytes

Age-dependent oocyte aneuploidy, a major cause of Down syndrome, is associated with declining sister chromatid cohesion in postnatal oocytes. Here we show that cohesion in postnatal mouse oocytes is regulated by Tex19.1. We show Tex19.1-/- oocytes have defects maintaining chiasmata, missegregate their chromosomes during meiosis, and transmit aneuploidies to the next generation. Furthermore, we show that mouse Tex19.1 inhibits N-end rule protein degradation mediated by its interacting partner UBR2, and that Ubr2 itself has a previously undescribed role in negatively regulating the acetylated SMC3 subpopulation of cohesin in mitotic somatic cells. Lastly, we show that acetylated SMC3 is associated with meiotic chromosome axes in mouse oocytes, and that this population of cohesin is specifically depleted in the absence of Tex19.1. These findings indicate that Tex19.1 regulates UBR protein activity to maintain acetylated SMC3 and sister chromatid cohesion in postnatal oocytes and prevent aneuploidy from arising in the female germline.

Sex Classification by Resting State Brain Connectivity

A large amount of brain imaging research has focused on group studies delineating differences between males and females with respect to both cognitive performance as well as structural and functional brain organization. To supplement existing findings, the present study employed a machine learning approach to assess how accurately participants' sex can be classified based on spatially specific resting state (RS) brain connectivity, using 2 samples from the Human Connectome Project (n1 = 434, n2 = 310) and 1 fully independent sample from the 1000BRAINS study (n = 941). The classifier, which was trained on 1 sample and tested on the other 2, was able to reliably classify sex, both within sample and across independent samples, differing both with respect to imaging parameters and sample characteristics. Brain regions displaying highest sex classification accuracies were mainly located along the cingulate cortex, medial and lateral frontal cortex, temporoparietal regions, insula, and precuneus. These areas were stable across samples and match well with previously described sex differences in functional brain organization. While our data show a clear link between sex and regionally specific brain connectivity, they do not support a clear-cut dimorphism in functional brain organization that is driven by sex alone.

Novel genetic loci affecting facial shape variation in humans

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.