Quantitative fate mapping: A general framework for analyzing progenitor state dynamics via retrospective lineage barcoding.

Natural and induced somatic mutations that accumulate in the genome during development record the phylogenetic relationships of cells; whether these lineage barcodes capture the complex dynamics of progenitor states remains unclear. We introduce quantitative fate mapping, an approach to reconstruct the hierarchy, commitment times, population sizes, and commitment biases of intermediate progenitor states during development based on a time-scaled phylogeny of their descendants. To reconstruct time-scaled phylogenies from lineage barcodes, we introduce Phylotime, a scalable maximum likelihood clustering approach based on a general barcoding mutagenesis model. We validate these approaches using realistic in silico and in vitro barcoding experiments. We further establish criteria for the number of cells that must be analyzed for robust quantitative fate mapping and a progenitor state coverage statistic to assess the robustness. This work demonstrates how lineage barcodes, natural or synthetic, enable analyzing progenitor fate and dynamics long after embryonic development in any organism.

NF-κB subunits direct kinetically distinct transcriptional cascades in antigen receptor-activated B cells.

The nuclear factor kappa B (NF-κB) family of transcription factors orchestrates signal-induced gene expression in diverse cell types. Cellular responses to NF-κB activation are regulated at the level of cell and signal specificity, as well as differential use of family members (subunit specificity). Here we used time-dependent multi-omics to investigate the selective functions of Rel and RelA, two closely related NF-κB proteins, in primary B lymphocytes activated via the B cell receptor. Despite large numbers of shared binding sites genome wide, Rel and RelA directed kinetically distinct cascades of gene expression in activated B cells. Single-cell RNA sequencing revealed marked heterogeneity of Rel- and RelA-specific responses, and sequential binding of these factors was not a major mechanism of protracted transcription. Moreover, nuclear co-expression of Rel and RelA led to functional antagonism between the factors. By rigorously identifying the target genes of each NF-κB subunit, these studies provide insights into exclusive functions of Rel and RelA in immunity and cancer.

TCF-1-Centered Transcriptional Network Drives an Effector versus Exhausted CD8 T Cell-Fate Decision.

TCF-1 is a key transcription factor in progenitor exhausted CD8 T cells (Tex). Moreover, this Tex cell subset mediates responses to PD-1 checkpoint pathway blockade. However, the role of the transcription factor TCF-1 in early fate decisions and initial generation of Tex cells is unclear. Single-cell RNA sequencing (scRNA-seq) and lineage tracing identified a TCF-1+Ly108+PD-1+ CD8 T cell population that seeds development of mature Tex cells early during chronic infection. TCF-1 mediated the bifurcation between divergent fates, repressing development of terminal KLRG1Hi effectors while fostering KLRG1Lo Tex precursor cells, and PD-1 stabilized this TCF-1+ Tex precursor cell pool. TCF-1 mediated a T-bet-to-Eomes transcription factor transition in Tex precursors by promoting Eomes expression and drove c-Myb expression that controlled Bcl-2 and survival. These data define a role for TCF-1 in early-fate-bifurcation-driving Tex precursor cells and also identify PD-1 as a protector of this early TCF-1 subset.

Sequential Enhancer Sequestration Dysregulates Recombination Center Formation at the IgH Locus.

Immunoglobulin heavy-chain (IgH) genes are assembled by DNA rearrangements that juxtapose a variable (VH), a diversity (DH), and a joining (JH) gene segment. Here, we report that in the absence of intergenic control region 1 (IGCR1), the intronic enhancer (Eµ) associates with the next available CTCF binding site located close to VH81X via putative heterotypic interactions involving YY1 and CTCF. The alternate Eµ/VH81X loop leads to formation of a distorted recombination center and altered DH rearrangements and disrupts chromosome conformation that favors distal VH recombination. Cumulatively, these features drive highly skewed, Eµ-dependent recombination of VH81X. Sequential deletion of CTCF binding regions on IGCR1-deleted alleles suggests that they influence recombination of single proximal VH gene segments. Our observations demonstrate that Eµ interacts differently with IGCR1- or VH-associated CTCF binding sites and thereby identify distinct roles for insulator-like elements in directing enhancer activity.

