Large Set Of Genes Research Articles

BdERECTA controls vasculature patterning and phloem-xylem organization in Brachypodium distachyon

BackgroundThe vascular system of plants consists of two main tissue types, xylem and phloem. These tissues are organized into vascular bundles that are arranged into a complex network running through the plant that is essential for the viability of land plants. Despite their obvious importance, the genes involved in the organization of vascular tissues remain poorly understood in grasses.ResultsWe studied in detail the vascular network in stems from the model grass Brachypodium distachyon (Brachypodium) and identified a large set of genes differentially expressed in vascular bundles versus parenchyma tissues. To decipher the underlying molecular mechanisms of vascularization in grasses, we conducted a forward genetic screen for abnormal vasculature. We identified a mutation that severely affected the organization of vascular tissues. This mutant displayed defects in anastomosis of the vascular network and uncommon amphivasal vascular bundles. The causal mutation is a premature stop codon in ERECTA, a LRR receptor-like serine/threonine-protein kinase. Mutations in this gene are pleiotropic indicating that it serves multiple roles during plant development. This mutant also displayed changes in cell wall composition, gene expression and hormone homeostasis.ConclusionIn summary, ERECTA has a pleiotropic role in Brachypodium. We propose a major role of ERECTA in vasculature anastomosis and vascular tissue organization in Brachypodium.

Genes on the circular code alphabet

The X motifs, motifs from the circular code X, are enriched in the (protein coding) genes of bacteria, archaea, eukaryotes, plasmids and viruses, moreover, in the minimal gene set belonging to the three domains of life, as well as in tRNA and rRNA sequences. They allow to retrieve, maintain and synchronize the reading frame in genes, and contribute to the regulation of gene expression. These results lead here to a theoretical study of genes based on the circular code alphabet. A new occurrence relation of the circular code X under the hypothesis of an equiprobable (balanced) strand pairing is given. Surprisingly, a statistical analysis of a large set of bacterial genes retrieves this relation on the circular code alphabet, but not on the DNA alphabet. Furthermore, the circular code X has the strongest balanced circular code pairing among 216 maximal C3 self-complementary trinucleotide circular codes, a new property of this circular code X. As an application of this theory, different tRNAs studied on the circular code alphabet reveal an unexpected stem structure. Thus, the circular code X would have constructed a coding stem in tRNAs as an outline of the future gene structure and the future DNA double helix.

Formative pluripotent stem cells show features of epiblast cells poised for gastrulation

The pluripotency of mammalian early and late epiblast could be recapitulated by naïve embryonic stem cells (ESCs) and primed epiblast stem cells (EpiSCs), respectively. However, these two states of pluripotency may not be sufficient to reflect the full complexity and developmental potency of the epiblast during mammalian early development. Here we report the establishment of self-renewing formative pluripotent stem cells (fPSCs) which manifest features of epiblast cells poised for gastrulation. fPSCs can be established from different mouse ESCs, pre-/early-gastrula epiblasts and induced PSCs. Similar to pre-/early-gastrula epiblasts, fPSCs show the transcriptomic features of formative pluripotency, which are distinct from naïve ESCs and primed EpiSCs. fPSCs show the unique epigenetic states of E6.5 epiblast, including the super-bivalency of a large set of developmental genes. Just like epiblast cells immediately before gastrulation, fPSCs can efficiently differentiate into three germ layers and primordial germ cells (PGCs) in vitro. Thus, fPSCs highlight the feasibility of using PSCs to explore the development of mammalian epiblast.

The Krüppel-like factor Cabut has cell cycle regulatory properties similar to E2F1

Using a gain-of-function screen in Drosophila, we identified the Krüppel-like factor Cabut (Cbt) as a positive regulator of cell cycle gene expression and cell proliferation. Enforced cbt expression is sufficient to induce an extra cell division in the differentiating fly wing or eye, and also promotes intestinal stem cell divisions in the adult gut. Although inappropriate cell proliferation also results from forced expression of the E2f1 transcription factor or its target, Cyclin E, Cbt does not increase E2F1 or Cyclin E activity. Instead, Cbt regulates a large set of E2F1 target genes independently of E2F1, and our data suggest that Cbt acts via distinct binding sites in target gene promoters. Although Cbt was not required for cell proliferation during wing or eye development, Cbt is required for normal intestinal stem cell divisions in the midgut, which expresses E2F1 at relatively low levels. The E2F1-like functions of Cbt identify a distinct mechanism for cell cycle regulation that may be important in certain normal cell cycles, or in cells that cycle inappropriately, such as cancer cells.

