Integrative analysis of gene expression and DNA methylation to identify biomarkers of non-genital warts induced by low-risk human papillomaviruses infection

Background

Human papillomaviruses have been shown to dysregulate the gene expression and DNA methylation profiles of their host cells over the course of infection. However, there is a lack of information on the impact of low-risk HPV infection and wart formation on host cell's expression and methylation patterns. Therefore, the objective of this study is to analyse the genome and methylome of common warts using an integrative approach.

Methods

In the present study, gene expression (GSE136347) and methylation (GSE213888) datasets of common warts were obtained from the GEO database. Identification of the differentially expressed and differentially methylated genes was carried out using the RnBeads R package and the edgeR Bioconductor package. Next, functional annotation of the identified genes was obtained using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). Network construction and analyses of the gene-gene, protein-protein, and signaling interactions of the differentially expressed and differentially methylated genes was performed using the GeneMANIA web interface, the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database, and the Signaling Network Open Resource 2.0 (SIGNOR 2.0), respectively. Lastly, significant hub genes were identified using the Cytoscape application CytoHubba.

Results

A total of 276 genes were identified as differentially expressed and differentially methylated in common warts, with 52% being upregulated and hypermethylated. Functional enrichment analysis identified extracellular components as the most enriched annotations, while network analyses identified ELN, ITGB1, TIMP1, MMP2, LGALS3, COL1A1 and ANPEP as significant hub genes.

Conclusions

To the best knowledge of the authors, this is the first integrative study to be carried out on non-genital warts induced by low-risk HPV types. Future studies are required to re-validate the findings in larger populations using alternative approaches.

Pathway Analysis of Patients with Severe Acute Respiratory Syndrome

BACKGROUND

Coronaviruses are emerging threats for human health, as demonstrated by the ongoing coronavirus disease 2019 (COVID-19) pandemic that is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). SARS-CoV-2 is closely related to SARS-CoV-1, which was the cause of the 2002-2004 SARS outbreak, but SARS-CoV-1 has been the subject of a relatively limited number of studies. Understanding the potential pathways and molecular targets of SARS-CoV-1 will contribute to current drug repurposing strategies by helping to predict potential drug-disease associations.

METHODS

A microarray dataset, GSE1739, of 10 SARS patients and 4 healthy controls was downloaded from NCBI's GEO repository, and differential expression was identified using NCBI's GEO2R software. Pathway and enrichment analysis of the differentially expressed genes was carried out using Ingenuity Pathway Analysis and Gene Set Enrichment Analysis, respectively.

RESULTS

Our findings show that the drugs dexamethasone, filgrastim, interferon alfacon-1, and levodopa were among the most significant upstream regulators of differential gene expression in SARS patients, while neutrophil degranulation was the most significantly enriched pathway.

CONCLUSION

An enhanced understanding of the pathways and molecular targets of SARS-CoV-1 in humans will contribute to current and future drug repurposing strategies, which are an essential tool to combat rapidly emerging health threats.

Comparing metabolic profiles between female endurance athletes and non-athletes reveals differences in androgen and corticosteroid levels

Endurance training is associated with physiological changes in elite athletes, but little is known about female-specific effects of endurance training. Despite the significant rise in female sports participation, findings from studies performed on male athletes are largely extrapolated to females without taking into consideration sex-specific differences in metabolism. Subsequently, this study aimed to investigate the steroid hormone profiles of elite female endurance athletes in comparison with their non-athletic counterparts. Untargeted metabolomics-based mass spectroscopy combined with ultra-high-performance liquid chromatography was performed on serum samples from 51 elite female endurance athletes and 197 non-athletic females. The results showed that, compared to non-athletic females, certain androgen, pregnenolone, and progestin steroid hormones were reduced in elite female endurance athletes, while corticosteroids were elevated. The most significantly altered steroid hormones were 5alpha-androstan-3alpha,17alpha-diol monosulfate (FDR = 1.90 × 10^-05), androstenediol (3alpha, 17alpha) monosulfate (FDR = 2.93 × 10^-04), and cortisol (FDR = 2.93 × 10^-04). Conclusively, the present study suggests that elite female endurance athletes have a unique steroid hormone profile with implications on their general health and performance.

Pathway analysis of smoking-induced changes in buccal mucosal gene expression

