The effect of X-linked dosage compensation on complex trait variation

Role of X chromosome and dosage-compensation mechanisms in complex trait genetics

The X chromosome (chrX) is often excluded from genome-wide association studies due to its unique biology complicating the analysis and interpretation of genetic data. Consequently, the influence of chrX on human complex traits remains debated. Here, we systematically assessed the relevance of chrX and the effect of its biology on complex traits by analyzing 48 quantitative traits in 343,695 individuals in UK Biobank with replication in 412,181 individuals from FinnGen. We show that, in the general population, chrX contributes to complex trait heritability at a rate of 3% of the autosomal heritability, consistent with the amount of genetic variation observed in chrX. We find that a pronounced male bias in chrX heritability supports the presence of near-complete dosage compensation between sexes through X chromosome inactivation (XCI). However, we also find subtle yet plausible evidence of escape from XCI contributing to human height. Assuming full XCI, the observed chrX contribution to complex trait heritability in both sexes is greater than expected given the presence of only a single active copy of chrX, mirroring potential dosage compensation between chrX and the autosomes. We find this enhanced contribution attributable to systematically larger active allele effects from chrX compared to autosomes in both sexes, independent of allele frequency and variant deleteriousness. Together, these findings support a model in which the two dosage-compensation mechanisms work in concert to balance the influence of chrX across the population while preserving sex-specific differences at a manageable level. Overall, our study advocates for more comprehensive locus discovery efforts in chrX.

X-Linked Gene Dosage and SOX2 Act as Key Roadblocks for Human Germ Cell Specification in Klinefelter Syndrome

Klinefelter syndrome (KS), characterized by the presence of at least one extra X-chromosome, is a common cause of male infertility. However, the mechanism underlying the failure of germline specification is not well studied. Intriguingly, the differentiation efficiency of female human pluripotent stem cells (hPSCs) is often lower than that of male. This study investigates how X-linked gene dosage affects human primordial germ cell-like cells (hPGCLCs) specification in both healthy and diseased conditions. This work reveals that X-linked genes play a multifaceted role against the fate competency to hPGCLCs, with escape genes IGSF1 and CHRDL1 inhibiting the TGF-beta/Activin A and BMP pathways, respectively. Notably, this work identifies a previously unrecognized role of SOX2, upregulated by the escape gene USP9X, elucidating a species-specific function in the mammalian germline. The USP9X-SOX2 regulatory axis profoundly influenced cellular metabolism, mitochondrial morphology, and progenitor competence in hPGCLCs specification. Furthermore, the inability to downregulate SOX2 and upregulate SOX17 in response to BMP signaling impedes downstream gene activation due to motif binding competition. These findings shed novel insights into the human germline specification by elucidating the divergent roles of SOX2 versus SOX17 in mammals, influenced by X-linked gene dosage effects. These results offer potential applications for improving the induction efficiency of hPGCLCs, facilitating disease mechanistic studies.

Mapping variants in thyroid hormone transporter MCT8 to disease severity by genomic, phenotypic, functional, structural and deep learning integration

Predicting and quantifying phenotypic consequences of genetic variants in rare disorders is a major challenge, particularly pertinent for 'actionable' genes such as thyroid hormone transporter MCT8 (encoded by the X-linked SLC16A2 gene), where loss-of-function (LoF) variants cause a rare neurodevelopmental and (treatable) metabolic disorder in males. The combination of deep phenotyping data with functional and computational tests and with outcomes in population cohorts, enabled us to: (i) identify the genetic aetiology of divergent clinical phenotypes of MCT8 deficiency with genotype-phenotype relationships present across survival and 24 out of 32 disease features; (ii) demonstrate a mild phenocopy in ~400,000 individuals with common genetic variants in MCT8; (iii) assess therapeutic effectiveness, which did not differ among LoF-categories; (iv) advance structural insights in normal and mutated MCT8 by delineating seven critical functional domains; (v) create a pathogenicity-severity MCT8 variant classifier that accurately predicted pathogenicity (AUC:0.91) and severity (AUC:0.86) for 8151 variants. Our information-dense mapping provides a generalizable approach to advance multiple dimensions of rare genetic disorders.

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.

Genetic parameters and genome-wide association studies including the X chromosome for various reproduction and semen quality traits in Nellore cattle

BACKGROUND

The profitability of the beef industry is directly influenced by the fertility rate and reproductive performance of both males and females, which can be improved through selective breeding. When performing genomic analyses, genetic markers located on the X chromosome have been commonly ignored despite the X chromosome being one of the largest chromosomes in the cattle genome. Therefore, the primary objectives of this study were to: (1) estimate variance components and genetic parameters for eighteen male and five female fertility and reproductive traits in Nellore cattle including X chromosome markers in the analyses; and (2) perform genome-wide association studies and functional genomic analyses to better understand the genetic background of male and female fertility and reproductive performance traits in Nellore cattle.

RESULTS

The percentage of the total direct heritability (h2total) explained by the X chromosome markers (h2x) ranged from 3 to 32% (average: 16.4%) and from 9 to 67% (average: 25.61%) for female reproductive performance and male fertility traits, respectively. Among the traits related to breeding soundness evaluation, the overall bull and semen evaluation and semen quality traits accounted for the highest proportion of h2x relative to h2total with an average of 39.5% and 38.75%, respectively. The total number of significant genomic markers per trait ranged from 7 (seminal vesicle width) to 43 (total major defects). The number of significant markers located on the X chromosome ranged from zero to five. A total of 683, 252, 694, 382, 61, and 77 genes overlapped with the genomic regions identified for traits related to female reproductive performance, semen quality, semen morphology, semen defects, overall bulls' fertility evaluation, and overall semen evaluation traits, respectively. The key candidate genes located on the X chromosome are PRR32, STK26, TMSB4X, TLR7, PRPS2, SMS, SMARCA1, UTP14A, and BCORL1. The main gene ontology terms identified are "Oocyte Meiosis", "Progesterone Mediated Oocyte Maturation", "Thermogenesis", "Sperm Flagellum", and "Innate Immune Response".

CONCLUSIONS

Our findings indicate the key role of genes located on the X chromosome on the phenotypic variability of male and female reproduction and fertility traits in Nellore cattle. Breeding programs aiming to improve these traits should consider adding the information from X chromosome markers in their genomic analyses.

Inter-chromosomal contacts demarcate genome topology along a spatial gradient