H3.3-K27M drives neural stem cell-specific gliomagenesis in a human iPSC-derived model.

Diffuse intrinsic pontine glioma (DIPG) is an aggressive childhood tumor of the brainstem with currently no curative treatment available. The vast majority of DIPGs carry a histone H3 mutation leading to a lysine 27-to-methionine exchange (H3K27M). We engineered human induced pluripotent stem cells (iPSCs) to carry an inducible H3.3-K27M allele in the endogenous locus and studied the effects of the mutation in different disease-relevant neural cell types. H3.3-K27M upregulated bivalent promoter-associated developmental genes, producing diverse outcomes in different cell types. While being fatal for iPSCs, H3.3-K27M increased proliferation in neural stem cells (NSCs) and to a lesser extent in oligodendrocyte progenitor cells (OPCs). Only NSCs gave rise to tumors upon induction of H3.3-K27M and TP53 inactivation in an orthotopic xenograft model recapitulating human DIPGs. In NSCs, H3.3-K27M leads to maintained expression of stemness and proliferative genes and a premature activation of OPC programs that together may cause tumor initiation.

A multiplex CRISPR interference tool for virulence gene interrogation in Legionella pneumophila

Catalytically inactive dCas9 imposes transcriptional gene repression by sterically precluding RNA polymerase activity at a given gene to which it was directed by CRISPR (cr)RNAs. This gene silencing technology, known as CRISPR interference (CRISPRi), has been employed in various bacterial species to interrogate genes, mostly individually or in pairs. Here, we developed a multiplex CRISPRi platform in the pathogen Legionella pneumophila capable of silencing up to ten genes simultaneously. Constraints on precursor-crRNA expression were overcome by combining a strong promoter with a boxA element upstream of a CRISPR array. Using crRNAs directed against virulence protein-encoding genes, we demonstrated that CRISPRi is fully functional not only during growth in axenic media, but also during macrophage infection, and that gene depletion by CRISPRi recapitulated the growth defect of deletion strains. By altering the position of crRNA-encoding spacers within the CRISPR array, our platform achieved the gradual depletion of targets that was mirrored by the severity in phenotypes. Multiplex CRISPRi thus holds great promise for probing large sets of genes in bulk in order to decipher virulence strategies of L. pneumophila and other bacterial pathogens.

XGSEA: CROSS-species gene set enrichment analysis via domain adaptation.

Gene set enrichment analysis (GSEA) has been widely used to identify gene sets with statistically significant difference between cases and controls against a large gene set. GSEA needs both phenotype labels and expression of genes. However, gene expression are assessed more often for model organisms than minor species. Also, importantly gene expression are not measured well under specific conditions for human, due to high risk of direct experiments, such as non-approved treatment or gene knockout, and then often substituted by mouse. Thus, predicting enrichment significance (on a phenotype) of a given gene set of a species (target, say human), by using gene expression measured under the same phenotype of the other species (source, say mouse) is a vital and challenging problem, which we call CROSS-species gene set enrichment problem (XGSEP). For XGSEP, we propose the CROSS-species gene set enrichment analysis (XGSEA), with three steps of: (1) running GSEA for a source species to obtain enrichment scores and $p$-values of source gene sets; (2) representing the relation between source and target gene sets by domain adaptation; and (3) using regression to predict $p$-values of target gene sets, based on the representation in (2). We extensively validated the XGSEA by using five regression and one classification measurements on four real data sets under various settings, proving that the XGSEA significantly outperformed three baseline methods in most cases. A case study of identifying important human pathways for T -cell dysfunction and reprogramming from mouse ATAC-Seq data further confirmed the reliability of the XGSEA. Source code of the XGSEA is available through https://github.com/LiminLi-xjtu/XGSEA.

TEAD1 protects against necroptosis in postmitotic cardiomyocytes through regulation of nuclear DNA-encoded mitochondrial genes.

