Corrigendum: Identifying predictors of clinical outcomes using the projection-predictive feature selection-a proof of concept on the example of Crohn's disease.

[This corrects the article DOI: 10.3389/fped.2023.1170563.].

Identifying predictors of clinical outcomes using the projection-predictive feature selection-a proof of concept on the example of Crohn's disease.

Objectives: Several clinical disease activity indices (DAIs) have been developed to noninvasively assess mucosal healing in pediatric Crohn's disease (CD). However, their clinical application can be complex. Therefore, we present a new way to identify the most informative biomarkers for mucosal inflammation from current markers in use and, based on this, how to obtain an easy-to-use DAI for clinical practice. A further aim of our proof-of-concept study is to demonstrate how the performance of such a new DAI can be compared to that of existing DAIs. Methods: The data of two independent study cohorts, with 167 visits from 109 children and adolescents with CD, were evaluated retrospectively. A variable selection based on a Bayesian ordinal regression model was applied to select clinical or standard laboratory parameters as predictors, using an endoscopic outcome. The predictive performance of the resulting model was compared to that of existing pediatric DAIs. Results: With our proof-of-concept dataset, the resulting model included C-reactive protein (CRP) and fecal calprotectin (FC) as predictors. In general, our model performed better than the existing DAIs. To show how our Bayesian approach can be applied in practice, we developed a web application for predicting disease activity for a new CD patient or visit. Conclusions: Our work serves as a proof-of-concept, showing that the statistical methods used here can identify biomarkers relevant for the prediction of a clinical outcome. In our case, a small number of biomarkers is sufficient, which, together with the web interface, facilitates the clinical application. However, the retrospective nature of our study, the rather small amount of data, and the lack of an external validation cohort do not allow us to consider our results as the establishment of a novel DAI for pediatric CD. This needs to be done with the help of a prospective study with more data and an external validation cohort in the future.

Large Duodenal Hematoma Causing an Ileus after an Endoscopic Duodenal Biopsy in a 6-Year-Old Child: A Case Report.

Intramural duodenal hematoma (IDH) in children is a rare complication after esophagogastroduodenoscopy. It is commonly described in patients with additional disorders or risk factors, such as coagulopathy. We present a case of a previously healthy 6-year-old boy with a large obstructing intramural duodenal hematoma and concomitant pancreatitis after an elective esophagogastroduodenoscopy. The patient presented with typical symptoms of an IDH, such as abdominal pain and distension, nausea and vomiting. IDH was diagnosed using ultrasound and magnetic resonance imaging examination. Conservative management with gastric decompression using a nasogastric feeding tube, bowel rest, total parenteral nutrition and analgesia was performed. After three weeks, the patient was discharged from the hospital without any complaints. Interventional management of IDH in pediatric patients with a lack of response to conservative therapy or complicating IDH should be discussed in an interdisciplinary team.

StreaMD: Advanced analysis of molecular dynamics using R.

Gromacs is one of the most popular molecular simulation suites currently available. In this contribution we present streaMD, the first interface between Gromacs trajectory files and the statistical language R. The amount of data created due to ever increasing computational power renders fast and efficient analysis of trajectories into a challenge. Especially as standard approaches such as root-mean square fluctuations and the like provide only limited physical insight. In our streaMD package integration of the Gromacs I/O libraries with advanced, graph-based analysis methods as the java library Stream leads to both: improved speed and analysis depth. We benchmark our results and highlight the applicability of the package by an interesting problem in RNA design, namely the interaction of tetracycline with an aptamer. © 2018 Wiley Periodicals, Inc.

StreAM-[Formula: see text]: algorithms for analyzing coarse grained RNA dynamics based on Markov models of connectivity-graphs.

