Discrete-state models identify pathway specific B cell states across diseases and infections at single-cell resolution.

Oxygen (O2) regulated pathways modulate B cell activation, migration and proliferation during infection, vaccination, and other diseases. Modeling these pathways in health and disease is critical to understand B cell states and ways to mediate them. To characterize B cells by their activation of O2 regulated pathways we develop pathway specific discrete state models using previously published single-cell RNA-sequencing (scRNA-seq) datasets from isolated B cells. Specifically, Single Cell Boolean Omics Network Invariant-Time Analysis (scBONITA) was used to infer logic gates for known pathway topologies. The simplest inferred set of logic gates that maximized the number of "OR" interactions between genes was used to simulate B cell networks involved in oxygen sensing until they reached steady network states (attractors). By focusing on the attractors that best represented sequenced cells, we identified genes critical in determining pathway specific cellular states that corresponded to diseased and healthy B cell phenotypes. Specifically, we investigate the transendothelial migration, regulation of actin cytoskeleton, HIF1A, and Citrate Cycle pathways. Our analysis revealed attractors that resembled the state of B cell exhaustion in HIV+ patients as well as attractors that promoted anerobic metabolism, angiogenesis, and tumorigenesis in breast cancer patients, which were eliminated after neoadjuvant chemotherapy (NACT). Finally, we investigated the attractors to which the Azimuth-annotated B cells mapped and found that attractors resembling B cells from HIV+ patients encompassed a significantly larger number of atypical memory B cells than HIV- attractors. Meanwhile, attractors resembling B cells from breast cancer patients post NACT encompassed a reduced number of atypical memory B cells compared to pre-NACT attractors.

Pillars of biology: Boolean modeling of gene-regulatory networks.

Boolean modeling is a mathematical modeling framework used for defining and studying gene-regulatory networks (GRNs). It serves as a means to develop mechanistic models, offering insights into the trajectories and dynamic properties of GRNs. In this review, I delve into seminal papers published in the Journal of Theoretical Biology that have spearheaded this field. Additionally, I explore the application of these modeling methods in the current era of data-intensive science.

The C. elegans Myc-family of transcription factors coordinate a dynamic adaptive response to dietary restriction.

Dietary restriction (DR), the process of decreasing overall food consumption over an extended period of time, has been shown to increase longevity across evolutionarily diverse species and delay the onset of age-associated diseases in humans. In Caenorhabditis elegans, the Myc-family transcription factors (TFs) MXL-2 (Mlx) and MML-1 (MondoA/ChREBP), which function as obligate heterodimers, and PHA-4 (orthologous to forkhead box transcription factor A) are both necessary for the full physiological benefits of DR. However, the adaptive transcriptional response to DR and the role of MML-1::MXL-2 and PHA-4 remains elusive. We identified the transcriptional signature of C. elegans DR, using the eat-2 genetic model, and demonstrate broad changes in metabolic gene expression in eat-2 DR animals, which requires both mxl-2 and pha-4. While the requirement for these factors in DR gene expression overlaps, we found many of the DR genes exhibit an opposing change in relative gene expression in eat-2;mxl-2 animals compared to wild-type, which was not observed in eat-2 animals with pha-4 loss. We further show functional deficiencies of the mxl-2 loss in DR outside of lifespan, as eat-2;mxl-2 animals exhibit substantially smaller brood sizes and lay a proportion of dead eggs, indicating that MML-1::MXL-2 has a role in maintaining the balance between resource allocation to the soma and to reproduction under conditions of chronic food scarcity. While eat-2 animals do not show a significantly different metabolic rate compared to wild-type, we also find that loss of mxl-2 in DR does not affect the rate of oxygen consumption in young animals. The gene expression signature of eat-2 mutant animals is consistent with optimization of energy utilization and resource allocation, rather than induction of canonical gene expression changes associated with acute metabolic stress -such as induction of autophagy after TORC1 inhibition. Consistently, eat-2 animals are not substantially resistant to stress, providing further support to the idea that chronic DR may benefit healthspan and lifespan through efficient use of limited resources rather than broad upregulation of stress responses, and also indicates that MML-1::MXL-2 and PHA-4 may have different roles in promotion of benefits in response to different pro-longevity stimuli.

