XWAS: A Software Toolset for Genetic Data Analysis and Association Studies of the X Chromosome

Methodological opportunities in genomic data analysis to advance health equity

The causes and consequences of inequities in genomic research and medicine are complex and widespread. However, it is widely acknowledged that underrepresentation of diverse populations in human genetics research risks exacerbating existing health disparities. Efforts to improve diversity are ongoing, but an often-overlooked source of inequity is the choice of analytical methods used to process, analyse and interpret genomic data. This choice can influence all areas of genomic research, from genome-wide association studies and polygenic score development to variant prioritization and functional genomics. New statistical and machine learning techniques to understand, quantify and correct for the impact of biases in genomic data are emerging within the wider genomic research and genomic medicine ecosystems. At this crucial time point, it is important to clarify where improvements in methods and practices can, or cannot, have a role in improving equity in genomics. Here, we review existing approaches to promote equity and fairness in statistical analysis for genomics, and propose future methodological developments that are likely to yield the most impact for equity.

Multi-ancestry GWAS reveals loci linked to human variation in LINE-1- and Alu-insertion numbers

LINE-1 (L1) and Alu are two families of transposable elements (TEs) occupying ~17% and ~11% of the human genome, respectively. Though only a small fraction of L1 copies is able to produce the machinery to mobilize autonomously, Alu and degenerate L1s can hijack their functional machinery and mobilize in trans. The expression and subsequent mobilization of L1 and Alu can exert pathological effects on their hosts. These features have made them promising focus subjects in studies of aging where they can become active. However, mechanisms regulating TE activity are incompletely characterized, especially in diverse human populations. To address these gaps, we leveraged genomic data from the 1000 Genomes Project to carry out a trans-ethnic GWAS of L1/Alu insertion singletons. These are rare, recently acquired insertions observed in only one person and which we used as proxies for variation in L1/Alu insertion numbers. Our approach identified SNVs in genomic regions containing genes with potential and known TE regulatory properties, and it enriched for SNVs in regions containing known regulators of L1 expression. Moreover, we identified reference TE copies and structural variants that associated with L1/Alu singletons, suggesting their potential contribution to TE insertion number variation. Finally, a transcriptional analysis of lymphoblastoid cells highlighted potential cell cycle alterations in a subset of samples harboring L1/Alu singletons. Collectively, our results suggest that known TE regulatory mechanisms may be active in diverse human populations, expand the list of loci implicated in TE insertion number variability, and reinforce links between TEs and disease.

Multi-ancestry GWAS reveals loci linked to human variation in LINE-1- and Alu-insertion numbers

LINE-1 (L1) and Alu are two families of transposable elements (TEs) occupying ~17% and ~11% of the human genome, respectively. Though only a small fraction of L1 copies is able to produce the machinery to mobilize autonomously, Alu and degenerate L1s can hijack their functional machinery and mobilize in trans. The expression and subsequent mobilization of L1 and Alu can exert pathological effects on their hosts. These features have made them promising focus subjects in studies of aging where they can become active. However, mechanisms regulating TE activity are incompletely characterized, especially in diverse human populations. To address these gaps, we leveraged genomic data from the 1000 Genomes Project to carry out a trans-ethnic GWAS of L1/Alu insertion singletons. These are rare, recently acquired insertions observed in only one person and which we used as proxies for variation in L1/Alu insertion numbers. Our approach identified SNVs in genomic regions containing genes with potential and known TE regulatory properties, and it enriched for SNVs in regions containing known regulators of L1 expression. Moreover, we identified reference TE copies and structural variants that associated with L1/Alu singletons, suggesting their potential contribution to TE insertion number variation. Finally, a transcriptional analysis of lymphoblastoid cells highlighted potential cell cycle alterations in a subset of samples harboring L1/Alu singletons. Collectively, our results suggest that known TE regulatory mechanisms may be active in diverse human populations, expand the list of loci implicated in TE insertion number variability, and reinforce links between TEs and disease.

The X chromosome's influences on the human brain

Genes on the X chromosome are extensively expressed in the human brain. However, little is known for the X chromosome's impact on the brain anatomy, microstructure, and functional networks. We examined 1045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X chromosome interactions while proposing an atlas outlining dosage compensation for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we found unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Large-scale audiometric phenotyping identifies distinct genes and pathways involved in hearing loss subtypes

Age-related hearing loss affects one-third of the population over 65 years. However, the diverse pathologies underlying these heterogenous phenotypes complicate genetic studies. To overcome challenges associated with accurate phenotyping for older adults with hearing loss, we applied computational phenotyping approaches based on audiometrically measured hearing loss. This novel phenotyping strategy uncovered distinct genetic variants associated with sensory and metabolic hearing loss. Sex-stratified analyses of these sexually dimorphic hearing loss phenotypes revealed a novel locus of relevance to sensory hearing loss in males, but not females. Enrichment analyses revealed that genes involved in frontotemporal dementia were implicated in metabolic hearing loss, while genes relating to sensory processing of sound by hair cells were implicated in sensory hearing loss. Our study has enhanced our understanding of these two distinct hearing loss phenotypes, representing the first step in the development of more precise treatments for these pathologically distinct hearing loss phenotypes.

Chromosome X-wide common variant association study in autism spectrum disorder

Autism spectrum disorder (ASD) displays a notable male bias in prevalence. Research into rare (<0.1) genetic variants on the X chromosome has implicated over 20 genes in ASD pathogenesis, such as MECP2, DDX3X, and DMD. The "female protective effect" in ASD suggests that females may require a higher genetic burden to manifest symptoms similar to those in males, yet the mechanisms remain unclear. Despite technological advances in genomics, the complexity of the biological nature of sex chromosomes leaves them underrepresented in genome-wide studies. Here, we conducted an X-chromosome-wide association study (XWAS) using whole-genome sequencing data from 6,873 individuals with ASD (82% males) across Autism Speaks MSSNG, Simons Simplex Collection (SSC), and Simons Powering Autism Research (SPARK), alongside 8,981 population controls (43% males). We analyzed 418,652 X chromosome variants, identifying 59 associated with ASD (p values 7.9 × 10-6 to 1.51 × 10-5), surpassing Bonferroni-corrected thresholds. Key findings include significant regions on Xp22.2 (lead SNP rs12687599, p = 3.57 × 10-7) harboring ASB9/ASB11 and another encompassing DDX53 and the PTCHD1-AS long non-coding RNA (lead SNP rs5926125, p = 9.47 × 10-6). When mapping genes within 10 kb of the 59 most significantly associated SNPs, 91 genes were found, 17 of which yielded association with ASD (GRPR, AP1S2, DDX53, HDAC8, PCDH19, PTCHD1, PCDH11X, PTCHD1-AS, DMD, SYAP1, CNKSR2, GLRA2, OFD1, CDKL5, GPRASP2, NXF5, and SH3KBP1). FGF13 emerged as an X-linked ASD candidate gene, highlighted by sex-specific differences in minor allele frequencies. These results reveal significant insights into X chromosome biology in ASD, confirming and nominating genes and pathways for further investigation.