Background

Cigarette smoking is the leading preventable cause of death worldwide, and it is the most common cause of oral cancers. This study aims to provide a deeper understanding of the molecular pathways in the oral cavity that are altered by exposure to cigarette smoke.

Methods

The gene expression dataset (accession number GSE8987, GPL96) of buccal mucosa samples from smokers (n = 5) and never smokers (n = 5) was downloaded from The National Center for Biotechnology Information's (NCBI) Gene Expression Omnibus (GEO) repository. Differential expression was ascertained via NCBI's GEO2R software, and Ingenuity Pathway Analysis (IPA) software was used to perform a pathway analysis.

Results

A total of 459 genes were found to be significantly differentially expressed in smoker buccal mucosa (p  < 0.05). A total of 261 genes were over-expressed while 198 genes were under-expressed. The top canonical pathways predicted by IPA were nitric oxide and reactive oxygen production at macrophages, macrophages/fibroblasts and endothelial cells in rheumatoid arthritis, and thyroid cancer pathways. The IPA upstream analysis predicted that the TP53, APP, SMAD3, and TNF proteins as well as dexamethasone drug would be top transcriptional regulators.

Conclusions

IPA highlighted critical pathways of carcinogenesis, mainly nitric oxide and reactive oxygen production at macrophages, and confirmed widespread injury in the buccal mucosa due to exposure to cigarette smoke. Our findings suggest that cigarette smoking significantly impacts gene pathways in the buccal mucosa and may highlight potential targets for treating the effects of cigarette smoking.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43042-022-00268-y.

Differential gene expression analysis of ankylosing spondylitis shows deregulation of the HLA-DRB, HLA-DQB, ITM2A, and CTLA4 genes

Background

Ankylosing spondylitis (AS) is a rare inflammatory disorder affecting the spinal joints. Although we know some of the genetic factors that are associated with the disease, the molecular basis of this illness has not yet been fully elucidated, and the genes involved in AS pathogenesis have not been entirely identified. The current study aimed at constructing a gene network that may serve as an AS gene signature and biomarker, both of which will help in disease diagnosis and the identification of therapeutic targets. Previously published gene expression profiles of 16 AS patients and 16 gender- and age-matched controls that were profiled on the Illumina HumanHT-12 V3.0 Expression BeadChip platform were mined. Patients were Portuguese, 21 to 64 years old, were diagnosed based on the modified New York criteria, and had Bath Ankylosing Spondylitis Disease Activity Index scores > 4 and Bath Ankylosing Spondylitis Functional Index scores > 4. All patients were receiving only NSAIDs and/or sulphasalazine. Functional enrichment and pathway analysis were performed to create an interaction network of differentially expressed genes.

Results

ITM2A, ICOS, VSIG10L, CD59, TRAC, and CTLA-4 were among the significantly differentially expressed genes in AS, but the most significantly downregulated genes were the HLA-DRB6, HLA-DRB5, HLA-DRB4, HLA-DRB3, HLA-DRB1, HLA-DQB1, ITM2A, and CTLA-4 genes. The genes in this study were mostly associated with the regulation of the immune system processes, parts of cell membrane, and signaling related to T cell receptor and antigen receptor, in addition to some overlaps related to the IL2 STAT signaling, as well as the androgen response. The most significantly over-represented pathways in the data set were associated with the “RUNX1 and FOXP3 which control the development of regulatory T lymphocytes (Tregs)” and the “GABA receptor activation” pathways.

Conclusions

Comprehensive gene analysis of differentially expressed genes in AS reveals a significant gene network that is involved in a multitude of important immune and inflammatory pathways. These pathways and networks might serve as biomarkers for AS and can potentially help in diagnosing the disease and identifying future targets for treatment.

Global gene methylation profiling of common warts caused by human papillomaviruses infection

