ComBatLS: A location- and scale-preserving method for multi-site image harmonization.

Recent work has leveraged massive datasets and advanced harmonization methods to construct normative models of neuroanatomical features and benchmark individuals' morphology. However, current harmonization tools do not preserve the effects of biological covariates including sex and age on features' variances; this failure may induce error in normative scores, particularly when such factors are distributed unequally across sites. Here, we introduce a new extension of the popular ComBat harmonization method, ComBatLS, that preserves biological variance in features' locations and scales. We use UK Biobank data to show that ComBatLS robustly replicates individuals' normative scores better than other ComBat methods when subjects are assigned to sex-imbalanced synthetic "sites". Additionally, we demonstrate that ComBatLS significantly reduces sex biases in normative scores compared to traditional methods. Finally, we show that ComBatLS successfully harmonizes consortium data collected across over 50 studies. R implementation of ComBatLS is available at https://github.com/andy1764/ComBatFamily.

An axis of genetic heterogeneity in autism is indexed by age at diagnosis and is associated with varying developmental and mental health profiles.

There is growing recognition that earliest signs of autism need not clearly manifest in the first three years of life. To what extent is this variation in developmental trajectories associated with age at autism diagnosis? Does the genetic profile of autism vary with age at autism diagnosis? Using longitudinal data from four birth cohorts, we demonstrate that two different trajectories of socio-emotional behaviours are associated with age at diagnosis. We further demonstrate that the age at autism diagnosis is partly heritable (h2 SNP = 0.12, s.e.m = 0.01), and is associated with two moderately correlated (rg = 0.38, s.e.m = 0.07) autism polygenic factors. One of these factors is associated with earlier diagnosis of autism, lower social and communication abilities in early childhood. The second factor is associated with later autism diagnosis, increased socio-emotional difficulties in adolescence, and has moderate to high positive genetic correlations with Attention-Deficit/Hyperactivity Disorder, mental health conditions, and trauma. Overall, our research identifies an axis of heterogeneity in autism, indexed by age at diagnosis, which partly explains heterogeneity in autism and the profiles of co-occurring neurodevelopmental and mental health profiles. Our findings have important implications for how we conceptualise autism and provide one model to explain some of the diversity within autism.

Child maltreatment as a transdiagnostic risk factor for the externalizing dimension: a Mendelian randomization study.

Observational studies suggest that child maltreatment increases the risk of externalizing spectrum disorders such as attention deficit hyperactivity disorder (ADHD), conduct disorder (CD), antisocial personality disorder (ASPD), and substance use disorder (SUD). Yet, only few of such associations have been investigated by approaches that provide strong evidence for causation, such as Mendelian Randomization (MR). Establishing causal inference is essential given the growing recognition of gene-environment correlations, which can confound observational research in the context of childhood maltreatment. Evaluating causality between child maltreatment and the externalizing phenotypes, we used genome-wide association study (GWAS) summary data for child maltreatment (143,473 participants), ADHD (20,183 cases; 35,191 controls), CD (451 cases; 256,859 controls), ASPD (381 cases; 252,877 controls), alcohol use disorder (AUD; 13,422 cases; 244,533 controls), opioid use disorder (OUD; 775 cases; 255,921 controls), and cannabinoid use disorder (CUD; 14,080 cases; 343,726 controls). We also generated a latent variable 'common externalizing factor' (EXT) using genomic structural equation modeling. Genetically predicted childhood maltreatment was consistently associated with ADHD (odds ratio [OR], 10.09; 95%-CI, 4.76-21.40; P = 1.63 × 10-09), AUD (OR, 3.72; 95%-CI, 1.85-7.52; P = 2.42 × 10-04), and the EXT (OR, 2.64; 95%-CI, 1.52-4.60; P = 5.80 × 10-04) across the different analyses and pleiotropy-robust methods. A subsequent GWAS on childhood maltreatment and the externalizing dimension from Externalizing Consortium (EXT-CON) confirmed these results. Two of the top five genes with the strongest associations in EXT GWAS, CADM2 and SEMA6D, are also ranked among the top 10 in the EXT-CON. The present results confirm the existence of a common externalizing factor and an increasing vulnerability caused by child maltreatment, with crucial implications for prevention. However, the partly diverging results also indicate that specific influences impact individual phenotypes separately.

Using Organoids to Model Sex Differences in the Human Brain.