Identifying X-chromosome variants associated with age-related macular degeneration

PURPOSE

In genome-wide association studies (GWAS), X chromosome (ChrX) variants are often not investigated. Sex-specific effects and ChrX-specific quality control (QC) are needed to examine these effects. Previous GWAS identified 52 autosomal variants associated with age-related macular degeneration (AMD) via the International AMD Genomics Consortium (IAMDGC), but did not analyze ChrX. Therefore¸ our goal was to investigate ChrX variants for association with AMD.

METHODS

We genotyped 29 629 non-Hispanic White (NHW) individuals (M/F:10404/18865; AMD12,087/14723) via a custom chip and imputed after ChrX-specific QC (XWAS 3.0) using the Michigan Imputation Server. Imputation generated 1 221 623 variants on ChrX. Age, informative PCs, and subphenotypes were covariates for logistic association analyses with Fisher's correction. Gene/pathway analyses were performed with VEGAS, GSEASNP, ICSNPathway, DAVID, and mirPath.

RESULTS

Logistic association on NHW individuals with sex correction identified variants in/near the genes SLITRK4, ARHGAP6, FGF13 and DMD associated with AMD (P < 1 × 10-6,Fisher's combined-corrected). Association testing of the subphenotypes of choroidal neovascularization and geographic atrophy (GA), identified variants in DMD associated with GA (P < 1 × 10-6, Fisher's combined-corrected). Via gene-based analysis with VEGAS, several genes were associated with AMD (P < 0.05, both truncated tail strength/truncated product P) including SLITRK4 and BHLHB9. Pathway analysis using GSEASNP and DAVID identified genes associated with nervous system development (FDR: P:0.02), and blood coagulation (FDR: P:0.03). Variants in the region of a microRNA (miR) were associated with AMD (P < 0.05, truncated tail strength/truncated product P). Via DIANA mirPath analysis, downstream targets of miRs showed association with brain disorders and fatty acid elongation (P < 0.05). A long noncoding RNA on ChrX near the DMD locus was also associated with AMD (P = 4 × 10-7). Epistatic analysis (t-statistic) for a quantitative trait of AMD vs control including covariates found a suggestive association in the XG gene (P = 2 × 10^-5).

CONCLUSIONS

Analysis of ChrX variation identifies several potential new locifor AMD risk and these variants nominate novel AMD pathways. Further analysis is needed to refine these results and to understand their biological significance and relationship with AMD development in worldwide populations.

Genetics of hypertension-related sex differences and hypertensive disorders of pregnancy

Background: Hypertension and hypertensive disorders of pregnancy (HDP) cause a significant burden of disease on societies and individuals by increasing cardiovascular disease risk. Environmental risk factors alone do not explain the observed sexual dimorphism in lifetime blood pressure (BP) trajectories nor inter-individual variation in HDP risk.

Methods: In this short review, we focus on the genetics of hypertension-related sex differences and HDP and discuss the importance of genetics utilization for sex-specific hypertension risk prediction.

Results: Population and twin studies estimate that 28-66% of variation in BP levels and HDP is explained by genetic variation, while genomic wide association studies suggest that BP traits and HDP partly share a common genetic background. Moreover, environmental and epigenetic regulation of these genes differ by sex and oestrogen receptors in particular are shown to convey cardio- and vasculoprotective effects through epigenetic regulation of DNA. The majority of known genetic variation in hypertension and HDP is polygenic. Polygenic risk scores for BP display stronger associations with hypertension risk in women than in men and are associated with sex-specific age of hypertension onset. Monogenic forms of hypertension are rare and mostly present equally in both sexes.

Conclusion: Despite recent genetic discoveries providing new insights into HDP and sex differences in BP traits, further research is needed to elucidate the underlying biology. Emphasis should be placed on demonstrating the added clinical value of these genetic discoveries, which may eventually facilitate genomics-based personalized treatments for hypertension and HDP.

Genetic analysis of cognitive preservation in the midwestern Amish reveals a novel locus on chromosome 2

INTRODUCTION

Alzheimer's disease (AD) remains a debilitating condition with limited treatments and additional therapeutic targets needed. Identifying AD protective genetic loci may identify new targets and accelerate identification of therapeutic treatments. We examined a founder population to identify loci associated with cognitive preservation into advanced age.

METHODS

Genome-wide association and linkage analyses were performed on 946 examined and sampled Amish individuals, aged 76-95, who were either cognitively unimpaired (CU) or impaired (CI).

RESULTS

A total of 12 single nucleotide polymorphisms (SNPs) demonstrated suggestive association (P ≤ 5 × 10-4) with cognitive preservation. Genetic linkage analyses identified > 100 significant (logarithm of the odds [LOD] ≥ 3.3) SNPs, some which overlapped with the association results. Only one locus on chromosome 2 retained significance across multiple analyses.

DISCUSSION

A novel significant result for cognitive preservation on chromosome 2 includes the genes LRRTM4 and CTNNA2. Additionally, the lead SNP, rs1402906, impacts the POU3F2 transcription factor binding affinity, which regulates LRRTM4 and CTNNA2.

HIGHLIGHTS

GWAS and linkage identified over 100 loci associated with cognitive preservation. One locus on Chromosome 2 retained significance over multiple analyses. Predicted TFBSs near rs1402906 regulate genes associated with neurocognition.

