A Note on Ising Network Analysis with Missing Data.

The Ising model has become a popular psychometric model for analyzing item response data. The statistical inference of the Ising model is typically carried out via a pseudo-likelihood, as the standard likelihood approach suffers from a high computational cost when there are many variables (i.e., items). Unfortunately, the presence of missing values can hinder the use of pseudo-likelihood, and a listwise deletion approach for missing data treatment may introduce a substantial bias into the estimation and sometimes yield misleading interpretations. This paper proposes a conditional Bayesian framework for Ising network analysis with missing data, which integrates a pseudo-likelihood approach with iterative data imputation. An asymptotic theory is established for the method. Furthermore, a computationally efficient Pólya-Gamma data augmentation procedure is proposed to streamline the sampling of model parameters. The method's performance is shown through simulations and a real-world application to data on major depressive and generalized anxiety disorders from the National Epidemiological Survey on Alcohol and Related Conditions (NESARC).

Establishing Virtual Bioequivalence and Clinically Relevant Specifications for Omeprazole Enteric-Coated Capsules by Incorporating Dissolution Data in PBPK Modeling.

Currently, Biopharmaceutics Classification System (BCS) classes I and III are the only biological exemptions of immediate-release solid oral dosage forms eligible for regulatory approval. However, through virtual bioequivalence (VBE) studies, BCS class II drugs may qualify for biological exemptions if reliable and validated modeling is used. Here, we sought to establish physiologically based pharmacokinetic (PBPK) models, in vitro-in vivo relationship (IVIVR), and VBE models for enteric-coated omeprazole capsules, to establish a clinically-relevant dissolution specification (CRDS) for screening BE and non-BE batches, and to ultimately develop evaluation criteria for generic omeprazole enteric-coated capsules. To establish omeprazole's IVIVR based on the PBPK model, we explored its in vitro dissolution conditions and then combined in vitro dissolution profile studies with in vivo clinical trials. The predicted omeprazole pharmacokinetics (PK) profiles and parameters closely matched the observed PK data. Based on the VBE results, the bioequivalence study of omeprazole enteric-coated capsules required at least 48 healthy Chinese subjects. Based on the CRDS, the capsules' in vitro dissolution should not be < 28%-54%, < 52%, or < 80% after two, three, and six hours, respectively. Failure to meet these dissolution criteria may result in non-bioequivalence. Here, PBPK modeling and IVIVR methods were used to bridge the in vitro dissolution of the drug with in vivo PK to establish the BE safety space of omeprazole enteric-coated capsules. The strategy used in this study can be applied in BE studies of other BCS II generics to obtain biological exemptions and accelerate drug development.

Evo-devo applied to sleep research: an approach whose time has come.

Sleep occurs in all animals but its amount, form, and timing vary considerably between species and between individuals. Currently, little is known about the basis for these differences, in part, because we lack a complete understanding of the brain circuitry controlling sleep-wake states and markers for the cell types which can identify similar circuits across phylogeny. Here, I explain the utility of an "Evo-devo" approach for comparative studies of sleep regulation and function as well as for sleep medicine. This approach focuses on the regulation of evolutionary ancient transcription factors which act as master controllers of cell-type specification. Studying these developmental transcription factor cascades can identify novel cell clusters which control sleep and wakefulness, reveal the mechanisms which control differences in sleep timing, amount, and expression, and identify the timepoint in evolution when different sleep-wake control neurons appeared. Spatial transcriptomic studies, which identify cell clusters based on transcription factor expression, will greatly aid this approach. Conserved developmental pathways regulate sleep in mice, Drosophila, and C. elegans. Members of the LIM Homeobox (Lhx) gene family control the specification of sleep and circadian neurons in the forebrain and hypothalamus. Increased Lhx9 activity may account for increased orexin/hypocretin neurons and reduced sleep in Mexican cavefish. Other transcription factor families specify sleep-wake circuits in the brainstem, hypothalamus, and basal forebrain. The expression of transcription factors allows the generation of specific cell types for transplantation approaches. Furthermore, mutations in developmental transcription factors are linked to variation in sleep duration in humans, risk for restless legs syndrome, and sleep-disordered breathing. This paper is part of the "Genetic and other molecular underpinnings of sleep, sleep disorders, and circadian rhythms including translational approaches" collection.

The molecular and cellular hematopoietic stem cell specification niche.