Sex differences are widespread during neurodevelopment and play a role in neuropsychiatric conditions such as autism, which is more prevalent in males than females. In humans, males have been shown to have larger brain volumes than females with development of the hippocampus and amygdala showing prominent sex differences. Mechanistically, sex steroids and sex chromosomes drive these differences in brain development, which seem to peak during prenatal and pubertal stages. Animal models have played a crucial role in understanding sex differences, but the study of human sex differences requires an experimental model that can recapitulate complex genetic traits. To fill this gap, human induced pluripotent stem cell-derived brain organoids are now being used to study how complex genetic traits influence prenatal brain development. For example, brain organoids from individuals with autism and individuals with X chromosome-linked Rett syndrome and fragile X syndrome have revealed prenatal differences in cell proliferation, a measure of brain volume differences, and excitatory-inhibitory imbalances. Brain organoids have also revealed increased neurogenesis of excitatory neurons due to androgens. However, despite growing interest in using brain organoids, several key challenges remain that affect its validity as a model system. In this review, we discuss how sex steroids and the sex chromosomes each contribute to sex differences in brain development. Then, we examine the role of X chromosome inactivation as a factor that drives sex differences. Finally, we discuss the combined challenges of modeling X chromosome inactivation and limitations of brain organoids that need to be taken into consideration when studying sex differences.

Sex differences are a contributing factor in neuropsychiatric conditions such as autism, which is more prevalent in males. Sex differences occur through interactions between sex steroid hormones such as estrogen and testosterone and sex chromosomes (chrX and chrY). Human stem cell­derived brain organoids are laboratory models that mimic brain development. For example, in individuals with neurodevelopmental conditions, brain organoids have revealed an imbalance of neuron populations compared with neurotypical individuals. In this review, we discuss sex steroid and sex chromosome influences on brain development and challenges of this model that need to be taken into account when studying sex differences.

Integrating genomic variants and developmental milestones to predict cognitive and adaptive outcomes in autistic children.

Although the first signs of autism are often observed as early as 18-36 months of age, there is a broad uncertainty regarding future development, and clinicians lack predictive tools to identify those who will later be diagnosed with co-occurring intellectual disability (ID). Here, we developed predictive models of ID in autistic children (n=5,633 from three cohorts), integrating different classes of genetic variants alongside developmental milestones. The integrated model yielded an AUC ROC=0.65, with this predictive performance cross-validated and generalised across cohorts. Positive predictive values reached up to 55%, accurately identifying 10% of ID cases. The ability to stratify the probabilities of ID using genetic variants was up to twofold greater in individuals with delayed milestones compared to those with typical development. These findings underscore the potential of models in neurodevelopmental medicine that integrate genomics and clinical observations to predict outcomes and target interventions.

Classification of Suicide Attempt Risk Using Environmental and Lifestyle Factors in 3 Large Youth Cohorts.

Importance: Suicide is the third-leading cause of death among US adolescents. Environmental and lifestyle factors influence suicidal behavior and can inform risk classification, yet quantifying and incorporating them in risk assessment presents a significant challenge for reproducibility and clinical translation. Objective: To quantify the aggregate contribution of environmental and lifestyle factors to youth suicide attempt risk classification. Design, Setting, and Participants: This was a cohort study in 3 youth samples: 2 national longitudinal cohorts from the US and the UK and 1 clinical cohort from a tertiary pediatric US hospital. An exposome-wide association study (ExWAS) approach was used to identify risk and protective factors and compute aggregate exposomic scores. Logistic regression models were applied to test associations and model fit of exposomic scores with suicide attempts in independent data. Youth from the Adolescent Brain Cognitive Development (ABCD) study, the UK Millennium Cohort Study (MCS), and the Children's Hospital of Philadelphia emergency department (CHOP-ED) were included in the study. Exposures: A single-weighted exposomic score that sums significant risk and protective environmental/lifestyle factors. Main Outcome and Measure: Self-reported suicide attempt. Results: A total of 40 364 youth were included in this analysis: 11 564 from the ABCD study (3 waves of assessment; mean [SD] age, 12.0 [0.7] years; 6034 male [52.2%]; 344 attempted suicide [3.0%]; 1154 environmental/lifestyle factors were included in the ABCD study), 9000 from the MCS cohort (mean [SD] age, 17.2 [0.3] years; 4593 female [51.0%]; 661 attempted suicide [7.3%]; 2864 environmental/lifestyle factors were included in the MCS cohort), and 19 800 from the CHOP-ED cohort (mean [SD] age, 15.3 [1.5] years; 12 937 female [65.3%]; 2051 attempted suicide [10.4%]; 36 environmental/lifestyle factors were included in the CHOP-ED cohort). In the ABCD discovery subsample, ExWAS identified 99 risk and protective exposures significantly associated with suicide attempt. A single weighted exposomic score that sums significant risk and protective exposures was associated with suicide attempt in an independent ABCD testing subsample (odds ratio [OR], 2.2; 95% CI, 2.0-2.6; P < .001) and explained 17.6% of the variance (based on regression pseudo-R2) in suicide attempt over and above that explained by age, sex, race, and ethnicity (2.8%) and by family history of suicide (6.3%). Findings were consistent in the MCS and CHOP-ED cohorts (explaining 22.6% and 19.3% of the variance in suicide attempt, respectively) despite clinical, demographic, and exposure differences. In all cohorts, compared with youth at the median quintile of the exposomic score, youth at the top fifth quintile were substantially more likely to have made a suicide attempt (OR, 4.3; 95% CI, 2.6-7.2 in the ABCD study; OR, 3.8; 95% CI, 2.7-5.3 in the MCS cohort; OR, 5.8; 95% CI, 4.7-7.1 in the CHOP-ED cohort). Conclusions and Relevance: Results suggest that exposomic scores of suicide attempt provided a generalizable method for risk classification that can be applied in diverse samples from clinical or population settings.

