Association between ADAM33 Single-Nucleotide Polymorphisms and Treatment Response to Inhaled Corticosteroids and a Long-Acting Beta-Agonist in Asthma.

ADAM33 has been linked to airway structural changes in patients with asthma, leading to airway hyperresponsiveness, narrowing, and ultimately poor treatment responsiveness. This study aimed to evaluate the genetic association of ADAM33 SNPs with asthma, disease severity, and treatment responsiveness to ICS+LABA in the South Indian population. In this case-control study (486 controls and 503 cases), we performed genotyping using MassArray for six SNPs of ADAM33, namely rs2280091, rs2787094, rs3918396, rs67044, rs2853209, and rs3918392. We studied the association with asthma and treatment responsiveness to ICS+LABA, using genotype, allele frequency distribution, and haplotype analysis. A significant clinical finding of the study was that certain patients in the disease severity group (moderate and mild) showed poor or no improvement after a three-month follow-up of regular ICS+LABA therapy. Of the studied ADAM33 SNPs, rs2853209 showed an association with asthma. The further analysis of asthma patients according to disease severity suggested an association between moderate disease and the minor allele "T" for rs2853209. The homozygous minor allele of SNP rs2787094 was found to be associated with poorer lung function and the least lung-function improvement after three months of ICS+LABA therapy. The haplotype analysis of six SNPs showed a significant association between the rs2853209 and rs3918396 blocks and asthma. ADAM33 gene polymorphism has clinical relevance in terms of disease association and response to treatment. SNP rs2853209 seemed most relevant to asthma, and SNP rs2787094 could be a genetic marker for predicting response to ICS+LABA therapy in the study population.

Prioritizing Candidate Genes for Type 2 Diabetes Mellitus using Integrated Network and Pathway Analysis.

Background: Type 2 Diabetes Mellitus (T2DM) has emerged as a major threat to global health that fosters life-threatening clinical complications, taking a huge toll on our society. More than 65 million Indians suffer from T2DM, making it one of the leading causes of death. T2DM and associated complications have to be constantly monitored and managed which reduces the overall quality of life and increases socioeconomic burden. Therefore, it is crucial to develop specific treatment and management strategies. In order to achieve this, it is essential to understand the underlying genetic causes and molecular mechanisms. Methods: Integrated gene network and ontology analyses facilitate prioritization of plausible candidate genes for T2DM and also aid in understanding their mechanistic pathways. In this study, T2DM-associated genes were subjected to sequential interaction network and gene set enrichment analysis. High ranking network clusters were derived and their interrelation with pathways was assessed. Results: About 23 significant candidate genes were prioritized from 615 T2DM-associated genes which were overrepresented in pathways related to insulin resistance, type 2 diabetes, signaling cascades such as insulin receptor signaling pathway, PI3K signaling, IGFR signaling pathway, ERBB signaling pathway, MAPK signaling pathway and their regulatory mechanisms. Conclusion: Of these, two tyrosine kinase receptor genes-EGFR and IGF1R were identified as common nodes and can be considered to be significant candidate genes in T2DM.

Differential expression of the Mycobacterium tuberculosis heat shock protein genes in response to drug-induced stress.

Heat shock proteins are essential in maintaining cellular protein function, especially during stress. Their influence in managing drug-induced stress in Tuberculosis is not clearly understood. AIMS: Study the expression of select genes of the DnaK/ClpB chaperone network to evaluate their role in stress response in Mycobacterium tuberculosis clinical isolates during exposure to Isoniazid (INH) and Rifampicin (RIF). METHODS: Sanger sequencing to detect drug-resistant mutations followed by Drug Susceptibility Testing and Minimum Inhibitory Concentration determination. Culturing the bacilli in vitro, exposed to 1/4, 1/2 and 1 × MIC, and RNA quantification of dnaK, dnaJ1, grpE and clpB genes by using Real-time PCR. RESULTS: Susceptible isolates showed marginal down-regulation of two genes for INH, whereas all genes under-expressed against RIF. INH-resistant isolates had distinct expression profiles for inhA-15 and katG315 mutants. RIF-resistant bacilli did not have significant differential expression. MDR isolate showed up-regulation of all the four genes, with two genes over-expressing (≥4-fold). CONCLUSIONS: We observed characteristic gene expression profiles for each isolate in response to lethal and sub-lethal doses of INH and RIF. This provides insight into the role of DnaK/ClpB chaperone network in managing drug-induced stress and facilitating resistance. Further, the knowledge could provide targets for new drugs and augmenters.

