Transcriptomics of developing embryos and organs: A raising tool for evo-devo

Deciphering the evolutionary development of the "Chinese lantern" within Solanaceae

MAIN CONCLUSION

The key genetic variation underlying the evo-devo of ICS in Solanaceae may be further pinpointed using an integrated strategy of forward and reverse genetics studies under the framework of phylogeny. The calyx of Physalis remains persistent throughout fruit development. Post-flowering, the fruiting calyx is inflated rapidly to encapsulate the berry, giving rise to a "Chinese lantern" structure called inflated calyx syndrome (ICS). It is unclear how this novelty arises. Over the past 2 decades, the role of MADS-box genes in the evolutionary development (evo-devo) of ICS has mainly been investigated within Solanaceae. In this review, we analyze the main achievements, challenges, and new progress. ICS acts as a source for fruit development, provides a microenvironment to protect fruit development, and assists in long-distance fruit dispersal. ICS is a typical post-floral trait, and the onset of its development is triggered by specific developmental signals that coincide with fertilization. These signals can be replaced by exogenous gibberellin and cytokinin application. MPF2-like heterotopic expression and MBP21-like loss have been proposed to be two essential evolutionary events for ICS origin, and manipulating the related MADS-box genes has been shown to affect the ICS size, sepal organ identity, and/or male fertility, but not completely disrupt ICS. Therefore, the core genes or key links in the ICS biosynthesis pathways may have undergone secondary mutations during evolution, or they have not yet been pinpointed. Recently, we have made some encouraging progress in acquiring lantern mutants in Physalis floridana. In addition to technological innovation, we propose an integrated strategy to further analyze the evo-devo mechanisms of ICS in Solanaceae using forward and reverse genetics studies under the framework of phylogeny.

Transcriptome analysis of East Asian common octopus, Octopus sinensis, paralarvae

BACKGROUND

The genes involved in cephalopod development and their association with hatching and survival during early life stages have been extensively studied. However, few studies have investigated the paralarvae transcriptome of the East Asian common octopus (Octopus sinen sis).

OBJECTIVE

This study aimed to identify the genes related to embryonic development and hatching in O. sinensis using RNA sequencing (RNA-seq) and verify the genes most relevant to different embryonic stages.

METHODS

RNA samples from hatched and 25 days post-hatching (dph) O. sinensis paralarvae were used to construct cDNA libraries. Clean reads from individual samples were aligned to the reference O. sinensis database to identify the differentially expressed genes (DEGs) between the 0- and 25-dph paralarvae libraries. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to supplement the RNA-seq data for embryogenic developmental stages.

RESULTS

A total of 12,597 transcripts were annotated and 5,468 DEGs were identified between the 0- and 25-dph O. sinensis paralarvae, including 2,715 upregulated and 2,753 downregulated transcripts in the 25-dph paralarvae. Several key DEGs were related to transmembrane transport, lipid biosynthesis, monooxygenase activity, lipid transport, neuropeptide signaling, transcription regulation, and protein-cysteine S-palmitoyltransferase activity during the post-hatching development of O. sinensis paralarvae. RT-qPCR analysis further revealed that SLC5A3A, ABCC12, and NPC1 transcripts in 20 and/or 30 days post-fertilization (dpf) embryos were significantly higher (p < 0.05) than those in 10-dpf embryos.

CONCLUSION

Transcriptome profiles provide molecular targets to understand the embryonic development, hatching, and survival of O. sinensis paralarvae, and enhance octopus production.

Sex-biased gene expression across mammalian organ development and evolution

Sexually dimorphic traits are common among mammals and are specified during development through the deployment of sex-specific genetic programs. Because little is known about these programs, we investigated them using a resource of gene expression profiles in males and females throughout the development of five organs in five mammals (human, mouse, rat, rabbit, and opossum) and a bird (chicken). We found that sex-biased gene expression varied considerably across organs and species and was often cell-type specific. Sex differences increased abruptly around sexual maturity instead of increasing gradually during organ development. Finally, sex-biased gene expression evolved rapidly at the gene level, with differences between organs in the evolutionary mechanisms used, but more slowly at the cellular level, with the same cell types being sexually dimorphic across species.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.

Developmental Gene Expression Differences between Humans and Mammalian Models

