Sequence-to-expression approach to identify etiological non-coding DNA variations in P53 and cMYC-driven diseases.

Disease risk prediction based on genomic sequence and transcriptional profile can improve disease screening and prevention. Despite identifying many disease-associated DNA variants, distinguishing deleterious non-coding DNA variations remains poor for most common diseases. In this study, we designed in vitro experiments to uncover the significance of occupancy and competitive binding between P53 and cMYC on common target genes. Analyzing publicly available ChIP-seq data for P53 and cMYC in embryonic stem cells showed that ~344-366 regions are co-occupied, and on average, two cis-overlapping motifs (CisOMs) per region were identified, suggesting that co-occupancy is evolutionarily conserved. Using U2OS and Raji cells untreated and treated with doxorubicin to increase P53 protein level while potentially reducing cMYC level, ChIP-seq analysis illustrated that around 16 to 922 genomic regions were co-occupied by P53 and cMYC, and substitutions of cMYC signals by P53 were detected post doxorubicin treatment. Around 187 expressed genes near co-occupied regions were altered at mRNA level according to RNA-seq data analysis. We utilized a computational motif-matching approach to illustrate that changes in predicted P53 binding affinity in CisOMs of co-occupied elements significantly correlate with alterations in reporter gene expression. We performed a similar analysis using SNPs mapped in CisOMs for P53 and cMYC from ChIP-seq data, and expression of target genes from GTEx portal. We found significant correlation between change in cMYC-motif binding affinity in CisOMs and altered expression. Our study brings us closer to developing a generally applicable approach to filter etiological non-coding variations associated with common diseases.

TWIST1 interacts with ß/δ-catenins during neural tube development and regulates fate transition in cranial neural crest cells.

Cell fate determination is a necessary and tightly regulated process for producing different cell types and structures during development. Cranial neural crest cells (CNCCs) are unique to vertebrate embryos and emerge from the neural plate borders into multiple cell lineages that differentiate into bone, cartilage, neurons and glial cells. We have previously reported that Irf6 genetically interacts with Twist1 during CNCC-derived tissue formation. Here, we have investigated the mechanistic role of Twist1 and Irf6 at early stages of craniofacial development. Our data indicate that TWIST1 is expressed in endocytic vesicles at the apical surface and interacts with ß/δ-catenins during neural tube closure, and Irf6 is involved in defining neural fold borders by restricting AP2α expression. Twist1 suppresses Irf6 and other epithelial genes in CNCCs during the epithelial-to-mesenchymal transition (EMT) process and cell migration. Conversely, a loss of Twist1 leads to a sustained expression of epithelial and cell adhesion markers in migratory CNCCs. Disruption of TWIST1 phosphorylation in vivo leads to epidermal blebbing, edema, neural tube defects and CNCC-derived structural abnormalities. Altogether, this study describes a previously uncharacterized function of mammalian Twist1 and Irf6 in the neural tube and CNCCs, and provides new target genes for Twist1 that are involved in cytoskeletal remodeling.

Molecular mechanisms of TWIST1-regulated transcription in EMT and cancer metastasis.

TWIST1 induces epithelial-to-mesenchymal transition (EMT) to drive cancer metastasis. It is yet unclear what determines TWIST1 functions to activate or repress transcription. We found that the TWIST1 N-terminus antagonizes TWIST1-regulated gene expression, cancer growth and metastasis. TWIST1 interacts with both the NuRD complex and the NuA4/TIP60 complex (TIP60-Com) via its N-terminus. Non-acetylated TWIST1-K73/76 selectively interacts with and recruits NuRD to repress epithelial target gene transcription. Diacetylated TWIST1-acK73/76 binds BRD8, a component of TIP60-Com that also binds histone H4-acK5/8, to recruit TIP60-Com to activate mesenchymal target genes and MYC. Knockdown of BRD8 abolishes TWIST1 and TIP60-Com interaction and TIP60-Com recruitment to TWIST1-activated genes, resulting in decreasing TWIST1-activated target gene expression and cancer metastasis. Both TWIST1/NuRD and TWIST1/TIP60-Com complexes are required for TWIST1 to promote EMT, proliferation, and metastasis at full capacity. Therefore, the diacetylation status of TWIST1-K73/76 dictates whether TWIST1 interacts either with NuRD to repress epithelial genes, or with TIP60-Com to activate mesenchymal genes and MYC. Since BRD8 is essential for TWIST1-acK73/76 and TIP60-Com interaction, targeting BRD8 could be a means to inhibit TWIST1-activated gene expression.

