Immunogenetics associated with severe coccidioidomycosis.

Disseminated coccidioidomycosis (DCM) is caused by Coccidioides, pathogenic fungi endemic to the southwestern United States and Mexico. Illness occurs in approximately 30% of those infected, less than 1% of whom develop disseminated disease. To address why some individuals allow dissemination, we enrolled patients with DCM and performed whole-exome sequencing. In an exploratory set of 67 patients with DCM, 2 had haploinsufficient STAT3 mutations, and defects in ß-glucan sensing and response were seen in 34 of 67 cases. Damaging CLEC7A and PLCG2 variants were associated with impaired production of ß-glucan-stimulated TNF-α from PBMCs compared with healthy controls. Using ancestry-matched controls, damaging CLEC7A and PLCG2 variants were overrepresented in DCM, including CLEC7A Y238* and PLCG2 R268W. A validation cohort of 111 patients with DCM confirmed the PLCG2 R268W, CLEC7A I223S, and CLEC7A Y238* variants. Stimulation with a DECTIN-1 agonist induced DUOX1/DUOXA1-derived hydrogen peroxide [H2O2] in transfected cells. Heterozygous DUOX1 or DUOXA1 variants that impaired H2O2 production were overrepresented in discovery and validation cohorts. Patients with DCM have impaired ß-glucan sensing or response affecting TNF-α and H2O2 production. Impaired Coccidioides recognition and decreased cellular response are associated with disseminated coccidioidomycosis.

Epithelial cell responses to rhinovirus identify an early-life-onset asthma phenotype in adults.

BACKGROUND: The study of pathogenic mechanisms in adult asthma is often marred by a lack of precise information about the natural history of the disease. Children who have persistent wheezing (PW) during the first 6 years of life and whose symptoms start before age 3 years (PW+) are much more likely to have wheezing illnesses due to rhinovirus (RV) in infancy and to have asthma into adult life than are those who do not have PW (PW-). OBJECTIVE: Our aim was to determine whether nasal epithelial cells from PW+ asthmatic adults as compared with cells from PW- asthmatic adults show distinct biomechanistic processes activated by RV exposure. METHODS: Air-liquid interface cultures derived from nasal epithelial cells of 36-year old participants with active asthma with and without a history of PW in childhood (10 PW+ participants and 20 PW- participants) from the Tucson Children's Respiratory Study were challenged with a human RV-A strain (RV-A16) or control, and their RNA was sequenced. RESULTS: A total of 35 differentially expressed genes involved in extracellular remodeling and angiogenesis distinguished the PW+ group from the PW- group at baseline and after RV-A stimulation. Notably, 22 transcriptomic pathways showed PW-by-RV interactions; the pathways were invariably overactivated in PW+ patients, and were involved in Toll-like receptor- and cytokine-mediated responses, remodeling, and angiogenic processes. CONCLUSIONS: Asthmatic adults with a history of persistent wheeze in the first 6 years of life have specific biomolecular alterations in response to RV-A that are not present in patients without such a history. Targeting these mechanisms may slow the progression of asthma in these patients.

'Single-subject studies'-derived analyses unveil altered biomechanisms between very small cohorts: implications for rare diseases.

MOTIVATION: Identifying altered transcripts between very small human cohorts is particularly challenging and is compounded by the low accrual rate of human subjects in rare diseases or sub-stratified common disorders. Yet, single-subject studies (S3) can compare paired transcriptome samples drawn from the same patient under two conditions (e.g. treated versus pre-treatment) and suggest patient-specific responsive biomechanisms based on the overrepresentation of functionally defined gene sets. These improve statistical power by: (i) reducing the total features tested and (ii) relaxing the requirement of within-cohort uniformity at the transcript level. We propose Inter-N-of-1, a novel method, to identify meaningful differences between very small cohorts by using the effect size of 'single-subject-study'-derived responsive biological mechanisms. RESULTS: In each subject, Inter-N-of-1 requires applying previously published S3-type N-of-1-pathways MixEnrich to two paired samples (e.g. diseased versus unaffected tissues) for determining patient-specific enriched genes sets: Odds Ratios (S3-OR) and S3-variance using Gene Ontology Biological Processes. To evaluate small cohorts, we calculated the precision and recall of Inter-N-of-1 and that of a control method (GLM+EGS) when comparing two cohorts of decreasing sizes (from 20 versus 20 to 2 versus 2) in a comprehensive six-parameter simulation and in a proof-of-concept clinical dataset. In simulations, the Inter-N-of-1 median precision and recall are > 90% and >75% in cohorts of 3 versus 3 distinct subjects (regardless of the parameter values), whereas conventional methods outperform Inter-N-of-1 at sample sizes 9 versus 9 and larger. Similar results were obtained in the clinical proof-of-concept dataset. AVAILABILITY AND IMPLEMENTATION: R software is available at Lussierlab.net/BSSD.

