RESUMEN
Heparan sulfate (HS) proteoglycans are important regulators of cellular responses to soluble mediators such as chemokines, cytokines and growth factors. We profiled changes in expression of genes encoding HS core proteins, biosynthesis enzymes and modifiers during macrophage polarisation, and found that the most highly regulated gene was Sulf2, an extracellular HS 6-O-sulfatase that was markedly downregulated in response to pro-inflammatory stimuli. We then generated Sulf2+/- bone marrow chimeric mice and examined inflammatory responses in antigen-induced arthritis, as a model of rheumatoid arthritis. Resolution of inflammation was impaired in myeloid Sulf2+/- chimeras, with elevated joint swelling and increased abundance of pro-arthritic Th17 cells in synovial tissue. Transcriptomic and in vitro analyses indicated that Sulf2 deficiency increased type I interferon signaling in bone marrow-derived macrophages, leading to elevated expression of the Th17-inducing cytokine IL6. This establishes that dynamic remodeling of HS by Sulf2 limits type I interferon signaling in macrophages, and so protects against Th17-driven pathology.
Asunto(s)
Macrófagos , Ratones Endogámicos C57BL , Transducción de Señal , Células Th17 , Animales , Células Th17/inmunología , Células Th17/metabolismo , Ratones , Macrófagos/metabolismo , Macrófagos/inmunología , Sulfatasas/metabolismo , Sulfatasas/genética , Sulfotransferasas/metabolismo , Sulfotransferasas/genética , Células Mieloides/metabolismo , Células Mieloides/inmunología , Artritis Experimental/inmunología , Artritis Experimental/patología , Artritis Experimental/metabolismo , Artritis Reumatoide/inmunología , Artritis Reumatoide/metabolismo , Artritis Reumatoide/patología , Inflamación/metabolismo , Inflamación/patología , Ratones Noqueados , Interleucina-6/metabolismo , Interleucina-6/genética , Heparitina Sulfato/metabolismoRESUMEN
Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 600 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.
RESUMEN
Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 599 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.
RESUMEN
BACKGROUND: Genome sequence variants affecting complex traits (quantitative trait loci, QTL) are enriched in functional regions of the genome, such as those marked by certain histone modifications. These variants are believed to influence gene expression. However, due to the linkage disequilibrium among nearby variants, pinpointing the precise location of QTL is challenging. We aimed to identify allele-specific binding (ASB) QTL (asbQTL) that cause variation in the level of histone modification, as measured by the height of peaks assayed by ChIP-seq (chromatin immunoprecipitation sequencing). We identified DNA sequences that predict the difference between alleles in ChIP-seq peak height in H3K4me3 and H3K27ac histone modifications in the mammary glands of cows. RESULTS: We used a gapped k-mer support vector machine, a novel best linear unbiased prediction model, and a multiple linear regression model that combines the other two approaches to predict variant impacts on peak height. For each method, a subset of 1000 sites with the highest magnitude of predicted ASB was considered as candidate asbQTL. The accuracy of this prediction was measured by the proportion where the predicted direction matched the observed direction. Prediction accuracy ranged between 0.59 and 0.74, suggesting that these 1000 sites are enriched for asbQTL. Using independent data, we investigated functional enrichment in the candidate asbQTL set and three control groups, including non-causal ASB sites, non-ASB variants under a peak, and SNPs (single nucleotide polymorphisms) not under a peak. For H3K4me3, a higher proportion of the candidate asbQTL were confirmed as ASB when compared to the non-causal ASB sites (P < 0.01). However, these candidate asbQTL did not enrich for the other annotations, including expression QTL (eQTL), allele-specific expression QTL (aseQTL) and sites conserved across mammals (P > 0.05). CONCLUSIONS: We identified putatively causal sites for asbQTL using the DNA sequence surrounding these sites. Our results suggest that many sites influencing histone modifications may not directly affect gene expression. However, it is important to acknowledge that distinguishing between putative causal ASB sites and other non-causal ASB sites in high linkage disequilibrium with the causal sites regarding their impact on gene expression may be challenging due to limitations in statistical power.
Asunto(s)
Alelos , Secuenciación de Inmunoprecipitación de Cromatina , Histonas , Sitios de Carácter Cuantitativo , Animales , Bovinos/genética , Histonas/genética , Histonas/metabolismo , Secuenciación de Inmunoprecipitación de Cromatina/métodos , Polimorfismo de Nucleótido Simple , Código de Histonas , Desequilibrio de Ligamiento , Anotación de Secuencia Molecular , FemeninoRESUMEN
Most genetic variants associated with fertility in mammals fall in non-coding regions of the genome and it is unclear how these variants affect fertility. Here we use genome-wide association summary statistics for Heifer puberty (pubertal or not at 600 days) from 27,707 Bos indicus, Bos taurus and crossbred cattle; multi-trait GWAS signals from 2119 indicine cattle for four fertility traits, including days to calving, age at first calving, pregnancy status, and foetus age in weeks (assessed by rectal palpation of the foetus); and expression quantitative trait locus for whole blood from 489 indicine cattle, to identify 87 putatively functional genes affecting cattle fertility. Our analysis reveals a significant overlap between the set of cattle and previously reported human fertility-related genes, impling the existence of a shared pool of genes that regulate fertility in mammals. These findings are crucial for developing approaches to improve fertility in cattle and potentially other mammals.
Asunto(s)
Fertilidad , Estudio de Asociación del Genoma Completo , Sitios de Carácter Cuantitativo , Animales , Bovinos/genética , Fertilidad/genética , Estudio de Asociación del Genoma Completo/veterinaria , Femenino , Polimorfismo de Nucleótido SimpleRESUMEN
The causes of temporal fluctuations in adult traits are poorly understood. Here, we investigate the genetic determinants of within-person trait variability of 8 repeatedly measured anthropometric traits in 50,117 individuals from the UK Biobank. We found that within-person (non-directional) variability had a SNP-based heritability of 2-5% for height, sitting height, body mass index (BMI) and weight (P ≤ 2.4 × 10-3). We also analysed longitudinal trait change and show a loss of both average height and weight beyond about 70 years of age. A variant tracking the Alzheimer's risk APOE- E 4 allele (rs429358) was significantly associated with weight loss ( ß = -0.047 kg per yr, s.e. 0.007, P = 2.2 × 10-11), and using 2-sample Mendelian Randomisation we detected a relationship consistent with causality between decreased lumbar spine bone mineral density and height loss (bxy = 0.011, s.e. 0.003, P = 3.5 × 10-4). Finally, population-level variance quantitative trait loci (vQTL) were consistent with within-person variability for several traits, indicating an overlap between trait variability assessed at the population or individual level. Our findings help elucidate the genetic influence on trait-change within an individual and highlight disease risks associated with these changes.
Asunto(s)
Apolipoproteínas E , Estatura , Índice de Masa Corporal , Polimorfismo de Nucleótido Simple , Sitios de Carácter Cuantitativo , Adulto , Anciano , Femenino , Humanos , Masculino , Persona de Mediana Edad , Alelos , Enfermedad de Alzheimer/genética , Antropometría , Apolipoproteínas E/genética , Estatura/genética , Peso Corporal/genética , Densidad Ósea/genética , Estudio de Asociación del Genoma Completo , Estudios Longitudinales , Vértebras Lumbares , Análisis de la Aleatorización Mendeliana , Biobanco del Reino Unido , Reino UnidoRESUMEN
We develop a method, SBayesRC, that integrates genome-wide association study (GWAS) summary statistics with functional genomic annotations to improve polygenic prediction of complex traits. Our method is scalable to whole-genome variant analysis and refines signals from functional annotations by allowing them to affect both causal variant probability and causal effect distribution. We analyze 50 complex traits and diseases using â¼7 million common single-nucleotide polymorphisms (SNPs) and 96 annotations. SBayesRC improves prediction accuracy by 14% in European ancestry and up to 34% in cross-ancestry prediction compared to the baseline method SBayesR, which does not use annotations, and outperforms other methods, including LDpred2, LDpred-funct, MegaPRS, PolyPred-S and PRS-CSx. Investigation of factors affecting prediction accuracy identifies a significant interaction between SNP density and annotation information, suggesting whole-genome sequence variants with annotations may further improve prediction. Functional partitioning analysis highlights a major contribution of evolutionary constrained regions to prediction accuracy and the largest per-SNP contribution from nonsynonymous SNPs.
Asunto(s)
Estudio de Asociación del Genoma Completo , Anotación de Secuencia Molecular , Herencia Multifactorial , Polimorfismo de Nucleótido Simple , Herencia Multifactorial/genética , Estudio de Asociación del Genoma Completo/métodos , Humanos , Anotación de Secuencia Molecular/métodos , Genómica/métodos , Genoma Humano , Modelos GenéticosRESUMEN
BACKGROUND: Bovine lactoferrin (Lf) is an iron absorbing whey protein with antibacterial, antiviral, and antifungal activity. Lactoferrin is economically valuable and has an extremely variable concentration in milk, partly driven by environmental influences such as milking frequency, involution, or mastitis. A significant genetic influence has also been previously observed to regulate lactoferrin content in milk. Here, we conducted genetic mapping of lactoferrin protein concentration in conjunction with RNA-seq, ChIP-seq, and ATAC-seq data to pinpoint candidate causative variants that regulate lactoferrin concentrations in milk. RESULTS: We identified a highly-significant lactoferrin protein quantitative trait locus (pQTL), as well as a cis lactotransferrin (LTF) expression QTL (cis-eQTL) mapping to the LTF locus. Using ChIP-seq and ATAC-seq datasets representing lactating mammary tissue samples, we also report a number of regions where the openness of chromatin is under genetic influence. Several of these also show highly significant QTL with genetic signatures similar to those highlighted through pQTL and eQTL analysis. By performing correlation analysis between these QTL, we revealed an ATAC-seq peak in the putative promotor region of LTF, that highlights a set of 115 high-frequency variants that are potentially responsible for these effects. One of the 115 variants (rs110000337), which maps within the ATAC-seq peak, was predicted to alter binding sites of transcription factors known to be involved in lactation-related pathways. CONCLUSIONS: Here, we report a regulatory haplotype of 115 variants with conspicuously large impacts on milk lactoferrin concentration. These findings could enable the selection of animals for high-producing specialist herds.
Asunto(s)
Lactancia , Lactoferrina , Leche , Animales , Femenino , Haplotipos , Lactancia/genética , Lactoferrina/genética , Lactoferrina/análisis , Lactoferrina/metabolismo , Leche/química , Leche/metabolismo , BovinosRESUMEN
Many quantitative trait loci (QTLs) are in non-coding regions. Therefore, QTLs are assumed to affect gene regulation. Gene expression and RNA splicing are primary steps of transcription, so DNA variants changing gene expression (eVariants) or RNA splicing (sVariants) are expected to significantly affect phenotypes. We quantify the contribution of eVariants and sVariants detected from 16 tissues (n = 4,725) to 37 traits of â¼120,000 cattle (average magnitude of genetic correlation between traits = 0.13). Analyzed in Bayesian mixture models, averaged across 37 traits, cis and trans eVariants and sVariants detected from 16 tissues jointly explain 69.2% (SE = 0.5%) of heritability, 44% more than expected from the same number of random variants. This 69.2% includes an average of 24% from trans e-/sVariants (14% more than expected). Averaged across 56 lipidomic traits, multi-tissue cis and trans e-/sVariants also explain 71.5% (SE = 0.3%) of heritability, demonstrating the essential role of proximal and distal regulatory variants in shaping mammalian phenotypes.
RESUMEN
Allele specific expression (ASE) is widespread in many species including cows. Therefore, regulatory regions which control gene expression should show cis-regulatory variation which mirrors this differential expression within the animal. ChIP-seq peaks for histone modifications and transcription factors measure activity at functional regions and the height of some peaks have been shown to correlate across tissues with the expression of particular genes, suggesting these peaks are putative regulatory regions. In this study we identified ASE in the bovine genome in multiple tissues and investigated whether ChIP-seq peaks for four histone modifications and the transcription factor CTCF show allele specific binding (ASB) differences in the same tissues. We then investigate whether peak height and gene expression, which correlates across tissues, also correlates within the animal by investigating whether the direction of ASB in putative regulatory regions, mirrors that of the ASE in the genes they are putatively regulating. We found that ASE and ASB were widespread in the bovine genome but vary in extent between tissues. However, even when the height of a peak was positively correlated across tissues with expression of an exon, ASE of the exon and ASB of the peak were in the same direction only half the time. A likely explanation for this finding is that the correlations between peak height and exon expression do not indicate that the height of the peak causes the extent of exon expression, at least in some cases.
Asunto(s)
Código de Histonas , Factores de Transcripción , Femenino , Animales , Bovinos , Código de Histonas/genética , Factores de Transcripción/genética , Alelos , Secuenciación de Inmunoprecipitación de Cromatina , Exones/genéticaRESUMEN
Invariant natural killer T (iNKT) cells act at the interface between lipid metabolism and immunity because of their restriction to lipid antigens presented on CD1d by antigen-presenting cells (APCs). How foreign lipid antigens are delivered to APCs remains elusive. Since lipoproteins routinely bind glycosylceramides structurally similar to lipid antigens, we hypothesized that circulating lipoproteins form complexes with foreign lipid antigens. In this study, we used 2-color fluorescence correlation spectroscopy to show, for the first time to our knowledge, stable complex formation of lipid antigens α-galactosylceramide (αGalCer), isoglobotrihexosylceramide, and OCH, a sphingosine-truncated analog of αGalCer, with VLDL and/or LDL in vitro and in vivo. We demonstrate LDL receptor-mediated (LDLR-mediated) uptake of lipoprotein-αGalCer complexes by APCs, leading to potent complex-mediated activation of iNKT cells in vitro and in vivo. Finally, LDLR-mutant PBMCs of patients with familial hypercholesterolemia showed impaired activation and proliferation of iNKT cells upon stimulation, underscoring the relevance of lipoproteins as a lipid antigen delivery system in humans. Taken together, circulating lipoproteins form complexes with lipid antigens to facilitate their transport and uptake by APCs, leading to enhanced iNKT cell activation. This study thereby reveals a potentially novel mechanism of lipid antigen delivery to APCs and provides further insight into the immunological capacities of circulating lipoproteins.
Asunto(s)
Células T Asesinas Naturales , Humanos , Células Presentadoras de Antígenos , Lipoproteínas/metabolismoRESUMEN
Studies have demonstrated that structural variants (SV) play a substantial role in the evolution of species and have an impact on Mendelian traits in the genome. However, unlike small variants (< 50 bp), it has been challenging to accurately identify and genotype SV at the population scale using short-read sequencing. Long-read sequencing technologies are becoming competitively priced and can address several of the disadvantages of short-read sequencing for the discovery and genotyping of SV. In livestock species, analysis of SV at the population scale still faces challenges due to the lack of resources, high costs, technological barriers, and computational limitations. In this review, we summarize recent progress in the characterization of SV in the major livestock species, the obstacles that still need to be overcome, as well as the future directions in this growing field. It seems timely that research communities pool resources to build global population-scale long-read sequencing consortiums for the major livestock species for which the application of genomic tools has become cost-effective.
Asunto(s)
Genómica , Ganado , Animales , Ganado/genética , Genotipo , FenotipoRESUMEN
BACKGROUND: Causal variants for complex traits, such as eQTL are often found in non-coding regions of the genome, where they are hypothesised to influence phenotypes by regulating gene expression. Many regulatory regions are marked by histone modifications, which can be assayed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). Sequence reads from ChIP-seq form peaks at putative regulatory regions, which may reflect the amount of regulatory activity at this region. Therefore, eQTL which are also associated with differences in histone modifications are excellent candidate causal variants. RESULTS: We assayed the histone modifications H3K4Me3, H3K4Me1 and H3K27ac and mRNA in the mammary gland of up to 400 animals. We identified QTL for peak height (histone QTL), exon expression (eeQTL), allele specific expression (aseQTL) and allele specific binding (asbQTL). By intersecting these results, we identify variants which may influence gene expression by altering regulatory regions of the genome, and may be causal variants for other traits. Lastly, we find that these variants are found in putative transcription factor binding sites, identifying a mechanism for the effect of many eQTL. CONCLUSIONS: We find that allele specific and traditional QTL analysis often identify the same genetic variants and provide evidence that many eQTL are regulatory variants which alter activity at regulatory regions of the bovine genome. Our work provides methodological and biological updates on how regulatory mechanisms interplay at multi-omics levels.
Asunto(s)
Código de Histonas , Multiómica , Bovinos/genética , Animales , Variación Genética , Expresión GénicaRESUMEN
BACKGROUND: Sharing individual phenotype and genotype data between countries is complex and fraught with potential errors, while sharing summary statistics of genome-wide association studies (GWAS) is relatively straightforward, and thus would be especially useful for traits that are expensive or difficult-to-measure, such as feed efficiency. Here we examined: (1) the sharing of individual cow data from international partners; and (2) the use of sequence variants selected from GWAS of international cow data to evaluate the accuracy of genomic estimated breeding values (GEBV) for residual feed intake (RFI) in Australian cows. RESULTS: GEBV for RFI were estimated using genomic best linear unbiased prediction (GBLUP) with 50k or high-density single nucleotide polymorphisms (SNPs), from a training population of 3797 individuals in univariate to trivariate analyses where the three traits were RFI phenotypes calculated using 584 Australian lactating cows (AUSc), 824 growing heifers (AUSh), and 2526 international lactating cows (OVE). Accuracies of GEBV in AUSc were evaluated by either cohort-by-birth-year or fourfold random cross-validations. GEBV of AUSc were also predicted using only the AUS training population with a weighted genomic relationship matrix constructed with SNPs from the 50k array and sequence variants selected from a meta-GWAS that included only international datasets. The genomic heritabilities estimated using the AUSc, OVE and AUSh datasets were moderate, ranging from 0.20 to 0.36. The genetic correlations (rg) of traits between heifers and cows ranged from 0.30 to 0.95 but were associated with large standard errors. The mean accuracies of GEBV in Australian cows were up to 0.32 and almost doubled when either overseas cows, or both overseas cows and AUS heifers were included in the training population. They also increased when selected sequence variants were combined with 50k SNPs, but with a smaller relative increase. CONCLUSIONS: The accuracy of RFI GEBV increased when international data were used or when selected sequence variants were combined with 50k SNP array data. This suggests that if direct sharing of data is not feasible, a meta-analysis of summary GWAS statistics could provide selected SNPs for custom panels to use in genomic selection programs. However, since this finding is based on a small cross-validation study, confirmation through a larger study is recommended.
Asunto(s)
Bovinos , Lactancia , Animales , Australia , Bovinos/genética , Femenino , Estudio de Asociación del Genoma Completo , Genómica , Genotipo , Fenotipo , Polimorfismo de Nucleótido SimpleRESUMEN
The quantitative geneticist W. G. ('Bill') Hill, awardee of the 2018 Darwin Medal of the Royal Society and the 2019 Mendel Medal of the Genetics Society (United Kingdom), died on 17 December 2021 at the age of 81 years. Here, we pay tribute to his multiple key scientific contributions, which span population and evolutionary genetics, animal and plant breeding and human genetics. We discuss his theoretical research on the role of linkage disequilibrium (LD) and mutational variance in the response to selection, the origin of the widely used LD metric r2 in genomic association studies, the genetic architecture of complex traits, the quantification of the variation in realized relationships given a pedigree relationship and much more. We demonstrate that basic theoretical research in quantitative and statistical genetics has led to profound insights into the genetics and evolution of complex traits and made predictions that were subsequently empirically validated, often decades later.
Asunto(s)
Genoma , Fitomejoramiento , Animales , Estudio de Asociación del Genoma Completo , Genómica , Humanos , Desequilibrio de LigamientoRESUMEN
Bayesian methods, such as BayesR, for predicting the genetic value or risk of individuals from their genotypes, such as Single Nucleotide Polymorphisms (SNP), are often implemented using a Markov Chain Monte Carlo (MCMC) process. However, the generation of Markov chains is computationally slow. We introduce a form of blocked Gibbs sampling for estimating SNP effects from Markov chains that greatly reduces computational time by sampling each SNP effect iteratively n-times from conditional block posteriors. Subsequent iteration over all blocks m-times produces chains of length m × n. We use this strategy to solve large-scale genomic prediction and fine-mapping problems using the Bayesian MCMC mixed-effects genetic model, BayesR3. We validate the method using simulated data, followed by analysis of empirical dairy cattle data using high dimension milk mid infra-red spectra data as an example of "omics" data and show its use to increase the precision of mapping variants affecting milk, fat, and protein yields relative to a univariate analysis of milk, fat, and protein.
Asunto(s)
Genoma , Genómica , Animales , Teorema de Bayes , Bovinos , Genómica/métodos , Cadenas de Markov , FenotipoRESUMEN
AIMS: Inflammation is a key factor in atherosclerosis. The transcription factor interferon regulatory factor-5 (IRF5) drives macrophages towards a pro-inflammatory state. We investigated the role of IRF5 in human atherosclerosis and plaque stability. METHODS AND RESULTS: Bulk RNA sequencing from the Carotid Plaque Imaging Project biobank were used to mine associations between major macrophage associated genes and transcription factors and human symptomatic carotid disease. Immunohistochemistry, proximity extension assays, and Helios cytometry by time of flight (CyTOF) were used for validation. The effect of IRF5 deficiency on carotid plaque phenotype and rupture in ApoE-/- mice was studied in an inducible model of plaque rupture. Interferon regulatory factor-5 and ITGAX/CD11c were identified as the macrophage associated genes with the strongest associations with symptomatic carotid disease. Expression of IRF5 and ITGAX/CD11c correlated with the vulnerability index, pro-inflammatory plaque cytokine levels, necrotic core area, and with each other. Macrophages were the predominant CD11c-expressing immune cells in the plaque by CyTOF and immunohistochemistry. Interferon regulatory factor-5 immunopositive areas were predominantly found within CD11c+ areas with a predilection for the shoulder region, the area of the human plaque most prone to rupture. Accordingly, an inducible plaque rupture model of ApoE-/-Irf5-/- mice had significantly lower frequencies of carotid plaque ruptures, smaller necrotic cores, and less CD11c+ macrophages than their IRF5-competent counterparts. CONCLUSION: Using complementary evidence from data from human carotid endarterectomies and a murine model of inducible rupture of carotid artery plaque in IRF5-deficient mice, we demonstrate a mechanistic link between the pro-inflammatory transcription factor IRF5, macrophage phenotype, plaque inflammation, and its vulnerability to rupture.