Mucosa-associated specific bacterial species disrupt the intestinal epithelial barrier in the autism phenome.

Gut-Brain Axis provides a bidirectional communicational route, an imbalance of which can have pathophysiological consequences. Differential gut microbiome studies have become a frontier in autism research, affecting 85% of autistic children. The present study aims to understand how gut microbiota of autism subjects differ from their neurotypical counterparts. This study would help to identify the abundance of bacterial signature species in autism and their associated metabolites. 16S rRNA metagenomic sequence datasets of 30 out of 206 autism subjects were selected from the American Gut Project Archive. First, the taxonomic assignment was inferred by similarity-based methods using the Quantitative Insights into Microbial Ecology (QIIME) toolkit. Next, species abundance was characterized, and a co-occurrence network was built to infer species interaction using measures of diversity. Thirdly, statistical parameters were incorporated to validate the findings. Finally, the identification of metabolites associated with these bacterial signature species connects with biological processes in the host through pathway analysis. Gut microbiome data revealed Akkermansia sp. and Faecalibacterium prausnitzii to be statistically lower in abundance in autistic children than their neurotypical peers with a five and two-fold decrease, respectively. While Prevotella sp. and Sutterella sp. showed a five and a two-fold increase in cases, respectively. The constructed pathway revealed succinate and butyrate as the significant metabolites for the bacterial signature species identified. The present study throws light on the role of mucosa-associated bacterial species: Veillonella sp., Prevotella sp., Akkermansia sp., Sutterella sp., Faecalibacterium prausnitzii, Lactobacillus sp., which can act as diagnostic criteria for detection of gut dysbiosis in autism.

De novo assembly, annotation and gene expression profiles of gonads of Cytorace-3, a hybrid lineage of Drosophila nasuta nasuta and D. n. albomicans.

Cytorace-3 is a laboratory evolved hybrid lineage of Drosophila nasuta nasuta males and Drosophila nasuta albomicans females currently passing ~850 generations. To assess interracial hybridization effects on gene expression in Cytorace-3 we profiled the transcriptomes of mature ovaries and testes by employing Illumina sequencing technology and de novo transcriptome assembling strategies. We found 26% of the ovarian, and 14% of testis genes to be differentially expressed in Cytorace-3 relative to the expressed genes in the parental gonadal transcriptomes. About 5% of genes exhibited additive gene expression pattern in the ovary and 3% in the testis, while the remaining genes were misexpressed in Cytorace-3. Nearly 772 of these misexpressed genes in the ovary and 413 in the testis were either over- or under-dominant. Genes following D. n. nasuta dominance was twice (270 genes) than D. n. albomicans dominance (133 genes) in the ovary. In contrast, only 105 genes showed D. n. nasuta dominance and 207 showed D. n. albomicans dominance in testis transcriptome. Of the six expression inheritance patterns, conserved inheritance pattern was predominant for both ovary (73%) and testis (85%) in Cytorace-3. This study is the first to provide an overview of the expression divergence and inheritance patterns of the transcriptomes in an independently evolving distinct hybrid lineage of Drosophila. This recorded expression divergence in Cytorace-3 surpasses that between parental lineages illustrating the strong impact of hybridization driving rapid gene expression changes.

Integrated Network and Gene Ontology Analysis Identifies Key Genes and Pathways for Coronary Artery Diseases.

BACKGROUND: The prevalence of Coronary Artery Disease (CAD) in developing countries is on the rise, owing to rapidly changing lifestyle. Therefore, it is imperative that the underlying genetic and molecular mechanisms be understood to develop specific treatment strategies. Comprehensive disease network and Gene Ontology (GO) studies aid in prioritizing potential candidate genes for CAD and also give insights into gene function by establishing gene and disease pathway relationships. METHODS: In the present study, CAD-associated genes were collated from different data sources and protein-protein interaction network was constructed using STRING. Highly interconnected network clusters were inferred and GO analysis was performed. RESULTS: Interrelation between genes and pathways were analyzed on ClueGO and 38 candidates were identified from 1475 CAD-associated genes, which were significantly enriched in CAD-related pathways such as metabolism and regulation of lipid molecules, platelet activation, macrophage derived foam cell differentiation, and blood coagulation and fibrin clot formation. DISCUSSION: Integrated network and ontology analysis enables biomarker prioritization for common complex diseases such as CAD. Experimental validation and future studies on the prioritized genes may reveal valuable insights into CAD development mechanism and targeted treatment strategies.

