Spatial and single-cell profiling of the metabolome, transcriptome and epigenome of the aging mouse liver.

Tissues within an organism and even cell types within a tissue can age with different velocities. However, it is unclear whether cells of one type experience different aging trajectories within a tissue depending on their spatial location. Here, we used spatial transcriptomics in combination with single-cell ATAC-seq and RNA-seq, lipidomics and functional assays to address how cells in the male murine liver are affected by age-related changes in the microenvironment. Integration of the datasets revealed zonation-specific and age-related changes in metabolic states, the epigenome and transcriptome. The epigenome changed in a zonation-dependent manner and functionally, periportal hepatocytes were characterized by decreased mitochondrial fitness, whereas pericentral hepatocytes accumulated large lipid droplets. Together, we provide evidence that changing microenvironments within a tissue exert strong influences on their resident cells that can shape epigenetic, metabolic and phenotypic outputs.

Preliminary results from the EMoLung clinical study showing early lung cancer detection by the LC score.

BACKGROUND: Lung cancer (LC) causes more deaths worldwide than any other cancer type. Despite advances in therapeutic strategies, the fatality rate of LC cases remains high (95%) since the majority of patients are diagnosed at late stages when patient prognosis is poor. Analysis of the International Association for the Study of Lung Cancer (IASLC) database indicates that early diagnosis is significantly associated with favorable outcome. However, since symptoms of LC at early stages are unspecific and resemble those of benign pathologies, current diagnostic approaches are mostly initiated at advanced LC stages. METHODS: We developed a LC diagnosis test based on the analysis of distinct RNA isoforms expressed from the GATA6 and NKX2-1 gene loci, which are detected in exhaled breath condensates (EBCs). Levels of these transcript isoforms in EBCs were combined to calculate a diagnostic score (the LC score). In the present study, we aimed to confirm the applicability of the LC score for the diagnosis of early stage LC under clinical settings. Thus, we evaluated EBCs from patients with early stage, resectable non-small cell lung cancer (NSCLC), who were prospectively enrolled in the EMoLung study at three sites in Germany. RESULTS: LC score-based classification of EBCs confirmed its performance under clinical conditions, achieving a sensitivity of 95.7%, 91.3% and 84.6% for LC detection at stages I, II and III, respectively. CONCLUSIONS: The LC score is an accurate and non-invasive option for early LC diagnosis and a valuable complement to LC screening procedures based on computed tomography.

Group size planning for breedings of gene-modified mice and other organisms following Mendelian inheritance.

Colony management of gene-modified animals is time-consuming, costly and affected by random events related to Mendelian genetics, fertility and litter size. Careful planning is mandatory to ensure successful outcomes using the least number of animals, hence adhering to the 3R principles of animal welfare. Here we have developed an R package, accessible also through an interactive public website, that optimizes breeding design by providing information about the optimal number of breedings needed to obtain defined breeding outcomes, taking into account specific species, strain, or line properties and success probability. Our software also enables breeding planning for balanced male-to-female ratio or single-sex experiments. We show that, for single-sex designs, the necessary number of breedings is at least doubled compared to the use of all born animals. While the presented tool provides preset parameters for the laboratory mouse, it can be readily used for any other species.

Non-canonical integrin signaling activates EGFR and RAS-MAPK-ERK signaling in small cell lung cancer.

Background: Small cell lung cancer (SCLC) is an extremely aggressive cancer type with a patient median survival of 6-12 months. Epidermal growth factor (EGF) signaling plays an important role in triggering SCLC. In addition, growth factor-dependent signals and alpha-, beta-integrin (ITGA, ITGB) heterodimer receptors functionally cooperate and integrate their signaling pathways. However, the precise role of integrins in EGF receptor (EGFR) activation in SCLC remains elusive. Methods: We analyzed human precision-cut lung slices (hPCLS), retrospectively collected human lung tissue samples and cell lines by classical methods of molecular biology and biochemistry. In addition, we performed RNA-sequencing-based transcriptomic analysis in human lung cancer cells and human lung tissue samples, as well as high-resolution mass spectrometric analysis of the protein cargo from extracellular vesicles (EVs) that were isolated from human lung cancer cells. Results: Our results demonstrate that non-canonical ITGB2 signaling activates EGFR and RAS/MAPK/ERK signaling in SCLC. Further, we identified a novel SCLC gene expression signature consisting of 93 transcripts that were induced by ITGB2, which may be used for stratification of SCLC patients and prognosis prediction of LC patients. We also found a cell-cell communication mechanism based on EVs containing ITGB2, which were secreted by SCLC cells and induced in control human lung tissue RAS/MAPK/ERK signaling and SCLC markers. Conclusions: We uncovered a mechanism of ITGB2-mediated EGFR activation in SCLC that explains EGFR-inhibitor resistance independently of EGFR mutations, suggesting the development of therapies targeting ITGB2 for patients with this extremely aggressive lung cancer type.