Dose-response effect of Montelukast on post-extraction dental socket repair and skeletal phenotype of mice.

Bone metabolism and repair are directly regulated by arachidonic acid metabolites. At present, we analyzed the dose-response effects of a selective cysteinyl leukotriene receptor type-1 antagonist during bone repair after tooth extraction and on non-injured skeleton. Sixty-three 129 Sv/Ev male mice composed the groups: C-Control (saline solution); MTK2-2 mg/Kg of Montelukast (MTK) and MTK4-4 mg/Kg of MTK, daily administered by mouth throughout all experimental periods set at 7, 14, and 21 days post-operative. Dental sockets were analyzed by computed microtomography (microCT), histopathology, and immunohistochemistry. Femurs, L5 vertebra and organs were also removed for observation. Blood was collected for plasma bone and liver markers. Histopathology and microCT analysis revealed early socket repair of MTK2 and MTK4 animals, with significant increased BV/TV at days 14 and 21 compared to C. Higher plasma calcium was detected at days 7 and 21 in MTK4 in comparison to C, while phosphate was significantly increased in MTK2 in the same periods in comparison to C and MTK4. No significant differences were found regarding plasma ALP and TRAP, neither for local TRAP and Runx2 immunolabeling at the healing sockets. Organs did not present histological abnormalities. Increased AST levels have been detected in distinct groups and periods. In general, femur phenotype was improved in MTK treated animals. Collectively, MTK promoted early bone formation after tooth extraction and increased bone quality of femurs and vertebra in a time-dose-dependent manner, and should be considered as an alternative therapy when improved post-extraction socket repair or skeleton preservation is required.

The TFAP2A-IRF6-GRHL3 genetic pathway is conserved in neurulation.

Mutations in IRF6, TFAP2A and GRHL3 cause orofacial clefting syndromes in humans. However, Tfap2a and Grhl3 are also required for neurulation in mice. Here, we found that homeostasis of Irf6 is also required for development of the neural tube and associated structures. Over-expression of Irf6 caused exencephaly, a rostral neural tube defect, through suppression of Tfap2a and Grhl3 expression. Conversely, loss of Irf6 function caused a curly tail and coincided with a reduction of Tfap2a and Grhl3 expression in tail tissues. To test whether Irf6 function in neurulation was conserved, we sequenced samples obtained from human cases of spina bifida and anencephaly. We found two likely disease-causing variants in two samples from patients with spina bifida. Overall, these data suggest that the Tfap2a-Irf6-Grhl3 genetic pathway is shared by two embryologically distinct morphogenetic events that previously were considered independent during mammalian development. In addition, these data suggest new candidates to delineate the genetic architecture of neural tube defects and new therapeutic targets to prevent this common birth defect.

Conservation of Epithelial-to-Mesenchymal Transition Process in Neural Crest Cells and Metastatic Cancer.

Epithelial to mesenchymal transition (EMT) is a highly conserved cellular process in several species, from worms to humans. EMT plays a fundamental role in early embryogenesis, wound healing, and cancer metastasis. For neural crest cell (NCC) development, EMT typically results in forming a migratory and potent cell population that generates a wide variety of cell and tissue, including cartilage, bone, connective tissue, endocrine cells, neurons, and glia amongst many others. The degree of conservation between the signaling pathways that regulate EMT during development and metastatic cancer (MC) has not been fully established, despite ample studies. This systematic review and meta-analysis dissects the major signaling pathways involved in EMT of NCC development and MC to unravel the similarities and differences. While the FGF, TGFß/BMP, SHH, and NOTCH pathways have been rigorously investigated in both systems, the EGF, IGF, HIPPO, Factor Receptor Superfamily, and their intracellular signaling cascades need to be the focus of future NCC studies. In general, meta-analyses of the associated signaling pathways show a significant number of overlapping genes (particularly ligands, transcription regulators, and targeted cadherins) involved in each signaling pathway of both systems without stratification by body segments and cancer type. Lack of stratification makes it difficult to meaningfully evaluate the intracellular downstream effectors of each signaling pathway. Finally, pediatric neuroblastoma and melanoma are NCC-derived malignancies, which emphasize the importance of uncovering the EMT events that convert NCC into treatment-resistant malignant cells.

