Twin-based Mendelian Randomization Analyses Highlight Smoking's Effects on Blood DNA Methylation, with Putative Reverse Causation.

Epigenome-wide association studies (EWAS) aim to identify differentially methylated loci associated with complex traits and disorders. EWAS of cigarette smoking shows some of the most widespread DNA methylation (DNAm) associations in blood. However, traditional EWAS cannot differentiate between causation and confounding, leading to ambiguity in etiological interpretations. Here, we apply an integrated approach combining Mendelian Randomization and twin-based Direction-of-Causation analyses (MR-DoC) to examine causality underlying smoking-associated blood DNAm changes in the Netherlands Twin Register (N=2577). Evidence across models suggests that current smoking's causal effects on DNAm likely drive many of the previous EWAS findings, implicating functional pathways relevant to several adverse health outcomes of smoking, including hemopoiesis, cell- and neuro-development, and immune regulation. Additionally, we find evidence of potential reverse causal influences at some DNAm sites, with 17 of these sites enriched for gene regulatory functional elements in the brain. The top three sites with evidence of DNAm's effects on smoking annotate to genes involved in G protein-coupled receptor signaling (GNG7, RGS3) and innate immune response (SLC15A4), elucidating potential biological risk factors for smoking. This study highlights the utility of integrating genotypic and DNAm measures in twin cohorts to clarify the causal relationships between health behaviors and blood DNAm.

Brain-age prediction: Systematic evaluation of site effects, and sample age range and size.