Contribution of autosomal rare and de novo variants to sex differences in autism.

Autism is four times more prevalent in males than females. To study whether this reflects a difference in genetic predisposition attributed to autosomal rare variants, we evaluated the sex differences in effect size of damaging protein-truncating and missense variants on autism predisposition in 47,061 autistic individuals, then compared effect sizes between individuals with and without cognitive impairment or motor delay. Although these variants mediated differential likelihood of autism with versus without motor or cognitive impairment, their effect sizes on the liability scale did not differ significantly by sex exome-wide or in genes sex-differentially expressed in the cortex. Although de novo mutations were enriched in genes with male-biased expression in the fetal cortex, the liability they conferred did not differ significantly from other genes with similar loss-of-function intolerance and sex-averaged cortical expression. In summary, autosomal rare coding variants confer similar liability for autism in females and males.

Polygenic scores for autism are associated with neurite density in adults and children from the general population.

Genetic variants linked to autism are thought to change cognition and behaviour by altering the structure and function of the brain. Although a substantial body of literature has identified structural brain differences in autism, it is unknown whether autism-associated common genetic variants are linked to changes in cortical macro- and micro-structure. We investigated this using neuroimaging and genetic data from adults (UK Biobank, N = 31,748) and children (ABCD, N = 4,928). Using polygenic scores and genetic correlations we observe a robust negative association between common variants for autism and a magnetic resonance imaging derived phenotype for neurite density (intracellular volume fraction) in the general population. This result is consistent across both children and adults, in both the cortex and in white matter tracts, and confirmed using polygenic scores and genetic correlations. There were no sex differences in this association. Mendelian randomisation analyses provide no evidence for a causal relationship between autism and intracellular volume fraction, although this should be revisited using better powered instruments. Overall, this study provides evidence for shared common variant genetics between autism and cortical neurite density.

Genetics of gambling disorder and related phenotypes: The potential uses of polygenic and multifactorial risk models to enable early detection and improve clinical outcomes.

Gambling Disorder (GD) is an impactful behavioural addiction for which there appear to be underpinning genetic contributors. Twin studies show significant GD heritability results and intergenerational transmission show high rates of transmission. Recent developments in polygenic and multifactorial risk prediction modelling provide promising opportunities to enable early identification and intervention for at risk individuals. People with GD often have significant delays in diagnosis and subsequent help-seeking that can compromise their recovery. In this paper we advocate for more research into the utility of polygenic and multifactorial risk modelling in GD research and treatment programs and rigorous evaluation of its costs and benefits.

The genetic relationships between brain structure and schizophrenia.