Identifying the molecular programs underlying human organ development and how they differ from model species is key for understanding human health and disease. Developmental gene expression profiles provide a window into the genes underlying organ development and a direct means to compare them across species. We use a transcriptomic resource covering the development of seven organs to characterize the temporal profiles of human genes associated with distinct disease classes and to determine, for each human gene, the similarity of its spatiotemporal expression with its orthologs in rhesus macaque, mouse, rat, and rabbit. We find clear associations between spatiotemporal profiles and the phenotypic manifestations of diseases. We also find that half of human genes differ from their mouse orthologs in their temporal trajectories in at least one of the organs. These include more than 200 genes associated with brain, heart, and liver disease for which mouse models should undergo extra scrutiny.

RNA-Seq Analysis Reveals Hub Genes Involved in Chicken Intramuscular Fat and Abdominal Fat Deposition During Development

Fat traits are important in the chicken industry where there is a desire for high intramuscular fat (IMF) and low abdominal fat. However, there is limited knowledge on the relationship between the dynamic status of gene expression and the body fat deposition in chicken. Transcriptome data were obtained from breast muscle and abdominal fat of female chickens from nine developmental stages (from embryonic day 12 to hatched day 180). In total, 8,545 genes in breast muscle and 6,824 genes in abdominal fat were identified as developmentally dynamic genes. Weighted correlation network analysis was used to identify gene modules and the hub genes. Twenty-one hub genes were identified, e.g., ENSGALG00000041996, which represents a candidate for high IMF, and CREB3L1, which relates to low abdominal fat weight. The transcript factor L3MBTL1 and the transcript factor cofactors TNIP1, HAT1, and BEND6 related to both high breast muscle IMF and low abdominal fat weight. Our results provide a resource of developmental transcriptome profiles in chicken breast muscle and abdominal fat. The candidate genes can be used in the selection for increased IMF content and/or a decrease in abdominal fat weight which would contribute to the improvement of these traits.

Gene expression trajectories during male and female reproductive development in balsam poplar (Populus balsamifera L.)

Plant reproductive development from the first appearance of reproductively committed axes through to floral maturation requires massive and rapid remarshalling of gene expression. In dioecious species such as poplar this is further complicated by divergent male and female developmental programs. We used seven time points in male and female balsam poplar (Populus balsamifera L.) buds and catkins representing the full annual flowering cycle, to elucidate the effects of time and sex on gene expression during reproductive development. Time (developmental stage) is dominant in patterning gene expression with the effect of sex nested within this. Here, we find (1) evidence for five successive waves of alterations to the chromatin landscape which may be important in setting the overall reproductive trajectory, regardless of sex. (2) Each individual developmental stage is further characterized by marked sex-differential gene expression. (3) Consistent sexually differentiated gene expression regardless of developmental stage reveal candidates for high-level regulators of sex and include the female-specific poplar ARR17 homologue. There is also consistent male-biased expression of the MADS-box genes PISTILLATA and APETALA3. Our work provides insights into expression trajectories shaping reproductive development, its potential underlying mechanisms, and sex-specific translation of the genome information into reproductive structures in balsam poplar.

Gene expression across mammalian organ development

The evolution of gene expression in mammalian organ development remains largely uncharacterized. Here we report the transcriptomes of seven organs (cerebrum, cerebellum, heart, kidney, liver, ovary and testis) across developmental time points from early organogenesis to adulthood for human, rhesus macaque, mouse, rat, rabbit, opossum and chicken. Comparisons of gene expression patterns identified correspondences of developmental stages across species, and differences in the timing of key events during the development of the gonads. We found that the breadth of gene expression and the extent of purifying selection gradually decrease during development, whereas the amount of positive selection and expression of new genes increase. We identified differences in the temporal trajectories of expression of individual genes across species, with brain tissues showing the smallest percentage of trajectory changes, and the liver and testis showing the largest. Our work provides a resource of developmental transcriptomes of seven organs across seven species, and comparative analyses that characterize the development and evolution of mammalian organs.

Perspectives and applications of machine learning for evolutionary developmental biology

Evolutionary Developmental Biology (Evo-Devo) is an ever-expanding field that aims to understand how development was modulated by the evolutionary process. In this sense, "omic" studies emerged as a powerful ally to unravel the molecular mechanisms underlying development. In this scenario, bioinformatics tools become necessary to analyze the growing amount of information. Among computational approaches, machine learning stands out as a promising field to generate knowledge and trace new research perspectives for bioinformatics. In this review, we aim to expose the current advances of machine learning applied to evolution and development. We draw clear perspectives and argue how evolution impacted machine learning techniques.

