Quality control and analytic best practices for testing genetic models of sex differences in large populations.

Phenotypic sex-based differences exist for many complex traits. In other cases, phenotypes may be similar, but underlying biology may vary. Thus, sex-aware genetic analyses are becoming increasingly important for understanding the mechanisms driving these differences. To this end, we provide a guide outlining the current best practices for testing various models of sex-dependent genetic effects in complex traits and disease conditions, noting that this is an evolving field. Insights from sex-aware analyses will not only teach us about the biology of complex traits but also aid in achieving the goals of precision medicine and health equity for all.

Sex modifies the effect of genetic risk scores for polycystic ovary syndrome on metabolic phenotypes.

Females with polycystic ovary syndrome (PCOS), the most common endocrine disorder in women, have an increased risk of developing cardiometabolic disorders such as insulin resistance, obesity, and type 2 diabetes (T2D). While only diagnosable in females, males with a family history of PCOS can also exhibit a poor cardiometabolic profile. Therefore, we aimed to elucidate the role of sex in the cardiometabolic comorbidities observed in PCOS by conducting bidirectional genetic risk score analyses in both sexes. We first conducted a phenome-wide association study (PheWAS) using PCOS polygenic risk scores (PCOSPRS) to identify potential pleiotropic effects of PCOSPRS across 1,380 medical conditions recorded in the Vanderbilt University Medical Center electronic health record (EHR) database, in females and males. After adjusting for age and genetic ancestry, we found that European (EUR)-ancestry males with higher PCOSPRS were significantly more likely to develop hypertensive diseases than females at the same level of genetic risk. We performed the same analysis in an African (AFR)-ancestry population, but observed no significant associations, likely due to poor trans-ancestry performance of the PRS. Based on observed significant associations in the EUR-ancestry population, we then tested whether the PRS for comorbid conditions (e.g., T2D, body mass index (BMI), hypertension, etc.) also increased the odds of a PCOS diagnosis. Only BMIPRS and T2DPRS were significantly associated with a PCOS diagnosis in EUR-ancestry females. We then further adjusted the T2DPRS for measured BMI and BMIresidual (regressed on the BMIPRS and enriched for the environmental contribution to BMI). Results demonstrated that genetically regulated BMI primarily accounted for the relationship between T2DPRS and PCOS. Overall, our findings show that the genetic architecture of PCOS has distinct sex differences in associations with genetically correlated cardiometabolic traits. It is possible that the cardiometabolic comorbidities observed in PCOS are primarily explained by their shared genetic risk factors, which can be further influenced by biological variables including sex and BMI.

The role of sex in the genomics of human complex traits.

Nearly all human complex traits and disease phenotypes exhibit some degree of sex differences, including differences in prevalence, age of onset, severity or disease progression. Until recently, the underlying genetic mechanisms of such sex differences have been largely unexplored. Advances in genomic technologies and analytical approaches are now enabling a deeper investigation into the effect of sex on human health traits. In this Review, we discuss recent insights into the genetic models and mechanisms that lead to sex differences in complex traits. This knowledge is critical for developing deeper insight into the fundamental biology of sex differences and disease processes, thus facilitating precision medicine.

Author Correction: The role of sex in the genomics of human complex traits.

In Box 4 of the originally published article, the text describing the Miami plot in part c of the figure contained some minor errors. Specifically, during pre-publication revision of the article, the authors updated the illustrative Miami plot (shown in figure part c) from that of reference 80 (Randall et al. PLoS Genet. (2013)) to a more recent study in reference 82 (Winkler et al. PLoS Genet. (2015)). The box text has now been updated to reflect that change. In paragraph 2, the trait has been updated from "hip circumference adjusted for body mass index" to "waist-to-hip ratio adjusted for body mass index (under 50 years old)" and it has been clarified that female GWAS data are shown on the top half of the plot with male data at the bottom. The original two citations of reference 80 in the Box 4 text have been updated to reference 82. Finally, a typographical artefact was corrected on the Y axis of the Miami plot, whereby the labels '14' and '16' in the top half of the plot we originally both shown as '12'. None of these corrections alter the overall illustrative point that genetic architectures for traits can differ between males and females, which was a conclusion of both reference 80 and reference 82.

