The dynamic behavior of chromatophores marks the transition from bands to spots in leopard geckos.

Reptilian skin coloration is spectacular and diverse, yet little is known about the ontogenetic processes that govern its establishment and the molecular signaling pathways that determine it. Here, we focus on the development of the banded pattern of leopard gecko hatchlings and the transition to black spots in the adult. With our histological analyses, we show that iridophores are present in the white and yellow bands of the hatchling and they gradually perish in the adult skin. Furthermore, we demonstrate that melanophores can autonomously form spots in the absence of the other chromatophores both on the regenerated skin of the tail and on the dorsal skin of the Mack Super Snow (MSS) leopard geckos. This color morph is characterized by uniform black coloration in hatchlings and black spots in adulthood; we establish that their skin is devoid of xanthophores and iridophores at both stages. Our genetic analyses identified a 13-nucleotide deletion in the PAX7 transcription factor of MSS geckos, affecting its protein coding sequence. With our single-cell transcriptomics analysis of embryonic skin, we confirm that PAX7 is expressed in iridophores and xanthophores, suggesting that it plays a key role in the differentiation of both chromatophores. Our in situ hybridizations on whole-mount embryos document the dynamics of the skin pattern formation and how it is impacted in the PAX7 mutants. We hypothesize that the melanophores-iridophores interactions give rise to the banded pattern of the hatchlings and black spot formation is an intrinsic capacity of melanophores in the postembryonic skin.

Single-cell transcriptional profile of CD34+ hematopoietic progenitor cells from del(5q) myelodysplastic syndromes and impact of lenalidomide.

While myelodysplastic syndromes with del(5q) (del(5q) MDS) comprises a well-defined hematological subgroup, the molecular basis underlying its origin remains unknown. Using single cell RNA-seq (scRNA-seq) on CD34+ progenitors from del(5q) MDS patients, we have identified cells harboring the deletion, characterizing the transcriptional impact of this genetic insult on disease pathogenesis and treatment response. Interestingly, both del(5q) and non-del(5q) cells present similar transcriptional lesions, indicating that all cells, and not only those harboring the deletion, may contribute to aberrant hematopoietic differentiation. However, gene regulatory network (GRN) analyses reveal a group of regulons showing aberrant activity that could trigger altered hematopoiesis exclusively in del(5q) cells, pointing to a more prominent role of these cells in disease phenotype. In del(5q) MDS patients achieving hematological response upon lenalidomide treatment, the drug reverts several transcriptional alterations in both del(5q) and non-del(5q) cells, but other lesions remain, which may be responsible for potential future relapses. Moreover, lack of hematological response is associated with the inability of lenalidomide to reverse transcriptional alterations. Collectively, this study reveals transcriptional alterations that could contribute to the pathogenesis and treatment response of del(5q) MDS.

Natural Hydrogels Support Kidney Organoid Generation and Promote In Vitro Angiogenesis.

To date, strategies aiming to modulate cell to extracellular matrix (ECM) interactions during organoid derivation remain largely unexplored. Here renal decellularized ECM (dECM) hydrogels are fabricated from porcine and human renal cortex as biomaterials to enrich cell-to-ECM crosstalk during the onset of kidney organoid differentiation from human pluripotent stem cells (hPSCs). Renal dECM-derived hydrogels are used in combination with hPSC-derived renal progenitor cells to define new approaches for 2D and 3D kidney organoid differentiation, demonstrating that in the presence of these biomaterials the resulting kidney organoids exhibit renal differentiation features and the formation of an endogenous vascular component. Based on these observations, a new method to produce kidney organoids with vascular-like structures is achieved through the assembly of hPSC-derived endothelial-like organoids with kidney organoids in 3D. Major readouts of kidney differentiation and renal cell morphology are assessed exploiting these culture platforms as new models of nephrogenesis. Overall, this work shows that exploiting cell-to-ECM interactions during the onset of kidney differentiation from hPSCs facilitates and optimizes current approaches for kidney organoid derivation thereby increasing the utility of these unique cell culture platforms for personalized medicine.

Generation of heart and vascular system in rodents by blastocyst complementation.

Generating organs from stem cells through blastocyst complementation is a promising approach to meet the clinical need for transplants. In order to generate rejection-free organs, complementation of both parenchymal and vascular cells must be achieved, as endothelial cells play a key role in graft rejection. Here, we used a lineage-specific cell ablation system to produce mouse embryos unable to form both the cardiac and vascular systems. By mouse intraspecies blastocyst complementation, we rescued heart and vascular system development separately and in combination, obtaining complemented hearts with cardiomyocytes and endothelial cells of exogenous origin. Complemented chimeras were viable and reached adult stage, showing normal cardiac function and no signs of histopathological defects in the heart. Furthermore, we implemented the cell ablation system for rat-to-mouse blastocyst complementation, obtaining xenogeneic hearts whose cardiomyocytes were completely of rat origin. These results represent an advance in the experimentation towards the in vivo generation of transplantable organs.