NTR 2.0: a rationally engineered prodrug-converting enzyme with substantially enhanced efficacy for targeted cell ablation.

Transgenic expression of bacterial nitroreductase (NTR) enzymes sensitizes eukaryotic cells to prodrugs such as metronidazole (MTZ), enabling selective cell-ablation paradigms that have expanded studies of cell function and regeneration in vertebrates. However, first-generation NTRs required confoundingly toxic prodrug treatments to achieve effective cell ablation, and some cell types have proven resistant. Here we used rational engineering and cross-species screening to develop an NTR variant, NTR 2.0, which exhibits ~100-fold improvement in MTZ-mediated cell-specific ablation efficacy, eliminating the need for near-toxic prodrug treatment regimens. NTR 2.0 therefore enables sustained cell-loss paradigms and ablation of previously resistant cell types. These properties permit enhanced interrogations of cell function, extended challenges to the regenerative capacities of discrete stem cell niches, and novel modeling of chronic degenerative diseases. Accordingly, we have created a series of bipartite transgenic reporter/effector resources to facilitate dissemination of NTR 2.0 to the research community.

Mobile interspersed repeats are major structural variants in the human genome.

Characterizing structural variants in the human genome is of great importance, but a genome wide analysis to detect interspersed repeats has not been done. Thus, the degree to which mobile DNAs contribute to genetic diversity, heritable disease, and oncogenesis remains speculative. We perform transposon insertion profiling by microarray (TIP-chip) to map human L1(Ta) retrotransposons (LINE-1 s) genome-wide. This identified numerous novel human L1(Ta) insertional polymorphisms with highly variant allelic frequencies. We also explored TIP-chip's usefulness to identify candidate alleles associated with different phenotypes in clinical cohorts. Our data suggest that the occurrence of new insertions is twice as high as previously estimated, and that these repeats are under-recognized as sources of human genomic and phenotypic diversity. We have just begun to probe the universe of human L1(Ta) polymorphisms, and as TIP-chip is applied to other insertions such as Alu SINEs, it will expand the catalog of genomic variants even further.

DNA methylation entropy is associated with DNA sequence features and developmental epigenetic divergence.

Epigenetic information defines tissue identity and is largely inherited in development through DNA methylation. While studied mostly for mean differences, methylation also encodes stochastic change, defined as entropy in information theory. Analyzing allele-specific methylation in 49 human tissue sample datasets, we find that methylation entropy is associated with specific DNA binding motifs, regulatory DNA, and CpG density. Then applying information theory to 42 mouse embryo methylation datasets, we find that the contribution of methylation entropy to time- and tissue-specific patterns of development is comparable to the contribution of methylation mean, and methylation entropy is associated with sequence and chromatin features conserved with human. Moreover, methylation entropy is directly related to gene expression variability in development, suggesting a role for epigenetic entropy in developmental plasticity.

Anti-PF4 antibodies associated with disease severity in COVID-19.

Severe COVID-19 is characterized by a prothrombotic state associated with thrombocytopenia, with microvascular thrombosis being almost invariably present in the lung and other organs at postmortem examination. We evaluated the presence of antibodies to platelet factor 4 (PF4)-polyanion complexes using a clinically validated immunoassay in 100 hospitalized patients with COVID-19 with moderate or severe disease (World Health Organization score, 4 to 10), 25 patients with acute COVID-19 visiting the emergency department, and 65 convalescent individuals. Anti-PF4 antibodies were detected in 95 of 100 hospitalized patients with COVID-19 (95.0%) irrespective of prior heparin treatment, with a mean optical density value of 0.871 ± 0.405 SD (range, 0.177 to 2.706). In contrast, patients hospitalized for severe acute respiratory disease unrelated to COVID-19 had markedly lower levels of the antibodies. In a high proportion of patients with COVID-19, levels of all three immunoglobulin (Ig) isotypes tested (IgG, IgM, and IgA) were simultaneously elevated. Antibody levels were higher in male than in female patients and higher in African Americans and Hispanics than in White patients. Anti-PF4 antibody levels were correlated with the maximum disease severity score and with significant reductions in circulating platelet counts during hospitalization. In individuals convalescent from COVID-19, the antibody levels returned to near-normal values. Sera from patients with COVID-19 induced higher levels of platelet activation than did sera from healthy blood donors, but the results were not correlated with the levels of anti-PF4 antibodies. These results demonstrate that the vast majority of patients with severe COVID-19 develop anti-PF4 antibodies, which may play a role in the clinical complications of COVID-19.