Body mass index stratified meta-analysis of genome-wide association studies of polycystic ovary syndrome in women of European ancestry.

BACKGROUND: Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder with a substantial genetic component. However, the clinical manifestations of PCOS are heterogeneous with notable differences between lean and obese women, implying a different pathophysiology manifesting in differential body mass index (BMI). We performed a meta-analysis of genome-wide association study (GWAS) data from six well-characterised cohorts, using a case-control study design stratified by BMI, aiming to identify genetic variants associated with lean and overweight/obese PCOS subtypes. RESULTS: The study comprised 254,588 women (5,937 cases and 248,651 controls) from individual studies performed in Australia, Estonia, Finland, the Netherlands and United States of America, and separated according to three BMI stratifications (lean, overweight and obese). Genome-wide association analyses were performed for each stratification within each cohort, with the data for each BMI group meta-analysed using METAL software. Almost half of the total study population (47%, n = 119,584) were of lean BMI (≤ 25 kg/m2). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10-12) and rs2228260 within XBP1 (P = 3.68 × 10-8). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10-6). There were no significant loci observed for the overweight or obese sub-strata when analysed separately, however, when these strata were combined, an association signal led by rs569675099 within DENND1A reached genome-wide significance (P = 3.22 × 10-9) and a gene-based association was identified with ERBB4 (P = 1.59 × 10-6). Nineteen of 28 signals identified in previous GWAS, were replicated with consistent allelic effect in the lean stratum. There were less replicated signals in the overweight and obese groups, and only 4 SNPs were replicated in each of the three BMI strata. CONCLUSIONS: Genetic variation at the XBP1, LINC02905 and ERBB4 loci were associated with PCOS within unique BMI strata, while DENND1A demonstrated associations across multiple strata, providing evidence of both distinct and shared genetic features between lean and overweight/obese PCOS-affected women. This study demonstrated that PCOS-affected women with contrasting body weight are not only phenotypically distinct but also show variation in genetic architecture; lean PCOS women typically display elevated gonadotrophin ratios, lower insulin resistance, higher androgen levels, including adrenal androgens, and more favourable lipid profiles. Overall, these findings add to the growing body of evidence supporting a genetic basis for PCOS as well as differences in genetic patterns relevant to PCOS BMI-subtype.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Phenomewide Association Study of Health Outcomes Associated With the Genetic Correlates of 25 Hydroxyvitamin D Concentration and Vitamin D Binding Protein Concentration.

While it is known that vitamin D deficiency is associated with adverse bone outcomes, it remains unclear whether low vitamin D status may increase the risk of a wider range of health outcomes. We had the opportunity to explore the association between common genetic variants associated with both 25 hydroxyvitamin D (25OHD) and the vitamin D binding protein (DBP, encoded by the GC gene) with a comprehensive range of health disorders and laboratory tests in a large academic medical center. We used summary statistics for 25OHD and DBP to generate polygenic scores (PGS) for 66,482 participants with primarily European ancestry and 13,285 participants with primarily African ancestry from the Vanderbilt University Medical Center Biobank (BioVU). We examined the predictive properties of PGS25OHD, and two scores related to DBP concentration with respect to 1322 health-related phenotypes and 315 laboratory-measured phenotypes from electronic health records. In those with European ancestry: (a) the PGS25OHD and PGSDBP scores, and individual SNPs rs4588 and rs7041 were associated with both 25OHD concentration and 1,25 dihydroxyvitamin D concentrations; (b) higher PGS25OHD was associated with decreased concentrations of triglycerides and cholesterol, and reduced risks of vitamin D deficiency, disorders of lipid metabolism, and diabetes. In general, the findings for the African ancestry group were consistent with findings from the European ancestry analyses. Our study confirms the utility of PGS and two key variants within the GC gene (rs4588 and rs7041) to predict the risk of vitamin D deficiency in clinical settings and highlights the shared biology between vitamin D-related genetic pathways a range of health outcomes.

Nonsynonymous single-nucleotide polymorphisms in the G6PC2 gene affect protein expression, enzyme activity, and fasting blood glucose.