Reliable genotyping of recombinant genomes using a robust hidden Markov model.

Meiotic recombination is an essential mechanism during sexual reproduction and includes the exchange of chromosome segments between homologous chromosomes. New allelic combinations are transmitted to the new generation, introducing novel genetic variation in the offspring genomes. With the improvement of high-throughput whole-genome sequencing technologies, large numbers of recombinant individuals can now be sequenced with low sequencing depth at low costs, necessitating computational methods for reconstructing their haplotypes. The main challenge is the uncertainty in haplotype calling that arises from the low information content of a single genomic position. Straightforward sliding window-based approaches are difficult to tune and fail to place recombination breakpoints precisely. Hidden Markov model (HMM)-based approaches, on the other hand, tend to over-segment the genome. Here, we present RTIGER, an HMM-based model that exploits in a mathematically precise way the fact that true chromosome segments typically have a certain minimum length. We further separate the task of identifying the correct haplotype sequence from the accurate placement of haplotype borders, thereby maximizing the accuracy of border positions. By comparing segmentations based on simulated data with known underlying haplotypes, we highlight the reasons for RTIGER outperforming traditional segmentation approaches. We then analyze the meiotic recombination pattern of segregants of 2 Arabidopsis (Arabidopsis thaliana) accessions and a previously described hyper-recombining mutant. RTIGER is available as an R package with an efficient Julia implementation of the core algorithm.

Investigating non-inferiority or equivalence in time-to-event data under non-proportional hazards.

The classical approach to analyze time-to-event data, e.g. in clinical trials, is to fit Kaplan-Meier curves yielding the treatment effect as the hazard ratio between treatment groups. Afterwards, a log-rank test is commonly performed to investigate whether there is a difference in survival or, depending on additional covariates, a Cox proportional hazard model is used. However, in numerous trials these approaches fail due to the presence of non-proportional hazards, resulting in difficulties of interpreting the hazard ratio and a loss of power. When considering equivalence or non-inferiority trials, the commonly performed log-rank based tests are similarly affected by a violation of this assumption. Here we propose a parametric framework to assess equivalence or non-inferiority for survival data. We derive pointwise confidence bands for both, the hazard ratio and the difference of the survival curves. Further we propose a test procedure addressing non-inferiority and equivalence by directly comparing the survival functions at certain time points or over an entire range of time. Once the model's suitability is proven the method provides a noticeable power benefit, irrespectively of the shape of the hazard ratio. On the other hand, model selection should be carried out carefully as misspecification may cause type I error inflation in some situations. We investigate the robustness and demonstrate the advantages and disadvantages of the proposed methods by means of a simulation study. Finally, we demonstrate the validity of the methods by a clinical trial example.

IGF1R expression by adult oligodendrocytes is not required in the steady-state but supports neuroinflammation.

In the central nervous system (CNS), insulin-like growth factor 1 (IGF-1) regulates myelination by oligodendrocyte (ODC) precursor cells and shows anti-apoptotic properties in neuronal cells in different in vitro and in vivo systems. Previous work also suggests that IGF-1 protects ODCs from cell death and enhances remyelination in models of toxin-induced and autoimmune demyelination. However, since evidence remains controversial, the therapeutic potential of IGF-1 in demyelinating CNS conditions is unclear. To finally shed light on the function of IGF1-signaling for ODCs, we deleted insulin-like growth factor 1 receptor (IGF1R) specifically in mature ODCs of the mouse. We found that ODC survival and myelin status were unaffected by the absence of IGF1R until 15 months of age, indicating that IGF-1 signaling does not play a major role in post-mitotic ODCs during homeostasis. Notably, the absence of IGF1R did neither affect ODC survival nor myelin status upon cuprizone intoxication or induction of experimental autoimmune encephalomyelitis (EAE), models for toxic and autoimmune demyelination, respectively. Surprisingly, however, the absence of IGF1R from ODCs protected against clinical neuroinflammation in the EAE model. Together, our data indicate that IGF-1 signaling is not required for the function and survival of mature ODCs in steady-state and disease.