The Hippo signaling effector, TEAD1 plays an essential role in cardiovascular development. However, a role for TEAD1 in postmitotic cardiomyocytes (CMs) remains incompletely understood. Herein we reported that TEAD1 is required for postmitotic CM survival. We found that adult mice with ubiquitous or CM-specific loss of Tead1 present with a rapid lethality due to an acute-onset dilated cardiomyopathy. Surprisingly, deletion of Tead1 activated the necroptotic pathway and induced massive cardiomyocyte necroptosis, but not apoptosis. In contrast to apoptosis, necroptosis is a pro-inflammatory form of cell death and consistent with this, dramatically higher levels of markers of activated macrophages and pro-inflammatory cytokines were observed in the hearts of Tead1 knockout mice. Blocking necroptosis by administration of necrostatin-1 rescued Tead1 deletion-induced heart failure. Mechanistically, genome-wide transcriptome and ChIP-seq analysis revealed that in adult hearts, Tead1 directly activates a large set of nuclear DNA-encoded mitochondrial genes required for assembly of the electron transfer complex and the production of ATP. Loss of Tead1 expression in adult CMs increased mitochondrial reactive oxygen species, disrupted the structure of mitochondria, reduced complex I-IV driven oxygen consumption and ATP levels, resulting in the activation of necroptosis. This study identifies an unexpected paradigm in which TEAD1 is essential for postmitotic CM survival by maintaining the expression of nuclear DNA-encoded mitochondrial genes required for ATP synthesis.

Pairing of segmentation clock genes drives robust pattern formation.

Gene expression is an inherently stochastic process 1,2; however, organismal development and homeostasis require that cells spatiotemporally coordinate the expression of large sets of genes. Coexpressed gene pairs in metazoans often reside in the same chromosomal neighborhood, with gene pairs representing 10% - 50% of all genes depending on species 3–6. As shared upstream regulators can ensure correlated gene expression, the selective advantage of maintaining adjacent gene pairs remains unknown 6. Here, using two linked zebrafish segmentation clock genes, her1 and her7, and combining single-cell transcript counting, genetic engineering, real-time imaging and computational modeling, we reveal that gene pairing boosts correlated transcription and provides phenotypic robustness for developmental pattern formation. Our results demonstrate that disrupting gene pairing disrupts oscillations and segmentation, identifying the selective pressure retaining correlated transcription to sustain a robust and rapid developmental clock. We anticipate that these findings will inspire investigating advantages of gene pairing in other systems and engineering precise synthetic clocks in embryos and organoids.

Fate of antibiotic resistant E. coli and antibiotic resistance genes during full scale conventional and advanced anaerobic digestion of sewage sludge.

Antibiotic resistant bacteria (ARB) and their genes (ARGs) have become recognised as significant emerging environmental pollutants. ARB and ARGs in sewage sludge can be transmitted back to humans via the food chain when sludge is recycled to agricultural land, making sludge treatment key to control the release of ARB and ARGs to the environment. This study investigated the fate of antibiotic resistant Escherichia coli and a large set of antibiotic resistance genes (ARGs) during full scale anaerobic digestion (AD) of sewage sludge at two U.K. wastewater treatment plants and evaluated the impact of thermal hydrolysis (TH) pre-treatment on their abundance and diversity. Absolute abundance of 13 ARGs and the Class I integron gene intI1 was calculated using single gene quantitative (q) PCR. High through-put qPCR analysis was also used to determine the relative abundance of 370 ARGs and mobile genetic elements (MGEs). Results revealed that TH reduced the absolute abundance of all ARGs tested and intI1 by 10-12,000 fold. After subsequent AD, a rebound effect was seen in many ARGs. The fate of ARGs during AD without pre-treatment was variable. Relative abundance of most ARGs and MGEs decreased or fluctuated, with the exception of macrolide resistance genes, which were enriched at both plants, and tetracyline and glycopeptide resistance genes which were enriched in the plant employing TH. Diversity of ARGs and MGEs decreased in both plants during sludge treatment. Principal coordinates analysis revealed that ARGs are clearly distinguished according to treatment step, whereas MGEs in digested sludge cluster according to site. This study provides a comprehensive within-digestor analysis of the fate of ARGs, MGEs and antibiotic resistant E. coli and highlights the effectiveness of AD, particularly when TH is used as a pre-treatment, at reducing the abundance of most ARGs and MGEs in sludgeand preventing their release into the environment.

