PheroxyPyrabenz and Carbopyrropyridin against major proteins of SARS CoV-2: a comprehensive in-silico molecular docking and dynamics simulation studies

The pandemic that started in 2020 left us with so much information about viruses and respiratory diseases, and the cause behind it was severe acute respiratory syndrome coronavirus-2 (SARS CoV-2). The world is still recovering, which costs so many economic and other indirect disasters; despite that, no medications are available on the market. Although the WHO approved a few vaccines on an emergency basis, the remarks and the reinfection chances are still under investigation, and a few pharmaceutical companies are also claiming that a few medications can be effective. However, there is no situation in control. SARS CoV-2 mutates and comes in different forms, making the situation unpredictable. In this study, we have screened the complete Asinex's BioDesign library, which contains 170,269 compounds, and shorted the data against the docking score that helps in the identification of 4-[5-(3-Ethoxy-4-hydroxyphenyl)-1-(2-hydroxyethyl)-1H-pyrazol-3-yl]-1, 2-benzenediol (PheroxyPyrabenz) and 1-[(3R,4R)-1-(5-Aminopentanoyl)-4-hydroxy-3-pyrrolidinyl]-1H-pyrrolo[2,3-b]pyridine-4-carboxamide (Carbopyrropyridin) as a significant drug candidate that can work against the multiple proteins of the SARS CoV-2 resulting in seizing the complete biological process of the virus. Further, the study extended to Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and molecular dynamics (MD) simulation of both the compounds with their complexity. The complete workflow of the study has shown satisfactory results, and both drug candidates can potentially stop the hunt for drugs against this virus after its experimental validation. Further, we checked both compounds' absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties, showing case-proof validatory results.Communicated by Ramaswamy H. Sarma.

Identification of novel interacts partners of ADAR1 enzyme mediating the oncogenic process in aggressive breast cancer

Triple-negative breast cancer (TNBC) subtype is characterized by aggressive clinical behavior and poor prognosis patient outcomes. Here, we show that ADAR1 is more abundantly expressed in infiltrating breast cancer (BC) tumors than in benign tumors. Further, ADAR1 protein expression is higher in aggressive BC cells (MDA-MB-231). Moreover, we identify a novel interacting partners proteins list with ADAR1 in MDA-MB-231, using immunoprecipitation assay and mass spectrometry. Using iLoop, a protein-protein interaction prediction server based on structural features, five proteins with high iloop scores were discovered: Histone H2A.V, Kynureninase (KYNU), 40S ribosomal protein SA, Complement C4-A, and Nebulin (ranged between 0.6 and 0.8). In silico analysis showed that invasive ductal carcinomas had the highest level of KYNU gene expression than the other classifications (p < 0.0001). Moreover, KYNU mRNA expression was shown to be considerably higher in TNBC patients (p < 0.0001) and associated with poor patient outcomes with a high-risk value. Importantly, we found an interaction between ADAR1 and KYNU in the more aggressive BC cells. Altogether, these results propose a new ADAR-KYNU interaction as potential therapeutic targeted therapy in aggressive BC.

The diagnostic yield of CGH and WES in neurodevelopmental disorders

Background

Neurodevelopmental disorders are a group of conditions characterized by developmental delays leading to abnormal brain functions. The methods of diagnosis and treatment of these conditions are complicated, and their treatment involves a combination of various forms of therapy. In recent years, the development of high-resolution technologies has played an important role in revealing the microdeletions, microduplications, and single-nucleotide variants of the chromosomes and how they are linked to the development of neurodevelopmental disorders. The wide implementation and application of molecular methodologies have started to shed light on the functional importance of using the appropriate methods in detecting these genetic variations that are categorized as either pathogenic or benign. The study aimed to compare the diagnostic yield of comparative hybridization (CGH) and whole exome sequencing (WES) in neurodevelopmental disorders among children attending the King Abdullah Specialist Children Hospital, Riyadh, Saudi Arabia.

Methods

