Search

Abstract Phosphorylation plays a crucial role in the regulation of many fundamental cellular processes. Phosphorylation levels are increased in many cancer cells where they may promote changes in mitochondrial homeostasis. Proteomic studies on various types of cancer identified 17 phosphorylation sites within the human ATP-dependent protease Lon, which degrades misfolded, unassembled and oxidatively damaged proteins in mitochondria. Most of these sites were found in Lon’s N-terminal (NTD) and ATPase domains, though little is known about the effects on their function. By combining the biochemical and cryo-electron microscopy studies, we show the effect of Tyr186 and Tyr394 phosphorylations in Lon’s NTD, which greatly reduce all Lon activities without affecting its ability to bind substrates or perturbing its tertiary structure. A substantial reduction in Lon’s activities is also observed in the presence of polyphosphate, whose amount significantly increases in cancer cells. Our study thus provides an insight into the possible fine-tuning of Lon activities in human diseases, which highlights Lon’s importance in maintaining proteostasis in mitochondria.

In the bloodstream of mammalian hosts, African trypanosomes face the challenge of protecting their invariant surface receptors from immune detection. This crucial role is fulfilled by a dense, glycosylated protein layer composed of variant surface glycoproteins (VSGs), which undergo antigenic variation and provide a physical barrier that shields the underlying invariant surface glycoproteins (ISGs). The protective shield's limited permeability comes at the cost of restricted access to the extracellular host environment, raising questions regarding the specific function of the ISG repertoire. In this study, we employ an integrative structural biology approach to show that intrinsically disordered membrane-proximal regions are a common feature of members of the ISG super-family, conferring the ability to switch between compact and elongated conformers. While the folded, membrane-distal ectodomain is buried within the VSG layer for compact conformers, their elongated counterparts would enable the extension beyond it. This dynamic behavior enables ISGs to maintain a low immunogenic footprint while still allowing them to engage with the host environment when necessary. Our findings add further evidence to a dynamic molecular organization of trypanosome surface antigens wherein intrinsic disorder underpins the characteristics of a highly flexible ISG proteome to circumvent the constraints imposed by the VSG coat.

Abstract African Trypanosomes have developed elaborate mechanisms to escape the adaptive immune response, but little is known about complement evasion particularly at the early stage of infection. Here we show that ISG65 of the human-infective parasite Trypanosoma brucei gambiense is a receptor for human complement factor C3 and its activation fragments and that it takes over a role in selective inhibition of the alternative pathway C5 convertase and thus abrogation of the terminal pathway. No deposition of C4b, as part of the classical and lectin pathway convertases, was detected on trypanosomes. We present the cryo-electron microscopy (EM) structures of native C3 and C3b in complex with ISG65 which reveal a set of modes of complement interaction. Based on these findings, we propose a model for receptor-ligand interactions as they occur at the plasma membrane of blood-stage trypanosomes and may facilitate innate immune escape of the parasite.

The surface proteins of parasitic protozoa mediate functions essential to survival within a host, including nutrient accumulation, environmental sensing and immune evasion. Several receptors involved in nutrient uptake and defence from the innate immune response have been described in African trypanosomes and, together with antigenic variation, contribute towards persistence within vertebrate hosts. Significantly, a superfamily of invariant surface glycoproteins (ISGs) populates the trypanosome surface, one of which, ISG75, is implicated in uptake of the century-old drug suramin. By CRISPR/Cas9 knockout and biophysical analysis, we show here that ISG75 directly binds suramin and mediates uptake of additional naphthol-related compounds, making ISG75 a conduit for entry of at least one structural class of trypanocidal compounds. However, ISG75 null cells present only modest attenuation of suramin sensitivity, have unaltered viability in vivo and in vitro and no alteration to suramin-invoked proteome responses. While ISG75 is demonstrated as a valid suramin cell entry pathway, we suggest the presence of additional mechanisms for suramin accumulation, further demonstrating the complexity of trypanosomatid drug interactions and potential for evolution of resistance.

Abstract Identification of a protein minimal fragment amenable to crystallisation can be time- and labour intensive especially if large amounts are required and the protein has a complex fold and functionally important post-translational modifications. In addition, a lack of homologues and structural information can further complicate the design of a minimal expression construct. Recombinant expression in E. coli promises high yields, low costs and fast turnover times, but falls short for many extracellular, eukaryotic proteins. Eukaryotic expression systems provide an alternative but are costly, slow and require special handling and equipment. Using a member of a structurally uncharacterized, eukaryotic receptor family as an example we employ hydrogen–deuterium exchange mass spectrometry (HDX-MS) guided construct design in conjunction with truncation scanning and targeted expression host switching to identify a minimal expression construct that can be produced with high yields and moderate costs.