BACKGROUND: In this work, we present a new coarse grained representation of RNA dynamics. It is based on adjacency matrices and their interactions patterns obtained from molecular dynamics simulations. RNA molecules are well-suited for this representation due to their composition which is mainly modular and assessable by the secondary structure alone. These interactions can be represented as adjacency matrices of k nucleotides. Based on those, we define transitions between states as changes in the adjacency matrices which form Markovian dynamics. The intense computational demand for deriving the transition probability matrices prompted us to develop StreAM-[Formula: see text], a stream-based algorithm for generating such Markov models of k-vertex adjacency matrices representing the RNA. RESULTS: We benchmark StreAM-[Formula: see text] (a) for random and RNA unit sphere dynamic graphs (b) for the robustness of our method against different parameters. Moreover, we address a riboswitch design problem by applying StreAM-[Formula: see text] on six long term molecular dynamics simulation of a synthetic tetracycline dependent riboswitch (500 ns) in combination with five different antibiotics. CONCLUSIONS: The proposed algorithm performs well on large simulated as well as real world dynamic graphs. Additionally, StreAM-[Formula: see text] provides insights into nucleotide based RNA dynamics in comparison to conventional metrics like the root-mean square fluctuation. In the light of experimental data our results show important design opportunities for the riboswitch.

Compile- and run-time approaches for the selection of efficient data structures for dynamic graph analysis.

Graphs are used to model a wide range of systems from different disciplines including social network analysis, biology, and big data processing. When analyzing these constantly changing dynamic graphs at a high frequency, performance is the main concern. Depending on the graph size and structure, update frequency, and read accesses of the analysis, the use of different data structures can yield great performance variations. Even for expert programmers, it is not always obvious, which data structure is the best choice for a given scenario. In previous work, we presented an approach for handling the selection of the most efficient data structures automatically using a compile-time approach well-suited for constant workloads. We extend this work with a measurement study of seven data structures and use the results to fit actual cost estimation functions. In addition, we evaluate our approach for the computations of seven different graph metrics. In analyses of real-world dynamic graphs with a constant workload, our approach achieves a speedup of up to 5.4× compared to basic data structure configurations. Such a compile-time based approach cannot yield optimal results when the behavior of the system changes later and the workload becomes non-constant. To close this gap we present a run-time approach which provides live profiling and facilitates automatic exchanges of data structures during execution. We analyze the performance of this approach using an artificial, non-constant workload where our approach achieves speedups of up to 7.3× compared to basic configurations.

Glucocorticoid receptor binds half sites as a monomer and regulates specific target genes.

BACKGROUND: Glucocorticoid receptor (GR) is a hormone-activated, DNA-binding transcriptional regulatory factor that controls inflammation, metabolism, stress responses, and other physiological processes. In vitro, GR binds as an inverted dimer to a motif consisting of two imperfectly palindromic 6 bp half sites separated by 3 bp spacers. In vivo, GR employs different patterns of functional surfaces of GR to regulate different target genes. The relationships between GR genomic binding and functional surface utilization have not been defined. RESULTS: We find that A477T, a GR mutant that disrupts the dimerization interface, differs from wild-type GRα in binding and regulation of target genes. Genomic regions strongly occupied by A477T are enriched for a novel half site motif. In vitro, GRα binds half sites as a monomer. Through the overlap between GRα- and A477T-bound regions, we identify GRα-bound regions containing only half sites. We further identify GR target genes linked with half sites and not with the full motif. CONCLUSIONS: Genomic regions bound by GR differ in underlying DNA sequence motifs and in the GR functional surfaces employed for regulation. Identification of GR binding regions that selectively utilize particular GR surfaces may discriminate sub-motifs, including the half site motif, that favor those surfaces. This approach may contribute to predictive models for GR activity and therapy.

Hic-5 is a transcription coregulator that acts before and/or after glucocorticoid receptor genome occupancy in a gene-selective manner.

