Start codon variant in LAG3 is associated with decreased LAG-3 expression and increased risk of autoimmune thyroid disease.

Autoimmune thyroid disease (AITD) is a common autoimmune disease. In a GWAS meta-analysis of 110,945 cases and 1,084,290 controls, 290 sequence variants at 225 loci are associated with AITD. Of these variants, 115 are previously unreported. Multiomics analysis yields 235 candidate genes outside the MHC-region and the findings highlight the importance of genes involved in T-cell regulation. A rare 5'-UTR variant (rs781745126-T, MAF = 0.13% in Iceland) in LAG3 has the largest effect (OR = 3.42, P = 2.2 × 10-16) and generates a novel start codon for an open reading frame upstream of the canonical protein translation initiation site. rs781745126-T reduces mRNA and surface expression of the inhibitory immune checkpoint LAG-3 co-receptor on activated lymphocyte subsets and halves LAG-3 levels in plasma among heterozygotes. All three homozygous carriers of rs781745126-T have AITD, of whom one also has two other T-cell mediated diseases, that is vitiligo and type 1 diabetes. rs781745126-T associates nominally with vitiligo (OR = 5.1, P = 6.5 × 10-3) but not with type 1 diabetes. Thus, the effect of rs781745126-T is akin to drugs that inhibit LAG-3, which unleash immune responses and can have thyroid dysfunction and vitiligo as adverse events. This illustrates how a multiomics approach can reveal potential drug targets and safety concerns.

Sequence variants affecting voice pitch in humans.

The genetic basis of the human vocal system is largely unknown, as are the sequence variants that give rise to individual differences in voice and speech. Here, we couple data on diversity in the sequence of the genome with voice and vowel acoustics in speech recordings from 12,901 Icelanders. We show how voice pitch and vowel acoustics vary across the life span and correlate with anthropometric, physiological, and cognitive traits. We found that voice pitch and vowel acoustics have a heritable component and discovered correlated common variants in ABCC9 that associate with voice pitch. The ABCC9 variants also associate with adrenal gene expression and cardiovascular traits. By showing that voice and vowel acoustics are influenced by genetics, we have taken important steps toward understanding the genetics and evolution of the human vocal system.

Genetic variants associated with syncope implicate neural and autonomic processes.

AIMS: Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications. METHODS AND RESULTS: This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders. CONCLUSION: The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope.

Multiomics analysis of rheumatoid arthritis yields sequence variants that have large effects on risk of the seropositive subset.

OBJECTIVES: To find causal genes for rheumatoid arthritis (RA) and its seropositive (RF and/or ACPA positive) and seronegative subsets. METHODS: We performed a genome-wide association study (GWAS) of 31 313 RA cases (68% seropositive) and ~1 million controls from Northwestern Europe. We searched for causal genes outside the HLA-locus through effect on coding, mRNA expression in several tissues and/or levels of plasma proteins (SomaScan) and did network analysis (Qiagen). RESULTS: We found 25 sequence variants for RA overall, 33 for seropositive and 2 for seronegative RA, altogether 37 sequence variants at 34 non-HLA loci, of which 15 are novel. Genomic, transcriptomic and proteomic analysis of these yielded 25 causal genes in seropositive RA and additional two overall. Most encode proteins in the network of interferon-alpha/beta and IL-12/23 that signal through the JAK/STAT-pathway. Highlighting those with largest effect on seropositive RA, a rare missense variant in STAT4 (rs140675301-A) that is independent of reported non-coding STAT4-variants, increases the risk of seropositive RA 2.27-fold (p=2.1×10-9), more than the rs2476601-A missense variant in PTPN22 (OR=1.59, p=1.3×10-160). STAT4 rs140675301-A replaces hydrophilic glutamic acid with hydrophobic valine (Glu128Val) in a conserved, surface-exposed loop. A stop-mutation (rs76428106-C) in FLT3 increases seropositive RA risk (OR=1.35, p=6.6×10-11). Independent missense variants in TYK2 (rs34536443-C, rs12720356-C, rs35018800-A, latter two novel) associate with decreased risk of seropositive RA (ORs=0.63-0.87, p=10-9-10-27) and decreased plasma levels of interferon-alpha/beta receptor 1 that signals through TYK2/JAK1/STAT4. CONCLUSION: Sequence variants pointing to causal genes in the JAK/STAT pathway have largest effect on seropositive RA, while associations with seronegative RA remain scarce.

