Rab32 and Rab38 genes in chordate pigmentation: an evolutionary perspective

Gene expression in notochord and nuclei pulposi: a study of gene families across the chordate phylum

The transition from notochord to vertebral column is a crucial milestone in chordate evolution and in prenatal development of all vertebrates. As ossification of the vertebral bodies proceeds, involutions of residual notochord cells into the intervertebral discs form the nuclei pulposi, shock-absorbing structures that confer flexibility to the spine. Numerous studies have outlined the developmental and evolutionary relationship between notochord and nuclei pulposi. However, the knowledge of the similarities and differences in the genetic repertoires of these two structures remains limited, also because comparative studies of notochord and nuclei pulposi across chordates are complicated by the gene/genome duplication events that led to extant vertebrates. Here we show the results of a pilot study aimed at bridging the information on these two structures. We have followed in different vertebrates the evolutionary trajectory of notochord genes identified in the invertebrate chordate Ciona, and we have evaluated the extent of conservation of their expression in notochord cells. Our results have uncovered evolutionarily conserved markers of both notochord development and aging/degeneration of the nuclei pulposi.

Transcriptome Reveals Multi Pigmentation Genes Affecting Dorsoventral Pattern in Avian Body

Certain animals exhibit a special dorsoventral pattern with a lighter ventral side compared to the dorsal one and this phenomenon was preserved in the long-term evolution process. Birds also retain this trait. Recently, Inaba et al. (2019) found that ASIP (agouti signal protein) regulated interconversion between different melanocyte types leads to dorsal stripe pattern, which may partly explain the birds' dorsoventral plumage color difference. In this study, we used the embryo samples of LBM (light brown mottling) ducks (Anas platyrhynchos) with white ventral and dark dorsal body parts to investigate the mechanism of dorsoventral color variation. Firstly, melanin deposition process of duck embryos was investigated. The result indicated that E13 and E16 were the active stages of melanin synthesis. Moreover, the melanin deposition on the dorsum of LBM ducks was higher than that on the ventral side throughout. Then, RNA-seq was conducted for the dorsal and ventral skin tissues from E7 (early), E13 (middle) and E19 (late) of LBM ducks. Expression pattern analysis showed that the mRNA expression of most melanin synthesis related genes were at the highest level at E13, which was consistent with the section analysis. A correlation was found between melanogenesis pathway and dorsoventral color difference by co-expression analysis. In the DEG (differentially expressed gene) analysis, we added the dorsal skin transcriptome of embryonic white and black duck of same subspecies (Anas platyrhynchos domestica) for horizontal comparison. The results showed that 8 melanogenesis related genes (TYR, TYRP1, MLANA, RAB38, OCA2, TSPAN10, MC1R, and MSLN) were the common DEGs (Differential expressed genes) in the comparisons of body parts and breeds suggesting that the underlying molecular regulatory mechanism of dorsoventral plumage color difference may be similar to that of albino and melanic duck, which were caused by the different expression of multiple genes in melanin synthesis pathway. In addition, the molecular regulation of melanin synthesis pathway in the dorsal and ventral side of LBM ducks was analyzed. In this pathway, ASIP, MC1R, TYR, and TYRP1 have differential mRNA expression. ASIP, as an upstream gene in this pathway, was likely to play a decisive role in determining the dorsoventral plumage pattern.

Genome analysis of American minks reveals link of mutations in Ras-related protein-38 gene to Moyle brown coat phenotype

Over 35 fur colours have been described in American mink (Neovison vison), only six of which have been previously linked to specific genes. Moyle fur colour belongs to a wide group of brownish colours that are highly similar to each other, which complicates selection and breeding procedures. We performed whole genome sequencing for two American minks with Moyle (m/m) and Violet (a/a m/m /p/p) phenotypes. We identified two frame-shift mutations in the gene encoding Ras-related protein-38 (RAB38), which regulates the trafficking of tyrosinase-containing vesicles to maturing melanosomes. The results highlight the role of RAB38 in the biogenesis of melanosomes and melanin and the genetic mechanism contributing to hair colour variety and intensity. These data are also useful for tracking economically valuable fur traits in mink breeding programmes.