Infection with the human papillomaviruses (HPV) often involves the epigenetic modification of the host genome. Despite its prevalence among the population, host genome methylation in HPV-induced warts is not clearly understood. In this study, genome-wide methylation profiling was carried out on paired healthy skin and wart samples in order to investigate the effects that benign HPV infection has on gene methylation status. To overcome this gap in knowledge, paired wart (n = 12) and normal skin (n = 12) samples were obtained from Arab males in order to perform DNA extraction and subsequent genome-wide methylation profiling on the Infinium Methylation EPIC Bead Chip microarray. Analysis of differential methylation revealed a clear pattern of discrimination between the wart and normal skin samples. In warts, the most differentially methylated (DM) genes included long non-coding RNAs (AC005884, AL049646.2, AC126121.2, AP001790.1, and AC107959.3), microRNAs (MIR374B, MIR596, MIR1255B1, MIR26B, and MIR196A2),snoRNAs (SNORD114-22, SNORD70, and SNORD114-31), pseudogenes (AC069366.1, RNU4ATAC11P, AC120057.1, NANOGP3, AC106038.2, TPT1P2, SDC4P, PKMP3, and VN2R3P), and protein-coding genes (AREG, GJB2, C12orf71, AC020909.2, S100A8, ZBED2, FABP7, and CYSLTR1). In addition, pathway analysis revealed that, among the most differentially methylated genes, STAT5A, RARA, MEF2D, MAP3K8, and THRA were the common regulators. It can be observed that HPV-induced warts involve a clear and unique epigenetic alteration to the host genome.

Bat-Borne Coronaviruses in Jordan and Saudi Arabia: A Threat to Public Health?

Emerging infectious diseases are of great concern to public health, as highlighted by the ongoing coronavirus disease 2019 (COVID-19) pandemic. Such diseases are of particular danger during mass gathering and mass influx events, as large crowds of people in close proximity to each other creates optimal opportunities for disease transmission. The Hashemite Kingdom of Jordan and the Kingdom of Saudi Arabia are two countries that have witnessed mass gatherings due to the arrival of Syrian refugees and the annual Hajj season. The mass migration of people not only brings exotic diseases to these regions but also brings new diseases back to their own countries, e.g., the outbreak of MERS in South Korea. Many emerging pathogens originate in bats, and more than 30 bat species have been identified in these two countries. Some of those bat species are known to carry viruses that cause deadly diseases in other parts of the world, such as the rabies virus and coronaviruses. However, little is known about bats and the pathogens they carry in Jordan and Saudi Arabia. Here, the importance of enhanced surveillance of bat-borne infections in Jordan and Saudi Arabia is emphasized, promoting the awareness of bat-borne diseases among the general public and building up infrastructure and capability to fill the gaps in public health preparedness to prevent future pandemics.

Genome-wide identification of methylated CpG sites in nongenital cutaneous warts

BACKGROUND

Low-risk HPV infection has not been the subject of epigenetic investigation. The present study was carried out in order to investigate the methylation status of CpG sites in non-genital cutaneous warts.

METHODS

Genomic DNA was extracted from 24 paired epidermal samples of warts and normal skin. DNA samples were bisulfite converted and underwent genome-wide methylation profiling using the Infinium MethylationEPIC BeadChip Kit.

RESULTS

From a total of 844,234 CpG sites, 56,960 and 43,040 CpG sites were found to be hypo- and hypermethylated, respectively, in non-genital cutaneous warts. The most differentially methylated CpG sites in warts were located within the C10orf26, FAM83H-AS1, ZNF644, LINC00702, GSAP, STAT5A, HDAC4, NCALD, and EXOC4 genes.

CONCLUSION

Non-genital cutaneous warts exhibit a unique CpG methylation signature.

HSF3 and Hsp70 Expression during Post-Hatch Cold Stress in Broiler Chickens Subjected to Embryonic Thermal Manipulation

Decades of selective breeding for commercial purposes have rendered the broiler chicken (Gallus gallus domesticus) highly susceptible to heat and cold stress. A multitude of studies have documented the effects of thermal manipulation (TM) on broiler thermotolerance during periods of post-hatch heat stress, but very few have focused on the effect of TM on a broiler’s ability to withstand cold stress. Therefore, the primary objective of the current study is to determine the effects of TM on the acquisition of thermotolerance in broilers via their expression of the stress-associated 70 kilodalton heat shock protein (Hsp70) gene and heat shock factor 3 (HSF3) gene. Briefly, Hubbard broiler embryos were subject to TM by increasing the incubation temperature to 39 °C and 65% relative humidity (RH) for 18 h daily, from embryonic days (ED) 10 to 18. Broilers were then exposed to cold stress by decreasing the room temperature to 16 °C during post-hatch days 32 to 37. After thermal challenge, broilers were euthanized and hepatic and splenic tissues were collected. Our results showed that TM decreased the hatchability rate and body temperature but improved the body weight gain. TM generally decreased the hepatic expression but did not change the splenic expression of HSF3 during cold stress. In contrast, both hepatic and splenic Hsp70 expression decreased during cold stress. The results of the present study may suggest that TM significantly affects a broiler’s genetic response to cold stress.