Genome-Wide Association Study in African Americans with Acute Respiratory Distress Syndrome Identifies the Selectin P Ligand Gene as a Risk Factor.

RATIONALE: Genetic factors are involved in acute respiratory distress syndrome (ARDS) susceptibility. Identification of novel candidate genes associated with increased risk and severity will improve our understanding of ARDS pathophysiology and enhance efforts to develop novel preventive and therapeutic approaches. OBJECTIVES: To identify genetic susceptibility targets for ARDS. METHODS: A genome-wide association study was performed on 232 African American patients with ARDS and 162 at-risk control subjects. The Identify Candidate Causal SNPs and Pathways platform was used to infer the association of known gene sets with the top prioritized intragenic SNPs. Preclinical validation of SELPLG (selectin P ligand gene) was performed using mouse models of LPS- and ventilator-induced lung injury. Exonic variation within SELPLG distinguishing patients with ARDS from sepsis control subjects was confirmed in an independent cohort. MEASUREMENTS AND MAIN RESULTS: Pathway prioritization analysis identified a nonsynonymous coding SNP (rs2228315) within SELPLG, encoding P-selectin glycoprotein ligand 1, to be associated with increased susceptibility. In an independent cohort, two exonic SELPLG SNPs were significantly associated with ARDS susceptibility. Additional support for SELPLG as an ARDS candidate gene was derived from preclinical ARDS models where SELPLG gene expression in lung tissues was significantly increased in both ventilator-induced (twofold increase) and LPS-induced (5.7-fold increase) murine lung injury models compared with controls. Furthermore, Selplg-/- mice exhibited significantly reduced LPS-induced inflammatory lung injury compared with wild-type C57/B6 mice. Finally, an antibody that neutralizes P-selectin glycoprotein ligand 1 significantly attenuated LPS-induced lung inflammation. CONCLUSIONS: These findings identify SELPLG as a novel ARDS susceptibility gene among individuals of European and African descent.

Novel disease syndromes unveiled by integrative multiscale network analysis of diseases sharing molecular effectors and comorbidities.

BACKGROUND: Forty-two percent of patients experience disease comorbidity, contributing substantially to mortality rates and increased healthcare costs. Yet, the possibility of underlying shared mechanisms for diseases remains not well established, and few studies have confirmed their molecular predictions with clinical datasets. METHODS: In this work, we integrated genome-wide association study (GWAS) associating diseases and single nucleotide polymorphisms (SNPs) with transcript regulatory activity from expression quantitative trait loci (eQTL). This allowed novel mechanistic insights for noncoding and intergenic regions. We then analyzed pairs of SNPs across diseases to identify shared molecular effectors robust to multiple test correction (False Discovery Rate FDReRNA < 0.05). We hypothesized that disease pairs found to be molecularly convergent would also be significantly overrepresented among comorbidities in clinical datasets. To assess our hypothesis, we used clinical claims datasets from the Healthcare Cost and Utilization Project (HCUP) and calculated significant disease comorbidities (FDRcomorbidity < 0.05). We finally verified if disease pairs resulting molecularly convergent were also statistically comorbid more than by chance using the Fisher's Exact Test. RESULTS: Our approach integrates: (i) 6175 SNPs associated with 238 diseases from ~ 1000 GWAS, (ii) eQTL associations from 19 tissues, and (iii) claims data for 35 million patients from HCUP. Logistic regression (controlled for age, gender, and race) identified comorbidities in HCUP, while enrichment analyses identified cis- and trans-eQTL downstream effectors of GWAS-identified variants. Among ~ 16,000 combinations of diseases, 398 disease-pairs were prioritized by both convergent eQTL-genetics (RNA overlap enrichment, FDReRNA < 0.05) and clinical comorbidities (OR > 1.5, FDRcomorbidity < 0.05). Case studies of comorbidities illustrate specific convergent noncoding regulatory elements. An intergenic architecture of disease comorbidity was unveiled due to GWAS and eQTL-derived convergent mechanisms between distinct diseases being overrepresented among observed comorbidities in clinical datasets (OR = 8.6, p-value = 6.4 × 10- 5 FET). CONCLUSIONS: These comorbid diseases with convergent eQTL genetic mechanisms suggest clinical syndromes. While it took over a decade to confirm the genetic underpinning of the metabolic syndrome, this study is likely highlighting hundreds of new ones. Further, this knowledge may improve the clinical management of comorbidities with precision and shed light on novel approaches of drug repositioning or SNP-guided precision molecular therapy inclusive of intergenic risks.