The Cis-Regulatory Code for Kelch-like 21/30 Specific Expression in Ciona robusta Sensory Organs

The tunicate Ciona robusta is an emerging model system to study the evolution of the nervous system. Due to their small embryos and compact genomes, tunicates, like Ciona robusta, have great potential to comprehend genetic circuitry underlying cell specific gene repertoire, among different neuronal cells. Their simple larvae possess a sensory vesicle comprising two pigmented sensory organs, the ocellus and the otolith. We focused here on Klhl21/30, a gene belonging to Kelch family, that, in Ciona robusta, starts to be expressed in pigmented cell precursors, becoming specifically maintained in the otolith precursor during embryogenesis. Evolutionary analyses demonstrated the conservation of Klhl21/30 in all the chordates. Cis-regulatory analyses and CRISPR/Cas9 mutagenesis of potential upstream factors, revealed that Klhl21/30 expression is controlled by the combined action of three transcription factors, Mitf, Dmrt, and Msx, which are downstream of FGF signaling. The central role of Mitf is consistent with its function as a fundamental regulator of vertebrate pigment cell development. Moreover, our results unraveled a new function for Dmrt and Msx as transcriptional co-activators in the context of the Ciona otolith.

Rimbp, a New Marker for the Nervous System of the Tunicate Ciona robusta

Establishment of presynaptic mechanisms by proteins that regulate neurotransmitter release in the presynaptic active zone is considered a fundamental step in animal evolution. Rab3 interacting molecule-binding proteins (Rimbps) are crucial components of the presynaptic active zone and key players in calcium homeostasis. Although Rimbp involvement in these dynamics has been described in distantly related models such as fly and human, the role of this family in most invertebrates remains obscure. To fill this gap, we defined the evolutionary history of Rimbp family in animals, from sponges to mammals. We report, for the first time, the expression of the two isoforms of the unique Rimbp family member in Ciona robusta in distinct domains of the larval nervous system. We identify intronic enhancers that are able to drive expression in different nervous system territories partially corresponding to Rimbp endogenous expression. The analysis of gene expression patterns and the identification of regulatory elements of Rimbp will positively impact our understanding of this family of genes in the context of Ciona embryogenesis.

Combined deficiency of RAB32 and RAB38 in the mouse mimics Hermansky-Pudlak syndrome and critically impairs thrombosis

The biogenesis of lysosome related organelles is defective in Hermansky-Pudlak syndrome (HPS), a disorder characterized by oculocutaneous albinism and platelet dense granule (DG) defects. The first animal model of HPS was the fawn-hooded rat, harboring a spontaneous mutation inactivating the small guanosine triphosphatase Rab38 This leads to coat color dilution associated with the absence of DGs and lung morphological defects. Another RAB38 mutant, the cht mouse, has normal DGs, which has raised controversy about the role of RAB38 in DG biogenesis. We show here that murine and human, but not rat, platelets also express the closely related RAB32. To elucidate the parts played by RAB32 and RAB38 in the biogenesis of DGs in vivo and their effects on platelet functions, we generated mice inactivated for Rab32, Rab38, and both genes. Single Rab38 inactivation mimicked cht mice, whereas single Rab32 inactivation had no effect in DGs, coat color, or lung morphology. By contrast, Rab32/38 double inactivation mimicked severe HPS, with strong coat and eye pigment dilution, some enlarged lung multilamellar bodies associated with a decrease in the number of DGs. These organelles were morphologically abnormal, decreased in number, and devoid of 5-hydroxytryptamine content. In line with the storage pool defect, platelet activation was affected, resulting in severely impaired thrombus growth and prolongation of the bleeding time. Overall, our study demonstrates the absence of impact of RAB38 or RAB32 single deficiency in platelet biogenesis and function resulting from full redundancy, and characterized a new mouse model mimicking HPS devoid of DG content.

Hyperpigmented melanistic skin lesions of smallmouth bass Micropterus dolomieu from the Chesapeake Bay watershed