Ligand activation and DNA-binding dictate the outcome of glucocorticoid receptor (GR)-mediated transcriptional regulation by inducing diverse receptor conformations that interact differentially with coregulators. GR recruits many coregulators via the well-characterized AF2 interaction surface in the GR ligand-binding domain, but Lin11, Isl-1, Mec-3 (LIM) domain coregulator Hic-5 (TGFB1I1) binds to the relatively uncharacterized tau2 activation domain in the hinge region of GR. Requirement of hydrogen peroxide-inducible clone-5 (Hic-5) for glucocorticoid-regulated gene expression was defined by Hic-5 depletion and global gene-expression analysis. Hic-5 depletion selectively affected both activation and repression of GR target genes, and Hic-5 served as an on/off switch for glucocorticoid regulation of many genes. For some hormone-induced genes, Hic-5 facilitated recruitment of Mediator complex. In contrast, many genes were not regulated by glucocorticoid until Hic-5 was depleted. On these genes Hic-5 prevented GR occupancy and chromatin remodeling and thereby inhibited their hormone-dependent regulation. Transcription factor binding to genomic sites is highly variable among different cell types; Hic-5 represents an alternative mechanism for regulating transcription factor-binding site selection that could apply both within a given cell type and among different cell types. Thus, Hic-5 is a versatile coregulator that acts by multiple gene-specific mechanisms that influence genomic occupancy of GR as well transcription complex assembly.

Defects in the C. elegans acyl-CoA synthase, acs-3, and nuclear hormone receptor, nhr-25, cause sensitivity to distinct, but overlapping stresses.

Metazoan transcription factors control distinct networks of genes in specific tissues, yet understanding how these networks are integrated into physiology, development, and homeostasis remains challenging. Inactivation of the nuclear hormone receptor nhr-25 ameliorates developmental and metabolic phenotypes associated with loss of function of an acyl-CoA synthetase gene, acs-3. ACS-3 activity prevents aberrantly high NHR-25 activity. Here, we investigated this relationship further by examining gene expression patterns following acs-3 and nhr-25 inactivation. Unexpectedly, we found that the acs-3 mutation or nhr-25 RNAi resulted in similar transcriptomes with enrichment in innate immunity and stress response gene expression. Mutants of either gene exhibited distinct sensitivities to pathogens and environmental stresses. Only nhr-25 was required for wild-type levels of resistance to the bacterial pathogen P. aeruginosa and only acs-3 was required for wild-type levels of resistance to osmotic stress and the oxidative stress generator, juglone. Inactivation of either acs-3 or nhr-25 compromised lifespan and resistance to the fungal pathogen D. coniospora. Double mutants exhibited more severe defects in the lifespan and P. aeruginosa assays, but were similar to the single mutants in other assays. Finally, acs-3 mutants displayed defects in their epidermal surface barrier, potentially accounting for the observed sensitivities. Together, these data indicate that inactivation of either acs-3 or nhr-25 causes stress sensitivity and increased expression of innate immunity/stress genes, most likely by different mechanisms. Elevated expression of these immune/stress genes appears to abrogate the transcriptional signatures relevant to metabolism and development.

The glucocorticoid receptor dimer interface allosterically transmits sequence-specific DNA signals.

Glucocorticoid receptor (GR) binds to genomic response elements and regulates gene transcription with cell and gene specificity. Within a response element, the precise sequence to which the receptor binds has been implicated in directing its structure and activity. Here, we use NMR chemical-shift difference mapping to show that nonspecific interactions with bases at particular positions in the binding sequence, such as those of the 'spacer', affect the conformation of distinct regions of the rat GR DNA-binding domain. These regions include the DNA-binding surface, the 'lever arm' and the dimerization interface, suggesting an allosteric pathway that signals between the DNA-binding sequence and the associated dimer partner. Disrupting this pathway by mutating the dimer interface alters sequence-specific conformations, DNA-binding kinetics and transcriptional activity. Our study demonstrates that GR dimer partners collaborate to read DNA shape and to direct sequence-specific gene activity.

Stalled spliceosomes are a signal for RNAi-mediated genome defense.

Using the yeast Cryptococcus neoformans, we describe a mechanism by which transposons are initially targeted for RNAi-mediated genome defense. We show that intron-containing mRNA precursors template siRNA synthesis. We identify a Spliceosome-Coupled And Nuclear RNAi (SCANR) complex required for siRNA synthesis and demonstrate that it physically associates with the spliceosome. We find that RNAi target transcripts are distinguished by suboptimal introns and abnormally high occupancy on spliceosomes. Functional investigations demonstrate that the stalling of mRNA precursors on spliceosomes is required for siRNA accumulation. Lariat debranching enzyme is also necessary for siRNA production, suggesting a requirement for processing of stalled splicing intermediates. We propose that recognition of mRNA precursors by the SCANR complex is in kinetic competition with splicing, thereby promoting siRNA production from transposon transcripts stalled on spliceosomes. Disparity in the strength of expression signals encoded by transposons versus host genes offers an avenue for the evolution of genome defense.