Non-homologous chromosomal contacts (NHCCs) between different chromosomes participate considerably in gene and genome regulation. Due to analytical challenges, NHCCs are currently considered as singular, stochastic events, and their extent and fundamental principles across cell types remain controversial. We develop a supervised and unsupervised learning algorithm, termed Signature, to call NHCCs in Hi-C datasets to advance our understanding of genome topology. Signature reveals 40,282 NHCCs and their properties across 62 Hi-C datasets of 53 diploid human cell types. Genomic regions of NHCCs are gene-dense, highly expressed, and harbor genes for cell-specific and sex-specific functions. Extensive inter-telomeric and inter-centromeric clustering occurs across cell types [Rabl's configuration] and 61 NHCCs are consistently found at the nuclear speckles. These constitutive 'anchor loci' facilitate an axis of genome activity whilst cell-type-specific NHCCs act in discrete hubs. Our results suggest that non-random chromosome positioning is supported by constitutive NHCCs that shape genome topology along an off-centered spatial gradient of genome activity.

Role of the X Chromosome in Alzheimer Disease Genetics

Importance

The X chromosome has remained enigmatic in Alzheimer disease (AD), yet it makes up 5% of the genome and carries a high proportion of genes expressed in the brain, making it particularly appealing as a potential source of unexplored genetic variation in AD.

Objectives

To perform the first large-scale X chromosome-wide association study (XWAS) of AD.

Design, Setting, and Participants

This was a meta-analysis of genetic association studies in case-control, family-based, population-based, and longitudinal AD-related cohorts from the US Alzheimer's Disease Genetics Consortium, the Alzheimer's Disease Sequencing Project, the UK Biobank, the Finnish health registry, and the US Million Veterans Program. Risk of AD was evaluated through case-control logistic regression analyses. Data were analyzed between January 2023 and March 2024. Genetic data available from high-density single-nucleotide variant microarrays and whole-genome sequencing and summary statistics for multitissue expression and protein quantitative trait loci available from published studies were included, enabling follow-up genetic colocalization analyses. A total of 1 629 863 eligible participants were selected from referred and volunteer samples, 477 596 of whom were excluded for analysis exclusion criteria. The number of participants who declined to participate in original studies was not available.

Main Outcome and Measures

Risk of AD, reported as odds ratios (ORs) with 95% CIs. Associations were considered at X chromosome-wide (P < 1 × 10-5) and genome-wide (P < 5 × 10-8) significance. Primary analyses are nonstratified, while secondary analyses evaluate sex-stratified effects.

Results

Analyses included 1 152 284 participants of non-Hispanic White, European ancestry (664 403 [57.7%] female and 487 881 [42.3%] male), including 138 558 individuals with AD. Six independent genetic loci passed X chromosome-wide significance, with 4 showing support for links between the genetic signal for AD and expression of nearby genes in brain and nonbrain tissues. One of these 4 loci passed conservative genome-wide significance, with its lead variant centered on an intron of SLC9A7 (OR, 1.03; 95% CI, 1.02-1.04) and colocalization analyses prioritizing both the SLC9A7 and nearby CHST7 genes. Of these 6 loci, 4 displayed evidence for escape from X chromosome inactivation with regard to AD risk.

Conclusion and Relevance

This large-scale XWAS of AD identified the novel SLC9A7 locus. SLC9A7 regulates pH homeostasis in Golgi secretory compartments and is anticipated to have downstream effects on amyloid β accumulation. Overall, this study advances our knowledge of AD genetics and may provide novel biological drug targets. The results further provide initial insights into elucidating the role of the X chromosome in sex-based differences in AD.

Genetic architecture of telomere length in 462,666 UK Biobank whole-genome sequences

Telomeres protect chromosome ends from damage and their length is linked with human disease and aging. We developed a joint telomere length metric, combining quantitative PCR and whole-genome sequencing measurements from 462,666 UK Biobank participants. This metric increased SNP heritability, suggesting that it better captures genetic regulation of telomere length. Exome-wide rare-variant and gene-level collapsing association studies identified 64 variants and 30 genes significantly associated with telomere length, including allelic series in ACD and RTEL1. Notably, 16% of these genes are known drivers of clonal hematopoiesis-an age-related somatic mosaicism associated with myeloid cancers and several nonmalignant diseases. Somatic variant analyses revealed gene-specific associations with telomere length, including lengthened telomeres in individuals with large SRSF2-mutant clones, compared with shortened telomeres in individuals with clonal expansions driven by other genes. Collectively, our findings demonstrate the impact of rare variants on telomere length, with larger effects observed among genes also associated with clonal hematopoiesis.

Quantification of escape from X chromosome inactivation with single-cell omics data reveals heterogeneity across cell types and tissues

Several X-linked genes escape from X chromosome inactivation (XCI), while differences in escape across cell types and tissues are still poorly characterized. Here, we developed scLinaX for directly quantifying relative gene expression from the inactivated X chromosome with droplet-based single-cell RNA sequencing (scRNA-seq) data. The scLinaX and differentially expressed gene analyses with large-scale blood scRNA-seq datasets consistently identified the stronger escape in lymphocytes than in myeloid cells. An extension of scLinaX to a 10x multiome dataset (scLinaX-multi) suggested a stronger escape in lymphocytes than in myeloid cells at the chromatin-accessibility level. The scLinaX analysis of human multiple-organ scRNA-seq datasets also identified the relatively strong degree of escape from XCI in lymphoid tissues and lymphocytes. Finally, effect size comparisons of genome-wide association studies between sexes suggested the underlying impact of escape on the genotype-phenotype association. Overall, scLinaX and the quantified escape catalog identified the heterogeneity of escape across cell types and tissues.

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.

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.

The Mechanisms of Long Non-coding RNA-XIST in Ischemic Stroke: Insights into Functional Roles and Therapeutic Potential

Ischemic stroke, which occurs due to the occlusion of cerebral arteries, is a common type of stroke. Recent research has highlighted the important role of long non-coding RNAs (lncRNAs) in the development of cerebrovascular diseases, specifically ischemic stroke. Understanding the functional roles of lncRNAs in ischemic stroke is crucial, given their potential contribution to the disease pathology. One noteworthy lncRNA is X-inactive specific transcript (XIST), which exhibits downregulation during the early stages of ischemic stroke and subsequent upregulation in later stages. XIST exert its influence on the development of ischemic stroke through interactions with multiple miRNAs and transcription factors. These interactions play a significant role in the pathogenesis of the condition. In this review, we have provided a comprehensive summary of the functional roles of XIST in ischemic stroke. By investigating the involvement of XIST in the disease process, we aim to enhance our understanding of the mechanisms underlying ischemic stroke and potentially identify novel therapeutic targets.

Sex-biased gene expression during neural differentiation of human embryonic stem cells

Sex differences in the developing human brain are primarily attributed to hormonal influence. Recently however, genetic differences and their impact on the developing nervous system have attracted increased attention. To understand genetically driven sexual dimorphisms in neurodevelopment, we investigated genome-wide gene expression in an in vitro differentiation model of male and female human embryonic stem cell lines (hESC), independent of the effects of human sex hormones. Four male and four female-derived hESC lines were differentiated into a population of mixed neurons over 37 days. Differential gene expression and gene set enrichment analyses were conducted on bulk RNA sequencing data. While similar differentiation tendencies in all cell lines demonstrated the robustness and reproducibility of our differentiation protocol, we found sex-biased gene expression already in undifferentiated ESCs at day 0, but most profoundly after 37 days of differentiation. Male and female cell lines exhibited sex-biased expression of genes involved in neurodevelopment, suggesting that sex influences the differentiation trajectory. Interestingly, the highest contribution to sex differences was found to arise from the male transcriptome, involving both Y chromosome and autosomal genes. We propose 13 sex-biased candidate genes (10 upregulated in male cell lines and 3 in female lines) that are likely to affect neuronal development. Additionally, we confirmed gene dosage compensation of X/Y homologs escaping X chromosome inactivation through their Y homologs and identified a significant overexpression of the Y-linked demethylase UTY and KDM5D in male hESC during neuron development, confirming previous results in neural stem cells. Our results suggest that genetic sex differences affect neuronal differentiation trajectories, which could ultimately contribute to sex biases during human brain development.

The Role of X Chromosome in Alzheimer's Disease Genetics

Importance

The X chromosome has remained enigmatic in Alzheimer's disease (AD), yet it makes up 5% of the genome and carries a high proportion of genes expressed in the brain, making it particularly appealing as a potential source of unexplored genetic variation in AD.

Objectives

Perform the first large-scale X chromosome-wide association study (XWAS) of AD. Primary analyses are non-stratified, while secondary analyses evaluate sex-stratified effects.

Design

Meta-analysis of genetic association studies in case-control, family-based, population-based, and longitudinal AD-related cohorts from the US Alzheimer's Disease Genetics Consortium (ADGC) and Alzheimer's Disease Sequencing Project (ADSP), the UK Biobank (UKB), the Finnish health registry (FinnGen), and the US Million Veterans Program (MVP). Risk for AD evaluated through case-control logistic regression analyses. Data were analyzed between January 2023 and March 2024.

Setting

Genetic data available from high-density single-nucleotide polymorphism (SNP) microarrays and whole-genome sequencing (WGS). Summary statistics for multi-tissue expression and protein quantitative trait loci (QTL) available from published studies, enabling follow-up genetic colocalization analyses.

Participants

1,629,863 eligible participants were selected from referred and volunteer samples, of which 477,596 were excluded for analysis exclusion criteria. Number of participants who declined to participate in original studies was not available.

Main Outcome and Measures

Risk for AD (odds ratio; OR) with 95% confidence intervals (CI). Associations were considered at X-chromosome-wide (P-value<1e-5) and genome-wide (P-value<5e-8) significance.

Results

Analyses included 1,152,284 non-Hispanic White European ancestry subjects (57.3% females), including 138,558 cases. 6 independent genetic loci passed X-chromosome-wide significance, with 4 showing support for causal links between the genetic signal for AD and expression of nearby genes in brain and non-brain tissues. One of these 4 loci passed conservative genome-wide significance, with its lead variant centered on an intron of SLC9A7 (OR=1.054, 95%-CI=[1.035, 1.075]) and colocalization analyses prioritizing both the SLC9A7 and nearby CHST7 genes.

Conclusion and Relevance

We performed the first large-scale XWAS of AD and identified the novel SLC9A7 locus. SLC9A7 regulates pH homeostasis in Golgi secretory compartments and is anticipated to have downstream effects on amyloid beta accumulation. Overall, this study significantly advances our knowledge of AD genetics and may provide novel biological drug targets.

Better together against genetic heterogeneity: A sex-combined joint main and interaction analysis of 290 quantitative traits in the UK Biobank

Genetic effects can be sex-specific, particularly for traits such as testosterone, a sex hormone. While sex-stratified analysis provides easily interpretable sex-specific effect size estimates, the presence of sex-differences in SNP effect implies a SNP×sex interaction. This suggests the usage of the often overlooked joint test, testing for an SNP's main and SNP×sex interaction effects simultaneously. Notably, even without individual-level data, the joint test statistic can be derived from sex-stratified summary statistics through an omnibus meta-analysis. Utilizing the available sex-stratified summary statistics of the UK Biobank, we performed such omnibus meta-analyses for 290 quantitative traits. Results revealed that this approach is robust to genetic effect heterogeneity and can outperform the traditional sex-stratified or sex-combined main effect-only tests. Therefore, we advocate using the omnibus meta-analysis that captures both the main and interaction effects. Subsequent sex-stratified analysis should be conducted for sex-specific effect size estimation and interpretation.

Polygenic response of sex chromosomes to sexual antagonism

Sexual antagonism occurs when males and females differ in their phenotypic fitness optima but are constrained in their evolution to these optima because of their shared genome. The sex chromosomes, which have distinct evolutionary "interests" relative to the autosomes, are theorized to play an important role in sexually antagonistic conflict. However, the evolutionary responses of sex chromosomes and autosomes have usually been considered independently, that is, via contrasting the response of a gene located on either an X chromosome or an autosome. Here, we study the coevolutionary response of the X chromosome and autosomes to sexually antagonistic selection acting on a polygenic phenotype. We model a phenotype initially under stabilizing selection around a single optimum, followed by a sudden divergence of the male and female optima. We find that, in the absence of dosage compensation, the X chromosome promotes evolution toward the female optimum, inducing coevolutionary male-biased responses on the autosomes. Dosage compensation obscures the female-biased interests of the X, causing it to contribute equally to male and female phenotypic change. We further demonstrate that fluctuations in an adaptive landscape can generate prolonged intragenomic conflict and accentuate the differential responses of the X and autosomes to this conflict.

The correlates of neonatal complement component 3 and 4 protein concentrations with a focus on psychiatric and autoimmune disorders

Complement components have been linked to schizophrenia and autoimmune disorders. We examined the association between neonatal circulating C3 and C4 protein concentrations in 68,768 neonates and the risk of six mental disorders. We completed genome-wide association studies (GWASs) for C3 and C4 and applied the summary statistics in Mendelian randomization and phenome-wide association studies related to mental and autoimmune disorders. The GWASs for C3 and C4 protein concentrations identified 15 and 36 independent loci, respectively. We found no associations between neonatal C3 and C4 concentrations and mental disorders in the total sample (both sexes combined); however, post-hoc analyses found that a higher C3 concentration was associated with a reduced risk of schizophrenia in females. Mendelian randomization based on C4 summary statistics found an altered risk of five types of autoimmune disorders. Our study adds to our understanding of the associations between C3 and C4 concentrations and subsequent mental and autoimmune disorders.

Unravelling the genetic variability of host resilience to endo- and ectoparasites in Nellore commercial herds

BACKGROUND

Host resilience (HR) to parasites can affect the performance of animals. Therefore, the aim of this study was to present a detailed investigation of the genetic mechanisms of HR to ticks (TICK), gastrointestinal nematodes (GIN), and Eimeria spp. (EIM) in Nellore cattle that were raised under natural infestation and a prophylactic parasite control strategy. In our study, HR was defined as the slope coefficient of body weight (BW) when TICK, GIN, and EIM burdens were used as environmental gradients in random regression models. In total, 1712 animals were evaluated at five measurement events (ME) at an average age of 331, 385, 443, 498, and 555 days, which generated 7307 body weight (BW) records. Of the 1712 animals, 1075 genotyped animals were used in genome-wide association studies to identify genomic regions associated with HR.

RESULTS

Posterior means of the heritability estimates for BW ranged from 0.09 to 0.54 across parasites and ME. The single nucleotide polymorphism (SNP)-derived heritability for BW at each ME ranged from a low (0.09 at ME.331) to a moderate value (0.23 at ME.555). Those estimates show that genetic progress can be achieved for BW through selection. Both genetic and genomic associations between BW and HR to TICK, GIN, and EIM confirmed that parasite infestation impacted the performance of animals. Selection for BW under an environment with a controlled parasite burden is an alternative to improve both, BW and HR. There was no impact of age of measurement on the estimates of genetic variance for HR. Five quantitative trait loci (QTL) were associated with HR to EIM but none with HR to TICK and to GIN. These QTL contain genes that were previously shown to be associated with the production of antibody modulators and chemokines that are released in the intestinal epithelium.

CONCLUSIONS

Selection for BW under natural infestation and controlled parasite burden, via prophylactic parasite control, contributes to the identification of animals that are resilient to nematodes and Eimeria ssp. Although we verified that sufficient genetic variation existed for HR, we did not find any genes associated with mechanisms that could justify the expression of HR to TICK and GIN.

Molecular bases of comorbidities: present and future perspectives

Co-occurrence of diseases decreases patient quality of life, complicates treatment choices, and increases mortality. Analyses of electronic health records present a complex scenario of comorbidity relationships that vary by age, sex, and cohort under study. The study of similarities between diseases using 'omics data, such as genes altered in diseases, gene expression, proteome, and microbiome, are fundamental to uncovering the origin of, and potential treatment for, comorbidities. Recent studies have produced a first generation of genetic interpretations for as much as 46% of the comorbidities described in large cohorts. Integrating different sources of molecular information and using artificial intelligence (AI) methods are promising approaches for the study of comorbidities. They may help to improve the treatment of comorbidities, including the potential repositioning of drugs.

Boosting the power of genome-wide association studies within and across ancestries by using polygenic scores

Genome-wide association studies (GWASs) have been mostly conducted in populations of European ancestry, which currently limits the transferability of their findings to other populations. Here, we show, through theory, simulations and applications to real data, that adjustment of GWAS analyses for polygenic scores (PGSs) increases the statistical power for discovery across all ancestries. We applied this method to analyze seven traits available in three large biobanks with participants of East Asian ancestry (n = 340,000 in total) and report 139 additional associations across traits. We also present a two-stage meta-analysis strategy whereby, in contributing cohorts, a PGS-adjusted GWAS is rerun using PGSs derived from a first round of a standard meta-analysis. On average, across traits, this approach yields a 1.26-fold increase in the number of detected associations (range 1.07- to 1.76-fold increase). Altogether, our study demonstrates the value of using PGSs to increase the power of GWASs in underrepresented populations and promotes such an analytical strategy for future GWAS meta-analyses.

Chromosomal inversion polymorphisms shape human brain morphology

The impact of chromosomal inversions on human brain morphology remains underexplored. We studied 35 common inversions classified from genotypes of 33,018 adults with European ancestry. The inversions at 2p22.3, 16p11.2, and 17q21.31 reach genome-wide significance, followed by 8p23.1 and 6p21.33, in their association with cortical and subcortical morphology. The 17q21.31, 8p23.1, and 16p11.2 regions comprise the LRRC37, OR7E, and NPIP duplicated gene families. We find the 17q21.31 MAPT inversion region, known for harboring neurological risk, to be the most salient locus among common variants for shaping and patterning the cortex. Overall, we observe the inverted orientations decreasing brain size, with the exception that the 2p22.3 inversion is associated with increased subcortical volume and the 8p23.1 inversion is associated with increased motor cortex. These significant inversions are in the genomic hotspots of neuropsychiatric loci. Our findings are generalizable to 3,472 children and demonstrate inversions as essential genetic variation to understand human brain phenotypes.

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.

Sex and gender in neurodevelopmental conditions

Health-related conditions often differ qualitatively or quantitatively between individuals of different birth-assigned sexes and gender identities, and/or with different gendered experiences, requiring tailored care. Studying the moderating and mediating effects of sex-related and gender-related factors on impairment, disability, wellbeing and health is of paramount importance especially for neurodivergent individuals, who are diagnosed with neurodevelopmental conditions with uneven sex/gender distributions. Researchers have become aware of the myriad influences that sex-related and gender-related variables have on the manifestations of neurodevelopmental conditions, and contemporary work has begun to investigate the mechanisms through which these effects are mediated. Here we describe topical concepts of sex and gender science, summarize current knowledge, and discuss research and clinical challenges related to autism, attention-deficit/hyperactivity disorder and other neurodevelopmental conditions. We consider sex and gender in the context of epidemiology, behavioural phenotypes, neurobiology, genetics, endocrinology and neighbouring disciplines. The available evidence supports the view that sex and gender are important contributors to the biological and behavioural variability in neurodevelopmental conditions. Methodological caveats such as frequent conflation of sex and gender constructs, inappropriate measurement of these constructs and under-representation of specific demographic groups (for example, female and gender minority individuals and people with intellectual disabilities) limit the translational potential of research so far. Future research and clinical implementation should integrate sex and gender into next-generation diagnostics, mechanistic investigations and support practices.

Physiological and molecular sex differences in human skeletal muscle in response to exercise training

Sex differences in exercise physiology, such as substrate metabolism and skeletal muscle fatigability, stem from inherent biological factors, including endogenous hormones and genetics. Studies investigating exercise physiology frequently include only males or do not take sex differences into consideration. Although there is still an underrepresentation of female participants in exercise research, existing studies have identified sex differences in physiological and molecular responses to exercise training. The observed sex differences in exercise physiology are underpinned by the sex chromosome complement, sex hormones and, on a molecular level, the epigenome and transcriptome. Future research in the field should aim to include both sexes, control for menstrual cycle factors, conduct large-scale and ethnically diverse studies, conduct meta-analyses to consolidate findings from various studies, leverage unique cohorts (such as post-menopausal, transgender, and those with sex chromosome abnormalities), as well as integrate tissue and cell-specific -omics data. This knowledge is essential for developing deeper insight into sex-specific physiological responses to exercise training, thus directing future exercise physiology studies and practical application.

Long noncoding RNA XIST: Mechanisms for X chromosome inactivation, roles in sex-biased diseases, and therapeutic opportunities

Sexual dimorphism has been reported in various human diseases including autoimmune diseases, neurological diseases, pulmonary arterial hypertension, and some types of cancers, although the underlying mechanisms remain poorly understood. The long noncoding RNA (lncRNA) X-inactive specific transcript (XIST) is involved in X chromosome inactivation (XCI) in female placental mammals, a process that ensures the balanced expression dosage of X-linked genes between sexes. XIST is abnormally expressed in many sex-biased diseases. In addition, escape from XIST-mediated XCI and skewed XCI also contribute to sex-biased diseases. Therefore, its expression or modification can be regarded as a biomarker for the diagnosis and prognosis of many sex-biased diseases. Genetic manipulation of XIST expression can inhibit the progression of some of these diseases in animal models, and therefore XIST has been proposed as a potential therapeutic target. In this manuscript, we summarize the current knowledge about the mechanisms for XIST-mediated XCI and the roles of XIST in sex-biased diseases, and discuss potential therapeutic strategies targeting XIST.

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.

Inheritance of hormonal stress response and temperament in infant rhesus macaques (Macaca Mulatta): Nonadditive and sex-specific effects

OBJECTIVE

Early life interindividual variation in hypothalamic-pituitary-adrenal (HPA) reactivity to stress is predictive of later life psychological and physical well-being, including the development of many pathological syndromes that are often sex-biased. A complex and interactive set of environmental and genetic causes for such variation has been implicated by previous studies, though little attention has been paid to nonadditive effects (e.g. dominance, X-linked) or sex-specific genetic effects.

METHOD

We used a large pedigreed sample of captive 3-4 months old infant rhesus macaques (N = 2,661, 54% female) to fit univariate and multivariate linear mixed quantitative genetic models for four longitudinal blood cortisol samples and three reliable ratings of infant temperament (nervousness, gentleness, confidence) during a mother-infant separation protocol.

RESULTS

Each trait had a moderate narrow-sense heritability (h², 0.26-0.46), but dominance effects caused the first two cortisol samples to have much larger broad-sense heritabilities (H², 0.57 and 0.77). We found no evidence for X-linked variance or common maternal environment variance. There was a sex difference in heritability of the first cortisol sample (h_f² < h_m²), suggesting differing genetic architecture of perception of maternal separation and relocation during infancy. Otherwise, genetic covariance matrices for the sexes were very similar. Genetic correlations between cortisol levels and temperament were weak (< |0.4|) but stronger than residual or phenotypic correlations.

CONCLUSIONS

HPA reactivity and temperament had a primarily additive genetic basis in infant macaques, but there were important complexities to the genetic architecture of including genetic dominance and sex differences in heritability at this early life stage. (PsycInfo Database Record (c) 2022 APA, all rights reserved).

Sex differences in deleterious mutational effects in Drosophila melanogaster: combining quantitative and population genetic insights

Fitness effects of deleterious mutations can differ between females and males due to: (i) sex differences in the strength of purifying selection; and (ii) sex differences in ploidy. Although sex differences in fitness effects have important broader implications (e.g., for the evolution of sex and lifespan), few studies have quantified their scope. Those that have belong to one of two distinct empirical traditions: (i) quantitative genetics, which focusses on multi-locus genetic variances in each sex, but is largely agnostic about their genetic basis; and (ii) molecular population genetics, which focusses on comparing autosomal and X-linked polymorphism, but is poorly suited for inferring contemporary sex differences. Here, we combine both traditions to present a comprehensive analysis of female and male adult reproductive fitness among 202 outbred, laboratory-adapted, hemiclonal genomes of Drosophila melanogaster. While we find no clear evidence for sex differences in the strength of purifying selection, sex differences in ploidy generate multiple signals of enhanced purifying selection for X-linked loci. These signals are present in quantitative genetic metrics-i.e., a disproportionate contribution of the X to male (but not female) fitness variation-and population genetic metrics-i.e., steeper regressions of an allele's average fitness effect on its frequency, and proportionally less nonsynonymous polymorphism on the X than autosomes. Fitting our data to models for both sets of metrics, we infer that deleterious alleles are partially recessive. Given the often-large gap between quantitative and population genetic estimates of evolutionary parameters, our study showcases the benefits of combining genomic and fitness data when estimating such parameters.

Skeletal muscle methylome and transcriptome integration reveals profound sex differences related to muscle function and substrate metabolism

Nearly all human complex traits and diseases exhibit some degree of sex differences, with epigenetics being one of the main contributing factors. Various tissues display sex differences in DNA methylation; however, this has not yet been explored in skeletal muscle, despite skeletal muscle being among the tissues with the most transcriptomic sex differences. For the first time, we investigated the effect of sex on autosomal DNA methylation in human skeletal muscle across three independent cohorts (Gene SMART, FUSION, and GSE38291) using a meta-analysis approach, totalling 369 human muscle samples (222 males and 147 females), and integrated this with known sex-biased transcriptomics. We found 10,240 differentially methylated regions (DMRs) at FDR < 0.005, 94% of which were hypomethylated in males, and gene set enrichment analysis revealed that differentially methylated genes were involved in muscle contraction and substrate metabolism. We then investigated biological factors underlying DNA methylation sex differences and found that circulating hormones were not associated with differential methylation at sex-biased DNA methylation loci; however, these sex-specific loci were enriched for binding sites of hormone-related transcription factors (with top TFs including androgen (AR), estrogen (ESR1), and glucocorticoid (NR3C1) receptors). Fibre type proportions were associated with differential methylation across the genome, as well as across 16% of sex-biased DNA methylation loci (FDR < 0.005). Integration of DNA methylomic results with transcriptomic data from the GTEx database and the FUSION cohort revealed 326 autosomal genes that display sex differences at both the epigenome and transcriptome levels. Importantly, transcriptional sex-biased genes were overrepresented among epigenetic sex-biased genes (p value = 4.6e−13), suggesting differential DNA methylation and gene expression between male and female muscle are functionally linked. Finally, we validated expression of three genes with large effect sizes (FOXO3A, ALDH1A1, and GGT7) in the Gene SMART cohort with qPCR. GGT7, involved in antioxidant metabolism, displays male-biased expression as well as lower methylation in males across the three cohorts. In conclusion, we uncovered 8420 genes that exhibit DNA methylation differences between males and females in human skeletal muscle that may modulate mechanisms controlling muscle metabolism and health.

Sex-Specific Association of the X Chromosome With Cognitive Change and Tau Pathology in Aging and Alzheimer Disease

Importance

The X chromosome represents 5% of the human genome in women and men, and its influence on cognitive aging and Alzheimer disease (AD) is largely unknown.

Objective

To determine whether the X chromosome is associated with sex-specific cognitive change and tau pathology in aging and AD.

Design, Setting, Participants

This study examined differential gene expression profiling of the X chromosome from an RNA sequencing data set of the dorsolateral prefrontal cortex obtained from autopsied, elderly individuals enrolled in the Religious Orders Study and Rush Memory and Aging Project joint cohorts. Samples were collected from the cohort study with enrollment from 1994 to 2017. Data were last analyzed in May 2021.

Main Outcomes and Measures

The main analysis examined whether X chromosome gene expression measured by RNA sequencing of the dorsolateral prefrontal cortex was associated with cognitive change during aging and AD, independent of AD pathology and at the transcriptome-wide level in women and men. Whether X chromosome gene expression was associated with neurofibrillary tangle burden, a measure of tau pathology that influences cognition, in women and men was also explored.

Results

Samples for RNA sequencing of the dorsolateral prefrontal cortex were obtained from 508 individuals (mean [SD] age at death, 88.4 [6.6] years; 315 [62.0%] were female; 197 [38.8%] had clinical diagnosis of AD at death; 293 [58.2%] had pathological diagnosis of AD at death) enrolled in the Religious Orders Study and Rush Memory and Aging Project joint cohorts and were followed up annually for a mean (SD) of 6.3 (3.9) years. X chromosome gene expression (29 genes), adjusted for age at death, education, and AD pathology, was significantly associated with cognitive change at the genome-wide level in women but not men. In the majority of identified X genes (19 genes), increased expression was associated with slower cognitive decline in women. In contrast with cognition, X chromosome gene expression (3 genes), adjusted for age at death and education, was associated with neuropathological tau burden at the genome-wide level in men but not women.

Conclusions and Relevance

In this study, the X chromosome was associated with cognitive trajectories and neuropathological tau burden in aging and AD in a sex-specific manner. This is important because specific X chromosome factors could contribute risk or resilience to biological pathways of aging and AD in women, men, or both.

Inferring genes that escape X-Chromosome inactivation reveals important contribution of variable escape genes to sex-biased diseases

The X Chromosome plays an important role in human development and disease. However, functional genomic and disease association studies of X genes greatly lag behind autosomal gene studies, in part owing to the unique biology of X-Chromosome inactivation (XCI). Because of XCI, most genes are only expressed from one allele. Yet, ∼30% of X genes “escape” XCI and are transcribed from both alleles, many only in a proportion of the population. Such interindividual differences are likely to be disease relevant, particularly for sex-biased disorders. To understand the functional biology for X-linked genes, we developed X-Chromosome inactivation for RNA-seq (XCIR), a novel approach to identify escape genes using bulk RNA-seq data. Our method, available as an R package, is more powerful than alternative approaches and is computationally efficient to handle large population-scale data sets. Using annotated XCI states, we examined the contribution of X-linked genes to the disease heritability in the United Kingdom Biobank data set. We show that escape and variable escape genes explain the largest proportion of X heritability, which is in large part attributable to X genes with Y homology. Finally, we investigated the role of each XCI state in sex-biased diseases and found that although XY homologous gene pairs have a larger overall effect size, enrichment for variable escape genes is significantly increased in female-biased diseases. Our results, for the first time, quantitate the importance of variable escape genes for the etiology of sex-biased disease, and our pipeline allows analysis of larger data sets for a broad range of phenotypes.

The genetic architecture of primary biliary cholangitis

Primary biliary cholangitis (PBC) is a rare autoimmune disease of the liver affecting the small bile ducts. From a genetic point of view, PBC is a complex trait and several genetic and environmental factors have been called in action to explain its etiopathogenesis. Similarly to other complex traits, PBC has benefited from the introduction of genome-wide association studies (GWAS), which identified many variants predisposing or protecting toward the development of the disease. While a progressive endeavour toward the characterization of candidate loci and downstream pathways is currently ongoing, there is still a relatively large portion of heritability of PBC to be revealed. In addition, genetic variation behind progression of the disease and therapeutic response are mostly to be investigated yet. This review outlines the state-of-the-art regarding the genetic architecture of PBC and provides some hints for future investigations, focusing on the study of gene-gene interactions, the application of whole-genome sequencing techniques, and the investigation of X chromosome that can be helpful to cover the missing heritability gap in PBC.

Effect of genomic X-chromosome regions on Nelore bull fertility

Scrotal circumference (SC) is a commonly used trait related to sexual precocity in bulls. Genome-wide association studies have uncovered a lot of genes related to this trait, however, only those present on autosomes. The inclusion of the second biggest chromosome (BTAX) can improve the knowledge of the genetic architecture of this trait. In this study, we performed a weighted, single-step, genome-wide association study using a 777 k BovineHD BeadChip (IllumHD) to analyze the association between SNPs and SC in Brazilian Nelore cattle. Phenotypes from 79,300 males and 3263 genotypes (2017 from females and 1246 from males)-(39,367 SNPs markers located at ChrX) were used. We identified eight regions on chromosome X that displayed important associations with SC. The results showed that together the genomic windows explained 28.52% of the genetic variance for the examined trait. Genes with potential functions in reproduction and fertility regulation were highlighted as candidates for sexual precocity rates in Nelore cattle (AFF2 and PJA1). Moreover, we found 10 genes that had not previously been identified as being associated with sexual precocity traits in cattle. These findings will further advance our understanding of the genetic architecture, considering mainly the presence of the chromosome X, for indicine cattle reproductive traits, being useful in the context of genomic prediction in beef cattle.

Genotype-by-environment interactions for reproduction, body composition, and growth traits in maternal-line pigs based on single-step genomic reaction norms

BACKGROUND

There is an increasing need to account for genotype-by-environment (G × E) interactions in livestock breeding programs to improve productivity and animal welfare across environmental and management conditions. This is even more relevant for pigs because selection occurs in high-health nucleus farms, while commercial pigs are raised in more challenging environments. In this study, we used single-step homoscedastic and heteroscedastic genomic reaction norm models (RNM) to evaluate G × E interactions in Large White pigs, including 8686 genotyped animals, for reproduction (total number of piglets born, TNB; total number of piglets born alive, NBA; total number of piglets weaned, NW), growth (weaning weight, WW; off-test weight, OW), and body composition (ultrasound muscle depth, MD; ultrasound backfat thickness, BF) traits. Genetic parameter estimation and single-step genome-wide association studies (ssGWAS) were performed for each trait.

RESULTS

The average performance of contemporary groups (CG) was estimated and used as environmental gradient in the reaction norm analyses. We found that the need to consider heterogeneous residual variance in RNM models was trait dependent. Based on estimates of variance components of the RNM slope and of genetic correlations across environmental gradients, G × E interactions clearly existed for TNB and NBA, existed for WW but were of smaller magnitude, and were not detected for NW, OW, MD, and BF. Based on estimates of the genetic variance explained by the markers in sliding genomic windows in ssGWAS, several genomic regions were associated with the RNM slope for TNB, NBA, and WW, indicating specific biological mechanisms underlying environmental sensitivity, and dozens of novel candidate genes were identified. Our results also provided strong evidence that the X chromosome contributed to the intercept and slope of RNM for litter size traits in pigs.

CONCLUSIONS

We provide a comprehensive description of G × E interactions in Large White pigs for economically-relevant traits and identified important genomic regions and candidate genes associated with GxE interactions on several autosomes and the X chromosome. Implementation of these findings will contribute to more accurate genomic estimates of breeding values by considering G × E interactions, in order to genetically improve the environmental robustness of maternal-line pigs.

Biological Function of Long Non-coding RNA (LncRNA) Xist

Long non-coding RNAs (lncRNAs) regulate gene expression in a variety of ways at epigenetic, chromatin remodeling, transcriptional, and translational levels. Accumulating evidence suggests that lncRNA X-inactive specific transcript (lncRNA Xist) serves as an important regulator of cell growth and development. Despites its original roles in X-chromosome dosage compensation, lncRNA Xist also participates in the development of tumor and other human diseases by functioning as a competing endogenous RNA (ceRNA). In this review, we comprehensively summarized recent progress in understanding the cellular functions of lncRNA Xist in mammalian cells and discussed current knowledge regarding the ceRNA network of lncRNA Xist in various diseases. Long non-coding RNAs (lncRNAs) are transcripts that are more than 200 nt in length and without an apparent protein-coding capacity (Furlan and Rougeulle, 2016; Maduro et al., 2016). These RNAs are believed to be transcribed by the approximately 98-99% non-coding regions of the human genome (Derrien et al., 2012; Fu, 2014; Montalbano et al., 2017; Slack and Chinnaiyan, 2019), as well as a large variety of genomic regions, such as exonic, tronic, and intergenic regions. Hence, lncRNAs are also divided into eight categories: Intergenic lncRNAs, Intronic lncRNAs, Enhancer lncRNAs, Promoter lncRNAs, Natural antisense/sense lncRNAs, Small nucleolar RNA-ended lncRNAs (sno-lncRNAs), Bidirectional lncRNAs, and non-poly(A) lncRNAs (Ma et al., 2013; Devaux et al., 2015; St Laurent et al., 2015; Chen, 2016; Quinn and Chang, 2016; Richard and Eichhorn, 2018; Connerty et al., 2020). A range of evidence has suggested that lncRNAs function as key regulators in crucial cellular functions, including proliferation, differentiation, apoptosis, migration, and invasion, by regulating the expression level of target genes via epigenomic, transcriptional, or post-transcriptional approaches (Cao et al., 2018). Moreover, lncRNAs detected in body fluids were also believed to serve as potential biomarkers for the diagnosis, prognosis, and monitoring of disease progression, and act as novel and potential drug targets for therapeutic exploitation in human disease (Jiang W. et al., 2018; Zhou et al., 2019a). Long non-coding RNA X-inactive specific transcript (lncRNA Xist) are a set of 15,000-20,000 nt sequences localized in the X chromosome inactivation center (XIC) of chromosome Xq13.2 (Brown et al., 1992; Debrand et al., 1998; Kay, 1998; Lee et al., 2013; da Rocha and Heard, 2017; Yang Z. et al., 2018; Brockdorff, 2019). Previous studies have indicated that lncRNA Xist regulate X chromosome inactivation (XCI), resulting in the inheritable silencing of one of the X-chromosomes during female cell development. Also, it serves a vital regulatory function in the whole spectrum of human disease (notably cancer) and can be used as a novel diagnostic and prognostic biomarker and as a potential therapeutic target for human disease in the clinic (Liu et al., 2018b; Deng et al., 2019; Dinescu et al., 2019; Mutzel and Schulz, 2020; Patrat et al., 2020; Wang et al., 2020a). In particular, lncRNA Xist have been demonstrated to be involved in the development of multiple types of tumors including brain tumor, Leukemia, lung cancer, breast cancer, and liver cancer, with the prominent examples outlined in Table 1. It was also believed that lncRNA Xist (Chaligne and Heard, 2014; Yang Z. et al., 2018) contributed to other diseases, such as pulmonary fibrosis, inflammation, neuropathic pain, cardiomyocyte hypertrophy, and osteoarthritis chondrocytes, and more specific details can be found in Table 2. This review summarizes the current knowledge on the regulatory mechanisms of lncRNA Xist on both chromosome dosage compensation and pathogenesis (especially cancer) processes, with a focus on the regulatory network of lncRNA Xist in human disease.

Patterns of recent natural selection on genetic loci associated with sexually differentiated human body size and shape phenotypes

Levels of sex differences for human body size and shape phenotypes are hypothesized to have adaptively reduced following the agricultural transition as part of an evolutionary response to relatively more equal divisions of labor and new technology adoption. In this study, we tested this hypothesis by studying genetic variants associated with five sexually differentiated human phenotypes: height, body mass, hip circumference, body fat percentage, and waist circumference. We first analyzed genome-wide association (GWAS) results for UK Biobank individuals (~194,000 females and ~167,000 males) to identify a total of 114,199 single nucleotide polymorphisms (SNPs) significantly associated with at least one of the studied phenotypes in females, males, or both sexes (P<5x10^-8). From these loci we then identified 3,016 SNPs (2.6%) with significant differences in the strength of association between the female- and male-specific GWAS results at a low false-discovery rate (FDR<0.001). Genes with known roles in sexual differentiation are significantly enriched for co-localization with one or more of these SNPs versus SNPs associated with the phenotypes generally but not with sex differences (2.73-fold enrichment; permutation test; P = 0.0041). We also confirmed that the identified variants are disproportionately associated with greater phenotype effect sizes in the sex with the stronger association value. We then used the singleton density score statistic, which quantifies recent (within the last ~3,000 years; post-agriculture adoption in Britain) changes in the frequencies of alleles underlying polygenic traits, to identify a signature of recent positive selection on alleles associated with greater body fat percentage in females (permutation test; P = 0.0038; FDR = 0.0380), directionally opposite to that predicted by the sex differences reduction hypothesis. Otherwise, we found no evidence of positive selection for sex difference-associated alleles for any other trait. Overall, our results challenge the longstanding hypothesis that sex differences adaptively decreased following subsistence transitions from hunting and gathering to agriculture.

There is uncertainty regarding the evolutionary history of human sex differences for quantitative body size and shape phenotypes. In this study we identified thousands of genetic loci that differentially impact body size and shape trait variation between females and males using a large sample of UK Biobank individuals. After confirming the biological plausibility of these loci, we used a population genomics approach to study the recent (over the past ~3,000 years) evolutionary histories of these loci in this population. We observed significant increases in the frequencies of alleles associated with greater body fat percentage in females. This result is contradictory to longstanding hypotheses that sex differences have adaptively decreased following subsistence transitions from hunting and gathering to agriculture.

Effect of the X chromosome in genomic evaluations of reproductive traits in beef cattle

The aim of this study was to evaluate whether there are predictive advantages for breeding values with inclusion of X chromosome genomic markers for reproductive (occurrence of early pregnancy - P16 and age at first calving - AFC) and andrological (scrotal circumference -SC) variables in beef cattle. There were 3263 genotypes of females and males evaluated. There were breeding value estimates for SC, AFC and P16 considering two scenarios: 1) only autosomal markers or 2) autosomal and X chromosome markers. To evaluate effects of inclusion of X chromosome markers on selection, responses to selection were compared including or not including genomic marker information from the X chromosome. There were greater heritability estimates for SC (0.40 and 0.31), AFC (0.11 and 0.09) and P16 (0.43 and 0.38) when analyses included, compared with not including, genomic marker information from the X chromosome. When selection is based on results from analyses that did not include information for the X chromosome, there was about a 7 % lesser mean genomic breeding value for the SC traits for selected animals. For P16, there was an approximate 4% lesser breeding value without inclusion of genomic marker information from the X chromosome, while this inclusion did not have as great an effect on the breeding value for AFC. There was an average predictive correlation of 0.79, 0.98 and 0.84 for SC, AFC and P16, respectively. These estimates indicate inclusion of the X chromosome genomic marker information in the analysis can improve prediction of genomic breeding values, especially for SC.