Structural neuroimaging data have been used to compute an estimate of the biological age of the brain (brain-age) which has been associated with other biologically and behaviorally meaningful measures of brain development and aging. The ongoing research interest in brain-age has highlighted the need for robust and publicly available brain-age models pre-trained on data from large samples of healthy individuals. To address this need we have previously released a developmental brain-age model. Here we expand this work to develop, empirically validate, and disseminate a pre-trained brain-age model to cover most of the human lifespan. To achieve this, we selected the best-performing model after systematically examining the impact of seven site harmonization strategies, age range, and sample size on brain-age prediction in a discovery sample of brain morphometric measures from 35,683 healthy individuals (age range: 5-90 years; 53.59% female). The pre-trained models were tested for cross-dataset generalizability in an independent sample comprising 2101 healthy individuals (age range: 8-80 years; 55.35% female) and for longitudinal consistency in a further sample comprising 377 healthy individuals (age range: 9-25 years; 49.87% female). This empirical examination yielded the following findings: (1) the accuracy of age prediction from morphometry data was higher when no site harmonization was applied; (2) dividing the discovery sample into two age-bins (5-40 and 40-90 years) provided a better balance between model accuracy and explained age variance than other alternatives; (3) model accuracy for brain-age prediction plateaued at a sample size exceeding 1600 participants. These findings have been incorporated into CentileBrain (https://centilebrain.org/#/brainAGE2), an open-science, web-based platform for individualized neuroimaging metrics.

Exposome-Wide Association Study of Body Mass Index Using a Novel Meta-Analytical Approach for Random Forest Models.

BACKGROUND: Overweight and obesity impose a considerable individual and social burden, and the urban environments might encompass factors that contribute to obesity. Nevertheless, there is a scarcity of research that takes into account the simultaneous interaction of multiple environmental factors. OBJECTIVES: Our objective was to perform an exposome-wide association study of body mass index (BMI) in a multicohort setting of 15 studies. METHODS: Studies were affiliated with the Dutch Geoscience and Health Cohort Consortium (GECCO), had different population sizes (688-141,825), and covered the entire Netherlands. Ten studies contained general population samples, others focused on specific populations including people with diabetes or impaired hearing. BMI was calculated from self-reported or measured height and weight. Associations with 69 residential neighborhood environmental factors (air pollution, noise, temperature, neighborhood socioeconomic and demographic factors, food environment, drivability, and walkability) were explored. Random forest (RF) regression addressed potential nonlinear and nonadditive associations. In the absence of formal methods for multimodel inference for RF, a rank aggregation-based meta-analytic strategy was used to summarize the results across the studies. RESULTS: Six exposures were associated with BMI: five indicating neighborhood economic or social environments (average home values, percentage of high-income residents, average income, livability score, share of single residents) and one indicating the physical activity environment (walkability in 5-km buffer area). Living in high-income neighborhoods and neighborhoods with higher livability scores was associated with lower BMI. Nonlinear associations were observed with neighborhood home values in all studies. Lower neighborhood home values were associated with higher BMI scores but only for values up to €300,000. The directions of associations were less consistent for walkability and share of single residents. DISCUSSION: Rank aggregation made it possible to flexibly combine the results from various studies, although between-study heterogeneity could not be estimated quantitatively based on RF models. Neighborhood social, economic, and physical environments had the strongest associations with BMI. https://doi.org/10.1289/EHP13393.

Breast Cancer Polygenic Risk Score Validation and Effects of Variable Imputation.

Breast cancer (BC) is a complex disease affecting one in eight women in the USA. Advances in population genomics have led to the development of polygenic risk scores (PRSs) with the potential to augment current risk models, but replication is often limited. We evaluated 2 robust PRSs with 313 and 3820 SNPs and the effects of multiple genotype imputation replications in BC cases and control populations. Biological samples from BC cases and cancer-free controls were drawn from three European ancestry cohorts. Genotyping on the Illumina Global Screening Array was followed by stringent quality control measures and 20 genotype imputation replications. A total of 468 unrelated cases and 4337 controls were scored, revealing significant differences in mean PRS percentiles between cases and controls (p < 0.001) for both SNP sets (313-SNP PRS: 52.81 and 48.07; 3820-SNP PRS: 55.45 and 49.81), with receiver operating characteristic curve analysis showing area under the curve values of 0.596 and 0.603 for the 313-SNP and 3820-SNP PRS, respectively. PRS fluctuations (from ~2-3% up to 9%) emerged across imputation iterations. Our study robustly reaffirms the predictive capacity of PRSs for BC by replicating their performance in an independent BC population and showcases the need to average imputed scores for reliable outcomes.

A Multivariate Genome-Wide Association Study Reveals Neural Correlates and Common Biological Mechanisms of Psychopathology Spectra.

There is considerable comorbidity across externalizing and internalizing behavior dimensions of psychopathology. We applied genomic structural equation modeling (gSEM) to genome-wide association study (GWAS) summary statistics to evaluate the factor structure of externalizing and internalizing psychopathology across 16 traits and disorders among European-ancestry individuals (n's = 16,400 to 1,074,629). We conducted GWAS on factors derived from well-fitting models. Downstream analyses served to identify biological mechanisms, explore drug repurposing targets, estimate genetic overlap between the externalizing and internalizing spectra, and evaluate causal effects of psychopathology liability on physical health. Both a correlated factors model, comprising two factors of externalizing and internalizing risk, and a higher-order single-factor model of genetic effects contributing to both spectra demonstrated acceptable fit. GWAS identified 409 lead single nucleotide polymorphisms (SNPs) associated with externalizing and 85 lead SNPs associated with internalizing, while the second-order GWAS identified 256 lead SNPs contributing to broad psychopathology risk. In bivariate causal mixture models, nearly all externalizing and internalizing causal variants overlapped, despite a genetic correlation of only 0.37 (SE = 0.02) between them. Externalizing genes showed cell-type specific expression in GABAergic, cortical, and hippocampal neurons, and internalizing genes were associated with reduced subcallosal cortical volume, providing insight into the neurobiological underpinnings of psychopathology. Genetic liability for externalizing, internalizing, and broad psychopathology exerted causal effects on pain, general health, cardiovascular diseases, and chronic illnesses. These findings underscore the complex genetic architecture of psychopathology, identify potential biological pathways for the externalizing and internalizing spectra, and highlight the physical health burden of psychiatric comorbidity.

Genome-wide analyses identify 21 infertility loci and over 400 reproductive hormone loci across the allele frequency spectrum.

Genome-wide association studies (GWASs) may help inform treatments for infertility, whose causes remain unknown in many cases. Here we present GWAS meta-analyses across six cohorts for male and female infertility in up to 41,200 cases and 687,005 controls. We identified 21 genetic risk loci for infertility (P≤5E-08), of which 12 have not been reported for any reproductive condition. We found positive genetic correlations between endometriosis and all-cause female infertility (rg=0.585, P=8.98E-14), and between polycystic ovary syndrome and anovulatory infertility (rg=0.403, P=2.16E-03). The evolutionary persistence of female infertility-risk alleles in EBAG9 may be explained by recent directional selection. We additionally identified up to 269 genetic loci associated with follicle-stimulating hormone (FSH), luteinising hormone, oestradiol, and testosterone through sex-specific GWAS meta-analyses (N=6,095-246,862). While hormone-associated variants near FSHB and ARL14EP colocalised with signals for anovulatory infertility, we found no rg between female infertility and reproductive hormones (P>0.05). Exome sequencing analyses in the UK Biobank (N=197,340) revealed that women carrying testosterone-lowering rare variants in GPC2 were at higher risk of infertility (OR=2.63, P=1.25E-03). Taken together, our results suggest that while individual genes associated with hormone regulation may be relevant for fertility, there is limited genetic evidence for correlation between reproductive hormones and infertility at the population level. We provide the first comprehensive view of the genetic architecture of infertility across multiple diagnostic criteria in men and women, and characterise its relationship to other health conditions.

Enhanced resolution profiling in twins reveals differential methylation signatures of type 2 diabetes with links to its complications.

BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression. METHODS: We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites. FINDINGS: We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2). INTERPRETATION: The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications. FUNDING: Funding acknowledgements for each cohort can be found in the Supplementary Note.

Genome-Wide Association Study of Obsessive-Compulsive Symptoms including 33,943 individuals from the general population.

While 1-2% of individuals meet the criteria for a clinical diagnosis of obsessive-compulsive disorder (OCD), many more (~13-38%) experience subclinical obsessive-compulsive symptoms (OCS) during their life. To characterize the genetic underpinnings of OCS and its genetic relationship to OCD, we conducted the largest genome-wide association study (GWAS) meta-analysis of parent- or self-reported OCS to date (N = 33,943 with complete phenotypic and genome-wide data), combining the results from seven large-scale population-based cohorts from Sweden, the Netherlands, England, and Canada (including six twin cohorts and one cohort of unrelated individuals). We found no genome-wide significant associations at the single-nucleotide polymorphism (SNP) or gene-level, but a polygenic risk score (PRS) based on the OCD GWAS previously published by the Psychiatric Genetics Consortium (PGC-OCD) was significantly associated with OCS (Pfixed = 3.06 × 10-5). Also, one curated gene set (Mootha Gluconeogenesis) reached Bonferroni-corrected significance (Ngenes = 28, Beta = 0.79, SE = 0.16, Pbon = 0.008). Expression of genes in this set is high at sites of insulin mediated glucose disposal. Dysregulated insulin signaling in the etiology of OCS has been suggested by a previous study describing a genetic overlap of OCS with insulin signaling-related traits in children and adolescents. We report a SNP heritability of 4.1% (P = 0.0044) in the meta-analyzed GWAS, and heritability estimates based on the twin cohorts of 33-43%. Genetic correlation analysis showed that OCS were most strongly associated with OCD (rG = 0.72, p = 0.0007) among all tested psychiatric disorders (N = 11). Of all 97 tested phenotypes, 24 showed a significant genetic correlation with OCS, and 66 traits showed concordant directions of effect with OCS and OCD. OCS have a significant polygenic contribution and share genetic risk with diagnosed OCD, supporting the hypothesis that OCD represents the extreme end of widely distributed OCS in the population.

Genomics Research of Lifetime Depression in the Netherlands: The BIObanks Netherlands Internet Collaboration (BIONIC) Project.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

Are Genetic and Environmental Risk Factors for Psychopathology Amplified in Children with Below-Average Intelligence? A Population-Based Twin Study.

There is a negative association between intelligence and psychopathology. We analyzed data on intelligence and psychopathology to assess this association in seven-year-old Dutch twin pairs (ranging from 616 to 14,150 depending on the phenotype) and estimated the degree to which genetic and environmental factors common to intelligence and psychopathology explain the association. Secondly, we examined whether genetic and environmental effects on psychopathology are moderated by intelligence. We found that intelligence, as assessed by psychometric IQ tests, correlated negatively with childhood psychopathology, as assessed by the DSM-oriented scales of the Child Behavior Check List (CBCL). The correlations ranged between - .09 and - .15 and were mainly explained by common genetic factors. Intelligence moderated genetic and environmental effects on anxiety and negative affect, but not those on ADHD, ODD, and autism. The heritability of anxiety and negative affect was greatest in individuals with below-average intelligence. We discuss mechanisms through which this effect could arise, and we end with some recommendations for future research.

Normative modelling of brain morphometry across the lifespan with CentileBrain: algorithm benchmarking and model optimisation.

The value of normative models in research and clinical practice relies on their robustness and a systematic comparison of different modelling algorithms and parameters; however, this has not been done to date. We aimed to identify the optimal approach for normative modelling of brain morphometric data through systematic empirical benchmarking, by quantifying the accuracy of different algorithms and identifying parameters that optimised model performance. We developed this framework with regional morphometric data from 37 407 healthy individuals (53% female and 47% male; aged 3-90 years) from 87 datasets from Europe, Australia, the USA, South Africa, and east Asia following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The multivariate fractional polynomial regression (MFPR) emerged as the preferred algorithm, optimised with non-linear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3000 study participants. This model can inform about the biological and behavioural implications of deviations from typical age-related neuroanatomical changes and support future study designs. The model and scripts described here are freely available through CentileBrain.

Using Instrumental Variables to Measure Causation over Time in Cross-Lagged Panel Models.

Cross-lagged panel models (CLPMs) are commonly used to estimate causal influences between two variables with repeated assessments. The lagged effects in a CLPM depend on the time interval between assessments, eventually becoming undetectable at longer intervals. To address this limitation, we incorporate instrumental variables (IVs) into the CLPM with two study waves and two variables. Doing so enables estimation of both the lagged (i.e., "distal") effects and the bidirectional cross-sectional (i.e., "proximal") effects at each wave. The distal effects reflect Granger-causal influences across time, which decay with increasing time intervals. The proximal effects capture causal influences that accrue over time and can help infer causality when the distal effects become undetectable at longer intervals. Significant proximal effects, with a negligible distal effect, would imply that the time interval is too long to estimate a lagged effect at that time interval using the standard CLPM. Through simulations and an empirical application, we demonstrate the impact of time intervals on causal inference in the CLPM and present modeling strategies to detect causal influences regardless of the time interval in a study. Furthermore, to motivate empirical applications of the proposed model, we highlight the utility and limitations of using genetic variables as IVs in large-scale panel studies.

Family History of Twinning and Fertility Traits in Nigerian Mothers of Dizygotic Twins.

Familial twinning and fertility traits were investigated in Nigerian mothers of dizygotic (DZ) twins (MoDZT; N = 972) and controls (N = 525) who responded to our person-to-person interview, which included questions on pregnancy history and family history of DZ twinning. Controls were defined as women who are not twins themselves and do not have twins in their first-degree relatives. Over 95% of the participants were Yoruba. We found that Nigerian MoDZT had an average of 4.0 (±2.6) pairs of twins among their relatives, and of these, the prevalence of DZ twins was significantly higher than that of monozygotic (MZ) twins (45.9% vs. 25.8%). Controls had an average of 0.5 (±0.4) pairs, and over 95% of the controls had no twins in their relatives. These results suggest genetic influences on DZ twinning in Nigerians. MoDZT were significantly younger in their mean age at first child, and had higher parity than controls, suggesting increased fertility in MoDZT. As compared to mothers with a single set of twins, mothers (N = 130) with multiple sets had significantly more twins among their relatives (5.4 pairs vs. 3.7 pairs) and had their first twins at a younger age (28.4 vs. 30.7 years), indicating that mothers with multiple sets of twins might have higher genetic propensity for twinning associated with earlier age at twin pregnancy. Our findings argue for genomewide association studies for DZ twinning in Nigerians, and may help to develop intervention strategies to overcome infertility/subfertility problems.

A multi-omics data analysis workflow packaged as a FAIR Digital Object.

BACKGROUND: Applying good data management and FAIR (Findable, Accessible, Interoperable, and Reusable) data principles in research projects can help disentangle knowledge discovery, study result reproducibility, and data reuse in future studies. Based on the concepts of the original FAIR principles for research data, FAIR principles for research software were recently proposed. FAIR Digital Objects enable discovery and reuse of Research Objects, including computational workflows for both humans and machines. Practical examples can help promote the adoption of FAIR practices for computational workflows in the research community. We developed a multi-omics data analysis workflow implementing FAIR practices to share it as a FAIR Digital Object. FINDINGS: We conducted a case study investigating shared patterns between multi-omics data and childhood externalizing behavior. The analysis workflow was implemented as a modular pipeline in the workflow manager Nextflow, including containers with software dependencies. We adhered to software development practices like version control, documentation, and licensing. Finally, the workflow was described with rich semantic metadata, packaged as a Research Object Crate, and shared via WorkflowHub. CONCLUSIONS: Along with the packaged multi-omics data analysis workflow, we share our experiences adopting various FAIR practices and creating a FAIR Digital Object. We hope our experiences can help other researchers who develop omics data analysis workflows to turn FAIR principles into practice.

Heritability of adult picky eating in the Netherlands.

Adult picky eating (APE), the rejection of familiar and unfamiliar foods leading to a diet with limited variety, is an understudied phenomenon which can have both physical and psychological negative consequences. The aetiology of individual differences in APE is understudied, although there is reason to believe that it is partly heritable. Therefore, we aimed to estimate the heritability of APE with data from the Netherlands Twin Register (n = 8016) with classical genetic structural equation modelling. In order to use these data, we firstly investigated whether a Food Preference Questionnaire (FPQ) could measure APE with a pre-registered prestudy. Adult participants (n = 414) filled in online questionnaires, including a FPQ and measures related to APE. Spearman's rho correlation quantified the relationship between different elements of the Dutch FPQ and different scores on measures of APE. Results of the prestudy showed that the mean liking score on the FPQ could be used to measure APE (ρ > .50). This measure was then used in the main study to estimate the heritability of APE. Results showed that broad-sense heritability for APE is 49 % (additive genetic effects 14 % (95 % CI [00, 38]) + dominance genetic effects 35 % (95 % CI [11, 52]), while the remaining variance is explained by unique environmental factors. Future studies may focus on uncovering the specific genetic and unique environmental factors that play a role in APE.

Evaluating the causal relationship between educational attainment and mental health.

We investigate the causal relationship between educational attainment (EA) and mental health using two research designs. First, we compare the relationship between EA and 18 psychiatric diagnoses within sibship in Dutch national registry data (N=1.7 million), thereby controlling for unmeasured familial factors. Second, we apply two-sample Mendelian Randomization, which uses genetic variants related to EA or psychiatric diagnosis as instrumental variables, to test whether there is a causal relation in either direction. Our results suggest that lower levels of EA causally increase the risk of MDD, ADHD, alcohol dependence, GAD and PTSD diagnoses. We also find evidence of a causal effect of ADHD on EA. For schizophrenia, anorexia nervosa, OCD, and bipolar disorder, results were inconsistent across the different approaches, highlighting the importance of using multiple research designs to understand complex relationships such as between EA and mental health.

Integrative multi-omics analysis of genomic, epigenomic, and metabolomics data leads to new insights for Attention-Deficit/Hyperactivity Disorder.

The evolving field of multi-omics combines data and provides methods for simultaneous analysis across several omics levels. Here, we integrated genomics (transmitted and non-transmitted polygenic scores [PGSs]), epigenomics, and metabolomics data in a multi-omics framework to identify biomarkers for Attention-Deficit/Hyperactivity Disorder (ADHD) and investigated the connections among the three omics levels. We first trained single- and next multi-omics models to differentiate between cases and controls in 596 twins (cases = 14.8%) from the Netherlands Twin Register (NTR) demonstrating reasonable in-sample prediction through cross-validation. The multi-omics model selected 30 PGSs, 143 CpGs, and 90 metabolites. We confirmed previous associations of ADHD with glucocorticoid exposure and the transmembrane protein family TMEM, show that the DNA methylation of the MAD1L1 gene associated with ADHD has a relation with parental smoking behavior, and present novel findings including associations between indirect genetic effects and CpGs of the STAP2 gene. However, out-of-sample prediction in NTR participants (N = 258, cases = 14.3%) and in a clinical sample (N = 145, cases = 51%) did not perform well (range misclassification was [0.40, 0.57]). The results highlighted connections between omics levels, with the strongest connections between non-transmitted PGSs, CpGs, and amino acid levels and show that multi-omics designs considering interrelated omics levels can help unravel the complex biology underlying ADHD.

New insights into the (epi)genetics of twinning.

Spontaneous dizygotic (DZ) twins, i.e. twins conceived without the use of ARTs, run in families and their prevalence varies widely around the globe. In contrast, monozygotic (MZ) twins occur at a constant rate across time and geographical regions and, with some rare exceptions, do not cluster in families. The leading hypothesis for MZ twins, which arise when a zygote splits during preimplantation stages of development, is random occurrence. We have found the first series of genes underlying the liability of being the mother of DZ twins and have shown that being an MZ twin is strongly associated with a stable DNA methylation signature in child and adult somatic tissues. Because identical twins keep this molecular signature across the lifespan, this discovery opens up completely new possibilities for the retrospective diagnosis of whether a person is an MZ twin whose co-twin may have vanished in the early stages of pregnancy. Here, we summarize the gene finding results for mothers of DZ twins based on genetic association studies followed by meta-analysis, and further present the striking epigenetic results for MZ twins.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers.

BACKGROUND: Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure. METHODS: Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference. RESULTS: For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness. CONCLUSIONS: We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Normative Modeling of Brain Morphometry Across the Lifespan Using CentileBrain: Algorithm Benchmarking and Model Optimization.

We present an empirically benchmarked framework for sex-specific normative modeling of brain morphometry that can inform about the biological and behavioral significance of deviations from typical age-related neuroanatomical changes and support future study designs. This framework was developed using regional morphometric data from 37,407 healthy individuals (53% female; aged 3-90 years) following a comparative evaluation of eight algorithms and multiple covariate combinations pertaining to image acquisition and quality, parcellation software versions, global neuroimaging measures, and longitudinal stability. The Multivariate Factorial Polynomial Regression (MFPR) emerged as the preferred algorithm optimized using nonlinear polynomials for age and linear effects of global measures as covariates. The MFPR models showed excellent accuracy across the lifespan and within distinct age-bins, and longitudinal stability over a 2-year period. The performance of all MFPR models plateaued at sample sizes exceeding 3,000 study participants. The model and scripts described here are freely available through CentileBrain (https://centilebrain.org/).