Single subject transcriptome analysis to identify functionally signed gene set or pathway activity.

Analysis of single-subject transcriptome response data is an unmet need of precision medicine, made challenging by the high dimension, dynamic nature and difficulty in extracting meaningful signals from biological or stochastic noise. We have proposed a method for single subject analysis that uses a mixture model for transcript fold-change clustering from isogenically paired samples, followed by integration of these distributions with Gene Ontology Biological Processes (GO-BP) to reduce dimension and identify functional attributes. We then extended these methods to develop functional signing metrics for gene set process regulation by incorporating biological repressor relationships encoded in GO-BP as negatively_regulates edges. Results revealed reproducible and biologically meaningful signals from analysis of a single subject's response, opening the door to future transcriptomic studies where subject and resource availability are currently limiting. We used inbred mouse strains fed different diets to provide isogenic biological replicates, permitting rigorous validation of our method. We compared significant genotype-specific GO-BP term results for overlap and rank order across three replicate pairs per genotype, and cross-methods to reference standards (limma+FET, SAM+FET, and GSEA). All single-subject analytics findings were robust and highly reproducible (median area under the ROC curve=0.96, n=24 genotypes × 3 replicates), providing confidence and validation of this approach for analyses in single subjects. R code is available online at http://www.lussiergroup.org/publications/PathwayActivity.

Identification of Jak-STAT signaling involvement in sarcoidosis severity via a novel microRNA-regulated peripheral blood mononuclear cell gene signature.

Sarcoidosis is a granulomatous lung disorder of unknown cause. The majority of individuals with sarcoidosis spontaneously achieve full remission (uncomplicated sarcoidosis), however, ~20% of sarcoidosis-affected individuals experience progressive lung disease or cardiac and nervous system involvement (complicated sarcoidosis). We investigated peripheral blood mononuclear cell (PBMC) microRNA and protein-coding gene expression data from healthy controls and patients with uncomplicated or complicated sarcoidosis. We identified 46 microRNAs and 1,559 genes that were differentially expressed across a continuum of sarcoidosis severity (healthy control → uncomplicated sarcoidosis → complicated sarcoidosis). A total of 19 microRNA-mRNA regulatory pairs were identified within these deregulated microRNAs and mRNAs, which consisted of 17 unique protein-coding genes yielding a 17-gene signature. Pathway analysis of the 17-gene signature revealed Jak-STAT signaling pathway as the most significantly represented pathway. A severity score was assigned to each patient based on the expression of the 17-gene signature and a significant increasing trend in the severity score was observed from healthy control, to uncomplicated sarcoidosis, and finally to complicated sarcoidosis. In addition, this microRNA-regulated gene signature differentiates sarcoidosis patients from healthy controls in independent validation cohorts. Our study suggests that PBMC gene expression is useful in diagnosis of sarcoidosis.

Complex genetics of pulmonary diseases: lessons from genome-wide association studies and next-generation sequencing.

The advent of high-throughput technologies has provided exceptional assistance for lung scientists to discover novel genetic variants underlying the development and progression of complex lung diseases. However, the discovered variants thus far do not explain much of the estimated heritability of complex lung diseases. Here, we review the literature of successfully used genome-wide association studies (GWASs) and identified the polymorphisms that reproducibly underpin the susceptibility to various noncancerous complex lung diseases or affect therapeutic responses. We also discuss the inherent limitations of GWAS approaches and how the use of next-generation sequencing technologies has furthered our understanding about the genetic determinants of these diseases. Next, we describe the contribution of the metagenomics to understand the interactions of the airways microbiome with lung diseases. We then highlight the urgent need for new integrative genomics-phenomics methods to more effectively interrogate and understand multiple downstream "omics" (eg, chromatin modification patterns). Finally, we address the scarcity of genetic studies addressing under-represented populations such as African Americans and Hispanics.