Selective Constraints on Coding Sequences of Nervous System Genes Are a Major Determinant of Duplicate Gene Retention in Vertebrates

The evolutionary history of vertebrates is marked by three ancient whole-genome duplications: two successive rounds in the ancestor of vertebrates, and a third one specific to teleost fishes. Biased loss of most duplicates enriched the genome for specific genes, such as slow evolving genes, but this selective retention process is not well understood. To understand what drives the long-term preservation of duplicate genes, we characterized duplicated genes in terms of their expression patterns. We used a new method of expression enrichment analysis, TopAnat, applied to in situ hybridization data from thousands of genes from zebrafish and mouse. We showed that the presence of expression in the nervous system is a good predictor of a higher rate of retention of duplicate genes after whole-genome duplication. Further analyses suggest that purifying selection against the toxic effects of misfolded or misinteracting proteins, which is particularly strong in nonrenewing neural tissues, likely constrains the evolution of coding sequences of nervous system genes, leading indirectly to the preservation of duplicate genes after whole-genome duplication. Whole-genome duplications thus greatly contributed to the expansion of the toolkit of genes available for the evolution of profound novelties of the nervous system at the base of the vertebrate radiation.

Sex differentiation in grayling (Salmonidae) goes through an all-male stage and is delayed in genetic males who instead grow faster

Fish populations can be threatened by distorted sex ratios that arise during sex differentiation. Here we describe sex differentiation in a wild grayling (Thymallus thymallus) population that suffers from distorted sex ratios. We verified that sex determination is linked to the sex determining locus (sdY) of salmonids. This allowed us to study sex-specific gene expression and gonadal development. Sex-specific gene expression could be observed during embryogenesis and was strong around hatching. About half of the fish showed immature testes around eleven weeks after fertilization. This phenotype was mostly replaced by the "testis-to-ovary" or "ovaries" phenotypes during development. The gonads of the remaining fish stayed undifferentiated until six months after fertilization. Genetic sexing revealed that fish with undifferentiated gonads were all males, who grew larger than the genetic females during the observational period. Only 12% of the genetic males showed testicular tissue six months after fertilization. We conclude that sex differentiation starts before hatching, goes through an all-male stage for both sexes (which represents a rare case of "undifferentiated" gonochoristic species that usually go through an all-female stage), and is delayed in males. During these juvenile stages males grow faster than females instead of developing their gonads.

Development and Evolution through the Lens of Global Gene Regulation

Evolution and development are two inherently intertwined processes. As the embryo develops it does so in ways that both reflect past constraints and bias the future evolution of the species. While research exploiting this insight typically studies individual genes, transcriptomic analyses have sparked a new wave of discoveries. In this opinion piece, I review the evidence arising from transcriptomics on the topics of the evolution of germ layers, the phylotypic stage, and developmental constraints. The spatiotemporal pattern of gene expression across germ layers provides evidence that the endoderm was the first germ layer to evolve. Comparing transcriptome dynamics throughout developmental time across distant species reveals a mid-developmental transition under strong developmental constraints. These studies highlight the efficiency of exploratory data analysis using computational tools and comparative approaches for discovery.

Optimal alpha reduces error rates in gene expression studies: a meta-analysis approach

Background

Transcriptomic approaches (microarray and RNA-seq) have been a tremendous advance for molecular science in all disciplines, but they have made interpretation of hypothesis testing more difficult because of the large number of comparisons that are done within an experiment. The result has been a proliferation of techniques aimed at solving the multiple comparisons problem, techniques that have focused primarily on minimizing Type I error with little or no concern about concomitant increases in Type II errors. We have previously proposed a novel approach for setting statistical thresholds with applications for high throughput omics-data, optimal α, which minimizes the probability of making either error (i.e. Type I or II) and eliminates the need for post-hoc adjustments.

Results

A meta-analysis of 242 microarray studies extracted from the peer-reviewed literature found that current practices for setting statistical thresholds led to very high Type II error rates. Further, we demonstrate that applying the optimal α approach results in error rates as low or lower than error rates obtained when using (i) no post-hoc adjustment, (ii) a Bonferroni adjustment and (iii) a false discovery rate (FDR) adjustment which is widely used in transcriptome studies.

Conclusions

