Rationale: IL-33 is a proinflammatory cytokine thought to play a role in the pathogenesis of asthma and chronic obstructive pulmonary disease (COPD). A recent clinical trial using an anti-IL-33... Show moreRationale: IL-33 is a proinflammatory cytokine thought to play a role in the pathogenesis of asthma and chronic obstructive pulmonary disease (COPD). A recent clinical trial using an anti-IL-33 antibody showed a reduction in exacerbation and improved lung function in ex-smokers but not current smokers with COPD.Objectives: This study aimed to understand the effects of smoking status on IL-33.Methods: We investigated the association of smoking status with the level of gene expression of IL-33 in the airways in eight independent transcriptomic studies of lung airways. Additionally, we performed Western blot analysis and immunohistochemistry for IL-33 in lung tissue to assess protein levels.Measurements and Main Results: Across the bulk RNA-sequencing datasets, IL-33 gene expression and its signaling pathway were significantly lower in current versus former or never-smokers and increased upon smoking cessation (P < 0.05). Single-cell sequencing showed that IL-33 is predominantly expressed in resting basal epithelial cells and decreases during the differentiation process triggered by smoke exposure. We also found a higher transitioning of this cellular subpopulation into amore differentiated cell type during chronic smoking, potentially driving the reduction of IL-33. Protein analysis demonstrated lower IL-33 levels in lung tissue from current versus former smokers with COPD and a lower proportion of IL-33-positive basal cells in current versus ex-smoking controls.Conclusions: We provide strong evidence that cigarette smoke leads to an overall reduction in IL-33 expression in transcriptomic and protein level, and this may be due to the decrease in resting basal cells. Together, these findings may explain the clinical observation that a recent antibody-based anti-IL-33 treatment is more effective in former than current smokers with COPD. Show less
Cortical cell loss is a core feature of Huntington’s disease (HD), beginning many years before clinical motor diagnosis, during the premanifest stage. However, it is unclear how genetic topography... Show moreCortical cell loss is a core feature of Huntington’s disease (HD), beginning many years before clinical motor diagnosis, during the premanifest stage. However, it is unclear how genetic topography relates to cortical cell loss. Here, we explore the biological processes and cell types underlying this relationship and validate these using cell-specific post-mortem data.Eighty premanifest participants on average 15 years from disease onset and 71 controls were included. Using volumetric and diffusion MRI we extracted HD-specific whole brain maps where lower grey matter volume and higher grey matter mean diffusivity, relative to controls, were used as proxies of cortical cell loss. These maps were combined with gene expression data from the Allen Human Brain Atlas (AHBA) to investigate the biological processes relating genetic topography and cortical cell loss.Cortical cell loss was positively correlated with the expression of developmental genes (i.e. higher expression correlated with greater atrophy and increased diffusivity) and negatively correlated with the expression of synaptic and metabolic genes that have been implicated in neurodegeneration. These findings were consistent for diffusion MRI and volumetric HD-specific brain maps.As wild-type huntingtin is known to play a role in neurodevelopment, we explored the association between wild-type huntingtin (HTT) expression and developmental gene expression across the AHBA. Co-expression network analyses in 134 human brains free of neurodegenerative disorders were also performed. HTT expression was correlated with the expression of genes involved in neurodevelopment while co-expression network analyses also revealed that HTT expression was associated with developmental biological processes.Expression weighted cell-type enrichment (EWCE) analyses were used to explore which specific cell types were associated with HD cortical cell loss and these associations were validated using cell specific single nucleus RNAseq (snRNAseq) data from post-mortem HD brains.The developmental transcriptomic profile of cortical cell loss in preHD was enriched in astrocytes and endothelial cells, while the neurodegenerative transcriptomic profile was enriched for neuronal and microglial cells. Astrocyte-specific genes differentially expressed in HD post-mortem brains relative to controls using snRNAseq were enriched in the developmental transcriptomic profile, while neuronal and microglial-specific genes were enriched in the neurodegenerative transcriptomic profile.Our findings suggest that cortical cell loss in preHD may arise from dual pathological processes, emerging as a consequence of neurodevelopmental changes, at the beginning of life, followed by neurodegeneration in adulthood, targeting areas with reduced expression of synaptic and metabolic genes. These events result in age-related cell death across multiple brain cell types. Show less
Objectives: Autoantibody responses increase years before the onset of inflammatory arthritis (IA) and are stable during transitioning from clinically suspect arthralgia (CSA) to IA. Cytokine and... Show moreObjectives: Autoantibody responses increase years before the onset of inflammatory arthritis (IA) and are stable during transitioning from clinically suspect arthralgia (CSA) to IA. Cytokine and chemokine levels also increase years before IA onset. However, the course in the at-risk stage of CSA during progression to disease or non-progression is unknown. To increase the understanding of processes mediating disease development, we studied the course of cytokine, chemokine and related receptors gene expression in CSA patients during progression to IA and in CSA patients who ultimately did not develop IA. Methods: Whole-blood RNA expression of 37 inflammatory cytokines, chemokines and related receptors was determined by dual-colour reverse transcription multiplex ligation-dependent probe amplification in paired samples of CSA patients at CSA onset and either at IA development or after 24 months without IA development. ACPA-positive and ACPA-negative CSA patients developing IA were compared at CSA onset and during progression to IA. Generalised estimating equations tested changes over time. A false discovery rate approach was applied. Results: None of the cytokine/chemokine genes significantly changed in expression between CSA onset and IA development. In CSA patients without IA development, G-CSF expression decreased (P = 0.001), whereas CCR6 and TNIP1 expression increased (P < 0.001 and P = 0.002, respectively) over a 2 year period. Expression levels in ACPA-positive and ACPA-negative CSA patients who developed IA were similar. Conclusion: Whole-blood gene expression of assessed cytokines, chemokines and related receptors did not change significantly from CSA to IA development. This suggests that changes in expression of these molecules may not be related to the final process of developing chronicity and may have occurred preceding CSA onset. Changes in gene expression in CSA patients without IA development may provide clues for processes related to resolution. Show less
Patients with metabolic syndrome are often prescribed statins to prevent the development of cardiovascular disease. Conversely, data on their effects on non-alcoholic steatohepatitis (NASH) are... Show morePatients with metabolic syndrome are often prescribed statins to prevent the development of cardiovascular disease. Conversely, data on their effects on non-alcoholic steatohepatitis (NASH) are lacking. We evaluated these effects by feeding APOE*3-Leiden mice a Western-type diet (WTD) with or without atorvastatin to induce NASH and hepatic fibrosis. Besides the well-known plasma cholesterol lowering (−30%) and anti-atherogenic effects (severe lesion size −48%), atorvastatin significantly reduced hepatic steatosis (−22%), the number of aggregated inflammatory cells in the liver (−80%) and hepatic fibrosis (−92%) compared to WTD-fed mice. Furthermore, atorvastatin-treated mice showed less immunohistochemically stained areas of inflammation markers. Atorvastatin prevented accumulation of free cholesterol in the form of cholesterol crystals (−78%). Cholesterol crystals are potent inducers of the NLRP3 inflammasome pathway and atorvastatin prevented its activation, which resulted in reduced expression of the pro-inflammatory cytokines interleukin (IL)-1β (−61%) and IL-18 (−26%). Transcriptome analysis confirmed strong reducing effects of atorvastatin on inflammatory mediators, including NLRP3, NFκB and TLR4. The present study demonstrates that atorvastatin reduces hepatic steatosis, inflammation and fibrosis and prevents cholesterol crystal formation, thereby precluding NLRP3 inflammasome activation. This may render atorvastatin treatment as an attractive approach to reduce NAFLD and prevent progression into NASH in dyslipidemic patients. Show less
Background: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation.... Show moreBackground: Thyroid hormones play a key role in differentiation and metabolism and are known regulators of gene expression through both genomic and epigenetic processes including DNA methylation. The aim of this study was to examine associations between thyroid hormones and DNA methylation.Methods: We carried out a fixed-effect meta-analysis of epigenome-wide association study (EWAS) of blood DNA methylation sites from 8 cohorts from the ThyroidOmics Consortium, incorporating up to 7073 participants of both European and African ancestry, implementing a discovery and replication stage. Statistical analyses were conducted using normalized beta CpG values as dependent and log-transformed thyrotropin (TSH), free thyroxine, and free triiodothyronine levels, respectively, as independent variable in a linear model. The replicated findings were correlated with gene expression levels in whole blood and tested for causal influence of TSH and free thyroxine by two-sample Mendelian randomization (MR).Results: Epigenome-wide significant associations (p-value <1.1E-7) of three CpGs for free thyroxine, five for free triiodothyronine, and two for TSH concentrations were discovered and replicated (combined p-values = 1.5E-9 to 4.3E-28). The associations included CpG sites annotated to KLF9 (cg00049440) and DOT1L (cg04173586) that overlap with all three traits, consistent with hypothalamic-pituitary-thyroid axis physiology. Significant associations were also found for CpGs in FKBP5 for free thyroxine, and at CSNK1D/LINCO1970 and LRRC8D for free triiodothyronine. MR analyses supported a causal effect of thyroid status on DNA methylation of KLF9. DNA methylation of cg00049440 in KLF9 was inversely correlated with KLF9 gene expression in blood. The CpG at CSNK1D/LINC01970 overlapped with thyroid hormone receptor alpha binding peaks in liver cells. The total additive heritability of the methylation levels of the six significant CpG sites was between 25% and 57%. Significant methylation QTLs were identified for CpGs at KLF9, FKBP5, LRRC8D, and CSNK1D/LINC01970.Conclusions: We report novel associations between TSH, thyroid hormones, and blood-based DNA methylation. This study advances our understanding of thyroid hormone action particularly related to KLF9 and serves as a proof-of-concept that integrations of EWAS with other -omics data can provide a valuable tool for unraveling thyroid hormone signaling in humans by complementing and feeding classical in vitro and animal studies. Show less
The glucocorticoid stress hormones affect brain function via high-affinity mineralocorticoid receptors (MRs) and lower-affinity glucocorticoid receptors (GRs). MR and GR not only differ in affinity... Show moreThe glucocorticoid stress hormones affect brain function via high-affinity mineralocorticoid receptors (MRs) and lower-affinity glucocorticoid receptors (GRs). MR and GR not only differ in affinity for ligands, but also have distinct, sometimes opposite, actions on neuronal excitability and other cellular and higher-order parameters related to cerebral function. GR and MR messenger RNA (mRNA) levels are often used as a proxy for the responsiveness to glucocorticoids, assuming proportionality between mRNA and protein levels. This may be especially relevant for the MR, which because of its high affinity is already largely occupied at low basal (trough) hormone levels. Here we explore how GR and MR mRNA levels are associated with the expression of a shared target gene, glucocorticoid-induced leucine zipper (GILZ, coded by Tsc22d3) with basal and elevated levels of corticosterone in male mice, using in situ hybridization. Depending on the hippocampal subfield and the corticosterone levels, mRNA levels of MR rather than GR mostly correlated with GILZ mRNA in the hippocampus and hypothalamus at the bulk tissue level. At the individual cell level, these correlations were much weaker. Using publicly available single-cell RNA sequencing data, we again observed that MR and GR mRNA levels were only weakly correlated with target gene expression in glutamatergic and GABAergic neurons. We conclude that MR mRNA levels can be limiting for receptor action, but many other cell-specific and region-specific factors ultimately determine corticosteroid receptor action. Altogether, our results argue for caution while interpreting the consequences of changed receptor expression for the response to glucocorticoids. Show less
Roelands, J.; Ploeg, M. van der; Ijsselsteijn, M.E.; Dang, H.; Boonstra, J.J.; Hardwick, J.C.H.; ... ; Miranda, N.F.C.C. de 2022
ObjectiveBiological insights into the stepwise development and progression of colorectal cancer (CRC) are imperative to develop tailored approaches for early detection and optimal clinical... Show moreObjectiveBiological insights into the stepwise development and progression of colorectal cancer (CRC) are imperative to develop tailored approaches for early detection and optimal clinical management of this disease. Here, we aimed to dissect the transcriptional and immunologic alterations that accompany malignant transformation in CRC and to identify clinically relevant biomarkers through spatial profiling of pT1 CRC samples. DesignWe employed digital spatial profiling (GeoMx) on eight pT1 CRCs to study gene expression in the epithelial and stromal segments across regions of distinct histology, including normal mucosa, low-grade and high-grade dysplasia and cancer. Consecutive histology sections were profiled by imaging mass cytometry to reveal immune contextures. Finally, publicly available single-cell RNA-sequencing data was analysed to determine the cellular origin of relevant transcripts. ResultsComparison of gene expression between regions within pT1 CRC samples identified differentially expressed genes in the epithelium (n=1394 genes) and the stromal segments (n=1145 genes) across distinct histologies. Pathway analysis identified an early onset of inflammatory responses during malignant transformation, typified by upregulation of gene signatures such as innate immune sensing. We detected increased infiltration of myeloid cells and a shift in macrophage populations from pro-inflammatory HLA-DR(+)CD204(-) macrophages to HLA-DR(-)CD204(+) immune-suppressive subsets from normal tissue through dysplasia to cancer, accompanied by the upregulation of the CD47/SIRP alpha 'don't eat me signal'. ConclusionSpatial profiling revealed the molecular and immunological landscape of CRC tumourigenesis at early disease stage. We identified biomarkers with strong association with disease progression as well as targetable immune processes that are exploitable in a clinical setting. Show less
Rathnayake, S.N.H.; Ditz, B.; Nijnatten, J. van; Sadaf, T.; Hansbro, P.M.; Brandsma, C.A.; ... ; Faiz, A. 2022
Background and ObjectiveSmoking disturbs the bronchial-mucus-barrier. This study assesses the cellular composition and gene expression shifts of the bronchial-mucus-barrier with smoking to... Show moreBackground and ObjectiveSmoking disturbs the bronchial-mucus-barrier. This study assesses the cellular composition and gene expression shifts of the bronchial-mucus-barrier with smoking to understand the mechanism of mucosal damage by cigarette smoke exposure. We explore whether single-cell-RNA-sequencing (scRNA-seq) based cellular deconvolution (CD) can predict cell-type composition in RNA-seq data. MethodsRNA-seq data of bronchial biopsies from three cohorts were analysed using CD. The cohorts included 56 participants with chronic obstructive pulmonary disease [COPD] (38 smokers; 18 ex-smokers), 77 participants without COPD (40 never-smokers; 37 smokers) and 16 participants who stopped smoking for 1 year (11 COPD and 5 non-COPD-smokers). Differential gene expression was used to investigate gene expression shifts. The CD-derived goblet cell ratios were validated by correlating with staining-derived goblet cell ratios from the COPD cohort. Statistics were done in the R software (false discovery rate p-value < 0.05). ResultsBoth CD methods indicate a shift in bronchial-mucus-barrier cell composition towards goblet cells in COPD and non-COPD-smokers compared to ex- and never-smokers. It shows that the effect was reversible within a year of smoking cessation. A reduction of ciliated and basal cells was observed with current smoking, which resolved following smoking cessation. The expression of mucin and sodium channel (ENaC) genes, but not chloride channel genes, were altered in COPD and current smokers compared to never smokers or ex-smokers. The goblet cell-derived staining scores correlate with CD-derived goblet cell ratios. ConclusionSmoking alters bronchial-mucus-barrier cell composition, transcriptome and increases mucus production. This effect is partly reversible within a year of smoking cessation. CD methodology can predict goblet-cell percentages from RNA-seq. Show less
Thonglim, A.; Bortolami, G.; Delzon, S.; Larter, M.; Offringa, R.; Keurentjes, J.J.B.; ... ; Lens, F. 2022
The synergy between drought-responsive traits across different organs is crucial in the whole-plant mechanism influencing drought resilience. These organ interactions, however, are poorly... Show moreThe synergy between drought-responsive traits across different organs is crucial in the whole-plant mechanism influencing drought resilience. These organ interactions, however, are poorly understood, limiting our understanding of drought response strategies at the whole-plant level. Therefore, we need more integrative studies, especially on herbaceous species that represent many important food crops but remain underexplored in their drought response. We investigated inflorescence stems and rosette leaves of six Arabidopsis thaliana genotypes with contrasting drought tolerance, and combined anatomical observations with hydraulic measurements and gene expression studies to assess differences in drought response. The soc1ful double mutant was the most drought-tolerant genotype based on its synergistic combination of low stomatal conductance, largest stomatal safety margin, more stable leaf water potential during non-watering, reduced transcript levels of drought stress marker genes, and reduced loss of chlorophyll content in leaves, in combination with stems showing the highest embolism resistance, most pronounced lignification, and thickest intervessel pit membranes. In contrast, the most sensitive Cvi ecotype shows the opposite extreme of the same set of traits. The remaining four genotypes show variations in this drought syndrome. Our results reveal that anatomical, ecophysiological, and molecular adaptations across organs are intertwined, and multiple (differentially combined) strategies can be applied to acquire a certain level of drought tolerance. Show less
Context: South Asian individuals are more prone to develop type 2 diabetes (T2D) coinciding with earlier complications than Europids. While inflammation plays a central role in the development and... Show moreContext: South Asian individuals are more prone to develop type 2 diabetes (T2D) coinciding with earlier complications than Europids. While inflammation plays a central role in the development and progression of T2D, this factor is still underexplored in South Asians. Objective: This work aimed to study whether circulating messenger RNA (mRNA) transcripts of immune genes are different between South Asian compared with Europid patients with T2D. Methods: A secondary analysis was conducted of 2 randomized controlled trials of Dutch South Asian (n = 45; age: 55 +/- 10 years, body mass index [BMI]: 29 +/- 4 kg/m(2)) and Dutch Europid (n = 44; age: 60 +/- 7 years, BMI: 32 +/- 4 kg/m(2)) patients with T2D. Main outcome measures included mRNA transcripts of 182 immune genes (microfluidic quantitative polymerase chain reaction; Fluidigm Inc) in fasted whole-blood, ingenuity pathway analyses (Qiagen). Results: South Asians, compared to Europids, had higher mRNA levels of B-cell markers (CD19, CD79A, CD79B, CR2, CXCR5, IGHD, MS4A1, PAX5; all fold change > 1.3, false discovery rate [FDR] < 0.008) and interferon (IFN)-signaling genes (CD274, GBP1, GBP2, GBP5, FCGR1A/B/CP, IFI16, IFIT3, IFITM1, IFITM3, TAP1; all FC > 1.2, FDR < 0.05). In South Asians, the IFN signaling pathway was the top canonical pathway (z score 2.6; P < .001) and this was accompanied by higher plasma IFN-gamma levels (FC = 1.5, FDR = 0.01). Notably, the ethnic difference in gene expression was larger for women (20/182 [11%]) than men (2/182 [1%]). Conclusion: South Asian patients with T2D show a more activated IFN-signaling pathway compared to Europid patients with T2D, which is more pronounced in women than men. We speculate that a more activated IFN-signaling pathway may contribute to the more rapid progression of T2D in South Asian compared with Europid individuals. Show less
Gabel, G.; Northoff, B.H.; Balboa, A.; Becirovic-Agic, M.; Petri, M.; Busch, A.; ... ; Lindeman, J.H.N. 2021
Background While numerous interventions effectively interfered with abdominal aortic aneurysm (AAA) formation/progression in preclinical models, none of the successes translated into clinical... Show moreBackground While numerous interventions effectively interfered with abdominal aortic aneurysm (AAA) formation/progression in preclinical models, none of the successes translated into clinical success. Hence, a systematic exploration of parallel and divergent processes in clinical AAA disease and its 2 primary models (the porcine pancreatic elastase and angiotensin-II infusion [AngII] murine model) was performed to identify mechanisms relevant for aneurysm disease. Methods and Results This study combines Movat staining and pathway analysis for histological and genomic comparisons between clinical disease and its models. The impact of a notable genomic signal for metabolic reprogramming was tested in a rescue trial (AngII model) evaluating the impact of 1-(4-pyridinyl)-3-(2-quinolinyl)-2-propen-1-one (PFK15)-mediated interference with main glycolytic switch PFKFB3. Histological evaluation characterized the AngII model as a dissection model that is accompanied by adventitial fibrosis. The porcine pancreatic elastase model showed a transient inflammatory response and aortic dilatation, followed by stabilization and fibrosis. Normalization of the genomic responses at day 14 confirmed the self-limiting nature of the porcine pancreatic elastase model. Clear parallel genomic responses with activated adaptive immune responses, and particularly strong signals for metabolic switching were observed in human AAA and the AngII model. Rescue intervention with the glycolysis inhibitor PFK15 in the AngII model showed that interference with the glycolytic switching quenches aneurysm formation. Conclusions Despite clear morphological contrasts, remarkable genomic parallels exist for clinical AAA disease and the AngII model. The metabolic response appears causatively involved in AAA progression and provides a novel therapeutic target. The clear transient genomic response classifies the porcine pancreatic elastase model as a disease initiation model. Show less
Gliddon, H.D.; Kaforou, M.; Alikian, M.; Habgood-Coote, D.; Zhou, C.X.; Oni, T.; ... ; ILULU Consortium 2021
Recently, host whole blood gene expression signatures have been identified for diagnosis of tuberculosis (TB). Absolute quantification of the concentrations of signature transcripts in blood have... Show moreRecently, host whole blood gene expression signatures have been identified for diagnosis of tuberculosis (TB). Absolute quantification of the concentrations of signature transcripts in blood have not been reported, but would facilitate diagnostic test development. To identify minimal transcript signatures, we applied a transcript selection procedure to microarray data from African adults comprising 536 patients with TB, other diseases (OD) and latent TB (LTBI), divided into training and test sets. Signatures were further investigated using reverse transcriptase (RT)-digital PCR (dPCR). A four-transcript signature (GBP6, TMCC1, PRDM1, and ARG1) measured using RT-dPCR distinguished TB patients from those with OD (area under the curve (AUC) 93.8% (CI95% 82.2-100%). A three-transcript signature (FCGR1A, ZNF296, and C1QB) differentiated TB from LTBI (AUC 97.3%, CI95%: 93.3-100%), regardless of HIV. These signatures have been validated across platforms and across samples offering strong, quantitative support for their use as diagnostic biomarkers for TB. Show less
One of the essential functions of microglia is to continuously sense changes in their environment and adapt to those changes. For this purpose, they use a set of genes termed the sensome. This... Show moreOne of the essential functions of microglia is to continuously sense changes in their environment and adapt to those changes. For this purpose, they use a set of genes termed the sensome. This sensome is comprised of the most abundantly expressed receptors on the surface of microglia. In this study, we updated previously identified mouse microglial sensome by incorporating an additional published RNAseq dataset into the data-analysis pipeline. We also identified members of the human microglial sensome using two independent human microglia RNAseq data sources. Using both the mouse and human microglia sensomes, we identified a key set of genes conserved between the mouse and human microglial sensomes as well as some differences between the species. We found a key set of 57 genes to be conserved in both mouse and human microglial sensomes. We define these genes as the "microglia core sensome". We then analyzed expression of genes in this core sensome in five different datasets from two neurodegenerative disease models at various stages of the diseases and found that, overall, changes in the level of expression of microglial sensome genes are specific to the disease or condition studied. Our results highlight the relevance of data generated in mice for understanding the biology of human microglia, but also stress the importance of species-specific gene sets for the investigation of diseases involving microglia. Defining this microglial specific core sensome may help identify pathological changes in microglia in humans and mouse models of human disease. Show less
Lemmers, R.J.L.F.; Vliet, P.J. van der; Blatnik, A.; Balog, J.; Zidar, J.; Henderson, D.; ... ; Maarel, S.M. van der 2021
Background: Facioscapulohumeral dystrophy (FSHD) is an inherited muscular dystrophy clinically characterised by muscle weakness starting with the facial and upper extremity muscles. A disease model... Show moreBackground: Facioscapulohumeral dystrophy (FSHD) is an inherited muscular dystrophy clinically characterised by muscle weakness starting with the facial and upper extremity muscles. A disease model has been developed that postulates that failure in somatic repression of the transcription factor DUX4 embedded in the D4Z4 repeat on chromosome 4q causes FSHD. However, due to the position of the D4Z4 repeat close to the telomere and the complex genetic and epigenetic aetiology of FSHD, there is ongoing debate about the transcriptional deregulation of closely linked genes and their involvement in FSHD. Method: Detailed genetic characterisation and gene expression analysis of patients with clinically confirmed FSHD and control individuals. Results: Identification of two FSHD families in which the disease is caused by repeat contraction and DUX4 expression from chromosome 10 due to a de novo D4Z4 repeat exchange between chromosomes 4 and 10. We show that the genetic lesion causal to FSHD in these families is physically separated from other candidate genes on chromosome 4. We demonstrate that muscle cell cultures from affected family members exhibit the characteristic molecular features of FSHD, including DUX4 and DUX4 target gene expression, without showing evidence for transcriptional deregulation of other chromosome 4-specific candidate genes. Conclusion: This study shows that in rare situations, FSHD can occur on chromosome 10 due to an interchromosomal rearrangement with the FSHD locus on chromosome 4q. These findings provide further evidence that DUX4 derepression is the dominant disease pathway for FSHD. Hence, therapeutic strategies should focus on DUX4 as the primary target. Show less
Gray, L.G.; Mills, J.D.; Curry-Hyde, A.; Devore, S.; Friedman, D.; Thom, M.; ... ; Janitz, M. 2020
Circular RNAs (circRNAs) regulate mRNA translation by binding to microRNAs (miRNAs), and their expression is altered in diverse disorders, including cancer, cardiovascular disease, and Parkinson's... Show moreCircular RNAs (circRNAs) regulate mRNA translation by binding to microRNAs (miRNAs), and their expression is altered in diverse disorders, including cancer, cardiovascular disease, and Parkinson's disease. Here, we compare circRNA expression patterns in the temporal cortex and hippocampus of patients with pharmacoresistant mesial temporal lobe epilepsy (MTLE) and healthy controls. Nine circRNAs showed significant differential expression, including circRNA-HOMER1, which is expressed in synapses. Further, we identified miRNA binding sites within the sequences of differentially expressed (DE) circRNAs; expression levels of mRNAs correlated with changes in complementary miRNAs. Gene set enrichment analysis of mRNA targets revealed functions in heterocyclic compound binding, regulation of transcription, and signal transduction, which maintain the structure and function of hippocampal neurons. The circRNA-miRNA-mRNA interaction networks illuminate the molecular changes in MTLE, which may be pathogenic or an effect of the disease or treatments and suggests that DE circRNAs and associated miRNAs may be novel therapeutic targets. Show less
Mechanical stress determines bone mass and structure. It is not known whether mechanical loading affects expression of bone regulatory genes in a combined deficiency of estrogen and vitamin D. We... Show moreMechanical stress determines bone mass and structure. It is not known whether mechanical loading affects expression of bone regulatory genes in a combined deficiency of estrogen and vitamin D. We studied the effect of mechanical loading on the messenger RNA (mRNA) expression of bone regulatory genes during vitamin D and/or estrogen deficiency. We performed a single bout in vivo axial loading with 14 N peak load, 2 Hz frequency and 360 cycles in right ulnae of nineteen weeks old female control Wistar rats with or without ovariectomy (OVX), vitamin D deficiency and the combination of OVX and vitamin D deficiency (N = 10/group). Total bone RNA was isolated 6 hours after loading, and mRNA expression was detected ofMepe, Fgf23, Dmp1, Phex, Sost, Col1a1, Cyp27b1, Vdr, andEsr1. Serum levels of 25(OH)D, 1,25(OH)(2)D and estradiol were also measured at this time point. The effect of loading, vitamin D and estrogen deficiency and their interaction on bone gene expression was tested using a mixed effect model analysis. Mechanical loading significantly increased the mRNA expression ofMepe, andSost, whereas it decreased the mRNA expression ofFgf23andEsr1. Mechanical loading showed a significant interaction with vitamin D deficiency with regard to mRNA expression ofVdrandEsr1. Mechanical loading affected gene expression ofMepe, Fgf23, Sost, andEsr1independently of vitamin D or estrogen, indicating that mechanical loading may affect bone turnover even during vitamin D deficiency and after menopause. Show less
Ouwerkerk, A.F. van; Bosada, F.M.; Liu, J.; Zhang, J.; Duijvenboden, K. van; Chaffin, M.; ... ; Christoffels, V.M. 2020
Rationale: Genome-wide association studies have identified a large number of common variants (single-nucleotide polymorphisms) associated with atrial fibrillation (AF). These variants are located... Show moreRationale: Genome-wide association studies have identified a large number of common variants (single-nucleotide polymorphisms) associated with atrial fibrillation (AF). These variants are located mainly in noncoding regions of the genome and likely include variants that modulate the function of transcriptional regulatory elements (REs) such as enhancers. However, the actual REs modulated by variants and the target genes of such REs remain to be identified. Thus, the biological mechanisms by which genetic variation promotes AF has thus far remained largely unexplored. Objective: To identify REs in genome-wide association study loci that are influenced by AF-associated variants. Methods and Results: We screened 2.45 Mbp of human genomic DNA containing 12 strongly AF-associated loci for RE activity using self-transcribing active regulatory region sequencing and a recently generated monoclonal line of conditionally immortalized rat atrial myocytes. We identified 444 potential REs, 55 of which contain AF-associated variants (P<10(-8)). Subsequently, using an adaptation of the self-transcribing active regulatory region sequencing approach, we identified 24 variant REs with allele-specific regulatory activity. By mining available chromatin conformation data, the possible target genes of these REs were mapped. To define the physiological function and target genes of such REs, we deleted the orthologue of an RE containing noncoding variants in theHcn4(potassium/sodium hyperpolarization-activated cyclic nucleotide-gated channel 4) locus of the mouse genome. Mice heterozygous for the RE deletion showed bradycardia, sinus node dysfunction, and selective loss ofHcn4expression. Conclusions: We have identified REs at multiple genetic loci for AF and found that loss of an RE at theHCN4locus results in sinus node dysfunction and reduced gene expression. Our approach can be broadly applied to facilitate the identification of human disease-relevant REs and target genes at cardiovascular genome-wide association studies loci. Show less
George, L.; Taylor, A.R.; Esteve-Codina, A.; Artigas, M.S.; Thun, G.A.; Bates, S.; ... ; EvA Study Team 2020
Background: Whether the clinical or pathophysiologic significance of the "treatable trait" high blood eosinophil count in COPD is the same as for asthma remains controversial. We sought to... Show moreBackground: Whether the clinical or pathophysiologic significance of the "treatable trait" high blood eosinophil count in COPD is the same as for asthma remains controversial. We sought to determine the relationship between the blood eosinophil count, clinical characteristics and gene expression from bronchial brushings in COPD and asthma.Methods: Subjects were recruited into a COPD (emphysema versus airway disease [EvA]) or asthma cohort (Unbiased BIOmarkers in PREDiction of respiratory disease outcomes, U-BIOPRED). We determined gene expression using RNAseq in EvA (n = 283) and Affymetrix microarrays in U-BIOPRED (n = 85). We ran linear regression analysis of the bronchial brushings transcriptional signal versus blood eosinophil counts as well as differential expression using a blood eosinophil > 200 cells/mu L as a cut-off. The false discovery rate was controlled at 1% (with continuous values) and 5% (with dichotomized values).Results: There were no differences in age, gender, lung function, exercise capacity and quantitative computed tomography between eosinophilic versus noneosinophilic COPD cases. Total serum IgE was increased in eosinophilic asthma and COPD. In EvA, there were 12 genes with a statistically significant positive association with the linear blood eosinophil count, whereas in U-BIOPRED, 1197 genes showed significant associations (266 positive and 931 negative). The transcriptome showed little overlap between genes and pathways associated with blood eosinophil counts in asthma versus COPD. Only CST1 was common to eosinophilic asthma and COPD and was replicated in independent cohorts.Conclusion: Despite shared "treatable traits" between asthma and COPD, the molecular mechanisms underlying these clinical entities are predominately different. Show less
Chatterjee, S.; Kang, S.D.; Alam, S.; Salzberg, A.C.; Milici, J.; Burg, S.H. van der; ... ; Meyers, C. 2019
Epidemiological data confirm a much higher incidence of high-risk human papillomavirus 16 (HPV16)-mediated carcinogenesis of the cervical epithelium than for other target sites. In order to... Show moreEpidemiological data confirm a much higher incidence of high-risk human papillomavirus 16 (HPV16)-mediated carcinogenesis of the cervical epithelium than for other target sites. In order to elucidate tissue-specific responses to virus infection, we compared gene expression changes induced by productive HPV16 infection of cervical, foreskin, and tonsil organotypic rafts. These rafts closely mimic persistent HPV16 infection, long before carcinogenesis sets in. The total number of gene expression changes varied considerably across the tissue types, with only 32 genes being regulated in common. Among them, we confirmed the Kelch-like family protein KLHL35 and the laminin-5 complex to be upregulated and downregulated, respectively, in all the three tissues. HPV16 infection induces upregulation of genes involved in cell cycle control, cell division, mitosis, DNA replication, and DNA damage repair in all the three tissues, indicative of a hyperproliferative environment. In the cervical and tonsil epithelium, we observe significant downregulation of genes involved in epidermis development, keratinocyte differentiation, and extracellular matrix organization. On the other hand, in HPV16-positive foreskin (HPV16 foreskin) tissue, several genes involved in interferon-mediated innate immunity, cytokine signaling, and cellular defenses were downregulated. Furthermore, pathway analysis and experimental validations identified important cellular pathways like STAT1 and transforming growth factor beta (TGF-beta) to be differentially regulated among the three tissue types. The differential modulation of important cellular pathways like TGF-beta 1 and STAT1 can explain the sensitivity of tissues to HPV cancer progression.IMPORTANCE Although the high-risk human papillomavirus 16 infects anogenital and oropharyngeal sites, the cervical epithelium has a unique vulnerability to progression of cancer. Host responses during persistent infection and preneoplastic stages can shape the outcome of cancer progression in a tissue-dependent manner. Our study for the first time reports differential regulation of critical cellular functions and signaling pathways during productive HPV16 infection of cervical, foreskin, and tonsil tissues. While the virus induces hyperproliferation in infected cells, it downregulates epithelial differentiation, epidermal development, and innate immune responses, according to the tissue type. Modulation of these biological functions can determine virus fitness and pathogenesis and illuminate key cellular mechanisms that the virus employs to establish persistence and finally initiate disease progression. Show less