Metabolomic profiles during early childhood and risk of food allergies and asthma in multiethnic children from a prospective birth cohort.

BACKGROUND: There are increasing numbers of metabolomic studies in food allergy (FA) and asthma, which, however, are predominantly limited by cross-sectional designs, small sample size, and being conducted in European populations. OBJECTIVE: We sought to identify metabolites unique to and shared by children with FA and/or asthma in a racially diverse prospective birth cohort, the Boston Birth Cohort. METHODS: Mass spectrometry-based untargeted metabolomic profiling was performed using venous plasma collected in early childhood (n = 811). FA was diagnosed according to clinical symptoms consistent with an acute hypersensitivity reaction at food ingestion and food specific-IgE > 0.35 kU/L. Asthma was defined on the basis of physician diagnosis. Generalized estimating equations were applied to analyze metabolomic associations with FA and asthma, adjusting for potential confounders. RESULTS: During a mean ± standard deviation follow-up of 11.8 ± 5.2 years from birth, 78 children developed FA and 171 developed asthma. Androgenic and pregnenolone steroids were significantly associated with a lower risk of FA, especially for egg allergy. N,N,N-trimethyl-5-aminovalerate (odds ratio [OR] = 0.65, 95% confidence interval [CI] = 0.48-0.87), and 1-oleoyl-2-arachidonoyl-sn-glycero-3-phosphoinositol (OR = 0.77; 95% CI = 0.66-0.90) were inversely associated with FA risk. Orotidine (OR = 4.73; 95% CI = 2.2-10.2) and 4-cholesten-3-one (OR = 0.52; 95% CI = 0.35-0.77) were the top 2 metabolites associated with risk of asthma, although they had no association with FA. In comparison, children with both FA and asthma exhibited an altered metabolomic profile that aligned with that of FA, including altered levels of lipids and steroids. CONCLUSION: In this US multiethnic prospective birth cohort, unique and shared alterations in plasma metabolites during early childhood were associated with risk of developing FA and/or asthma. These findings await further validation.

Emergency myelopoiesis distinguishes multisystem inflammatory syndrome in children from pediatric severe COVID-19.

BACKGROUND: Multisystem inflammatory syndrome in children (MIS-C) is a hyperinflammatory condition caused by recent SARS-CoV-2 infection, but the underlying immunological mechanisms driving this distinct syndrome are unknown. METHODS: We utilized high dimensional flow cytometry, cell-free (cf) DNA, and cytokine and chemokine profiling to identify mechanisms of critical illness distinguishing MIS-C from severe acute COVID-19 (SAC). RESULTS: Compared to SAC, MIS-C patients demonstrated profound innate immune cell death and features of emergency myelopoiesis (EM), an understudied phenomenon observed in severe inflammation. EM signatures were characterized by fewer mature myeloid cells in the periphery and decreased expression of HLA-DR and CD86 on antigen presenting cells. IL-27, a cytokine known to drive hematopoietic stem cells towards EM, was increased in MIS-C, and correlated with immature cell signatures in MIS-C. Upon recovery, EM signatures decreased, and IL-27 plasma levels returned to normal levels. Despite profound lymphopenia, we report a lack of cfDNA released by adaptive immune cells and increased CCR7 expression on T cells indicative of egress out of peripheral blood. CONCLUSIONS: Immune cell signatures of EM combined with elevated innate immune cell-derived cfDNA levels distinguish MIS-C from SAC in children and provide mechanistic insight into dysregulated immunity contributing towards MIS-C, offering potential diagnostic and therapeutic targets.

The BAF chromatin complex component SMARCC1 does not mediate GLI transcriptional repression of Hedgehog target genes in limb buds.

