Through-Ice Acoustic Communication for Ocean Worlds Exploration.

Subsurface exploration of ice-covered planets and moons presents communications challenges because of the need to communicate through kilometers of ice. The objective of this task is to develop the capability to wirelessly communicate through kilometers of ice and thus complement the potentially failure-prone tethers deployed behind an ice-penetrating probe on Ocean Worlds. In this paper, the preliminary work on the development of wireless deep-ice communication is presented and discussed. The communication test and acoustic attenuation measurements in ice have been made by embedding acoustic transceivers in glacial ice at the Matanuska Glacier, Anchorage, Alaska. Field test results show that acoustic communication is viable through ice, demonstrating the transmission of data and image files in the 13-18 kHz band over 100 m. The results suggest that communication over many kilometers of ice thickness could be feasible by employing reduced transmitting frequencies around 1 kHz, though future work is needed to better constrain the likely acoustic attenuation properties through a refrozen borehole.

Generation and characterization of an immunodeficient mouse model of mucopolysaccharidosis type II.

Mucopolysaccharidosis type II (Hunter syndrome, MPS II) is an inherited X-linked recessive disease caused by deficiency of iduronate-2-sulfatase (IDS), resulting in the accumulation of the glycosaminoglycans (GAG) heparan and dermatan sulfates. Mouse models of MPS II have been used in several reports to study disease pathology and to conduct preclinical studies for current and next generation therapies. Here, we report the generation and characterization of an immunodeficient mouse model of MPS II, where CRISPR/Cas9 was employed to knock out a portion of the murine IDS gene on the NOD/SCID/Il2rγ (NSG) immunodeficient background. IDS-/- NSG mice lacked detectable IDS activity in plasma and all analyzed tissues and exhibited elevated levels of GAGs in those same tissues and in the urine. Histopathology revealed vacuolized cells in both the periphery and CNS of NSG-MPS II mice. This model recapitulates skeletal disease manifestations, such as increased zygomatic arch diameter and decreased femur length. Neurocognitive deficits in spatial memory and learning were also observed in the NSG-MPS II model. We anticipate that this new immunodeficient model will be appropriate for preclinical studies involving xenotransplantation of human cell products intended for the treatment of MPS II.

A novel process for H&E, immunofluorescence, and imaging mass cytometry on a single slide with a concise analytics pipeline.

Imaging mass cytometry (IMC) is a powerful spatial technology that utilizes cytometry time of flight to acquire multiplexed image datasets with up to 40 markers, via metal-tagged antibodies. Recent advances in IMC have led to the inclusion of RNAScope probes and multiple new analysis pipelines have led to faster analyses and better results. However, IMC still suffers from lower resolution (1 µm2 pixels) and relatively small regions of interest (ROIs) (<2 mm2 ) compared to other, light-based microscope technologies. Capturing higher-resolution images on serial sections causes great difficulty when attempting to align cells and structures across serial sections, especially when observing smaller cell types and structures. Therefore, we demonstrate the combination of H&E and multiplex immunofluorescence imaging, for much higher resolution of the structural and cellular compartments found throughout the entire tissue section, with the high-dimensionality of IMC for specific ROIs on a single slide. Additionally, we demonstrate a simple and effective open-source cell segmentation and IMC analysis pipeline with previously published and freely available software.

Contribution of the innate and adaptive immune systems to aortic dilation in murine mucopolysaccharidosis type I.