Simultaneous Detection of Drug-resistant Mutations in Mycobacterium tuberculosis and Determining their Role through In Silico Docking.

A molecular method for diagnosis of drug-resistant Tuberculosis is Multiplex allele-specific PCR (MAS-PCR), which is more time-efficient, and its accuracy is studied using DNA sequencing. Also, understanding the role of mutations, when translated to protein, in causing resistance helps in better drug designing. AIM: To study MAS-PCR in the detection of drug resistance in comparison to DNA sequencing in Mycobacterium tuberculosis, and understand the mechanism of interaction of drugs with mutant proteins. METHODS: MAS-PCR was used for the detection of drug-resistant mutations and validation was done through DNA sequencing. MAS-PCR targeted four genes, iniA for the drug Ethambutol, rpsL and rrs for Streptomycin, and gyrA for Fluoroquinolone resistance, respectively. Further, the sequence data was analysed and modeled for in silico docking to study the effect on the interaction of the anti-TB drug molecule with the target protein. RESULTS: We identified drug-resistant mutations in four out of 95 isolates with one of them carrying a mutation at codon iniA501, two at gyrA94, and one for both iniA501 and gyrA94 using MASPCR. DNA sequencing confirmed drug-resistant mutations in only two isolates, whereas two others had mutation adjacent to the target allele. Drug-protein docking showed Estimated Free Energy of Binding to be higher for Fluoroquinolone binding with GyrA D94V mutant. Both wild and mutant IniA interact with EMB but had no significant effect on binding energy. CONCLUSION: DNA sequencing-based drug resistance detection of TB is more accurate than MASPCR. Understanding the role of mutations in influencing the drug-protein interaction will help in designing effective drug alternatives.

Allopatric sibling species pair Drosophila nasuta nasuta and Drosophila nasuta albomicans exhibit expression divergence in ovarian transcriptomes.

Drosophila nasuta nasuta and Drosophila nasuta albomicans represent cross fertile members of the immigrans species group of Drosophila with an allopatric mode of distribution exhibiting characteristic novelties. Illumina sequencing technology and de novo transcriptome assembling strategies were used for the current study. The analysis revealed 8% of the transcriptome to be differentially expressed between the ovaries of these two species, of which 9% was related to female reproduction. The majority of the differentially expressed genes were enriched for genetic information processing pathways, biosynthesis, and metabolism-related pathways. SNPs in D. n. albomicans transcriptome was double in number than in D. n. nasuta and only 5% of these SNPs were fixed. Ka/Ks ratios indicated the lineages were under strong purifying selection. The genes which are differentially expressed are evolving at a similar rate as that of genes with conserved expression. Thus, the current findings provide useful insights on the expression dynamics during incipient species divergence of D. n. nasuta and D. n. albomicans since their divergence time of ~ 0.5 million years.

Integrated Functional Analysis Implicates Syndromic and Rare Copy Number Variation Genes as Prominent Molecular Players in Pathogenesis of Autism Spectrum Disorders.

Autism Spectrum Disorders (ASD) are caused by disrupted neurodevelopment leading to socio-communication and behavioural abnormalities. Although genetic anomalies like Copy Number Variations (CNV) have been implicated in ASD, their overall genomic landscape and pathogenicity remain elusive. Therefore, we created a CNV map for ASD using 9337 cases and 5650 controls from SFARI database, statistically marked genomic regions with high and low frequencies of CNVs (i.e., common and rare CNV regions respectively), performed gene function enrichment for CNV genes, built functional networks, pathways and examined their expression in brain tissues. Information thus obtained were cumulatively integrated using a weighted scoring strategy to rank CNV regions by their neuro-functional attributes. Subsequently, we mapped 105 genic CNV regions across 20 chromosomes. They encompassed 537 genes, of which only 59 (11%) genes were identified with Single Nucleotide Variations (SNV) in ASD subjects through sequencing and functional studies, which indicated that diverse sets of genes were affected by CNVs and SNVs in ASD. Overall, syndromic CNV regions displayed the most prominent neuronal functions. While common CNV regions were found in loci 15q11.2, 16p11.2, 22q11.21, 15q13.2-13.3, rare CNV regions in loci 4p16.3, 9q34.3, 7q11.23, 17p11.2 contributed significantly to protein interaction networks and were highly expressed in brain. Enriched CNV genes were clustered in six functional categories with either direct roles in neurodevelopment or auxiliary roles like cellular signalling via MAPK pathway, cytoskeletal organization and transport or immune regulation. Mechanisms through which these molecular systems could independently or in combination trigger an ASD phenotype were predicted.