Transcriptional responses to the Hedgehog (HH) signaling pathway are primarily modulated by GLI repression in the mouse limb. Previous studies suggested a role for the BAF chromatin remodeling complex in mediating GLI repression. Consistent with this possibility, the core BAF complex protein SMARCC1 is present at most active limb enhancers including the majority of GLI enhancers. However, in contrast to GLI repression which reduces chromatin accessibility, SMARCC1 maintains chromatin accessibility at most enhancers, including those bound by GLI. Moreover, SMARCC1 binding at GLI-regulated enhancers occurs independently of GLI3. Consistent with previous studies, some individual GLI target genes are mis-regulated in Smarcc1 conditional knockouts, though most GLI target genes are unaffected. Moreover, SMARCC1 is not necessary for mediating constitutive GLI repression in HH mutant limb buds. We conclude that SMARCC1 does not mediate GLI3 repression, which we propose utilizes alternative chromatin remodeling complexes.

EDClust: an EM-MM hybrid method for cell clustering in multiple-subject single-cell RNA sequencing.

MOTIVATION: Single-cell RNA sequencing (scRNA-seq) has revolutionized biological research by enabling the measurement of transcriptomic profiles at the single-cell level. With the increasing application of scRNA-seq in larger-scale studies, the problem of appropriately clustering cells emerges when the scRNA-seq data are from multiple subjects. One challenge is the subject-specific variation; systematic heterogeneity from multiple subjects may have a significant impact on clustering accuracy. Existing methods seeking to address such effects suffer from several limitations. RESULTS: We develop a novel statistical method, EDClust, for multi-subject scRNA-seq cell clustering. EDClust models the sequence read counts by a mixture of Dirichlet-multinomial distributions and explicitly accounts for cell-type heterogeneity, subject heterogeneity and clustering uncertainty. An EM-MM hybrid algorithm is derived for maximizing the data likelihood and clustering the cells. We perform a series of simulation studies to evaluate the proposed method and demonstrate the outstanding performance of EDClust. Comprehensive benchmarking on four real scRNA-seq datasets with various tissue types and species demonstrates the substantial accuracy improvement of EDClust compared to existing methods. AVAILABILITY AND IMPLEMENTATION: The R package is freely available at https://github.com/weix21/EDClust. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Contrasting Association of Maternal Plasma Biomarkers of Smoking and 1-Carbon Micronutrients with Offspring DNA Methylation: Evidence of Aryl Hydrocarbon Receptor Repressor Gene-Smoking-Folate Interaction.

BACKGROUND: Maternal prenatal smoking is known to alter offspring DNA methylation (DNAm). However, there are no effective interventions to mitigate smoking-induced DNAm alteration. OBJECTIVES: This study investigated whether 1-carbon nutrients (folate, vitamins B6, and B12) can protect against prenatal smoking-induced offspring DNAm alterations in the aryl hydrocarbon receptor repressor (AHRR) (cg05575921), GFI1 (cg09935388), and CYP1A1 (cg05549655) genes. METHODS: This study included mother-newborn dyads from a racially diverse US birth cohort. The cord blood DNAm at the above 3 sites were derived from a previous study using the Illumina Infinium MethylationEPIC BeadChip. Maternal smoking was assessed by self-report and plasma biomarkers (hydroxycotinine and cotinine). Maternal plasma folate, and vitamins B6 and B12 concentrations were obtained shortly after delivery. Linear regressions, Bayesian kernel machine regression, and quantile g-computation were applied to test the study hypothesis by adjusting for covariables and multiple testing. RESULTS: The study included 834 mother-newborn dyads (16.7% of newborns exposed to maternal smoking). DNAm at cg05575921 (AHRR) and at cg09935388 (GFI1) was inversely associated with maternal smoking biomarkers in a dose-response fashion (all P < 7.01 × 10-13). In contrast, cg05549655 (CYP1A1) was positively associated with maternal smoking biomarkers (P < 2.4 × 10-6). Folate concentrations only affected DNAm levels at cg05575921 (AHRR, P = 0.014). Regression analyses showed that compared with offspring with low hydroxycotinine exposure (<0.494) and adequate maternal folate concentrations (quartiles 2-4), an offspring with high hydroxycotinine exposure (≥0.494) and low folate concentrations (quartile 1) had a significant reduction in DNAm at cg05575921 (M-value, ß ± SE = -0.801 ± 0.117, P = 1.44 × 10-11), whereas adequate folate concentrations could cut smoking-induced hypomethylation by almost half. Exposure mixture models further supported the protective role of adequate folate concentrations against smoking-induced aryl hydrocarbon receptor repressor (AHRR) hypomethylation. CONCLUSIONS: This study found that adequate maternal folate can attenuate maternal smoking-induced offspring AHRR cg05575921 hypomethylation, which has been previously linked to a range of pediatric and adult diseases.