Genetic factors contributing to extensive variability of sex-specific hepatic gene expression in Diversity Outbred mice

Sex-specific transcription characterizes hundreds of genes in mouse liver, many implicated in sex-differential drug and lipid metabolism and disease susceptibility. While the regulation of liver sex differences by growth hormone-activated STAT5 is well established, little is known about autosomal genetic factors regulating the sex-specific liver transcriptome. Here we show, using genotyping and expression data from a large population of Diversity Outbred mice, that genetic factors work in tandem with growth hormone to control the individual variability of hundreds of sex-biased genes, including many long non-coding RNA genes. Significant associations between single nucleotide polymorphisms and sex-specific gene expression were identified as expression quantitative trait loci (eQTLs), many of which showed strong sex-dependent associations. Remarkably, autosomal genetic modifiers of sex-specific genes were found to account for more than 200 instances of gain or loss of sex-specificity across eight Diversity Outbred mouse founder strains. Sex-biased STAT5 binding sites and open chromatin regions with strain-specific variants were significantly enriched at eQTL regions regulating correspondingly sex-specific genes, supporting the proposed functional regulatory nature of the eQTL regions identified. Binding of the male-biased, growth hormone-regulated repressor BCL6 was most highly enriched at trans-eQTL regions controlling female-specific genes. Co-regulated gene clusters defined by overlapping eQTLs included sets of highly correlated genes from different chromosomes, further supporting trans-eQTL action. These findings elucidate how an unexpectedly large number of autosomal factors work in tandem with growth hormone signaling pathways to regulate the individual variability associated with sex differences in liver metabolism and disease.

