Mapping short tandem repeats for liver gene expression traits helps prioritize potential causal variants for complex traits in pigs

Background Short tandem repeats (STRs) were recently found to have significant impacts on gene expression and diseases in humans, but their roles on gene expression and complex traits in pigs remain unexplored. This study investigates the effects of STRs on gene expression in liver tissues based on the whole-genome sequences and RNA-Seq data of a discovery cohort of 260 F6 individuals and a validation population of 296 F7 individuals from a heterogeneous population generated from crosses among eight pig breeds. Results We identified 5203 and 5868 significantly expression STRs (eSTRs, FDR < 1%) in the F6 and F7 populations, respectively, most of which could be reciprocally validated (π1 = 0.92). The eSTRs explained 27.5% of the cis-heritability of gene expression traits on average. We further identified 235 and 298 fine-mapped STRs through the Bayesian fine-mapping approach in the F6 and F7 pigs, respectively, which were significantly enriched in intron, ATAC peak, compartment A and H3K4me3 regions. We identified 20 fine-mapped STRs located in 100 kb windows upstream and downstream of published complex trait-associated SNPs, which colocalized with epigenetic markers such as H3K27ac and ATAC peaks. These included eSTR of the CLPB, PGLS, PSMD6 and DHDH genes, which are linked with genome-wide association study (GWAS) SNPs for blood-related traits, leg conformation, growth-related traits, and meat quality traits, respectively. Conclusions This study provides insights into the effects of STRs on gene expression traits. The identified eSTRs are valuable resources for prioritizing causal STRs for complex traits in pigs.

SparsePro: an efficient genome-wide fine-mapping method integrating summary statistics and functional annotations

Identifying causal variants from genome-wide association studies (GWASs) is challenging due to widespread linkage disequilibrium (LD). Functional annotations of the genome may help prioritize variants that are biologically relevant and thus improve fine-mapping of GWAS results. However, classical fine-mapping methods have a high computational cost, particularly when the underlying genetic architecture and LD patterns are complex. Here, we propose a novel approach, SparsePro, to efficiently conduct genome-wide fine-mapping. Our method enjoys two major innovations: First, by creating a sparse low-dimensional projection of the high-dimensional genotype data, we enable a linear search of causal variants instead of a combinatorial search of causal configurations used in most existing methods; Second, we adopt a probabilistic framework with a highly efficient variational expectation-maximization algorithm to integrate statistical associations and functional priors. We evaluate SparsePro through extensive simulations using resources from the UK Biobank. Compared to state-of-the-art methods, SparsePro achieved more accurate and well-calibrated posterior inference with greatly reduced computation time. We demonstrate the utility of SparsePro by investigating the genetic architecture of five functional biomarkers of vital organs. We show that, compared to other methods, the causal variants identified by SparsePro are highly enriched for expression quantitative trait loci and explain a larger proportion of trait heritability. We also identify potential causal variants contributing to the genetically encoded coordination mechanisms between vital organs, and pinpoint target genes with potential pleiotropic effects. In summary, we have developed an efficient genome-wide fine-mapping method with the ability to integrate functional annotations. Our method may have wide utility in understanding the genetics of complex traits as well as in increasing the yield of functional follow-up studies of GWASs. SparsePro software is available on GitHub at https://github.com/zhwm/SparsePro.

Functional dissection of inherited non-coding variation influencing multiple myeloma risk

Thousands of non-coding variants have been associated with increased risk of human diseases, yet the causal variants and their mechanisms-of-action remain obscure. In an integrative study combining massively parallel reporter assays (MPRA), expression analyses (eQTL, meQTL, PCHiC) and chromatin accessibility analyses in primary cells (caQTL), we investigate 1,039 variants associated with multiple myeloma (MM). We demonstrate that MM susceptibility is mediated by gene-regulatory changes in plasma cells and B-cells, and identify putative causal variants at six risk loci (SMARCD3, WAC, ELL2, CDCA7L, CEP120, and PREX1). Notably, three of these variants co-localize with significant plasma cell caQTLs, signaling the presence of causal activity at these precise genomic positions in an endogenous chromosomal context in vivo. Our results provide a systematic functional dissection of risk loci for a hematologic malignancy.

Genetic Regulation of Cytokine Response in Patients with Acute Community-Acquired Pneumonia

