Prostaglandin signaling regulates renal multiciliated cell specification and maturation.

Multiciliated cells (MCCs) are specialized epithelia with apical bundles of motile cilia that direct fluid flow. MCC dysfunction is associated with human diseases of the respiratory, reproductive, and central nervous systems. Further, the appearance of renal MCCs has been cataloged in several kidney conditions, where their function is unknown. Despite their pivotal health importance, many aspects of MCC development remain poorly understood. Here, we utilized a chemical screen to identify molecules that affect MCC ontogeny in the zebrafish embryo kidney, and found prostaglandin signaling is essential both for renal MCC progenitor formation and terminal differentiation. Moreover, we show that prostaglandin activity is required downstream of the transcription factor ets variant 5a (etv5a) during MCC fate choice, where modulating prostaglandin E2 (PGE2) levels rescued MCC number. The discovery that prostaglandin signaling mediates renal MCC development has broad implications for other tissues, and could provide insight into a multitude of pathological states.

Enhanced genome assembly and a new official gene set for Tribolium castaneum.

BACKGROUND: The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. RESULTS: Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. CONCLUSIONS: The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

MicroRNA evolution, expression, and function during short germband development in Tribolium castaneum.

MicroRNAs are well-established players in the development of multicellular animals. Most of our understanding of microRNA function in arthropod development comes from studies in Drosophila. Despite their advantages as model systems, the long germband embryogenesis of fruit flies is an evolutionary derived state restricted to several holometabolous insect lineages. MicroRNA evolution and expression across development in animals exhibiting the ancestral and more widespread short germband mode of embryogenesis has not been characterized. We sequenced small RNA libraries of oocytes and successive intervals covering the embryonic development of the short germband model organism, Tribolium castaneum. We analyzed the evolution and temporal expression of the microRNA complement and sequenced libraries of total RNA to investigate the relationships with microRNA target expression. We show microRNA maternal loading and sequence-specific 3' end nontemplate oligoadenylation of maternally deposited microRNAs that is conserved between Tribolium and Drosophila. We further uncover large clusters encoding multiple paralogs from several Tribolium-specific microRNA families expressed during a narrow interval of time immediately after the activation of zygotic transcription. These novel microRNAs, together with several early expressed conserved microRNAs, target a significant number of maternally deposited transcripts. Comparison with Drosophila shows that microRNA-mediated maternal transcript targeting is a conserved process in insects, but the number and sequences of microRNAs involved have diverged. The expression of fast-evolving and species-specific microRNAs in the early blastoderm of T. castaneum is consistent with previous findings in Drosophila and shows that the unique permissiveness for microRNA innovation at this stage is a conserved phenomenon.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution.

BACKGROUND: The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum. RESULTS: We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication. CONCLUSIONS: Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

Fast-evolving microRNAs are highly expressed in the early embryo of Drosophila virilis.

MicroRNAs are short non-protein-coding RNAs that regulate gene expression at the post-transcriptional level and are essential for the embryonic development of multicellular animals. Comparative genome-scale analyses have revealed that metazoan evolution is accompanied by the continuous acquisition of novel microRNA genes. This suggests that novel microRNAs may promote innovation and diversity in development. We determined the evolutionary origins of extant Drosophila microRNAs and estimated the sequence divergence between the 130 orthologous microRNAs in Drosophila melanogaster and Drosophila virilis, separated by 63 million years of evolution. We then generated small RNA sequencing data sets covering D. virilis development and explored the relationship between microRNA conservation and expression in a developmental context. We find that late embryonic, larval, and adult stages are dominated by conserved microRNAs. This pattern, however, does not hold for the early embryo, where rapidly evolving microRNAs are uniquely present at high levels in both species. The group of fast-evolving microRNAs that are highly expressed in the early embryo belong to two Drosophilid lineage-specific clusters: mir-310 ∼ 313 and mir-309 ∼ 6. These clusters have particularly complex evolutionary histories of duplication, gain, and loss. Our analyses suggest that the early embryo is a more permissive environment for microRNA changes and innovations. Fast-evolving microRNAs, therefore, have the opportunity to become preferentially integrated in early developmental processes, and may impact the evolution of development. The relationship between microRNA conservation and expression throughout the development of Drosophila differs from that previously observed for protein-coding genes.

Clusters of microRNAs emerge by new hairpins in existing transcripts.