A retrospective study was conducted between 2015 and 2018 on 105 patients diagnosed with neurodevelopmental disorders through array-based CGH (Array-CGH) and WES.

Results

In a sample of 105 patients, 16% was the hit rate of copy number variations (CNVs). WES was requested for CNV-negative patients (n = 79), of which 30% was the hit rate of pathogenic or likely pathogenic single-nucleotide variants. There was a difference in the diagnostic yield between CGH (16%) and WES (30%).

Conclusion

WES was a better approach than Array-CGH to detect various DNA mutations or variants. Our findings could guide clinicians, researchers, and testing laboratories select the most cost-effective and appropriate approach for diagnosing their patients.

Congenital disorder of glycosylation with defective fucosylation 2 ( FCSK gene defect): The third report in the literature with a mild phenotype

BACKGROUND

Congenital disorders of glycosylation (CDG) are a group of heterogeneous disorders caused by abnormal lipid or protein glycosylation. Variants in the FCSK gene have been reported to cause CDG. Defective FCSK-induced CDG (FCSK-CDG) has only been reported previously in three unrelated children.

METHODS

In this study, we genetically and clinically examined a 3-year-old proband with resolved infantile spasms and normal development. Standard whole-exome sequencing (WES) and Sanger sequencing were performed to identify the functional impact of the variant.

RESULTS

WES revealed a rare biallelic missense variant (c.3013G>C; p.Val1005Leu) in FCSK. RT-qPCR showed a significant depletion in FCSK gene expression in the affected individual. Western blotting revealed reduced FCSK expression at the protein level compared to that in the control. Furthermore, 3D protein modeling suggested changes in the secondary structure, which might affect the overall FCSK protein function.

CONCLUSION

This study broadens the mutation and phenotypic spectrum of FCSK-associated developmental disorders.

Non-invasive prenatal testing for autosomal recessive disorders: A new promising approach

Background: In pregnant women at risk of autosomal recessive (AR) disorders, prenatal diagnosis of AR disorders primarily involves invasive procedures, such as chorionic villus sampling and amniocentesis.

Methods: We collected blood samples from four pregnant women in their first trimester who presented a risk of having a child with an AR disorder. Cell-free DNA (cfDNA) was extracted, amplified, and double-purified to reduce maternal DNA interference. Additionally, whole-genome amplification was performed for traces of residual purified cfDNA for utilization in subsequent applications.

Results: Based on our findings, we detected the fetal status with the family corresponding different genes, i.e., LZTR1, DVL2, HBB, RNASEH2B, and MYO7A, as homozygous affected, wild-type, and heterozygous carriers, respectively. Results were subsequently confirmed by prenatal amniocentesis. The results of AmpFLSTR™ Identifiler™ presented a distinct profile from the corresponding mother profile, thereby corroborating the result reflecting the genetic material of the fetus.

Conclusion: Herein, we detected AR disease mutations in the first trimester of pregnancy while surmounting limitations associated with maternal genetic material interference. Importantly, such detection strategies would allow the screening of pregnant women for common AR diseases, especially in highly consanguineous marriage populations. This technique would open avenues for the early detection and prevention of recessive diseases among the population.

Case Report: Bi-allelic missense variant in the desmocollin 3 gene causes hypotrichosis and recurrent skin vesicles

Background: Hypotrichosis with Recurrent Skin Vesicles (HYPTSV) is an extremely rare condition, having autosomal recessive inheritance. Here in we report a 4-years- old Saudi boy who presented with a history of recurrent skin blisters that are localized to the extremities and hypotrichosis since birth.

Methods: The present study describes a consanguineous Saudi family segregating HYPTSV in an autosomal recessive fashion. A single proband (II-1) exhibited features such as diffused non-scarring alopecia on the scalp, intraepidermal blister, post-inflammatory hyperpigmented macules, and follicular hyperkeratosis. DNA of the index was subjected to whole-genome sequencing (WGS). Furthermore, 3D protein modeling was performed for the mutated and normal protein.