Systematic review of hormonal and genetic factors involved in the nonsyndromic disorders of the lower jaw.

Mandibular disorders are among the most common birth defects in humans, yet the etiological factors are largely unknown. Most of the neonates affected by mandibular abnormalities have a sequence of secondary anomalies, including airway obstruction and feeding problems, that reduce the quality of life. In the event of lacking corrective surgeries, patients with mandibular congenital disorders suffer from additional lifelong problems such as sleep apnea and temporomandibular disorders, among others. The goal of this systematic review is to gather evidence on hormonal and genetic factors that are involved in signaling pathways and interactions that are potentially associated with the nonsyndromic mandibular disorders. We found that members of FGF and BMP pathways, including FGF8/10, FGFR2/3, BMP2/4/7, BMPR1A, ACVR1, and ACVR2A/B, have a prominent number of gene-gene interactions among all identified genes in this review. Gene ontology of the 154 genes showed that the functional gene sets are involved in all aspects of cellular processes and organogenesis. Some of the genes identified by the genome-wide association studies of common mandibular disorders are involved in skeletal formation and growth retardation based on animal models, suggesting a potential direct role as genetic risk factors in the common complex jaw disorders. Developmental Dynamics 248:162-172, 2019. © 2018 Wiley Periodicals, Inc.

A cleft lip and palate gene, Irf6, is involved in osteoblast differentiation of craniofacial bone.

BACKGROUND: Interferon regulatory factor 6 (IRF6) plays a critical role in embryonic tissue development, including differentiation of epithelial cells. Besides orofacial clefting due to haploinsufficiency of IRF6, recent human genetic studies indicated that mutations in IRF6 are linked to small mandible and digit abnormalities. The function of IRF6 has been well studied in oral epithelium; however, its role in craniofacial skeletal formation remains unknown. In this study, we investigated the role of Irf6 in craniofacial bone development using comparative analyses between wild-type (WT) and Irf6-null littermate mice. RESULTS: Immunostaining revealed the expression of IRF6 in hypertrophic chondrocytes, osteocytes, and bone matrix of craniofacial tissues. Histological analysis of Irf6-null mice showed a remarkable reduction in the number of lacunae, embedded osteocytes in matrices, and a reduction in mineralization during bone formation. These abnormalities may explain the decreased craniofacial bone density detected by micro-CT, loss of incisors, and mandibular bone abnormality of Irf6-null mice. To validate the autonomous role of IRF6 in bone, extracted primary osteoblasts from calvarial bone of WT and Irf6-null pups showed no effect on osteoblastic viability and proliferation. However, a reduction in mineralization was detected in Irf6-null cells. CONCLUSIONS: Altogether, these findings suggest an autonomous role of Irf6 in regulating bone differentiation and mineralization. Developmental Dynamics 248:221-232, 2019. © 2019 Wiley Periodicals, Inc.

The effect of non-coding DNA variations on P53 and cMYC competitive inhibition at cis-overlapping motifs.

Non-coding DNA variations play a critical role in increasing the risk for development of common complex diseases, and account for the majority of SNPs highly associated with cancer. However, it remains a challenge to identify etiologic variants and to predict their pathological effects on target gene expression for clinical purposes. Cis-overlapping motifs (COMs) are elements of enhancer regions that impact gene expression by enabling competitive binding and switching between transcription factors. Mutations within COMs are especially important when the involved transcription factors have opposing effects on gene regulation, like P53 tumor suppressor and cMYC proto-oncogene. In this study, genome-wide analysis of ChIP-seq data from human cancer and mouse embryonic cells identified a significant number of putative regulatory elements with signals for both P53 and cMYC. Each co-occupied element contains, on average, two COMs, and one common SNP every two COMs. Gene ontology of predicted target genes for COMs showed that the majority are involved in DNA damage, apoptosis, cell cycle regulation, and RNA processing. EMSA results showed that both cMYC and P53 bind to cis-overlapping motifs within a ChIP-seq co-occupied region in Chr12. In vitro functional analysis of selected co-occupied elements verified enhancer activity, and also showed that the occurrence of SNPs within three COMs significantly altered enhancer activity. We identified a list of COM-associated functional SNPs that are in close proximity to SNPs associated with common diseases in large population studies. These results suggest a potential molecular mechanism to identify etiologic regulatory mutations associated with common diseases.