G6PC2 encodes a glucose-6-phosphatase (G6Pase) catalytic subunit that modulates the sensitivity of insulin secretion to glucose and thereby regulates fasting blood glucose (FBG). A common single-nucleotide polymorphism (SNP) in G6PC2, rs560887 is an important determinant of human FBG variability. This SNP has a subtle effect on G6PC2 RNA splicing, which raises the question as to whether nonsynonymous SNPs with a major impact on G6PC2 stability or enzyme activity might have a broader disease/metabolic impact. Previous attempts to characterize such SNPs were limited by the very low inherent G6Pase activity and expression of G6PC2 protein in islet-derived cell lines. In this study, we describe the use of a plasmid vector that confers high G6PC2 protein expression in islet cells, allowing for a functional analysis of 22 nonsynonymous G6PC2 SNPs, 19 of which alter amino acids that are conserved in mouse G6PC2 and the human and mouse variants of the related G6PC1 isoform. We show that 16 of these SNPs markedly impair G6PC2 protein expression (>50% decrease). These SNPs have variable effects on the stability of human and mouse G6PC1, despite the high sequence homology between these isoforms. Four of the remaining six SNPs impaired G6PC2 enzyme activity. Electronic health record-derived phenotype analyses showed an association between high-impact SNPs and FBG, but not other diseases/metabolites. While homozygous G6pc2 deletion in mice increases the risk of hypoglycemia, these human data reveal no evidence that the beneficial use of partial G6PC2 inhibitors to lower FBG would be associated with unintended negative consequences.

Fate or coincidence: do COPD and major depression share genetic risk factors?

Major depressive disorder (MDD) is a common comorbidity in chronic obstructive pulmonary disease (COPD), affecting up to 57% of patients with COPD. Although the comorbidity of COPD and MDD is well established, the causal relationship between these two diseases is unclear. A large-scale electronic health record clinical biobank and genome-wide association study summary statistics for MDD and lung function traits were used to investigate potential shared underlying genetic susceptibility between COPD and MDD. Linkage disequilibrium score regression was used to estimate genetic correlation between phenotypes. Polygenic risk scores (PRS) for MDD and lung function traits were developed and used to perform a phenome-wide association study (PheWAS). Multi-trait-based conditional and joint analysis identified single-nucleotide polymorphisms (SNPs) influencing both lung function and MDD. We found genetic correlations between MDD and all lung function traits were small and not statistically significant. A PRS-MDD was significantly associated with an increased risk of COPD in a PheWAS [odds ratio (OR) = 1.12, 95% confidence interval (CI): 1.09-1.16] when adjusting for age, sex and genetic ancestry, but this relationship became attenuated when controlling for smoking history (OR = 1.08, 95% CI: 1.04-1.13). No significant associations were found between the lung function PRS and MDD. Multi-trait-based conditional and joint analysis identified three SNPs that may contribute to both traits, two of which were previously associated with mood disorders and COPD. Our findings suggest that the observed relationship between COPD and MDD may not be driven by a strong shared genetic architecture.

The genetic contribution to the comorbidity of depression and anxiety: a multi-site electronic health records study of almost 178 000 people.

BACKGROUND: Depression and anxiety are common and highly comorbid, and their comorbidity is associated with poorer outcomes posing clinical and public health concerns. We evaluated the polygenic contribution to comorbid depression and anxiety, and to each in isolation. METHODS: Diagnostic codes were extracted from electronic health records for four biobanks [N = 177 865 including 138 632 European (77.9%), 25 612 African (14.4%), and 13 621 Hispanic (7.7%) ancestry participants]. The outcome was a four-level variable representing the depression/anxiety diagnosis group: neither, depression-only, anxiety-only, and comorbid. Multinomial regression was used to test for association of depression and anxiety polygenic risk scores (PRSs) with the outcome while adjusting for principal components of ancestry. RESULTS: In total, 132 960 patients had neither diagnosis (74.8%), 16 092 depression-only (9.0%), 13 098 anxiety-only (7.4%), and 16 584 comorbid (9.3%). In the European meta-analysis across biobanks, both PRSs were higher in each diagnosis group compared to controls. Notably, depression-PRS (OR 1.20 per s.d. increase in PRS; 95% CI 1.18-1.23) and anxiety-PRS (OR 1.07; 95% CI 1.05-1.09) had the largest effect when the comorbid group was compared with controls. Furthermore, the depression-PRS was significantly higher in the comorbid group than the depression-only group (OR 1.09; 95% CI 1.06-1.12) and the anxiety-only group (OR 1.15; 95% CI 1.11-1.19) and was significantly higher in the depression-only group than the anxiety-only group (OR 1.06; 95% CI 1.02-1.09), showing a genetic risk gradient across the conditions and the comorbidity. CONCLUSIONS: This study suggests that depression and anxiety have partially independent genetic liabilities and the genetic vulnerabilities to depression and anxiety make distinct contributions to comorbid depression and anxiety.