Results: WGS revealed a novel bi-allelic missense variant (c.154G&gt;C; p. Val52Leu) in the DSC3 gene, which segregated perfectly using Sanger sequencing. In addition, 3D protein modeling revealed a substantial change in the mutated DSC3 protein as compared to the normal DSC3 protein.

Conclusion: This is the 3rd novel variant reported in the DSC3 gene associated with the HYPTSV phenotype. This report further strengthens the evidence that bi-allelic variants in the DSC3 cause severe HYPTSV in humans.

Novel homozygous pathogenic mitochondrial DNAJC19 variant in a patient with dilated cardiomyopathy and global developmental delay

BACKGROUND

Dilated cardiomyopathy with ataxia syndrome (DCMA) or 3-methylglutaconic aciduria type V is a rare global autosomal recessive mitochondrial syndrome that is clinically and genetically heterogeneous. It is characterized by early-onset dilated cardiomyopathy and increased urinary excretion of 3-methylglutaconic acid. As a result, some patients die due to cardiac failure, while others manifest with growth retardation, microcytic anemia, mild ataxia, and mild muscle weakness. DCMA is caused by variants in the DnaJ heat shock protein family (Hsp40) member C19 gene (DNAJC19), which plays an important role in mitochondrial protein import machinery in the inner mitochondrial membrane.

METHODS

We describe a single affected family member who presented with cardiomyopathy, global developmental delay, chest infection, seizures, elevated excretion of 3-methylglutaconic acid, and 3-methylglutaric acid in the urine.

RESULTS

Whole-exome sequencing followed by Sanger sequencing revealed a homozygous frameshift variant in the reading frame starting at codon 54 in exon 4 in the DNAJC19 gene (c.159del [Phe54Leufs*5]), which results in a stop codon four positions downstream. Quantitative gene expression analysis revealed that DNAJC19 mRNA expression in this patient was substantially reduced compared to the control.

CONCLUSIONS

We present a novel variant in the DNAJC19 gene that causes rare autosomal recessive mitochondrial 3-methylglutaconic aciduria type V. By comparing the current case with previously reported ones, we conclude that the disease is extremely heterogeneous for reasons that are still unknown.

Mutated VWA8 Is Associated With Developmental Delay, Microcephaly, and Scoliosis and Plays a Novel Role in Early Development and Skeletal Morphogenesis in Zebrafish

Von Willebrand A domain-containing protein 8 (VWA8), also named KIAA0564, is a poorly characterized, mitochondrial matrix-targeted protein having a putative ATPase activity. VWA8 is comprising of ATPase-associated domains and a VWFA domain associated with ATPase activity inside the cell. In the present study, we describe a large consanguineous family of Saudi origin segregating a complex developmental syndrome in an autosomal recessive fashion. All the affected individuals exhibited severe developmental disorders. DNA from three patients was subjected to whole-exome sequencing followed by Sanger sequencing. VWA8 knock-down zebrafish morpholinos were used to study the phenotypic effect of this gene on zebrafish development. A homozygous missense variant [c.947A > G; p.(Asp316Gly)] was identified in exon 8 of the VWA8 gene, which perfectly segregated with the disease phenotype. Using zebrafish morpholino, we observed delayed development at an early stage, lack of movement, light sensitivity, severe skeletal deformity such as scoliosis, and facial dysmorphism. This is the first homozygous variant identified in the VWA8 gene underlying global developmental delay, microcephaly, scoliosis, limbs, and cardiovascular malformations in humans. We provide genetic and molecular evidence using zebrafish morpholino for a homozygous variant in the VWA8 gene, associated with such a complex developmental syndrome in humans.

Herbal melanin induces interleukin-1β secretion and production by human THP-1 monocytes via Toll-like receptor 2 and p38 MAPK activation