Homeodomain-interacting protein kinase maintains neuronal homeostasis during normal Caenorhabditis elegans aging and systemically regulates longevity from serotonergic and GABAergic neurons.

Aging and the age-associated decline of the proteome is determined in part through neuronal control of evolutionarily conserved transcriptional effectors, which safeguard homeostasis under fluctuating metabolic and stress conditions by regulating an expansive proteostatic network. We have discovered the Caenorhabditis elegans homeodomain-interacting protein kinase (HPK-1) acts as a key transcriptional effector to preserve neuronal integrity, function, and proteostasis during aging. Loss of hpk-1 results in drastic dysregulation in expression of neuronal genes, including genes associated with neuronal aging. During normal aging hpk-1 expression increases throughout the nervous system more broadly than any other kinase. Within the aging nervous system, hpk-1 induction overlaps with key longevity transcription factors, which suggests that hpk-1 expression mitigates natural age-associated physiological decline. Consistently, pan-neuronal overexpression of hpk-1 extends longevity, preserves proteostasis both within and outside of the nervous system, and improves stress resistance. Neuronal HPK-1 improves proteostasis through kinase activity. HPK-1 functions cell non-autonomously within serotonergic and γ-aminobutyric acid (GABA)ergic neurons to improve proteostasis in distal tissues by specifically regulating distinct components of the proteostatic network. Increased serotonergic HPK-1 enhances the heat shock response and survival to acute stress. In contrast, GABAergic HPK-1 induces basal autophagy and extends longevity, which requires mxl-2 (MLX), hlh-30 (TFEB), and daf-16 (FOXO). Our work establishes hpk-1 as a key neuronal transcriptional regulator critical for preservation of neuronal function during aging. Further, these data provide novel insight as to how the nervous system partitions acute and chronic adaptive response pathways to delay aging by maintaining organismal homeostasis.

Proteins are essential for nearly every cellular process to sustain a healthy organism. A complex network of pathways and signalling molecules regulates the proteins so that they work correctly in a process known as proteostasis. As the body ages, this network can become damaged, which leads to the production of faulty proteins. Many proteins end up being misfolded ­ in other words, they are misshapen on the molecular level, which can be toxic for the cell. A build-up of such misfolded proteins is implicated in several neurological conditions, including Alzheimer's, Parkinson's and Huntington's disease. Cells have various ways to detect and respond to internal stressors, such as tissue or organ damage. For example, specific proteins in the nervous system can raise a 'central' alert when damage is detected, which then primes and coordinates the body's systems to respond in the peripheral cells and tissues. But exactly how this happens is still unclear. To find out more about the central coordination of stress responses, Lazaro-Pena et al. studied one such sensor protein, called HPK-1, in the roundworm C. elegans. They first overexpressed the protein in various tissues. This revealed that only when HPK-1 was overactive in nerve tissue, it protected proteins and prolonged the lifespan of the worms. An increased amount of HPK-1 improved the health span of the worms and older worms also moved better. However, genetically manipulated worms lacking HPK-1 in their nerve cells showed a faster decline in nervous system health as they aged, which could be reversed once HPK-1 was activated again. Lazaro-Pena et al. then measured the amount of HPK-1 in worms at different stages of their life. This showed that as the worms aged, the amount of HPK-1 increased in the nerve cells. The nerve cells in which HPK-1 levels increased overlapped with an increased expression of proteins associated with longevity. Moreover, when HPK-1 was overexpressed, it stimulated the release of other cell signals, which then triggered protective responses to prevent the misfolding and aggregation of proteins and to help degrade damaged proteins. This study shows for the first time that HPK-1 appears to play a protective role during normal ageing and that it may act as a key switch to stimulate other protective mechanisms. These findings may give rise to new insights into how the nervous system can coordinate many different stress responses, and ultimately delay ageing throughout the whole body.