Identification of functional variants for cleft lip with or without cleft palate in or near PAX7, FGFR2, and NOG by targeted sequencing of GWAS loci.

Although genome-wide association studies (GWASs) for nonsyndromic orofacial clefts have identified multiple strongly associated regions, the causal variants are unknown. To address this, we selected 13 regions from GWASs and other studies, performed targeted sequencing in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses. Within a cluster of strongly associated common variants near NOG, we found that one, rs227727, disrupts enhancer activity. We furthermore identified significant clusters of non-coding rare variants near NTN1 and NOG and found several rare coding variants likely to affect protein function, including four nonsense variants in ARHGAP29. We confirmed 48 de novo mutations and, based on best biological evidence available, chose two of these for functional assays. One mutation in PAX7 disrupted the DNA binding of the encoded transcription factor in an in vitro assay. The second, a non-coding mutation, disrupted the activity of a neural crest enhancer downstream of FGFR2 both in vitro and in vivo. This targeted sequencing study provides strong functional evidence implicating several specific variants as primary contributory risk alleles for nonsyndromic clefting in humans.

Call for action: how to improve use of patient-reported outcomes to guide clinical decision making in rheumatoid arthritis.

Current guidelines for the management of rheumatoid arthritis (RA) recommend early treatment and a treat-to-target goal of remission or low disease activity. Over the past decade, this approach has been extremely successful in reducing disease activity and joint damage in patients with RA. At the same time, however, overall patient perception of well-being appears to have decreased with respect to outcome measures considered important by patients themselves, such as pain, fatigue, physical function and quality of life. The timely and effective use of patient-reported outcomes (PROs) could encourage physicians to focus more on the impact of RA on patients and how patients are feeling. This in turn would facilitate shared decision making between patients and physicians, ultimately leading to a more patient-centered approach and improved patient care. Indeed, PROs provide information about individual patients that complements information provided by physical assessment and composite scores, and can also be used to guide patient care, such as determining whether a clinic visit is needed or whether treatment modifications are necessary. This is particularly important for patients who do not achieve the aspirational target of remission or low disease activity with pharmacological treatment. A number of validated PRO questionnaires are available, but how and which PROs should be incorporated into rheumatology clinical practice as part of the decision-making process is still controversial. Combining PROs with technology, such as computer adaptive tests, electronic PRO systems, web-based platforms and patient dashboards, could further aid PRO integration into daily rheumatology clinical practice.

An etiologic regulatory mutation in IRF6 with loss- and gain-of-function effects.

DNA variation in Interferon Regulatory Factor 6 (IRF6) causes Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate (CLP). However, an etiologic variant in IRF6 has been found in only 70% of VWS families. To test whether DNA variants in regulatory elements cause VWS, we sequenced three conserved elements near IRF6 in 70 VWS families that lack an etiologic mutation within IRF6 exons. A rare mutation (350dupA) was found in a conserved IRF6 enhancer element (MCS9.7) in a Brazilian family. The 350dupA mutation abrogated the binding of p63 and E47 transcription factors to cis-overlapping motifs, and significantly disrupted enhancer activity in human cell cultures. Moreover, using a transgenic assay in mice, the 350dupA mutation disrupted the activation of MCS9.7 enhancer element and led to failure of lacZ expression in all head and neck pharyngeal arches. Interestingly, disruption of the p63 Motif1 and/or E47 binding sites by nucleotide substitution did not fully recapitulate the effect of the 350dupA mutation. Rather, we recognized that the 350dupA created a CAAAGT motif, a binding site for Lef1 protein. We showed that Lef1 binds to the mutated site and that overexpression of Lef1/ß-Catenin chimeric protein repressed MCS9.7-350dupA enhancer activity. In conclusion, our data strongly suggest that 350dupA variant is an etiologic mutation in VWS patients and disrupts enhancer activity by a loss- and gain-of-function mechanism, and thus support the rationale for additional screening for regulatory mutations in patients with CLP.