BACKGROUND: Adult immunocompetent male C57Bl/6 mucopolysaccharidosis, type I (MPSI) mice develop aortic insufficiency (AI), dilated ascending aortas and decreased cardiac function, findings not observed in immune incompetent adult male NSG MPSI mice. We sought to determine why. METHODS: Cardiac ultrasound measurements of ascending aorta and left ventricular dimensions and Doppler interrogation for AI were performed in 6-month-old male B6 MPSI (N = 12), WT (N = 6), NSG MPSI (N = 8), NSG (N = 6) mice. Urinary glycosaminoglycans, RNA sequencing with quantitative PCR were performed and aortic pathology assessed by routine and immunohistochemical staining on subsets of murine aortas. RESULTS: Ascending aortic diameters were significantly greater, left ventricular function significantly decreased, and AI significantly more frequent in B6 MPSI mice compared to NSG MPSI mice (p < 0.0001, p = 0.008 and p = 0.02, respectively); NSG and B6 WT mice showed no changes. Urinary glycosaminoglycans were significantly greater in B6 and NSG MPSI mice and both were significantly elevated compared to WT controls (p = 0.003 and p < 0.0001, respectively). By RNA sequencing, all 11 components of the inflammasome pathway were upregulated in B6 MUT, but only Aim2 and Ctsb in NSG MUT mice and none in WT controls. Both B6 and NSG MUT mice demonstrated variably-severe intramural inflammation, vacuolated cells, elastin fragmentation and disarray, and intense glycosaminoglycans on histological staining. B6 MPSI mice demonstrated numerous medial MAC2+ macrophages and adventitial CD3+ T-cells while MAC2+ macrophages were sparse and CD3+ T-cells absent in NSG MPSI mice. CONCLUSIONS: Aortic dilation, AI and decreased cardiac function occur in immunocompetent B6 MPSI male mice but not in immune incompetent NSG MPSI mice, unrelated to GAG excretion, upregulation of Ctsb, or routine histologic appearance. Upregulation of all components of the inflammasome pathway in B6 MUT, but not NSG MUT mice, and abundant medial MAC2 and adventitial CD3 infiltrates in B6, but not NSG, MPSI aortas differentiated the two strains. These results suggest that the innate and adaptive immune systems play a role in these cardiac findings which may be relevant to human MPSI.

Autoantibody-positive healthy individuals with lower lupus risk display a unique immune endotype.

BACKGROUND: Autoimmune diseases comprise a spectrum of illnesses and are on the rise worldwide. Although antinuclear antibodies (ANAs) are detected in many autoimmune diseases, up to 20% of healthy women are ANA-positive (ANA+) and most will never develop clinical symptoms. Furthermore, disease transition is higher among ANA+ African Americans compared with ANA+ European Americans. OBJECTIVE: We sought to determine the immune features that might define and prevent transition to clinical autoimmunity in ANA+ healthy individuals. METHODS: We comprehensively phenotyped immune profiles of African Americans and European Americans who are ANA-negative (ANA-) healthy, ANA+ healthy, or have SLE using single cell mass cytometry, next-generation RNA-sequencing, multiplex cytokine profiling, and phospho-signaling analyses. RESULTS: We found that, compared with both ANA- and ANA+ healthy individuals, patients with SLE of both races displayed T-cell expansion and elevated expression of type I and II interferon pathways. We discovered a unique immune signature that suggests a suppressive immune phenotype and reduced CD11C+ autoimmunity-associated B cells in healthy ANA+ European Americans that is absent in their SLE or even healthy ANA- counterparts, or among African American cohorts. In contrast, ANA+ healthy African Americans exhibited elevated expression of T-cell activation markers and higher plasma levels of IL-6 than did healthy ANA+ European Americans. CONCLUSIONS: We propose that this novel immune signature identified in ANA+ healthy European Americans may protect them from T-cell expansion, heightened activation of interferon pathways, and disease transition.

Unique Sjögren's syndrome patient subsets defined by molecular features.

OBJECTIVE: To address heterogeneity complicating primary SS (pSS) clinical trials, research and care by characterizing and clustering patients by their molecular phenotypes. METHODS: pSS patients met American-European Consensus Group classification criteria and had at least one systemic manifestation and stimulated salivary flow of ⩾0.1 ml/min. Correlated transcriptional modules were derived from gene expression microarray data from blood (n = 47 with appropriate samples). Patients were clustered based on this molecular information using an unbiased random forest modelling approach. In addition, multiplex, bead-based assays and ELISAs were used to assess 30 serum cytokines, chemokines and soluble receptors. Eleven autoantibodies, including anti-Ro/SSA and anti-La/SSB, were measured by Bio-Rad Bioplex 2200. RESULTS: Transcriptional modules distinguished three clusters of pSS patients. Cluster 1 showed no significant elevation of IFN or inflammation modules. Cluster 2 showed strong IFN and inflammation modular network signatures, as well as high plasma protein levels of IP-10/CXCL10, MIG/CXCL9, BLyS (BAFF) and LIGHT. Cluster 3 samples exhibited moderately elevated IFN modules, but with suppressed inflammatory modules, increased IP-10/CXCL10 and B cell-attracting chemokine 1/CXCL13 and trends toward increased MIG/CXCL9, IL-1α, and IL-21. Anti-Ro/SSA and anti-La/SSB were present in all three clusters. CONCLUSION: Molecular profiles encompassing IFN, inflammation and other signatures can be used to separate patients with pSS into distinct clusters. In the future, such profiles may inform patient selection for clinical trials and guide treatment decisions.