Somitic positional information guides self-organized patterning of snake scales.

Two influential concepts in tissue patterning are Wolpert's positional information and Turing's self-organized reaction-diffusion (RD). The latter establishes the patterning of hair and feathers. Here, our morphological, genetic, and functional-by CRISPR-Cas9-mediated gene disruption-characterization of wild-type versus "scaleless" snakes reveals that the near-perfect hexagonal pattern of snake scales is established through interactions between RD in the skin and somitic positional information. First, we show that ventral scale development is guided by hypaxial somites and, second, that ventral scales and epaxial somites guide the sequential RD patterning of the dorsolateral scales. The RD intrinsic length scale evolved to match somite periodicity, ensuring the alignment of ribs and scales, both of which play a critical role in snake locomotion.

Loss of the matrix metalloproteinase-10 causes premature features of aging in satellite cells.

Aged muscles accumulate satellite cells with a striking decline response to damage. Although intrinsic defects in satellite cells themselves are the major contributors to aging-associated stem cell dysfunction, increasing evidence suggests that changes in the muscle-stem cell local microenvironment also contribute to aging. Here, we demonstrate that loss of the matrix metalloproteinase-10 (MMP-10) in young mice alters the composition of the muscle extracellular matrix (ECM), and specifically disrupts the extracellular matrix of the satellite cell niche. This situation causes premature features of aging in the satellite cells, contributing to their functional decline and a predisposition to enter senescence under proliferative pressure. Similarly, reduction of MMP-10 levels in young satellite cells from wild type animals induces a senescence response, while addition of the protease delays this program. Significantly, the effect of MMP-10 on satellite cell aging can be extended to another context of muscle wasting, muscular dystrophy. Systemic treatment of mdx dystrophic mice with MMP-10 prevents the muscle deterioration phenotype and reduces cellular damage in the satellite cells, which are normally under replicative pressure. Most importantly, MMP-10 conserves its protective effect in the satellite cell-derived myoblasts isolated from a Duchenne muscular dystrophy patient by decreasing the accumulation of damaged DNA. Hence, MMP-10 provides a previously unrecognized therapeutic opportunity to delay satellite cell aging and overcome satellite cell dysfunction in dystrophic muscles.

Gene therapy restores the transcriptional program of hematopoietic stem cells in Fanconi anemia.

Clinical trials have shown that lentiviral-mediated gene therapy can ameliorate bone marrow failure (BMF) in nonconditioned Fanconi anemia (FA) patients resulting from the proliferative advantage of corrected FA hematopoietic stem and progenitor cells (HSPC). However, it is not yet known if gene therapy can revert affected molecular pathways in diseased HSPC. Single-cell RNA sequencing was performed in chimeric populations of corrected and uncorrected HSPC co-existing in the BM of gene therapy-treated FA patients. Our study demonstrates that gene therapy reverts the transcriptional signature of FA HSPC, which then resemble the transcriptional program of healthy donor HSPC. This includes a down-regulated expression of TGF-ß and p21, typically up-regulated in FA HSPC, and upregulation of DNA damage response and telomere maintenance pathways. Our results show for the first time the potential of gene therapy to rescue defects in the HSPC transcriptional program from patients with inherited diseases; in this case, in FA characterized by BMF and cancer predisposition.

Revealing cell populations catching the early stages of human embryo development in naive pluripotent stem cell cultures.

Naive human pluripotent stem cells (hPSCs) are defined as the in vitro counterpart of the human preimplantation embryo's epiblast and are used as a model system to study developmental processes. In this study, we report the discovery and characterization of distinct cell populations coexisting with epiblast-like cells in 5iLAF naive human induced PSC (hiPSC) cultures. It is noteworthy that these populations closely resemble different cell types of the human embryo at early developmental stages. While epiblast-like cells represent the main cell population, interestingly we detect a cell population with gene and transposable element expression profile closely resembling the totipotent eight-cell (8C)-stage human embryo, and three cell populations analogous to trophectoderm cells at different stages of their maturation process: transition, early, and mature stages. Moreover, we reveal the presence of cells resembling primitive endoderm. Thus, 5iLAF naive hiPSC cultures provide an excellent opportunity to model the earliest events of human embryogenesis, from the 8C stage to the peri-implantation period.

One-Step In Vitro Generation of ETV2-Null Pig Embryos.