Background: Community-acquired pneumonia (CAP) is an acute disease condition with a high risk of rapid deteriorations. We analysed the influence of genetics on cytokine regulation to obtain a better understanding of patient’s heterogeneity. Methods: For up to N = 389 genotyped participants of the PRO-GRESS study of hospitalised CAP patients, we performed a genome-wide association study of ten cytokines IL-1β, IL-6, IL-8, IL-10, IL-12, MCP-1 (MCAF), MIP-1α (CCL3), VEGF, VCAM-1, and ICAM-1. Consecutive secondary analyses were performed to identify independent hits and corresponding causal variants. Results: 102 SNPs from 14 loci showed genome-wide significant associations with five of the cytokines. The most interesting associations were found at 6p21.1 for VEGF (p = 1.58 × 10−20), at 17q21.32 (p = 1.51 × 10−9) and at 10p12.1 (p = 2.76 × 10−9) for IL-1β, at 10p13 for MIP-1α (CCL3) (p = 2.28 × 10−9), and at 9q34.12 for IL-10 (p = 4.52 × 10−8). Functionally plausible genes could be assigned to the majority of loci including genes involved in cytokine secretion, granulocyte function, and cilial kinetics. Conclusion: This is the first context-specific genetic association study of blood cytokine concentrations in CAP patients revealing numerous biologically plausible candidate genes. Two of the loci were also associated with atherosclerosis with probable common or consecutive pathomechanisms.

An exploration of linkage fine-mapping on sequences from case-control studies

Linkage analysis maps genetic loci for a heritable trait by identifying genomic regions with excess relatedness among individuals with similar trait values. Analysis may be conducted on related individuals from families, or on samples of unrelated individuals from a population. For allelically heterogeneous traits, population-based linkage analysis can be more powerful than genotypic-association analysis. Here, we focus on linkage analysis in a population sample, but use sequences rather than individuals as our unit of observation. Earlier investigations of sequence-based linkage mapping relied on known sequence relatedness, whereas we infer relatedness from the sequence data. We propose two ways to associate similarity in relatedness of sequences with similarity in their trait values and compare the resulting linkage methods to two genotypic- association methods. We also introduce a procedure to label case sequences as potential carriers or non-carriers of causal variants after an association has been found. This post-hoc labeling of case sequences is based on inferred relatedness to other case sequences. Our simulation results indicate that methods based on sequence-relatedness improve localization and perform as well as genotypic-association methods for detecting rare causal variants. Sequence-based linkage analysis therefore has potential to fine-map allelically heterogeneous disease traits.

Genetic overlap between idiopathic pulmonary fibrosis and COVID−19

Genome-wide association studies (GWAS) of coronavirus disease 2019 (COVID-19) and idiopathic pulmonary fibrosis (IPF) have identified genetic loci associated with both traits, suggesting possible shared biological mechanisms. Using updated GWAS of COVID-19 and IPF, we evaluated the genetic overlap between these two diseases and identified four genetic loci (including one novel) with likely shared causal variants between severe COVID-19 and IPF. Although there was a positive genetic correlation between COVID-19 and IPF, two of these four shared genetic loci had an opposite direction of effect. IPF-associated genetic variants related to telomere dysfunction and spindle assembly showed no association with COVID-19 phenotypes. Together, these results suggest there are both shared and distinct biological processes driving IPF and severe COVID-19 phenotypes.

Effect sizes of causal variants for gene expression and complex traits differ between populations

Despite the growing number of genome-wide association studies (GWAS) for complex traits, it remains unclear whether effect sizes of causal genetic variants differ between populations. In principle, effect sizes of causal variants could differ between populations due to gene-by-gene or gene-by-environment interactions. However, comparing causal variant effect sizes is challenging: it is difficult to know which variants are causal, and comparisons of variant effect sizes are confounded by differences in linkage disequilibrium (LD) structure between ancestries. Here, we develop a method to assess causal variant effect size differences that overcomes these limitations. Specifically, we leverage the fact that segments of European ancestry shared between European-American and admixed African-American individuals have similar LD structure, allowing for unbiased comparisons of variant effect sizes in European ancestry segments. We apply our method to two types of traits: gene expression and low-density lipoprotein cholesterol (LDL-C). We find that causal variant effect sizes for gene expression are significantly different between European-Americans and African-Americans; for LDL-C, we observe a similar point estimate although this is not significant, likely due to lower statistical power. Cross-population differences in variant effect sizes highlight the role of genetic interactions in trait architecture and will contribute to the poor portability of polygenic scores across populations, reinforcing the importance of conducting GWAS on individuals of diverse ancestries and environments.

Summary statistics knockoff inference empowers identification of putative causal variants in genome-wide association studies