Parameter estimation for a point-source diffusion-decay morphogen model.

In this paper we present a novel method for finding unknown parameters for an unknown morphogen. We postulate the existence of an unknown morphogen in a given three-dimensional domain due to the spontaneous arrangement of a downstream species on the domain boundary for which data is known. Assuming a modified Helmholtz model for the morphogen and that it is produced from a single source in the domain, our method accurately estimates the source location and other model parameters. Notably, our method does not require the forward solution of the model to be computed which can often be a challenge for three-dimensional PDE model parameter fitting. Instead, an extension is made from the problem domain to an infinite domain and the analytic nature of the fundamental solution is exploited. We explore in this manuscript strategies for best conditioning the problem and rigorously explore the accuracy of the method on two test problems. Our tests focus on the effect of source location on accuracy but also the robustness of the algorithm to experimental noise.

Generating New Space-Filling Test Instances for Continuous Black-Box Optimization.

This article presents a method to generate diverse and challenging new test instances for continuous black-box optimization. Each instance is represented as a feature vector of exploratory landscape analysis measures. By projecting the features into a two-dimensional instance space, the location of existing test instances can be visualized, and their similarities and differences revealed. New instances are generated through genetic programming which evolves functions with controllable characteristics. Convergence to selected target points in the instance space is used to drive the evolutionary process, such that the new instances span the entire space more comprehensively. We demonstrate the method by generating two-dimensional functions to visualize its success, and ten-dimensional functions to test its scalability. We show that the method can recreate existing test functions when target points are co-located with existing functions, and can generate new functions with entirely different characteristics when target points are located in empty regions of the instance space. Moreover, we test the effectiveness of three state-of-the-art algorithms on the new set of instances. The results demonstrate that the new set is not only more diverse than a well-known benchmark set, but also more challenging for the tested algorithms. Hence, the method opens up a new avenue for developing test instances with controllable characteristics, necessary to expose the strengths and weaknesses of algorithms, and drive algorithm development.

Performance Analysis of Continuous Black-Box Optimization Algorithms via Footprints in Instance Space.

This article presents a method for the objective assessment of an algorithm's strengths and weaknesses. Instead of examining the performance of only one or more algorithms on a benchmark set, or generating custom problems that maximize the performance difference between two algorithms, our method quantifies both the nature of the test instances and the algorithm performance. Our aim is to gather information about possible phase transitions in performance, that is, the points in which a small change in problem structure produces algorithm failure. The method is based on the accurate estimation and characterization of the algorithm footprints, that is, the regions of instance space in which good or exceptional performance is expected from an algorithm. A footprint can be estimated for each algorithm and for the overall portfolio. Therefore, we select a set of features to generate a common instance space, which we validate by constructing a sufficiently accurate prediction model. We characterize the footprints by their area and density. Our method identifies complementary performance between algorithms, quantifies the common features of hard problems, and locates regions where a phase transition may lie.

Approximate Bayesian computation schemes for parameter inference of discrete stochastic models using simulated likelihood density.

BACKGROUND: Mathematical modeling is an important tool in systems biology to study the dynamic property of complex biological systems. However, one of the major challenges in systems biology is how to infer unknown parameters in mathematical models based on the experimental data sets, in particular, when the data are sparse and the regulatory network is stochastic. RESULTS: To address this issue, this work proposed a new algorithm to estimate parameters in stochastic models using simulated likelihood density in the framework of approximate Bayesian computation. Two stochastic models were used to demonstrate the efficiency and effectiveness of the proposed method. In addition, we designed another algorithm based on a novel objective function to measure the accuracy of stochastic simulations. CONCLUSIONS: Simulation results suggest that the usage of simulated likelihood density improves the accuracy of estimates substantially. When the error is measured at each observation time point individually, the estimated parameters have better accuracy than those obtained by a published method in which the error is measured using simulations over the entire observation time period.

A note on the relationship between turbulent coherent structures and phase correlation.

