FinnGen provides genetic insights from a well-phenotyped isolated population.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

High-Resolution Genotyping of Formalin-Fixed Tissue Accurately Estimates Polygenic Risk Scores in Human Diseases.

Formalin-fixed paraffin-embedded (FFPE) tissues stored in biobanks and pathology archives are a vast but underutilized source for molecular studies on different diseases. Beyond being the "gold standard" for preservation of diagnostic human tissues, FFPE samples retain similar genetic information as matching blood samples, which could make FFPE samples an ideal resource for genomic analysis. However, research on this resource has been hindered by the perception that DNA extracted from FFPE samples is of poor quality. Here, we show that germline disease-predisposing variants and polygenic risk scores (PRS) can be identified from FFPE normal tissue (FFPE-NT) DNA with high accuracy. We optimized the performance of FFPE-NT DNA on a genome-wide array containing 657,675 variants. Via a series of testing and validation phases, we established a protocol for FFPE-NT genotyping with results comparable with blood genotyping. The median call rate of FFPE-NT samples in the validation phase was 99.85% (range 98.26%-99.94%) and median concordance with matching blood samples was 99.79% (range 98.85%-99.9%). We also demonstrated that a rare pathogenic PALB2 genetic variant predisposing to cancer can be correctly identified in FFPE-NT samples. We further imputed the FFPE-NT genotype data and calculated the FFPE-NT genome-wide PRS in 3 diseases and 4 disease risk variables. In all cases, FFPE-NT and matching blood PRS were highly concordant (all Pearson's r > 0.95). The ability to precisely genotype FFPE-NT on a genome-wide array enables translational genomics applications of archived FFPE-NT samples with the possibility to link to corresponding phenotypes and longitudinal health data.

Missense mutations in small muscle protein X-linked (SMPX) cause distal myopathy with protein inclusions.

Using deep phenotyping and high-throughput sequencing, we have identified a novel type of distal myopathy caused by mutations in the Small muscle protein X-linked (SMPX) gene. Four different missense mutations were identified in ten patients from nine families in five different countries, suggesting that this disease could be prevalent in other populations as well. Haplotype analysis of patients with similar ancestry revealed two different founder mutations in Southern Europe and France, indicating that the prevalence in these populations may be higher. In our study all patients presented with highly similar clinical features: adult-onset, usually distal more than proximal limb muscle weakness, slowly progressing over decades with preserved walking. Lower limb muscle imaging showed a characteristic pattern of muscle involvement and fatty degeneration. Histopathological and electron microscopic analysis of patient muscle biopsies revealed myopathic findings with rimmed vacuoles and the presence of sarcoplasmic inclusions, some with amyloid-like characteristics. In silico predictions and subsequent cell culture studies showed that the missense mutations increase aggregation propensity of the SMPX protein. In cell culture studies, overexpressed SMPX localized to stress granules and slowed down their clearance.

Mutation update: the spectra of nebulin variants and associated myopathies.

A mutation update on the nebulin gene (NEB) is necessary because of recent developments in analysis methodology, the identification of increasing numbers and novel types of variants, and a widening in the spectrum of clinical and histological phenotypes associated with this gigantic, 183 exons containing gene. Recessive pathogenic variants in NEB are the major cause of nemaline myopathy (NM), one of the most common congenital myopathies. Moreover, pathogenic NEB variants have been identified in core-rod myopathy and in distal myopathies. In this update, we present the disease-causing variants in NEB in 159 families, 143 families with NM, and 16 families with NM-related myopathies. Eighty-eight families are presented here for the first time. We summarize 86 previously published and 126 unpublished variants identified in NEB. Furthermore, we have analyzed the NEB variants deposited in the Exome Variant Server (http://evs.gs.washington.edu/EVS/), identifying that pathogenic variants are a minor fraction of all coding variants (∼7%). This indicates that nebulin tolerates substantial changes in its amino acid sequence, providing an explanation as to why variants in such a large gene result in relatively rare disorders. Lastly, we discuss the difficulties of drawing reliable genotype-phenotype correlations in NEB-associated disease.

Abnormal actin binding of aberrant ß-tropomyosins is a molecular cause of muscle weakness in TPM2-related nemaline and cap myopathy.

NM (nemaline myopathy) is a rare genetic muscle disorder defined on the basis of muscle weakness and the presence of structural abnormalities in the muscle fibres, i.e. nemaline bodies. The related disorder cap myopathy is defined by cap-like structures located peripherally in the muscle fibres. Both disorders may be caused by mutations in the TPM2 gene encoding ß-Tm (tropomyosin). Tm controls muscle contraction by inhibiting actin-myosin interaction in a calcium-sensitive manner. In the present study, we have investigated the pathogenetic mechanisms underlying five disease-causing mutations in Tm. We show that four of the mutations cause changes in affinity for actin, which may cause muscle weakness in these patients, whereas two show defective Ca2+ activation of contractility. We have also mapped the amino acids altered by the mutation to regions important for actin binding and note that two of the mutations cause altered protein conformation, which could account for impaired actin affinity.

Pharmacoepigenetics of hypertension: genome-wide methylation analysis of responsiveness to four classes of antihypertensive drugs using a double-blind crossover study design.

Essential hypertension remains the leading risk factor of global disease burden, but its treatment goals are often not met. We investigated whether DNA methylation is associated with antihypertensive responses to a diuretic, a beta-blocker, a calcium channel blocker or an angiotensin receptor antagonist. In addition, since we previously showed an SNP at the transcription start site (TSS) of the catecholamine biosynthesis-related ACY3 gene to associate with blood pressure (BP) response to beta-blockers, we specifically analysed the association of methylation sites close to the ACY3 TSS with BP responses to beta-blockers. We conducted an epigenome-wide association study between leukocyte DNA methylation and BP responses to antihypertensive monotherapies in two hypertensive Finnish cohorts: the GENRES (https://clinicaltrials.gov/ct2/show/NCT03276598; amlodipine 5 mg, bisoprolol 5 mg, hydrochlorothiazide 25 mg, or losartan 50 mg daily) and the LIFE-Fin studies (https://clinicaltrials.gov/ct2/show/NCT00338260; atenolol 50 mg or losartan 50 mg daily). The monotherapy groups consisted of approximately 200 individuals each. We identified 64 methylation sites to suggestively associate (P < 1E-5) with either systolic or diastolic BP responses to a particular study drug in GENRES. These associations did not replicate in LIFE-Fin . Three methylation sites close to the ACY3 TSS were associated with systolic BP responses to bisoprolol in GENRES but not genome-wide significantly (P < 0.05). No robust associations between DNA methylation and BP responses to four different antihypertensive drugs were identified. However, the findings on the methylation sites close to the ACY3 TSS may support the role of ACY3 genetic and epigenetic variation in BP response to bisoprolol.

Dominant Distal Myopathy 3 (MPD3) Caused by a Deletion in the HNRNPA1 Gene.

BACKGROUND AND OBJECTIVES: To determine the genetic cause of the disease in the previously reported family with adult-onset autosomal dominant distal myopathy (myopathy, distal, 3; MPD3). METHODS: Continued clinical evaluation including muscle MRI and muscle pathology. A linkage analysis with single nucleotide polymorphism arrays and genome sequencing were used to identify the genetic defect, which was verified by Sanger sequencing. RNA sequencing was used to investigate the transcriptional effects of the identified genetic defect. RESULTS: Small hand muscles (intrinsic, thenar, and hypothenar) were first involved with spread to the lower legs and later proximal muscles. Dystrophic changes with rimmed vacuoles and cytoplasmic inclusions were observed in muscle biopsies at advanced stage. A single nucleotide polymorphism array confirmed the previous microsatellite-based linkage to 8p22-q11 and 12q13-q22. Genome sequencing of three affected family members combined with structural variant calling revealed a small heterozygous deletion of 160 base pairs spanning the second last exon 10 of the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1) gene, which is in the linked region on chromosome 12. Segregation of the mutation with the disease was confirmed by Sanger sequencing. RNA sequencing showed that the mutant allele produces a shorter mutant mRNA transcript compared with the wild-type allele. Immunofluorescence studies on muscle biopsies revealed small p62 and larger TDP-43 inclusions. DISCUSSION: A small exon 10 deletion in the gene HNRNPA1 was identified as the cause of MPD3 in this family. The new HNRNPA1-related phenotype, upper limb presenting distal myopathy, was thus confirmed, and the family displays the complexities of gene identification.

Common genetic variation of beta1- and beta2-adrenergic receptor and response to four classes of antihypertensive treatment.

Varying results have been reported on the association of beta-adrenergic receptor polymorphisms with blood pressure (BP) response to beta-blockers. We investigated the influence of ADRB1 Ser49Gly and Arg389Gly, and ADRB2 Gly16Arg and Glu27Gln polymorphisms on ambulatory BP response to bisoprolol and three other antihypertensive drug monotherapies in a placebo-controlled, double-blind, cross-over study with 233 moderately hypertensive men. ADRB1 Ser49Ser homozygotes tended to have a better ambulatory BP response to bisoprolol but the difference was statistically nonsignificant. ADRB1 Arg389Arg homozygotes did not show better BP response to bisoprolol than the other genotypes. There were no significant associations of ADRB2 polymorphisms with BP responses to any of the study drugs. The results from this controlled study in hypertensive men do not support clinical use of common polymorphisms in ADRB1 and ADRB2 in predicting BP responses to beta-blockers or to three other antihypertensive drugs.

Genomics of asthma, allergy and chronic rhinosinusitis: novel concepts and relevance in airway mucosa.

Genome wide association studies (GWASs) have revealed several airway disease-associated risk loci. Their role in the onset of asthma, allergic rhinitis (AR) or chronic rhinosinusitis (CRS), however, is not yet fully understood. The aim of this review is to evaluate the airway relevance of loci and genes identified in GWAS studies. GWASs were searched from databases, and a list of loci associating significantly (p < 10-8) with asthma, AR and CRS was created. This yielded a total of 267 significantly asthma/AR-associated loci from 31 GWASs. No significant CRS -associated loci were found in this search. A total of 170 protein coding genes were connected to these loci. Of these, 76/170 (44%) showed bronchial epithelial protein expression in stained microscopic figures of Human Protein Atlas (HPA), and 61/170 (36%) had a literature report of having airway epithelial function. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analyses were performed, and 19 functional protein categories were found as significantly (p < 0.05) enriched among these genes. These were related to cytokine production, cell activation and adaptive immune response, and all were strongly connected in network analysis. We also identified 15 protein pathways that were significantly (p < 0.05) enriched in these genes, related to T-helper cell differentiation, virus infection, JAK-STAT signaling pathway, and asthma. A third of GWAS-level risk loci genes of asthma or AR seemed to have airway epithelial functions according to our database and literature searches. In addition, many of the risk loci genes were immunity related. Some risk loci genes also related to metabolism, neuro-musculoskeletal or other functions. Functions overlapped and formed a strong network in our pathway analyses and are worth future studies of biomarker and therapeutics.

Heme oxygenase-1 repeat polymorphism in septic acute kidney injury.

Acute kidney injury (AKI) is a syndrome that frequently affects the critically ill. Recently, an increased number of dinucleotide repeats in the HMOX1 gene were reported to associate with development of AKI in cardiac surgery. We aimed to test the replicability of this finding in a Finnish cohort of critically ill septic patients. This multicenter study was part of the national FINNAKI study. We genotyped 300 patients with severe AKI (KDIGO 2 or 3) and 353 controls without AKI (KDIGO 0) for the guanine-thymine (GTn) repeat in the promoter region of the HMOX1 gene. The allele calling was based on the number of repeats, the cut off being 27 repeats in the S-L (short to long) classification, and 27 and 34 repeats for the S-M-L2 (short to medium to very long) classification. The plasma concentrations of heme oxygenase-1 (HO-1) enzyme were measured on admission. The allele distribution in our patients was similar to that published previously, with peaks at 23 and 30 repeats. The S-allele increases AKI risk. An adjusted OR was 1.30 for each S-allele in an additive genetic model (95% CI 1.01-1.66; p = 0.041). Alleles with a repeat number greater than 34 were significantly associated with lower HO-1 concentration (p<0.001). In septic patients, we report an association between a short repeat in HMOX1 and AKI risk.

Identification of a founder mutation in TPM3 in nemaline myopathy patients of Turkish origin.

To date, six genes are known to cause nemaline (rod) myopathy (NM), a rare congenital neuromuscular disorder. In an attempt to find a seventh gene, we performed linkage and subsequent sequence analyses in 12 Turkish families with recessive NM. We found homozygosity in two of the families at 1q12-21.2, a region encompassing the gamma-tropomyosin gene (TPM3) encoding slow skeletal muscle alpha-tropomyosin, a known NM gene. Sequencing revealed homozygous deletion of the first nucleotide of the last exon, c.913delA of TPM3 in both families. The mutation removes the last nucleotide before the stop codon, causing a frameshift and readthrough across the termination signal. The encoded alphaTm(slow) protein is predicted to be 73 amino acids longer than normal, and the extension to the protein is hypothesised to be unable to form a coiled coil. The resulting tropomyosin protein may therefore be non-functional. The affected children in both families were homozygous for the mutation, while the healthy parents were mutation carriers. Both of the patients in Family 1 had the severe form of NM, and also an unusual chest deformity. The affected children in Family 2 had the intermediate form of NM. Muscle biopsies showed type 1 (slow) fibres to be markedly smaller than type 2 (fast) fibres. Previously, there had been five reports, only, of NM caused by mutations in TPM3. The mutation reported here is the first deletion to be identified in TPM3, and it is likely to be a founder mutation in the Turkish population.

Laboratory tests as predictors of the antihypertensive effects of amlodipine, bisoprolol, hydrochlorothiazide and losartan in men: results from the randomized, double-blind, crossover GENRES Study.

OBJECTIVE: Individual blood pressure responses to antihypertensive therapy are difficult to predict. To improve optimization of antihypertensive therapy, we analyzed correlations of relevant laboratory tests with blood pressure responses to four antihypertensive monotherapies. METHODS: In the GENRES study, 208 Finnish men aged 35-60 years with moderate hypertension used amlodipine 5 mg, bisoprolol 5 mg, hydrochlorothiazide 25 mg and losartan 50 mg daily, each for 4 weeks as a monotherapy in a double-blind, randomized, placebo-controlled crossover study; that is, each subject received each type of monotherapy in a random order. The treatment periods were preceded and separated by 4-week placebo periods. Ambulatory 24-h and office blood pressure measurements were carried out after all study periods. Data from several biochemical tests were correlated to antihypertensive drug responses. RESULTS: Serum total calcium concentration was negatively correlated with blood pressure responses to amlodipine (P values 0.001-0.002). Plasma renin activity was positively correlated with blood pressure responses to losartan (P values 0.001-0.005) and bisoprolol (P values 0.03-0.17), and negatively with blood pressure responses to hydrochlorothiazide (P values 0.01-0.07). Daily urinary excretion of sodium was negatively correlated with ambulatory blood pressure responses to amlodipine (P values 0.001-0.01). CONCLUSIONS: In this carefully controlled study, marked individual variations in antihypertensive drug responsiveness were found to correlate to several baseline laboratory parameters. The negative correlation between serum calcium levels and amlodipine responses is intriguing and suggests an underlying mechanistic association. Collectively, our data imply that laboratory tests may have some value in prediction of the efficacy of various antihypertensive drug therapies, although great patient-to-patient variation remains an obstacle for exact predictive classification.

PRKCDBP (CAVIN3) and CRY2 associate with major depressive disorder.

BACKGROUND: Dysfunctions in the intrinsic clocks are suggested in patients with depressive disorders. The cryptochrome circadian clocks 1 and 2 (CRY1 and CRY2) proteins modulate circadian rhythms in a cell and influence emotional reactions and mood in an individual. The protein kinase C delta binding protein (PRKCDBP, or CAVIN3), similar to the serum deprivation response protein (SDPR, or CAVIN2), reduces metabolic stability of the PER2-CRY2 transcription factor complex that plays a role in the circadian rhythm synchronization. Our aim was to study SDPR, PRKCDBP, CRY1 and CRY2 genetic variants in depressive disorders. METHODS: The sample included 5910 Finnish individuals assessed with the Munich-Composite International Diagnostic Interview (M-CIDI) in year 2000. In year 2011, 3424 individuals were assessed again. After genotype quality control, there were 383 subjects with major depressive disorder, 166 with dysthymia, and 479 with depressive disorders (major depressive disorder, dysthymia or both), and 4154 healthy controls. A total of 48 single-nucleotide polymorphisms from SDPR, PRKCDBP, CRY1 and CRY2 genes were analyzed using logistic regression models controlling for age and gender. RESULTS: The earlier reported association of CRY2 variants with dysthymia was confirmed and extended to major depressive disorder (q<0.05). In addition, novel associations of PRKCDBP rs1488864 with depressive disorders (q=0.02) and with major depressive disorder in specific (q=0.007) were found. LIMITATIONS: The number of cases was moderate and coverage of PRKCDB was limited. CONCLUSIONS: CRY2 and PRKCDBP variants may be risk factors of major depressive disorder and provide information for diagnosis.

Replicated evidence for aminoacylase 3 and nephrin gene variations to predict antihypertensive drug responses.

AIM: To replicate the genome-wide associations of the antihypertensive effects of bisoprolol and losartan in GENRES, using the Finnish patients of LIFE study. PATIENTS & METHODS: We analyzed association of four SNPs with atenolol and three SNPs with losartan response in 927 Finnish LIFE patients (467 for atenolol and 460 for losartan). RESULTS: rs2514036, a variation at a transcription start site of ACY3, was associated with blood pressure response to atenolol in men in LIFE. Response to bisoprolol was correlated to baseline plasma levels of N-acetylphenylalanine and phenylalanine (ACY3 substrate and end product, respectively) in GENRES study. NPHS1 variation rs3814995 was associated with losartan effect in LIFE. CONCLUSION: We provide support for two pharmacogenomic markers for beta-blockers and angiotensin receptor antagonists.

Identification of 45 novel mutations in the nebulin gene associated with autosomal recessive nemaline myopathy.

Nemaline myopathy (NM) is a clinically and genetically heterogeneous disorder of skeletal muscle caused by mutations in at least five different genes encoding thin filament proteins of the striated muscle sarcomere. We have previously described 18 different mutations in the last 42 exons of the nebulin gene (NEB) in 18 families with NM. Here we report 45 novel NEB mutations detected by denaturing high-performance liquid chromatography (dHPLC) and sequence analysis of all 183 NEB exons in NM patients from 44 families. Altogether we have identified, including the deletion of exon 55 identified in the Ashkenazi Jewish population, 64 different mutations in NEB segregating with autosomal recessive NM in 55 families. The majority (55%) of the mutations in NEB are frameshift or nonsense mutations predicted to cause premature truncation of nebulin. Point mutations (25%) or deletions (3%) affecting conserved splice signals are predicted in the majority of cases to cause in-frame exon skipping, possibly leading to impaired nebulin-tropomyosin interaction along the thin filament. Patients in 18 families had one of nine missense mutations (14%) affecting conserved amino acids at or in the vicinity of actin or tropomyosin binding sites. In addition, we found the exon 55 deletion in four families. The majority of the patients (in 49/55 families) were shown to be compound heterozygous for two different mutations. The mutations were found in both constitutively and alternatively expressed exons throughout the NEB gene, and there were no obvious mutational hotspots. Patients with more severe clinical pictures tended to have mutations predicted to be more disruptive than patients with milder forms.

CRY1 and CRY2 genetic variants in seasonality: A longitudinal and cross-sectional study.

Cryptochromes are key components of the circadian clocks that generate and maintain seasonal variations. The aim of our study was to analyze the associations of CRY1 and CRY2 genetic variants with the problematicity of seasonal variations, and whether the problematicity of seasonal variations changed during the follow-up of 11 years. Altogether 21 CRY1 and 16 CRY2 single-nucleotide polymorphisms (SNPs) were genotyped and analyzed in 5910 individuals from a Finnish nationwide population-based sample who had filled in the self-report on the seasonal variations in mood and behavior in the year 2000. In the year 2011, 3356 of these individuals filled in the same self-report on the seasonal variations in mood and behavior. Regression models were used to test whether any of the SNPs associated with the problematicity of seasonal variations or with a change in the problematicity from 2000 to 2011. In the longitudinal analysis, CRY2 SNP rs61884508 was protective from worsening of problematicity of seasonal variations. In the cross-sectional analysis, CRY2 SNP rs72902437 showed evidence of association with problematicity of seasonal variations, as did SNP rs1554338 (in the MAPK8IP1 and downstream of CRY2).

SIRT1 Polymorphisms Associate with Seasonal Weight Variation, Depressive Disorders, and Diastolic Blood Pressure in the General Population.

SIRT1 polymorphisms have previously been associated with depressive and anxiety disorders. We aimed at confirming these earlier findings and extending the analyses to seasonal variations in mood and behavior. Three tag single-nucleotide polymorphisms (SNPs) were selected to capture the common variation in the SIRT1 gene. 5910 individuals (with blood sample, diagnostic interview, self-report of on seasonal changes in mood and behavior) were selected from a representative Finnish nationwide population-based sample. Logistic and linear regression models were used to analyze the associations between the SNPs and depressive and anxiety disorders, metabolic syndrome (EGIR criteria) and its components, and health examination measurements, Homeostasis Model Assessments, and diagnoses of type 2 and type 1 diabetes. SIRT1 rs2273773 showed evidence of association with seasonal variation in weight (C-allele, OR = 0.85, 95% CI = 0.76-0.95, p = 0.005). In addition, our study gave further support for the association of SIRT1 gene with depressive disorders (rs3758391) and diastolic blood pressure (rs2273773).

CRY1, CRY2 and PRKCDBP genetic variants in metabolic syndrome.

The circadian clock affects metabolic cycles, and there is a link between circadian clock genes and metabolic syndrome. Therefore, we wanted to investigate whether variants of the core circadian clock genes, cryptochrome circadian clocks 1 and 2 (CRY1 and CRY2), or those of protein kinase C, delta binding protein (PRKCDBP), which regulate the interactions and abundance of dimers of the period and cryptochrome proteins, are associated with metabolic syndrome or its components. The association of 48 single-nucleotide polymorphisms (SNPs) from CRY1, CRY2 and PRKCDBP genes with metabolic disorder or its components was analyzed in a sample of 5910 individuals. Genotyping was performed using the Sequenom MassARRAY system. SNPs and haplotypes were analyzed using linear or logistic regression with additive models controlling for age and sex. Continuous phenotypes were permuted 10,000 times. False discovery rate q-values were calculated to correct for multiple testing. Overall, CRY1 and CRY2 variants showed nominal association with the metabolic syndrome components, hypertension and triglyceride levels, and one CRY2 variant had an association with metabolic syndrome, although none of these associations yielded significant q-values. However, the haplotype analysis of these variants supported the association of CRY1 with arterial hypertension and elevated blood pressure. Further studies are warranted regarding the role of CRY1 in cardiovascular diseases.

Genetic Loci Associated with Allergic Sensitization in Lithuanians.

Allergic rhinitis (AR) is a common and complex disease. It is associated with environmental as well as genetic factors. Three recent genome-wide association studies (GWAS) reported altogether 47 single nucleotide polymorphisms (SNPs) associated with AR or allergic sensitization (AS) in Europeans and North Americans. Two follow up studies in Swedish and Chinese replicated 15 associations. In these studies individuals were selected based on the self-reported AR, or AR/AS diagnosed using blood IgE test or skin prick test (SPT), which were performed often without restriction to specific allergens. Here we performed third replication study in Lithuanians. We used SPT and carefully selected set of allergens prevalent in Lithuania, as well as Illumina Core Exome chip for SNP detection. We genotyped 270 SPT-positive individuals (137 Betulaceae -, 174 Poaceae-, 199 Artemisia-, 70 Helianthus-, 22 Alternaria-, 22 Cladosporium-, 140 mites-, 95 cat- and 97 dog dander-sensitive cases) and 162 SPT-negative controls. We found altogether 13 known SNPs associated with AS (p ≤0.05). Three SNPs were found in Lithuanians sensitive to several allergens, and 10 SNPs were found in Lithuanians sensitive to a certain allergen. For the first time, SNP rs7775228:C was associated with patient sensitivity to dog allergens (F_A=0,269, F_U=0.180, P=0.008). Thus, careful assessment of AS allowed us to detect known genetic variants associated with AS/AR in relatively small cohort of Lithuanians.