Pseudotime analysis reveals novel regulatory factors for multigenic onset and monogenic transition of odorant receptor expression.

During their maturation from horizontal basal stem cells, olfactory sensory neurons (OSNs) are known to select exactly one out of hundreds of olfactory receptors (ORs) and express it on their surface, a process called monogenic selection. Monogenic expression is preceded by a multigenic phase during which several OR genes are expressed in a single OSN. Here, we perform pseudotime analysis of a single cell RNA-Seq dataset of murine olfactory epithelium to precisely align the multigenic and monogenic expression phases with the cell types occurring during OSN differentiation. In combination with motif analysis of OR gene cluster-associated enhancer regions, we identify known and novel transcription (co-)factors (Ebf1, Lhx2, Ldb1, Fos and Ssbp2) and chromatin remodelers (Kdm1a, Eed and Zmynd8) associated with OR expression. The inferred temporal order of their activity suggests novel mechanisms contributing to multigenic OR expression and monogenic selection.

RNA transcription and degradation of Alu retrotransposons depends on sequence features and evolutionary history.

Alu elements are one of the most successful groups of RNA retrotransposons and make up 11% of the human genome with over 1 million individual loci. They are linked to genetic defects, increases in sequence diversity, and influence transcriptional activity. Still, their RNA metabolism is poorly understood yet. It is even unclear whether Alu elements are mostly transcribed by RNA Polymerase II or III. We have conducted a transcription shutoff experiment by α-amanitin and metabolic RNA labeling by 4-thiouridine combined with RNA fragmentation (TT-seq) and RNA-seq to shed further light on the origin and life cycle of Alu transcripts. We find that Alu RNAs are more stable than previously thought and seem to originate in part from RNA Polymerase II activity, as previous reports suggest. Their expression however seems to be independent of the transcriptional activity of adjacent genes. Furthermore, we have developed a novel statistical test for detecting the expression of quantitative trait loci in Alu elements that relies on the de Bruijn graph representation of all Alu sequences. It controls for both statistical significance and biological relevance using a tuned k-mer representation, discovering influential sequence features missed by regular motif search. In addition, we discover several point mutations using a generalized linear model, and motifs of interest, which also match transcription factor-binding motifs.

Unmasking lower gastrointestinal bleeding under administration of norepinephrine.

BACKGROUND: Bleeding in the gastrointestinal tract is common and transarterial embolization enables the clinician to control gastrointestinal bleeding. Contrast extravasation is a prerequisite for successful embolization. Provocative angiography is helpful in the detection of elusive bleeding. AIM: We performed a retrospective analysis of angiographic treatment in patients with lower gastrointestinal hemorrhage and initially negative angiographies, as well as the role of norepinephrine (NE) in unmasking bleeding. METHODS: We analyzed 41 patients with lower gastrointestinal bleeding after angiography who had undergone treatment over a period of 10 years. All patients had a positive shock index and needed intensive care. RESULTS: In three of four patients, angiography disclosed the site of bleeding when NE was used during the procedure for hemodynamic stabilization. CONCLUSION: We suggest that angiography performed after the administration of NE in unstable patients with gastrointestinal bleeding and an initially negative angiography has the potential to unmask bleeding sites for successful embolization. However, this statement must be confirmed in prospective studies.

Developmental HCN channelopathy results in decreased neural progenitor proliferation and microcephaly in mice.

The development of the cerebral cortex relies on the controlled division of neural stem and progenitor cells. The requirement for precise spatiotemporal control of proliferation and cell fate places a high demand on the cell division machinery, and defective cell division can cause microcephaly and other brain malformations. Cell-extrinsic and -intrinsic factors govern the capacity of cortical progenitors to produce large numbers of neurons and glia within a short developmental time window. In particular, ion channels shape the intrinsic biophysical properties of precursor cells and neurons and control their membrane potential throughout the cell cycle. We found that hyperpolarization-activated cyclic nucleotide-gated cation (HCN) channel subunits are expressed in mouse, rat, and human neural progenitors. Loss of HCN channel function in rat neural stem cells impaired their proliferation by affecting the cell-cycle progression, causing G1 accumulation and dysregulation of genes associated with human microcephaly. Transgene-mediated, dominant-negative loss of HCN channel function in the embryonic mouse telencephalon resulted in pronounced microcephaly. Together, our findings suggest a role for HCN channel subunits as a part of a general mechanism influencing cortical development in mammals.