Embryonic Thermal Manipulation Affects the Antioxidant Response to Post-Hatch Thermal Exposure in Broiler Chickens

Simple Summary

The broiler chicken is one of the most important livestock species in the world, as it occupies a major role in the modern human diet. Due to uneven artificial selection pressures, the broiler has increased in size over the past few decades at the expense of its ability to withstand oxidative damage, the latter of which is often a byproduct of thermal stress. In order to attenuate the effects of heat stress, thermal manipulation (TM), which involves changes in incubation temperature at certain points of embryonic development, is increasingly being presented as a way in which to improve broiler thermotolerance. Therefore, the objective of this study was to investigate how TM might affect broiler response to post-hatch thermal stress in the context of the genes that help combat oxidative damage, namely the catalase, NADPH oxidase 4 (NOX4), and superoxide dismutase 2 (SOD2) genes. Expression of all three aforementioned genes differed significantly between TM and control chickens after exposure to cold and heat stress. Conclusively, TM may act as a viable mode of preventative treatment for broilers at risk of thermally induced oxidative stress.

Abstract

Thermal stress is a major source of oxidative damage in the broiler chicken (Gallus gallus domesticus) due to the latter’s impaired metabolic function. While heat stress has been extensively studied in broilers, the effects of cold stress on broiler physiologic and oxidative function are still relatively unknown. The present study aimed to understand how thermal manipulation (TM) might affect a broiler’s oxidative response to post-hatch thermal stress in terms of the mRNA expression of the catalase, NADPH oxidase 4 (NOX4), and superoxide dismutase 2 (SOD2) genes. During embryonic days 10 to 18, TM was carried out by raising the temperature to 39 °C at 65% relative humidity for 18 h/day. To induce heat stress, room temperature was raised from 21 to 35 °C during post-hatch days (PD) 28 to 35, while cold stress was induced during PD 32 to 37 by lowering the room temperature from 21 to 16 °C. At the end of the thermal stress periods, a number of chickens were euthanized to extract hepatic and splenic tissue from the heat-stressed group and cardiac, hepatic, muscular, and splenic tissue from the cold-stressed group. Catalase, NOX4, and SOD2 expression in the heart, liver, and spleen were decreased in TM chickens compared to controls after both cold and heat stress. In contrast, the expression levels of these genes in the breast muscles of the TM group were increased or not affected. Moreover, TM chicks possessed an increased body weight (BW) and decreased cloacal temperature (T^C) compared to controls on PD 37. In addition, TM led to increased BW and lower T^C after both cold and heat stress. Conclusively, our findings suggest that TM has a significant effect on the oxidative function of thermally stressed broilers.

Epigenome-wide analysis of common warts reveals aberrant promoter methylation

Epigenetic alteration of host DNA is a common occurrence in both low- and high-risk human papillomavirus (HPV) infection. Although changes in promoter methylation have been widely studied in HPV-associated cancers, they have not been the subject of much investigation in HPV-induced warts, which are a temporary manifestation of HPV infection. The present study sought to examine the differences in promoter methylation between warts and normal skin. To achieve this, DNA was extracted from 24 paired wart and normal skin samples and inputted into the Infinium MethylationEPIC BeadChip microarray. Differential methylation analysis revealed a clear pattern of hyper- and hypomethylation in warts compared to normal skin, and the most differentially methylated promoters were found within the EIF3EP2, CYSLTR1, C10orf99, KRT6B, LAMA4, and H3F3B genes as well as the C9orf30 pseudogene. Moreover, pathway analysis showed that the H3F3A, CDKN1A, and MAPK13 genes were the most common regulators among the most differentially methylated promoters. Since the tissue samples were excised from active warts, however, this differential methylation could either be a cellular response to HPV infection or an HPV-driven process to establish the wart and/or promote disease progression. Conclusively, it is apparent that HPV infection alters the methylation status of certain genes to possibly initiate the formation of a wart and maintain its presence.