Hyperpigmented melanistic skin lesions (HPMLs) of smallmouth bass Micropterus dolomieu are observed in the Potomac and Susquehanna rivers, Chesapeake Bay watershed, USA. Routine, nonlethal population surveys were conducted at 8 sites on the mainstem Susquehanna River and 9 on the Juniata River, a tributary of the Susquehanna River, between 2012 and 2018, and the prevalence of HPMLs was documented. A total of 4078 smallmouth bass were collected from the mainstem Susquehanna River and 6478 from the Juniata River. Lesions were primarily seen in bass greater than 200 mm, and prevalence in the Susquehanna River (8%) was higher (p < 0.001) than in the Juniata River (2%). As part of ongoing fish health monitoring projects, smallmouth bass were collected at additional sites, primarily tributaries of the Susquehanna (n = 758) and Potomac (n = 545) rivers between 2013 and 2018. Prevalence in the Susquehanna River (13%) was higher (p < 0.001) than the Potomac (3%). Microscopically, HPMLs were characterized by an increased number of melanocytes in the epidermis or within the dermis and epidermis. RNAseq analyses of normal and melanistic skin identified 3 unique sequences in HPMLs. Two were unidentified and the third was a viral helicase (E1). Transcript abundance in 16 normal skin samples and 16 HPMLs showed upregulation of genes associated with melanogenesis and cell proliferation in HPMLs. The E1 transcript was detected in 12 of the 16 melanistic areas but in no samples from normal skin. Further research will be necessary to identify the putative new virus and determine its role in melanocyte proliferation.

Dark-rearing uncovers novel gene expression patterns in an obligate cave-dwelling fish

Extreme environments often result in the evolution of dramatic adaptive features. The Mexican tetra, Astyanax mexicanus, includes 30 different populations of cave-dwelling forms that live in perpetual darkness. As a consequence, many populations have evolved eye loss, reduced pigmentation, and amplification of nonvisual sensory systems. Closely-related surface-dwelling morphs demonstrate typical vision, pigmentation, and sensation. Transcriptomic assessments in this system have revealed important developmental changes associated with the cave morph, however, they have not accounted for photic rearing conditions. Prior studies reared individuals under a 12:12 hr light/dark (LD) cycle. Here, we reared cavefish under constant darkness (DD) for 5+ years. From these experimental individuals, we performed mRNA sequencing and compared gene expression of surface fish reared under LD conditions to cavefish reared under DD conditions to identify photic-dependent gene expression differences. Gene Ontology enrichment analyses revealed a number of previously underappreciated cave-associated changes impacting blood physiology and olfaction. We further evaluated the position of differentially expressed genes relative to QTL positions from prior studies and found several candidate genes associated with these ecologically relevant lighting conditions. In sum, this work highlights photic conditions as a key environmental factor impacting gene expression patterns in blind cave-dwelling fish.

The evolutionary landscape of the Rab family in chordates

Intracellular traffic amongst organelles represents a key feature for eukaryotes and is orchestrated principally by members of Rab family, the largest within Ras superfamily. Given that variations in Rab repertoire have been fundamental in animal diversification, we provided the most exhaustive survey regarding the Rab toolkit of chordates. Our findings reveal the existence of 42 metazoan conserved subfamilies exhibiting a univocal intron/exon structure preserved from cnidarians to vertebrates. Since the current view does not capture the Rab complexity, we propose a new Rab family classification in three distinct monophyletic clades. The Rab complement of chordates shows a dramatic diversification due to genome duplications and independent gene duplications and losses with sharp differences amongst cephalochordates, tunicates and gnathostome vertebrates. Strikingly, the analysis of the domain architecture of this family highlighted the existence of chimeric calcium-binding Rabs, which are animal novelties characterized by a complex evolutionary history in gnathostomes and whose role in cellular metabolism is obscure. This work provides novel insights in the knowledge of Rab family: our hypothesis is that chordates represent a hotspot of Rab variability, with many events of gene gains and losses impacting intracellular traffic capabilities. Our results help to elucidate the role of Rab members in the transport amongst endomembranes and shed light on intracellular traffic routes in vertebrates. Then, since the predominant role of Rabs in the molecular communication between different cellular districts, this study paves to way to comprehend inherited or acquired human disorders provoked by dysfunctions in Rab genes.

Ran GTPase, an eukaryotic gene novelty, is involved in amphioxus mitosis