Generation and Functional Analysis of Defective Viral Genomes during SARS-CoV-2 Infection.

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from transcriptome sequencing (RNA-seq) data sets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hot spots were identified for DVG recombination, and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single-cell RNA-seq analysis indicated the interferon (IFN) stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the next-generation sequencing (NGS) data set from a published cohort study and observed a significantly higher amount and frequency of DVG in symptomatic patients than those in asymptomatic patients. Finally, we observed exceptionally diverse DVG populations in one immunosuppressive patient up to 140 days after the first positive test of COVID-19, suggesting for the first time an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and into how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. IMPORTANCE Defective viral genomes (DVGs) are generated ubiquitously in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide the potential for them to be used in novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex, and this recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hot spots for nonhomologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide evidence to harness the immunostimulatory potential of DVGs in the development of a vaccine and antivirals for SARS-CoV-2.

Executable Network Models of Integrated Multiomics Data.

Multiomics profiling provides a holistic picture of a condition being examined and captures the complexity of signaling events, beginning from the original cause (environmental or genetic), to downstream functional changes at multiple molecular layers. Pathway enrichment analysis has been used with multiomics data sets to characterize signaling mechanisms. However, technical and biological variability between these layered data limit an integrative computational analyses. We present a Boolean network-based method, multiomics Boolean Omics Network Invariant-Time Analysis (mBONITA), to integrate omics data sets that quantify multiple molecular layers. mBONITA utilizes prior knowledge networks to perform topology-based pathway analysis. In addition, mBONITA identifies genes that are consistently modulated across molecular measurements by combining observed fold-changes and variance, with a measure of node (i.e., gene or protein) influence over signaling, and a measure of the strength of evidence for that gene across data sets. We used mBONITA to integrate multiomics data sets from RAMOS B cells treated with the immunosuppressant drug cyclosporine A under varying O2 tensions to identify pathways involved in hypoxia-mediated chemotaxis. We compare mBONITA's performance with 6 other pathway analysis methods designed for multiomics data and show that mBONITA identifies a set of pathways with evidence of modulation across all omics layers. mBONITA is freely available at https://github.com/Thakar-Lab/mBONITA.

Preparation of noninfectious scRNAseq samples from SARS-CoV-2-infected epithelial cells.

Coronavirus disease (COVID-19) is an infectious disease caused by the SARS coronavirus 2 (SARS-CoV-2) virus. Direct assessment, detection, and quantitative analysis using high throughput methods like single-cell RNA sequencing (scRNAseq) is imperative to understanding the host response to SARS-CoV-2. One barrier to studying SARS-CoV-2 in the laboratory setting is the requirement to process virus-infected cell cultures, and potentially infectious materials derived therefrom, under Biosafety Level 3 (BSL-3) containment. However, there are only 190 BSL3 laboratory facilities registered with the U.S. Federal Select Agent Program, as of 2020, and only a subset of these are outfitted with the equipment needed to perform high-throughput molecular assays. Here, we describe a method for preparing non-hazardous RNA samples from SARS-CoV-2 infected cells, that enables scRNAseq analyses to be conducted safely in a BSL2 facility-thereby making molecular assays of SARS-CoV-2 cells accessible to a much larger community of researchers. Briefly, we infected African green monkey kidney epithelial cells (Vero-E6) with SARS-CoV-2 for 96 hours, trypsin-dissociated the cells, and inactivated them with methanol-acetone in a single-cell suspension. Fixed cells were tested for the presence of infectious SARS-CoV-2 virions using the Tissue Culture Infectious Dose Assay (TCID50), and also tested for viability using flow cytometry. We then tested the dissociation and methanol-acetone inactivation method on primary human lung epithelial cells that had been differentiated on an air-liquid interface. Finally, we performed scRNAseq quality control analysis on the resulting cell populations to evaluate the effects of our virus inactivation and sample preparation protocol on the quality of the cDNA produced. We found that methanol-acetone inactivated SARS-CoV-2, fixed the lung epithelial cells, and could be used to obtain noninfectious, high-quality cDNA libraries. This methodology makes investigating SARS-CoV-2, and related high-containment RNA viruses at a single-cell level more accessible to an expanded community of researchers.