Mapping genomic regions for reproductive traits in beef cattle: Inclusion of the X chromosome

Although the second largest chromosome of the genome, the X chromosome is usually excluded from genome-wide association studies (GWAS). Considering the presence and importance of genes on this chromosome that are involved in reproduction, the aim of this study was to evaluate the effect of its inclusion in GWAS on reproductive traits (scrotal circumference [SC], early pregnancy [P16] and age at first calving [AFC]) in a Nelore herd. Genotype data from 3,263 animals with the above-mentioned phenotypes were used. The results showed an increase in the variances explained by the autosomal markers for all traits when the X chromosome was not included. For SC, there was an increase of more than 10% for the windows on chromosomes 2 and 6. For P16, the effect was increased by almost 20% for windows on chromosome 5. The same pattern was found for AFC, with an increase of more than 10% for the most important windows. The results indicate that the noninclusion of the X chromosome can overestimate the effects of autosomes on SC, P16 and AFC not only because of the additive effect of the X chromosome itself but also because of its epistatic effect on autosomal genes.

Theoretical and empirical comparisons of expected and realized relationships for the X-chromosome

Background

X-chromosomal loci present different inheritance patterns compared to autosomal loci and must be modeled accordingly. Sexual chromosomes are not systematically considered in whole-genome relationship matrices although rules based on genealogical or marker information have been derived. Loci on the X-chromosome could have a significant contribution to the additive genetic variance, in particular for some traits such as those related to reproduction. Thus, accounting for the X-chromosome relationship matrix might be informative to better understand the architecture of complex traits (e.g., by estimating the variance associated to this chromosome) and to improve their genomic prediction. For such applications, previous studies have shown the benefits of combining information from genotyped and ungenotyped individuals.