Genetically regulated multi-omics study for symptom clusters of posttraumatic stress disorder highlights pleiotropy with hematologic and cardio-metabolic traits.

Posttraumatic stress disorder (PTSD) is a psychiatric disorder that may arise in response to severe traumatic event and is diagnosed based on three main symptom clusters (reexperiencing, avoidance, and hyperarousal) per the Diagnostic Manual of Mental Disorders (version DSM-IV-TR). In this study, we characterized the biological heterogeneity of PTSD symptom clusters by performing a multi-omics investigation integrating genetically regulated gene, splicing, and protein expression in dorsolateral prefrontal cortex tissue within a sample of US veterans enrolled in the Million Veteran Program (N total = 186,689). We identified 30 genes in 19 regions across the three PTSD symptom clusters. We found nine genes to have cell-type specific expression, and over-representation of miRNA-families - miR-148, 30, and 8. Gene-drug target prioritization approach highlighted cyclooxygenase and acetylcholine compounds. Next, we tested molecular-profile based phenome-wide impact of identified genes with respect to 1678 phenotypes derived from the Electronic Health Records of the Vanderbilt University biorepository (N = 70,439). Lastly, we tested for local genetic correlation across PTSD symptom clusters which highlighted metabolic (e.g., obesity, diabetes, vascular health) and laboratory traits (e.g., neutrophil, eosinophil, tau protein, creatinine kinase). Overall, this study finds comprehensive genomic evidence including clinical and regulatory profiles between PTSD, hematologic and cardiometabolic traits, that support comorbidities observed in epidemiologic studies of PTSD.

Life is pain: Fibromyalgia as a nexus of multiple liability distributions.

Fibromyalgia is a complex disease of unclear etiology that is complicated by difficulties in diagnosis, treatment, and clinical heterogeneity. To clarify this etiology, healthcare-based data are leveraged to assess the influences on fibromyalgia in several domains. Prevalence is less than 1% of females in our population register data, and about 1/10th that in males. Fibromyalgia often presents with co-occurring conditions including back pain, rheumatoid arthritis, and anxiety. More comorbidities are identified with hospital-associated biobank data, falling into three broad categories of pain-related, autoimmune, and psychiatric disorders. Selecting representative phenotypes with published genome-wide association results for polygenic scoring, we confirm genetic predispositions to psychiatric, pain sensitivity, and autoimmune conditions show associations with fibromyalgia, although these may differ by ancestry group. We conduct a genome-wide association analysis of fibromyalgia in biobank samples, which did not result in any genome-wide significant loci; further studies with increased sample size are necessary to identify specific genetic effects on fibromyalgia. Overall, fibromyalgia appears to have strong clinical and likely genetic links to several disease categories, and could usefully be understood as a composite manifestation of these etiological sources.

Frequency of benign neutropenia among Black versus White individuals undergoing a bone marrow assessment.