Homeodomain-interacting protein kinase maintains neuronal homeostasis during normal Caenorhabditis elegans aging and systemically regulates longevity from serotonergic and GABAergic neurons.

Aging and the age-associated decline of the proteome is determined in part through neuronal control of evolutionarily conserved transcriptional effectors, which safeguard homeostasis under fluctuating metabolic and stress conditions by regulating an expansive proteostatic network. We have discovered the Caenorhabditis elegans h omeodomain-interacting p rotein k inase (HPK-1) acts as a key transcriptional effector to preserve neuronal integrity, function, and proteostasis during aging. Loss of hpk-1 results in drastic dysregulation in expression of neuronal genes, including genes associated with neuronal aging. During normal aging hpk-1 expression increases throughout the nervous system more broadly than any other kinase. Within the aging nervous system, hpk-1 induction overlaps with key longevity transcription factors, which suggests hpk-1 expression mitigates natural age-associated physiological decline. Consistently, pan-neuronal overexpression of hpk-1 extends longevity, preserves proteostasis both within and outside of the nervous system, and improves stress resistance. Neuronal HPK-1 improves proteostasis through kinase activity. HPK-1 functions cell non-autonomously within serotonergic and GABAergic neurons to improve proteostasis in distal tissues by specifically regulating distinct components of the proteostatic network. Increased serotonergic HPK-1 enhances the heat shock response and survival to acute stress. In contrast, GABAergic HPK-1 induces basal autophagy and extends longevity, which requires mxl-2 (MLX), hlh-30 (TFEB), and daf-16 (FOXO). Our work establishes hpk-1 as a key neuronal transcriptional regulator critical for preservation of neuronal function during aging. Further, these data provide novel insight as to how the nervous system partitions acute and chronic adaptive response pathways to delay aging by maintaining organismal homeostasis.

Executable models of immune signaling pathways in HIV-associated atherosclerosis.

Atherosclerosis (AS)-associated cardiovascular disease is an important cause of mortality in an aging population of people living with HIV (PLWH). This elevated risk has been attributed to viral infection, anti-retroviral therapy, chronic inflammation, and lifestyle factors. However, the rates at which PLWH develop AS vary even after controlling for length of infection, treatment duration, and for lifestyle factors. To investigate the molecular signaling underlying this variation, we sequenced 9368 peripheral blood mononuclear cells (PBMCs) from eight PLWH, four of whom have atherosclerosis (AS+). Additionally, a publicly available dataset of PBMCs from persons before and after HIV infection was used to investigate the effect of acute HIV infection. To characterize dysregulation of pathways rather than just measuring enrichment, we developed the single-cell Boolean Omics Network Invariant Time Analysis (scBONITA) algorithm. scBONITA infers executable dynamic pathway models and performs a perturbation analysis to identify high impact genes. These dynamic models are used for pathway analysis and to map sequenced cells to characteristic signaling states (attractor analysis). scBONITA revealed that lipid signaling regulates cell migration into the vascular endothelium in AS+ PLWH. Pathways implicated included AGE-RAGE and PI3K-AKT signaling in CD8+ T cells, and glucagon and cAMP signaling pathways in monocytes. Attractor analysis with scBONITA facilitated the pathway-based characterization of cellular states in CD8+ T cells and monocytes. In this manner, we identify critical cell-type specific molecular mechanisms underlying HIV-associated atherosclerosis using a novel computational method.

Baseline host determinants of robust human HIV-1 vaccine-induced immune responses: A meta-analysis of 26 vaccine regimens.