Each year, tens of thousands of people worldwide die of end-stage organ failure due to the limited availability of organs for use in transplantation. To meet this clinical demand, one of the last frontiers of regenerative medicine is the generation of humanized organs in pigs from pluripotent stem cells (PSCs) via blastocyst complementation. For this, organ-disabled pig models are needed. As endothelial cells (ECs) play a critical role in xenotransplantation rejection in every organ, we aimed to produce hematoendothelial-disabled pig embryos targeting the master transcription factor ETV2 via CRISPR-Cas9-mediated genome modification. In this study, we designed five different guide RNAs (gRNAs) against the DNA-binding domain of the porcine ETV2 gene, which were tested on porcine fibroblasts in vitro. Four out of five guides showed cleavage capacity and, subsequently, these four guides were microinjected individually as ribonucleoprotein complexes (RNPs) into one-cell-stage porcine embryos. Next, we combined the two gRNAs that showed the highest targeting efficiency and microinjected them at higher concentrations. Under these conditions, we significantly improved the rate of biallelic mutation. Hence, here, we describe an efficient one-step method for the generation of hematoendothelial-disabled pig embryos via CRISPR-Cas9 microinjection in zygotes. This model could be used in experimentation related to the in vivo generation of humanized organs.

A diabetic milieu increases ACE2 expression and cellular susceptibility to SARS-CoV-2 infections in human kidney organoids and patient cells.

It is not well understood why diabetic individuals are more prone to develop severe COVID-19. To this, we here established a human kidney organoid model promoting early hallmarks of diabetic kidney disease development. Upon SARS-CoV-2 infection, diabetic-like kidney organoids exhibited higher viral loads compared with their control counterparts. Genetic deletion of the angiotensin-converting enzyme 2 (ACE2) in kidney organoids under control or diabetic-like conditions prevented viral detection. Moreover, cells isolated from kidney biopsies from diabetic patients exhibited altered mitochondrial respiration and enhanced glycolysis, resulting in higher SARS-CoV-2 infections compared with non-diabetic cells. Conversely, the exposure of patient cells to dichloroacetate (DCA), an inhibitor of aerobic glycolysis, resulted in reduced SARS-CoV-2 infections. Our results provide insights into the identification of diabetic-induced metabolic programming in the kidney as a critical event increasing SARS-CoV-2 infection susceptibility, opening the door to the identification of new interventions in COVID-19 pathogenesis targeting energy metabolism.

Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles.

Reptiles exhibit a spectacular diversity of skin colors and patterns brought about by the interactions among three chromatophore types: black melanophores with melanin-packed melanosomes, red and yellow xanthophores with pteridine- and/or carotenoid-containing vesicles, and iridophores filled with light-reflecting platelets generating structural colors. Whereas the melanosome, the only color-producing endosome in mammals and birds, has been documented as a lysosome-related organelle, the maturation paths of xanthosomes and iridosomes are unknown. Here, we first use 10x Genomics linked-reads and optical mapping to assemble and annotate a nearly chromosome-quality genome of the corn snake Pantherophis guttatus The assembly is 1.71 Gb long, with an N50 of 16.8 Mb and L50 of 24. Second, we perform mapping-by-sequencing analyses and identify a 3.9-Mb genomic interval where the lavender variant resides. The lavender color morph in corn snakes is characterized by gray, rather than red, blotches on a pink, instead of orange, background. Third, our sequencing analyses reveal a single nucleotide polymorphism introducing a premature stop codon in the lysosomal trafficking regulator gene (LYST) that shortens the corresponding protein by 603 amino acids and removes evolutionary-conserved domains. Fourth, we use light and transmission electron microscopy comparative analyses of wild type versus lavender corn snakes and show that the color-producing endosomes of all chromatophores are substantially affected in the LYST mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

Impact of genome architecture on the functional activation and repression of Hox regulatory landscapes.

BACKGROUND: The spatial organization of the mammalian genome relies upon the formation of chromatin domains of various scales. At the level of gene regulation in cis, collections of enhancer sequences define large regulatory landscapes that usually match with the presence of topologically associating domains (TADs). These domains often contain ranges of enhancers displaying similar or related tissue specificity, suggesting that in some cases, such domains may act as coherent regulatory units, with a global on or off state. By using the HoxD gene cluster, which specifies the topology of the developing limbs via highly orchestrated regulation of gene expression, as a paradigm, we investigated how the arrangement of regulatory domains determines their activity and function. RESULTS: Proximal and distal cells in the developing limb express different levels of Hoxd genes, regulated by flanking 3' and 5' TADs, respectively. We characterized the effect of large genomic rearrangements affecting these two TADs, including their fusion into a single chromatin domain. We show that, within a single hybrid TAD, the activation of both proximal and distal limb enhancers globally occurred as when both TADs are intact. However, the activity of the 3' TAD in distal cells is generally increased in the fused TAD, when compared to wild type where it is silenced. Also, target gene activity in distal cells depends on whether or not these genes had previously responded to proximal enhancers, which determines the presence or absence of H3K27me3 marks. We also show that the polycomb repressive complex 2 is mainly recruited at the Hox gene cluster and can extend its coverage to far-cis regulatory sequences as long as confined to the neighboring TAD structure. CONCLUSIONS: We conclude that antagonistic limb proximal and distal enhancers can exert their specific effects when positioned into the same TAD and in the absence of their genuine target genes. We also conclude that removing these target genes reduced the coverage of a regulatory landscape by chromatin marks associated with silencing, which correlates with its prolonged activity in time.