eQTL networks unveil enriched mRNA master integrators downstream of complex disease-associated SNPs.

The causal and interplay mechanisms of Single Nucleotide Polymorphisms (SNPs) associated with complex diseases (complex disease SNPs) investigated in genome-wide association studies (GWAS) at the transcriptional level (mRNA) are poorly understood despite recent advancements such as discoveries reported in the Encyclopedia of DNA Elements (ENCODE) and Genotype-Tissue Expression (GTex). Protein interaction network analyses have successfully improved our understanding of both single gene diseases (Mendelian diseases) and complex diseases. Whether the mRNAs downstream of complex disease genes are central or peripheral in the genetic information flow relating DNA to mRNA remains unclear and may be disease-specific. Using expression Quantitative Trait Loci (eQTL) that provide DNA to mRNA associations and network centrality metrics, we hypothesize that we can unveil the systems properties of information flow between SNPs and the transcriptomes of complex diseases. We compare different conditions such as naïve SNP assignments and stringent linkage disequilibrium (LD) free assignments for transcripts to remove confounders from LD. Additionally, we compare the results from eQTL networks between lymphoblastoid cell lines and liver tissue. Empirical permutation resampling (p<0.001) and theoretic Mann-Whitney U test (p<10(-30)) statistics indicate that mRNAs corresponding to complex disease SNPs via eQTL associations are likely to be regulated by a larger number of SNPs than expected. We name this novel property mRNA hubness in eQTL networks, and further term mRNAs with high hubness as master integrators. mRNA master integrators receive and coordinate the perturbation signals from large numbers of polymorphisms and respond to the personal genetic architecture integratively. This genetic signal integration contrasts with the mechanism underlying some Mendelian diseases, where a genetic polymorphism affecting a single protein hub produces a divergent signal that affects a large number of downstream proteins. Indeed, we verify that this property is independent of the hubness in protein networks for which these mRNAs are transcribed. Our findings provide novel insights into the pleiotropy of mRNAs targeted by complex disease polymorphisms and the architecture of the information flow between the genetic polymorphisms and transcriptomes of complex diseases.

Accelerating precision biology and medicine with computational biology and bioinformatics.

A report on the 22nd Annual International Conference on Intelligent Systems for Molecular Biology, held in Boston, Massachusetts, USA, July 11-15, 2014.

Combining functional genomics strategies identifies modular heterogeneity of breast cancer intrinsic subtypes.

BACKGROUND: The discovery of breast cancer subtypes and subsequent development of treatments aimed at them has allowed for a great reduction in the mortality of breast cancer. But despite this progress, tumors with similar characteristics that belong to the same subtype continue to respond differently to the same treatment. Five subtypes of breast cancer, namely intrinsic subtypes, have been characterized to date based on their gene expression profiles. Among other characteristics, subtypes vary in their degree of intra-subtype heterogeneity. It is not clear, however, whether this heterogeneity is shared across all tumor traits. It is also unclear whether individual traits can be highly heterogeneous among a majority of homogeneous traits. RESULTS: We employ network theory to uncover gene modules and accordingly consider them as tumor traits, which capture shared biological processes among the subtypes. We then use the ß-diversity metric from ecology to quantify the heterogeneity in these gene modules. In doing so, we show that breast cancer heterogeneity is contained in gene modules and that this modular heterogeneity increases monotonically across the subtypes. We identify a core of two modules that are shared among all subtypes which contain nucleosome assembly and mammary morphogenesis genes, and a number of modules that are specific to subtypes. This modular heterogeneity, which increases with global heterogeneity, relates to tumor aggressiveness. Indeed, we observe that Luminal A, the most treatable of subtypes, has the lowest modular heterogeneity whereas the Basal-like subtype, which is among the hardest to treat, has the highest. Furthermore, our analysis shows that a higher degree of global heterogeneity does not imply higher heterogeneity for all modules, as Luminal B shows the highest heterogeneity for core modules. CONCLUSIONS: Overall, modular heterogeneity provides a framework with which to dissect cancer heterogeneity and better understand its underpinnings, thereby ultimately advancing our knowledge towards a more effective personalized cancer therapy.

Evaluation of liver diseases in Iranian patients with primary antibody deficiencies.