BACKGROUND: The identification of baseline host determinants that associate with robust HIV-1 vaccine-induced immune responses could aid HIV-1 vaccine development. We aimed to assess both the collective and relative performance of baseline characteristics in classifying individual participants in nine different Phase 1-2 HIV-1 vaccine clinical trials (26 vaccine regimens, conducted in Africa and in the Americas) as High HIV-1 vaccine responders. METHODS: This was a meta-analysis of individual participant data, with studies chosen based on participant-level (vs. study-level summary) data availability within the HIV-1 Vaccine Trials Network. We assessed the performance of 25 baseline characteristics (demographics, safety haematological measurements, vital signs, assay background measurements) and estimated the relative importance of each characteristic in classifying 831 participants as High (defined as within the top 25th percentile among positive responders or above the assay upper limit of quantification) versus Non-High responders. Immune response outcomes included HIV-1-specific serum IgG binding antibodies and Env-specific CD4+ T-cell responses assessed two weeks post-last dose, all measured at central HVTN laboratories. Three variable importance approaches based on SuperLearner ensemble machine learning were considered. FINDINGS: Overall, 30.1%, 50.5%, 36.2%, and 13.9% of participants were categorized as High responders for gp120 IgG, gp140 IgG, gp41 IgG, and Env-specific CD4+ T-cell vaccine-induced responses, respectively. When including all baseline characteristics, moderate performance was achieved for the classification of High responder status for the binding antibody responses, with cross-validated areas under the ROC curve (CV-AUC) of 0.72 (95% CI: 0.68, 0.76) for gp120 IgG, 0.73 (0.69, 0.76) for gp140 IgG, and 0.67 (95% CI: 0.63, 0.72) for gp41 IgG. In contrast, the collection of all baseline characteristics yielded little improvement over chance for predicting High Env-specific CD4+ T-cell responses [CV-AUC: 0.53 (0.48, 0.58)]. While estimated variable importance patterns differed across the three approaches, female sex assigned at birth, lower height, and higher total white blood cell count emerged as significant predictors of High responder status across multiple immune response outcomes using Approach 1. Of these three baseline variables, total white blood cell count ranked highly across all three approaches for predicting vaccine-induced gp41 and gp140 High responder status. INTERPRETATION: The identified features should be studied further in pursuit of intervention strategies to improve vaccine responses and may be adjusted for in analyses of immune response data to enhance statistical power. FUNDING: National Institute of Allergy and Infectious Diseases (UM1AI068635 to YH, UM1AI068614 to GDT, UM1AI068618 to MJM, and UM1 AI069511 to MCK), the Duke CFAR P30 AI064518 to GDT, and National Institute of Dental and Craniofacial Research (R01DE027245 to JJK). This work was also supported by the Bill and Melinda Gates Foundation. The content is solely the responsibility of the authors and does not necessarily represent the official views of any of the funding sources.

Generation and functional analysis of defective viral genomes during SARS-CoV-2 infection.

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from RNA-seq datasets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hotspots were identified for DVG recombination and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single cell RNA-seq analysis indicated the IFN stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the NGS dataset from a published cohort study and observed significantly higher DVG amount and frequency in symptomatic patients than that in asymptomatic patients. Finally, we observed unusually high DVG frequency in one immunosuppressive patient up to 140 days after admitted to hospital due to COVID-19, first-time suggesting an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. Importance: Defective viral genomes (DVGs) are ubiquitously generated in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide them the potential for novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex and the recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hotspots for non-homologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide the evidence to harness DVGs’ immunostimulatory potential in the development of vaccine and antivirals for SARS-CoV-2.

Cyclosporine A Modulates LSP1 Protein Levels in Human B Cells to Attenuate B Cell Migration at Low O2 Levels.