Genetic linkage may result in the expression of multiple products from a polycistronic transcript, under the control of a single promoter. In animals, protein-coding polycistronic transcripts are rare. However, microRNAs are frequently clustered in the genomes of animals, and these clusters are often transcribed as a single unit. The evolution of microRNA clusters has been the subject of much speculation, and a selective advantage of clusters of functionally related microRNAs is often proposed. However, the origin of microRNA clusters has not been so far explored. Here, we study the evolution of microRNA clusters in Drosophila melanogaster. We observed that the majority of microRNA clusters arose by the de novo formation of new microRNA-like hairpins in existing microRNA transcripts. Some clusters also emerged by tandem duplication of a single microRNA. Comparative genomics show that these clusters are unlikely to split or undergo rearrangements. We did not find any instances of clusters appearing by rearrangement of pre-existing microRNA genes. We propose a model for microRNA cluster evolution in which selection over one of the microRNAs in the cluster interferes with the evolution of the other linked microRNAs. Our analysis suggests that the study of microRNAs and small RNAs must consider linkage associations.

Structure, evolution and function of the bi-directionally transcribed iab-4/iab-8 microRNA locus in arthropods.

In Drosophila melanogaster, the iab-4/iab-8 locus encodes bi-directionally transcribed microRNAs that regulate the function of flanking Hox transcription factors. We show that bi-directional transcription, temporal and spatial expression patterns and Hox regulatory function of the iab-4/iab-8 locus are conserved between fly and the beetle Tribolium castaneum. Computational predictions suggest iab-4 and iab-8 microRNAs can target common sites, and cell-culture assays confirm that iab-4 and iab-8 function overlaps on Hox target sites in both fly and beetle. However, we observe key differences in the way Hox genes are targeted. For instance, abd-A transcripts are targeted only by iab-8 in Drosophila, whereas both iab-4 and iab-8 bind to Tribolium abd-A. Our evolutionary and functional characterization of a bi-directionally transcribed microRNA establishes the iab-4/iab-8 system as a model for understanding how multiple products from sense and antisense microRNAs target common sites.

MicroRNA evolution by arm switching.

MicroRNAs (miRNAs) modulate transcript stability and translation. Functional mature miRNAs are processed from one or both arms of the hairpin precursor. The miR-100/10 family has undergone three independent evolutionary events that have switched the arm from which the functional miRNA is processed. The dominant miR-10 sequences in the insects Drosophila melanogaster and Tribolium castaneum are processed from opposite arms. However, the duplex produced by Dicer cleavage has an identical sequence in fly and beetle. Expression of the Tribolium miR-10 sequence in Drosophila S2 cells recapitulates the native beetle pattern. Thus, arm usage is encoded in the primary miRNA sequence, but outside the mature miRNA duplex. We show that the predicted messenger RNA targets and inferred function of sequences from opposite arms differ significantly. Arm switching is likely to be general, and provides a fundamental mechanism to evolve the function of a miRNA locus and target gene network.

miR-9a regulates levels of both rhomboid mRNA and protein in the early Drosophila melanogaster embryo.

MicroRNAs can have subtle and combinatorial effects on the levels of the targets and pathways they act on. Studying the consequences of a single microRNA knockout often proves difficult as many such knockouts exhibit phenotypes only under stress conditions. This has often led to the hypothesis that microRNAs buffer the effects of intrinsic and environmental stochasticity on gene expression. Observing and understanding this buffering effect entails quantitative analysis of microRNA and target expression in single cells. To this end, we have employed single-molecule fluorescence in situ hybridization, immunofluorescence, and high-resolution confocal microscopy to investigate the effects of miR-9a loss on the expression of the serine-protease Rhomboid in Drosophila melanogaster early embryos. Our single-cell quantitative approach shows that spatially, the rhomboid mRNA pattern is identical in WT and miR-9a knockout embryos. However, we find that the number of mRNA molecules per cell is higher when miR-9a is absent, and the level and temporal accumulation of rhomboid protein shows a more dramatic increase in the miR-9a knockout. Specifically, we see accumulation of rhomboid protein in miR-9a mutants by stage 5, much earlier than in WT. The data, therefore, show that miR-9a functions in the regulation of rhomboid mRNA and protein levels. While further work is required to establish whether rhomboid is a direct target of miR-9 in Drosophila, our results further establish the miR-9 family microRNAs as conserved regulators of timing in neurogenic processes. This study shows the power of single-cell quantification as an experimental tool to study phenotypic consequences of microRNA mis-regulation.

The embryonic transcriptome of Parhyale hawaiensis reveals different dynamics of microRNAs and mRNAs during the maternal-zygotic transition.

