No abstract

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.

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.

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.

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. Electronic supplementary material The online version of this article (10.1186/s12915-019-0677-x) contains supplementary material, which is available to authorized users.

The SIB Swiss Institute of Bioinformatics (www.isb-sib.ch) provides world-class bioinformatics databases, software tools, services and training to the international life science community in academia and industry. These solutions allow life scientists to turn the exponentially growing amount of data into knowledge. Here, we provide an overview of SIB's resources and competence areas, with a strong focus on curated databases and SIB's most popular and widely used resources. In particular, SIB's Bioinformatics resource portal ExPASy features over 150 resources, including UniProtKB/Swiss-Prot, ENZYME, PROSITE, neXtProt, STRING, UniCarbKB, SugarBindDB, SwissRegulon, EPD, arrayMap, Bgee, SWISS-MODEL Repository, OMA, OrthoDB and other databases, which are briefly described in this article.

Albinism and leucism are phenotypes resulting from impaired melanin pigmentation in the skin and skin appendages. However, melanin pigmentation of eyes remains unaffected in leucism. Here, using transmission electron microscopy, we show that the leucistic morph of the Texas rat snake (Pantherophis obsoletus lindheimeri) lacks both melanophores and xanthophores in its skin and exhibits a uniform ivory white color generated by iridophores and collagen fibers. In addition, we sequenced the full genome of a leucistic individual and obtained a highly-contiguous near-chromosome quality assembly of 1.69 Gb with an N50 of 14.5 Mb and an L50 of 29 sequences. Using a candidate-gene approach, we then identify in the leucistic genome a single-nucleotide deletion that generates a frameshift and a premature termination codon in the melanocyte inducing transcription factor (MITF) gene. This mutation shortens the translated protein from 574 to 286 amino acids, removing the helix-loop-helix DNA-binding domain that is highly conserved among vertebrates. Genotyping leucistic animals of independent lineages showed that not all leucistic individuals carry this single-nucleotide deletion. Subsequent gene expression analyses reveal that all leucistic individuals that we analyzed exhibit a significantly decreased expression of MITF. We thus suggest that mutations affecting the regulation and, in some cases, the coding sequence of MITF, the former probably predating the latter, could be associated with the leucistic phenotype in Texas rat snakes. MITF is involved in the development and survival of melanophores in vertebrates. In zebrafish, a classical model species for pigmentation that undergoes metamorphosis, larvae and adults of homozygous mitfa mutants lack melanophores, show an excess of iridophores and exhibit reduced yellow pigmentation. On the contrary, in the leucistic Texas rat snake, a non-metamorphic species, only iridophores persist. Our results suggest that fate determination of neural-crest derived melanophores and xanthophores, but not of iridophores, could require the expression of MITF during snake embryonic development.

SARS-CoV-2 infections lead to a high risk of hospitalization and mortality in diabetic patients. Why diabetic individuals are more prone to develop severe COVID-19 remains unclear. Here, we established a novel human kidney organoid model that mimics early hallmarks of diabetic nephropathy. High oscillatory glucose exposure resulted in metabolic changes, expansion of extracellular membrane components, gene expression changes determined by scRNAseq, and marked upregulation of angiotensin-converting enzyme 2 (ACE2). Upon SARS-CoV-2 infection, hyperglycemic conditions lead to markedly higher viral loads in kidney organoids compared to normoglycemia. Genetic deletion of ACE2, but not of the candidate receptor BSG/CD147, in kidney organoids demonstrated the essential role of ACE2 in SARS-CoV-2 infections and completely prevented SARS-CoV-2 infection in the diabetogenic microenvironment. These data introduce a novel organoid model for diabetic kidney disease and show that diabetic-induced ACE2 licenses the diabetic kidney to enhanced SARS-CoV-2 replication.