Coordinated migration of B cells within and between secondary lymphoid tissues is required for robust antibody responses to infection or vaccination. Secondary lymphoid tissues normally expose B cells to a low O2 (hypoxic) environment. Recently, we have shown that human B cell migration is modulated by an O2-dependent molecular switch, centrally controlled by the hypoxia-induced (transcription) factor-1α (HIF1A), which can be disrupted by the immunosuppressive calcineurin inhibitor, cyclosporine A (CyA). However, the mechanisms by which low O2 environments attenuate B cell migration remain poorly defined. Proteomics analysis has linked CXCR4 chemokine receptor signaling to cytoskeletal rearrangement. We now hypothesize that the pathways linking the O2 sensing molecular switch to chemokine receptor signaling and cytoskeletal rearrangement would likely contain phosphorylation events, which are typically missed in traditional transcriptomic and/or proteomic analyses. Hence, we have performed a comprehensive phosphoproteomics analysis of human B cells treated with CyA after engagement of the chemokine receptor CXCR4 with CXCL12. Statistical analysis of the separate and synergistic effects of CyA and CXCL12 revealed 116 proteins whose abundance is driven by a synergistic interaction between CyA and CXCL12. Further, we used our previously described algorithm BONITA to reveal a critical role for Lymphocyte Specific Protein 1 (LSP1) in cytoskeletal rearrangement. LSP1 is known to modulate neutrophil migration. Validating these modeling results, we show experimentally that LSP1 levels in B cells increase with low O2 exposure, and CyA treatment results in decreased LSP1 protein levels. This correlates with the increased chemotactic activity observed after CyA treatment. Lastly, we directly link LSP1 levels to chemotactic capacity, as shRNA knock-down of LSP1 results in significantly increased B cell chemotaxis at low O2 levels. These results directly link CyA to LSP1-dependent cytoskeletal regulation, demonstrating a previously unrecognized mechanism by which CyA modulates human B cell migration. Data are available via ProteomeXchange with identifier PXD036167.

Biomarkers of Development of Immunity and Allergic Diseases in Farming and Non-farming Lifestyle Infants: Design, Methods and 1 Year Outcomes in the "Zooming in to Old Order Mennonites" Birth Cohort Study.

Traditional farming lifestyle has been shown to be protective against asthma and allergic diseases. The individual factors that appear to be associated with this "farm-life effect" include consumption of unpasteurized farm milk and exposure to farm animals and stables. However, the biomarkers of the protective immunity and those associated with early development of allergic diseases in infancy remain unclear. The "Zooming in to Old Order Mennonites (ZOOM)" study was designed to assess the differences in the lifestyle and the development of the microbiome, systemic and mucosal immunity between infants born to traditional farming lifestyle at low risk for allergic diseases and those born to urban/suburban atopic families with a high risk for allergic diseases in order to identify biomarkers of development of allergic diseases in infancy. 190 mothers and their infants born to Old Order Mennonite population protected from or in Rochester families at high risk for allergic diseases were recruited before birth from the Finger Lakes Region of New York State. Questionnaires and samples are collected from mothers during pregnancy and after delivery and from infants at birth and at 1-2 weeks, 6 weeks, 6, 12, 18, and 24 months, with 3-, 4-, and 5-year follow-up ongoing. Samples collected include maternal blood, stool, saliva, nasal and skin swabs and urine during pregnancy; breast milk postnatally; infant blood, stool, saliva, nasal and skin swabs. Signs and symptoms of allergic diseases are assessed at every visit and serum specific IgE is measured at 1 and 2 years of age. Allergic diseases are diagnosed by clinical history, exam, and sensitization by skin prick test and/or serum specific IgE. By the end of the first year of life, the prevalence of food allergy and atopic dermatitis were higher in ROC infants compared to the rates observed in OOM infants as was the number of infants sensitized to foods. These studies of immune system development in a population protected from and in those at risk for allergic diseases will provide critical new knowledge about the development of the mucosal and systemic immunity and lay the groundwork for future studies of prevention of allergic diseases.

The Replica Set Method is a Robust, Accurate, and High-Throughput Approach for Assessing and Comparing Lifespan in C. elegans Experiments.