Chromosome X-Wide Common Variant Association Study (XWAS) in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) displays a notable male bias in prevalence. Research into rare (<0.1) genetic variants on the X chromosome has implicated over 20 genes in ASD pathogenesis, such as MECP2, DDX3X, and DMD. The "female protective effect" in ASD suggests that females may require a higher genetic burden to manifest similar symptoms as males, yet the mechanisms remain unclear. Despite technological advances in genomics, the complexity of the biological nature of sex chromosomes leave them underrepresented in genome-wide studies. Here, we conducted an X chromosome-wide association study (XWAS) using whole-genome sequencing data from 6,873 individuals with ASD (82% males) across Autism Speaks MSSNG, Simons Simplex Cohort SSC, and Simons Foundation Powering Autism Research SPARK, alongside 8,981 population controls (43% males). We analyzed 418,652 X-chromosome variants, identifying 59 associated with ASD (p-values 7.9×10-6 to 1.51×10-5), surpassing Bonferroni-corrected thresholds. Key findings include significant regions on chrXp22.2 (lead SNP=rs12687599, p=3.57×10-7) harboring ASB9/ASB11, and another encompassing DDX53/PTCHD1-AS long non-coding RNA (lead SNP=rs5926125, p=9.47×10-6). When mapping genes within 10kb of the 59 most significantly associated SNPs, 91 genes were found, 17 of which yielded association with ASD (GRPR, AP1S2, DDX53, HDAC8, PCDH19, PTCHD1, PCDH11X, PTCHD1-AS, DMD, SYAP1, CNKSR2, GLRA2, OFD1, CDKL5, GPRASP2, NXF5, SH3KBP1). FGF13 emerged as a novel X-linked ASD candidate gene, highlighted by sex-specific differences in minor allele frequencies. These results reveal significant new insights into X chromosome biology in ASD, confirming and nominating genes and pathways for further investigation.

Leveraging sex-genetic interactions to understand brain disorders: recent advances and current gaps

It is established that there are sex differences in terms of prevalence, age of onset, clinical manifestations, and response to treatment for a variety of brain disorders, including neurodevelopmental, psychiatric, and neurodegenerative disorders. Cohorts of increasing sample sizes with diverse data types collected, including genetic, transcriptomic and/or phenotypic data, are providing the building blocks to permit analytical designs to test for sex-biased genetic variant-trait associations, and for sex-biased transcriptional regulation. Such molecular assessments can contribute to our understanding of the manifested phenotypic differences between the sexes for brain disorders, offering the future possibility of delivering personalized therapy for females and males. With the intention of raising the profile of this field as a research priority, this review aims to shed light on the importance of investigating sex-genetic interactions for brain disorders, focusing on two areas: (i) variant-trait associations and (ii) transcriptomics (i.e. gene expression, transcript usage and regulation). We specifically discuss recent advances in the field, current gaps and provide considerations for future studies.

An eQTL-based approach reveals candidate regulators of LINE-1 RNA levels in lymphoblastoid cells

Long interspersed element 1 (LINE-1; L1) are a family of transposons that occupy ~17% of the human genome. Though a small number of L1 copies remain capable of autonomous transposition, the overwhelming majority of copies are degenerate and immobile. Nevertheless, both mobile and immobile L1s can exert pleiotropic effects (promoting genome instability, inflammation, or cellular senescence) on their hosts, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSD17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for regulators of transposon RNA levels.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Mechanisms of sex differences in Alzheimer's disease

Alzheimer's disease (AD) and the mechanisms underlying its etiology and progression are complex and multifactorial. The higher AD risk in women may serve as a clue to better understand these complicated processes. In this review, we examine aspects of AD that demonstrate sex-dependent effects and delve into the potential biological mechanisms responsible, compiling findings from advanced technologies such as single-cell RNA sequencing, metabolomics, and multi-omics analyses. We review evidence that sex hormones and sex chromosomes interact with various disease mechanisms during aging, encompassing inflammation, metabolism, and autophagy, leading to unique characteristics in disease progression between men and women.

Comparing the effect of imputation reference panel composition in four distinct Latin American cohorts

Genome-wide association studies have been useful in identifying genetic risk factors for various phenotypes. These studies rely on imputation and many existing panels are largely composed of individuals of European ancestry, resulting in lower levels of imputation quality in underrepresented populations. We aim to analyze how the composition of imputation reference panels affects imputation quality in four target Latin American cohorts. We compared imputation quality for chromosomes 7 and X when altering the imputation reference panel by: 1) increasing the number of Latin American individuals; 2) excluding either Latin American, African, or European individuals, or 3) increasing the Indigenous American (IA) admixture proportions of included Latin Americans. We found that increasing the number of Latin Americans in the reference panel improved imputation quality in the four populations; however, there were differences between chromosomes 7 and X in some cohorts. Excluding Latin Americans from analysis resulted in worse imputation quality in every cohort, while differential effects were seen when excluding Europeans and Africans between and within cohorts and between chromosomes 7 and X. Finally, increasing IA-like admixture proportions in the reference panel increased imputation quality at different levels in different populations. The difference in results between populations and chromosomes suggests that existing and future reference panels containing Latin American individuals are likely to perform differently in different Latin American populations.

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

The genetic architecture of youth anxiety: a study protocol

BACKGROUND

Anxiety disorders are the most common psychiatric problems among Canadian youth and typically have an onset in childhood or adolescence. They are characterized by high rates of relapse and chronicity, often resulting in substantial impairment across the lifespan. Genetic factors play an important role in the vulnerability toward anxiety disorders. However, genetic contribution to anxiety in youth is not well understood and can change across developmental stages. Large-scale genetic studies of youth are needed with detailed assessments of symptoms of anxiety disorders and their major comorbidities to inform early intervention or preventative strategies and suggest novel targets for therapeutics and personalization of care.

METHODS

The Genetic Architecture of Youth Anxiety (GAYA) study is a Pan-Canadian effort of clinical and genetic experts with specific recruitment sites in Calgary, Halifax, Hamilton, Toronto, and Vancouver. Youth aged 10-19 (n = 13,000) will be recruited from both clinical and community settings and will provide saliva samples, complete online questionnaires on demographics, symptoms of mental health concerns, and behavioural inhibition, and complete neurocognitive tasks. A subset of youth will be offered access to a self-managed Internet-based cognitive behavioral therapy resource. Analyses will focus on the identification of novel genetic risk loci for anxiety disorders in youth and assess how much of the genetic risk for anxiety disorders is unique or shared across the life span.

DISCUSSION

Results will substantially inform early intervention or preventative strategies and suggest novel targets for therapeutics and personalization of care. Given that the GAYA study will be the biggest genomic study of anxiety disorders in youth in Canada, this project will further foster collaborations nationally and across the world.

Genetic analysis of the X chromosome in people with Lewy body dementia nominates new risk loci

Sex influences the prevalence and symptoms of Lewy body dementia (LBD). However, genome-wide association studies typically focus on autosomal variants and exclude sex-specific risk factors. We addressed this gap by performing an X chromosome-wide association study using whole-genome sequence data from 2591 LBD cases and 4391 controls. We identified a significant risk locus within intron 1 of MAP3K15 (rs141773145, odds ratio = 2.42, 95% confidence interval = 1.65-3.56, p-value = 7.0 × 10-6) in female LBD cases conditioned for APOE ε4 dosage. The locus includes an enhancer region that regulates MAP3K15 expression in ganglionic eminence cells derived from primary cultured neurospheres. Rare variant burden testing showed differential enrichment of missense mutations in TEX13A in female LBD cases, that did not reach significance (p-value = 1.34 × 10-4). These findings support the sex-specific effects of genetic factors and a potential role of Alzheimer's-related risk for females with LBD.