Herbal melanin (HM), extracted from Nigella sativa, is known for its immunogenic properties through the modulation of cytokine production via Toll-like receptor (TLR)4. TLRs play a crucial role in the host defense through the regulation of innate and adaptive immune responses. However, the potential effect of HM on the production of interleukin-1β (IL-1β), the main immunoregulatory cytokine secreted by activated monocytes, has not been reported. The present study aimed to investigate the effects of HM on IL-1β secretion and production, detected by enzyme-linked immunosorbent assay, western blotting and mRNA expression monitored by reverse transcription-PCR, in human monocytes and a monocytic cell line, THP-1. Signaling pathways involved in the HM-induced IL-1β production was investigated in the THP-1 cells. It was shown that HM upregulated the IL-1β mRNA in the THP-1 cells and induced the secretion of IL-1β in the monocytes and THP-1 cells, in a dose-dependent manner, compared to the untreated cells. HM increased the protein expression of IL-1β, TLR2, the main receptor for IL-1β production, and activated p38 mitogen-activated protein kinase (MAPK), a key mediator for stress-induced IL-1β gene expression. The blockade of the p38 MAPK pathway, with the pharmacological inhibitor SB202190, and TLR2 receptor with a neutralization antibody, resulted in the decrease of HM-induced IL-1β production in THP-1 cells. The TLR4 receptor blockade also decreased HM-induced IL-1β production, but to a lesser extent than TLR2 blockade. In conclusion, the present study demonstrated that HM stimulates IL-1β production in monocytes and THP-1 cells, in a TLR2/p38 MAPK pathway-dependent manner, suggesting promising immunoregulatory potentials of HM against inflammatory-associated diseases.

Identification of the TTC26 Splice Variant in a Novel Complex Ciliopathy Syndrome with Biliary, Renal, Neurological, and Skeletal Manifestations

Ciliopathies constitute heterogeneous disorders that result from mutations in ciliary proteins. These proteins play an important role in the development of organs, physiology, and signaling pathways, and sequence variations in the genes encoding these proteins are associated with multisystem disorders. In this study, we describe a severe ciliopathy disorder that segregates in an autosomal recessive manner in a nonconsanguineous Saudi family. The proband exhibited features such as cholestasis, cystic dilatation of intrahepatic biliary ducts, diabetes insipidus, dysmorphic facial features, optic atrophy, pituitary hypoplasia, hydrocephalus, aqueductal stenosis, hyperextensible knee joints, bilateral knee dislocation, polydactyly, and syndactyly. Whole-genome sequencing and Sanger sequencing revealed a homozygous splice site variant (c.4-1G>C; NM_024926.3) in the tetratricopeptide repeat domain 26 (TTC26) gene located in chromosome 7q34, which cosegregated perfectly with the disease phenotype. qRT-PCR revealed a substantial decrease in the expression of the TTC26 gene as compared to the normal control, suggesting the pathogenicity of the identified variant. This report further strengthens the evidence that homozygous variants in the TTC26 gene cause severe ciliopathies with diverse phenotypes. We named this newly characterized condition as BRENS syndrome, which stands for biliary, renal, neurological, and skeletal features.

Next-Generation Sequencing-Based Pre-Implantation Genetic Testing for Aneuploidy (PGT-A): First Report from Saudi Arabia

Recently, high-throughput next-generation sequencing (NGS)-based preimplantation genetic testing for aneuploidies techniques came into use. This technique is essential for successful embryo transfer and accomplishing pregnancy, thus reducing the time and cost of additional cycles. In this study, we describe our first experience in introducing an NGS-based preimplantation genetic testing for aneuploidy (PGT-A) service using next-generation sequencing in King Abdulaziz Medical City located in Riyadh, Saudi Arabia. Our main goal was to report the successful implementation of this new technology in clinical practice and highlight the factors that may affect the results. In total, 200 blastomere biopsies were obtained from 36 in vitro fertilization (IVF) cycles from Saudi couples suffering from prolonged infertility or recurrent embryo transfer failure. NGS-based PGT-A was performed in all embryos. The results were analyzed in five age groups, showing that aneuploidy rates increased with maternal age. Moreover, the results also showed that complex abnormal embryos with (2-5) aneuploidy are the most common type of embryos. Additionally, our data showed that chromosome 16-related abnormality was the most frequent abnormality detected among all reported abnormalities. In conclusion, our study suggests that NGS-based PGT-A is an applicable and reliable technique for routine-based embryo screening, especially for couples suffering from recurrent miscarriages or multiple embryo transfer failures.