The advent of feeding based RNAi in Caenorhabditis elegans led to an era of gene discovery in aging research. Hundreds of gerogenes were discovered, and many are evolutionarily conserved, raising the exciting possibility that the underlying genetic basis for healthy aging in higher vertebrates could be quickly deciphered. Yet, the majority of putative gerogenes have still only been cursorily characterized, highlighting the need for high-throughput, quantitative assessments of changes in aging. A widely used surrogate measure of aging is lifespan. The traditional way to measure mortality in C. elegans tracks the deaths of individual animals over time within a relatively small population. This traditional method provides straightforward, direct measurements of median and maximum lifespan for the sampled population. However, this method is time consuming, often underpowered, and involves repeated handling of a set of animals over time, which in turn can introduce contamination or possibly damage increasingly fragile, aged animals. We have previously developed an alternative "Replica Set" methodology, which minimizes handling and increases throughput by at least an order of magnitude. The Replica Set method allows changes in lifespan to be measured for over one hundred feeding-based RNAi clones by one investigator in a single experiment- facilitating the generation of large quantitative phenotypic datasets, a prerequisite for development of biological models at a systems level. Here, we demonstrate through analysis of lifespan experiments simulated in silico that the Replica Set method is at least as precise and accurate as the traditional method in evaluating and estimating lifespan, and requires many fewer total animal observations across the course of an experiment. Furthermore, we show that the traditional approach to lifespan experiments is more vulnerable than the Replica Set method to experimental and measurement error. We find no compromise in statistical power for Replica Set experiments, even for moderate effect sizes, or when simulated experimental errors are introduced. We compare and contrast the statistical analysis of data generated by the two approaches, and highlight pitfalls common with the traditional methodology. Collectively, our analysis provides a standard of measure for each method across comparable parameters, which will be invaluable in both experimental design and evaluation of published data for lifespan studies.

WikiNetworks: translating manually created biological pathways for topological analysis.

SUMMARY: WikiPathways is a database of 2979 biological pathways across 31 species created using the drawing software PathVisio. Many of these pathways are not directly usable for network-based topological analyses due to differences in curation styles and drawings. We developed the WikiNetworks package to standardize and construct directed networks by combining geometric information and manual annotations from WikiPathways. WikiNetworks performs significantly better than existing tools. This enables the use of high-quality WikiPathways resource for network-based topological analysis of high-throughput data. AVAILABILITY AND IMPLEMENTATION: WikiNetworks is written in Python3 and is available on github.com/Thakar-Lab/wikinetworks and on PyPI. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A systems genomics approach uncovers molecular associates of RSV severity.

Respiratory syncytial virus (RSV) infection results in millions of hospitalizations and thousands of deaths each year. Variations in the adaptive and innate immune response appear to be associated with RSV severity. To investigate the host response to RSV infection in infants, we performed a systems-level study of RSV pathophysiology, incorporating high-throughput measurements of the peripheral innate and adaptive immune systems and the airway epithelium and microbiota. We implemented a novel multi-omic data integration method based on multilayered principal component analysis, penalized regression, and feature weight back-propagation, which enabled us to identify cellular pathways associated with RSV severity. In both airway and immune cells, we found an association between RSV severity and activation of pathways controlling Th17 and acute phase response signaling, as well as inhibition of B cell receptor signaling. Dysregulation of both the humoral and mucosal response to RSV may play a critical role in determining illness severity.

Traditional Farming Lifestyle in Old Older Mennonites Modulates Human Milk Composition.