Results

In this paper, we start by presenting rules to compute a genomic relationship matrix (GRM) for the X-chromosome (G^X) without making any assumption on dosage compensation, and based on coding of gene content with 0/1 for males and 0/1/2 for females. This coding adjusts naturally to previously derived pedigree-based relationships (S) for the X-chromosome. When needed, we propose to accommodate and estimate dosage compensation and genetic heterogeneity across sexes via multiple trait models. Using a Holstein dairy cattle dataset, including males and females, we then empirically illustrate that realized relationships (G^X) matches expectations (S). However, G^X presents high deviations from S. G^X has also a lower dimensionality compared to the autosomal GRM. In particular, individuals are frequently identical along the entire chromosome. Finally, we confirm that the heritability of gene content for markers on the X-chromosome that are estimated by using S is 1, further demonstrating that S and G^X can be combined. For the pseudo-autosomal region, we demonstrate that the expected relationships vary according to position because of the sex-gradient. We end by presenting the rules to construct the 'H matrix’ by combining both relationship matrices.

Conclusions

This work shows theoretically and empirically that a pedigree-based relationship matrix built with rules specifically developed for the X-chromosome (S) matches the realized GRM for the X-chromosome. Therefore, applications that combine expected relationships and genotypes for markers on the X-chromosome should use S and G^X.