Salivarian trypanosomes comprise a group of extracellular anthroponotic and zoonotic parasites. The only sustainable method for global control of these infection is through vaccination of livestock animals. Despite multiple reports describing promising laboratory results, no single field-applicable solution has been successful so far. Conventionally, vaccine research focusses mostly on exposed immunogenic antigens, or the structural molecular knowledge of surface exposed invariant immunogens. Unfortunately, extracellular parasites (or parasites with extracellular life stages) have devised efficient defense systems against host antibody attacks, so they can deal with the mammalian humoral immune response. In the case of trypanosomes, it appears that these mechanisms have been perfected, leading to vaccine failure in natural hosts. Here, we provide two examples of potential vaccine candidates that, despite being immunogenic and accessible to the immune system, failed to induce a functionally protective memory response. First, trypanosomal enolase was tested as a vaccine candidate, as it was recently characterized as a highly conserved enzyme that is readily recognized during infection by the host antibody response. Secondly, we re-addressed a vaccine approach towards the Invariant Surface Glycoprotein ISG75, and showed that despite being highly immunogenic, trypanosomes can avoid anti-ISG75 mediated parasitemia control.

The major autolysins (Atl) of Staphylococcus epidermidis and S. aureus play an important role in cell separation, and their mutants are also attenuated in virulence. Therefore, autolysins represent a promising target for the development of new types of antibiotics. Here, we report the high-resolution structure of the catalytically active amidase domain AmiE (amidase S. epidermidis) from the major autolysin of S. epidermidis. This is the first protein structure with an amidase-like fold from a bacterium with a gram-positive cell wall architecture. AmiE adopts a globular fold, with several alpha-helices surrounding a central beta-sheet. Sequence comparison reveals a cluster of conserved amino acids that define a putative binding site with a buried zinc ion. Mutations of key residues in the putative active site result in loss of activity, enabling us to propose a catalytic mechanism. We also identified and synthesized muramyltripeptide, the minimal peptidoglycan fragment that can be used as a substrate by the enzyme. Molecular docking and digestion assays with muramyltripeptide derivatives allow us to identify key determinants of ligand binding. This results in a plausible model of interaction of this ligand not only for AmiE, but also for other PGN-hydrolases that share the same fold. As AmiE active-site mutations also show a severe growth defect, our findings provide an excellent platform for the design of specific inhibitors that target staphylococcal cell separation and can thereby prevent growth of this pathogen.

Summary: In African American adults, the strongest genetic predictor of pneumonia appears to be the A allele of rs334, a variant in the β-globin gene, which in homozygous form causes sickle cell disease (SCD). No comparable studies have been done in African American children. We performed genome-wide association analyses of 482 African American children with documented pneumonia and 2,048 African American control individuals using genotypes imputed from two reference panels: 1000 Genomes (1KG) (which contains rs334) and TOPMed (does not contain rs334). Using 1KG imputed genotypes, the most significant variant was rs334 (A allele; odds ratio [OR] = 2.76; 95% CI, 2.21–3.74; p = 5.9 × 10−19); using TOPMed imputed genotypes the most significant variant was rs2226952, found in the β-globin locus control region (G allele; OR = 2.14; 95% CI, 1.78–2.57; p = 5.1 × 10−16). After conditioning on rs334, the most strongly associated variant in the β-globin locus, rs33930165 (T allele, 1KG: OR = 4.09; 95% CI, 2.29–7.29; p = 1.7 × 10−6; TOPMed: OR = 3.58; 95% CI, 2.18–5.90; p = 4.7 × 10−7), which as a compound heterozygote with rs334 A allele, can cause SCD. To compare the power of different sample sets we developed a way to estimate the power of sample sets with different sample sizes, genotype arrays, and imputation platforms. Our results suggest that, in African American children, the strongest genetic determinants of pneumonia are those that increase the risk of SCD.

Abstract To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Abstract We assess performance and limitations of polygenic risk scores (PRSs) for multiple blood pressure (BP) phenotypes in diverse population groups. We compare “clumping-and-thresholding” (PRSice2) and LD-based (LDPred2) methods to construct PRSs from each of multiple GWAS, as well as multi-PRS approaches that sum PRSs with and without weights, including PRS-CSx. We use datasets from the MGB Biobank, TOPMed study, UK biobank, and from All of Us to train, assess, and validate PRSs in groups defined by self-reported race/ethnic background (Asian, Black, Hispanic/Latino, and White). For both SBP and DBP, the PRS-CSx based PRS, constructed as a weighted sum of PRSs developed from multiple independent GWAS, perform best across all race/ethnic backgrounds. Stratified analysis in All of Us shows that PRSs are better predictive of BP in females compared to males, individuals without obesity, and middle-aged (40-60 years) compared to older and younger individuals.