Genome-Wide Tiling Array Analysis of HPV-Induced Warts Reveals Aberrant Methylation of Protein-Coding and Non-Coding Regions

The human papillomaviruses (HPV) are a group of double-stranded DNA viruses that exhibit an exclusive tropism for squamous epithelia. HPV can either be low- or high-risk depending on its ability to cause benign lesions or cancer, respectively. Unsurprisingly, the majority of epigenetic research has focused on the high-risk HPV types, neglecting the low-risk types in the process. Therefore, the main objective of this study is to better understand the epigenetics of wart formation by investigating the differences in methylation between HPV-induced cutaneous warts and normal skin. A number of clear and very significant differences in methylation patterns were found between cutaneous warts and normal skin. Around 55% of the top-ranking 100 differentially methylated genes in warts were protein coding, including the EXOC4, KCNU, RTN1, LGI1, IRF2, and NRG1 genes. Additionally, non-coding RNA genes, such as the AZIN1-AS1, LINC02008, and MGC27382 genes, constituted 11% of the top-ranking 100 differentially methylated genes. Warts exhibited a unique pattern of methylation that is a possible explanation for their transient nature. Since the genetics of cutaneous wart formation are not completely known, the findings of the present study could contribute to a better understanding of how HPV infection modulates host methylation to give rise to warts in the skin.

Gene Expression Profiling of MicroRNAs in HPV-Induced Warts and Normal Skin

: Infection with the human papillomavirus (HPV) is a common occurrence among the global population, with millions of new cases emerging on an annual basis. Dysregulated microRNA (miRNA) expression is increasingly being identified to play a role in a number of different diseases, especially in the context of high-risk HPV infection. The present study investigated the miRNA expression profiles of warts induced by low-risk HPV. In warts, miR-27b, miR-24-1, miR-3654, miR-647, and miR-1914 were downregulated while miR-612 was upregulated compared to normal skin. Using miRTargetLink Human, experimentally supported evidence was obtained showing that miR-27b targeted the vascular endothelial growth factor C (VEGFC) and CAMP-responsive element binding protein 1 (CREB1) genes. The VEGFC and CREB1 genes have been reported to be involved in tumorigenesis and wart formation, respectively. Similarly, the oxidized low-density lipoprotein receptor 1 (OLR1) gene, which plays an important role in the humoral immunity of the skin, and the plexin D1 (PLXND1) gene, which is highly expressed in tumor vasculature, were both found to be common targets of miR-27b, miR-1914, and miR-612.

Genome-Wide CpG Island Methylation Profiles of Cutaneous Skin with and without HPV Infection

HPV infection is one of the most commonly transmitted diseases among the global population. While it can be asymptomatic, non-genital HPV infection often gives rise to cutaneous warts, which are benign growths arising from the epidermal layer of the skin. This study aimed to produce a global analysis of the ways in which cutaneous wart formation affected the CpG island methylome. The Infinium MethylationEPIC BeadChip microarray was utilized in order to quantitatively interrogate CpG island methylation in genomic DNA extracted from 24 paired wart and normal skin samples. Differential methylation analysis was carried out by means of assigning a combined rank score using RnBeads. The 1000 top-ranking CpG islands were then subject to Locus Overlap Analysis (LOLA) for enrichment of genomic ranges, while signaling pathway analysis was carried out on the top 100 differentially methylated CpG islands. Differential methylation analysis illustrated that the most differentially methylated CpG islands in warts lay within the ITGB5, DTNB, RBFOX3, SLC6A9, and C2orf27A genes. In addition, the most enriched genomic region sets in warts were Sheffield's tissue-clustered DNase hypersensitive sites, ENCODE's segmentation and transcription factor binding sites, codex sites, and the epigenome sites from cistrome. Lastly, signaling pathway analysis showed that the GRB2, GNB1, NTRK1, AXIN1, and SKI genes were the most common regulators of the genes associated with the top 100 most differentially methylated CpG islands in warts. Our study shows that HPV-induced cutaneous warts have a clear CpG island methylation profile that sets them apart from normal skin. Such a finding could account for the temporary nature of warts and the capacity for individuals to undergo clinical remission.