Effects of sex and gender on the etiologies and presentation of select internalizing psychopathologies

The internalizing spectrum encompasses a subset of psychopathologies characterized by emotional liability, anhedonia, anxiousness, distress, and fear, and includes, among others, diagnoses of major depressive disorder (MDD), generalized anxiety disorder (GAD), and posttraumatic stress disorder (PTSD). In this review, we describe the vast body of work highlighting a role for sex and gender in the environment, symptom onset, genetic liability, and disorder progression and comorbidities of MDD, GAD, and PTSD. We also point the reader to different language used in diverse fields to describe sexual and gender minorities that may complicate the interpretation of emerging literature from the social sciences, psychiatric and psychological sciences, and genetics. Finally, we identify several gaps in knowledge that we hope serve as launch-points for expanding the scope of psychiatric studies beyond binarized sex-stratification. Despite being under-represented in genomics studies, placing emphasis on inclusion of sexual and gender diverse participants in these works will hopefully improve our understanding of disorder etiology using genetics as one tool to inform how biology (e.g., hormone concentration) and environmental variables (e.g., exposure to traumatic events) contribute to differences in symptom onset, pattern, and long-term trajectory.

Genetic analysis of cognitive preservation in the midwestern Amish reveals a novel locus on chromosome 2

INTRODUCTION

Alzheimer disease (AD) remains a debilitating condition with limited treatments and additional therapeutic targets needed. Identifying AD protective genetic loci may identify new targets and accelerate identification of therapeutic treatments. We examined a founder population to identify loci associated with cognitive preservation into advanced age.

METHODS

Genome-wide association and linkage analyses were performed on 946 examined and sampled Amish individuals, aged 76-95, who were either cognitively unimpaired (CU) or impaired (CI).

RESULTS

12 SNPs demonstrated suggestive association (P≤5×10-4) with cognitive preservation. Genetic linkage analyses identified >100 significant (LOD≥3.3) SNPs, some which overlapped with the association results. Only one locus on chromosome 2 retained significance across multiple analyses.

DISCUSSION

A novel significant result for cognitive preservation on chromosome 2 includes the genes LRRTM4 and CTNNA2. Additionally, the lead SNP, rs1402906, impacts the POU3F2 transcription factor binding affinity, which regulates LRRTM4 and CTNNA2.

Genome-wide association study of antisocial personality disorder diagnostic criteria provides evidence for shared risk factors across disorders

INTRODUCTION

While progress has been made in determining the genetic basis of antisocial behaviour, little progress has been made for antisocial personality disorder (ASPD), a condition that often co-occurs with other psychiatric conditions including substance use disorders, attention deficit hyperactivity disorder (ADHD), and anxiety disorders. This study aims to improve the understanding of the genetic risk for ASPD and its relationship with other disorders and traits.

METHODS

We conducted a genome-wide association study (GWAS) of the number of ASPD diagnostic criteria data from 3217 alcohol-dependent participants recruited in the UK (UCL, N = 644) and the USA (Yale-Penn, N = 2573).

RESULTS

We identified rs9806493, a chromosome 15 variant, that showed a genome-wide significant association ( Z -score = -5.501, P = 3.77 × 10 -8 ) with ASPD criteria. rs9806493 is an eQTL for SLCO3A1 (Solute Carrier Organic Anion Transporter Family Member 3A1), a ubiquitously expressed gene with strong expression in brain regions that include the anterior cingulate and frontal cortices. Polygenic risk score analysis identified positive correlations between ASPD and smoking, ADHD, depression traits, and posttraumatic stress disorder. Negative correlations were observed between ASPD PRS and alcohol intake frequency, reproductive traits, and level of educational attainment.

CONCLUSION

This study provides evidence for an association between ASPD risk and SLCO3A1 and provides insight into the genetic architecture and pleiotropic associations of ASPD.

An eQTL-based Approach Reveals Candidate Regulators of LINE-1 RNA Levels in Lymphoblastoid Cells

Long interspersed element 1 (L1) are a family of autonomous, actively mobile transposons that occupy ~17% of the human genome. A number of pleiotropic effects induced by L1 (promoting genome instability, inflammation, or cellular senescence) have been observed, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSDB17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for transposon transcriptional regulators.

Identification of candidate novel production variants on the Bos taurus chromosome X

Chromosome X is often excluded from bovine genetic studies due to complications caused by the sex specific nature of the chromosome. As chromosome X is the second largest cattle chromosome and makes up approximately 6% of the female genome, finding ways to include chromosome X in dairy genetic studies is important. Using female animals and treating chromosome X as an autosome, we performed X chromosome inclusive genome-wide association studies in the selective breeding environment of the New Zealand dairy industry, aiming to identify chromosome X variants associated with milk production traits. We report on the findings of these genome-wide association studies and their potential effect within the dairy industry. We identify missense mutations in the MOSPD1 and CCDC160 genes that are associated with decreased milk volume and protein production and increased fat production. Both of these mutations are exonic SNP that are more prevalent in the Jersey breed than in Holstein-Friesians. Of the 2 candidates proposed it is likely that only one is causal, though we have not been able to identify which is more likely.

Sex-Dimorphic Analyses Identify Novel and Sex-Specific Genetic Associations in Inflammatory Bowel Disease

BACKGROUND

Sex is an integral variable often overlooked in complex disease genetics. Differences between sexes have been reported in natural history, disease complications, and age of onset in inflammatory bowel disease (IBD). While association studies have identified >230 IBD loci, there have been a limited number of studies investigating sex differences underlying these genetic associations.

METHODS

We report the first investigation of sex-dimorphic associations via meta-analysis of a sex-stratified association study (34 579 IBD cases, 39 125 controls). In addition, we performed chromosome (chr) X-specific analyses, considering models of X inactivation (XCI) and XCI escape. Demographic and clinical characteristics were also compared between sexes.

RESULTS

We identified significant differences between sexes for disease location and perianal complication in Crohn's disease and disease extent in ulcerative colitis. We observed genome-wide-significant sex-dimorphic associations (P < 5 × 10-8) at loci not previously reported in large-scale IBD genetic studies, including at chr9q22, CARMIL1, and UBASH3A. We identified variants in known IBD loci, including in chr2p15 and within the major histocompatibility complex on chr6, exhibiting sex-specific patterns of association (P < 5 × 10-7 in one sex only). We identified 3 chrX associations with IBD, including a novel Crohn's disease susceptibility locus at Xp22.

CONCLUSIONS

