Knowing how the human brain is shaped by migration and admixture is a critical step in studying human evolution [1, 2], as well as in preventing the bias of hidden population structure in brain research [3, 4]. Yet, the neuroanatomical differences engendered by population history are still poorly understood. Most of the inference relies on craniometric measurements, because morphology of the brain is presumed to be the neurocranium's main shaping force before bones are fused and ossified [5]. Although studies have shown that the shape variations of cranial bones are consistent with population history [6-8], it is unknown how much human ancestry information is retained by the human cortical surface. In our group's previous study, we found that area measures of cortical surface and total brain volumes of individuals of European descent in the United States correlate significantly with their ancestral geographic locations in Europe [9]. Here, we demonstrate that the three-dimensional geometry of cortical surface is highly predictive of individuals' genetic ancestry in West Africa, Europe, East Asia, and America, even though their genetic background has been shaped by multiple waves of migratory and admixture events. The geometry of the cortical surface contains richer information about ancestry than the areal variability of the cortical surface, independent of total brain volumes. Besides explaining more ancestry variance than other brain imaging measurements, the 3D geometry of the cortical surface further characterizes distinct regional patterns in the folding and gyrification of the human brain associated with each ancestral lineage.

Little is known about how genetic variation contributes to neuroanatomical variability, and whether particular genomic regions comprising genes or evolutionarily conserved elements are enriched for effects that influence brain morphology. Here, we examine brain imaging and single-nucleotide polymorphisms (SNPs) data from ∼2,700 individuals. We show that a substantial proportion of variation in cortical surface area is explained by additive effects of SNPs dispersed throughout the genome, with a larger heritable effect for visual and auditory sensory and insular cortices (h(2)∼0.45). Genome-wide SNPs collectively account for, on average, about half of twin heritability across cortical regions (N=466 twins). We find enriched genetic effects in or near genes. We also observe that SNPs in evolutionarily more conserved regions contributed significantly to the heritability of cortical surface area, particularly, for medial and temporal cortical regions. SNPs in less conserved regions contributed more to occipital and dorsolateral prefrontal cortices.

Additive genetic variance (VA ) and total genetic variance (VG ) are core concepts in biomedical, evolutionary and production-biology genetics. What determines the large variation in reported VA /VG ratios from line-cross experiments is not well understood. Here we report how the VA /VG ratio, and thus the ratio between narrow and broad sense heritability (h(2) /H(2) ), varies as a function of the regulatory architecture underlying genotype-to-phenotype (GP) maps. We studied five dynamic models (of the cAMP pathway, the glycolysis, the circadian rhythms, the cell cycle, and heart cell dynamics). We assumed genetic variation to be reflected in model parameters and extracted phenotypes summarizing the system dynamics. Even when imposing purely linear genotype to parameter maps and no environmental variation, we observed quite low VA /VG ratios. In particular, systems with positive feedback and cyclic dynamics gave more non-monotone genotype-phenotype maps and much lower VA /VG ratios than those without. The results show that some regulatory architectures consistently maintain a transparent genotype-to-phenotype relationship, whereas other architectures generate more subtle patterns. Our approach can be used to elucidate these relationships across a whole range of biological systems in a systematic fashion.

Chronic abuse of heroin leads to long-lasting and complicated cognitive impairment. Dopamine receptors are critically involved in the impulsive drug-driven behavior and the altered attention, processing speed, and mental flexibility that are associated with higher relapse rates. However, the effects of the different dopamine receptors and their possible involvement in heroin-induced cognitive impairment remain unclear.