Background: In addition to farming exposures in childhood, maternal farming exposures provide strong protection against allergic disease in their children; however, the effect of farming lifestyle on human milk (HM) composition is unknown. Objective: This study aims to characterize the maternal immune effects of Old Order Mennonite (OOM) traditional farming lifestyle when compared with Rochester (ROC) families at higher risk for asthma and allergic diseases using HM as a proxy. Methods: HM samples collected at median 2 months of lactation from 52 OOM and 29 ROC mothers were assayed for IgA1 and IgA2 antibodies, cytokines, endotoxin, HM oligosaccharides (HMOs), and targeted fatty acid (FA) metabolites. Development of early childhood atopic diseases in children by 3 years of age was assessed. In addition to group comparisons, systems level network analysis was performed to identify communities of multiple HM factors in ROC and OOM lifestyle. Results: HM contains IgA1 and IgA2 antibodies broadly recognizing food, inhalant, and bacterial antigens. OOM HM has significantly higher levels of IgA to peanut, ovalbumin, dust mites, and Streptococcus equii as well TGF-ß2, and IFN-λ3. A strong correlation occurred between maternal antibiotic use and levels of several HMOs. Path-based analysis of HMOs shows lower activity in the path involving lactoneohexaose (LNH) in the OOM as well as higher levels of lacto-N-neotetraose (LNnT) and two long-chain FAs C-18OH (stearic acid) and C-23OH (tricosanoic acid) compared with Rochester HM. OOM and Rochester milk formed five different clusters, e.g., butyrate production was associated with Prevotellaceae, Veillonellaceae, and Micrococcaceae cluster. Development of atopic disease in early childhood was more common in Rochester and associated with lower levels of total IgA, IgA2 to dust mite, as well as of TSLP. Conclusion: Traditional, agrarian lifestyle, and antibiotic use are strong regulators of maternally derived immune and metabolic factors, which may have downstream implications for postnatal developmental programming of infant's gut microbiome and immune system.

SysMod: the ISCB community for data-driven computational modelling and multi-scale analysis of biological systems.

Computational models of biological systems can exploit a broad range of rapidly developing approaches, including novel experimental approaches, bioinformatics data analysis, emerging modelling paradigms, data standards and algorithms. A discussion about the most recent advances among experts from various domains is crucial to foster data-driven computational modelling and its growing use in assessing and predicting the behaviour of biological systems. Intending to encourage the development of tools, approaches and predictive models, and to deepen our understanding of biological systems, the Community of Special Interest (COSI) was launched in Computational Modelling of Biological Systems (SysMod) in 2016. SysMod's main activity is an annual meeting at the Intelligent Systems for Molecular Biology (ISMB) conference, which brings together computer scientists, biologists, mathematicians, engineers, computational and systems biologists. In the five years since its inception, SysMod has evolved into a dynamic and expanding community, as the increasing number of contributions and participants illustrate. SysMod maintains several online resources to facilitate interaction among the community members, including an online forum, a calendar of relevant meetings and a YouTube channel with talks and lectures of interest for the modelling community. For more than half a decade, the growing interest in computational systems modelling and multi-scale data integration has inspired and supported the SysMod community. Its members get progressively more involved and actively contribute to the annual COSI meeting and several related community workshops and meetings, focusing on specific topics, including particular techniques for computational modelling or standardisation efforts.

Infant gut microbiome is enriched with Bifidobacterium longum ssp. infantis in Old Order Mennonites with traditional farming lifestyle.

BACKGROUND: Growing up on traditional, single-family farms is associated with protection against asthma in school age, but the mechanisms against early manifestations of atopic disease are largely unknown. We sought determine the gut microbiome and metabolome composition in rural Old Order Mennonite (OOM) infants at low risk and Rochester, NY urban/suburban infants at high risk for atopic diseases. METHODS: In a cohort of 65 OOM and 39 Rochester mother-infant pairs, 101 infant stool and 61 human milk samples were assessed by 16S rRNA gene sequencing for microbiome composition and qPCR to quantify Bifidobacterium spp. and B. longum ssp. infantis (B. infantis), a consumer of human milk oligosaccharides (HMOs). Fatty acids (FAs) were analyzed in 34 stool and human 24 milk samples. Diagnoses and symptoms of atopic diseases by 3 years of age were assessed by telephone. RESULTS: At a median age of 2 months, stool was enriched with Bifidobacteriaceae, Clostridiaceae, and Aerococcaceae in the OOM compared with Rochester infants. B. infantis was more abundant (p < .001) and prevalent, detected in 70% of OOM compared with 21% of Rochester infants (p < .001). Stool colonized with B. infantis had higher levels of lactate and several medium- to long/odd-chain FAs. In contrast, paired human milk was enriched with a distinct set of FAs including butyrate. Atopic diseases were reported in 6.5% of OOM and 35% of Rochester children (p < .001). CONCLUSION: A high rate of B. infantis colonization, similar to that seen in developing countries, is found in the OOM at low risk for atopic diseases.