Association of G6PD variants with hemoglobin A1c and impact on diabetes diagnosis in East Asian individuals

OBJECTIVE

Hemoglobin A1c (HbA1c) accuracy is important for diabetes diagnosis and estimation of overall glycemia. The G6PD-Asahi variant which causes glucose-6-phosphate dehydrogenase (G6PD) deficiency has been shown to lower HbA1c independently of glycemia in African ancestry populations. As different G6PD variants occur in Asian ancestry, we sought to identify Asian-specific G6PD variants associated with HbA1c.

RESEARCH DESIGN AND METHODS

In eight Asian population-based cohorts, we performed imputation on the X chromosome using the 1000 Genomes reference panel and tested for association with HbA1c (10 005 East Asians and 2051 South Asians). Results were meta-analyzed across studies. We compared the proportion of individuals classified as having diabetes/pre-diabetes by fasting glucose ≥100 mg/dL or HbA1c ≥5.7% units among carriers and non-carriers of HbA1c-associated variants.

RESULTS

The strongest association was a missense variant (G6PD-Canton, rs72554665, minor allele frequency=2.2%, effect in men=-0.76% unit, 95% CI -0.88 to -0.64, p=1.25×10^-27, n=2844). Conditional analyses identified a secondary distinct signal, missense variant (G6PD-Kaiping, rs72554664, minor allele frequency=1.6%, effect in men=-1.12 % unit, 95% CI -1.32 to -0.92, p=3.12×10^-15, p_{conditional_Canton}=7.57×10^-11). Adjusting for glucose did not attenuate their effects. The proportion of individuals with fasting glucose ≥100 mg/dL did not differ by carrier status of G6PD-Canton (p=0.21). Whereas the proportion of individuals with HbA1c ≥5.7% units was lower in carriers (5%) compared with non-carriers of G6PD-Canton (30%, p=0.03).