Health-related quality of life outcomes with prasugrel among medically managed non-ST-segment elevation acute coronary syndrome patients: Insights from the Targeted Platelet Inhibition to Clarify the Optimal Strategy to Medically Manage Acute Coronary Syndromes (TRILOGY ACS) trial.

BACKGROUND: Few studies have assessed treatment effects on health-related quality of life (HRQoL) in patients with acute coronary syndrome (ACS) treated without revascularization. The TRILOGY ACS trial randomized patients with ACS to either prasugrel or clopidogrel therapy plus aspirin. Outcomes showed a complex pattern suggestive of late benefits with respect to repeat clinical events and benefits confined to patients who underwent angiography. Here, we examine the HRQoL correlates of these patterns. METHODS: HRQoL was measured at baseline and 3, 12, and 24 months or end of study (EOS) in 7243 patients aged <75 years using the EuroQol 3-level, group 5-dimension index (EQ-5D). Linear mixed effects models for repeated measures were used to examine treatment differences in HRQoL overall, stratified by angiography status, and among patients who did and did not have non-fatal events. RESULTS: No baseline differences in HRQoL were seen between patients randomized to prasugrel (n=3620) or clopidogrel (n=3623). At 24 months, remaining patients assigned to prasugrel (n=1450) vs. clopidogrel (n=1443) had higher EQ-5D index scores (86.4 vs. 84.9, P=.01). Mixed effects models found no difference in EQ-5D scores among prasugrel and clopidogrel patients overall across subgroups stratified by angiography status. However, among patients with non-fatal clinical events, patients on clopidogrel reported a larger decrement in HRQoL than patients on prasugrel (79.5±18.1 vs. 80.6±18.0; P=.02). CONCLUSIONS: Overall, there was no difference in HRQoL outcomes among patients receiving prasugrel vs. clopidogrel. However, the differential effects of the treatments among patients with non-fatal events require further investigation.

Analysis of the Relationship Between Micrognathia and Cleft Palate: A Systematic Review.

Objective To gather data from relevant experimental and observational studies to determine the relationship between micrognathia and cleft palate. The goal is to raise awareness and motivate clinicians to consider the cause and effect relationship when confronted with patients with cleft palate, even if there is no clearly noticeable mandibular abnormality. Design Several electronic databases were systematically examined to find articles for this review, using search terms including "cleft palate," "micrognathia," "tongue," and "airway obstruction." PubMed was the source of all the articles chosen to be included. Exclusion criteria included case reports, articles focused on treatment options, and articles only tangentially related to cleft palate and/or micrognathia. Results A total of 930 articles were screened for relevance, and 82 articles were chosen for further analysis. Evidence gathered in this review includes a variety of etiological factors that are causative or associated with both micrognathia and cleft palate. Observational studies relating the two abnormalities are also included. Much of the included literature recognizes a cause-and-effect relationship between micrognathia and cleft palate. Conclusion On the basis of the published data, we suggest that micrognathia does induce cleft palate in humans and animals. With knowledge of this causative relationship, clinicians should consider the importance of gathering cephalometric data on the mandibles and tongues of patients presenting with isolated cleft palate to determine whether they have micrognathia as well. With more data, patterns may emerge that could give insight into the complex etiology of nonsyndromic cleft palate.

Effects of a supraphysiological dose of testosterone cypionate on salivary gland function in adult male Wistar rats.