Genes Contributing to the Unique Biology and Intrinsic Antibiotic Resistance of Enterococcus faecalis.

The enterococci, which are among the leading causes of multidrug-resistant (MDR) hospital infection, are notable for their environmental ruggedness, which extends to intrinsic antibiotic resistance. To identify genes that confer this unique property, we used Tn-seq to comprehensively explore the genome of MDR Enterococcus faecalis strain MMH594 for genes important for growth in nutrient-containing medium and with low-level antibiotic challenge. As expected, a large core of genes for DNA replication, expression, and central metabolism, shared with other bacteria, are intolerant to transposon disruption. However, genes were identified that are important to E. faecalis that are either absent from or unimportant for Staphylococcus aureus and Streptococcus pneumoniae fitness when similarly tested. Further, 217 genes were identified that when challenged by sub-MIC antibiotic levels exhibited reduced tolerance to transposon disruption, including those previously shown to contribute to intrinsic resistance, and others not previously ascribed this role. E. faecalis is one of the few Gram-positive bacteria experimentally shown to possess a functional Entner-Doudoroff pathway for carbon metabolism, a pathway that contributes to stress tolerance in other microbes. Through functional genomics and network analysis we defined the unusual structure of this pathway in E. faecalis and assessed its importance. These approaches also identified toxin-antitoxin and related systems that are unique and active in E. faecalis Finally, we identified genes that are absent in the closest nonenterococcal relatives, the vagococci, and that contribute importantly to fitness with and without antibiotic selection, advancing an understanding of the unique biology of enterococci.IMPORTANCE Enterococci are leading causes of antibiotic-resistant infection transmitted in hospitals. The intrinsic hardiness of these organisms allows them to survive disinfection practices and then proliferate in the gastrointestinal tracts of antibiotic-treated patients. The objective of this study was to identify the underlying genetic basis for its unusual hardiness. Using a functional genomic approach, we identified traits and pathways of general importance for enterococcal survival and growth that distinguish them from closely related pathogens as well as ancestrally related species. We further identified unique traits that enable them to survive antibiotic challenge, revealing a large set of genes that contribute to intrinsic antibiotic resistance and a smaller set of uniquely important genes that are rare outside enterococci.

The molecular mechanisms that determine different degrees of polyphagy in the Bemisia tabaci species complex.