CONCLUSIONS

We identified two G6PD variants in East Asian men associated with non-glycemic lowering of HbA1c. Carriers of these variants are more likely to be underdiagnosed for diabetes or pre-diabetes than non-carriers if screened by HbA1c without confirmation by direct glucose measurements.

Statistical methods for testing X chromosome variant associations: application to sex-specific characteristics of bipolar disorder

BACKGROUND

Bipolar disorder (BD) affects both sexes, but important sex differences exist with respect to its symptoms and comorbidities. For example, rapid cycling (RC) is more prevalent in females, and alcohol use disorder (AUD) is more prevalent in males. We hypothesize that X chromosome variants may be associated with sex-specific characteristics of BD. Few studies have explored the role of the X chromosome in BD, which is complicated by X chromosome inactivation (XCI). This process achieves "dosage compensation" for many X chromosome genes by silencing one of the two copies in females, and most statistical methods either ignore that XCI occurs or falsely assume that one copy is inactivated at all loci. We introduce new statistical methods that do not make these assumptions.

METHODS

We investigated this hypothesis in 1001 BD patients from the Genetic Association Information Network (GAIN) and 957 BD patients from the Mayo Clinic Bipolar Disorder Biobank. We examined the association of over 14,000 X chromosome single nucleotide polymorphisms (SNPs) with sex-associated BD traits using two statistical approaches that account for whether a SNP may be undergoing or escaping XCI. In the "XCI-informed approach," we fit a sex-adjusted logistic regression model assuming additive genetic effects where we coded the SNP either assuming one copy is expressed or two copies are expressed based on prior knowledge about which regions are inactivated. In the "XCI-robust approach," we fit a logistic regression model with sex, SNP, and SNP-sex interaction effects that is flexible to whether the region is inactivated or escaping XCI.