Parhyale hawaiensis has emerged as the crustacean model of choice due to its tractability, ease of imaging, sequenced genome, and development of CRISPR/Cas9 genome editing tools. However, transcriptomic datasets spanning embryonic development are lacking, and there is almost no annotation of non-protein-coding RNAs, including microRNAs. We have sequenced microRNAs, together with mRNAs and long non-coding RNAs, in Parhyale using paired size-selected RNA-seq libraries at seven time-points covering important transitions in embryonic development. Focussing on microRNAs, we annotate 175 loci in Parhyale, 88 of which have no known homologs. We use these data to annotate the microRNAome of 37 crustacean genomes, and suggest a core crustacean microRNA set of around 61 sequence families. We examine the dynamic expression of microRNAs and mRNAs during the maternal-zygotic transition. Our data suggest that zygotic genome activation occurs in two waves in Parhyale with microRNAs transcribed almost exclusively in the second wave. Contrary to findings in other arthropods, we do not predict a general role for microRNAs in clearing maternal transcripts. These data significantly expand the available transcriptomics resources for Parhyale, and facilitate its use as a model organism for the study of small RNAs in processes ranging from embryonic development to regeneration.

Age-Related Alterations in Macrophage Distribution and Function Are Associated With Delayed Cutaneous Wound Healing.

Ageing-related delays and dysregulated inflammation in wound healing are well-documented in both human and animal models. However, cellular and molecular changes underlying this impairment in healing progression are not fully understood. In this study, we characterised ageing-associated changes to macrophages in wounds of young and aged mice and investigated transcriptomic differences that may impact the progression of wound healing. Full-thickness wounds created on the dorsum of C57BL/6J young and aged mice were excised on Days 3 and 7 post-wounding for analysis by immunohistochemistry, flow cytometry, and RNA sequencing. Our data revealed that macrophages were significantly reduced in aged wounds in comparison to young. Functional transcriptomic analyses showed that macrophages from aged wounds exhibited significantly reduced expression of cell cycle, DNA replication, and repair pathway genes. Furthermore, we uncovered an elevated pro-inflammatory gene expression program in the aged macrophages correlated with poor inflammation resolution and excessive tissue damage observed in aged wounds. Altogether, our work provides insights into how poorly healing aged wounds are phenotypically defined by the presence of macrophages with reduced proliferative capacity and an exacerbated inflammatory response, both of which are pathways that can be targeted to improve healing in the elderly.

Evolution of the ventral midline in insect embryos.

The ventral midline is a source of signals that pattern the nerve cord of insect embryos. In dipterans such as the fruitfly Drosophila melanogaster (D. mel.) and the mosquito Anopheles gambiae (A. gam.), the midline is narrow and spans just 1-2 cells. However, in the honeybee, Apis mellifera (A. mel.), the ventral midline is broad and encompasses 5-6 cells. slit and other midline-patterning genes display a corresponding expansion in expression. Evidence is presented that this difference is due to divergent cis regulation of the single-minded (sim) gene, which encodes a bHLH-PAS transcription factor essential for midline differentiation. sim is regulated by a combination of Notch signaling and a Twist (Twi) activator gradient in D. mel., but it is activated solely by Twi in A. mel. We suggest that the Twi-only mode of regulation--and the broad ventral midline--represents the ancestral form of CNS patterning in Holometabolous insects.

smiFISH and embryo segmentation for single-cell multi-gene RNA quantification in arthropods.

Recently, advances in fluorescent in-situ hybridization techniques and in imaging technology have enabled visualization and counting of individual RNA molecules in single cells. This has greatly enhanced the resolution in our understanding of transcriptional processes. Here, we adapt a recently published smiFISH protocol (single-molecule inexpensive fluorescent in-situ hybridization) to whole embryos across a range of arthropod model species, and also to non-embryonic tissues. Using multiple fluorophores with distinct spectra and white light laser confocal imaging, we simultaneously detect and separate single RNAs from up to eight different genes in a whole embryo. We also combine smiFISH with cell membrane immunofluorescence, and present an imaging and analysis pipeline for 3D cell segmentation and single-cell RNA counting in whole blastoderm embryos. Finally, using whole embryo single-cell RNA count data, we propose two alternative single-cell variability measures to the commonly used Fano factor, and compare the capacity of these three measures to address different aspects of single-cell expression variability.

Single-cell visualization of mir-9a and Senseless co-expression during Drosophila melanogaster embryonic and larval peripheral nervous system development.