A Step-by-Step Guide to Assemble a Reptilian Genome.

Multiple technologies and software are now available facilitating the de novo sequencing and assembly of any vertebrate genome. Yet the quality of most available sequenced genomes is substantially poorer than that of the golden standard in the field: the human genome. Here, we present a step-by-step protocol for the successful sequencing and assembly of a high-quality snake genome that can be applied to any other reptilian or avian species. We combine the great sequencing depth and accuracy of short reads with the use of different insert size libraries for extended scaffolding followed by optical mapping. We show that this procedure improved the corn snake scaffold N50 from 3.7 kbp to 1.4 Mbp, currently making it one of the snake genomes with the longest scaffolds. Short guidelines are also given on the extraction of long DNA molecules from reptilian blood and the necessary modifications in DNA extraction protocols. This chapter is accompanied by a website ( www.reptilomics.org/stepbystep.html ), where we provide links to the suggested software, examples of input and output files, and running parameters.

Reptilian Transcriptomes v2.0: An Extensive Resource for Sauropsida Genomics and Transcriptomics.

Despite the availability of deep-sequencing techniques, genomic and transcriptomic data remain unevenly distributed across phylogenetic groups. For example, reptiles are poorly represented in sequence databases, hindering functional evolutionary and developmental studies in these lineages substantially more diverse than mammals. In addition, different studies use different assembly and annotation protocols, inhibiting meaningful comparisons. Here, we present the "Reptilian Transcriptomes Database 2.0," which provides extensive annotation of transcriptomes and genomes from species covering the major reptilian lineages. To this end, we sequenced normalized complementary DNA libraries of multiple adult tissues and various embryonic stages of the leopard gecko and the corn snake and gathered published reptilian sequence data sets from representatives of the four extant orders of reptiles: Squamata (snakes and lizards), the tuatara, crocodiles, and turtles. The LANE runner 2.0 software was implemented to annotate all assemblies within a single integrated pipeline. We show that this approach increases the annotation completeness of the assembled transcriptomes/genomes. We then built large concatenated protein alignments of single-copy genes and inferred phylogenetic trees that support the positions of turtles and the tuatara as sister groups of Archosauria and Squamata, respectively. The Reptilian Transcriptomes Database 2.0 resource will be updated to include selected new data sets as they become available, thus making it a reference for differential expression studies, comparative genomics and transcriptomics, linkage mapping, molecular ecology, and phylogenomic analyses involving reptiles. The database is available at www.reptilian-transcriptomes.org and can be enquired using a wwwblast server installed at the University of Geneva.

The genome sequence of the corn snake (Pantherophis guttatus), a valuable resource for EvoDevo studies in squamates.

Squamates (snakes and lizards) exhibit a striking variety of phenotypes, with little known on their generative mechanisms. Studies aiming to understand the genetic basis of this wide diversity in morphology, physiology and ecology will greatly benefit from whole genome sequencing initiatives, as they provide the foundation for comparative analyses and improve our understanding of the evolution, development and diversification of traits. Here, we present the first draft genome of the corn snake Pantherophis guttatus, an oviparous snake that we promote as a particularly appropriate model species for evolutionary developmental studies in squamates. We sequenced 100-base paired-end reads from multiple individuals of a single family (parents and offspring) that produced a genome assembly of 1.53 gigabases (Gb), roughly covering 75% of the expected total genome size, and 297,768 scaffolds >1 Kb. We were able to fully retrieve 86, and partially another 106, of the 248 CEGMA core genes, indicating that a high genome completeness was achieved, even though the assembly is fragmented. Using MAKER2, we annotated 10,917 genes with high confidence (Annotation Edit Distance (AED)<1) and an additional 5,263 predicted genes matched with the species' transcriptome. Numerous colour and colour pattern morphs exist in P. guttatus, making it an ideal model to study the genetic determinism, development, and evolution of adaptive colour traits in reptiles. Using our draft genome and a Single-Nucleotide Polymorphism (SNP) calling approach, we confirmed the interval with the causative mutation for the amelanistic phenotype, a result supported by a parallel exome-based study.