Recent advances in genome sequencing and imputation technologies provide an exciting opportunity to comprehensively study the contribution of genetic variants to complex phenotypes. However, our ability to translate genetic discoveries into mechanistic insights remains limited at this point. In this paper, we propose an efficient knockoff-based method, GhostKnockoff, for genome-wide association studies (GWAS) that leads to improved power and ability to prioritize putative causal variants relative to conventional GWAS approaches. The method requires only Z-scores from conventional GWAS and hence can be easily applied to enhance existing and future studies. The method can also be applied to meta-analysis of multiple GWAS allowing for arbitrary sample overlap. We demonstrate its performance using empirical simulations and two applications: (1) analysis of 1,403 binary phenotypes from the UK Biobank data in 408,961 samples of European ancestry, and (2) a meta-analysis for Alzheimer’s disease (AD) comprising nine overlapping large-scale GWAS, whole-exome and whole-genome sequencing studies. The UK Biobank analysis demonstrates superior performance of the proposed method compared to conventional GWAS in both statistical power (2.05-fold more discoveries) and localization of putative causal variants at each locus (46% less proxy variants due to linkage disequilibrium). The AD meta-analysis identified 55 risk loci (including 31 new loci) with ~70% of the proximal genes at these loci showing suggestive signal in downstream single-cell transcriptomic analyses. Our results demonstrate that GhostKnockoff can identify putatively functional variants with weaker statistical effects that are missed by conventional association tests.

The adducin saga: Pleiotropic genomic targets for precision medicine in human hypertension; vascular, renal, and cognitive diseases

Hypertension is a leading risk factor for stroke, heart disease, chronic kidney disease, vascular cognitive impairment, and Alzheimer's disease. Previous genetic studies have nominated hundreds of genes linked to hypertension and renal and cognitive diseases. Some have been advanced as candidate genes by showing that they can alter blood pressure or renal and cerebral vascular function in knockout animals; however, final validation of the causal variants and underlying mechanisms have remained elusive. This review chronicles 40 years of work, from the initial identification of adducin (ADD) as an ACTIN-binding protein suggested to increase blood pressure in Milan hypertensive rats, to the discovery of a mutation in ADD1 as a candidate gene for hypertension in rats that were subsequently linked to hypertension in man. More recently, a recessive K572Q mutation in ADD3 was identified in Fawn-Hooded Hypertensive (FHH) and Milan Normotensive (MNS) rats that develop renal disease, which is absent in resistant strains. ADD3 dimerizes with ADD1 to form functional ADD protein. The mutation in ADD3 disrupts a critical ACTIN-binding site necessary for its interactions with actin and spectrin to regulate the cytoskeleton. Studies using Add3 knockout and transgenic strains, as well as a genetic complementation study in FHH and MNS rats, confirmed that the K572Q mutation in ADD3 plays a causal role in altering the myogenic response and autoregulation of renal and cerebral blood flow, resulting in increased susceptibility to hypertension-induced renal disease and cerebral vascular and cognitive dysfunction.

Genomic analysis of childhood hearing loss in the Yoruba population of Nigeria

Although variant alleles of hundreds of genes are associated with sensorineural deafness in children, the genes and alleles involved remain largely unknown in the Sub-Saharan regions of Africa. We ascertained 56 small families mainly of Yoruba ethno-lingual ancestry in or near Ibadan, Nigeria, that had at least one individual with nonsyndromic, severe-to-profound, prelingual-onset, bilateral hearing loss not attributed to nongenetic factors. We performed a combination of exome and Sanger sequencing analyses to evaluate both nuclear and mitochondrial genomes. No biallelic pathogenic variants were identified in GJB2, a common cause of deafness in many populations. Potential causative variants were identified in genes associated with nonsyndromic hearing loss (CIB2, COL11A1, ILDR1, MYO15A, TMPRSS3, and WFS1), nonsyndromic hearing loss or Usher syndrome (CDH23, MYO7A, PCDH15, and USH2A), and other syndromic forms of hearing loss (CHD7, OPA1, and SPTLC1). Several rare mitochondrial variants, including m.1555A>G, were detected in the gene MT-RNR1 but not in control Yoruba samples. Overall, 20 (33%) of 60 independent cases of hearing loss in this cohort of families were associated with likely causal variants in genes reported to underlie deafness in other populations. None of these likely causal variants were present in more than one family, most were detected as compound heterozygotes, and 77% had not been previously associated with hearing loss. These results indicate an unusually high level of genetic heterogeneity of hearing loss in Ibadan, Nigeria and point to challenges for molecular genetic screening, counseling, and early intervention in this population.