Various definitions of coherent structures exist in turbulence research, but a common assumption is that coherent structures have correlated spectral phases. As a result, randomization of phases is believed, generally, to remove coherent structures from the measured data. Here, we reexamine these assumptions using atmospheric turbulence measurements. Small-scale coherent structures are detected in the usual way using the wavelet transform. A considerable percentage of the detected structures are not phase correlated, although some of them are clearly organized in space and time. At larger scales, structures have even higher degree of spatiotemporal coherence but are also associated with weak phase correlation. A series of specific examples are shown to demonstrate this. These results warn about the vague terminology and assumptions around coherent structures, particularly for complex real-world turbulence.

Comparative dose effectiveness of intravenous and intrathecal AAV9.CB7.hIDS, RGX-121, in mucopolysaccharidosis type II mice.

Mucopolysaccharidosis type II (MPS II) is an X-linked recessive lysosomal disease caused by iduronate-2-sulfatase (IDS) deficiency, leading to accumulation of glycosaminoglycans (GAGs) and the emergence of progressive disease. Enzyme replacement therapy is the only currently approved treatment, but it leaves neurological disease unaddressed. Cerebrospinal fluid (CSF)-directed administration of AAV9.CB7.hIDS (RGX-121) is an alternative treatment strategy, but it is unknown if this approach will affect both neurologic and systemic manifestations. We compared the effectiveness of intrathecal (i.t.) and intravenous (i.v.) routes of administration (ROAs) at a range of vector doses in a mouse model of MPS II. While lower doses were completely ineffective, a total dose of 1 × 109 gc resulted in appreciable IDS activity levels in plasma but not tissues. Total doses of 1 × 1010 and 1 × 1011 gc by either ROA resulted in supraphysiological plasma IDS activity, substantial IDS activity levels and GAG reduction in nearly all tissues, and normalized zygomatic arch diameter. In the brain, a dose of 1 × 1011 gc i.t. achieved the highest IDS activity levels and the greatest reduction in GAG content, and it prevented neurocognitive deficiency. We conclude that a dose of 1 × 1010 gc normalized metabolic and skeletal outcomes, while neurologic improvement required a dose of 1 × 1011 gc, thereby suggesting the prospect of a similar direct benefit in humans.

Strategies in Infertile Azoospermic Patients with Negative Microdissection Testicular Sperm Extraction Surgery.

Non-obstructive azoospermia is reported to affect 1 in 100 men, and despite advances in surgical practice, the succesful sperm retrieval rate for microdissection testicular sperm extraction surgery (mTESE) is only 46%. This article reviews the potential causes for mTESE failure and provides a management strategy to guide the clinicians on how to treat this challenging cohort of patients.

Genetic load in incomplete lupus erythematosus.

OBJECTIVE: Patients with incomplete lupus erythematosus (ILE) have lupus features but insufficient criteria for SLE classification. Some patients with ILE transition to SLE, but most avoid major organ involvement. This study tested whether the milder disease course in ILE is influenced by reduced SLE risk allele genetic load. METHODS: We calculated the genetic load based on 99 SLE-associated risk alleles in European American patients with SLE (≥4 American College of Rheumatology (ACR) 1997 criteria, n=170), patients with ILE (3 ACR 1997 criteria, n=169), a subset of patients with ILE not meeting Systemic Lupus International Collaborating Clinics (SLICC) classification (ILESLICC, n=119) and healthy controls (n=133). Unweighted genetic loads were calculated as the total sum of risk alleles for each individual, while weighted genetic loads were defined as the sum of risk alleles multiplied by the natural logarithm of the previously published OR of each risk allele for SLE susceptibility. RESULTS: The median unweighted and weighted SLE risk allele genetic load was significantly greater in patients with ILE (unweighted: 81, p value=0.01; weighted: 16.3, p value=0.001) and patients with SLE (80, p value=0.02; 16.29, p value=0.0006) compared with healthy controls (78, 15.76). Patients with ILESLICC trended towards an increased genetic load, although not statistically significant (unweighted: 80, p value=0.14; weighted: 16.05, p value=0.07). However, the median genetic load did not significantly differ between ILE and SLE, and genetic load did not differentiate patients with ILE and SLE (area under the curve=0.51, p=0.78) by receiver operator characteristic analysis. CONCLUSIONS: Patients with ILE and SLE have comparable genetic loads of SLE risk loci, suggesting similar genetic predispositions between these conditions. Phenotypical differences between SLE and ILE may instead be influenced by ILE-specific variants and gene-environment interactions.