Ran (ras-related nuclear protein) is a small GTPase belonging to the RAS superfamily that is specialized in nuclear trafficking. Through different accessory proteins, Ran plays key roles in several processes including nuclear import-export, mitotic progression and spindle assembly. Consequently, Ran dysfunction has been linked to several human pathologies. This work illustrates the high degree of amino acid conservation of Ran orthologues across evolution, reflected in its conserved role in nuclear trafficking. Moreover, we studied the evolutionary scenario of the pre-metazoan genetic linkage between Ran and Stx, and we hypothesized that chromosomal proximity of these two genes across metazoans could be related to a regulatory logic or a functional linkage. We studied, for the first time, Ran expression during amphioxus development and reported its presence in the neural vesicle, mouth, gill slits and gut corresponding to body regions involved in active cell division.

Mechanisms of developmental toxicity in zebrafish embryos (Danio rerio) induced by boscalid

Boscalid has been widely used for controlling various plant diseases. It is one of the most frequently detected pesticides in main coastal estuaries in California, with concentrations as high as 36μg/L, but its ecotoxicology information is scarce. To assess the aquatic risk of boscalid, acute toxicity and sub-lethal developmental toxicity toward zebrafish embryos were determined in the present study. In the acute toxicity test, a series of toxic symptoms of embryos were observed, including abnormal spontaneous movement, slow heartbeat, yolk sac oedema, pericardial oedema, spine deformation and hatching inhibition, and 96-h-LC50 (50% lethal concentration) of boscalid toward zebrafish embryos was 2.65 (2.506-2.848)mg/L. From the results of the sub-lethal developmental toxicity test, boscalid was confirmed to have a great impact on development mechanisms of zebrafish embryos. Cell apoptosis in embryos was induced by boscalid with upregulation of genes in the cell apoptosis and an increase of capspase-3 and caspase-9 activity in the present study. Lipid metabolism was affected in embryos due to changes in gene expression and the contents of total triacylglyceride and cholesterol. Melanin synthesis and deposition was caused in embryos due to alterations in related gene expression. Overall, changes in cell apoptosis, lipid metabolism and melanin synthesis and deposition might be responsible for developmental toxicity of boscalid to zebrafish embryos.

Combination of novel and public RNA-seq datasets to generate an mRNA expression atlas for the domestic chicken

Background

The domestic chicken (Gallus gallus) is widely used as a model in developmental biology and is also an important livestock species. We describe a novel approach to data integration to generate an mRNA expression atlas for the chicken spanning major tissue types and developmental stages, using a diverse range of publicly-archived RNA-seq datasets and new data derived from immune cells and tissues.

Results

Randomly down-sampling RNA-seq datasets to a common depth and quantifying expression against a reference transcriptome using the mRNA quantitation tool Kallisto ensured that disparate datasets explored comparable transcriptomic space. The network analysis tool Graphia was used to extract clusters of co-expressed genes from the resulting expression atlas, many of which were tissue or cell-type restricted, contained transcription factors that have previously been implicated in their regulation, or were otherwise associated with biological processes, such as the cell cycle. The atlas provides a resource for the functional annotation of genes that currently have only a locus ID. We cross-referenced the RNA-seq atlas to a publicly available embryonic Cap Analysis of Gene Expression (CAGE) dataset to infer the developmental time course of organ systems, and to identify a signature of the expansion of tissue macrophage populations during development.

Conclusion

Expression profiles obtained from public RNA-seq datasets – despite being generated by different laboratories using different methodologies – can be made comparable to each other. This meta-analytic approach to RNA-seq can be extended with new datasets from novel tissues, and is applicable to any species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4972-7) contains supplementary material, which is available to authorized users.

Rab32/38 and the xenophagic restriction of intracellular bacteria replication

Rab GTPases' subversion by intracellular pathogens during infection has been extensively documented. Recent findings have implicated a key intracellular bacterial restriction/containment function for Rab32/38 in Salmonella species in macrophages and Listeria monocytogenes in dendritic cells. Rab32/38 aids the phagolysosome maturation, and mediates a parallel xenophagy mechanism by engaging prohibitins.