Gata6 potently initiates reprograming of pluripotent and differentiated cells to extraembryonic endoderm stem cells.

Transcription factor-mediated reprograming is a powerful method to study cell fate changes. In this study, we demonstrate that the transcription factor Gata6 can initiate reprograming of multiple cell types to induced extraembryonic endoderm stem (iXEN) cells. Intriguingly, Gata6 is sufficient to drive iXEN cells from mouse pluripotent cells and differentiated neural cells. Furthermore, GATA6 induction in human embryonic stem (hES) cells also down-regulates pluripotency gene expression and up-regulates extraembryonic endoderm (ExEn) genes, revealing a conserved function in mediating this cell fate switch. Profiling transcriptional changes following Gata6 induction in mES cells reveals step-wise pluripotency factor disengagement, with initial repression of Nanog and Esrrb, then Sox2, and finally Oct4, alongside step-wise activation of ExEn genes. Chromatin immunoprecipitation and subsequent high-throughput sequencing analysis shows Gata6 enrichment near pluripotency and endoderm genes, suggesting that Gata6 functions as both a direct repressor and activator. Together, this demonstrates that Gata6 is a versatile and potent reprograming factor that can act alone to drive a cell fate switch from diverse cell types.

Regional Gene Expression in the Retina, Optic Nerve Head, and Optic Nerve of Mice with Optic Nerve Crush and Experimental Glaucoma.

A major risk factor for glaucomatous optic neuropathy is the level of intraocular pressure (IOP), which can lead to retinal ganglion cell axon injury and cell death. The optic nerve has a rostral unmyelinated portion at the optic nerve head followed by a caudal myelinated region. The unmyelinated region is differentially susceptible to IOP-induced damage in rodent models and human glaucoma. While several studies have analyzed gene expression changes in the mouse optic nerve following optic nerve injury, few were designed to consider the regional gene expression differences that exist between these distinct areas. We performed bulk RNA-sequencing on the retina and separately micro-dissected unmyelinated and myelinated optic nerve regions from naïve C57BL/6 mice, mice after optic nerve crush, and mice with microbead-induced experimental glaucoma (total = 36). Gene expression patterns in the naïve unmyelinated optic nerve showed significant enrichment of the Wnt, Hippo, PI3K-Akt, and transforming growth factor ß pathways, as well as extracellular matrix-receptor and cell membrane signaling pathways, compared to the myelinated optic nerve and retina. Gene expression changes induced by both injuries were more extensive in the myelinated optic nerve than the unmyelinated region, and greater after nerve crush than glaucoma. Changes present three and fourteen days after injury largely subsided by six weeks. Gene markers of reactive astrocytes did not consistently differ between injury states. Overall, the transcriptomic phenotype of the mouse unmyelinated optic nerve was significantly different from immediately adjacent tissues, likely dominated by expression in astrocytes, whose junctional complexes are inherently important in responding to IOP elevation.

Differential Cytokine Signatures of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and Influenza Infection Highlight Key Differences in Pathobiology.