The abusive use of anabolic androgenic steroids has become a serious health problem worldwide, but its effects on oral health are still poorly understood. Therefore, the objective of this study was to evaluate the effects of a supraphysiological dose of testosterone cypionate (TC) on salivary biochemical, histomorphology, immunohistochemistry, and redox state parameters of parotid and submandibular glands. Twenty male Wistar rats, 12 weeks old, were divided into two groups (n=10/group): a control group and TC group, which received a dose of 20 mg/kg, once a week, for 6 weeks. Post treatment, the saliva and glands were collected. A supraphysiological dose of TC increased plasma and salivary testosterone concentrations. Although TC did not alter salivary flow, pH, and buffering capacity, the treatment increased the salivary secretion of total protein and reduced amylase, calcium, phosphate, and potassium. TC reduced the connective tissue area in the parotid gland and acinar area of the submandibular gland, while increasing the granular convoluted tubule area in the submandibular gland. Proliferating cell nuclear antigen was higher in the acinar cells of the submandibular glands from the TC group. Moreover, TC increased concentrations of total oxidant capacity and damaged lipids in both salivary glands, while total antioxidant activity and uric acid were lower in the submandibular gland, and reduced glutathione was higher in both glands. Superoxide dismutase, catalase, and glutathione peroxidase activities were higher in the parotid gland, while only glutathione peroxidase activity was lower in the submandibular gland of the TC group. In conclusion, TC abuse may be a potential factor for dysfunction of the parotid and submandibular glands, becoming a risk factor for the oral and systemic health of users.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6.

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

MCS9.7 enhancer activity is highly, but not completely, associated with expression of Irf6 and p63.

BACKGROUND: DNA variation in Interferon Regulatory Factor 6 (IRF6) contributes risk for orofacial clefting, including a common DNA variant rs642961. This DNA variant is located in a multi-species conserved sequence that is 9.7 kb upstream from the IRF6 transcriptional start site (MCS9.7). The MCS9.7 element was shown to possess enhancer activity that mimicked the expression of endogenous Irf6 at embryonic day 11.5 in transient transgenic embryos, and also contains a p63 binding site that transactivates IRF6 expression. To analyze whether the MCS9.7 enhancer is sufficient to drive IRF6 expression, we generated stable transgenic murine lines that carry a MCS9.7-lacZ transgene. We hypothesized that MCS9.7 was sufficient to recapitulate the endogenous expression of Irf6 at other time-points during embryonic development. RESULTS: We observed that MCS9.7 activity recapitulated endogenous Irf6 expression in most tissues, but not in the medial edge epithelium (MEE) at E14.5, when Irf6 expression was high during secondary palatal fusion. Also, while MCS9.7 activity and Irf6 expression were associated with p63 expression, we observed MCS9.7 activity and Irf6 expression in periderm, although p63 was absent. CONCLUSION: These data suggest that MCS9.7 enhancer activity is not sufficient to recapitulate IRF6 expression, and that p63 expression is not always necessary nor sufficient for transactivation of IRF6.

Sequence-to-expression approach to identify etiological non-coding DNA variations in P53 and cMYC-driven diseases.

Background and methods: Disease risk prediction based on DNA sequence and transcriptional profile can improve disease screening, prevention, and potential therapeutic approaches by revealing contributing genetic factors and altered regulatory networks. Despite identifying many disease-associated DNA variants through genome-wide association studies, distinguishing deleterious non-coding DNA variations remains poor for most common diseases. We previously reported that non-coding variations disrupting cis-overlapping motifs (CisOMs) of opposing transcription factors significantly affect enhancer activity. We designed in vitro experiments to uncover the significance of the co-occupancy and competitive binding and inhibition between P53 and cMYC on common target gene expression. Results: Analyzing publicly available ChIP-seq data for P53 and cMYC in human embryonic stem cells and mouse embryonic cells showed that ~ 344-366 genomic regions are co-occupied by P53 and cMYC. We identified, on average, two CisOMs per region, suggesting that co-occupancy is evolutionarily conserved in vertebrates. Our data showed that treating U2OS cells with doxorubicin increased P53 protein level while reducing cMYC level. In contrast, no change in protein levels was observed in Raji cells. ChIP-seq analysis illustrated that 16-922 genomic regions were co-occupied by P53 and cMYC before and after treatment, and substitutions of cMYC signals by P53 were detected after doxorubicin treatment in U2OS. Around 187 expressed genes near co-occupied regions were altered at mRNA level according to RNA-seq data. We utilized a computational motif-matching approach to determine that changes in predicted P53 binding affinity by DNA variations in CisOMs of co-occupied elements significantly correlate with alterations in reporter gene expression. We performed a similar analysis using SNPs mapped in CisOMs for P53 and cMYC from ChIP-seq data in U2OS and Raji, and expression of target genes from the GTEx portal. Conclusions: We found a significant correlation between change in motif-predicted cMYC binding affinity by SNPs in CisOMs and altered gene expression. Our study brings us closer to developing a generally applicable approach to filter etiological non-coding variations associated with P53 and cMYC-dependent diseases.