The Drosophila melanogaster peripheral nervous system (PNS) comprises the sensory organs that allow the fly to detect environmental factors such as temperature and pressure. PNS development is a highly specified process where each sensilla originates from a single sensory organ precursor (SOP) cell. One of the major genetic orchestrators of PNS development is Senseless, which encodes a zinc finger transcription factor (Sens). Sens is both necessary and sufficient for SOP differentiation. Senseless expression and SOP number are regulated by the microRNA miR-9a. However, the reciprocal dynamics of Senseless and miR-9a are still obscure. By coupling single-molecule FISH with immunofluorescence, we are able to visualize transcription of the mir-9a locus and expression of Sens simultaneously. During embryogenesis, we show that the expression of mir-9a in SOP cells is rapidly lost as Senseless expression increases. However, this mutually exclusive expression pattern is not observed in the third instar imaginal wing disk, where some Senseless-expressing cells show active sites of mir-9a transcription. These data challenge and extend previous models of Senseless regulation and show complex co-expression dynamics between mir-9a and Senseless. The differences in this dynamic relationship between embryonic and larval PNS development suggest a possible switch in miR-9a function. Our work brings single-cell resolution to the understanding of dynamic regulation of PNS development by Senseless and miR-9a.

Silencing of an abdominal Hox gene during early development is correlated with limb development in a crustacean trunk.

We tested whether Artemia abd-A could repress limbs in Drosophila embryos, and found that although abd-A transcripts were produced, ABD-A protein was not. Similarly, developing Artemia epidermal cells showed expression of abd-A transcripts without accumulation of ABD-A protein. This finding in Artemia reveals a new variation in Hox gene function that is associated with morphological evolution. In this case, a HOX protein expression pattern is completely absent during early development, although the HOX protein is expressed at later stages in the central nervous system in a "homeotic-like" pattern. The combination of an absence of ABD-A protein expression in the Artemia limb primordia and the weak repressive function of Artemia UBX protein on the limb-promoting gene Dll are likely to be two reasons why homonomous limbs develop throughout the entire Artemia trunk.

The Transcription Factor-microRNA Regulatory Network during hESC-chondrogenesis.

Human embryonic stem cells (ESCs) offer a promising therapeutic approach for osteoarthritis (OA). The unlimited source of cells capable of differentiating to chondrocytes has potential for repairing damaged cartilage or to generate disease models via gene editing. However their use is limited by the efficiency of chondrogenic differentiation. An improved understanding of the transcriptional and post-transcriptional regulation of chondrogenesis will enable us to improve hESC chondrogenic differentiation protocols. Small RNA-seq and whole transcriptome sequencing was performed on distinct stages of hESC-directed chondrogenesis. This revealed significant changes in the expression of several microRNAs including upregulation of known cartilage associated microRNAs and those transcribed from the Hox complexes, and the downregulation of pluripotency associated microRNAs. Integration of miRomes and transcriptomes generated during hESC-directed chondrogenesis identified key functionally related clusters of co-expressed microRNAs and protein coding genes, associated with pluripotency, primitive streak, limb development and extracellular matrix. Analysis identified regulators of hESC-directed chondrogenesis such as miR-29c-3p with 10 of its established targets identified as co-regulated 'ECM organisation' genes and miR-22-3p which is highly co-expressed with ECM genes and may regulate these genes indirectly by targeting the chondrogenic regulators SP1 and HDAC4. We identified several upregulated transcription factors including HOXA9/A10/D13 involved in limb patterning and RELA, JUN and NFAT5, which have targets enriched with ECM associated genes. We have developed an unbiased approach for integrating transcriptome and miRome using protein-protein interactions, transcription factor regulation and miRNA target interactions and identified key regulatory networks prominent in hESC chondrogenesis.

Visualization of trans-homolog enhancer-promoter interactions at the Abd-B Hox locus in the Drosophila embryo.

The Hox gene Abdominal-B (Abd-B) controls the morphogenesis of posterior abdominal segments in Drosophila. Expression is regulated by a series of 3' enhancers that are themselves transcribed. RNA FISH was used to visualize nascent transcripts associated with coding and noncoding regions of Abd-B in developing embryos. Confocal imaging suggests that distal enhancers often loop to the Abd-B promoter region. Surprisingly, enhancers located on one chromosome frequently associate with the Abd-B transcription unit located on the other homolog. These trans-homolog interactions can be interpreted as the direct visualization of a genetic phenomenon known as transvection, whereby certain mutations in Abd-B can be rescued in trans by the other copy of the gene. A 10 kb sequence in the 3' flanking region mediates tight pairing of Abd-B alleles, thereby facilitating trans looping of distal enhancers. Such trans-homolog interactions might be a common mechanism of gene regulation in higher metazoans.