BACKGROUND: Several inflammatory cytokines are upregulated in severe coronavirus disease 2019 (COVID-19). We compared cytokines in COVID-19 versus influenza to define differentiating features of the inflammatory response to these pathogens and their association with severe disease. Because elevated body mass index (BMI) is a known risk factor for severe COVID-19, we examined the relationship of BMI to cytokines associated with severe disease. METHODS: Thirty-seven cytokines and chemokines were measured in plasma from 135 patients with COVID-19, 57 patients with influenza, and 30 healthy controls. Controlling for BMI, age, and sex, differences in cytokines between groups were determined by linear regression and random forest prediction was used to determine the cytokines most important in distinguishing severe COVID-19 and influenza. Mediation analysis was used to identify cytokines that mediate the effect of BMI and age on disease severity. RESULTS: Interleukin-18 (IL-18), IL-1ß, IL-6, and tumor necrosis factor-α (TNF-α) were significantly increased in COVID-19 versus influenza patients, whereas granulocyte macrophage colony-stimulating factor, interferon-γ (IFN-γ), IFN-λ1, IL-10, IL-15, and monocyte chemoattractant protein 2 were significantly elevated in the influenza group. In subgroup analysis based on disease severity, IL-18, IL-6, and TNF-α were elevated in severe COVID-19, but not in severe influenza. Random forest analysis identified high IL-6 and low IFN-λ1 levels as the most distinct between severe COVID-19 and severe influenza. Finally, IL-1RA was identified as a potential mediator of the effects of BMI on COVID-19 severity. CONCLUSIONS: These findings point to activation of fundamentally different innate immune pathways in severe acute respiratory syndrome coronavirus 2 and influenza infection, and emphasize drivers of severe COVID-19 to focus both mechanistic and therapeutic investigations.

Fetal lipidome and incident risk of food allergy: A prospective birth cohort study.

BACKGROUND: Lipids are proposed to be important in developing adaptive immunity and allergy. However, studies to date reported inconsistent findings. OBJECTIVE: To examine newborn lipidome (a comprehensive profiling of circulating lipid metabolites) on child's risk of developing food allergy (FA). The maternal-cord joint effects of lipid metabolites on FA development were also investigated. METHODS: This study included 647 mother-child pairs from the Boston Birth Cohort and analyzed 202 lipid metabolites in cord plasma profiled by liquid chromatography tandem mass spectrometry. FA was defined based on standard clinical criteria. Logistic regression was applied to examine the relationships between individual metabolites and risk of FA. RESULTS: Of the 647 children, 61 developed FA. Cord triacylglycerols of long carbon chains and multiple double bonds were significantly associated with decreased risk of FA. These associations were comparable across strata of pertinent maternal and child covariates, and were independent of maternal triacylglycerols when assessed simultaneously. Besides, cord and maternal triacylglycerols had an additive effect in association with risk of FA: Children having high (≥Median) C56:8 triacylglycerol levels in both cord and maternal plasma were at the lowest risk of developing FA (OR = 0.24, 95% CI = 0.10-0.56, p = .001), compared to those having low levels in both cord and maternal plasma. CONCLUSION: This is the first birth cohort study to link altered cord plasma lipidome with future risk of development FA during childhood. It calls for further investigation on triacylglycerols of long carbon chains and multiple double bonds as potential novel predictive biomarkers and therapeutic targets for FA.

Distinct Cytokine and Chemokine Dysregulation in Hospitalized Children With Acute Coronavirus Disease 2019 and Multisystem Inflammatory Syndrome With Similar Levels of Nasopharyngeal Severe Acute Respiratory Syndrome Coronavirus 2 Shedding.

BACKGROUND: Multisystem inflammatory syndrome in children (MIS-C) is a severe clinical phenotype of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection that remains poorly understood. METHODS: Hospitalized children <18 years of age with suspected coronavirus disease 2019 (COVID-19) (N = 53) were recruited into a prospective cohort study; 32 had confirmed COVID-19, with 16 meeting the US Centers for Disease Control criteria for MIS-C. Differences in nasopharyngeal viral ribonucleic acid (RNA) levels, SARS-CoV-2 seropositivity, and cytokine/chemokine profiles were examined, including after adjustments for age and sex. RESULTS: The median ages for those with and without MIS-C were 8.7 years (interquartile range [IQR], 5.5-13.9) and 2.2 years (IQR, 1.1-10.5), respectively (P = .18), and nasopharyngeal levels of SARS-CoV-2 RNA did not differ significantly between the 2 groups (median 63 848.25 copies/mL versus 307.1 copies/mL, P = .66); 75% of those with MIS-C were antibody positive compared with 44% without (P = .026). Levels of 14 of 37 cytokines/chemokines (interleukin [IL]-1RA, IL-2RA, IL-6, IL-8, tumor necrosis factor-α, IL-10, IL-15, IL-18, monocyte chemoattractant protein [MCP]-1, IP-10, macrophage-inflammatory protein [MIP]-1α, MCP-2, MIP-1ß, eotaxin) were significantly higher in children with MIS-C compared to those without, irrespective of age or sex (false discovery rate <0.05; P < .05). CONCLUSIONS: The distinct pattern of heightened cytokine/chemokine dysregulation observed with MIS-C, compared with acute COVID-19, occurs across the pediatric age spectrum and with similar levels of nasopharyngeal SARS-CoV-2 RNA.