We conclude that optimal α can reduce error rates associated with transcripts in both microarray and RNA-seq experiments, but point out that improved statistical techniques alone cannot solve the problems associated with high throughput datasets – these approaches need to be coupled with improved experimental design that considers larger sample sizes and/or greater study replication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1728-3) contains supplementary material, which is available to authorized users.

Cross-kingdom comparison of the developmental hourglass

The developmental hourglass model has its foundations in classic anatomical studies by von Baer and Haeckel. In this context, even the conservation of animal body plans has been explained by evolutionary constraints acting on mid-embryogenic development. Recent studies have shown that developmental hourglass patterns also exist on the transcriptomic level, mirroring the corresponding morphological patterns. The identification of similar patterns in embryonic, post-embryonic, and life cycle spanning transcriptomes in plant and fungus development, however, contradict the notion of a direct coupling between morphological and molecular patterns. To explain the existence of hourglass patterns across kingdoms and developmental processes, we propose the organizational checkpoint model that integrates the developmental hourglass model into a framework of transcriptome switches.

Transcriptomic signatures shaped by cell proportions shed light on comparative developmental biology

Background

Comparative transcriptomics can answer many questions in developmental and evolutionary developmental biology. Most transcriptomic studies start by showing global patterns of variation in transcriptomes that differ between species or organs through developmental time. However, little is known about the kinds of expression differences that shape these patterns.

Results

We compared transcriptomes during the development of two morphologically distinct serial organs, the upper and lower first molars of the mouse. We found that these two types of teeth largely share the same gene expression dynamics but that three major transcriptomic signatures distinguish them, all of which are shaped by differences in the relative abundance of different cell types. First, lower/upper molar differences are maintained throughout morphogenesis and stem from differences in the relative abundance of mesenchyme and from constant differences in gene expression within tissues. Second, there are clear time-shift differences in the transcriptomes of the two molars related to cusp tissue abundance. Third, the transcriptomes differ most during early-mid crown morphogenesis, corresponding to exaggerated morphogenetic processes in the upper molar involving fewer mitotic cells but more migrating cells. From these findings, we formulate hypotheses about the mechanisms enabling the two molars to reach different phenotypes. We also successfully applied our approach to forelimb and hindlimb development.

Conclusions

Gene expression in a complex tissue reflects not only transcriptional regulation but also abundance of different cell types. This knowledge provides valuable insights into the cellular processes underpinning differences in organ development. Our approach should be applicable to most comparative developmental contexts.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1157-7) contains supplementary material, which is available to authorized users.

Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy

BACKGROUND

Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy.

RESULTS

Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism.

CONCLUSIONS

This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.

Dinosaurs, Chameleons, Humans, and Evo-Devo Path: Linking Étienne Geoffroy's Teratology, Waddington's Homeorhesis, Alberch's Logic of "Monsters," and Goldschmidt Hopeful "Monsters"

Since the rise of evo-devo (evolutionary developmental biology) in the 1980s, few authors have attempted to combine the increasing knowledge obtained from the study of model organisms and human medicine with data from comparative anatomy and evolutionary biology in order to investigate the links between development, pathology, and macroevolution. Fortunately, this situation is slowly changing, with a renewed interest in evolutionary developmental pathology (evo-devo-path) in the past decades, as evidenced by the idea to publish this special, and very timely, issue on "Developmental Evolution in Biomedical Research." As all of us have recently been involved, independently, in works related in some way or another with evolution and developmental anomalies, we decided to join our different perspectives and backgrounds in the present contribution for this special issue. Specifically, we provide a brief historical account on the study of the links between evolution, development, and pathologies, followed by a review of the recent work done by each of us, and then by a general discussion on the broader developmental and macroevolutionary implications of our studies and works recently done by other authors. Our primary aims are to highlight the strength of studying developmental anomalies within an evolutionary framework to understand morphological diversity and disease by connecting the recent work done by us and others with the research done and broader ideas proposed by authors such as Étienne Geoffroy Saint-Hilaire, Waddington, Goldschmidt, Gould, and Per Alberch, among many others to pave the way for further and much needed work regarding abnormal development and macroevolution.

Using gene networks in EvoDevo analyses