Sequence Variants in TAAR5 and Other Loci Affect Human Odor Perception and Naming.

Olfactory receptor (OR) genes in humans form a special class characterized by unusually high DNA sequence diversity, which should give rise to differences in perception and behavior. In the largest genome-wide association study to date based on olfactory testing, we investigated odor perception and naming with smell tasks performed by 9,122 Icelanders, with replication in a separate sample of 2,204 individuals. We discovered an association between a low-frequency missense variant in TAAR5 and reduced intensity rating of fish odor containing trimethylamine (p.Ser95Pro, pcombined = 5.6 × 10-15). We demonstrate that TAAR5 genotype affects aversion to fish odor, reflected by linguistic descriptions of the odor and pleasantness ratings. We also discovered common sequence variants in two canonical olfactory receptor loci that associate with increased intensity and naming of licorice odor (trans-anethole: lead variant p.Lys233Asn in OR6C70, pcombined = 8.8 × 10-16 and pcombined = 1.4 × 10-9) and enhanced naming of cinnamon (trans-cinnamaldehyde; intergenic variant rs317787-T, pcombined = 5.0 × 10-17). Together, our results show that TAAR5 genotype variation influences human odor responses and highlight that sequence diversity in canonical OR genes can lead to enhanced olfactory ability, in contrast to the view that greater tolerance for mutations in the human OR repertoire leads to diminished function.

Structural modelling of the DNAJB6 oligomeric chaperone shows a peptide-binding cleft lined with conserved S/T-residues at the dimer interface.

The remarkably efficient suppression of amyloid fibril formation by the DNAJB6 chaperone is dependent on a set of conserved S/T-residues and an oligomeric structure, features unusual among DNAJ chaperones. We explored the structure of DNAJB6 using a combination of structural methods. Lysine-specific crosslinking mass spectrometry provided distance constraints to select a homology model of the DNAJB6 monomer, which was subsequently used in crosslink-assisted docking to generate a dimer model. A peptide-binding cleft lined with S/T-residues is formed at the monomer-monomer interface. Mixed isotope crosslinking showed that the oligomers are dynamic entities that exchange subunits. The purified protein is well folded, soluble and composed of oligomers with a varying number of subunits according to small-angle X-ray scattering (SAXS). Elongated particles (160 × 120 Å) were detected by electron microscopy and single particle reconstruction resulted in a density map of 20 Å resolution into which the DNAJB6 dimers fit. The structure of the oligomer and the S/T-rich region is of great importance for the understanding of the function of DNAJB6 and how it can bind aggregation-prone peptides and prevent amyloid diseases.

Chaperone-client interactions between Hsp21 and client proteins monitored in solution by small angle X-ray scattering and captured by crosslinking mass spectrometry.

The small heat shock protein (sHsp) chaperones are important for stress survival, yet the molecular details of how they interact with client proteins are not understood. All sHsps share a folded middle domain to which is appended flexible N- and C-terminal regions varying in length and sequence between different sHsps which, in different ways for different sHsps, mediate recognition of client proteins. In plants there is a chloroplast-localized sHsp, Hsp21, and a structural model suggests that Hsp21 has a dodecameric arrangement with six N-terminal arms located on the outside of the dodecamer and six inwardly-facing. Here, we investigated the interactions between Hsp21 and thermosensitive model substrate client proteins in solution, by small-angle X-ray scattering (SAXS) and crosslinking mass spectrometry. The chaperone-client complexes were monitored and the Rg -values were found to increase continuously during 20 min at 45°, which could reflect binding of partially unfolded clients to the flexible N-terminal arms of the Hsp21 dodecamer. No such increase in Rg -values was observed with a mutational variant of Hsp21, which is mainly dimeric and has reduced chaperone activity. Crosslinking data suggest that the chaperone-client interactions involve the N-terminal region in Hsp21 and only certain parts in the client proteins. These parts are peripheral structural elements presumably the first to unfold under destabilizing conditions. We propose that the flexible and hydrophobic N-terminal arms of Hsp21 can trap and refold early-unfolding intermediates with or without dodecamer dissociation.

The chloroplast-localized small heat shock protein Hsp21 associates with the thylakoid membranes in heat-stressed plants.