Genetic risks for schizophrenia are theoretically mediated by genetic effects on brain structure but it has been unclear which genes are associated with both schizophrenia and cortical phenotypes. We accessed genome-wide association studies (GWAS) of schizophrenia (N = 69,369 cases; 236,642 controls), and of three magnetic resonance imaging (MRI) metrics (surface area, cortical thickness, neurite density index) measured at 180 cortical areas (N = 36,843, UK Biobank). Using Hi-C-coupled MAGMA, 61 genes were significantly associated with both schizophrenia and one or more MRI metrics. Whole genome analysis with partial least squares demonstrated significant genetic covariation between schizophrenia and area or thickness of most cortical regions. Genetic similarity between cortical areas was strongly coupled to their phenotypic covariance, and genetic covariation between schizophrenia and brain phenotypes was strongest in the hubs of structural covariance networks. Pleiotropically associated genes were enriched for neurodevelopmental processes and positionally concentrated in chromosomes 3p21, 17q21 and 11p11. Mendelian randomization analysis indicated that genetically determined variation in a posterior cingulate cortical area could be causal for schizophrenia. Parallel analyses of GWAS on bipolar disorder, Alzheimer's disease and height showed that pleiotropic association with MRI metrics was stronger for schizophrenia compared to other disorders.

Childhood abuse v. neglect and risk for major psychiatric disorders.

BACKGROUND: Childhood maltreatment (CM) is a strong risk factor for psychiatric disorders but serves in its current definitions as an umbrella for various fundamentally different childhood experiences. As first step toward a more refined analysis of the impact of CM, our objective is to revisit the relation of abuse and neglect, major subtypes of CM, with symptoms across disorders. METHODS: Three longitudinal studies of major depressive disorder (MDD, N = 1240), bipolar disorder (BD, N = 1339), and schizophrenia (SCZ, N = 577), each including controls (N = 881), were analyzed. Multivariate regression models were used to examine the relation between exposure to abuse, neglect, or their combination to the odds for MDD, BD, SCZ, and symptoms across disorders. Bidirectional Mendelian randomization (MR) was used to probe causality, using genetic instruments of abuse and neglect derived from UK Biobank data (N = 143 473). RESULTS: Abuse was the stronger risk factor for SCZ (OR 3.51, 95% CI 2.17-5.67) and neglect for BD (OR 2.69, 95% CI 2.09-3.46). Combined CM was related to increased risk exceeding additive effects of abuse and neglect for MDD (RERI = 1.4) and BD (RERI = 1.1). Across disorders, abuse was associated with hallucinations (OR 2.16, 95% CI 1.55-3.01) and suicide attempts (OR 2.16, 95% CI 1.55-3.01) whereas neglect was associated with agitation (OR 1.24, 95% CI 1.02-1.51) and reduced need for sleep (OR 1.64, 95% CI 1.08-2.48). MR analyses were consistent with a bidirectional causal effect of abuse with SCZ (IVWforward = 0.13, 95% CI 0.01-0.24). CONCLUSIONS: Childhood abuse and neglect are associated with different risks to psychiatric symptoms and disorders. Unraveling the origin of these differences may advance understanding of disease etiology and ultimately facilitate development of improved personalized treatment strategies.

The molecular genetic landscape of human brain size variation.

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

Two human brain systems micro-structurally associated with obesity.

The relationship between obesity and human brain structure is incompletely understood. Using diffusion-weighted MRI from ∼30,000 UK Biobank participants, we test the hypothesis that obesity (waist-to-hip ratio, WHR) is associated with regional differences in two micro-structural MRI metrics: isotropic volume fraction (ISOVF), an index of free water, and intra-cellular volume fraction (ICVF), an index of neurite density. We observed significant associations with obesity in two coupled but distinct brain systems: a prefrontal/temporal/striatal system associated with ISOVF and a medial temporal/occipital/striatal system associated with ICVF. The ISOVF~WHR system colocated with expression of genes enriched for innate immune functions, decreased glial density, and high mu opioid (MOR) and other neurotransmitter receptor density. Conversely, the ICVF~WHR system co-located with expression of genes enriched for G-protein coupled receptors and decreased density of MOR and other receptors. To test whether these distinct brain phenotypes might differ in terms of their underlying shared genetics or relationship to maps of the inflammatory marker C-reactive Protein (CRP), we estimated the genetic correlations between WHR and ISOVF (rg = 0.026, P = 0.36) and ICVF (rg = 0.112, P < 9×10-4) as well as comparing correlations between WHR maps and equivalent CRP maps for ISOVF and ICVF (P<0.05). These correlational results are consistent with a two-way mechanistic model whereby genetically determined differences in neurite density in the medial temporal system may contribute to obesity, whereas water content in the prefrontal system could reflect a consequence of obesity mediated by innate immune system activation.