Healthy individuals in the United States identified as having Black race have lower neutrophil counts, on average, than individuals identified as having White race, which could result in more negative diagnostic evaluations for neutropenia. To test this hypothesis, the proportion of evaluations where the final diagnosis was clinically insignificant neutropenia for Black and White individuals who underwent an evaluation by a haematologist that included a bone marrow (BM) biopsy to investigate neutropenia was assessed. 172 individuals without prior haematological diagnoses who underwent a haematological evaluation to investigate neutropenia. Individuals diagnosed with clinically insignificant neutropenia between Black and White individuals were compared using a propensity-score-adjusted logistic regression. Of 172 individuals, 42 (24%) were classified as Black race, 86 (50%) were males, and the 79 (46%) were over 18 years old. A BM biopsy did not identify pathology in 95% (40 of 42) of Black individuals and 68% (89 of 130) of White Individuals. Black individuals (25 of 42 [60%]) received a final diagnosis of clinically insignificant neutropenia, compared to White individuals (12 of 130 [9%]) (adjusted odds ratio =7.9, 95% CI: 3.1 - 21.1). We conclude that black individuals were more likely to receive a diagnosis of clinically insignificant neutropenia after haematological assessment.

Dense phenotyping from electronic health records enables machine learning-based prediction of preterm birth.

BACKGROUND: Identifying pregnancies at risk for preterm birth, one of the leading causes of worldwide infant mortality, has the potential to improve prenatal care. However, we lack broadly applicable methods to accurately predict preterm birth risk. The dense longitudinal information present in electronic health records (EHRs) is enabling scalable and cost-efficient risk modeling of many diseases, but EHR resources have been largely untapped in the study of pregnancy. METHODS: Here, we apply machine learning to diverse data from EHRs with 35,282 deliveries to predict singleton preterm birth. RESULTS: We find that machine learning models based on billing codes alone can predict preterm birth risk at various gestational ages (e.g., ROC-AUC = 0.75, PR-AUC = 0.40 at 28 weeks of gestation) and outperform comparable models trained using known risk factors (e.g., ROC-AUC = 0.65, PR-AUC = 0.25 at 28 weeks). Examining the patterns learned by the model reveals it stratifies deliveries into interpretable groups, including high-risk preterm birth subtypes enriched for distinct comorbidities. Our machine learning approach also predicts preterm birth subtypes (spontaneous vs. indicated), mode of delivery, and recurrent preterm birth. Finally, we demonstrate the portability of our approach by showing that the prediction models maintain their accuracy on a large, independent cohort (5978 deliveries) from a different healthcare system. CONCLUSIONS: By leveraging rich phenotypic and genetic features derived from EHRs, we suggest that machine learning algorithms have great potential to improve medical care during pregnancy. However, further work is needed before these models can be applied in clinical settings.

LGBTQ+ Perspectives on Conducting Genomic Research on Sexual Orientation and Gender Identity.

We conducted in-depth, semi-structured interviews with LGBTQ+-identified individuals (n = 31) to explore the range of LGBTQ+ perspectives on genomic research using either sexual orientation or gender identity (SOGI) data. Most interviewees presumed that research would confirm genetic contributions to sexual orientation and gender identity. Primary hopes for such confirmation included validating LGBTQ+ identities, improved access to and quality of healthcare and other resources, and increased acceptance in familial, socio-cultural, and political environments. Areas of concern included threats of pathologizing and medicalizing LGBTQ+ identities and experiences, undermining reproductive rights, gatekeeping of health or social systems, and malicious testing or misuse of genetic results, particularly for LGBTQ+ youth. Overall, interviewees were divided on the acceptability of genomic research investigating genetic contributions to sexual orientation and gender identity. Participants emphasized researchers' ethical obligations to LGBTQ+ individuals and endorsed engagement with LGBTQ+ communities throughout all aspects of genomic research using SOGI data.

Genome-wide association study of problematic opioid prescription use in 132,113 23andMe research participants of European ancestry.