Pancytopenia, Recurrent Infection, Poor Wound Healing, Heterotopia of the Brain Probably Associated with A Candidate Novel de Novo CDC42 Gene Defect: Expanding the Molecular and Phenotypic Spectrum

CDC42 (cell division cycle protein 42) belongs to the Rho GTPase family that is known to control the signaling axis that regulates several cellular functions, including cell cycle progression, migration, and proliferation. However, the functional characterization of the CDC42 gene in mammalian physiology remains largely unclear. Here, we report the genetic and functional characterization of a non-consanguineous Saudi family with a single affected individual. Clinical examinations revealed poor wound healing, heterotopia of the brain, pancytopenia, and recurrent infections. Whole exome sequencing revealed a de novo missense variant (c.101C > A, p.Pro34Gln) in the CDC42 gene. The functional assays revealed a substantial reduction in the growth and motility of the patient cells as compared to the normal cells control. Homology three-dimensional (3-D) modeling of CDC42 revealed that the Pro34 is important for the proper protein secondary structure. In conclusion, we report a candidate disease-causing variant, which requires further confirmation for the etiology of CDC42 pathogenesis. This represents the first case from the Saudi population. The current study adds to the spectrum of mutations in the CDC42 gene that might help in genetic counseling and contributes to the CDC42-related genetic and functional characterization. However, further studies into the molecular mechanisms that are involved are needed in order to determine the role of the CDC42 gene associated with aberrant cell migration and immune response.

Next Generation Sequencing Based Non-invasive Prenatal Testing (NIPT): First Report From Saudi Arabia

Background: Non-invasive prenatal testing (NIPT) for aneuploidy in pregnant women screening has been recently established in Saudi Arabia. We aim from this study to report our experience in the implementation of this new technology in clinical practice and to assess factors influencing cell-free fetal (cffDNA) fraction and successful NIPT reporting.

Methods: In total, 200 pregnant women were subjected to the NIPT test using standard methods. Next-generation sequencing (NGS) was used to analyze cffDNA in maternal plasma.

Results: Out of the 200 NIPT cases, the average age of pregnant women was 35 ± 6 years (range: 21–48 years). The average cffDNA fraction of reported cases was 13.72% (range: 3–31%). Out of these 200 cases, 187 (93.5%) were at low risk, while 13 (6.5%) cases revealed high risk for aneuploidy. Among these chromosomal abnormalities, 7 (3.5%) cases of Down’s syndrome, 5 (2.5%) Edwards’ Syndrome, and only 1 case of (0.5%) Patau’s syndrome was observed. Out of the 13 high-risk cases, 2 (15.3%) were found in women below the age of 30.

Conclusion: This is the first study reporting the successful implementation of an in-house NIPT screening service in Saudi Arabia. Our data showed high accuracy and sensitivity to detect high-risk cases indicating the usefulness of such a technique as an alternative to invasive testing and (hopefully) will change the common screening practice for pregnant women in Saudi Arabia.

EMC10 homozygous variant identified in a family with global developmental delay, mild intellectual disability, and speech delay