Hematopoietic stem cells (HSCs) are a population of tissue-specific stem cells that reside in the bone marrow of adult mammals where they self-renew and continuously regenerate the adult hematopoietic lineages over the life of the individual. Prominence as a stem cell model and clinical usefulness has driven interest in understanding the physiological processes that lead to specification of HSCs during embryonic development. High efficiency directed differentiation of HSCs by instruction of defined progenitor cells using sequentially defined instructive molecules and conditions remains impossible, indicating that comprehensive knowledge of the complete set of precursor intermediate identities and required inductive inputs remains incompletely understood. Recently, interest in the molecular and cellular microenvironment where HSCs are specified from endothelial precursors-the "specification niche"-has increased. Here we review recent progress in understanding these niche spaces across vertebrate phyla, as well as how a better characterization of the origin and molecular phenotypes of the niche cell populations has helped inform and complicate previous understanding of signaling required for HSC emergence and maturation.

PHLDA1-PRDM1 mediates the effect of lentiviral vectors on fate-determination of human retinal progenitor cells.

Lentiviral vectors have markedly enhanced gene therapy efficiency in treating congenital diseases, but their long-term safety remains controversial. Most gene therapies for congenital eye diseases need to be carried out at early ages, yet the assessment of related risks to ocular development posed by lentiviral vectors is challenging. Utilizing single-cell transcriptomic profiling on human retinal organoids, this study explored the impact of lentiviral vectors on the retinal development and found that lentiviral vectors can cause retinal precursor cells to shift toward photoreceptor fate through the up-regulation of key fate-determining genes such as PRDM1. Further investigation demonstrated that the intron and intergenic region of PRDM1 was bound by PHLDA1, which was also up-regulated by lentiviral vectors exposure. Importantly, knockdown of PHLDA1 successfully suppressed the lentivirus-induced differentiation bias of photoreceptor cells. The findings also suggest that while lentiviral vectors may disrupt the fate determination of retinal precursor cells, posing risks in early-stage retinal gene therapy, these risks could potentially be reduced by inhibiting the PHLDA1-PRDM1 axis.

Mapping MAVE data for use in human genomics applications.

The large-scale experimental measures of variant functional assays submitted to MaveDB have the potential to provide key information for resolving variants of uncertain significance, but the reporting of results relative to assayed sequence hinders their downstream utility. The Atlas of Variant Effects Alliance mapped multiplexed assays of variant effect data to human reference sequences, creating a robust set of machine-readable homology mappings. This method processed approximately 2.5 million protein and genomic variants in MaveDB, successfully mapping 98.61% of examined variants and disseminating data to resources such as the UCSC Genome Browser and Ensembl Variant Effect Predictor.

Neural Tissue-Like, not Supraphysiological, Electrical Conductivity Stimulates Neuronal Lineage Specification through Calcium Signaling and Epigenetic Modification.

Electrical conductivity is a pivotal biophysical factor for neural interfaces, though optimal values remain controversial due to challenges isolating this cue. To address this issue, conductive substrates made of carbon nanotubes and graphene oxide nanoribbons, exhibiting a spectrum of conductivities from 0.02 to 3.2 S m-1, while controlling other surface properties is designed. The focus is to ascertain whether varying conductivity in isolation has any discernable impact on neural lineage specification. Remarkably, neural-tissue-like low conductivity (0.02-0.1 S m-1) prompted neural stem/progenitor cells to exhibit a greater propensity toward neuronal lineage specification (neurons and oligodendrocytes, not astrocytes) compared to high supraphysiological conductivity (3.2 S m-1). High conductivity instigated the apoptotic process, characterized by increased apoptotic fraction and decreased neurogenic morphological features, primarily due to calcium overload. Conversely, cells exposed to physiological conductivity displayed epigenetic changes, specifically increased chromatin openness with H3acetylation (H3ac) and neurogenic-transcription-factor activation, along with a more balanced intracellular calcium response. The pharmacological inhibition of H3ac further supported the idea that such epigenetic changes might play a key role in driving neuronal specification in response to neural-tissue-like, not supraphysiological, conductive cues. These findings underscore the necessity of optimal conductivity when designing neural interfaces and scaffolds to stimulate neuronal differentiation and facilitate the repair process.

Tlx3 mediates neuronal differentiation and proper condensation of the developing trigeminal ganglion.