AbstractPopulation genetics has adapted as technological advances in next-generation sequencing have resulted in an exponential increase of genetic data. A common approach to efficiently analyze genetic variation present in large sequencing data is through the allele frequency spectrum, defined as the distribution of allele frequencies in a sample. While the frequency spectrum serves to summarize patterns of genetic variation, it implicitly assumes mutation types (A→C vs C→T) as interchangeable. However, mutations of different types arise and spread due to spatial and temporal variation in forces such as mutation rate and biased gene conversion that result in heterogeneity in the distribution of allele frequencies across sites. In this work, we explore the impact of this simplification on multiple aspects of population genetic modeling. As a site’s mutation rate is strongly affected by flanking nucleotides, we defined a mutation subtype by the base pair change and adjacent nucleotides (e.g. AAA→ATA) and systematically assessed the heterogeneity in the frequency spectrum across 96 distinct 3-mer mutation subtypes using n

Summary: Biobanks of linked clinical patient histories and biological samples are an efficient strategy to generate large cohorts for modern genetics research. Biobank recruitment varies by factors such as geographic catchment and sampling strategy, which affect biobank demographics and research utility. Here, we describe the Michigan Genomics Initiative (MGI), a single-health-system biobank currently consisting of >91,000 participants recruited primarily during surgical encounters at Michigan Medicine. The surgical enrollment results in a biobank enriched for many diseases and ideally suited for a disease genetics cohort. Compared with the much larger population-based UK Biobank, MGI has higher prevalence for nearly all diagnosis-code-based phenotypes and larger absolute case counts for many phenotypes. Genome-wide association study (GWAS) results replicate known findings, thereby validating the genetic and clinical data. Our results illustrate that opportunistic biobank sampling within single health systems provides a unique and complementary resource for exploring the genetics of complex diseases.

Summary: Polygenic risk scores (PRSs) have been widely explored in precision medicine. However, few studies have thoroughly investigated their best practices in global populations across different diseases. We here utilized data from Global Biobank Meta-analysis Initiative (GBMI) to explore methodological considerations and PRS performance in 9 different biobanks for 14 disease endpoints. Specifically, we constructed PRSs using pruning and thresholding (P + T) and PRS-continuous shrinkage (CS). For both methods, using a European-based linkage disequilibrium (LD) reference panel resulted in comparable or higher prediction accuracy compared with several other non-European-based panels. PRS-CS overall outperformed the classic P + T method, especially for endpoints with higher SNP-based heritability. Notably, prediction accuracy is heterogeneous across endpoints, biobanks, and ancestries, especially for asthma, which has known variation in disease prevalence across populations. Overall, we provide lessons for PRS construction, evaluation, and interpretation using GBMI resources and highlight the importance of best practices for PRS in the biobank-scale genomics era.

Fibromuscular dysplasia is a cardiovascular disease affecting mostly women with a mostly unknown genetic basis. Here the authors have performed a genome-wide association meta-analysis of Fibromuscular dysplasia to identify genetic loci, some of which are shared with common cardiovascular disease and traits.

Drugs targeting cardiovascular disease (CVD) can have negative consequences for liver function. Here, the authors combine genome wide analyses on 69,479 individuals to identify loss-of-function variants with beneficial effects on CVD-related traits without negative impacts on liver function.

Abstract Background The spectrum of somatic single-nucleotide variants in cancer genomes often reflects the signatures of multiple distinct mutational processes, which can provide clinically actionable insights into cancer etiology. Existing software tools for identifying and evaluating these mutational signatures do not scale to analyze large datasets containing thousands of individuals or millions of variants. Results We introduce Helmsman, a program designed to perform mutation signature analysis on arbitrarily large sequencing datasets. Helmsman is up to 300 times faster than existing software. Helmsman’s memory usage is independent of the number of variants, resulting in a small enough memory footprint to analyze datasets that would otherwise exceed the memory limitations of other programs. Conclusions Helmsman is a computationally efficient tool that enables users to evaluate mutational signatures in massive sequencing datasets that are otherwise intractable with existing software. Helmsman is freely available at https://github.com/carjed/helmsman.

Abstract Background Disentangling the etiology of common, complex diseases is a major challenge in genetic research. For bipolar disorder (BD), several genome-wide association studies (GWAS) have been performed. Similar to other complex disorders, major breakthroughs in explaining the high heritability of BD through GWAS have remained elusive. To overcome this dilemma, genetic research into BD, has embraced a variety of strategies such as the formation of large consortia to increase sample size and sequencing approaches. Here we advocate a complementary approach making use of already existing GWAS data: a novel data mining procedure to identify yet undetected genotype–phenotype relationships. We adapted association rule mining, a data mining technique traditionally used in retail market research, to identify frequent and characteristic genotype patterns showing strong associations to phenotype clusters. We applied this strategy to three independent GWAS datasets from 2835 phenotypically characterized patients with BD. In a discovery step, 20,882 candidate association rules were extracted. Results Two of these rules—one associated with eating disorder and the other with anxiety—remained significant in an independent dataset after robust correction for multiple testing. Both showed considerable effect sizes (odds ratio ~ 3.4 and 3.0, respectively) and support previously reported molecular biological findings. Conclusion Our approach detected novel specific genotype–phenotype relationships in BD that were missed by standard analyses like GWAS. While we developed and applied our method within the context of BD gene discovery, it may facilitate identifying highly specific genotype–phenotype relationships in subsets of genome-wide data sets of other complex phenotype with similar epidemiological properties and challenges to gene discovery efforts.