INTRODUCTION: Patients with primary antibody deficiency (PAD) can complicate with liver disease. This study was performed in order to study the prevalence and causes of hepatobiliary diseases in Iranian patients with PAD. MATERIAL AND METHODS: Sixty-two patients with PAD were followed-up and signs and symptoms of liver disease were recorded. All patients were screened for hepatitis C virus (HCV-RNA) and those patients with any sign of liver disease or gastrointestinal complaints were tested for Cryptosporidium parvum. RESULTS: Clinical evidences of liver disease, including hepatomegaly, were documented in 22 patients (35.5%). Eight patients (13%) had clinical and/or laboratory criteria of chronic liver disease. Only one patient was HCV-RNA positive; he had stigmata of chronic liver disease and pathologic evidence of chronic active hepatitis with cirrhosis. Cryptosporidium parvum test was positive for one patient with hyper-IgM syndrome. In liver biopsy of patients with liver involvement, one had histological findings related to sclerosing cholangitis, and five had mild to moderate chronic active hepatitis with unknown reason. CONCLUSIONS: Chronic active hepatitis is the most common pathologic feature of liver injury in Iranian patients with PAD. Liver disease in PAD usually accompanies with other organ involvements and could increase the mortality of PAD. Whether this high rate of liver disease with unknown origin (75%) is the result of an unidentified hepatotropic virus or other mechanisms such as autoimmunity, is currently difficult to understand.

Presence of Idiopathic Thrombocytopenic Purpura and autoimmune hemolytic anemia in the patients with common variable immunodeficiency.

Common Variable Immunodeficiency (CVID) is a heterogeneous group of disorders characterized by hypogammaglobulinemia and an increased susceptibility to recurrent infections as well as autoimmunity and malignancies. Idiopathic Thrombocytopenic Purpura (ITP) and Autoimmune Hemolytic Anemia (AIHA) are two autoimmune disorders which may be seen in association with CVID. Among 85 CVID patients, seven cases had ITP and/or AIHA (8%). Four of these patients had one or more episodes of ITP, one patient had AIHA, and two patients had both ITP and AIHA (Evans syndrome). Almost, all patients experienced chronic and recurrent infections mostly in respiratory and gastrointestinal systems during the course of the disease. Among the seven patients, five presented their underlying disease with recurrent respiratory and/or gastrointestinal tract infections, while in two remaining patients, CVID was presented with ITP. Three patients died until now; two because of hepatic failure and one due to pulmonary hemorrhage. As CVID is prone to autoimmune disorders, it should be considered as a differential diagnosis of adult-onset ITP and possibly in children. Chronic and recurrent ITP, especially in the presence of propensity to respiratory and gastrointestinal infections mandate the evaluation for an underlying immune dysregulation such as CVID.

Genotype-phenotype correlation in Bruton's tyrosine kinase deficiency.

Bruton's tyrosine kinase (Btk) belongs to the Tec family of nonreceptor protein tyrosine kinases. Mutations in the BTK gene cause X-linked agammaglobulinemia (XLA); a primary immunodeficiency disorder in human. No clear genotype-phenotype correlation has been established in XLA so far. To determine how differently mutations in BTK affect the severity of the disease and if BTK promoter polymorphic variant or intron 1 polymorphic variant in Tec, a cytoplasmic tyrosine kinase that might substitute for Btk, could contribute to the clinical phenotype, we analyzed the clinical and molecular findings in a cohort of XLA patients. Polymorphisms in BTK promoter and TEC intron 1 regions include substitutions of C>T (rs2071219) and T>C (rs2664019), respectively. Btk expression was evaluated by means of western immunoblotting and fluorescence-activated cell sorter analysis. Mutations were categorized as mild or severe and patients were evaluated for the clinical severity of disease. On the basis of the results, severe genotypes do not necessarily lead to severe phenotypes. More over, in a considerable number of patients with mild phenotype we showed a severe mutation with a tendency toward C substitution in the polymorphic site on TEC intron 1.

Immunologic evaluation of patients with recurrent ear, nose, and throat infections.