These analyses identified novel IBD loci, in addition to characterizing sex-specific patterns of associations underlying sex-dimorphic associations. By elucidating the role of sex in IBD genetics, our study will help enhance our understanding of the differences between the sexes in IBD biology and underscores a need to move beyond conventional sex-combined analyses to appreciate the genetic architecture of IBD more comprehensively.

eXclusionarY: 10 years later, where are the sex chromosomes in GWASs?

10 years ago, a detailed analysis showed that only 33% of genome-wide association study (GWAS) results included the X chromosome. Multiple recommendations were made to combat such exclusion. Here, we re-surveyed the research landscape to determine whether these earlier recommendations had been translated. Unfortunately, among the genome-wide summary statistics reported in 2021 in the NHGRI-EBI GWAS Catalog, only 25% provided results for the X chromosome and 3% for the Y chromosome, suggesting that the exclusion phenomenon not only persists but has also expanded into an exclusionary problem. Normalizing by physical length of the chromosome, the average number of studies published through November 2022 with genome-wide-significant findings on the X chromosome is ∼1 study/Mb. By contrast, it ranges from ∼6 to ∼16 studies/Mb for chromosomes 4 and 19, respectively. Compared with the autosomal growth rate of ∼0.086 studies/Mb/year over the last decade, studies of the X chromosome grew at less than one-seventh that rate, only ∼0.012 studies/Mb/year. Among the studies that reported significant associations on the X chromosome, we noted extreme heterogeneities in data analysis and reporting of results, suggesting the need for clear guidelines. Unsurprisingly, among the 430 scores sampled from the PolyGenic Score Catalog, 0% contained weights for sex chromosomal SNPs. To overcome the dearth of sex chromosome analyses, we provide five sets of recommendations and future directions. Finally, until the sex chromosomes are included in a whole-genome study, instead of GWASs, we propose such studies would more properly be referred to as "AWASs," meaning "autosome-wide scans."

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.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

An Efficient Bayesian Method for Estimating the Degree of the Skewness of X Chromosome Inactivation Based on the Mixture of General Pedigrees and Unrelated Females

Skewed X chromosome inactivation (XCI-S) has been reported to be associated with some X-linked diseases. Several methods have been proposed to estimate the degree of XCI-S (denoted as γ) for quantitative and qualitative traits based on unrelated females. However, there is no method available for estimating γ based on general pedigrees. Therefore, in this paper, we propose a Bayesian method to obtain the point estimate and the credible interval of γ based on the mixture of general pedigrees and unrelated females (called mixed data for brevity), which is also suitable for only general pedigrees. We consider the truncated normal prior and the uniform prior for γ. Further, we apply the eigenvalue decomposition and Cholesky decomposition to our proposed methods to accelerate the computation speed. We conduct extensive simulation studies to compare the performances of our proposed methods and two existing Bayesian methods which are only applicable to unrelated females. The simulation results show that the incorporation of general pedigrees can improve the efficiency of the point estimation and the precision and the accuracy of the interval estimation of γ. Finally, we apply the proposed methods to the Minnesota Center for Twin and Family Research data for their practical use.

Sex-Specific Genetic and Transcriptomic Liability to Neuroticism

BACKGROUND

The presentation, etiology, and relative risk of psychiatric disorders are strongly influenced by biological sex. Neuroticism is a transdiagnostic feature of psychiatric disorders displaying prominent sex differences. We performed genome-wide association studies of neuroticism separately in males and females to identify sex-specific genetic and transcriptomic profiles.

METHODS

Neuroticism scores were derived from the Eysenck Personality Inventory Neuroticism scale. Genome-wide association studies were performed in 145,669 females and 129,229 males from the UK Biobank considering autosomal and X chromosomal variation. Two-sided z tests were used to test for sex-specific effects of discovered loci, genetic correlates (n = 673 traits), tissue and gene transcriptomic profiles, and polygenic associations across health outcomes in the Vanderbilt University Biobank (39,692 females and 31,268 males).

RESULTS

The single nucleotide polymorphism heritability of neuroticism was not statistically different between males (h2 = 10.6%) and females (h2 = 11.85%). Four female-specific (rs10736549-CNTN5, rs6507056-ASXL3, rs2087182-MMS22L, and rs72995548-HSPB2) and 2 male-specific (rs10507274-MED13L and rs7984597) neuroticism risk loci reached genome-wide significance. Male- and female-specific neuroticism polygenic scores were most significantly associated with mood disorders (males: odds ratio = 1.11, p = 1.40 × 10-9; females: odds ratio = 1.14, p = 6.05 × 10-22). They also associated with sex-specific laboratory measurements related to erythrocyte count, distribution, and hemoglobin concentration. Gene expression variation in the pituitary was enriched for neuroticism loci in males (male: b = 0.026, p = .002), and genetically regulated transcriptomic changes highlighted the effect of SHISHA9, TEX26, and NCOA6.

CONCLUSIONS

Through a comprehensive assessment of genetic risk for neuroticism and the associated biological processes, this study identified several molecular pathways that can partially explain the known sex differences in neurotic symptoms and their psychiatric comorbidities.

Why does the X chromosome lag behind autosomes in GWAS findings?

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

Altered gene expression in excitatory neurons is associated with Alzheimer's disease and its higher incidence in women

Introduction

Alzheimer's disease (AD) is a neurodegenerative disorder involving interactions between different cell types in the brain. Previous single-cell and bulk expression Alzheimer's studies have reported conflicting findings about the key cell types and cellular pathways whose expression is primarily altered in this disease. We re-analyzed these data in a uniform, coherent manner aiming to resolve and extend past findings. Our analysis sheds light on the observation that females have higher AD incidence than males.

Methods

We re-analyzed three single-cell transcriptomics datasets. We used the software Model-based Analysis of Single-cell Transcriptomics (MAST) to seek differentially expressed genes comparing AD cases to matched controls for both sexes together and each sex separately. We used the GOrilla software to search for enriched pathways among the differentially expressed genes. Motivated by the male/female difference in incidence, we studied genes on the X-chromosome, focusing on genes in the pseudoautosomal region (PAR) and on genes that are heterogeneous across individuals or tissues for X-inactivation. We validated findings by analyzing bulk AD datasets from the cortex in the Gene Expression Omnibus.

Results

Our results resolve a contradiction in the literature, showing that by comparing AD patients to unaffected controls, excitatory neurons have more differentially expressed genes than do other cell types. Synaptic transmission and related pathways are altered in a sex-specific analysis of excitatory neurons. PAR genes and X-chromosome heterogeneous genes, including, for example, BEX1 and ELK1, may contribute to the difference in sex incidence of Alzheimer's disease. GRIN1, stood out as an overexpressed autosomal gene in cases versus controls in all three single-cell datasets and as a functional candidate gene contributing to pathways upregulated in cases.