In recent years, several genes have been implicated in the variable disease presentation of global developmental delay (GDD) and intellectual disability (ID). The endoplasmic reticulum membrane protein complex (EMC) family is known to be involved in GDD and ID. Homozygous variants of EMC1 are associated with GDD, scoliosis, and cerebellar atrophy, indicating the relevance of this pathway for neurogenetic disorders. EMC10 is a bone marrow‐derived angiogenic growth factor that plays an important role in infarct vascularization and promoting tissue repair. However, this gene has not been previously associated with human disease. Herein, we describe a Saudi family with two individuals segregating a recessive neurodevelopmental disorder. Both of the affected individuals showed mild ID, speech delay, and GDD. Whole‐exome sequencing (WES) and Sanger sequencing were performed to identify candidate genes. Further, to elucidate the functional effects of the variant, quantitative real‐time PCR (RT‐qPCR)‐based expression analysis was performed. WES revealed a homozygous splice acceptor site variant (c.679‐1G>A) in EMC10 (chromosome 19q13.33) that segregated perfectly within the family. RT‐qPCR showed a substantial decrease in the relative EMC10 gene expression in the patients, indicating the pathogenicity of the identified variant. For the first time in the literature, the EMC10 gene variant was associated with mild ID, speech delay, and GDD. Thus, this gene plays a key role in developmental milestones, with the potential to cause neurodevelopmental disorders in humans.

Mutated RAP1GDS1 causes a new syndrome of dysmorphic feature, intellectual disability & speech delay

Background

RAP1GDS1 (RAP1, GTP‐GDP dissociation stimulator 1), also known as SmgGDS, is a guanine nucleotide exchange factor (GEF) that regulates small GTPases, including, RHOA, RAC1, and KRAS. RAP1GDS1 was shown to be highly expressed in different tissue types including the brain. However, mutations in the RAP1GDS1 gene associated with human diseases have not previously been reported.

Methods

We report on four affected individuals, presenting intellectual disability, global developmental delay (GDD), and hypotonia. The probands’ DNA was subjected to whole‐genome sequencing, revealing a homozygous splice acceptor site mutation in the RAP1GDS1 gene (1444‐1G > A). Sanger sequencing was performed to confirm the segregation of the variant in two Saudi families. The possible aberrant splicing in the patients’ RNA was investigated using RT‐PCR and changes in mRNA expression of the patients were confirmed using qRT‐PCR.

Results

The identified splice variant was found to segregate within the two families. RT‐PCR showed that the mutation affected RAP1GDS1 gene splicing, resulting in the production of aberrant transcripts in the affected individuals. Quantitative gene expression analysis demonstrated that the RAP1GDS1 mRNA expression in all the probands was significantly decreased compared to that of the control, and Sanger sequencing of the probands’ cDNA revealed skipping of exon 13, further strengthening the pathogenicity of this variant.

Conclusion

We are the first to report the mutation of the RAP1GDS1 gene as a potential cause of GDD and hypotonia. However, further investigations into the molecular mechanisms involved are required to confirm the role of RAP1GDS1 gene in causing GDD and hypotonia.

MYT1L mutation in a patient causes intellectual disability and early onset of obesity: a case report and review of the literature

Background Obesity has become one of the greatest health risks worldwide. Recently, there was an explosion of information regarding the role of the central nervous system (CNS) in the development of monogenic and syndromic obesity. Case presentation Over the last decade, terminal and interstitial submicroscopic deletions of copy number variants (CNVs) in 2p25.3 and single nucleotide variants (SNVs) in myelin transcription factor 1 like (MYT1L) were detected by genome-wide array analysis and whole exome sequencing (WES) in patients with a nonspecific clinical phenotype that commonly includes intellectual disability (ID), early onset of obesity and speech delay. Here, we report the first Saudi female patient with mild to moderate ID, early onset of obesity and speech delay associated with a de novo pathogenic SNV in the MYT1L gene (c. 1585G>A [Gly529Arg]), which causes an amino acid change from Gly to Arg at position 529 that leads to mental retardation, autosomal dominant 39.

X chromosome dosage and presence of SRY shape sex-specific differences in DNA methylation at an autosomal region in human cells

BACKGROUND

Sexual dimorphism in DNA methylation levels is a recurrent epigenetic feature in different human cell types and has been implicated in predisposition to disease, such as psychiatric and autoimmune disorders. To elucidate the genetic origins of sex-specific DNA methylation, we examined DNA methylation levels in fibroblast cell lines and blood cells from individuals with different combinations of sex chromosome complements and sex phenotypes focusing on a single autosomal region--the differentially methylated region (DMR) in the promoter of the zona pellucida binding protein 2 (ZPBP2) as a reporter.