Haplotype analysis of ADAM33 polymorphisms in asthma: A pilot study.

Background & objectives: ADAM33 is implicated as a potentially strong candidate gene for asthma and bronchial hyper-responsiveness. Many polymorphisms of ADAM33 have been studied along with ADAM33 expression in various cells of the lungs. Haplotype analysis also showed association with asthma in different populations across the world. Therefore, the aim of this study was to perform a comprehensive screening of ADAM33 polymorphisms in adult patients with asthma. Methods: Thirty five polymorphisms of ADAM33 were genotyped in 55 patients with asthma and 53 controls. The association of single nucleotide polymorphisms (SNPs) and haplotypes with phenotypes of asthma was analysed. Results: The genotype, minor allele frequency, odds ratio and Hardy-Weinberg equilibrium did not show any significant difference among cases and controls. No association was found between SNPs of ADAM33 with the severity of asthma. Correlation analysis of ADAM33 SNPs to the phenotypes, based on clinical variables and allergen sensitization, did not show significant difference. Haplotype analysis showed that rs2280090 and rs2280091 were associated with asthma in the patient group. Interpretation & conclusions: Haplotype analysis showed an association of the two SNP variations with asthma. These SNPs lead to amino acid change and are prone to phosphorylation, which may affect expression levels and protein function of ADAM33 and asthma susceptibility.

Detection of First-Line Drug Resistance Mutations and Drug-Protein Interaction Dynamics from Tuberculosis Patients in South India.

BACKGROUND: Diagnosis of drug-resistant tuberculosis predominantly relies on culture-based drug susceptibility testing, which take weeks to produce a result and a more time-efficient alternative method is multiplex allele-specific PCR (MAS-PCR). Also, understanding the role of mutations in causing resistance helps better drug designing. AIMS: To evaluate the ability of MAS-PCR in the detection of drug resistance and to understand the mechanism of interaction of drugs with mutant proteins in Mycobacterium tuberculosis. METHODS: Detection of drug-resistant mutations using MAS-PCR and validation through DNA sequencing. MAS-PCR targeted five loci on three genes, katG 315 and inhA -15 for the drug isoniazid (INH), and rpoB 516, 526, and 531 for rifampicin (RIF). Furthermore, the sequence data were analyzed to study the effect on interaction of the anti-TB drug molecule with the target protein using in silico docking. RESULTS: We identified drug-resistant mutations in 8 out of 114 isolates with 2 of them as multidrug-resistant TB using MAS-PCR. DNA sequencing confirmed only six of these, recording a sensitivity of 85.7% and specificity of 99.3% for MAS-PCR. Molecular docking showed estimated free energy of binding (ΔG) being higher for RIF binding with RpoB S531L mutant. Codon 315 in KatG does not directly interact with INH but blocks the drug access to active site. CONCLUSIONS: We propose DNA sequencing-based drug resistance detection for TB, which is more accurate than MAS-PCR. Understanding the action of resistant mutations in disrupting the normal drug-protein interaction aids in designing effective drug alternatives.

Identifying the risk of producing aneuploids using meiotic recombination genes as biomarkers: A copy number variation approach.

BACKGROUND & OBJECTIVES: Aneuploids are the most common chromosomal abnormality in liveborns and are usually the result of non-disjunction (NDJ) in meiosis. Copy number variations (CNVs) are large structural variations affecting the human genome. CNVs influence critical genes involved in causing NDJ by altering their copy number which affects the clinical outcome. In this study influence of CNVs on critical meiotic recombination was examined using new computational technologies to assess their role in causing aneuploidy. METHODS: This investigation was based on the analysis of 12 random normal populations consisting of 1714 individuals for aneuploid causing genes under CNV effect. To examine the effect of CNVs on genes causing aneuploidy, meiotic recombination genes were analyzed using EnrichR, WebGestalt and Ingenuity Pathway Analysis (IPA). RESULTS: Forty three NDJ genes were found under CNV burden; IPA (Ingenuity Pathway Analysis) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis of CNV in meiotic recombination genes revealed a significant role of breast cancer gene 1, amyloid protein precursor, mitogen-activated protein kinase and nerve growth factor as key molecular players involved in causing aneuploidy. Interaction between these genes with other CNV-overlapping genes involved in cell cycle, recombination and meiosis might lead to increased incidences of aneuploidy. INTERPRETATION & CONCLUSIONS: The findings of this study implied that the effect of CNVs on normal genome contributed in amplifying the occurrences of chromosomal aneuploidies. The normal individuals consisting of variations in the susceptible genes causing aneuploids in the population remain undetected until the disorder genes express in the succeeding generations.