Discussion

Taken together, these results point to a potential linkage between two longstanding questions concerning AD pathogenesis, involving which cell type is the most important and why females have a higher incidence than males.

Highlights

By reanalyzing three, published, single-cell RNAseq datasets, we resolved a contradiction in the literature and showed that when comparing AD patients to unaffected controls, excitatory neurons have more differentially expressed genes than do other cell types.Further analysis of the published single-cell datasets showed that synaptic transmission and related pathways are altered in a sex-specific analysis of excitatory neurons.Combining analysis of single-cell datasets and publicly available bulk transcriptomics datasets revealed that X-chromosome genes, such as BEX1, ELK1, and USP11, whose X-inactivation status is heterogeneous may contribute to the higher incidence in females of Alzheimer's disease.

A systematic review of analytical methods used in genetic association analysis of the X-chromosome

Genetic association studies have been very successful at elucidating the genetic background of many complex diseases/traits. However, the X-chromosome is often neglected in these studies because of technical difficulties and the fact that most tools only utilize genetic data from autosomes. In this review, we aim to provide an overview of different practical approaches that are followed to incorporate the X-chromosome in association analysis, such as Genome-Wide Association Studies and Expression Quantitative Trait Loci Analysis. In general, the choice of which test statistics is most appropriate will depend on three main criteria: (1) the underlying X-inactivation model, (2) if Hardy-Weinberg equilibrium holds and sex-specific allele frequencies are expected and (3) whether adjustment for confounding variables is required. All in all, it is recommended that a combination of different association tests should be used for the analysis of X-chromosome.

Illuminating the Mechanisms Underlying Sex Differences in Cardiovascular Disease

Sex is a key risk factor for many types of cardiovascular disease. It is imperative to understand the mechanisms underlying sex differences to devise optimal preventive and therapeutic approaches for all individuals. Both biological sex (determined by sex chromosomes and gonadal hormones) and gender (social and cultural behaviors associated with femininity or masculinity) influence differences between men and women in disease susceptibility and pathology. Here, we focus on the application of experimental mouse models that elucidate the influence of 2 components of biological sex-sex chromosome complement (XX or XY) and gonad type (ovaries or testes). These models have revealed that in addition to well-known effects of gonadal hormones, sex chromosome complement influences cardiovascular risk factors, such as plasma cholesterol levels and adiposity, as well as the development of atherosclerosis and pulmonary hypertension. One mechanism by which sex chromosome dosage influences cardiometabolic traits is through sex-biased expression of X chromosome genes that escape X inactivation. These include chromatin-modifying enzymes that regulate gene expression throughout the genome. The identification of factors that determine sex-biased gene expression and cardiometabolic traits will expand our mechanistic understanding of cardiovascular disease processes and provide insight into sex differences that remain throughout the lifespan as gonadal hormone levels alter with age.

BEXCIS: Bayesian methods for estimating the degree of the skewness of X chromosome inactivation

BACKGROUND

X chromosome inactivation (XCI) is an epigenetic phenomenon that one of two X chromosomes in females is transcriptionally silenced during early embryonic development. Skewed XCI has been reported to be associated with some X-linked diseases. There have been several methods measuring the degree of the skewness of XCI. However, these methods may still have several limitations.

RESULTS

We propose a Bayesian method to obtain the point estimate and the credible interval of the degree of XCI skewing by incorporating its prior information of being between 0 and 2. We consider a normal prior and a uniform prior for it (respectively denoted by BN and BU). We also propose a penalized point estimate based on the penalized Fieller's method and derive the corresponding confidence interval. Simulation results demonstrate that the BN and BU methods can solve the problems of extreme point estimates, noninformative intervals, empty sets and discontinuous intervals. The BN method generally outperforms other methods with the lowest mean squared error in the point estimation, and well controls the coverage probability with the smallest median and the least variation of the interval width in the interval estimation. We apply all the methods to the Graves' disease data and the Minnesota Center for Twin and Family Research data, and find that SNP rs3827440 in the Graves' disease data may undergo skewed XCI towards the allele C.

CONCLUSIONS

We recommend the BN method for measuring the degree of the skewness of XCI in practice. The R package BEXCIS is publicly available at https://github.com/Wen-YiYu/BEXCIS .

Representing sex chromosomes in genome assemblies

Sex chromosomes have evolved hundreds of independent times across eukaryotes. As genome sequencing, assembly, and scaffolding techniques rapidly improve, it is now feasible to build fully phased sex chromosome assemblies. Despite technological advances enabling phased assembly of whole chromosomes, there are currently no standards for representing sex chromosomes when publicly releasing a genome. Furthermore, most computational analysis tools are unable to efficiently investigate their unique biology relative to autosomes. We discuss a diversity of sex chromosome systems and consider the challenges of representing sex chromosome pairs in genome assemblies. By addressing these issues now as technologies for full phasing of chromosomal assemblies are maturing, we can collectively ensure that future genome analysis toolkits can be broadly applied to all eukaryotes with diverse types of sex chromosome systems. Here we provide best practice guidelines for presenting a genome assembly that contains sex chromosomes. These guidelines can also be applied to other non-recombining genomic regions, such as S-loci in plants and mating-type loci in fungi and algae.

XCMAX4: A Robust X Chromosomal Genetic Association Test Accounting for Covariates

Although the X chromosome accounts for about 5% of the human genes, it is routinely excluded from genome-wide association studies probably due to its unique structure and complex biological patterns. While some statistical methods have been proposed for testing the association between X chromosomal markers and diseases, very a few of them can adjust for covariates. Unfortunately, those methods that can incorporate covariates either need to specify an X chromosome inactivation model or require the permutation procedure to compute the p value. In this article, we proposed a novel analytic approach based on logistic regression that allows for covariates and does not need to specify the underlying X chromosome inactivation pattern. Simulation studies showed that our proposed method controls the size well and has robust performance in power across various practical scenarios. We applied the proposed method to analyze Graves’ disease data to show its usefulness in practice.

Gene-Based Methods for Estimating the Degree of the Skewness of X Chromosome Inactivation