RESULTS

Our data show that the presence of the sex determining region Y (SRY) was associated with lower methylation levels, whereas higher X chromosome dosage in the absence of SRY led to an increase in DNA methylation levels at the ZPBP2 DMR. We mapped the X-linked modifier of DNA methylation to the long arm of chromosome X (Xq13-q21) and tested the impact of mutations in the ATRX and RLIM genes, located in this region, on methylation levels. Neither ATRX nor RLIM mutations influenced ZPBP2 methylation in female carriers.

CONCLUSIONS

We conclude that sex-specific methylation differences at the autosomal locus result from interaction between a Y-linked factor SRY and at least one X-linked factor that acts in a dose-dependent manner.

Local genotype influences DNA methylation at two asthma-associated regions, 5q31 and 17q21, in a founder effect population

BACKGROUND

Two asthma-associated regions 17q12-q21 and 5q31.1 harbour genes that show strong effect of genotype on expression levels. DNA methylation has an important role in gene regulation; therefore, we examined DNA methylation at promoters of 12 genes from 5q31 and 17q12-q21 regions. Our goal was to determine whether DNA methylation was associated with predisposition to asthma and whether such a relationship was independent from genetic association.

METHODS

Using sodium bisulfite sequencing and pyrosequencing methylation assays, we examined the effect of genotype on DNA methylation in peripheral blood cells from individuals from the Saguenay-Lac-Saint-Jean asthma familial collection and lymphoblastoid cell lines.

RESULTS

The local genotype influenced methylation levels of solute carrier family 22 (organic 3 cation/carnitine transporter) member 5 (SLC22A5), zona pellucida binding protein 2 (ZPBP2) and gasdermin A (GSDMA) promoter regions. The genotype had a dominant effect on ZPBP2 and GSDMA methylation with lower methylation levels in individuals that carry the asthma-predisposing alleles. Males also had lower methylation at the ZPBP2 promoter than females. We did not observe an effect of asthma status that would be independent of the genotype and the sex effects in the GSDMA, ZPBP2 and SLC22A5 regions; however, GSDMA and ZPBP2 data were suggestive of interaction between asthma and methylation levels in females and SLC22A5 in males.

CONCLUSIONS

The local genotype influences methylation levels at SLC22A5 and ZPBP2 promoters independently of the asthma status. Further studies are necessary to confirm the relationship between GSDMA-ZPBP2 and SLC22A5 methylation and asthma in females and males separately.

Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines

Phenotypic variation results from variation in gene expression, which is modulated by genetic and/or epigenetic factors. To understand the molecular basis of human disease, interaction between genetic and epigenetic factors needs to be taken into account. The asthma-associated region 17q12-q21 harbors three genes, the zona pellucida binding protein 2 (ZPBP2), gasdermin B (GSDMB) and ORM1-like 3 (ORMDL3), that show allele-specific differences in expression levels in lymphoblastoid cell lines (LCLs) and CD4+ T cells. Here, we report a molecular dissection of allele-specific transcriptional regulation of the genes within the chromosomal region 17q12-q21 combining in vitro transfection, formaldehyde-assisted isolation of regulatory elements, chromatin immunoprecipitation and DNA methylation assays in LCLs. We found that a single nucleotide polymorphism rs4795397 influences the activity of ZPBP2 promoter in vitro in an allele-dependent fashion, and also leads to nucleosome repositioning on the asthma-associated allele. However, variable methylation of exon 1 of ZPBP2 masks the strong genetic effect on ZPBP2 promoter activity in LCLs. In contrast, the ORMDL3 promoter is fully unmethylated, which allows detection of genetic effects on its transcription. We conclude that the cis-regulatory effects on 17q12-q21 gene expression result from interaction between several regulatory polymorphisms and epigenetic factors within the cis-regulatory haplotype region.