The trigeminal ganglion, the largest of the vertebrate cranial ganglia, is comprised of sensory neurons that relay sensations of pain, touch, and temperature to the brain. These neurons are derived from two embryonic cell types, the neural crest and ectodermal placodes, whose interactions are critical for proper ganglion formation. While the T-cell leukemia homeobox 3 (Tlx3) gene is known to be expressed in placodally-derived sensory neurons and necessary for their differentiation, little was known about Tlx3 expression and/or function in the neural crest-derived component of the developing trigeminal ganglion. By combining lineage labeling with in situ hybridization in the chick embryo, we show that neural crest-derived cells that contribute to the cranial trigeminal ganglion express Tlx3 at a time point that coincides with the onset of ganglion condensation. Importantly, loss of Tlx3 function in vivo diminishes the overall size and abundance of neurons within the trigeminal ganglion. Conversely, ectopic expression of Tlx3 in migrating cranial neural crest results in their premature neuronal differentiation. Taken together, our results demonstrate a critical role for Tlx3 in neural crest-derived cells during chick trigeminal gangliogenesis.

The multiverse of universes: A tutorial to plan, execute and interpret multiverses analyses using the R package multiverse.

Even when guided by strong theories and sound methods, researchers must often choose a singular course of action from multiple viable alternatives. Regardless of the choice, it, along with all other choices made during the research process, individually and collectively affects study results, often in unpredictable ways. The inability to disentangle how much of an observed effect is attributable to the phenomenon of interest, and how much is attributable to what have come to be known as researcher degrees of freedom (RDF), slows theoretical progress and stymies practical implementation. However, if one could examine the results from a particular set of RDF (known as a universe) against a systematically and comprehensively determined background of alternative viable universes (known as a multiverse), then the effects of RDF can be directly examined to provide greater context and clarity to future researchers, and greater confidence in the recommendations to practitioners. This tutorial demonstrates a means to map result variability directly and efficiently, and empirically investigate RDF impact on conclusions via multiverse analysis. Using the R package multiverse, we outline best practices in planning, executing and interpreting of multiverse analyses.

Ebf1 is a mouse deafness gene and deletion causes a dramatic increase in hair cells and support cells of the organ of Corti.

Following up on our previous observation that early B cell factor (EBF) sites are enriched in open chromatin of the developing sensory epithelium of the mouse cochlea, we investigated the effect of deletion of Ebf1 on inner ear development. We used a Cre driver to delete Ebf1 at the otocyst stage prior to development of the cochlea. We examined the cochlea at postnatal day (P) 1 and found that the sensory epithelium had doubled in size but the length of the cochlear duct was unaffected. We also found that deletion of Ebf1 led to ectopic sensory patches in the Kölliker's organ. Innervation of the developing organ of Corti was disrupted with no obvious spiral bundles. The ectopic patches were also innervated. All the extra hair cells (HCs) within the sensory epithelium and Kölliker's organ contained mechanoelectrical transduction channels as indicated by rapid uptake of FM1-43. The excessive numbers of HCs were still present in the adult Ebf1 conditional knockout (cKO) animal. The animals had no detectable auditory brainstem response (ABR) suggesting that this gene is essential for hearing development.

Revisiting trophectoderm-inner cell mass lineage segregation in the mammalian preimplantation embryo.

In the first days of life, cells of the mammalian embryo segregate into two distinct lineages, trophectoderm and inner cell mass. Unlike nonmammalian species, mammalian development does not proceed from predetermined factors in the oocyte. Rather, asymmetries arise de novo in the early embryo incorporating cues from cell position, contractility, polarity, and cell-cell contacts. Molecular heterogeneities, including transcripts and non-coding RNAs, have now been characterized as early as the 2-cell stage. However, it's debated whether these early heterogeneities bias cells toward one fate or the other or whether lineage identity arises stochastically at the 16-cell stage. This review summarizes what is known about early blastomere asymmetries and our understanding of lineage allocation in the context of historical models. Preimplantation development is reviewed coupled with what is known about changes in morphology, contractility, and transcription factor networks. The addition of single-cell atlases of human embryos has begun to reveal key differences between human and mouse, including the timing of events and core transcription factors. Furthermore, the recent generation of blastoid models will provide valuable tools to test and understand fate determinants. Lastly, new techniques are reviewed, which may better synthesize existing knowledge with emerging data sets and reconcile models with the regulative capacity unique to the mammalian embryo.

Physiologically Based Biopharmaceutics Modeling for Gefapixant IR Formulation Development and Defining the Bioequivalence Dissolution Safe Space.