Skewed X chromosome inactivation (XCI-S) has been reported to be associated with some X-linked diseases, and currently several methods have been proposed to estimate the degree of the XCI-S (denoted as γ) for a single locus. However, no method has been available to estimate γ for genes. Therefore, in this paper, we first propose the point estimate and the penalized point estimate of γ for genes, and then derive its confidence intervals based on the Fieller’s and penalized Fieller’s methods, respectively. Further, we consider the constraint condition of γ∈[0, 2] and propose the Bayesian methods to obtain the point estimates and the credible intervals of γ, where a truncated normal prior and a uniform prior are respectively used (denoted as GBN and GBU). The simulation results show that the Bayesian methods can avoid the extreme point estimates (0 or 2), the empty sets, the noninformative intervals ([0, 2]) and the discontinuous intervals to occur. GBN performs best in both the point estimation and the interval estimation. Finally, we apply the proposed methods to the Minnesota Center for Twin and Family Research data for their practical use. In summary, in practical applications, we recommend using GBN to estimate γ of genes.

Major sex differences in allele frequencies for X chromosomal variants in both the 1000 Genomes Project and gnomAD

An unexpectedly high proportion of SNPs on the X chromosome in the 1000 Genomes Project phase 3 data were identified with significant sex differences in minor allele frequencies (sdMAF). sdMAF persisted for many of these SNPs in the recently released high coverage whole genome sequence of the 1000 Genomes Project that was aligned to GRCh38, and it was consistent between the five super-populations. Among the 245,825 common (MAF>5%) biallelic X-chromosomal SNPs in the phase 3 data presumed to be of high quality, 2,039 have genome-wide significant sdMAF (p-value <5e-8). sdMAF varied by location: non-pseudo-autosomal region (NPR) = 0.83%, pseudo-autosomal regions (PAR1) = 0.29%, PAR2 = 13.1%, and X-transposed region (XTR)/PAR3 = 0.85% of SNPs had sdMAF, and they were clustered at the NPR-PAR boundaries, among others. sdMAF at the NPR-PAR boundaries are biologically expected due to sex-linkage, but have generally been ignored in association studies. For comparison, similar analyses found only 6, 1 and 0 SNPs with significant sdMAF on chromosomes 1, 7 and 22, respectively. Similar sdMAF results for the X chromosome were obtained from the high coverage whole genome sequence data from gnomAD V 3.1.2 for both the non-Finnish European and African/African American samples. Future X chromosome analyses need to take sdMAF into account.

The Role of Epigenetics in Primary Biliary Cholangitis

Primary Biliary Cholangitis (PBC) is a rare autoimmune disease of the liver, affecting mostly females. There is evidence that epigenetic changes have a pathogenic role in PBC. Epigenetic modifications are related to methylation of CpG DNA islands, post-translational modifications of histone proteins, and non-coding RNAs. In PBC, there are data showing a dysregulation of all these levels, especially in immune cells. In addition, epigenetics seems to be involved in complex phenomena such as X monosomy or abnormalities in the process of X chromosome inactivation, which have been reported in PBC and appear to influence its sex imbalance and pathogenesis. We review here historical data on epigenetic modifications in PBC, present new data, and discuss possible links among X-chromosome abnormalities at a genetic and epigenetic level, PBC pathogenesis, and PBC sex imbalance.

Sex-linked genetic mechanisms and atrial fibrillation risk

Atrial fibrillation (AF) is a cardiac condition characterised by an irregular heartbeat, atrial pathology and an elevated downstream risk of thrombosis and heart failure, as well as neurological sequelae including stroke and dementia. The prevalence and presentation of, risk factors for, and therapeutic responses to, AF differ by sex, and this sex bias may be partially explained in terms of genetics. Here, we consider four sex-linked genetic mechanisms that may influence sex-biased phenotypes related to AF and provide examples of each: X-linked gene dosage, X-linked genomic imprinting, sex-biased autosomal gene expression, and male-limited Y-linked gene expression. We highlight novel candidate risk genes and pathways that warrant further investigation in clinical and preclinical studies. Understanding the biological basis of sex differences in AF should allow better prediction of disease risk, identification of novel risk/protective factors, and the development of more effective sex-tailored interventions.

Genetic, Molecular, and Cellular Determinants of Sex-Specific Cardiovascular Traits

Despite the well-known sex dimorphism in cardiovascular disease traits, the exact genetic, molecular, and cellular underpinnings of these differences are not well understood. A growing body of evidence currently points at the links between cardiovascular disease traits and the genome, epigenome, transcriptome, and metabolome. However, the sex-specific differences in these links remain largely unstudied due to challenges in bioinformatic methods, inadequate statistical power, analytic costs, and paucity of valid experimental models. This review article provides an overview of the literature on sex differences in genetic architecture, heritability, epigenetic changes, transcriptomic signatures, and metabolomic profiles in relation to cardiovascular disease traits. We also review the literature on the associations between sex hormones and cardiovascular disease traits and discuss the potential mechanisms underlying these associations, focusing on human studies.

Sex Differences in Neurodegeneration: The Role of the Immune System in Humans

Growing evidence supports significant involvement of immune dysfunction in the etiology of neurodegenerative diseases, several of which also display prominent sex differences across prevalence, pathology, and symptomology. In this review, we summarize evidence from human studies of established and recent findings of sex differences in multiple sclerosis, Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis and discuss how sex-specific central nervous system innate immune activity could contribute to downstream sex differences in these diseases. We examine human genomic and transcriptomics studies in each neurodegenerative disease through the lens of sex differences in the neuroimmune system and highlight the importance of stratifying sex in clinical and translational research studies. Finally, we discuss the limitations of the existing studies and outline recommendations for further advancing sex-based analyses to uncover novel disease mechanisms that could ultimately help treat both sexes.