An intricate web of regulatory relationships between DNA, RNA, proteins and metabolites regulates how organisms achieve form and function. Genome sequencing combined with computational methods has allowed us to look at diverse readouts and generate a comprehensive framework for how molecules generate morphological phenotypes. RNAseq has evolved and proved useful for identifying links between transcription factor activity and transcript abundance, and for the generation of transcriptomes in non-model species through de novo assembly. Gene coexpression networks, combined with differential correlation analysis, offer the most versatile gene interaction exploratory tools, using gene expression patterns to determine potential associations and modularity of gene interactions and the differential associations between networks. Networks incorporating different types of biological data can help reduce the complexity of the entirety of biological information into discrete and interpretable pieces of data that can be the starting point for hypothesis construction and testing.

TROM: A Testing-Based Method for Finding Transcriptomic Similarity of Biological Samples

Comparative transcriptomics has gained increasing popularity in genomic research thanks to the development of high-throughput technologies including microarray and next-generation RNA sequencing that have generated numerous transcriptomic data. An important question is to understand the conservation and divergence of biological processes in different species. We propose a testing-based method TROM (Transcriptome Overlap Measure) for comparing transcriptomes within or between different species, and provide a different perspective, in contrast to traditional correlation analyses, about capturing transcriptomic similarity. Specifically, the TROM method focuses on identifying associated genes that capture molecular characteristics of biological samples, and subsequently comparing the biological samples by testing the overlap of their associated genes. We use simulation and real data studies to demonstrate that TROM is more powerful in identifying similar transcriptomes and more robust to stochastic gene expression noise than Pearson and Spearman correlations. We apply TROM to compare the developmental stages of six Drosophila species, C. elegans, S. purpuratus, D. rerio and mouse liver, and find interesting correspondence patterns that imply conserved gene expression programs in the development of these species. The TROM method is available as an R package on CRAN (https://cran.r-project.org/package=TROM) with manuals and source codes available at http://www.stat.ucla.edu/~jingyi.li/software-and-data/trom.html.

The analysis of Gene Regulatory Networks in plant evo-devo

We provide an overview of methods and workflows that can be used to investigate the topologies of Gene Regulatory Networks (GRNs) in the context of plant evolutionary-developmental (evo-devo) biology. Many of the species that occupy key positions in plant phylogeny are poorly adapted as laboratory models and so we focus here on techniques that can be efficiently applied to both model and non-model species of interest to plant evo-devo. We outline methods that can be used to describe gene expression patterns and also to elucidate the transcriptional, post-transcriptional, and epigenetic regulatory mechanisms underlying these patterns, in any plant species with a sequenced genome. We furthermore describe how the technique of Protein Resurrection can be used to confirm inferences on ancestral GRNs and also to provide otherwise-inaccessible points of reference in evolutionary histories by exploiting paralogues generated in gene and whole genome duplication events. Finally, we argue for the better integration of molecular data with information from paleobotanical, paleoecological, and paleogeographical studies to provide the fullest possible picture of the processes that have shaped the evolution of plant development.

On the evolutionary relationship between chondrocytes and osteoblasts

Vertebrates are the only animals that produce bone, but the molecular genetic basis for this evolutionary novelty remains obscure. Here, we synthesize information from traditional evolutionary and modern molecular genetic studies in order to generate a working hypothesis on the evolution of the gene regulatory network (GRN) underlying bone formation. Since transcription factors are often core components of GRNs (i.e., kernels), we focus our analyses on Sox9 and Runx2. Our argument centers on three skeletal tissues that comprise the majority of the vertebrate skeleton: immature cartilage, mature cartilage, and bone. Immature cartilage is produced during early stages of cartilage differentiation and can persist into adulthood, whereas mature cartilage undergoes additional stages of differentiation, including hypertrophy and mineralization. Functionally, histologically, and embryologically, these three skeletal tissues are very similar, yet unique, suggesting that one might have evolved from another. Traditional studies of the fossil record, comparative anatomy and embryology demonstrate clearly that immature cartilage evolved before mature cartilage or bone. Modern molecular approaches show that the GRNs regulating differentiation of these three skeletal cell fates are similar, yet unique, just like the functional and histological features of the tissues themselves. Intriguingly, the Sox9 GRN driving cartilage formation appears to be dominant to the Runx2 GRN of bone. Emphasizing an embryological and evolutionary transcriptomic view, we hypothesize that the Runx2 GRN underlying bone formation was co-opted from mature cartilage. We discuss how modern molecular genetic experiments, such as comparative transcriptomics, can test this hypothesis directly, meanwhile permitting levels of constraint and adaptation to be evaluated quantitatively. Therefore, comparative transcriptomics may revolutionize understanding of not only the clade-specific evolution of skeletal cells, but also the generation of evolutionary novelties, providing a modern paradigm for the evolutionary process.

Using the phenoscape knowledgebase to relate genetic perturbations to phenotypic evolution

The abundance of phenotypic diversity among species can enrich our knowledge of development and genetics beyond the limits of variation that can be observed in model organisms. The Phenoscape Knowledgebase (KB) is designed to enable exploration and discovery of phenotypic variation among species. Because phenotypes in the KB are annotated using standard ontologies, evolutionary phenotypes can be compared with phenotypes from genetic perturbations in model organisms. To illustrate the power of this approach, we review the use of the KB to find taxa showing evolutionary variation similar to that of a query gene. Matches are made between the full set of phenotypes described for a gene and an evolutionary profile, the latter of which is defined as the set of phenotypes that are variable among the daughters of any node on the taxonomic tree. Phenoscape's semantic similarity interface allows the user to assess the statistical significance of each match and flags matches that may only result from differences in annotation coverage between genetic and evolutionary studies. Tools such as this will help meet the challenge of relating the growing volume of genetic knowledge in model organisms to the diversity of phenotypes in nature. The Phenoscape KB is available at http://kb.phenoscape.org.

BrAD-seq: Breath Adapter Directional sequencing: a streamlined, ultra-simple and fast library preparation protocol for strand specific mRNA library construction

Next Generation Sequencing (NGS) is driving rapid advancement in biological understanding and RNA-sequencing (RNA-seq) has become an indispensable tool for biology and medicine. There is a growing need for access to these technologies although preparation of NGS libraries remains a bottleneck to wider adoption. Here we report a novel method for the production of strand specific RNA-seq libraries utilizing the terminal breathing of double-stranded cDNA to capture and incorporate a sequencing adapter. Breath Adapter Directional sequencing (BrAD-seq) reduces sample handling and requires far fewer enzymatic steps than most available methods to produce high quality strand-specific RNA-seq libraries. The method we present is optimized for 3-prime Digital Gene Expression (DGE) libraries and can easily extend to full transcript coverage shotgun (SHO) type strand-specific libraries and is modularized to accommodate a diversity of RNA and DNA input materials. BrAD-seq offers a highly streamlined and inexpensive option for RNA-seq libraries.

What to compare and how: Comparative transcriptomics for Evo-Devo

Evolutionary developmental biology has grown historically from the capacity to relate patterns of evolution in anatomy to patterns of evolution of expression of specific genes, whether between very distantly related species, or very closely related species or populations. Scaling up such studies by taking advantage of modern transcriptomics brings promising improvements, allowing us to estimate the overall impact and molecular mechanisms of convergence, constraint or innovation in anatomy and development. But it also presents major challenges, including the computational definitions of anatomical homology and of organ function, the criteria for the comparison of developmental stages, the annotation of transcriptomics data to proper anatomical and developmental terms, and the statistical methods to compare transcriptomic data between species to highlight significant conservation or changes. In this article, we review these challenges, and the ongoing efforts to address them, which are emerging from bioinformatics work on ontologies, evolutionary statistics, and data curation, with a focus on their implementation in the context of the development of our database Bgee (http://bgee.org). J. Exp. Zool. (Mol. Dev. Evol.) 324B: 372–382, 2015. © 2015 The Authors. J. Exp. Zool. (Mol. Dev. Evol.) published by Wiley Periodicals, Inc.

BrAD-seq: Breath Adapter Directional sequencing: a streamlined, ultra-simple and fast library preparation protocol for strand specific mRNA library construction

Next Generation Sequencing (NGS) is driving rapid advancement in biological understanding and RNA-sequencing (RNA-seq) has become an indispensable tool for biology and medicine. There is a growing need for access to these technologies although preparation of NGS libraries remains a bottleneck to wider adoption. Here we report a novel method for the production of strand specific RNA-seq libraries utilizing the terminal breathing of double-stranded cDNA to capture and incorporate a sequencing adapter. Breath Adapter Directional sequencing (BrAD-seq) reduces sample handling and requires far fewer enzymatic steps than most available methods to produce high quality strand-specific RNA-seq libraries. The method we present is optimized for 3-prime Digital Gene Expression (DGE) libraries and can easily extend to full transcript coverage shotgun (SHO) type strand-specific libraries and is modularized to accommodate a diversity of RNA and DNA input materials. BrAD-seq offers a highly streamlined and inexpensive option for RNA-seq libraries.