Bone Mineral Density (BMD) is a highly heritable trait, but genome-wide association studies have identified few genetic risk factors. Epidemiological studies suggest associations between BMD and several traits and diseases, but the nature of the suggestive comorbidity is still unknown. We used a novel genetic pleiotropy-informed conditional False Discovery Rate (FDR) method to identify single nucleotide polymorphisms (SNPs) associated with BMD by leveraging cardiovascular disease (CVD) associated disorders and metabolic traits. By conditioning on SNPs associated with the CVD-related phenotypes, type 1 diabetes, type 2 diabetes, systolic blood pressure, diastolic blood pressure, high density lipoprotein, low density lipoprotein, triglycerides and waist hip ratio, we identified 65 novel independent BMD loci (26 with femoral neck BMD and 47 with lumbar spine BMD) at conditional FDR < 0.01. Many of the loci were confirmed in genetic expression studies. Genes validated at the mRNA levels were characteristic for the osteoblast/osteocyte lineage, Wnt signaling pathway and bone metabolism. The results provide new insight into genetic mechanisms of variability in BMD, and a better understanding of the genetic underpinnings of clinical comorbidity.

Most of the genetic architecture of schizophrenia (SCZ) has not yet been identified. Here, we apply a novel statistical algorithm called Covariate-Modulated Mixture Modeling (CM3), which incorporates auxiliary information (heterozygosity, total linkage disequilibrium, genomic annotations, pleiotropy) for each single nucleotide polymorphism (SNP) to enable more accurate estimation of replication probabilities, conditional on the observed test statistic ("z-score") of the SNP. We use a multiple logistic regression on z-scores to combine information from auxiliary information to derive a "relative enrichment score" for each SNP. For each stratum of these relative enrichment scores, we obtain nonparametric estimates of posterior expected test statistics and replication probabilities as a function of discovery z-scores, using a resampling-based approach that repeatedly and randomly partitions meta-analysis sub-studies into training and replication samples. We fit a scale mixture of two Gaussians model to each stratum, obtaining parameter estimates that minimize the sum of squared differences of the scale-mixture model with the stratified nonparametric estimates. We apply this approach to the recent genome-wide association study (GWAS) of SCZ (n = 82,315), obtaining a good fit between the model-based and observed effect sizes and replication probabilities. We observed that SNPs with low enrichment scores replicate with a lower probability than SNPs with high enrichment scores even when both they are genome-wide significant (p < 5x10-8). There were 693 and 219 independent loci with model-based replication rates ≥80% and ≥90%, respectively. Compared to analyses not incorporating relative enrichment scores, CM3 increased out-of-sample yield for SNPs that replicate at a given rate. This demonstrates that replication probabilities can be more accurately estimated using prior enrichment information with CM3.

Recently, the insular cortex (IC) was identified as part of the neuronal circuit responsible for the reward expectations in cue-triggered behaviours. Moreover, there are evidences that connections between the IC and the ventral striatum, particularly with the nucleus accumbens (NAc), may mediate the retrieval and performance of actions based on incentive memory. However, the precise role of the IC-NAc connections in cue-related drug-seeking behaviour remains unclear. We used the morphine-induced conditioned place preference (CPP) paradigm to assess the formation and relapse of cue-related drug-seeking. cFos immunostaining was used to determine the activation of the brain regions. Chemogenetic and optogenetic methods were used to manipulate the activity of IC-to-NAc projection neurons. The result showed that neurons in IC and NAc core but not NAc shell were activated following cue-induced morphine-seeking behaviour. Negligible effect of inhibition of IC-to-NAc core projection (IC→NAc core) on morphine CPP expression, whereas chemogenetic inactivation of this projection potently blocked the reinstatement of expressed morphine CPP. Furthermore, optogenetic inhibition of glutamatergic IC→NAc core inputs significant suppressed the CPP reinstatement without significant effect on CPP expression. We demonstrated here, for the first time, that IC→NAc core glutamatergic projection is required for the reinstatement of cue-associated drug seeking behaviour in mice. The present study provide insights into modulations of relapse of cue-associated drug-seeking behaviour following repeated overexposure to opioids in humans.

