Continent-wide genomic analysis of the African buffalo (Syncerus caffer).

The African buffalo (Syncerus caffer) is a wild bovid with a historical distribution across much of sub-Saharan Africa. Genomic analysis can provide insights into the evolutionary history of the species, and the key selective pressures shaping populations, including assessment of population level differentiation, population fragmentation, and population genetic structure. In this study we generated the highest quality de novo genome assembly (2.65 Gb, scaffold N50 69.17 Mb) of African buffalo to date, and sequenced a further 195 genomes from across the species distribution. Principal component and admixture analyses provided little support for the currently described four subspecies. Estimating Effective Migration Surfaces analysis suggested that geographical barriers have played a significant role in shaping gene flow and the population structure. Estimated effective population sizes indicated a substantial drop occurring in all populations 5-10,000 years ago, coinciding with the increase in human populations. Finally, signatures of selection were enriched for key genes associated with the immune response, suggesting infectious disease exert a substantial selective pressure upon the African buffalo. These findings have important implications for understanding bovid evolution, buffalo conservation and population management.

Genomic Reference Resource for African Cattle: Genome Sequences and High-Density Array Variants.

The diversity in genome resources is fundamental to designing genomic strategies for local breed improvement and utilisation. These resources also support gene discovery and enhance our understanding of the mechanisms of resilience with applications beyond local breeds. Here, we report the genome sequences of 555 cattle (208 of which comprise new data) and high-density (HD) array genotyping of 1,082 samples (537 new samples) from indigenous African cattle populations. The new sequences have an average genome coverage of ~30X, three times higher than the average (~10X) of the over 300 sequences already in the public domain. Following variant quality checks, we identified approximately 32.3 million sequence variants and 661,943 HD autosomal variants mapped to the Bos taurus reference genome (ARS-UCD1.2). The new datasets were generated as part of the Centre for Tropical Livestock Genetics and Health (CTLGH) Genomic Reference Resource for African Cattle (GRRFAC) initiative, which aspires to facilitate the generation of this livestock resource and hopes for its utilisation for complete indigenous breed characterisation and sustainable global livestock improvement.

Investigating the Role of Brain Natriuretic Peptide (BNP) and N-Terminal-proBNP in Thrombosis and Acute Ischemic Stroke Etiology.