Neurologic Recovery in MPS I and MPS II Mice by AAV9-Mediated Gene Transfer to the CNS After the Development of Cognitive Dysfunction.

The mucopolysaccharidoses (MPS) are a group of recessively inherited conditions caused by deficiency of lysosomal enzymes essential to the catabolism of glycosaminoglycans (GAG). MPS I is caused by deficiency of the lysosomal enzyme alpha-L-iduronidase (IDUA), while MPS II is caused by a lack of iduronate-2-sulfatase (IDS). Lack of these enzymes leads to early mortality and morbidity, often including neurological deficits. Enzyme replacement therapy has markedly improved the quality of life for MPS I and MPS II affected individuals but is not effective in addressing neurologic manifestations. For MPS I, hematopoietic stem cell transplant has shown effectiveness in mitigating the progression of neurologic disease when carried out in early in life, but neurologic function is not restored in patients transplanted later in life. For both MPS I and II, gene therapy has been shown to prevent neurologic deficits in affected mice when administered early, but the effectiveness of treatment after the onset of neurologic disease manifestations has not been characterized. To test if neurocognitive function can be recovered in older animals, human IDUA or IDS-encoding AAV9 vector was administered by intracerebroventricular injection into MPS I and MPS II mice, respectively, after the development of neurologic deficit. Vector sequences were distributed throughout the brains of treated animals, associated with high levels of enzyme activity and normalized GAG storage. Two months after vector infusion, treated mice exhibited spatial navigation and learning skills that were normalized, that is, indistinguishable from those of normal unaffected mice, and significantly improved compared to untreated, affected animals. We conclude that cognitive function was restored by AAV9-mediated, central nervous system (CNS)-directed gene transfer in the murine models of MPS I and MPS II, suggesting that gene transfer may result in neurodevelopment improvements in severe MPS I and MPS II when carried out after the onset of cognitive decline.

Empirical evidence of climate change and urbanization impacts on warming stream temperatures.

Climate change and urbanization threaten streams and the biodiversity that rely upon them worldwide. Emissions of greenhouse gases are causing air and sea surface temperatures to increase, and even small areas of urbanization are degrading stream biodiversity, water quality and hydrology. However, empirical evidence of how increasing air temperatures and urbanization together affect stream temperatures over time and their relative influence on stream temperatures is limited. This study quantifies changes in stream temperatures in a region in South-East Australia with an urban-agricultural-forest landcover gradient and where increasing air temperatures have been observed. Using Random Forest models we identify air temperature and urbanization drive increasing stream temperatures and that their combined effects are larger than their individual effects occurring alone. Furthermore, we identify potential mitigation measures useful for waterway managers and policy makers. The results show that both local and global solutions are needed to reduce future increases to stream temperature.

Non-invasive intravenous administration of AAV9 transducing iduronate sulfatase leads to global metabolic correction and prevention of neurologic deficits in a mouse model of Hunter syndrome.

Hunter syndrome is a rare x-linked recessive genetic disorder that affects lysosomal metabolism due to deficiency of iduronate-2-sulfatase (IDS), with subsequent accumulation of glycosaminoglycans heparan and dermatan sulfates (GAG). Enzyme replacement therapy is the only FDA-approved remedy and is an expensive life-time treatment that alleviates some symptoms of the disease without neurocognitive benefit. We previously reported successful treatment in a mouse model of mucopolysaccharidosis type II (MPS II) using adeno-associated viral vector serotype 9 encoding human IDS (AAV9.hIDS) via intracerebroventricular injection. As a less invasive and more straightforward procedure, here we report intravenously administered AAV9.hIDS in a mouse model of MPS II. In animals administered 1.5 × 1012 vg of AAV9.hIDS at 2 months of age, we observed supraphysiological levels of IDS enzyme activity in the circulation (up to 9100-fold higher than wild-type), in the tested peripheral organs (up to 560-fold higher than wild-type), but only 4% to 50% of wild type levels in the CNS. GAG levels were normalized on both sides of the blood-brain-barrier (BBB) in most of tissues tested. Despite low levels of the IDS observed in the CNS, this treatment prevented neurocognitive decline as shown by testing in the Barnes maze and by fear conditioning. This study demonstrates that a single dose of IV-administered AAV9.hIDS may be an effective and non-invasive procedure to treat MPS II that benefits both sides of the BBB, with implications for potential use of IV-administered AAV9 for other neuronopathic lysosomal diseases.