PURPOSE: In this study, we aimed to study the frequency of possible underlying immunodeficiency responsible for susceptibility to ear, nose, and throat (ENT) infection. MATERIALS AND METHODS: One hundred three (72 males and 31 females) consecutive children and adult patients with history of recurrent or chronic ENT infections, referred by otolaryngologists to the Department of Allergy and Clinical Immunology, Children's Medical Center, Tehran University of Medical Sciences (Tehran, Iran), were enrolled to the study from March 2003 to March 2006. For each patient, demographic information and medical histories of any ENT infections were collected by reviewing the patient's records. We measured immunoglobulin isotype concentrations and immunoglobulin (Ig) G subclasses by nephelometry and enzyme-linked immunosorbent assay methods, respectively. Of 103 patients, 75 received unconjugated pneumococcus polyvalent vaccine, and blood samples were taken before and 21 days after vaccination. Specific antibodies against whole pneumococcal antigens were measured using enzyme-linked immunosorbent assay method. Existence of bronchiectasis was confirmed in each patient using high resolution computed tomography scan. RESULTS: Among 103 patients, 17 (16.5%) patients were diagnosed to have defects in antibody-mediated immunity including 6 patients with immunoglobulin class deficiency (2 common variable deficiency and 4 IgA deficiency), 3 with IgG subclass deficiency (2 IgG2 and 1 IgG3), and 8 with specific antibody deficiency against polysaccharide antigens. In our series, bronchiectasis was detected in 5 cases associated with primary immunodeficiency. CONCLUSIONS: Long-standing history of ENT infections could be an alarm for ENT infections associated with primary antibody deficiency.

Progressive multifocal leukoencephalopathy in purine nucleoside phosphorylase deficiency.

Progressive multifocal leukoencephalopathy is a demyelinating disease caused by JC virus, an opportunistic infection of the central nervous system. Although the majority of cases are infected with the human immunodeficiency virus (HIV), other immunocompromised patients are also at risk. Purine nucleoside phosphorylase is an enzyme in the purine salvage pathway that reversibly converts inosine to hypoxanthine and guanosine to guanine. Purine nucleoside phosphorylase deficiency is a combined immunodeficiency with a profound cellular defect. Neurologic abnormalities are salient features of this syndrome. We describe for the first time a patient with this rare disorder presented with progressive multifocal leukoencephalopathy.

Mortality and morbidity in common variable immunodeficiency.

Common variable immunodeficiency (CVID) is a heterogeneous group of disorders, characterized by hypogammaglobulinemia, defective specific-antibody production resulting in recurrent bacterial infections. Delay in diagnosis and inadequate treatment result in increased irreversible complications and mortality. To determine persistent morbidities, mortality rate and survival in Iranian patients with CVID, hospital records of 72 (39 males and 33 females) diagnosed CVID patients were reviewed. Probabilities of survival after diagnosis of CVID were estimated from Kaplan-Meier life tables. Studied patients were enrolled over a 20-year period (1984-2005). The most commonly observed complication was bronchiectasis (24 cases), followed by splenomegaly, intestinal villous atrophy (11 cases), and failure to thrive (10 cases). Post-diagnosis survival was estimated as 65% for the first 6.5 years, which remains the same until 14 years after diagnosis when the survival curve drops to nearly 45%. The mortality rate among patients who had no regular visits and did not receive periodical IVIG was more remarkable when compared with those who had been followed up timely (p-value = 0.001). The most common cause of death was respiratory failure. Based on our observation, it can be highlighted that all patients with CVID, even under regular immunoglobulin replacement, need close monitoring for early detection of complications and introduction of appropriate management.

Clinical, immunological and molecular characteristics of 37 Iranian patients with X-linked agammaglobulinemia.

BACKGROUND: X-linked agammaglobulinemia (XLA) is a hereditary immunodeficiency characterized by an early onset of recurrent bacterial infections, a profound deficiency of all immunoglobulin isotypes and a markedly reduced number of peripheral B lymphocytes. Eighty-five percent of the patients with this phenotype have mutations in Bruton's tyrosine kinase (BTK) gene. METHODS: To provide an informative outlook of clinical and immunological manifestations of XLA in Iran, 37 Iranian male patients with an age range of 1-34 years, followed over a period of 25 years, were studied. Twenty-four of the 37 patients were screened for BTK gene mutation using PCR-SSCP followed by direct sequencing. BTK protein expression assay was done by flow cytometry in 9 families. RESULTS: All patients first presented with infectious diseases, the most common of which were respiratory tract infections. Eighteen different mutations were identified, 13 of which were novel: IVS1+5G>C, 1896G>A, 349delA, 1618C>T, 1783T>C, 2084A>G, 1346delT, 1351delGAG, 587A>G, IVS14-1G>A, IVS3+2T>C, 1482G>A, 1975C>A. CONCLUSION: The fact that we found a great number of novel mutations in a relatively limited number of patients underlines the heterogeneity of BTK mutations in the Iranian population. The large number of new mutations indicates that extended studies in this region would be rewarding.