The RA-BE-REAL Multinational, Prospective, Observational Study in Patients with Rheumatoid Arthritis Receiving Baricitinib, Targeted Synthetic, or Biologic Disease-Modifying Therapies: a 6-Month Interim Analysis.

INTRODUCTION: RA-BE-REAL has the overall aim of defining a profile of patients with rheumatoid arthritis (RA) starting baricitinib or any other targeted synthetic (ts) or any biologic (b) disease-modifying antirheumatic drug (DMARD) for the first time, and the primary objective of estimating time until discontinuation from any cause (excluding sustained response) of the initial treatment. METHODS: RA-BE-REAL is an ongoing, prospective, observational, 36-month study in patients with RA initiating treatment with baricitinib (cohort A) or any other tsDMARD or any bDMARD (cohort B) for the first time. The primary objective is to assess the time until treatment discontinuation from any cause (excluding sustained response) at 24 months, (i.e., the rate of discontinuation of initial baricitinib or ts/bDMARD). Patient profiles of each cohort are described and compared. Post-baseline data are descriptively analyzed. This manuscript presents baseline and interim (6-month) outcomes for European patients with RA participating in the global RA-BE-REAL study. RESULTS: Data from 1074 patients (cohort A: 509; cohort B: 565) were analyzed. For cohorts A and B, respectively, the 6-month cumulative incidence (95% confidence interval) of treatment discontinuation was 16.5 (12.9-21.1) and 23.3 (19.1-28.2), and the proportions of patients achieving remission were 25.6% and 18.5%. At baseline, mean patient age was 59.1 and 57.0 years (p = 0.010) and mean disease duration was 10.0 and 8.9 years (p = 0.047), respectively. The proportions of patients exposed to ts/bDMARDs at any time before study entry were 51.9% and 39.1%, and the proportions of patients initiated on monotherapy were 50.9% and 31.2%, respectively. CONCLUSION: In real-world settings, patients with RA initiating treatment with baricitinib were older and had longer disease duration than those initiating treatment with any other tsDMARD or any bDMARD. Initial descriptive data regarding treatment discontinuation (including reasons for discontinuation), effectiveness, and treatment patterns will be enriched as the study progresses.

Assessment of the toxic effects of levetiracetam on biochemical, functional, and redox parameters of salivary glands in male Wistar rats.

Levetiracetam (LEV) is an anticonvulsant for epilepsy. The toxic effects of this medication in tissues have been associated with redox state imbalance, which can lead to salivary gland dysfunction. Therefore, the current work investigated the effects of LEV on the biochemical, functional, and redox parameters of the parotid and submandibular glands in rats. For this, male Wistar rats (Rattus norvegicus albinus) were randomly divided into 3 groups (n = 10/group): Control (0.9% saline solution), LEV100 (100 mg/kg), and LEV300 (300 mg/kg). After 21 consecutive days of intragastric gavage treatments, pilocarpine stimulated saliva secretion was collected for salivary biochemical analysis. The extracted salivary glands were utilized for histomorphometry and redox state analyses. Our results showed that LEV300 increased plasma hepatotoxicity markers and reduced salivary amylase activity and the acinar surface area of the parotid gland. Total oxidant capacity and oxidative damage to lipids and proteins were higher in the parotid gland, while total antioxidant capacity and uric acid levels were reduced in the submandibular gland of the LEV100 group compared to Control. On the other hand, total oxidant capacity, oxidative damage to lipids and proteins, total antioxidant capacity, and uric acid levels were lower in both salivary glands of the LEV300 group compared to Control. Superoxide dismutase and glutathione peroxidase activities were lower in the salivary glands of treated animals compared to Control. In conclusion our data suggest that treatment with LEV represents a potentially toxic agent, that contributes to drug-induced salivary gland dysfunction.