Endoscopic hemostasis makes the difference: Angiographic treatment in patients with lower gastrointestinal bleeding.

BACKGROUND: The large majority of gastrointestinal bleedings subside on their own or after endoscopic treatment. However, a small number of these may pose a challenge in terms of therapy because the patients develop hemodynamic instability, and endoscopy does not achieve adequate hemostasis. Interventional radiology supplemented with catheter angiography (CA) and transarterial embolization have gained importance in recent times. AIM: To evaluate clinical predictors for angiography in patients with lower gastrointestinal bleeding (LGIB). METHODS: We compared two groups of patients in a retrospective analysis. One group had been treated for more than 10 years with CA for LGIB (n = 41). The control group had undergone non-endoscopic or endoscopic treatment for two years and been registered in a bleeding registry (n = 92). The differences between the two groups were analyzed using decision trees with the goal of defining clear rules for optimal treatment. RESULTS: Patients in the CA group had a higher shock index, a higher Glasgow-Blatchford bleeding score (GBS), lower serum hemoglobin levels, and more rarely achieved hemostasis in primary endoscopy. These patients needed more transfusions, had longer hospital stays, and had to undergo subsequent surgery more frequently (P < 0.001). CONCLUSION: Endoscopic hemostasis proved to be the crucial difference between the two patient groups. Primary endoscopic hemostasis, along with GBS and the number of transfusions, would permit a stratification of risks. After prospective confirmation of the present findings, the use of decision trees would permit the identification of patients at risk for subsequent diagnosis and treatment based on interventional radiology.

Computational Tumor Infiltration Phenotypes Enable the Spatial and Genomic Analysis of Immune Infiltration in Colorectal Cancer.

Cancer immunotherapy has led to significant therapeutic progress in the treatment of metastatic and formerly untreatable tumors. However, drug response rates are variable and often only a subgroup of patients will show durable response to a treatment. Biomarkers that help to select those patients that will benefit the most from immunotherapy are thus of crucial importance. Here, we aim to identify such biomarkers by investigating the tumor microenvironment, i.e., the interplay between different cell types like immune cells, stromal cells and malignant cells within the tumor and developed a computational method that determines spatial tumor infiltration phenotypes. Our method is based on spatial point pattern analysis of immunohistochemically stained colorectal cancer tumor tissue and accounts for the intra-tumor heterogeneity of immune infiltration. We show that, compared to base-line models, tumor infiltration phenotypes provide significant additional support for the prediction of established biomarkers in a colorectal cancer patient cohort (n = 80). Integration of tumor infiltration phenotypes with genetic and genomic data from the same patients furthermore revealed significant associations between spatial infiltration patterns and common mutations in colorectal cancer and gene expression signatures. Based on these associations, we computed novel gene signatures that allow one to predict spatial tumor infiltration patterns from gene expression data only and validated this approach in a separate dataset from the Cancer Genome Atlas.

FADD and Caspase-8 Regulate Gut Homeostasis and Inflammation by Controlling MLKL- and GSDMD-Mediated Death of Intestinal Epithelial Cells.

Pathways controlling intestinal epithelial cell (IEC) death regulate gut immune homeostasis and contribute to the pathogenesis of inflammatory bowel diseases. Here, we show that caspase-8 and its adapter FADD act in IECs to regulate intestinal inflammation downstream of Z-DNA binding protein 1 (ZBP1)- and tumor necrosis factor receptor-1 (TNFR1)-mediated receptor interacting protein kinase 1 (RIPK1) and RIPK3 signaling. Mice with IEC-specific FADD or caspase-8 deficiency developed colitis dependent on mixed lineage kinase-like (MLKL)-mediated epithelial cell necroptosis. However, MLKL deficiency fully prevented ileitis caused by epithelial caspase-8 ablation, but only partially ameliorated ileitis in mice lacking FADD in IECs. Our genetic studies revealed that caspase-8 and gasdermin-D (GSDMD) were both required for the development of MLKL-independent ileitis in mice with epithelial FADD deficiency. Therefore, FADD prevents intestinal inflammation downstream of ZBP1 and TNFR1 by inhibiting both MLKL-induced necroptosis and caspase-8-GSDMD-dependent pyroptosis-like death of epithelial cells.