People with obesity are at greater risk of cardiovascular diseases and metabolic conditions such as type 2 diabetes. More recently obesity has also been linked to changes in the brain that are associated with age-related dementia and cognitive decline. This includes a thinner cortex (the brain's outer layer) and lower volume of grey matter which is where cognitive processes, such as learning, take place. However, questions remain about how obesity and grey matter are connected. For instance, it is unclear whether the change in volume is due to there being fewer cells (and thus more water between them) or fewer connections between cells in these brain areas. It is also unknown whether the reduced volume of grey matter is a cause or consequence of obesity. To address these questions, Kitzbichler et al. analysed 30,000 MRI scans of the human brain which are stored in the UK Biobank. This revealed two characteristics in grey matter that were linked to obesity: higher amounts of water between cells in some areas, and a lower density of connections between neurons in others. The areas with higher levels of free water are known to have more glial cells which provide support to neurons. They also have more receptors that bind to fatty acids (which are often raised in people with obesity) and more receptors for molecules and cells involved in the immune response. In contrast, the areas with a lower density of connections between neurons usually were more closely associated with genetic risk factors associated with obesity, and fewer receptors involved in feeding, appetite and energy use. The findings of Kitzblicher et al. suggest that differences in the density of connections between neurons may contribute to obesity. High water content in grey matter, on the other hand, may be a consequence of obesity that occurs as a result of immune receptors becoming activated. This provides new insights in to how obesity and grey matter in the brain are connected.

The overlapping genetic architecture of psychiatric disorders and cortical brain structure.

Both psychiatric vulnerability and cortical structure are shaped by the cumulative effect of common genetic variants across the genome. However, the shared genetic underpinnings between psychiatric disorders and brain structural phenotypes, such as thickness and surface area of the cerebral cortex, remains elusive. In this study, we employed pleiotropy-informed conjunctional false discovery rate analysis to investigate shared loci across genome-wide association scans of regional cortical thickness, surface area, and seven psychiatric disorders in approximately 700,000 individuals of European ancestry. Aggregating regional measures, we identified 50 genetic loci shared between psychiatric disorders and surface area, as well as 26 genetic loci shared with cortical thickness. Risk alleles exhibited bidirectional effects on both cortical thickness and surface area, such that some risk alleles for each disorder increased regional brain size while other risk alleles decreased regional brain size. Due to bidirectional effects, in many cases we observed extensive pleiotropy between an imaging phenotype and a psychiatric disorder even in the absence of a significant genetic correlation between them. The impact of genetic risk for psychiatric disorders on regional brain structure did exhibit a consistent pattern across highly comorbid psychiatric disorders, with 80% of the genetic loci shared across multiple disorders displaying consistent directions of effect. Cortical patterning of genetic overlap revealed a hierarchical genetic architecture, with the association cortex and sensorimotor cortex representing two extremes of shared genetic influence on psychiatric disorders and brain structural variation. Integrating multi-scale functional annotations and transcriptomic profiles, we observed that shared genetic loci were enriched in active genomic regions, converged on neurobiological and metabolic pathways, and showed differential expression in postmortem brain tissue from individuals with psychiatric disorders. Cumulatively, these findings provide a significant advance in our understanding of the overlapping polygenic architecture between psychopathology and cortical brain structure.

Genetic insights into human cortical organization and development through genome-wide analyses of 2,347 neuroimaging phenotypes.

Our understanding of the genetics of the human cerebral cortex is limited both in terms of the diversity and the anatomical granularity of brain structural phenotypes. Here we conducted a genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging-derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,663 individuals and identified 4,349 experiment-wide significant loci. These phenotypes include cortical thickness, surface area, gray matter volume, measures of folding, neurite density and water diffusion. We identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with cortical expansion are associated with cephalic disorders. Finally, we identified complex interphenotype and inter-regional genetic relationships among the 13 phenotypes, reflecting the developmental differences among them. Together, these analyses identify distinct genetic organizational principles of the cortex and their correlates with neurodevelopment.

Robust estimation of cortical similarity networks from brain MRI.

Structural similarity is a growing focus for magnetic resonance imaging (MRI) of connectomes. Here we propose Morphometric INverse Divergence (MIND), a new method to estimate within-subject similarity between cortical areas based on the divergence between their multivariate distributions of multiple MRI features. Compared to the prior approach of morphometric similarity networks (MSNs) on n > 11,000 scans spanning three human datasets and one macaque dataset, MIND networks were more reliable, more consistent with cortical cytoarchitectonics and symmetry and more correlated with tract-tracing measures of axonal connectivity. MIND networks derived from human T1-weighted MRI were more sensitive to age-related changes than MSNs or networks derived by tractography of diffusion-weighted MRI. Gene co-expression between cortical areas was more strongly coupled to MIND networks than to MSNs or tractography. MIND network phenotypes were also more heritable, especially edges between structurally differentiated areas. MIND network analysis provides a biologically validated lens for cortical connectomics using readily available MRI data.

Phenotypic effects of genetic variants associated with autism.

While over 100 genes have been associated with autism, little is known about the prevalence of variants affecting them in individuals without a diagnosis of autism. Nor do we fully appreciate the phenotypic diversity beyond the formal autism diagnosis. Based on data from more than 13,000 individuals with autism and 210,000 undiagnosed individuals, we estimated the odds ratios for autism associated to rare loss-of-function (LoF) variants in 185 genes associated with autism, alongside 2,492 genes displaying intolerance to LoF variants. In contrast to autism-centric approaches, we investigated the correlates of these variants in individuals without a diagnosis of autism. We show that these variants are associated with a small but significant decrease in fluid intelligence, qualification level and income and an increase in metrics related to material deprivation. These effects were larger for autism-associated genes than in other LoF-intolerant genes. Using brain imaging data from 21,040 individuals from the UK Biobank, we could not detect significant differences in the overall brain anatomy between LoF carriers and non-carriers. Our results highlight the importance of studying the effect of the genetic variants beyond categorical diagnosis and the need for more research to understand the association between these variants and sociodemographic factors, to best support individuals carrying these variants.

The genetics of autism and steroid-related traits in prenatal and postnatal life.

Background: Autism likelihood is a largely heritable trait. Autism prevalence has a skewed sex ratio, with males being diagnosed more often than females. Steroid hormones play a mediating role in this, as indicated by studies of both prenatal biology and postnatal medical conditions in autistic men and women. It is currently unclear if the genetics of steroid regulation or production interact with the genetic liability for autism. Methods: To address this, two studies were conducted using publicly available datasets, which focused respectively on rare genetic variants linked to autism and neurodevelopmental conditions (study 1) and common genetic variants (study 2) for autism. In Study 1 an enrichment analysis was conducted, between autism-related genes (SFARI database) and genes that are differentially expressed (FDR<0.1) between male and female placentas, in 1st trimester chorionic villi samples of viable pregnancies (n=39). In Study 2 summary statistics of genome wide association studies (GWAS) were used to investigate the genetic correlation between autism and bioactive testosterone, estradiol and postnatal PlGF levels, as well as steroid-related conditions such as polycystic ovaries syndrome (PCOS), age of menarche, and androgenic alopecia. Genetic correlation was calculated based on LD Score regression and results were corrected for multiple testing with FDR. Results: In Study 1, there was significant enrichment of X-linked autism genes in male-biased placental genes, independently of gene length (n=5 genes, p<0.001). In Study 2, common genetic variance associated with autism did not correlate to the genetics for the postnatal levels of testosterone, estradiol or PlGF, but was associated with the genotypes associated with early age of menarche in females (b=-0.109, FDR-q=0.004) and protection from androgenic alopecia for males (b=-0.135, FDR-q=0.007). Conclusion: The rare genetic variants associated with autism appear to interact with placental sex differences, while the common genetic variants associated with autism appear to be involved in the regulation of steroid-related traits. These lines of evidence indicate that the likelihood for autism is partly linked to factors mediating physiological sex differences throughout development.

Identifying rare genetic variants in 21 highly multiplex autism families: the role of diagnosis and autistic traits.

Autism is a highly heritable, heterogeneous, neurodevelopmental condition. Large-scale genetic studies, predominantly focussing on simplex families and clinical diagnoses of autism have identified hundreds of genes associated with autism. Yet, the contribution of these classes of genes to multiplex families and autistic traits still warrants investigation. Here, we conducted whole-genome sequencing of 21 highly multiplex autism families, with at least three autistic individuals in each family, to prioritise genes associated with autism. Using a combination of both autistic traits and clinical diagnosis of autism, we identify rare variants in genes associated with autism, and related neurodevelopmental conditions in multiple families. We identify a modest excess of these variants in autistic individuals compared to individuals without an autism diagnosis. Finally, we identify a convergence of the genes identified in molecular pathways related to development and neurogenesis. In sum, our analysis provides initial evidence to demonstrate the value of integrating autism diagnosis and autistic traits to prioritise genes.