The need for biomarkers for acute ischemic stroke (AIS) to understand the mechanisms implicated in pathological clot formation is critical. The levels of the brain natriuretic peptides known as brain natriuretic peptide (BNP) and NT-proBNP have been shown to be increased in patients suffering from heart failure and other heart conditions. We measured their expression in AIS clots of cardioembolic (CE) and large artery atherosclerosis (LAA) etiology, evaluating their location inside the clots, aiming to uncover their possible role in thrombosis. We analyzed 80 thrombi from 80 AIS patients in the RESTORE registry of AIS clots, 40 of which were of CE and 40 of LAA etiology. The localization of BNP and NT-BNP, quantified using immunohistochemistry and immunofluorescence, in AIS-associated white blood cell subtypes was also investigated. We found a statistically significant positive correlation between BNP and NT-proBNP expression levels (Spearman's rho = 0.668 p < 0.0001 *). We did not observe any statistically significant difference between LAA and CE clots in BNP expression (0.66 [0.13-3.54]% vs. 0.53 [0.14-3.07]%, p = 0.923) or in NT-proBNP expression (0.29 [0.11-0.58]% vs. 0.18 [0.05-0.51]%, p = 0.119), although there was a trend of higher NT-proBNP expression in the LAA clots. It was noticeable that BNP was distributed throughout the thrombus and especially within platelet-rich regions. However, NT-proBNP colocalized with neutrophils, macrophages, and T-lymphocytes, suggesting its association with the thrombo-inflammatory process.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds.

BACKGROUND: The genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome. RESULTS: We identified several thousand variants that had significantly steep clines ('SCV') across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle. CONCLUSIONS: To our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Meta-analysis of heritability estimates and genome-wide association for tick-borne haemoparasites in African cattle.

The control of tick-borne haemoparasites in cattle largely relies on the use of acaricide drugs against the tick vectors, with some vaccination also being used against selected pathogens. These interventions can be difficult in Africa, where accessibility and cost of vaccines can be issues, and the increasing resistance of tick vectors to the widely used acaricides is a complication to disease control. A potential complementary control strategy could be the exploitation of any natural host genetic resistance to the pathogens. However, there are currently very few estimates of the extent of host resistance to tick-borne haemoparasites, and a significant contributing factor to this knowledge gap is likely to be the difficulty of collecting appropriate samples and data in the smallholder systems that predominate livestock production in low- and middle-income countries, particularly at scale. In this study, we have estimated the heritability for the presence/absence of several important haemoparasite species (including Anaplasma marginale, Babesia bigemina, Babesia bovis, and Ehrlichia ruminantium), as well as for relevant traits such as body weight and body condition score (BCS), in 1,694 cattle from four African countries (Burkina Faso, Ghana, Nigeria, and Tanzania). Heritability estimates within countries were mostly not significant, ranging from 0.05 to 0.84 across traits and countries, with standard errors between 0.07 and 0.91. However, the weighted mean of heritability estimates was moderate and significant for body weight and BCS (0.40 and 0.49, respectively), with significant heritabilities also observed for the presence of A. marginale (0.16) and E. ruminantium (0.19). In a meta-analysis of genome-wide association studies (GWAS) for these traits, two peaks were identified as reaching the suggestive significance threshold (p < 1.91 × 10-7 and p < 1.89 × 10-7, respectively): one on chromosome 24 for BCS and one on chromosome 8 for the E. ruminantium infection status. These findings indicate that there is likely to be a genetic basis that contributes to pathogen presence/absence for tick-borne haemoparasite species, which could potentially be exploited to improve cattle resistance in Africa to the economically important diseases caused by these pathogens.

The Bovine Pangenome Consortium: democratizing production and accessibility of genome assemblies for global cattle breeds and other bovine species.

The Bovine Pangenome Consortium (BPC) is an international collaboration dedicated to the assembly of cattle genomes to develop a more complete representation of cattle genomic diversity. The goal of the BPC is to provide genome assemblies and a community-agreed pangenome representation to replace breed-specific reference assemblies for cattle genomics. The BPC invites partners sharing our vision to participate in the production of these assemblies and the development of a common, community-approved, pangenome reference as a public resource for the research community ( https://bovinepangenome.github.io/ ). This community-driven resource will provide the context for comparison between studies and the future foundation for cattle genomic selection.

Profiling the immune epigenome across global cattle breeds.

BACKGROUND: Understanding the variation between well and poorly adapted cattle breeds to local environments and pathogens is essential for breeding cattle with improved climate and disease-resistant phenotypes. Although considerable progress has been made towards identifying genetic differences between breeds, variation at the epigenetic and chromatin levels remains poorly characterized. Here, we generate, sequence and analyse over 150 libraries at base-pair resolution to explore the dynamics of DNA methylation and chromatin accessibility of the bovine immune system across three distinct cattle lineages. RESULTS: We find extensive epigenetic divergence between the taurine and indicine cattle breeds across immune cell types, which is linked to the levels of local DNA sequence divergence between the two cattle sub-species. The unique cell type profiles enable the deconvolution of complex cellular mixtures using digital cytometry approaches. Finally, we show distinct sub-categories of CpG islands based on their chromatin and methylation profiles that discriminate between classes of distal and gene proximal islands linked to discrete transcriptional states. CONCLUSIONS: Our study provides a comprehensive resource of DNA methylation, chromatin accessibility and RNA expression profiles of three diverse cattle populations. The findings have important implications, from understanding how genetic editing across breeds, and consequently regulatory backgrounds, may have distinct impacts to designing effective cattle epigenome-wide association studies in non-European breeds.

Rare genetic variants underlie outlying levels of DNA methylation and gene-expression.

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on DNA methylation (DNAm) as effect sizes are expected to be larger for molecular traits compared with complex traits. Here, we investigate DNAm in healthy ageing populations-the Lothian Birth Cohorts of 1921 and 1936-and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic single nucleotide polymorphisms (SNPs) within 1 kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels toward the tails of the population distribution. This study demonstrates the role of rare SNPs in the phenotypic variation of DNAm and the effect of extreme levels of DNAm on gene expression.

Assessment of genotyping array performance for genome-wide association studies and imputation in African cattle.

BACKGROUND: In cattle, genome-wide association studies (GWAS) have largely focused on European or Asian breeds, using genotyping arrays that were primarily designed for European cattle. Because there is growing interest in performing GWAS in African breeds, we have assessed the performance of 23 commercial bovine genotyping arrays for capturing the diversity across African breeds and performing imputation. We used 409 whole-genome sequences (WGS) spanning global cattle breeds, and a real cohort of 2481 individuals (including African breeds) that were genotyped with the Illumina high-density (HD) array and the GeneSeek bovine 50 k array. RESULTS: We found that commercially available arrays were not effective in capturing variants that segregate among African indicine animals. Only 6% of these variants in high linkage disequilibrium (LD) (r2 > 0.8) were on the best performing arrays, which contrasts with the 17% and 25% in African and European taurine cattle, respectively. However, imputation from available HD arrays can successfully capture most variants (accuracies up to 0.93), mainly when using a global, not continent-specific, reference panel, which partially reflects the unusually high levels of admixture on the continent. When considering functional variants, the GGPF250 array performed best for tagging WGS variants and imputation. Finally, we show that imputation from low-density arrays can perform almost as well as HD arrays, if a two-stage imputation approach is adopted, i.e. first imputing to HD and then to WGS, which can potentially reduce the costs of GWAS. CONCLUSIONS: Our results show that the choice of an array should be based on a balance between the objective of the study and the breed/population considered, with the HD and BOS1 arrays being the best choice for both taurine and indicine breeds when performing GWAS, and the GGPF250 being preferable for fine-mapping studies. Moreover, our results suggest that there is no advantage to using the indicus-specific arrays for indicus breeds, regardless of the objective. Finally, we show that using a reference panel that better represents global bovine diversity improves imputation accuracy, particularly for non-European taurine populations.

The conservation of human functional variants and their effects across livestock species.

Despite the clear potential of livestock models of human functional variants to provide important insights into the biological mechanisms driving human diseases and traits, their use to date has been limited. Generating such models via genome editing is costly and time consuming, and it is unclear which variants will have conserved effects across species. In this study we address these issues by studying naturally occurring livestock models of human functional variants. We show that orthologues of over 1.6 million human variants are already segregating in domesticated mammalian species, including several hundred previously directly linked to human traits and diseases. Models of variants linked to particular phenotypes, including metabolomic disorders and height, are preferentially shared across species, meaning studying the genetic basis of these phenotypes is particularly tractable in livestock. Using machine learning we demonstrate it is possible to identify human variants that are more likely to have an existing livestock orthologue, and, importantly, we show that the effects of functional variants are often conserved in livestock, acting on orthologous genes with the same direction of effect. Consequently, this work demonstrates the substantial potential of naturally occurring livestock carriers of orthologues of human functional variants to disentangle their functional impacts.

Graphia: A platform for the graph-based visualisation and analysis of high dimensional data.

Graphia is an open-source platform created for the graph-based analysis of the huge amounts of quantitative and qualitative data currently being generated from the study of genomes, genes, proteins metabolites and cells. Core to Graphia's functionality is support for the calculation of correlation matrices from any tabular matrix of continuous or discrete values, whereupon the software is designed to rapidly visualise the often very large graphs that result in 2D or 3D space. Following graph construction, an extensive range of measurement algorithms, routines for graph transformation, and options for the visualisation of node and edge attributes are available, for graph exploration and analysis. Combined, these provide a powerful solution for the interpretation of high-dimensional data from many sources, or data already in the form of a network or equivalent adjacency matrix. Several use cases of Graphia are described, to showcase its wide range of applications in the analysis biological data. Graphia runs on all major desktop operating systems, is extensible through the deployment of plugins and is freely available to download from https://graphia.app/.

Next generation sequencing of transcribed genes in ruminant γδ T cell populations.

Bovine γδ T cells are distinguished by expression of WC1, hybrid pattern recognition receptors and co-receptors to the T cell receptor (TCR), or their absence. WC1 molecules bind pathogens and the ability of γδ T cells to respond to pathogens largely correlates with their expression of particular WC1 genes. Following activation, the TCR and WC1 molecules co-localize and knocking down WC1 abrogates the ability of WC1-expressing γδ T cells to respond to antigen. It is known that these two major populations, WC1+ and WC1-, differ in their TCR gene expression and previous studies showed other differences using semi-quantitative RT-PCR and serial analysis of gene expression. Differences in genes expressed would influence the functional outcome when WC1+ vs. WC1- γδ T cells respond to pathogens. To identify unique aspects of their transcriptome, here we performed RNA-Seq of flow cytometrically sorted bovine WC1+ and WC1- γδ T cells and compared them to all mononuclear cells in blood. The greatest differences in gene expression were found between γδ T cells and other mononuclear cells and included those involved in lymphocyte activation and effector processes. Only minor differences occurred between ex vivo WC1+ vs. WC1- γδ T cells with those gene products being involved in cell adhesion and chemotaxis. After culturing cells from primed animals with Leptospira antigens major difference in the transcriptome was evident, with over 600 genes significantly differentially expressed including those focused on cytokine signaling. Unexpectedly, antigen-responding and non-responding populations of WC1+ γδ T cells had few differences in their transcriptomes outside of cytotoxic factors although they had more WC1-1, WC1-2 and WC1-13 transcripts. Through differential gene expression we were able to define properties of ex vivo and stimulated WC1+ cells which will be useful in understanding their functional biology.

A locus conferring tolerance to Theileria infection in African cattle.

East Coast fever, a tick-borne cattle disease caused by the Theileria parva parasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Bos indicus (Boran) cattle demonstrating heritable tolerance to infection with T. parva (h2 = 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region on bovine chromosome 15 that is significantly associated with survival outcome following T. parva exposure. Testing this locus in an independent cohort of animals replicates this association with survival following T. parva infection. A stop gained variant in a paralogue of the FAF1 gene in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa's most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection by T. parva.

Using machine learning to detect the differential usage of novel gene isoforms.

BACKGROUND: Differential isoform usage is an important driver of inter-individual phenotypic diversity and is linked to various diseases and traits. However, accurately detecting the differential usage of different gene transcripts between groups can be difficult, in particular in less well annotated genomes where the spectrum of transcript isoforms is largely unknown. RESULTS: We investigated whether machine learning approaches can detect differential isoform usage based purely on the distribution of reads across a gene region. We illustrate that gradient boosting and elastic net approaches can successfully identify large numbers of genes showing potential differential isoform usage between Europeans and Africans, that are enriched among relevant biological pathways and significantly overlap those identified by previous approaches. We demonstrate that diversity at the 3' and 5' ends of genes are primary drivers of these differences between populations. CONCLUSION: Machine learning methods can effectively detect differential isoform usage from read fraction data, and can provide novel insights into the biological differences between groups.

S100b in acute ischemic stroke clots is a biomarker for post-thrombectomy intracranial hemorrhages.

Background and purpose: Post-thrombectomy intracranial hemorrhages (PTIH) are dangerous complications of acute ischemic stroke (AIS) following mechanical thrombectomy. We aimed to investigate if S100b levels in AIS clots removed by mechanical thrombectomy correlated to increased risk of PTIH. Methods: We analyzed 122 thrombi from 80 AIS patients in the RESTORE Registry of AIS clots, selecting an equal number of patients having been pre-treated or not with rtPA (40 each group). Within each subgroup, 20 patients had developed PTIH and 20 patients showed no signs of hemorrhage. Gross photos of each clot were taken and extracted clot area (ECA) was measured using ImageJ. Immunohistochemistry for S100b was performed and Orbit Image Analysis was used for quantification. Immunofluorescence was performed to investigate co-localization between S100b and T-lymphocytes, neutrophils and macrophages. Chi-square or Kruskal-Wallis test were used for statistical analysis. Results: PTIH was associated with higher S100b levels in clots (0.33 [0.08-0.85] vs. 0.07 [0.02-0.27] mm2, H1 = 6.021, P = 0.014*), but S100b levels were not significantly affected by acute thrombolytic treatment (P = 0.386). PTIH was also associated with patients having higher NIHSS at admission (20.0 [17.0-23.0] vs. 14.0 [10.5-19.0], H1 = 8.006, P = 0.005) and higher number of passes during thrombectomy (2 [1-4] vs. 1 [1-2.5], H1 = 5.995, P = 0.014*). S100b co-localized with neutrophils, macrophages and with T-lymphocytes in the clots. Conclusions: Higher S100b expression in AIS clots, higher NIHSS at admission and higher number of passes during thrombectomy are all associated with PTIH. Further investigation of S100b expression in AIS clots by neutrophils, macrophages and T-lymphocytes could provide insight into the role of S100b in thromboinflammation.

Inherited Tolerance in Cattle to the Apicomplexan Protozoan Theileria parva is Associated with Decreased Proliferation of Parasite-Infected Lymphocytes.

Theileria parva is the causative agent of East Coast fever and Corridor disease, which are fatal, economically important diseases of cattle in eastern, central and southern Africa. Improved methods of control of the diseases are urgently required. The parasite transforms host lymphocytes, resulting in a rapid, clonal expansion of infected cells. Resistance to the disease has long been reported in cattle from T. parva-endemic areas. We reveal here that first- and second-generation descendants of a single Bos indicus bull survived severe challenge with T. parva, (overall survival rate 57.3% compared to 8.7% for unrelated animals) in a series of five field studies. Tolerant cattle displayed a delayed and less severe parasitosis and febrile response than unrelated animals. The in vitro proliferation of cells from surviving cattle was much reduced compared to those from animals that succumbed to infection. Additionally, some pro-inflammatory cytokines such as IL1ß, IL6, TNFα or TGFß which are usually strongly expressed in susceptible animals and are known to regulate cell growth or motility, remain low in tolerant animals. This correlates with the reduced proliferation and less severe clinical reactions observed in tolerant cattle. The results show for the first time that the inherited tolerance to T. parva is associated with decreased proliferation of infected lymphocytes. The results are discussed in terms of whether the reduced proliferation is the result of a perturbation of the transformation mechanism induced in infected cells or is due to an innate immune response present in the tolerant cattle.

Clinical Evaluation of Corridor Disease in Bos indicus (Boran) Cattle Naturally Infected With Buffalo-Derived Theileria parva.

Corridor disease (CD) is a fatal condition of cattle caused by buffalo-derived Theileria parva. Unlike the related condition, East Coast fever, which results from infection with cattle-derived T. parva, CD has not been extensively studied. We describe in detail the clinical and laboratory findings in cattle naturally infected with buffalo-derived T. parva. Forty-six cattle were exposed to buffalo-derived T. parva under field conditions at the Ol Pejeta Conservancy, Kenya, between 2013 and 2018. The first signs of disease observed in all animals were nasal discharge (mean day of onset was 9 days post-exposure), enlarged lymph nodes (10 days post-exposure), and pyrexia (13.7 days post-exposure). Coughing and labored breathing were observed in more than 50% of animals (14 days post-exposure). Less commonly observed signs, corneal edema (22%) and diarrhea (11%), were observed later in the disease progression (19 days post-exposure). All infections were considered clinically severe, and 42 animals succumbed to infection. The mean time to death across all studies was 18.4 days. The mean time from onset of clinical signs to death was 9 days and from pyrexia to death was 4.8 days, indicating a relatively short duration of clinical illness. There were significant relationships between days to death and the days to first temperature (chi2 = 4.00, p = 0.046), and days to peak temperature (chi2 = 25.81, p = 0.001), animals with earlier onset pyrexia died sooner. These clinical indicators may be useful for assessing the severity of disease in the future. All infections were confirmed by the presence of macroschizonts in lymph node biopsies (mean time to parasitosis was 11 days). Piroplasms were detected in the blood of two animals (4%) and 20 (43%) animals seroconverted. In this study, we demonstrate the successful approach to an experimental field study for CD in cattle. We also describe the clinical progression of CD in naturally infected cattle, including the onset and severity of clinical signs and pathology. Laboratory diagnoses based on examination of blood samples are unreliable, and alternatives may not be available to cattle keepers. The rapid development of CD requires recognition of the clinical signs, which may be useful for early diagnosis of the disease and effective intervention for affected animals.

Comparison of V̇o2peak and V̇o2max at Different Sampling Intervals in Collegiate Wrestlers.

ABSTRACT: Sell, KM, Ghigiarelli, JJ, Prendergast, JM, Ciani, GJ, Martin, J, and Gonzalez, AM. Comparison of V̇o2peak and V̇o2max at different sampling intervals in collegiate wrestlers. J Strength Cond Res 35(10): 2915-2917, 2021-The purpose of this study was to determine the difference in the highest oxygen uptake (V̇o2peak) achieved during a maximal effort graded exercise test (GXT) in 20 NCAA Division I male wrestlers using breath-by-breath (BbB) values to the maximal uptake averaged across different time- and breath-based oxygen consumption sampling intervals (V̇o2max). Given the need for aerobic fitness and anaerobic power in wrestling, the accurate determination of V̇o2max is imperative if it is to be used to identify current aerobic fitness and consequently guide sport-specific training programs to address weaknesses in this area. Each subject completed a cycle ergometer GXT during which BbB data were collected via indirect calorimetry and V̇o2peak determined as the highest value. V̇o2max was considered as the average value of 3-s, 5-s, 10-s, 20-s, and 30-s sampling, and 3-b, 7-b, and 11-b sampling during the GXT. Results show that the BbB V̇o2peak was significantly higher than the 5-s, 10-s, 20-s, 30-s, and 11-b (p < 0.05). The 3-b V̇o2max was significantly higher than the 20-s and 30-s V̇o2max values (p < 0.05). The underestimation of V̇o2peak for each time-based interval sampling approach compared with BbB V̇o2peak is consistent with previous research, but the comparison of BbB data to breath-based interval sampling has not been widely addressed in prior research. The use of a 7-b sampling interval for the determination of V̇o2max may be a promising approach to minimize the systematic errors associated with BbB or less frequent sampling intervals, but future research is needed to further support its application with elite athletic populations such as those in the current study.

Sex differences in genetic architecture in the UK Biobank.

Males and females present differences in complex traits and in the risk of a wide array of diseases. Genotype by sex (GxS) interactions are thought to account for some of these differences. However, the extent and basis of GxS are poorly understood. In the present study, we provide insights into both the scope and the mechanism of GxS across the genome of about 450,000 individuals of European ancestry and 530 complex traits in the UK Biobank. We found small yet widespread differences in genetic architecture across traits. We also found that, in some cases, sex-agnostic analyses may be missing trait-associated loci and looked into possible improvements in the prediction of high-level phenotypes. Finally, we studied the potential functional role of the differences observed through sex-biased gene expression and gene-level analyses. Our results suggest the need to consider sex-aware analyses for future studies to shed light onto possible sex-specific molecular mechanisms.