The small heat shock protein (sHsp) chaperones are crucial for cell survival and can prevent aggregation of client proteins that partially unfold under destabilizing conditions. Most investigations on the chaperone activity of sHsps are based on a limited set of thermosensitive model substrate client proteins since the endogenous targets are often not known. There is a high diversity among sHsps with a single conserved ß-sandwich fold domain defining the family, the α-crystallin domain, whereas the N-terminal and C-terminal regions are highly variable in length and sequence among various sHsps and conserved only within orthologues. The endogenous targets are probably also varying among various sHsps, cellular compartments, cell type and organism. Here we have investigated Hsp21, a non-metazoan sHsp expressed in the chloroplasts in green plants which experience huge environmental fluctuations not least in temperature. We describe how Hsp21 can also interact with the chloroplast thylakoid membranes, both when isolated thylakoid membranes are incubated with Hsp21 protein and when plants are heat-stressed. The amount of Hsp21 associated with the thylakoid membranes was precisely determined by quantitative mass spectrometry after metabolic 15 N-isotope labeling of either recombinantly expressed and purified Hsp21 protein or intact Arabidopsis thaliana plants. We found that Hsp21 is among few proteins that become associated with the thylakoid membranes in heat-stressed plants, and that approximately two thirds of the pool of chloroplast Hsp21 is affected. We conclude that for a complete picture of the role of sHsps in plant stress resistance also their association with the membranes should be considered.

Structural model of dodecameric heat-shock protein Hsp21: Flexible N-terminal arms interact with client proteins while C-terminal tails maintain the dodecamer and chaperone activity.

Small heat-shock proteins (sHsps) prevent aggregation of thermosensitive client proteins in a first line of defense against cellular stress. The mechanisms by which they perform this function have been hard to define due to limited structural information; currently, there is only one high-resolution structure of a plant sHsp published, that of the cytosolic Hsp16.9. We took interest in Hsp21, a chloroplast-localized sHsp crucial for plant stress resistance, which has even longer N-terminal arms than Hsp16.9, with a functionally important and conserved methionine-rich motif. To provide a framework for investigating structure-function relationships of Hsp21 and understanding these sequence variations, we developed a structural model of Hsp21 based on homology modeling, cryo-EM, cross-linking mass spectrometry, NMR, and small-angle X-ray scattering. Our data suggest a dodecameric arrangement of two trimer-of-dimer discs stabilized by the C-terminal tails, possibly through tail-to-tail interactions between the discs, mediated through extended IXVXI motifs. Our model further suggests that six N-terminal arms are located on the outside of the dodecamer, accessible for interaction with client proteins, and distinct from previous undefined or inwardly facing arms. To test the importance of the IXVXI motif, we created the point mutant V181A, which, as expected, disrupts the Hsp21 dodecamer and decreases chaperone activity. Finally, our data emphasize that sHsp chaperone efficiency depends on oligomerization and that client interactions can occur both with and without oligomer dissociation. These results provide a generalizable workflow to explore sHsps, expand our understanding of sHsp structural motifs, and provide a testable Hsp21 structure model to inform future investigations.

Probing the transient interaction between the small heat-shock protein Hsp21 and a model substrate protein using crosslinking mass spectrometry.

Small heat-shock protein chaperones are important players in the protein quality control system of the cell, because they can immediately respond to partially unfolded proteins, thereby protecting the cell from harmful aggregates. The small heat-shock proteins can form large polydisperse oligomers that are exceptionally dynamic, which is implicated in their function of protecting substrate proteins from aggregation. Yet the mechanism of substrate recognition remains poorly understood, and little is known about what parts of the small heat-shock proteins interact with substrates and what parts of a partially unfolded substrate protein interact with the small heat-shock proteins. The transient nature of the interactions that prevent substrate aggregation rationalize probing this interaction by crosslinking mass spectrometry. Here, we used a workflow with lysine-specific crosslinking and offline nano-liquid chromatography matrix-assisted laser desorption/ionization tandem time-of-flight mass spectrometry to explore the interaction between the plant small heat-shock protein Hsp21 and a thermosensitive model substrate protein, malate dehydrogenase. The identified crosslinks point at an interaction between the disordered N-terminal region of Hsp21 and the C-terminal presumably unfolding part of the substrate protein.

Detection of crosslinks within and between proteins by LC-MALDI-TOFTOF and the software FINDX to reduce the MSMS-data to acquire for validation.