RESULTS

Using the "XCI-informed approach," which considers only the main effect of SNP and does not allow the SNP effect to differ by sex, no significant associations were identified for any of the phenotypes. Using the "XCI-robust approach," intergenic SNP rs5932307 was associated with BD (P = 8.3 × 10-8), with a stronger effect in females (odds ratio in males (ORM) = 1.13, odds ratio in females for a change of two allele copies (ORW2) = 3.86).

CONCLUSION

X chromosome association studies should employ methods which account for its unique biology. Future work is needed to validate the identified associations with BD, to formally assess the performance of both approaches under different true genetic architectures, and to apply these approaches to study sex differences in other conditions.

Impact of Rare and Common Genetic Variants on Diabetes Diagnosis by Hemoglobin A1c in Multi-Ancestry Cohorts: The Trans-Omics for Precision Medicine Program

Hemoglobin A1c (HbA1c) is widely used to diagnose diabetes and assess glycemic control in individuals with diabetes. However, nonglycemic determinants, including genetic variation, may influence how accurately HbA1c reflects underlying glycemia. Analyzing the NHLBI Trans-Omics for Precision Medicine (TOPMed) sequence data in 10,338 individuals from five studies and four ancestries (6,158 Europeans, 3,123 African-Americans, 650 Hispanics, and 407 East Asians), we confirmed five regions associated with HbA1c (GCK in Europeans and African-Americans, HK1 in Europeans and Hispanics, FN3K and/or FN3KRP in Europeans, and G6PD in African-Americans and Hispanics) and we identified an African-ancestry-specific low-frequency variant (rs1039215 in HBG2 and HBE1, minor allele frequency (MAF) = 0.03). The most associated G6PD variant (rs1050828-T, p.Val98Met, MAF = 12% in African-Americans, MAF = 2% in Hispanics) lowered HbA1c (-0.88% in hemizygous males, -0.34% in heterozygous females) and explained 23% of HbA1c variance in African-Americans and 4% in Hispanics. Additionally, we identified a rare distinct G6PD coding variant (rs76723693, p.Leu353Pro, MAF = 0.5%; -0.98% in hemizygous males, -0.46% in heterozygous females) and detected significant association with HbA1c when aggregating rare missense variants in G6PD. We observed similar magnitude and direction of effects for rs1039215 (HBG2) and rs76723693 (G6PD) in the two largest TOPMed African American cohorts, and we replicated the rs76723693 association in the UK Biobank African-ancestry participants. These variants in G6PD and HBG2 were monomorphic in the European and Asian samples. African or Hispanic ancestry individuals carrying G6PD variants may be underdiagnosed for diabetes when screened with HbA1c. Thus, assessment of these variants should be considered for incorporation into precision medicine approaches for diabetes diagnosis.

Conservation, acquisition, and functional impact of sex-biased gene expression in mammals

Sex differences abound in human health and disease, as they do in other mammals used as models. The extent to which sex differences are conserved at the molecular level across species and tissues is unknown. We surveyed sex differences in gene expression in human, macaque, mouse, rat, and dog, across 12 tissues. In each tissue, we identified hundreds of genes with conserved sex-biased expression-findings that, combined with genomic analyses of human height, explain ~12% of the difference in height between females and males. We surmise that conserved sex biases in expression of genes otherwise operating equivalently in females and males contribute to sex differences in traits. However, most sex-biased expression arose during the mammalian radiation, which suggests that careful attention to interspecies divergence is needed when modeling human sex differences.