Major depressive disorder (MDD) is a common illness accompanied by considerable morbidity, mortality, costs, and heightened risk of suicide. We conducted a genome-wide association meta-analysis based in 135,458 cases and 344,901 controls and identified 44 independent and significant loci. The genetic findings were associated with clinical features of major depression and implicated brain regions exhibiting anatomical differences in cases. Targets of antidepressant medications and genes involved in gene splicing were enriched for smaller association signal. We found important relationships of genetic risk for major depression with educational attainment, body mass, and schizophrenia: lower educational attainment and higher body mass were putatively causal, whereas major depression and schizophrenia reflected a partly shared biological etiology. All humans carry lesser or greater numbers of genetic risk factors for major depression. These findings help refine the basis of major depression and imply that a continuous measure of risk underlies the clinical phenotype.

The risk of posttraumatic stress disorder (PTSD) following trauma is heritable, but robust common variants have yet to be identified. In a multi-ethnic cohort including over 30,000 PTSD cases and 170,000 controls we conduct a genome-wide association study of PTSD. We demonstrate SNP-based heritability estimates of 5-20%, varying by sex. Three genome-wide significant loci are identified, 2 in European and 1 in African-ancestry analyses. Analyses stratified by sex implicate 3 additional loci in men. Along with other novel genes and non-coding RNAs, a Parkinson's disease gene involved in dopamine regulation, PARK2, is associated with PTSD. Finally, we demonstrate that polygenic risk for PTSD is significantly predictive of re-experiencing symptoms in the Million Veteran Program dataset, although specific loci did not replicate. These results demonstrate the role of genetic variation in the biology of risk for PTSD and highlight the necessity of conducting sex-stratified analyses and expanding GWAS beyond European ancestry populations.