Visualizing gene expression during zebrafish pronephros development and regeneration.

The vertebrate kidney is comprised of functional units known as nephrons. Defects in nephron development or activity are a common feature of kidney disease. Current medical treatments are unable to ameliorate the dire consequences of nephron deficit or injury. Although there have been tremendous advancements in our understanding of nephron ontogeny and the response to damage, many significant knowledge gaps still remain. The zebrafish embryo kidney, or pronephros, is an ideal model for many renal development and regeneration studies because it is comprised of nephrons that share conserved features with the nephron units that comprise the mammalian metanephric kidney. In this chapter, we provide an overview about the benefits of using the zebrafish pronephros to study the mechanisms underlying nephrogenesis as well as epithelial repair and regeneration. We subsequently detail methods for the spatiotemporal assessment of gene and protein expression in zebrafish embryos that can be used to extend the understanding of nephron development and disease, and thereby create new opportunities to identify therapeutic strategies for regenerative medicine.

Abundant expression of somatic transposon-derived piRNAs throughout Tribolium castaneum embryogenesis.

BACKGROUND: Piwi-interacting RNAs (piRNAs) are a class of short (~26-31-nucleotide) non-protein-coding RNAs expressed in the metazoan germline. The piRNA pathway in arthropods is best understood in the ovary of Drosophila melanogaster, where it acts to silence active transposable elements (TEs). Maternal loading of piRNAs in oocytes is further required for the inheritance of piRNA-mediated transposon defence. However, our understanding of the diversity, evolution and function of the piRNA complement beyond drosophilids is limited. The red flour beetle, Tribolium castaneum, is an emerging model organism separated from Drosophila by ~ 350 million years of evolution that displays a number of features ancestral to arthropods, including short germ embryogenesis. Here, we characterize the maternally deposited and zygotically expressed small RNA and mRNA complements throughout T. castaneum embryogenesis. RESULTS: We find that beetle oocytes and embryos of all stages are abundant in heterogeneous ~ 28-nucleotide RNAs. These small RNAs originate from discrete genomic loci enriched in TE sequences and display the molecular signatures of transposon-derived piRNAs. In addition to the maternally loaded primary piRNAs, Tribolium embryos produce secondary piRNAs by the cleavage of zygotically activated TE transcripts via the ping-pong mechanism. The two Tribolium piRNA pathway effector proteins, Tc-Piwi/Aub and Tc-Ago3, are also expressed throughout the soma of early embryos. CONCLUSIONS: Our results show that the piRNA pathway in Tribolium is not restricted to the germline, but also operates in the embryo and may act to antagonize zygotically activated transposons. Taken together, these data highlight a functional divergence of the piRNA pathway between insects.

Functional and Genetic Analysis of Spectraplakins in Drosophila.

The cytoskeleton is a dynamic network of filamentous protein polymers required for virtually all cellular processes. It consists of three major classes, filamentous actin (F-actin), intermediate filaments, and microtubules, all displaying characteristic structural properties, functions, cellular distributions, and sets of interacting regulatory proteins. One unique class of proteins, the spectraplakins, bind, regulate, and integrate the functions of all three classes of cytoskeleton proteins. Spectraplakins are giant, evolutionary conserved multidomain proteins (spanning up to 9000 aa) that are true members of the plakin, spectrin, and Gas2-like protein families. They have OMIM-listed disease links to epidermolysis bullosa and hereditary sensory and autonomic neuropathy. Their role in disease is likely underrepresented since studies in model animal systems have revealed critical roles in polarity, morphogenesis, differentiation and maintenance, migration, signaling, and intracellular trafficking in a variety of tissues. This enormous diversity of spectraplakin function is consistent with the numerous isoforms produced from single genomic loci that combine different sets of functional domains in distinct cellular contexts. To study the broad range of functions and complexity of these proteins, Drosophila is a powerful model. Thus, the fly spectraplakin Short stop (Shot) acts as an actin-microtubule linker and plays important roles in many developmental processes, which provide experimentally amenable and relevant contexts in which to study spectraplakin functions. For these studies, a versatile range of relevant experimental resources that facilitate genetics and transgenic approaches, highly refined genomics tools, and an impressive set of spectraplakin-specific genetic and molecular tools are readily available. Here, we use the example of Shot to illustrate how the various tools and strategies available for Drosophila can be employed to decipher and dissect cellular roles and molecular mechanisms of spectraplakins.