The growing prevalence of opioid use disorder (OUD) constitutes an urgent health crisis. Ample evidence indicates that risk for OUD is heritable. As a surrogate (or proxy) for OUD, we explored the genetic basis of using prescription opioids 'not as prescribed'. We hypothesized that misuse of opiates might be a heritable risk factor for OUD. To test this hypothesis, we performed a genome-wide association study (GWAS) of problematic opioid use (POU) in 23andMe research participants of European ancestry (N = 132,113; 21% cases). We identified two genome-wide significant loci (rs3791033, an intronic variant of KDM4A; rs640561, an intergenic variant near LRRIQ3). POU showed positive genetic correlations with the two largest available GWAS of OUD and opioid dependence (rg = 0.64, 0.80, respectively). We also identified numerous additional genetic correlations with POU, including alcohol dependence (rg = 0.74), smoking initiation (rg = 0.63), pain relief medication intake (rg = 0.49), major depressive disorder (rg = 0.44), chronic pain (rg = 0.42), insomnia (rg = 0.39), and loneliness (rg = 0.28). Although POU was positively genetically correlated with risk-taking (rg = 0.38), conditioning POU on risk-taking did not substantially alter the magnitude or direction of these genetic correlations, suggesting that POU does not simply reflect a genetic tendency towards risky behavior. Lastly, we performed phenome- and lab-wide association analyses, which uncovered additional phenotypes that were associated with POU, including respiratory failure, insomnia, ischemic heart disease, and metabolic and blood-related biomarkers. We conclude that opioid misuse can be measured in population-based cohorts and provides a cost-effective complementary strategy for understanding the genetic basis of OUD.

Genetic risk for major depressive disorder and loneliness in sex-specific associations with coronary artery disease.

Major depressive disorder (MDD) and loneliness are phenotypically and genetically correlated with coronary artery disease (CAD), but whether these associations are explained by pleiotropic genetic variants or shared comorbidities is unclear. To tease apart these scenarios, we first assessed the medical morbidity pattern associated with genetic risk factors for MDD and loneliness by conducting a phenome-wide association study in 18,385 European-ancestry individuals in the Vanderbilt University Medical Center biobank, BioVU. Polygenic scores for MDD and loneliness were developed for each person using previously published meta-GWAS summary statistics, and were tested for association with 882 clinical diagnoses ascertained via billing codes in electronic health records. We discovered strong associations with heart disease diagnoses, and next embarked on targeted analyses of CAD in 3893 cases and 4197 controls. We found odds ratios of 1.11 (95% CI, 1.04-1.18; P 8.43 × 10-4) and 1.13 (95% CI, 1.07-1.20; P 4.51 × 10-6) per 1-SD increase in the polygenic scores for MDD and loneliness, respectively. Results were similar in patients without psychiatric symptoms, and the increased risk persisted in females even after adjusting for multiple conventional risk factors and a polygenic score for CAD. In a final sensitivity analysis, we statistically adjusted for the genetic correlation between MDD and loneliness and re-computed polygenic scores. The polygenic score unique to loneliness remained associated with CAD (OR 1.09, 95% CI 1.03-1.15; P 0.002), while the polygenic score unique to MDD did not (OR 1.00, 95% CI 0.95-1.06; P 0.97). Our replication sample was the Atherosclerosis Risk in Communities (ARIC) cohort of 7197 European-ancestry participants (1598 incident CAD cases). In ARIC, polygenic scores for MDD and loneliness were associated with hazard ratios of 1.07 (95% CI, 0.99-1.14; P = 0.07) and 1.07 (1.01-1.15; P = 0.03), respectively, and we replicated findings from the BioVU sensitivity analyses. We conclude that genetic risk factors for MDD and loneliness act pleiotropically to increase CAD risk in females.

Functional seizures are associated with cerebrovascular disease and functional stroke is more common in patients with functional seizures than epileptic seizures.

PURPOSE: To characterize the relationship between functional seizures (FSe), cerebrovascular disease (CVD), and functional stroke. METHOD: A retrospective case-control study of 189 patients at a single large tertiary medical center. We performed a manual chart review of medical records of patients with FSe or epileptic seizures (ES), who also had ICD code evidence of CVD. The clinical characteristics of FSe, ES, CVD, and functional stroke were recorded. Logistic regression and Welch's t-tests were used to evaluate the differences between the FSe and ES groups. RESULTS: Cerebrovascular disease was confirmed in 58.7% and 87.6% of patients with FSe or ES through manual chart review. Stroke was significantly more common in patients with ES (76.29%) than FSe (43.48%) (p = 4.07 × 10-6). However, compared to nonepileptic controls FSe was associated with both CVD (p < 0.0019) and stroke (p < 6.62 × 10-10). Functional stroke was significantly more common in patients with FSe (39.13%) than patients with ES (4.12%) (p = 4.47 × 10-9). Compared to patients with ES, patients with FSe were younger (p = 0.00022), more likely to be female (p = 0.00040), and more likely to have comorbid mental health needs including anxiety (p = 1.06 × 10-6), PTSD or history of trauma (e.g., sexual abuse) (p = 1.06 × 10-13), and bipolar disorder (p = 0.0011). CONCLUSION: Our results confirm the initial observation of increased CVD in patients with FSe and further suggest that patients with FSe may be predisposed to developing another functional neurological disorder (FND) (i.e., functional stroke). We speculate that this may be due to shared risk factors and pathophysiological processes that are common to various manifestations of FND.