Lysine-specific chemical crosslinking in combination with mass spectrometry is emerging as a tool for the structural characterization of protein complexes and protein-protein interactions. After tryptic digestion of crosslinked proteins there are thousands of peptides amenable to MSMS, of which only very few are crosslinked peptides of interest. Here we describe how the advantage offered by off-line LC-MALDI-TOF/TOF mass spectrometry is exploited in a two-step workflow to focus the MSMS-acquisition on crosslinks mainly. In a first step, MS-data are acquired and all the peak list files from the LC-separated fractions are merged by the FINDX software and screened for presence of crosslinks which are recognized as isotope-labeled doublet peaks. Information on the isotope doublet peak mass and intensity can be used as search constraints to reduce the number of false positives that match randomly to the observed peak masses. Based on the MS-data a precursor ion inclusion list is generated and used in a second step, where a restricted number of MSMS-spectra are acquired for crosslink validation. The decoupling of MS and MSMS and the peptide sorting with FINDX based on MS-data has the advantage that MSMS can be restricted to and focused on crosslinks of Type 2, which are of highest biological interest but often lowest in abundance. The LC-MALDI TOF/TOF workflow here described is applicable to protein multisubunit complexes and using (14)N/(15)N mixed isotope strategy for the detection of inter-protein crosslinks within protein oligomers.

Thiol-exchange in DTSSP crosslinked peptides is proportional to cysteine content and precisely controlled in crosslink detection by two-step LC-MALDI MSMS.

The lysine-specific crosslinker 3,3'-dithiobis(sulfosuccinimidylpropionate) (DTSSP) is commonly used in the structural characterization of proteins by chemical crosslinking and mass spectrometry and we here describe an efficient two-step LC-MALDI-TOF/TOF procedure to detect crosslinked peptides. First MS data are acquired, and the properties of isotope-labeled DTSSP are used in data analysis to identify candidate crosslinks. MSMS data are then acquired for a restricted number of precursor ions per spot for final crosslink identification. We show that the thiol-catalyzed exchange between crosslinked peptides, which is due to the disulfide bond in DTSSP and known to possibly obscure data, can be precisely quantified using isotope-labeled DTSSP. Crosslinked peptides are recognized as 8 Da doublet peaks and a new isotopic peak with twice the intensity appears in the middle of the doublet as a consequence of the thiol-exchange. False-positive crosslinks, formed exclusively by thiol-exchange, yield a 1:2:1 isotope pattern, whereas true crosslinks, formed by two lysine residues within crosslinkable distance in the native protein structure, yield a 1:0:1 isotope pattern. Peaks with a 1:X:1 isotope pattern, where 0 < X < 2, can be trusted as true crosslinks, with a defined proportion of the signal [2X/(2 + X)] being noise from the thiol-exchange. The thiol-exchange was correlated with the protein cysteine content and was minimized by shortening the trypsin incubation time, and for two molecular chaperone proteins with known structure all crosslinks fitted well to the structure data. The thiol-exchange can thus be controlled and isotope-labeled DTSSP safely used to detect true crosslinks between lysine residues in proteins.

The activity of barley NADPH-dependent thioredoxin reductase C is independent of the oligomeric state of the protein: tetrameric structure determined by cryo-electron microscopy.

Thioredoxin and thioredoxin reductase can regulate cell metabolism through redox regulation of disulfide bridges or through removal of H(2)O(2). These two enzymatic functions are combined in NADPH-dependent thioredoxin reductase C (NTRC), which contains an N-terminal thioredoxin reductase domain fused with a C-terminal thioredoxin domain. Rice NTRC exists in different oligomeric states, depending on the absence or presence of its NADPH cofactor. It has been suggested that the different oligomeric states may have diverse activity. Thus, the redox status of the chloroplast could influence the oligomeric state of NTRC and thereby its activity. We have characterized the oligomeric states of NTRC from barley (Hordeum vulgare L.). This also includes a structural model of the tetrameric NTRC derived from cryo-electron microscopy and single-particle reconstruction. We conclude that the tetrameric NTRC is a dimeric arrangement of two NTRC homodimers. Unlike that of rice NTRC, the quaternary structure of barley NTRC complexes is unaffected by addition of NADPH. The activity of NTRC was tested with two different enzyme assays. The N-terminal part of NTRC was tested in a thioredoxin reductase assay. A peroxide sensitive Mg-protoporphyrin IX monomethyl ester (MPE) cyclase enzyme system of the chlorophyll biosynthetic pathway was used to test the catalytic ability of both the N- and C-terminal parts of NTRC. The different oligomeric assembly states do not exhibit significantly different activities. Thus, it appears that the activities are independent of the oligomeric state of barley NTRC.