Benefits of a factorial design focusing on inclusion of female and male animals in one experiment.

Disease occurrence, clinical manifestations, and outcomes differ between men and women. Yet, women and men are most of the time treated similarly, which is often based on experimental data over-representing one sex. Accounting for persisting sex bias in biomedical research is the misconception that the analysis of sex-specific effects would double sample size and costs. We designed an analysis to test the potential benefits of a factorial study design in the context of a study including male and female animals. We chose a 2 × 2 factorial design approach to study the effect of treatment, sex, and an interaction term of treatment and sex in a hypothetical situation. We calculated the sample sizes required to detect an effect of a given magnitude with sufficient power and under different experimental setups. We demonstrated that the inclusion of both sexes in experimental setups, without testing for sex effects, requires no or few additional animals in our scenarios. These experimental designs still allow for the exploration of sex effects at low cost. In a confirmatory instead of an exploratory design, we observed an increase in total sample sizes by 33%, at most. Since the complexities associated with this mathematical model require statistical expertise, we generated and provide a sample size calculator for planning factorial design experiments. For the inclusion of sex, a factorial design is advisable, and a sex-specific analysis can be performed without excessive additional effort. Our easy-to-use calculation tool provides help in designing studies with both sexes and addresses the current sex bias in preclinical studies. KEY MESSAGES: • Both sexes should be included into animal studies. • Exploratory study of sex effects can be conducted with no or small increase in animal number. • Confirmatory analysis of sex effects requires maximum 33% more animals per study. • Our calculation tool supports the design of studies with both sexes.

Clustering of samples with a tree-shaped dependence structure, with an application to microscopic time lapse imaging.

MOTIVATION: Recent imaging technologies allow for high-throughput tracking of cells as they migrate, divide, express fluorescent markers and change their morphology. The interpretation of these data requires unbiased, efficient statistical methods that model the dynamics of cell phenotypes. RESULTS: We introduce treeHFM, a probabilistic model which generalizes the theory of hidden Markov models to tree structured data. While accounting for the entire genealogy of a cell, treeHFM categorizes cells according to their primary phenotypic features. It models all relevant events in a cell's life, including cell division, and thereby enables the analysis of event order and cell fate heterogeneity. Simulations show higher accuracy in predicting correct state labels when modeling the more complex, tree-shaped dependency of samples over standard HMM modeling. Applying treeHFM to time lapse images of hematopoietic progenitor cell differentiation, we demonstrate that progenitor cells undergo a well-ordered sequence of differentiation events. AVAILABILITY AND IMPLEMENTATION: The treeHFM is implemented in C++. We provide wrapper functions for the programming languages R (CRAN package, https://CRAN.R-project.org/package=treeHFM) and Matlab (available at Mathworks Central, http://se.mathworks.com/matlabcentral/fileexchange/57575-treehfml). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MowJoe: a method for automated-high throughput dissected leaf phenotyping.

BACKGROUND: Accurate and automated phenotyping of leaf images is necessary for high throughput studies of leaf form like genome-wide association analysis and other forms of quantitative trait locus mapping. Dissected leaves (also referred to as compound) that are subdivided into individual units are an attractive system to study diversification of form. However, there are only few software tools for their automated analysis. Thus, high-throughput image processing algorithms are needed that can partition these leaves in their phenotypically relevant units and calculate morphological features based on these units. RESULTS: We have developed MowJoe, an image processing algorithm that dissects a dissected leaf into leaflets, petiolule, rachis and petioles. It employs image skeletonization to convert leaves into graphs, and thereafter applies algorithms operating on graph structures. This partitioning of a leaf allows the derivation of morphological features such as leaf size, or eccentricity of leaflets. Furthermore, MowJoe automatically places landmarks onto the terminal leaflet that can be used for further leaf shape analysis. It generates specific output files that can directly be imported into downstream shape analysis tools. We applied the algorithm to two accessions of Cardamine hirsuta and show that our features are able to robustly discriminate between these accessions. CONCLUSION: MowJoe is a tool for the semi-automated, quantitative high throughput shape analysis of dissected leaf images. It provides the statistical power for the detection of the genetic basis of quantitative morphological variations.

G-EYE advanced colonoscopy for improved polyp detection rates - a randomized tandem pilot study with different endoscopists.