MEDTEC Students against Coronavirus: Investigating the Role of Hemostatic Genes in the Predisposition to COVID-19 Severity

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiologic agent of the coronavirus disease 2019 (COVID-19) pandemic. Besides virus intrinsic characteristics, the host genetic makeup is predicted to account for the extreme clinical heterogeneity of the disease, which is characterized, among other manifestations, by a derangement of hemostasis associated with thromboembolic events. To date, large-scale studies confirmed that genetic predisposition plays a role in COVID-19 severity, pinpointing several susceptibility genes, often characterized by immunologic functions. With these premises, we performed an association study of common variants in 32 hemostatic genes with COVID-19 severity. We investigated 49,845 single-nucleotide polymorphism in a cohort of 332 Italian severe COVID-19 patients and 1668 controls from the general population. The study was conducted engaging a class of students attending the second year of the MEDTEC school (a six-year program, held in collaboration between Humanitas University and the Politecnico of Milan, allowing students to gain an MD in Medicine and a Bachelor's Degree in Biomedical Engineering). Thanks to their willingness to participate in the fight against the pandemic, we evidenced several suggestive hits (p < 0.001), involving the PROC, MTHFR, MTR, ADAMTS13, and THBS2 genes (top signal in PROC: chr2:127192625:G:A, OR = 2.23, 95%CI = 1.50-3.34, p = 8.77 × 10-5). The top signals in PROC, MTHFR, MTR, ADAMTS13 were instrumental for the construction of a polygenic risk score, whose distribution was significantly different between cases and controls (p = 1.62 × 10-8 for difference in median levels). Finally, a meta-analysis performed using data from the Regeneron database confirmed the contribution of the MTHFR variant chr1:11753033:G:A to the predisposition to severe COVID-19 (pooled OR = 1.21, 95%CI = 1.09-1.33, p = 4.34 × 10-14 in the weighted analysis).

Improving clinical outcomes through attention to sex and hormones in research

Biological sex, fluctuations in sex steroid hormones throughout life and gender as a social construct all influence every aspect of health and disease. Yet, for decades, most basic and clinical studies have included only male individuals. As modern health care moves towards personalized medicine, it is clear that considering sex and hormonal status in basic and clinical studies will bring precision to the development of novel therapeutics and treatment paradigms. To this end, funding, regulatory and policy agencies now require inclusion of female animals and women in basic and clinical studies. However, inclusion of female animals and women often does not mean that information regarding potential hormonal interactions with pharmacological treatments or clinical outcomes is available. All sex steroid hormones can interact with receptors for drug targets, metabolism and transport. Genetic variation in receptors or in enzymatic function might contribute to sex differences in therapeutic efficacy and adverse drug reactions. Outcomes from clinical trials are often not reported by sex, and, if the data are available, they are not translated into clinical practice guidelines. This Review will provide a historical perspective for the current state of research related to hormone trials and provide concrete strategies that, if implemented, will improve the health of all people.

The X factor: A robust and powerful approach to X-chromosome-inclusive whole-genome association studies

The X‐chromosome is often excluded from genome‐wide association studies because of analytical challenges. Some of the problems, such as the random, skewed, or no X‐inactivation model uncertainty, have been investigated. Other considerations have received little to no attention, such as the value in considering nonadditive and gene–sex interaction effects, and the inferential consequence of choosing different baseline alleles (i.e., the reference vs. the alternative allele). Here we propose a unified and flexible regression‐based association test for X‐chromosomal variants. We provide theoretical justifications for its robustness in the presence of various model uncertainties, as well as for its improved power when compared with the existing approaches under certain scenarios. For completeness, we also revisit the autosomes and show that the proposed framework leads to a more robust approach than the standard method. Finally, we provide supporting evidence by revisiting several published association studies. Supporting Information for this article are available online.

A statistical measure for the skewness of X chromosome inactivation for quantitative traits and its application to the MCTFR data

Background

X chromosome inactivation (XCI) is that one of two chromosomes in mammalian females is silenced during early development of embryos. There has been a statistical measure for the degree of the skewness of XCI for qualitative traits. However, no method is available for such task at quantitative trait loci.

Results

In this article, we extend the existing statistical measure for the skewness of XCI for qualitative traits, and the likelihood ratio, Fieller’s and delta methods for constructing the corresponding confidence intervals, and make them accommodate quantitative traits. The proposed measure is a ratio of two linear regression coefficients when association exists. Noting that XCI may cause variance heterogeneity of the traits across different genotypes in females, we obtain the point estimate and confidence intervals of the measure by incorporating such information. The hypothesis testing of the proposed methods is also investigated. We conduct extensive simulation studies to assess the performance of the proposed methods. Simulation results demonstrate that the median of the point estimates of the measure is very close to the pre-specified true value. The likelihood ratio and Fieller’s methods control the size well, and have the similar test power and accurate coverage probability, which perform better than the delta method. So far, we are not aware of any association study for the X-chromosomal loci in the Minnesota Center for Twin and Family Research data. So, we apply our proposed methods to these data for their practical use and find that only the rs792959 locus, which is simultaneously associated with the illicit drug composite score and behavioral disinhibition composite score, may undergo XCI skewing. However, this needs to be confirmed by molecular genetics.

Conclusions

We recommend the Fieller’s method in practical use because it is a non-iterative procedure and has the similar performance to the likelihood ratio method.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00978-z.

A Multiancestry Sex-Stratified Genome-Wide Association Study of Spontaneous Clearance of Hepatitis C Virus

BACKGROUND

Spontaneous clearance of acute hepatitis C virus (HCV) infection is more common in women than in men, independent of known risk factors.

METHODS

To identify sex-specific genetic loci, we studied 4423 HCV-infected individuals (2903 male, 1520 female) of European, African, and Hispanic ancestry. We performed autosomal, and X chromosome sex-stratified and combined association analyses in each ancestry group.

RESULTS

A male-specific region near the adenosine diphosphate-ribosylation factor-like 5B (ARL5B) gene was identified. Individuals with the C allele of rs76398191 were about 30% more likely to have chronic HCV infection than individuals with the T allele (OR, 0.69; P = 1.98 × 10-07), and this was not seen in females. The ARL5B gene encodes an interferon-stimulated gene that inhibits immune response to double-stranded RNA viruses. We also identified suggestive associations near septin 6 and ribosomal protein L39 genes on the X chromosome. In box sexes, allele G of rs12852885 was associated with a 40% increase in HCV clearance compared with the A allele (OR, 1.4; P = 2.46 × 10-06). Septin 6 facilitates HCV replication via interaction with the HCV NS5b protein, and ribosomal protein L39 acts as an HCV core interactor.

CONCLUSIONS

These novel gene associations support differential mechanisms of HCV clearance between the sexes and provide biological targets for treatment or vaccine development.

Testing and estimation of X-chromosome SNP effects: Impact of model assumptions

Interest in analyzing X chromosome single nucleotide polymorphisms (SNPs) is growing and several approaches have been proposed. Prior studies have compared power of different approaches, but bias and interpretation of coefficients have received less attention. We performed simulations to demonstrate the impact of X chromosome model assumptions on effect estimates. We investigated the coefficient biases of SNP and sex effects with commonly used models for X chromosome SNPs, including models with and without assumptions of X chromosome inactivation (XCI), and with and without SNP–sex interaction terms. Sex and SNP coefficient biases were observed when assumptions made about XCI and sex differences in SNP effect in the analysis model were inconsistent with the data‐generating model. However, including a SNP–sex interaction term often eliminated these biases. To illustrate these findings, estimates under different genetic model assumptions are compared and interpreted in a real data example. Models to analyze X chromosome SNPs make assumptions beyond those made in autosomal variant analysis. Assumptions made about X chromosome SNP effects should be stated clearly when reporting and interpreting X chromosome associations. Fitting models with SNP × Sex interaction terms can avoid reliance on assumptions, eliminating coefficient bias even in the absence of sex differences in SNP effect.