Germline mutation rate is a critical parameter in the study of genetics and evolution. Here, Carlson et al. infer fine-scale patterns of human mutation rate heterogeneity by analyzing ~36 million singleton variants from 3560 whole-genome sequences.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Background: Alzheimer’s disease (AD) and bipolar disorder (BIP) are complex traits influenced by numerous common genetic variants, most of which remain to be detected. Clinical and epidemiological evidence suggest that AD and BIP are related. However, it is not established if this relation is of genetic origin. Here, we applied statistical methods based on the conditional false discovery rate (FDR) framework to detect genetic overlap between AD and BIP and utilized this overlap to increase the power to identify common genetic variants associated with either or both traits.Methods: We obtained genome wide association studies data from the International Genomics of Alzheimer’s Project part 1 (17,008 AD cases and 37,154 controls) and the Psychiatric Genetic Consortium Bipolar Disorder Working Group (20,352 BIP cases and 31,358 controls). We used conditional QQ-plots to assess overlap in common genetic variants between AD and BIP. We exploited the genetic overlap to re-rank test-statistics for AD and BIP and improve detection of genetic variants using the conditional FDR framework.Results: Conditional QQ-plots demonstrated a polygenic overlap between AD and BIP. Using conditional FDR, we identified one novel genomic locus associated with AD, and nine novel loci associated with BIP. Further, we identified two novel loci jointly associated with AD and BIP implicating the MARK2 gene (lead SNP rs10792421, conjunctional FDR = 0.030, same direction of effect) and the VAC14 gene (lead SNP rs11649476, conjunctional FDR = 0.022, opposite direction of effect).Conclusion: We found polygenic overlap between AD and BIP and identified novel loci for each trait and two jointly associated loci. Further studies should examine if the shared loci implicating the MARK2 and VAC14 genes could explain parts of the shared and distinct features of AD and BIP.

Multivariate statistical techniques such as principal components analysis (PCA) and multidimensional scaling (MDS) have been widely used to summarize the structure of human genetic variation, often in easily visualized two-dimensional maps. Many recent studies have reported similarity between geographic maps of population locations and MDS or PCA maps of genetic variation inferred from single-nucleotide polymorphisms (SNPs). However, this similarity has been evident primarily in a qualitative sense; and, because different multivariate techniques and marker sets have been used in different studies, it has not been possible to formally compare genetic variation datasets in terms of their levels of similarity with geography. In this study, using genome-wide SNP data from 128 populations worldwide, we perform a systematic analysis to quantitatively evaluate the similarity of genes and geography in different geographic regions. For each of a series of regions, we apply a Procrustes analysis approach to find an optimal transformation that maximizes the similarity between PCA maps of genetic variation and geographic maps of population locations. We consider examples in Europe, Sub-Saharan Africa, Asia, East Asia, and Central/South Asia, as well as in a worldwide sample, finding that significant similarity between genes and geography exists in general at different geographic levels. The similarity is highest in our examples for Asia and, once highly distinctive populations have been removed, Sub-Saharan Africa. Our results provide a quantitative assessment of the geographic structure of human genetic variation worldwide, supporting the view that geography plays a strong role in giving rise to human population structure.

Although a highly heritable and disabling disease, bipolar disorder's (BD) genetic variants have been challenging to identify. We present new genotype data for 1,190 cases and 401 controls and perform a genome-wide association study including additional samples for a total of 2,191 cases and 1,434 controls. We do not detect genome-wide significant associations for individual loci; however, across all SNPs, we show an association between the power to detect effects calculated from a previous genome-wide association study and evidence for replication (P = 1.5×10(-7)). To demonstrate that this result is not likely to be a false positive, we analyze replication rates in a large meta-analysis of height and show that, in a large enough study, associations replicate as a function of power, approaching a linear relationship. Within BD, SNPs near exons exhibit a greater probability of replication, supporting an enrichment of reproducible associations near functional regions of genes. These results indicate that there is likely common genetic variation associated with BD near exons (±10 kb) that could be identified in larger studies and, further, provide a framework for assessing the potential for replication when combining results from multiple studies.