Global gene expression and chromatin accessibility of the peripheral nervous system in animal models of persistent pain.

BACKGROUND: Efforts to understand genetic variability involved in an individual's susceptibility to chronic pain support a role for upstream regulation by epigenetic mechanisms. METHODS: To examine the transcriptomic and epigenetic basis of chronic pain that resides in the peripheral nervous system, we used RNA-seq and ATAC-seq of the rat dorsal root ganglion (DRG) to identify novel molecular pathways associated with pain hypersensitivity in two well-studied persistent pain models induced by chronic constriction injury (CCI) of the sciatic nerve and intra-plantar injection of complete Freund's adjuvant (CFA) in rats. RESULTS: Our RNA-seq studies identify a variety of biological process related to synapse organization, membrane potential, transmembrane transport, and ion binding. Interestingly, genes that encode transcriptional regulators were disproportionately downregulated in both models. Our ATAC-seq data provide a comprehensive map of chromatin accessibility changes in the DRG. A total of 1123 regions showed changes in chromatin accessibility in one or both models when compared to the naïve and 31 shared differentially accessible regions (DAR)s. Functional annotation of the DARs identified disparate molecular functions enriched for each pain model which suggests that chromatin structure may be altered differently following sciatic nerve injury and hind paw inflammation. Motif analysis identified 17 DNA sequences known to bind transcription factors in the CCI DARs and 33 in the CFA DARs. Two motifs were significantly enriched in both models. CONCLUSIONS: Our improved understanding of the changes in chromatin accessibility that occur in chronic pain states may identify regulatory genomic elements that play essential roles in modulating gene expression in the DRG.

A prospective cohort study on the intersectionality of obesity, chronic disease, social factors, and incident risk of COVID-19 in US low-income minority middle-age mothers.

OBJECTIVE: Coronavirus disease 2019 (COVID-19) has disproportionally affected communities of color. We aimed to determine what factors are associated with COVID-19 testing and test positivity in an underrepresented, understudied, and underreported (U3) population of mothers. METHODS: This study included 2996 middle-aged mothers of the Boston Birth Cohort (a sample of predominantly urban, low-income, Black and Hispanic mothers) who were enrolled shortly after they gave birth and followed onward at the Boston Medical Center. COVID-19 testing and test positivity were defined by the SARS-CoV-2 nucleic acid test. Two-probit Heckman selection models were performed to identify factors associated with test positivity while accounting for potential selection associated with COVID testing. RESULTS: The mean (SD) age of study mothers was 41.9 (±7.7) years. In the sample, 1741 (58.1%) and 667 (22.3%) mothers were self-identified as Black and Hispanic, respectively. A total of 396 mothers had COVID-19 testing and of those, 95 mothers tested positive from March 2020 to February 2021. Among a multitude of factors examined, factors associated with the probability of being tested were obesity (RR = 1.27; 95% confidence interval (CI): 1.08-1.49); and presence of preexisting chronic medical conditions including hypertension, asthma, stroke, and other comorbidities (coronary heart disease, chronic kidney disease, and sickle cell disease) with a corresponding RR = 1.40 (95% CI: 1.23-1.60); 1.29 (95% CI: 1.11-1.50); 1.44 (95% CI: 1.23-1.68); and 1.37 (95% CI: 1.12-1.67), respectively. Factors associated with higher incident risk of a positive COVID-19 test were body mass index, birthplace outside of the USA, and being without a college-level education. CONCLUSIONS: This study demonstrated the intersectionality of obesity and social factors in modulating incident risk of COVID-19 in this sample of US Black and Hispanic middle-aged mothers. Methodologically, our findings underscore the importance of accounting for potential selection bias in COVID-19 testing in order to obtain unbiased estimates of COVID-19 infection.