High-resolution arrays reveal burden of copy number variations on Parkinson disease genes associated with increased disease risk in random cohorts.

BACKGROUND: Parkinson disease (PD) is a neurological disease responsible for a considerable rate of mortality and morbidity in the society. Since the symptoms of the disease appear much later than the actual onset of neuron degeneration, a majority of cases remain undiagnosed until the manifestation of the symptoms. OBJECTIVES: In order to investigate the existence of such susceptibility in the population, we analyzed Copy Number Variation (CNV) influences on PD genes in 1715 individuals from 12 different populations. RESULTS: Overall, 16 CNV-PD genes, 3 known to be causal and 13 associated, were found to be significantly enriched. PARK2, was under heavy burden with ~1% of the population containing CNV in the exonic region. The impact of these genes on the genome and disease pathway was analyzed using several genome analysis tools. Protein interaction network of CNV-PD genes revealed a complex interaction of molecules forming a major hub by the α-Synuclein, whose direct interactors, LRRK2, PARK2 and ATP13A2 are under CNV influence. CONCLUSIONS: We hypothesize that CNVs may not be the initiating event in the pathogenesis of PD and remain latent until additional secondary hits are acquired and also propose novel genes that may fall under the PD pathway which contribute in pathogenesis.

A newly evolved Drosophila Cytorace-9 shows trade-off between longevity and immune response.

Species with an efficient immune system would be at an advantage to evade pathogenic challenges and adapt to an ever changing ecological niche. The upkeep of immunity is a costly affair, thus trade-offs between immunity and other life history traits are expected. However, studies on the relation between immunity and life span have yielded paradoxical results. Drosophila Cytoraces, being at different stages of evolutionary divergence, provide an excellent experimental model system to study how evolving populations gain novel traits in the absence of selection. We found that in the absence of pathogenic infections, the Cytorace-9 flies lived longer than those of Cytorace-3. However, when these Cytoraces were challenged with different pathogenic microbes, the trend was opposite. After infection with pathogens, the long-lived Cytorace-9 survived worse than the short lived Cytorace-3, which can be attributed to a reduction in its immune response. This study provides evidence to support the existence of a trade-off between life span and immunity.

Copy Number Variation of UGT 2B Genes in Indian Families Using Whole Genome Scans.

Background and Objectives. Uridine diphospho-glucuronosyltransferase 2B (UGT2B) is a family of genes involved in metabolizing steroid hormones and several other xenobiotics. These UGT2B genes are highly polymorphic in nature and have distinct polymorphisms associated with specific regions around the globe. Copy number variations (CNVs) status of UGT2B17 in Indian population is not known and their disease associations have been inconclusive. It was therefore of interest to investigate the CNV profile of UGT2B genes. Methods. We investigated the presence of CNVs in UGT2B genes in 31 members from eight Indian families using Affymetrix Genome-Wide Human SNP Array 6.0 chip. Results. Our data revealed >50% of the study members carried CNVs in UGT2B genes, of which 76% showed deletion polymorphism. CNVs were observed more in UGT2B17 (76.4%) than in UGT2B15 (17.6%). Molecular network and pathway analysis found enrichment related to steroid metabolic process, carboxylesterase activity, and sequence specific DNA binding. Interpretation and Conclusion. We report the presence of UGT2B gene deletion and duplication polymorphisms in Indian families. Network analysis indicates the substitutive role of other possible genes in the UGT activity. The CNVs of UGT2B genes are very common in individuals indicating that the effect is neutral in causing any suspected diseases.

Type 2 diabetes mellitus disease risk genes identified by genome wide copy number variation scan in normal populations.