BACKGROUND AND AIMS: The most commonly missed polyps in colonoscopy are those located behind haustral folds. The G-EYE system is a standard colonoscope consisting of re-processable balloon at its distal tip. The G-EYE balloon improves the detection of polyps by straightening the haustral folds. In our back-to-back tandem study, we aimed to determine whether and to what extent the G-EYE system could reduce adenoma miss rates in screening colonoscopy. METHODS: Patients referred to colonoscopy were randomized into 2 groups. Group A underwent a standard colonoscopy (SC) followed by balloon colonoscopy (BC), and Group B underwent BC followed by SC. In this randomized tandem study, the investigator's level of training and the endoscopists themselves were changed after each withdrawal. Each endoscopist was blinded to the results of the first withdrawal. RESULTS: Fifty-eight patients were enrolled and randomized into 2 groups with similar baseline characteristics. Nine patients were excluded from the study. Twenty-five patients underwent SC followed by BC while 24 underwent BC followed by SC. The adenoma miss rate for SC was 41 %, with an additional detection rate of 69 % for BC (ratio 1.69). The overall miss rate for polyps was 60 % for SC, with an additional detection rate of 150 % for BC (ratio 2.5). Experienced investigators who used BC were able to identify an additional 7 polyps while inexperienced investigators. CONCLUSIONS: Although our results could not clearly confirm that BC improves adenoma detection, the investigator's experience appears to be a major determinant of the adenoma detection rate.

The Manchester Triage System (MTS): a score for emergency management of patients with acute gastrointestinal bleeding.

BACKGROUND: Suspected gastrointestinal (GI) bleeding is a common initial diagnosis in emergency departments. Despite existing endoscopic scores to estimate the risk of GI bleeding, the primary clinical assessment of urgency can remain challenging. The 5-step Manchester Triage System (MTS) is a validated score that is often applied for the initial assessment of patients presenting in emergency departments. METHODS: All computer-based records of patients who were admitted between January 2014 and December 2014 to our emergency department in a tertiary referral hospital were analyzed retrospectively. The aim of our retrospective analysis was to determine if patient triage using the MTS is associated with rates of endoscopy and with presence of active GI bleeding. RESULTS: In summary, 5689 patients with a GI condition were treated at our emergency department. Two hundred eighty-four patients (4.9 %) presented with suspected GI bleeding, and 165 patients (58 %) received endoscopic diagnostic. Endoscopic intervention for hemostasis was needed in 34 patients (21 %). In patients who underwent emergency endoscopy, triage into MTS categories with higher urgency was associated with higher rates of endoscopic confirmation of suspected GI bleeding (79 % of patients with MTS priority levels 1 or 2, 53 % in level 3 patients, and 40 % in levels 4 or 5 patients; p = 0.024). CONCLUSIONS: The MTS is an established tool for triage in emergency departments and could have a potential to guide early clinical decision-making with regards to urgency of endoscopic evaluation in patients with suspected GI bleeding.

GenoGAM: genome-wide generalized additive models for ChIP-Seq analysis.

MOTIVATION: Chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) is a widely used approach to study protein-DNA interactions. Often, the quantities of interest are the differential occupancies relative to controls, between genetic backgrounds, treatments, or combinations thereof. Current methods for differential occupancy of ChIP-Seq data rely however on binning or sliding window techniques, for which the choice of the window and bin sizes are subjective. RESULTS: Here, we present GenoGAM (Genome-wide Generalized Additive Model), which brings the well-established and flexible generalized additive models framework to genomic applications using a data parallelism strategy. We model ChIP-Seq read count frequencies as products of smooth functions along chromosomes. Smoothing parameters are objectively estimated from the data by cross-validation, eliminating ad hoc binning and windowing needed by current approaches. GenoGAM provides base-level and region-level significance testing for full factorial designs. Application to a ChIP-Seq dataset in yeast showed increased sensitivity over existing differential occupancy methods while controlling for type I error rate. By analyzing a set of DNA methylation data and illustrating an extension to a peak caller, we further demonstrate the potential of GenoGAM as a generic statistical modeling tool for genome-wide assays. AVAILABILITY AND IMPLEMENTATION: Software is available from Bioconductor: https://www.bioconductor.org/packages/release/bioc/html/GenoGAM.html . CONTACT: gagneur@in.tum.de. SUPPLEMENTARY INFORMATION: Supplementary information is available at Bioinformatics online.