The whitefly Bemisia tabaci is a closely related group of >35 cryptic species that feed on the phloem sap of a broad range of host plants. Species in the complex differ in their host‐range breadth, but the mechanisms involved remain poorly understood. We investigated, therefore, how six different B. tabaci species cope with the environmental unpredictability presented by a set of four common and novel host plants. Behavioral studies indicated large differences in performances on the four hosts and putative specialization of one of the species to cassava plants. Transcriptomic analyses revealed two main insights. First, a large set of genes involved in metabolism (>85%) showed differences in expression between the six species, and each species could be characterized by its own unique expression pattern of metabolic genes. However, within species, these genes were constitutively expressed, with a low level of environmental responsiveness (i.e., to host change). Second, within each species, sets of genes mainly associated with the super‐pathways “environmental information processing” and “organismal systems” responded to the host switching events. These included genes encoding for proteins involved in sugar homeostasis, signal transduction, membrane transport, and immune, endocrine, sensory and digestive responses. Our findings suggested that the six B. tabaci species can be divided into four performance/transcriptomic “Types” and that polyphagy can be achieved in multiple ways. However, polyphagy level is determined by the specific identity of the metabolic genes/pathways that are enriched and overexpressed in each species (the species' individual metabolic “tool kit”).

PyGNA: a unified framework for geneset network analysis

BackgroundGene and protein interaction experiments provide unique opportunities to study the molecular wiring of a cell. Integrating high-throughput functional genomics data with this information can help identifying networks associated with complex diseases and phenotypes.ResultsHere we introduce an integrated statistical framework to test network properties of single and multiple genesets under different interaction models. We implemented this framework as an open-source software, called Python Geneset Network Analysis (PyGNA). Our software is designed for easy integration into existing analysis pipelines and to generate high quality figures and reports. We also developed PyGNA to take advantage of multi-core systems to generate calibrated null distributions on large datasets. We then present the results of extensive benchmarking of the tests implemented in PyGNA and a use case inspired by RNA sequencing data analysis, showing how PyGNA can be easily integrated to study biological networks. PyGNA is available at http://github.com/stracquadaniolab/pygna and can be easily installed using the PyPi or Anaconda package managers, and Docker.ConclusionsWe present a tool for network-aware geneset analysis. PyGNA can either be readily used and easily integrated into existing high-performance data analysis pipelines or as a Python package to implement new tests and analyses. With the increasing availability of population-scale omic data, PyGNA provides a viable approach for large scale geneset network analysis.

HDAC3 Mediates the Inflammatory Response and LPS Tolerance in Human Monocytes and Macrophages.

Histone deacetylases (HDACs) are a group of enzymes that control histone deacetylation and bear potential to direct expression of large gene sets. We determined the effect of HDAC inhibitors (HDACi) on human monocytes and macrophages, with respect to their polarization, activation, and their capabilities of inducing endotoxin tolerance. To address the role for HDACs in macrophage polarization, we treated monocytes with HDAC3i, HDAC6i or pan-HDACi prior to polarization into M1 or M2 macrophages using IFNγ or IL-4 respectively. To study the HDAC inhibition effect on cytokine expression, macrophages were treated with HDACi prior to LPS-stimulation. TNFα, IL-6, and p40 were measured with ELISA, whereas modifications of Histone 3 and STAT1 were assessed using western blot. To address the role for HDAC3 in repeated LPS challenge induction, HDAC3i or HDAC3 siRNA was added to monocytes prior to incubation with IFNγ, which were then repeatedly challenged with LPS and analyzed by means of protein analyses and transcriptional profiling. Pan-HDACi and HDAC3i reduced cytokine secretion in monocytes and M1 macrophages, whereas HDAC6i yielded no such effect. Notably, neither pan-HDACi nor HDAC3i reduced cytokine secretion in M2 macrophages. In contrast to previous reports in mouse macrophages, HDAC3i did not affect macrophage polarization in human cells. Likewise, HDAC3 was not required for IFNγ signaling or IFNβ secretion. Cytokine and gene expression analyses confirmed that IFNγ-treated macrophages consistently develop a cytokine response after LPS repeated challenge, but pretreatment with HDAC3i or HDAC3 siRNA reinstates a state of tolerance reflected by general suppression of tolerizable genes, possibly through decreasing TLRs expression, and particularly TLR4/CD14. The development of endotoxin tolerance in macrophages is important to reduce exacerbated immune response and limit tissue damage. We conclude that HDAC3 is an attractive protein target to mediate macrophage reactivity and tolerance induction in inflammatory macrophages.

When the Search for Stemness Genes Meets the Skin Substitute Bioengineering Field: KLF4 Transcription Factor under the Light.

The transcription factor “Kruppel-like factor 4” (KLF4) is a central player in the field of pluripotent stem cell biology. In particular, it was put under the spotlight as one of the four factors of the cocktail originally described for reprogramming into induced pluripotent stem cells (iPSCs). In contrast, its possible functions in native tissue stem cells remain largely unexplored. We recently published that KLF4 is a regulator of “stemness” in human keratinocytes. We show that reducing the level of expression of this transcription factor by RNA interference or pharmacological repression promotes the ex vivo amplification and regenerative capacity of two types of cells of interest for cutaneous cell therapy: native keratinocyte stem and progenitor cells from adult epidermis, which have been used for more than three decades in skin graft bioengineering, and keratinocytes generated by the lineage-oriented differentiation of embryonic stem cells (ESCs), which have potential for the development of skin bio-bandages. At the mechanistic level, KLF4 repression alters the expression of a large set of genes involved in TGF-β1 and WNT signaling pathways. Major regulators of TGF-β bioavailability and different TGF-β receptors were targeted, notably modulating the ALK1/Smad1/5/9 axis. At a functional level, KLF4 repression produced an antagonist effect on TGF-β1-induced keratinocyte differentiation.

Cross-Sectional Transcriptional Analysis of the Aging Murine Heart.

Cardiovascular disease accounts for millions of deaths each year and is currently the leading cause of mortality worldwide. The aging process is clearly linked to cardiovascular disease, however, the exact relationship between aging and heart function is not fully understood. Furthermore, a holistic view of cardiac aging, linking features of early life development to changes observed in old age, has not been synthesized. Here, we re-purpose RNA-sequencing data previously-collected by our group, investigating gene expression differences between wild-type mice of different age groups that represent key developmental milestones in the murine lifespan. DESeq2's generalized linear model was applied with two hypothesis testing approaches to identify differentially-expressed (DE) genes, both between pairs of age groups and across mice of all ages. Pairwise comparisons identified genes associated with specific age transitions, while comparisons across all age groups identified a large set of genes associated with the aging process more broadly. An unsupervised machine learning approach was then applied to extract common expression patterns from this set of age-associated genes. Sets of genes with both linear and non-linear expression trajectories were identified, suggesting that aging not only involves the activation of gene expression programs unique to different age groups, but also the re-activation of gene expression programs from earlier ages. Overall, we present a comprehensive transcriptomic analysis of cardiac gene expression patterns across the entirety of the murine lifespan.

Association Mapping and Development of Marker-Assisted Selection Tools for the Resistance to White Pine Blister Rust in the Alberta Limber Pine Populations.

Since its introduction to North America in the early 1900s, white pine blister rust (WPBR) caused by the fungal pathogen Cronartium ribicola has resulted in substantial economic losses and ecological damage to native North American five-needle pine species. The high susceptibility and mortality of these species, including limber pine (Pinus flexilis), creates an urgent need for the development and deployment of resistant germplasm to support recovery of impacted populations. Extensive screening for genetic resistance to WPBR has been underway for decades in some species but has only started recently in limber pine using seed families collected from wild parental trees in the USA and Canada. This study was conducted to characterize Alberta limber pine seed families for WPBR resistance and to develop reliable molecular tools for marker-assisted selection (MAS). Open-pollinated seed families were evaluated for host reaction following controlled infection using C. ribicola basidiospores. Phenotypic segregation for presence/absence of stem symptoms was observed in four seed families. The segregation ratios of these families were consistent with expression of major gene resistance (MGR) controlled by a dominant R locus. Based on linkage disequilibrium (LD)-based association mapping used to detect single nucleotide polymorphism (SNP) markers associated with MGR against C. ribicola, MGR in these seed families appears to be controlled by Cr4 or other R genes in very close proximity to Cr4. These associated SNPs were located in genes involved in multiple molecular mechanisms potentially underlying limber pine MGR to C. ribicola, including NBS-LRR genes for recognition of C. ribicola effectors, signaling components, and a large set of defense-responsive genes with potential functions in plant effector-triggered immunity (ETI). Interactions of associated loci were identified for MGR selection in trees with complex genetic backgrounds. SNPs with tight Cr4-linkage were further converted to TaqMan assays to confirm their effectiveness as MAS tools. This work demonstrates the successful translation and deployment of molecular genetic knowledge into specific MAS tools that can be easily applied in a selection or breeding program to efficiently screen MGR against WPBR in Alberta limber pine populations.

Diet Modifies Pioglitazone's Influence on Hepatic PPARγ-Regulated Mitochondrial Gene Expression.

Pioglitazone (Pio) is a thiazolidinedione (TZD) insulin-sensitizing drug whose effects result predominantly from its modulation of the transcriptional activity of peroxisome proliferator-activated-receptor-gamma (PPARγ). Pio is used to treat human insulin-resistant diabetes and also frequently considered for treatment of nonalcoholic steatohepatitis (NASH). In both settings, Pio's beneficial effects are believed to result primarily from its actions on adipose PPARγ activity, which improves insulin sensitivity and reduces the delivery of fatty acids to the liver. Nevertheless, a recent clinical trial showed variable efficacy of Pio in human NASH. Hepatocytes also express PPARγ, and such expression increases with insulin resistance and in nonalcoholic fatty liver disease (NAFLD). Furthermore, mice that overexpress hepatocellular PPARγ and Pio-treated mice with extrahepatic PPARγ gene disruption develop features of NAFLD. Thus, Pio's direct impact on hepatocellular gene expression might also be a determinant of this drug's ultimate influence on insulin resistance and NAFLD. Previous studies have characterized Pio's PPARγ-dependent effects on hepatic expression of specific adipogenic, lipogenic, and other metabolic genes. However, such transcriptional regulation has not been comprehensively assessed. The studies reported here address that consideration by genome-wide comparisons of Pio's hepatic transcriptional effects in wildtype (WT) and liver-specific PPARγ-knockout (KO) mice given either control or high-fat (HFD) diets. The results identify a large set of hepatic genes for which Pio's liver PPARγ-dependent transcriptional effects are concordant with its effects on RXR-DNA binding in WT mice. These data also show that HFD modifies Pio's influence on a subset of such transcriptional regulation. Finally, our findings reveal a broader influence of Pio on PPARγ-dependent hepatic expression of nuclear genes encoding mitochondrial proteins than previously recognized. Taken together, these studies provide new insights about the tissue-specific mechanisms by which Pio affects hepatic gene expression and the broad scope of this drug's influence on such regulation.

Dissecting the epigenomic dynamics of human fetal germ cell development at single-cell resolution

Proper development of fetal germ cells (FGCs) is vital for the precise transmission of genetic and epigenetic information through generations. The transcriptional landscapes of human FGC development have been revealed; however, the epigenetic reprogramming process of FGCs remains elusive. Here, we profiled the genome-wide DNA methylation and chromatin accessibility of human FGCs at different phases as well as gonadal niche cells at single-cell resolution. First, we found that DNA methylation levels of FGCs changed in a temporal manner, whereas FGCs at different phases in the same embryo exhibited comparable DNA methylation levels and patterns. Second, we revealed the phase-specific chromatin accessibility signatures at the promoter regions of a large set of critical transcription factors and signaling pathway genes. We also identified potential distal regulatory elements including enhancers in FGCs. Third, compared with other hominid-specific retrotransposons, SVA_D might have a broad spectrum of binding capacity for transcription factors, including SOX15 and SOX17. Finally, using an in vitro culture system of human FGCs, we showed that the BMP signaling pathway promoted the cell proliferation of FGCs, and regulated the WNT signaling pathway by orchestrating the chromatin accessibility of its ligand genes. Our single-cell epigenomic atlas and functional assays provide valuable insights for understanding the strongly heterogeneous, unsynchronized, yet highly robust nature of human germ cell development.

STRESS-FORCES EMBRYONIC STEM CELLS TO INCREASE EXPRESSION OF HEPATITIS A AND HERPES SIMPLEX 1 VIRUS RECEPTORS AND TWO GENES NECESSARY FOR COVID-19 UPTAKE

Determine whether various stresses change the expression in embryonic stem cells (ESC) and their differentiated lineages, of a large set of viral receptors and susceptibility genes. Experimental. ESC were tested by bulk or single cell (sc) RNAseq after 72 hr exposures to 0-300 mM hyperosmotic sorbitol (with stemness-maintaining Leukemia Inhibitory Factor, LIF) to quantify stress-forced differentiation. Controls for normal stemness were LIF+ and normal differentiation were LIF-. RNA was isolated by RNAeasy lysis or 10XGenomics Dropseq, RNA quality was checked by Agilent TapeStation, cDNA was synthesized using Lexogen’s QuantSeq library kit, and barcoded, multiplexed and sequenced by Illumina NovaSeq 6000. Data were demultiplexed by CASAVA software and FC expression was compared between conditions. In triplicate experiments, significant fold change (FC) genes (FC≥2; FDR≤0.05, P<0.05) identified affected pathways. Proliferation or death were assayed by Hoechst staining or Trypan blue staining, respectively. Validating studies including qPCR. Using bulk RNAseq it was shown that VIM and TMPRSS2, two genes necessary for Covid-19 cell uptake are up-regulated significantly as measured by P-value and FDR. VIM is up-regulated 5.6 fold higher at 300 mM sorbitol compared with normal differentiation and TMPRSS2 is also up-regulated 3.1 fold higher by high sorbitol despite LIF presence. Havcr1 the hepatitis A virus receptor was up-regulated 55 fold over normal ESC differentiation with significance by p-Value and FDR. Pvrl1 the Herpes simplex 1 (HSV1) virus receptor was up-regulated 4.5 fold compared with ESC undergoing normal differentiation by LIF removal. These are proof-of-principle studies that suggest that some viruses have high susceptibility when ESC and their progeny differentiated lineages undergo hyperosmotic stress. HSV1 is known to cross the placenta, and recent data suggest that Covid-19 may cross the placenta and induce IgM detected in cord blood. We intend to extend these tests using Phthalates and cortisol and to test human embryonic stem cells for similar stress-enhanced susceptibility.