The Evolutionary History of Common Genetic Variants Influencing Human Cortical Surface Area.

Structural brain changes along the lineage leading to modern Homo sapiens contributed to our distinctive cognitive and social abilities. However, the evolutionarily relevant molecular variants impacting key aspects of neuroanatomy are largely unknown. Here, we integrate evolutionary annotations of the genome at diverse timescales with common variant associations from large-scale neuroimaging genetic screens. We find that alleles with evidence of recent positive polygenic selection over the past 2000-3000 years are associated with increased surface area (SA) of the entire cortex, as well as specific regions, including those involved in spoken language and visual processing. Therefore, polygenic selective pressures impact the structure of specific cortical areas even over relatively recent timescales. Moreover, common sequence variation within human gained enhancers active in the prenatal cortex is associated with postnatal global SA. We show that such variation modulates the function of a regulatory element of the developmentally relevant transcription factor HEY2 in human neural progenitor cells and is associated with structural changes in the inferior frontal cortex. These results indicate that non-coding genomic regions active during prenatal cortical development are involved in the evolution of human brain structure and identify novel regulatory elements and genes impacting modern human brain structure.

Association of Preinjury Medical Diagnoses With Pediatric Persistent Postconcussion Symptoms in Electronic Health Records.

OBJECTIVE: To identify risk factors and generate hypotheses for pediatric persistent postconcussion symptoms (PPCS). SETTING: A regional healthcare system in the Southeastern United States. PARTICIPANTS: An electronic health record-based algorithm was developed and validated to identify PPCS cases and controls from an institutional database of more than 2.8 million patients. PPCS cases (n = 274) were patients aged 5 to 18 years with PPCS-related diagnostic codes or with PPCS key words identified by natural language processing of clinical notes. Age, sex, and year of index event-matched controls (n = 1096) were patients with mild traumatic brain injury codes only. Patients with moderate or severe traumatic brain injury were excluded. All patients used our healthcare system at least 3 times 180 days before their injury. DESIGN: Case-control study. MAIN MEASURES: The outcome was algorithmic classification of PPCS. Exposures were all preinjury medical diagnoses assigned at least 180 days before the injury. RESULTS: Cases and controls both had a mean of more than 9 years of healthcare system use preinjury. Of 221 preinjury medical diagnoses, headache disorder was associated with PPCS after accounting for multiple testing (odds ratio [OR] = 2.9; 95% confidence interval [CI]: 1.6-5.0; P = 2.1e-4). Six diagnoses were associated with PPCS at a suggestive threshold for statistical significance (false discovery rate P < .10): gastritis/duodenitis (OR = 2.8; 95% CI: 1.6-5.1; P = 5.0e-4), sleep disorders (OR = 2.3; 95% CI: 1.4-3.7; P = 7.4e-4), abdominal pain (OR = 1.6; 95% CI: 1.2-2.2; P = 9.2e-4), chronic sinusitis (OR = 2.8; 95% CI: 1.5-5.2; P = 1.3e-3), congenital anomalies of the skin (OR = 2.9; 95% CI: 1.5-5.5; P = 1.9e-3), and chronic pharyngitis/nasopharyngitis (OR = 2.4; 95% CI: 1.4-4.3; P = 2.5e-3). CONCLUSIONS: These results support the strong association of preinjury headache disorders with PPCS. An association of PPCS with prior gastritis/duodenitis, sinusitis, and pharyngitis/nasopharyngitis suggests a role for chronic inflammation in PPCS pathophysiology and risk, although results could equally be attributable to a higher likelihood of somatization among PPCS cases. Identified risk factors should be investigated further and potentially considered during the management of pediatric mild traumatic brain injury cases.