Gefapixant is a weakly basic drug which has been formulated as an immediate release tablet for oral administration. A physiologically based biopharmaceutics model (PBBM) was developed based on gefapixant physicochemical properties and clinical pharmacokinetics to aid formulation selection, bioequivalence safe space assessment and dissolution specification settings. In vitro dissolution profiles of different free base and citrate salt formulations were used as an input to the model. The model was validated against the results of independent studies, which included a bioequivalence and a relative bioavailability study, as well as a human ADME study, all meeting acceptance criteria of prediction errors ≤ 20% for both Cmax and AUC.  PBBM was also applied to evaluate gastric pH-mediated drug-drug-interaction potential with co-administration of a proton pump inhibitor (PPI), omeprazole. Model results showed good agreement with clinical data in which omeprazole lowered gefapixant exposure for the free base formulation but did not significantly alter gefapixant pharmacokinetics for the citrate based commercial drug product. An extended virtual dissolution bioequivalence safe space was established.  Gefapixant drug product batches are anticipated to be bioequivalent with the clinical reference batch when their dissolution is > 80% in 60 minutes. PBBM established a wide dissolution bioequivalence space as part of assuring product quality.

Lineage specification in glioblastoma is regulated by METTL7B.

Glioblastomas are the most common malignant brain tumors in adults; they are highly aggressive and heterogeneous and show a high degree of plasticity. Here, we show that methyltransferase-like 7B (METTL7B) is an essential regulator of lineage specification in glioblastoma, with an impact on both tumor size and invasiveness. Single-cell transcriptomic analysis of these tumors and of cerebral organoids derived from expanded potential stem cells overexpressing METTL7B reveal a regulatory role for the gene in the neural stem cell-to-astrocyte differentiation trajectory. Mechanistically, METTL7B downregulates the expression of key neuronal differentiation players, including SALL2, via post-translational modifications of histone marks.

A model specification test for semiparametric nonignorable missing data modeling.

The instrumental variable approaches have been demonstrated effective for semiparametrically modeling the propensity function in analyzing data that may be missing not at random. A model specification test is considered for a class of parsimonious semiparametric propensity models. The test is constructed based on assessing an over-identification so as to detect possible incompatibility in the moment conditions when the model and/or instrumental variables are misspecified. Validity of the test under the null hypothesis is established; and its power is studied when the model is misspecified. A data analysis and simulations are presented to demonstrate the effectiveness of our methods.

The diversification of the shoot branching system: A quantitative and comparative perspective in meristem determinacy.

Reiterative shoot branching largely defines important yield components of crops and is essentially controlled by programs that direct the initiation, dormancy release, and differentiation of meristems in the axils of leaves. Here, we focus on meristem determinacy, defining the number of reiterations that shape the shoot architectures and exhibit enormous diversity in a wide range of species. The meristem determinacy per se is hierarchically complex and context-dependent for the successively emerged meristems, representing a crucial mechanism in shaping the complexity of the shoot branching. In addition, we have highlighted that two key components of axillary meristem developmental programs may have been co-opted in controlling flower/ear number of an axillary inflorescence in legumes/maize, hinting at the diversification of axillary-meristem-patterning programs in different lineages. This begs the question how axillary meristem patterning programs may have diversified during plant evolution and hence helped shape the rich variation in shoot branching systems.

Making Ramón y Cajal proud: Development of cell identity and diversity in the cerebral cortex.

Since the beautiful images of Santiago Ramón y Cajal provided a first glimpse into the immense diversity and complexity of cell types found in the cerebral cortex, neuroscience has been challenged and inspired to understand how these diverse cells are generated and how they interact with each other to orchestrate the development of this remarkable tissue. Some fundamental questions drive the field's quest to understand cortical development: what are the mechanistic principles that govern the emergence of neuronal diversity? How do extrinsic and intrinsic signals integrate with physical forces and activity to shape cell identity? How do the diverse populations of neurons and glia influence each other during development to guarantee proper integration and function? The advent of powerful new technologies to profile and perturb cortical development at unprecedented resolution and across a variety of modalities has offered a new opportunity to integrate past knowledge with brand new data. Here, we review some of this progress using cortical excitatory projection neurons as a system to draw out general principles of cell diversification and the role of cell-cell interactions during cortical development.

Metabolic regulator ERRγ governs gastric stem cell differentiation into acid-secreting parietal cells.