Obexelimab in Systemic Lupus Erythematosus With Exploration of Response Based on Gene Pathway Co-Expression Patterns: A Double-Blind, Randomized, Placebo-Controlled, Phase 2 Trial.

OBJECTIVE: Obexelimab is an investigational, bifunctional, noncytolytic monoclonal antibody that binds CD19 and FcyRIIb to inhibit B cells, plasmablasts, and plasma cells. This trial evaluated the efficacy and safety of obexelimab in the treatment of patients with systemic lupus erythematosus (SLE). METHODS: During screening, patients with active, non-organ-threatening SLE received corticosteroid injections to ameliorate symptoms while immunosuppressants were withdrawn (≤10 mg/day prednisone equivalent and ≤400 mg/day hydroxychloroquine allowed). Patients with improved disease activity were randomized 1:1 to obexelimab 5 mg/kg intravenously or placebo once every 2 weeks until week 32 or loss of improvement (LOI). RESULTS: In this study, 104 patients were randomized. Analysis of the primary endpoint, proportion of patients reaching week 32 without LOI, used an efficacy-evaluable (EE) population defined as patients who completed the study or withdrew for flare or treatment-related toxicity. This endpoint did not reach statistical significance: 21 of 50 obexelimab-treated patients (42.0%) versus 12 of 42 patients (28.6%) treated with a placebo (P = 0.183). Time to LOI was increased in obexelimab-treated patients versus patients treated with a placebo in the EE (hazard ratio [HR] 0.53, P = 0.025) and intention-to-treat (HR 0.59, P = 0.062) populations. In obexelimab-treated patients, B cells decreased approximately 50%, and trough concentration and inclusion in baseline gene expression clusters with high B cell pathway modules were associated with increased time to LOI. Obexelimab was associated with infusion reactions but was generally safe and well-tolerated. CONCLUSION: Although the primary endpoint was not reached, secondary analysis showed time to LOI was significantly increased in obexelimab-treated patients, and analysis of patient subsets defined by gene expression patterns at baseline suggests a responding subpopulation.

Ancestry-based differences in the immune phenotype are associated with lupus activity.

Systemic lupus erythematosus (SLE) affects 1 in 537 Black women, which is >2-fold more than White women. Black patients develop the disease at a younger age, have more severe symptoms, and have a greater chance of early mortality. We used a multiomics approach to uncover ancestry-associated immune alterations in patients with SLE and healthy controls that may contribute biologically to disease disparities. Cell composition, signaling, epigenetics, and proteomics were evaluated by mass cytometry; droplet-based single-cell transcriptomics and proteomics; and bead-based multiplex soluble mediator levels in plasma. We observed altered whole blood frequencies and enhanced activity in CD8+ T cells, B cells, monocytes, and DCs in Black patients with more active disease. Epigenetic modifications in CD8+ T cells (H3K27ac) could distinguish disease activity level in Black patients and differentiate Black from White patient samples. TLR3/4/7/8/9-related gene expression was elevated in immune cells from Black patients with SLE, and TLR7/8/9 and IFN-α phospho-signaling and cytokine responses were heightened even in immune cells from healthy Black control patients compared with White individuals. TLR stimulation of healthy immune cells recapitulated the ancestry-associated SLE immunophenotypes. This multiomic resource defines ancestry-associated immune phenotypes that differ between Black and White patients with SLE, which may influence the course and severity of SLE and other diseases.

Mathematical modelling of stem cell differentiation: the PU.1-GATA-1 interaction.

The transcription factors PU.1 and GATA-1 are known to be important in the development of blood progenitor cells. Specifically they are thought to regulate the differentiation of progenitor cells into the granulocyte/macrophage lineage and the erythrocyte/megakaryocite lineage. While several mathematical models have been proposed to investigate the interaction between the transcription factors in recent years, there is still debate about the nature of the progenitor state in the dynamical system, and whether the existing models adequately capture new knowledge about the interactions gleaned from experimental data. Further, the models utilise different formalisms to represent the genetic regulation, and it appears that the resulting dynamical system depends upon which formalism is adopted. In this paper we analyse the four existing models, and propose an alternative model which is shown to demonstrate a rich variety of dynamical systems behaviours found across the existing models, including both bistability and tristability required for modelling the undifferentiated progenitors.