Discovering genetic variants associated with human brain structures is an on-going effort. The ENIGMA consortium conducted genome-wide association studies (GWAS) with standard multi-study analytical methodology and identified several significant single nucleotide polymorphisms (SNPs). Here we employ a novel analytical approach that incorporates functional genome annotations (e.g., exon or 5'UTR), total linkage disequilibrium (LD) scores and heterozygosity to construct enrichment scores for improved identification of relevant SNPs. The method provides increased power to detect associated SNPs by estimating stratum-specific false discovery rate (FDR), where strata are classified according to enrichment scores. Applying this approach to the GWAS summary statistics of putamen volume in the ENIGMA cohort, a total of 15 independent significant SNPs were identified (conditional FDR < 0.05). In contrast, 4 SNPs were found based on standard GWAS analysis (P < 5 × 10-8). These 11 novel loci include GATAD2B, ASCC3, DSCAML1, and HELZ, which are previously implicated in various neural related phenotypes. The current findings demonstrate the boost in power with the annotation-informed FDR method, and provide insight into the genetic architecture of the putamen.

Primary open-angle glaucoma (POAG) is a leading cause of irreversible blindness worldwide, and elevated intraocular pressure (IOP) is a major risk factor. While IOP is mainly controlled by adjusting the outflow resistance in the trabecular meshwork (TM), drugs that act directly on the TM are rare. In this study, we discovered a novel compound and pathway that acts on the TM and decreases IOP by genomic, proteomic, and bioinformatic analyses of POAG-derived TMs and experimental validation. Overlapping differentially expressed genes of the TM between patients with POAG and normal controls from two independent gene expression profiles in public databases were analyzed and matched by using the Connectivity Map (CMap). Rottlerin was identified as a potential compound. Subsequent experiments confirmed that rottlerin reversed POAG phenotypes in vitro and that it decreased IOP and actin/extracellular matrix accumulation in vivo with no detectable ocular side effects. SwissTargetPrediction in combination with pathway analysis predicted that the effects of rottlerin may be mediated by activation of the Rap1 pathway. Finally, we confirmed that rottlerin upregulated Rap1 and the downstream PI3K/AKT pathway independent of the MAPK/ERK pathway in a dexamethasone-induced POAG cell model.

Higher cognitive functions are regarded as one of the main distinctive traits of humans. Evidence for the cognitive evolution of human beings is mainly based on fossil records of an expanding cranium and an increasing complexity of material culture artefacts. However, the molecular genetic factors involved in the evolution are still relatively unexplored. Here, we investigated whether genomic regions that underwent positive selection in humans after divergence from Neanderthals are enriched for genetic association with phenotypes related to cognitive functions. We used genome wide association data from a study of college completion (N = 111,114), one of educational attainment (N = 293,623) and two different studies of general cognitive ability (N = 269,867 and 53,949). We found nominally significant polygenic enrichment of associations with college completion (p = 0.025), educational attainment (p = 0.043) and general cognitive ability (p = 0.015 and 0.025, respectively), suggesting that variants influencing these phenotypes are more prevalent in evolutionarily salient regions. The enrichment remained significant after controlling for other known genetic enrichment factors, and for affiliation to genes highly expressed in the brain. These findings support the notion that phenotypes related to higher order cognitive skills typical of humans have a recent genetic component that originated after the separation of the human and Neanderthal lineages.

Development and aging of the cerebral cortex show similar topographic organization and are governed by the same genes. It is unclear whether the same is true for subcortical regions, which follow fundamentally different ontogenetic and phylogenetic principles. We tested the hypothesis that genetically governed neurodevelopmental processes can be traced throughout life by assessing to which degree brain regions that develop together continue to change together through life. Analyzing over 6000 longitudinal MRIs of the brain, we used graph theory to identify five clusters of coordinated development, indexed as patterns of correlated volumetric change in brain structures. The clusters tended to follow placement along the cranial axis in embryonic brain development, suggesting continuity from prenatal stages, and correlated with cognition. Across independent longitudinal datasets, we demonstrated that developmental clusters were conserved through life. Twin-based genetic correlations revealed distinct sets of genes governing change in each cluster. Single-nucleotide polymorphisms-based analyses of 38,127 cross-sectional MRIs showed a similar pattern of genetic volume-volume correlations. In conclusion, coordination of subcortical change adheres to fundamental principles of lifespan continuity and genetic organization.

Early metastasis of pancreatic cancer (PC) leads to high mortality, and the underlying mechanism of metastasis remains unclear. Tumor necrosis factor superfamily member 9 (TNFSF9) is associated with poor prognosis in PC. Here, we investigated the effect of TNFSF9 on PC proliferation and apoptosis, and focused on the effect of TNFSF9 on PC metastasis and its potential mechanism. We found that TNFSF9 promotes PC metastasis in vivo and in vitro, and may be partially dependent on the Wnt/Snail signaling pathway. In addition, TNFSF9 also regulates the release of cytokines IL-10 and transforming growth factor-β (TGF-β) in pancreatic cancer cells through Wnt signaling to induce the M2 polarization of macrophages and promote the migration of PC cells. Overall, our study found that TNFSF9 may directly promote PC metastasis or indirectly promote PC metastasis through macrophage M2 polarization. Our study provides a new costimulatory target for the treatment of PC.

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

A promise of genomics in precision medicine is to provide individualized genetic risk predictions. Polygenic risk scores (PRS), computed by aggregating effects from many genomic variants, have been developed as a useful tool in complex disease research. However, the application of PRS as a tool for predicting an individual's disease susceptibility in a clinical setting is challenging because PRS typically provide a relative measure of risk evaluated at the level of a group of people but not at individual level. Here, we introduce a machine-learning technique, Mondrian Cross-Conformal Prediction (MCCP), to estimate the confidence bounds of PRS-to-disease-risk prediction. MCCP can report disease status conditional probability value for each individual and give a prediction at a desired error level. Moreover, with a user-defined prediction error rate, MCCP can estimate the proportion of sample (coverage) with a correct prediction.

The nucleus accumbens (NAc) is the most promising target for drug use disorder treatment. Deep brain stimulation (DBS) of NAc is effective for drug use disorder treatment. However, the mechanisms by which DBS produces its therapeutic effects remain enigmatic. Here, we define a behavioral cutoff criterion to distinguish depressive-like behaviors and non-depressive-like behaviors in mice after morphine withdrawal. We identified a basolateral amygdala (BLA) to NAc D1 medium spiny neuron (MSN) pathway that controls depressive-like behaviors after morphine withdrawal. Furthermore, the paraventricular nucleus of thalamus (PVT) to NAc D2 MSN pathway controls naloxone-induced acute withdrawal symptoms. Optogenetically induced long-term potentiation with κ-opioid receptor (KOR) antagonism enhanced BLA to NAc D1 MSN signaling and also altered the excitation/inhibition balance of NAc D2 MSN signaling. We also verified that a new 50 Hz DBS protocol reversed morphine withdrawal-evoked abnormal plasticity in NAc. Importantly, this refined DBS treatment effectively alleviated naloxone-induced withdrawal symptoms and depressive-like behaviors and prevented stress-induced reinstatement. Taken together, the results demonstrated that input- and cell type-specific synaptic plasticity underlies morphine withdrawal, which may lead to novel targets for the treatment of opioid use disorder.

To restrain the chemical reaction at cathode interface of organic solar cells, two cathode interfacial materials are synthesized by connecting phenanthroline with carbolong unit. Consequently, the D18:L8-BO based organic solar cell with double-phenanthroline-carbolong achieve the highest efficiency of 18.2%. Double-phenanthroline-carbolong with larger steric hindrance and stronger electron-withdrawing property confirms to suppress the interfacial reaction with norfullerene acceptor, resulting the most stable device. Double-phenanthroline-carbolong based device can sustain 80% of its initial efficiency for 2170 h in dark N2 atmosphere, 96 h under 85 oC and keep 68% initial efficiency after been illuminated for 2200 h, which are significantly better than bathocuproin based devices. Moreover, superb interfacial stability of double-phenanthroline-carbolong cathode interface enables thermal posttreatment of organic sub-cell in perovskite/organic tandem solar cells and obtained a remarkable efficiency of 21.7% with excellent thermal stability, which indicates the potentially wide application of phenanthroline-carbolong materials for stable and efficient solar device fabrications.

Cortical asymmetry is a ubiquitous feature of brain organization that is subtly altered in some neurodevelopmental disorders, yet we lack knowledge of how its development proceeds across life in health. Achieving consensus on the precise cortical asymmetries in humans is necessary to uncover the developmental timing of asymmetry and the extent to which it arises through genetic and later influences in childhood. Here, we delineate population-level asymmetry in cortical thickness and surface area vertex-wise in seven datasets and chart asymmetry trajectories longitudinally across life (4-89 years; observations = 3937; 70% longitudinal). We find replicable asymmetry interrelationships, heritability maps, and test asymmetry associations in large-scale data. Cortical asymmetry was robust across datasets. Whereas areal asymmetry is predominantly stable across life, thickness asymmetry grows in childhood and peaks in early adulthood. Areal asymmetry is low-moderately heritable (max h2SNP ~19%) and correlates phenotypically and genetically in specific regions, indicating coordinated development of asymmetries partly through genes. In contrast, thickness asymmetry is globally interrelated across the cortex in a pattern suggesting highly left-lateralized individuals tend towards left-lateralization also in population-level right-asymmetric regions (and vice versa), and exhibits low or absent heritability. We find less areal asymmetry in the most consistently lateralized region in humans associates with subtly lower cognitive ability, and confirm small handedness and sex effects. Results suggest areal asymmetry is developmentally stable and arises early in life through genetic but mainly subject-specific stochastic effects, whereas childhood developmental growth shapes thickness asymmetry and may lead to directional variability of global thickness lateralization in the population.

The question of whether immune dysfunction contributes to risk of psychiatric disorders has long been a subject of interest. To assert this hypothesis a plethora of correlative evidence has been accumulated from the past decades; however, a variety of technical and practical obstacles impeded on a cause-effect interpretation of these data. With the advent of large-scale omics technology and advanced statistical models, particularly Mendelian randomization, new studies testing this old hypothesis are accruing. Here we synthesize these new findings from genomics and genetic causal inference studies on the role of immune dysfunction in major psychiatric disorders and reconcile these new data with pre-omics findings. By reconciling these evidences, we aim to identify key gaps and propose directions for future studies in the field.