Disruption of the abdominal-B promoter tethering element results in a loss of long-range enhancer-directed Hox gene expression in Drosophila.

There are many examples within gene complexes of transcriptional enhancers interacting with only a subset of target promoters. A number of molecular mechanisms including promoter competition, insulators and chromatin looping are thought to play a role in regulating these interactions. At the Drosophila bithorax complex (BX-C), the IAB5 enhancer specifically drives gene expression only from the Abdominal-B (Abd-B) promoter, even though the enhancer and promoter are 55 kb apart and are separated by at least three insulators. In previous studies, we discovered that a 255 bp cis-regulatory module, the promoter tethering element (PTE), located 5' of the Abd-B transcriptional start site is able to tether IAB5 to the Abd-B promoter in transgenic embryo assays. In this study we examine the functional role of the PTE at the endogenous BX-C using transposon-mediated mutagenesis. Disruption of the PTE by P element insertion results in a loss of enhancer-directed Abd-B expression during embryonic development and a homeotic transformation of abdominal segments. A partial deletion of the PTE and neighboring upstream genomic sequences by imprecise excision of the P element also results in a similar loss of Abd-B expression in embryos. These results demonstrate that the PTE is an essential component of the regulatory network at the BX-C and is required in vivo to mediate specific long-range enhancer-promoter interactions.

Functional evolution of cis-regulatory modules at a homeotic gene in Drosophila.

It is a long-held belief in evolutionary biology that the rate of molecular evolution for a given DNA sequence is inversely related to the level of functional constraint. This belief holds true for the protein-coding homeotic (Hox) genes originally discovered in Drosophila melanogaster. Expression of the Hox genes in Drosophila embryos is essential for body patterning and is controlled by an extensive array of cis-regulatory modules (CRMs). How the regulatory modules functionally evolve in different species is not clear. A comparison of the CRMs for the Abdominal-B gene from different Drosophila species reveals relatively low levels of overall sequence conservation. However, embryonic enhancer CRMs from other Drosophila species direct transgenic reporter gene expression in the same spatial and temporal patterns during development as their D. melanogaster orthologs. Bioinformatic analysis reveals the presence of short conserved sequences within defined CRMs, representing gap and pair-rule transcription factor binding sites. One predicted binding site for the gap transcription factor KRUPPEL in the IAB5 CRM was found to be altered in Superabdominal (Sab) mutations. In Sab mutant flies, the third abdominal segment is transformed into a copy of the fifth abdominal segment. A model for KRUPPEL-mediated repression at this binding site is presented. These findings challenge our current understanding of the relationship between sequence evolution at the molecular level and functional activity of a CRM. While the overall sequence conservation at Drosophila CRMs is not distinctive from neighboring genomic regions, functionally critical transcription factor binding sites within embryonic enhancer CRMs are highly conserved. These results have implications for understanding mechanisms of gene expression during embryonic development, enhancer function, and the molecular evolution of eukaryotic regulatory modules.

Non-genic transcription at the Drosophila bithorax complex functional activity of the dark matter of the genome.

Drosophila melanogaster is a powerful model system for the study of gene regulation due to its short generation time, high fertility and the availability of various genetic tools to manipulate the genome. Investigation into the regulation of homeotic genes and their role in embryonic patterning during development was pioneered in Drosophila. Recently, the molecular mechanisms responsible for regulating gene expression in the bithorax complex have been the focus of active study. Many of these studies have pointed to the importance of cis-regulatory modules, genetic sequences that direct the temporal and spatial patterns of gene expression over large genomic distances. Additional components of the regulatory code have emerged beyond the primary DNA sequence. In particular, non-genic transcription is an important mechanism for controlling gene expression either through direct transcriptional mechanisms that mediate dynamic epigenetic control of the chromatin environment or through functional activity of the RNA products.