Parietal cells (PCs) produce gastric acid to kill pathogens and aid digestion. Dysregulated PC census is common in disease, yet how PCs differentiate is unclear. Here, we identify the PC progenitors arising from isthmal stem cells, using mouse models and human gastric cells, and show that they preferentially express cell-metabolism regulator and orphan nuclear receptor Estrogen-related receptor gamma (Esrrg, encoding ERRγ). Esrrg expression facilitated the tracking of stepwise molecular, cellular, and ultrastructural stages of PC differentiation. EsrrgP2ACreERT2 lineage tracing revealed that Esrrg expression commits progenitors to differentiate into mature PCs. scRNA-seq indicated the earliest Esrrg+ PC progenitors preferentially express SMAD4 and SP1 transcriptional targets and the GTPases regulating acid-secretion signal transduction. As progenitors matured, ERRγ-dependent metabolic transcripts predominated. Organoid and mouse studies validated the requirement of ERRγ for PC differentiation. Our work chronicles stem cell differentiation along a single lineage in vivo and suggests ERRγ as a therapeutic target for PC-related disorders.

Within- and between-subject biological variation estimates for the enumeration of lymphocyte deep immunophenotyping and monocyte subsets.

OBJECTIVES: This study aimed to deliver biological variation (BV) estimates for 25 types of lymphocyte subpopulations subjected to deep immunophenotyping (memory T/B cells, regulatory T cells, etc.) and classical, intermediate, and nonclassical monocyte subsets based on the full spectrum flow cytometry (FS-FCM) and a Biological Variation Data Critical Appraisal Checklist (BIVAC) design. METHODS: Samples were collected biweekly from 60 healthy Chinese adults over 10 consecutive two-week periods. Each sample was measured in duplicate within a single run for lymphocyte deep immunophenotyping and monocyte subset determination using FS-FCM, including the percentage (%) and absolute count (cells/µL). After trend adjustment, a Bayesian model was applied to deliver the within-subject BV (CVI) and between-subject BV (CVG) estimates with 95 % credibility intervals. RESULTS: Enumeration (% and cells/µL) for 25 types of lymphocyte deep immunophenotyping and three types of monocyte subset percentages showed considerable variability in terms of CVI and CVG. CVI ranged from 4.23 to 47.47 %. Additionally, CVG ranged between 10.32 and 101.30 %, except for CD4+ effector memory T cells re-expressing CD45RA. No significant differences were found between males and females for CVI and CVG estimates. Nevertheless, the CVGs of PD-1+ T cells (%) may be higher in females than males. Based on the desired analytical performance specification, the maximum allowable imprecision immune parameter was the CD8+PD-1+ T cell (cells/µL), with 23.7 %. CONCLUSIONS: This is the first study delivering BV estimates for 25 types of lymphocyte subpopulations subjected to deep immunophenotyping, along with classical, intermediate, and nonclassical monocyte subsets, using FS-FCM and adhering to the BIVAC design.

Regulatory Considerations for Producing mRNA Vaccines for Clinical Trials.

The approval of clinical trials by the competent authorities requires comprehensive quality documentation on the new drug to be used on the clinical trial participant. In the EU, quality data is summarized as investigational medicinal product dossier (IMPD), in the United States, as investigational new drug (IND) application. For that, several preconditions concerning production, quality control, and assurance have to be fulfilled. Here, specific requirements related to mRNA vaccines are addressed on the basis of European standards.

Specification curve analysis to identify heterogeneity in risk factors for dementia: findings from the UK Biobank.

BACKGROUND: In 2020, the Lancet Commission identified 12 risk factors as priorities for prevention of dementia, and other studies identified APOE e4/e4 genotype and family history of Alzheimer's disease strongly associated with dementia outcomes; however, it is unclear how robust these relationships are across dementia subtypes and analytic scenarios. Specification curve analysis (SCA) is a new tool to probe how plausible analytical scenarios influence outcomes. METHODS: We evaluated the heterogeneity of odds ratios for 12 risk factors reported from the Lancet 2020 report and two additional strong associated non-modifiable factors (APOE e4/e4 genotype and family history of Alzheimer's disease) with dementia outcomes across 450,707 UK Biobank participants using SCA with 5357 specifications across dementia subtypes (outcomes) and analytic models (e.g., standard demographic covariates such as age or sex and/or 14 correlated risk factors). RESULTS: SCA revealed variable dementia risks by subtype and age, with associations for TBI and APOE e4/e4 robust to model specification; in contrast, diabetes showed fluctuating links with dementia subtypes. We found that unattributed dementia participants had similar risk factor profiles to participants with defined subtypes. CONCLUSIONS: We observed heterogeneity in the risk of dementia, and estimates of risk were influenced by the inclusion of a combination of other modifiable risk factors; non-modifiable demographic factors had a minimal role in analytic heterogeneity. Future studies should report multiple plausible analytic scenarios to test the robustness of their association. Considering these combinations of risk factors could be advantageous for the clinical development and evaluation of novel screening models for different types of dementia.