AIMS: To identify the role of copy number variations (CNVs) on disease risk genes and its effect on disease phenotypes in type 2 diabetes mellitus (T2DM) in 12 random populations using high throughput arrays. METHODS: CNV analysis was carried out on a total of 1715 individuals from 12 populations, from ArrayExpress Archive of the European Bioinformatics Institute along with our subjects using Affymetrix Genome Wide SNP 6.0 array. CNV effect on T2DM genes were analyzed using several bioinformatics tools and a molecular protein interaction network was constructed to identify the disease mechanism altered by the CNVs. RESULTS: Analysis showed 34.4% of the total population to be under CNV burden for T2DM, with 83 disease causal and associated genes being under CNV influence. Hotspots were identified on chromosomes 22, 12, 6, 19 and 11.Overlap studies with case cohorts revealed significant disease risk genes such as EGFR, E2F1, PPP1R3A, HLA and TSPAN8. CONCLUSIONS: CNVs play a significant role in predisposing T2DM in normal cohorts and contribute to the phenotypic effects. Thus, CNVs should be considered as one of the major contributors in predisposition of the disease.

Development of Rapid Automatized Naming (RAN) in Simultaneous Kannada-English Biliterate Children.

RAN tests were administered to 600 typically developing children, 60 each from grade level one through grade ten (30 boys and 30 girls), who learn two distinct languages, English and Kannada simultaneously from the very first grade. The overall results were in accordance with similar previous studies in English and other European languages. The developmental trajectories were similar across the languages to a large extent; but the results also showed some differences across languages with respect to synchrony between the measures and the overall naming speed. Though some of the differences could be ascribed to the bilingual/biliterate culture and language use, there are enough scopes for future researches to examine these issues.

Global patterns of large copy number variations in the human genome reveal complexity in chromosome organization.

Global patterns of copy number variations (CNVs) in chromosomes are required to understand the dynamics of genome organization and complexity. For this study, analysis was performed using the Affymetrix Genome-Wide Human SNP Array 6.0 chip and CytoScan High-Density arrays. We identified a total of 44 109 CNVs from 1715 genomes with a mean of 25 CNVs in an individual, which established the first drafts of population-specific CNV maps providing a rationale for prioritizing chromosomal regions. About 19 905 ancient CNVs were identified across all chromosomes and populations at varying frequencies. CNV count, and sometimes CNV size, contributed to the bulk CNV size of the chromosome. Population specific lengthening and shortening of chromosomal length was observed. Sex bias for CNV presence was largely dependent on ethnicity. Lower CNV inheritance rate was observed for India, compared to YRI and CEU. A total of 33 candidate CNV hotspots from 5382 copy number (CN) variable region (CNVR) clusters were identified. Population specific CNV distribution patterns in p and q arms disturbed the assumption that CNV counts in the p arm are less common compared to long arms, and the CNV occurrence and distribution in chromosomes is length independent. This study unraveled the force of independent evolutionary dynamics on genome organization and complexity across chromosomes and populations.

IL-10 and IL-17F Promoter Single Nucleotide Polymorphism and Asthma: A Case-Control Study in South India.

BACKGROUND: Several studies have assessed the association between IL-17F and IL-10 promoter polymorphisms and asthma, but the results were conflicting. Furthermore, few studies have evaluated the association of cytokine polymorphisms with asthma and its clinical phenotypes. OBJECTIVE: This study was conducted to evaluate the association of IL-10 (interleukin 10) and IL-17F (interleukin 17F) promoter polymorphisms (rs1800871, rs1800896 and rs1889570) with asthma and its clinical phenotypes including severity, atopic status, spirometric parameters, and response to treatment in south Indian population. A sub-study was conducted to assess cytokine levels in subjects with different gene variants. METHODS: IL-10 and IL-17F polymorphisms were genotyped in 419 asthmatic patients and 393 controls using Mass ARRAY. RESULTS: Our results showed an association between IL-10 SNPs and mild asthma. No association was found with any of three SNPs in moderate to severe asthma. Comparison of genotype distribution of IL-17F rs1887570 AA variant among atopic and non-atopic patients showed significant difference (p = 0.024). Correlation analysis of IL-10 and IL-17F SNPs to clinical variables showed a positive correlation between IL-17F rs1887570 AA and number of allergen sensitized (rs = 0.142, p = 0.004). Significant improvement in lung function was observed after 2 months of ICS (Inhaled corticosteroids) and LABA (long acting ß2 agonist) treatment in all subjects with no statistically significant difference among SNPs variants. Cytokines levels were similar in different SNP variants. CONCLUSION: We observed an association between IL-10 rs1800871 and rs1800896 SNPs and mild asthma, as well as IL-17F rs1887570 AA variant and number of allergens sensitized.