Gene duplication and the origins of morphological complexity in pancrustacean eyes, a genomic approach

Efficient CRISPR genome editing and integrative genomic analyses reveal the mosaicism of Cas-induced mutations and pleiotropic effects of scarlet gene in an emerging model system

Despite the revolutionary impacts of CRISPR-Cas gene editing systems, the effective and widespread use of CRISPR technologies in emerging model organisms still faces significant challenges. These include the inefficiency in generating heritable mutations at the organismal level, limited knowledge about the genomic consequences of gene editing, and an inadequate understanding of the inheritance patterns of CRISPR-Cas-induced mutations. This study addresses these issues by 1) developing an efficient microinjection delivery method for CRISPR editing in the microcrustacean Daphnia pulex; 2) assessing the editing efficiency of Cas9 and Cas12a nucleases, examining mutation inheritance patterns, and analyzing the local and global mutation spectrum in the scarlet mutants; and 3) investigating the transcriptomes of scarlet mutants to understand the pleiotropic effects of scarlet underlying their swimming behavior changes. Our reengineered CRISPR microinjection method results in efficient biallelic editing with both nucleases. While indels are dominant in Cas-induced mutations, a few on-site large deletions (>1kb) are observed, most likely caused by microhomology-mediated end joining repair. Knock-in of a stop codon cassette to the scarlet locus was successful, despite complex induced mutations surrounding the target site. Moreover, extensive germline mosaicism exists in some mutants, which unexpectedly produce different phenotypes/genotypes in their asexual progenies. Lastly, our transcriptomic analyses unveil significant gene expression changes associated with scarlet knock-out and altered swimming behavior in mutants, including several genes (e.g., NMDA1, ABAT, CNTNAP2) involved in human neurodegenerative diseases. This study expands our understanding of the dynamics of gene editing in the tractable model organism Daphnia and highlights its promising potential as a neurological disease model.

Origins and evolution of biological novelty

Understanding the origins and impacts of novel traits has been a perennial interest in many realms of ecology and evolutionary biology. Here, we build on previous evolutionary and philosophical treatments of this subject to encompass novelties across biological scales and eco-evolutionary perspectives. By defining novelties as new features at one biological scale that have emergent effects at other biological scales, we incorporate many forms of novelty that have previously been treated in isolation (such as novelty from genetic mutations, new developmental pathways, new morphological features, and new species). Our perspective is based on the fundamental idea that the emergence of a novelty, at any biological scale, depends on its environmental and genetic context. Through this lens, we outline a broad array of generative mechanisms underlying novelty and highlight how genomic tools are transforming our understanding of the origins of novelty. Lastly, we present several case studies to illustrate how novelties across biological scales and systems can be understood based on common mechanisms of change and their environmental and genetic contexts. Specifically, we highlight how gene duplication contributes to the evolution of new complex structures in visual systems; how genetic exchange in symbiosis alters functions of both host and symbiont, resulting in a novel organism; and how hybridisation between species can generate new species with new niches.

Deep Diversity: Extensive Variation in the Components of Complex Visual Systems across Animals

Understanding the molecular underpinnings of the evolution of complex (multi-part) systems is a fundamental topic in biology. One unanswered question is to what the extent do similar or different genes and regulatory interactions underlie similar complex systems across species? Animal eyes and phototransduction (light detection) are outstanding systems to investigate this question because some of the genetics underlying these traits are well characterized in model organisms. However, comparative studies using non-model organisms are also necessary to understand the diversity and evolution of these traits. Here, we compare the characteristics of photoreceptor cells, opsins, and phototransduction cascades in diverse taxa, with a particular focus on cnidarians. In contrast to the common theme of deep homology, whereby similar traits develop mainly using homologous genes, comparisons of visual systems, especially in non-model organisms, are beginning to highlight a "deep diversity" of underlying components, illustrating how variation can underlie similar complex systems across taxa. Although using candidate genes from model organisms across diversity was a good starting point to understand the evolution of complex systems, unbiased genome-wide comparisons and subsequent functional validation will be necessary to uncover unique genes that comprise the complex systems of non-model groups to better understand biodiversity and its evolution.

Nocturnal mosquito Cryptochrome 1 mediates greater electrophysiological and behavioral responses to blue light relative to diurnal mosquito Cryptochrome 1

Nocturnal Anopheles mosquitoes exhibit strong behavioral avoidance to blue-light while diurnal Aedes mosquitoes are behaviorally attracted to blue-light and a wide range of other wavelengths of light. To determine the molecular mechanism of these effects, we expressed light-sensing Anopheles gambiae (AgCRY1) and Aedes aegypti (AeCRY1) Cryptochrome 1 (CRY) genes under a crypGAL4-24 driver line in a mutant Drosophila genetic background lacking native functional CRY, then tested behavioral and electrophysiological effects of mosquito CRY expression relative to positive and negative CRY control conditions. Neither mosquito CRY stops the circadian clock as shown by robust circadian behavioral rhythmicity in constant darkness in flies expressing either AgCRY1 or AeCRY1. AgCRY1 and AeCRY1 both mediate acute increases in large ventral lateral neuronal firing rate evoked by 450 nm blue-light, corresponding to CRY's peak absorbance in its base state, indicating that both mosquito CRYs are functional, however, AgCRY1 mediates significantly stronger sustained electrophysiological light-evoked depolarization in response to blue-light relative to AeCRY1. In contrast, neither AgCRY1 nor AeCRY1 expression mediates measurable increases in large ventral lateral neuronal firing rates in response to 405 nm violet-light, the peak of the Rhodopsin-7 photoreceptor that is co-expressed in the large lateral ventral neurons. These results are consistent with the known action spectra of type 1 CRYs and lack of response in cry-null controls. AgCRY1 and AeCRY1 expressing flies show behavioral attraction to low intensity blue-light, but AgCRY1 expressing flies show behavioral avoidance to higher intensity blue-light. These results show that nocturnal and diurnal mosquito Cryptochrome 1 proteins mediate differential physiological and behavioral responses to blue-light that are consistent with species-specific mosquito behavior.

The genome of the black cutworm Agrotis ipsilon

The black cutworm (BCW), Agrotis ipsilon, is a worldwide polyphagous and underground pest that causes a high level of economic loss to a wide range of crops through the damage of roots. This species performs non-directed migration throughout East and Southeast Asia seasonally. Lack of a genome information has limited further studies on its unique biology and the development of novel management approaches. In this study, we present a 476 Mb de novo assembly of BCW, along with a consensus gene set of 14,801 protein-coding gene models. Quality controls show that both genome assembly and annotations are high-quality and mostly complete. We focus manual annotation and comparative genomics on gene families that related to the unique attributes of this species, such as nocturnality, long-distance migration, and host adaptation. We find that the BCW genome encodes a similar gene repertoire in various migration-related gene families to the diural migratory butterfly Danaus plexiipus, with additional copies of long wavelength opsin and two eye development-related genes. On the other hand, we find that the genomes of BCW and many other polyphagous lepidopterans encode many more gustatory receptor genes, particularly the lineage-specific expanded bitter receptor genes, than the mono- or oligo-phagous species, suggesting a common role of gustatory receptors (GRs) expansion in host range expansion. The availability of a BCW genome provides valuable resources to study the molecular mechanisms of non-directed migration in lepidopteran pests and to develop novel strategies to control migratory nocturnal pests.

Molecular Evolution of Phototransduction Pathway Genes in Nocturnal and Diurnal Fireflies (Coleoptera: Lampyridae)

Most organisms are dependent on sensory cues from their environment for survival and reproduction. Fireflies (Coleoptera: Lampyridae) represent an ideal system for studying sensory niche adaptation due to many species relying on bioluminescent communication; as well as a diversity of ecologies. Here; using transcriptomics; we examine the phototransduction pathway in this non-model organism; and provide some of the first evidence for positive selection in the phototransduction pathway beyond opsins in beetles. Evidence for gene duplications within Lampyridae are found in inactivation no afterpotential C and inactivation no afterpotential D. We also find strong support for positive selection in arrestin-2; inactivation no afterpotential D; and transient receptor potential-like; with weak support for positive selection in guanine nucleotide-binding protein G(q) subunit alpha and neither inactivation nor afterpotential C. Taken with other recent work in flies; butterflies; and moths; this represents an exciting new avenue of study as we seek to further understand diversification and constraint on the phototransduction pathway in light of organism ecology.

CRISPR/Cas9-mediated mutation reveals Pax6 is essential for development of the compound eye in Decapoda Exopalaemon carinicauda

The compound eye in crustaceans is a main eye type in the animal kingdom, knowledge about the mechanism to determine the development of compound eye is very limited. Paired box protein 6 (Pax6) is generally regarded as a master regulator for eye development. In the present study, a genome-based analysis of the Pax6 gene in the ridge tail white prawn Exopalaemon carinicauda was performed and two members of Pax6 homologs, named Ec-Eyeless (EcEy) and Ec-Twin of eyeless (EcToy) were identified. To understand the function of these two homologs of Pax6 gene in the prawn, the CRISPR/Cas9 genome editing technique was applied to generate EcEy and EcToy knock-out (KO) prawns and their phenotypes were analyzed. The surviving EcEy-KO embryos and larvae exhibited severe abnormal eye morphology, suggesting that EcEy is necessary for the compound eye development in prawn, while no mutant phenotype was found in EcToy-KO individuals. These findings highlighted the conservative role of Pax6 gene in the compound eye formation, and the functional differentiation between EcEy and EcToy gene may reveal a novel regulating mechanism of Pax6 on the compound eye development in the decapods. These data will provide important information for understanding the regulation mechanism for crustacean compound eye development.

Molecular evolution and expression of opsin genes in Hydra vulgaris

BACKGROUND

The evolution of opsin genes is of great interest because it can provide insight into the evolution of light detection and vision. An interesting group in which to study opsins is Cnidaria because it is a basal phylum sister to Bilateria with much visual diversity within the phylum. Hydra vulgaris (H. vulgaris) is a cnidarian with a plethora of genomic resources to characterize the opsin gene family. This eyeless cnidarian has a behavioral reaction to light, but it remains unknown which of its many opsins functions in light detection. Here, we used phylogenetics and RNA-seq to investigate the molecular evolution of opsin genes and their expression in H. vulgaris. We explored where opsin genes are located relative to each other in an improved genome assembly and where they belong in a cnidarian opsin phylogenetic tree. In addition, we used RNA-seq data from different tissues of the H. vulgaris adult body and different time points during regeneration and budding stages to gain insight into their potential functions.

RESULTS

We identified 45 opsin genes in H. vulgaris, many of which were located near each other suggesting evolution by tandem duplications. Our phylogenetic tree of cnidarian opsin genes supported previous claims that they are evolving by lineage-specific duplications. We identified two H. vulgaris genes (HvOpA1 and HvOpB1) that fall outside of the two commonly determined Hydra groups; these genes possibly have a function in nematocytes and mucous gland cells respectively. We also found opsin genes that have similar expression patterns to phototransduction genes in H. vulgaris. We propose a H. vulgaris phototransduction cascade that has components of both ciliary and rhabdomeric cascades.

CONCLUSIONS

This extensive study provides an in-depth look at the molecular evolution and expression of H. vulgaris opsin genes. The expression data that we have quantified can be used as a springboard for additional studies looking into the specific function of opsin genes in this species. Our phylogeny and expression data are valuable to investigations of opsin gene evolution and cnidarian biology.

Evolution of Phototransduction Genes in Lepidoptera

Vision is underpinned by phototransduction, a signaling cascade that converts light energy into an electrical signal. Among insects, phototransduction is best understood in Drosophila melanogaster. Comparison of D. melanogaster against three insect species found several phototransduction gene gains and losses, however, lepidopterans were not examined. Diurnal butterflies and nocturnal moths occupy different light environments and have distinct eye morphologies, which might impact the expression of their phototransduction genes. Here we investigated: 1) how phototransduction genes vary in gene gain or loss between D. melanogaster and Lepidoptera, and 2) variations in phototransduction genes between moths and butterflies. To test our prediction of phototransduction differences due to distinct visual ecologies, we used insect reference genomes, phylogenetics, and moth and butterfly head RNA-Seq and transcriptome data. As expected, most phototransduction genes were conserved between D. melanogaster and Lepidoptera, with some exceptions. Notably, we found two lepidopteran opsins lacking a D. melanogaster ortholog. Using antibodies we found that one of these opsins, a candidate retinochrome, which we refer to as unclassified opsin (UnRh), is expressed in the crystalline cone cells and the pigment cells of the butterfly, Heliconius melpomene. Our results also show that butterflies express similar amounts of trp and trpl channel mRNAs, whereas moths express ∼50× less trp, a potential adaptation to darkness. Our findings suggest that while many single-copy D. melanogaster phototransduction genes are conserved in lepidopterans, phototransduction gene expression differences exist between moths and butterflies that may be linked to their visual light environment.

Pluripotency and the origin of animal multicellularity

A widely held-but rarely tested-hypothesis for the origin of animals is that they evolved from a unicellular ancestor, with an apical cilium surrounded by a microvillar collar, that structurally resembled modern sponge choanocytes and choanoflagellates^1-4. Here we test this view of animal origins by comparing the transcriptomes, fates and behaviours of the three primary sponge cell types-choanocytes, pluripotent mesenchymal archaeocytes and epithelial pinacocytes-with choanoflagellates and other unicellular holozoans. Unexpectedly, we find that the transcriptome of sponge choanocytes is the least similar to the transcriptomes of choanoflagellates and is significantly enriched in genes unique to either animals or sponges alone. By contrast, pluripotent archaeocytes upregulate genes that control cell proliferation and gene expression, as in other metazoan stem cells and in the proliferating stages of two unicellular holozoans, including a colonial choanoflagellate. Choanocytes in the sponge Amphimedon queenslandica exist in a transient metastable state and readily transdifferentiate into archaeocytes, which can differentiate into a range of other cell types. These sponge cell-type conversions are similar to the temporal cell-state changes that occur in unicellular holozoans⁵. Together, these analyses argue against homology of sponge choanocytes and choanoflagellates, and the view that the first multicellular animals were simple balls of cells with limited capacity to differentiate. Instead, our results are consistent with the first animal cell being able to transition between multiple states in a manner similar to modern transdifferentiating and stem cells.

A likelihood ratio test for changes in homeolog expression bias

BACKGROUND

Gene duplications are a major source of raw material for evolution and a likely contributor to the diversity of life on earth. Duplicate genes (i.e., homeologs, in the case of a whole genome duplication) may retain their ancestral function, sub- or neofunctionalize, or be lost entirely. A primary way that duplicate genes evolve new functions is by altering their expression patterns. Comparing the expression patterns of duplicate genes gives clues as to whether any of these evolutionary processes have occurred.

RESULTS

We develop a likelihood ratio test for the analysis of the expression ratios of duplicate genes across two conditions (e.g., tissues). We demonstrate an application of this test by comparing homeolog expression patterns of 1448 homeologous gene pairs using RNA-seq data generated from leaves and petals of an allotetraploid monkeyflower (Mimulus luteus). We assess the sensitivity of this test to different levels of homeolog expression bias and compare the method to several alternatives.

CONCLUSIONS

The likelihood ratio test derived here is a direct, transparent, and easily implemented method for detecting changes in homeolog expression bias that outperforms alternative approaches. While our method was derived with homeolog analysis in mind, this method can be used to analyze changes in the ratio of expression levels between any two genes in any two conditions.

Phototransduction Gene Expression and Evolution in Cave and Surface Crayfishes

In the absence of light in caves, animals have repeatedly evolved reduced eyes and visual systems. Whether the underlying genetic components remain intact in blind species remains unanswered across taxa. The freshwater crayfish have evolved to live in caves multiple times throughout their history; therefore, this system provides an opportunity to probe the genetic patterns and processes underlying repeated vision loss. Using transcriptomic data from the eyes of 14 species of cave and surface crayfishes, we identify the expression of 17 genes putatively related to visual phototransduction. We find a similarly complete repertoire of phototransduction gene families expressed in cave and surface species, but that the expression levels of those transcripts are consistently lower in cave species. We find statistical support for episodic positive selection, increased and decreased selection strength in caves, depending on the gene family. Analyses of gene expression evolution suggest convergent and possibly adaptive downregulation of these genes across eye-reduction events. Our results reveal a combination of evolutionary processes acting on the sequences and gene expression levels of vision-related genes underlying the loss of vision in caves.

Rapid diversification of homothorax expression patterns after gene duplication in spiders

Background

Gene duplications provide genetic material for the evolution of new morphological and physiological features. One copy can preserve the original gene functions while the second copy may evolve new functions (neofunctionalisation). Gene duplications may thus provide new genes involved in evolutionary novelties.

Results

We have studied the duplicated homeobox gene homothorax (hth) in the spider species Parasteatoda tepidariorum and Pholcus phalangioides and have compared these data with previously published data from additional spider species. We show that the expression pattern of hth1 is highly conserved among spiders, consistent with the notion that this gene copy preserves the original hth functions. By contrast, hth2 has a markedly different expression profile especially in the prosomal appendages. The pattern in the pedipalps and legs consists of several segmental rings, suggesting a possible role of hth2 in limb joint development. Intriguingly, however, the hth2 pattern is much less conserved between the species than hth1 and shows a species specific pattern in each species investigated so far.

Conclusions

We hypothesise that the hth2 gene has gained a new patterning function after gene duplication, but has then undergone a second phase of diversification of its new role in the spider clade. The evolution of hth2 may thus provide an interesting example for a duplicated gene that has not only contributed to genetic diversity through neofunctionalisation, but beyond that has been able to escape evolutionary conservation after neofunctionalisation thus forming the basis for further genetic diversification.

Electronic supplementary material

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

First Insights into the Subterranean Crustacean Bathynellacea Transcriptome: Transcriptionally Reduced Opsin Repertoire and Evidence of Conserved Homeostasis Regulatory Mechanisms

Bathynellacea (Crustacea, Syncarida, Parabathynellidae) are subterranean aquatic crustaceans that typically inhabit freshwater interstitial spaces (e.g., groundwater) and are occasionally found in caves and even hot springs. In this study, we sequenced the whole transcriptome of Allobathynella bangokensis using RNA-seq. De novo sequence assembly produced 74,866 contigs including 28,934 BLAST hits. Overall, the gene sequences were most similar to those of the waterflea Daphnia pulex. In the A. bangokensis transcriptome, no opsin or related sequences were identified, and no contig aligned to the crustacean visual opsins and non-visual opsins (i.e. arthropsins, peropsins, and melaopsins), suggesting potential regressive adaptation to the dark environment. However, A. bangokensis expressed conserved gene family sets, such as heat shock proteins and those related to key innate immunity pathways and antioxidant defense systems, at the transcriptional level, suggesting that this species has evolved adaptations involving molecular mechanisms of homeostasis. The transcriptomic information of A. bangokensis will be useful for investigating molecular adaptations and response mechanisms to subterranean environmental conditions.

Structural differences and differential expression among rhabdomeric opsins reveal functional change after gene duplication in the bay scallop, Argopecten irradians (Pectinidae)

Background

Opsins are the only class of proteins used for light perception in image-forming eyes. Gene duplication and subsequent functional divergence of opsins have played an important role in expanding photoreceptive capabilities of organisms by altering what wavelengths of light are absorbed by photoreceptors (spectral tuning). However, new opsin copies may also acquire novel function or subdivide ancestral functions through changes to temporal, spatial or the level of gene expression. Here, we test how opsin gene copies diversify in function and evolutionary fate by characterizing four rhabdomeric (G_q-protein coupled) opsins in the scallop, Argopecten irradians, identified from tissue-specific transcriptomes.

Results

Under a phylogenetic analysis, we recovered a pattern consistent with two rounds of duplication that generated the genetic diversity of scallop G_q-opsins. We found strong support for differential expression of paralogous G_q-opsins across ocular and extra-ocular photosensitive tissues, suggesting that scallop G_q-opsins are used in different biological contexts due to molecular alternations outside and within the protein-coding regions. Finally, we used available protein models to predict which amino acid residues interact with the light-absorbing chromophore. Variation in these residues suggests that the four G_q-opsin paralogs absorb different wavelengths of light.

Conclusions

Our results uncover novel genetic and functional diversity in the light-sensing structures of the scallop, demonstrating the complicated nature of G_q-opsin diversification after gene duplication. Our results highlight a change in the nearly ubiquitous shadow response in molluscs to a narrowed functional specificity for visual processes in the eyed scallop. Our findings provide a starting point to study how gene duplication may coincide with eye evolution, and more specifically, different ways neofunctionalization of G_q-opsins may occur.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0823-9) contains supplementary material, which is available to authorized users.

A second visual rhodopsin gene, rh1-2, is expressed in zebrafish photoreceptors and found in other ray-finned fishes

Rhodopsin (rh1) is the visual pigment expressed in rod photoreceptors of vertebrates that is responsible for initiating the critical first step of dim-light vision. Rhodopsin is usually a single copy gene; however, we previously discovered a novel rhodopsin-like gene expressed in the zebrafish retina, rh1-2, which we identified as a functional photosensitive pigment that binds 11-cis retinal and activates in response to light. Here, we localized expression of rh1-2 in the zebrafish retina to a subset of peripheral photoreceptor cells, which indicates a partially overlapping expression pattern with rh1 We also expressed, purified and characterized Rh1-2, including investigation of the stability of the biologically active intermediate. Using fluorescence spectroscopy, we found the half-life of the rate of retinal release of Rh1-2 following photoactivation to be more similar to that of the visual pigment rhodopsin than to the non-visual pigment exo-rhodopsin (exorh), which releases retinal around 5 times faster. Phylogenetic and molecular evolutionary analyses show that rh1-2 has ancient origins within teleost fishes, is under similar selective pressure to rh1, and likely experienced a burst of positive selection following its duplication and divergence from rh1 These findings indicate that rh1-2 is another functional visual rhodopsin gene, which contradicts the prevailing notion that visual rhodopsin is primarily found as a single copy gene within ray-finned fishes. The reasons for retention of this duplicate gene, as well as possible functional consequences for the visual system, are discussed.

Brain evolution by brain pathway duplication

Understanding the mechanisms of evolution of brain pathways for complex behaviours is still in its infancy. Making further advances requires a deeper understanding of brain homologies, novelties and analogies. It also requires an understanding of how adaptive genetic modifications lead to restructuring of the brain. Recent advances in genomic and molecular biology techniques applied to brain research have provided exciting insights into how complex behaviours are shaped by selection of novel brain pathways and functions of the nervous system. Here, we review and further develop some insights to a new hypothesis on one mechanism that may contribute to nervous system evolution, in particular by brain pathway duplication. Like gene duplication, we propose that whole brain pathways can duplicate and the duplicated pathway diverge to take on new functions. We suggest that one mechanism of brain pathway duplication could be through gene duplication, although other mechanisms are possible. We focus on brain pathways for vocal learning and spoken language in song-learning birds and humans as example systems. This view presents a new framework for future research in our understanding of brain evolution and novel behavioural traits.

Genome-wide profiling of 24 hr diel rhythmicity in the water flea, Daphnia pulex: network analysis reveals rhythmic gene expression and enhances functional gene annotation

Background

Marine and freshwater zooplankton exhibit daily rhythmic patterns of behavior and physiology which may be regulated directly by the light:dark (LD) cycle and/or a molecular circadian clock. One of the best-studied zooplankton taxa, the freshwater crustacean Daphnia, has a 24 h diel vertical migration (DVM) behavior whereby the organism travels up and down through the water column daily. DVM plays a critical role in resource tracking and the behavioral avoidance of predators and damaging ultraviolet radiation. However, there is little information at the transcriptional level linking the expression patterns of genes to the rhythmic physiology/behavior of Daphnia.

Results

Here we analyzed genome-wide temporal transcriptional patterns from Daphnia pulex collected over a 44 h time period under a 12:12 LD cycle (diel) conditions using a cosine-fitting algorithm. We used a comprehensive network modeling and analysis approach to identify novel co-regulated rhythmic genes that have similar network topological properties and functional annotations as rhythmic genes identified by the cosine-fitting analyses. Furthermore, we used the network approach to predict with high accuracy novel gene-function associations, thus enhancing current functional annotations available for genes in this ecologically relevant model species. Our results reveal that genes in many functional groupings exhibit 24 h rhythms in their expression patterns under diel conditions. We highlight the rhythmic expression of immunity, oxidative detoxification, and sensory process genes. We discuss differences in the chronobiology of D. pulex from other well-characterized terrestrial arthropods.

Conclusions

This research adds to a growing body of literature suggesting the genetic mechanisms governing rhythmicity in crustaceans may be divergent from other arthropod lineages including insects. Lastly, these results highlight the power of using a network analysis approach to identify differential gene expression and provide novel functional annotation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2998-2) contains supplementary material, which is available to authorized users.

Towards Consensus Gene Ages

Correctly estimating the age of a gene or gene family is important for a variety of fields, including molecular evolution, comparative genomics, and phylogenetics, and increasingly for systems biology and disease genetics. However, most studies use only a point estimate of a gene’s age, neglecting the substantial uncertainty involved in this estimation. Here, we characterize this uncertainty by investigating the effect of algorithm choice on gene-age inference and calculate consensus gene ages with attendant error distributions for a variety of model eukaryotes. We use 13 orthology inference algorithms to create gene-age datasets and then characterize the error around each age-call on a per-gene and per-algorithm basis. Systematic error was found to be a large factor in estimating gene age, suggesting that simple consensus algorithms are not enough to give a reliable point estimate. We also found that different sources of error can affect downstream analyses, such as gene ontology enrichment. Our consensus gene-age datasets, with associated error terms, are made fully available at so that researchers can propagate this uncertainty through their analyses (geneages.org).

Two temporal functions of Glass: Ommatidium patterning and photoreceptor differentiation

Much progress has been made in elucidating the molecular networks required for specifying retinal cells, including photoreceptors, but the downstream mechanisms that maintain identity and regulate differentiation remain poorly understood. Here, we report that the transcription factor Glass has a dual role in establishing a functional Drosophila eye. Utilizing conditional rescue approaches, we confirm that persistent defects in ommatidium patterning combined with cell death correlate with the overall disruption of eye morphology in glass mutants. In addition, we reveal that Glass exhibits a separable role in regulating photoreceptor differentiation. In particular, we demonstrate the apparent loss of glass mutant photoreceptors is not only due to cell death but also a failure of the surviving photoreceptors to complete differentiation. Moreover, the late reintroduction of Glass in these developmentally stalled photoreceptors is capable of restoring differentiation in the absence of correct ommatidium patterning. Mechanistically, transcription profiling at the time of differentiation reveals that Glass is necessary for the expression of many genes implicated in differentiation, i.e. rhabdomere morphogenesis, phototransduction, and synaptogenesis. Specifically, we show Glass directly regulates the expression of Pph13, which encodes a transcription factor necessary for opsin expression and rhabdomere morphogenesis. Finally, we demonstrate the ability of Glass to choreograph photoreceptor differentiation is conserved between Drosophila and Tribolium, two holometabolous insects. Altogether, our work identifies a fundamental regulatory mechanism to generate the full complement of cells required for a functional rhabdomeric visual system and provides a critical framework to investigate the basis of differentiation and maintenance of photoreceptor identity.

Sexually dimorphic gene expression in the lateral eyes of Euphilomedes carcharodonta (Ostracoda, Pancrustacea)

Background

The evolution and development of sexual dimorphism illuminates a central question in biology: How do similar genomes produce different phenotypes? In an XX/XO system especially the state of a sexually dimorphic trait is determined by differences in gene expression, as there are no additional genetic loci in either sex. Here, we examine the XX/XO ostracod crustacean species Euphilomedes carcharodonta. This species exhibits radical sexual dimorphism of their lateral eyes, females have only a tiny simple lateral eye while males have elaborate ommatidial eyes.

Results

We find that males express three of nine eye-development gene homologs at significantly higher levels during juvenile eye development, compared to females. We also find that most eye-development genes examined are pleiotropic, with high expression levels during embryonic development as well as during juvenile eye development. Later, in adults, we find that phototransduction genes are expressed at higher levels in males than in females, as we might expect when comparing ommatidial to simple eyes.

Conclusions

We show here that expression changes of a handful of developmental genes may underlie the radical difference in a dimorphic character. This work gives an important point of comparison for studying eye evolution and development in the Pancrustacea.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0026-2) contains supplementary material, which is available to authorized users.

Transcriptome-Wide Differential Gene Expression in Bicyclus anynana Butterflies: Female Vision-Related Genes Are More Plastic

Vision is energetically costly to maintain. Consequently, over time many cave-adapted species downregulate the expression of vision genes or even lose their eyes and associated eye genes entirely. Alternatively, organisms that live in fluctuating environments, with different requirements for vision at different times, may evolve phenotypic plasticity for expression of vision genes. Here, we use a global transcriptomic and candidate gene approach to compare gene expression in the heads of a polyphenic butterfly. Bicyclus anynana have two seasonal forms that display sexual dimorphism and plasticity in eye morphology, and female-specific plasticity in opsin gene expression. Nonchoosy dry season females downregulate opsin expression, consistent with the high physiological cost of vision. To identify other genes associated with sexually dimorphic and seasonally plastic differences in vision, we analyzed RNA-sequencing data from whole head tissues. We identified two eye development genes (klarsicht and warts homologs) and an eye pigment biosynthesis gene (henna) differentially expressed between seasonal forms. By comparing sex-specific expression across seasonal forms, we found that klarsicht, warts, henna, and another eye development gene (domeless) were plastic in a female-specific manner. In a male-only analysis, white (w) was differentially expressed between seasonal forms. Reverse transcription polymerase chain reaction confirmed that warts and white are expressed in eyes only, whereas klarsicht, henna and domeless are expressed in both eyes and brain. We find that differential expression of eye development and eye pigment genes is associated with divergent eye phenotypes in B. anynana seasonal forms, and that there is a larger effect of season on female vision-related genes.

Applications of comparative evolution to human disease genetics

Direct comparison of human diseases with model phenotypes allows exploration of key areas of human biology which are often inaccessible for practical or ethical reasons. We review recent developments in comparative evolutionary approaches for finding models for genetic disease, including high-throughput generation of gene/phenotype relationship data, the linking of orthologous genes and phenotypes across species, and statistical methods for linking human diseases to model phenotypes.

Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative “parasitism genes.” Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.

Using phylogenetically-informed annotation (PIA) to search for light-interacting genes in transcriptomes from non-model organisms

Background

Tools for high throughput sequencing and de novo assembly make the analysis of transcriptomes (i.e. the suite of genes expressed in a tissue) feasible for almost any organism. Yet a challenge for biologists is that it can be difficult to assign identities to gene sequences, especially from non-model organisms. Phylogenetic analyses are one useful method for assigning identities to these sequences, but such methods tend to be time-consuming because of the need to re-calculate trees for every gene of interest and each time a new data set is analyzed. In response, we employed existing tools for phylogenetic analysis to produce a computationally efficient, tree-based approach for annotating transcriptomes or new genomes that we term Phylogenetically-Informed Annotation (PIA), which places uncharacterized genes into pre-calculated phylogenies of gene families.

Results

We generated maximum likelihood trees for 109 genes from a Light Interaction Toolkit (LIT), a collection of genes that underlie the function or development of light-interacting structures in metazoans. To do so, we searched protein sequences predicted from 29 fully-sequenced genomes and built trees using tools for phylogenetic analysis in the Osiris package of Galaxy (an open-source workflow management system). Next, to rapidly annotate transcriptomes from organisms that lack sequenced genomes, we repurposed a maximum likelihood-based Evolutionary Placement Algorithm (implemented in RAxML) to place sequences of potential LIT genes on to our pre-calculated gene trees. Finally, we implemented PIA in Galaxy and used it to search for LIT genes in 28 newly-sequenced transcriptomes from the light-interacting tissues of a range of cephalopod mollusks, arthropods, and cubozoan cnidarians. Our new trees for LIT genes are available on the Bitbucket public repository (http://bitbucket.org/osiris_phylogenetics/pia/) and we demonstrate PIA on a publicly-accessible web server (http://galaxy-dev.cnsi.ucsb.edu/pia/).

Conclusions

Our new trees for LIT genes will be a valuable resource for researchers studying the evolution of eyes or other light-interacting structures. We also introduce PIA, a high throughput method for using phylogenetic relationships to identify LIT genes in transcriptomes from non-model organisms. With simple modifications, our methods may be used to search for different sets of genes or to annotate data sets from taxa outside of Metazoa.

Electronic supplementary material

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

UniRef clusters: a comprehensive and scalable alternative for improving sequence similarity searches

MOTIVATION

UniRef databases provide full-scale clustering of UniProtKB sequences and are utilized for a broad range of applications, particularly similarity-based functional annotation. Non-redundancy and intra-cluster homogeneity in UniRef were recently improved by adding a sequence length overlap threshold. Our hypothesis is that these improvements would enhance the speed and sensitivity of similarity searches and improve the consistency of annotation within clusters.

RESULTS

Intra-cluster molecular function consistency was examined by analysis of Gene Ontology terms. Results show that UniRef clusters bring together proteins of identical molecular function in more than 97% of the clusters, implying that clusters are useful for annotation and can also be used to detect annotation inconsistencies. To examine coverage in similarity results, BLASTP searches against UniRef50 followed by expansion of the hit lists with cluster members demonstrated advantages compared with searches against UniProtKB sequences; the searches are concise (∼7 times shorter hit list before expansion), faster (∼6 times) and more sensitive in detection of remote similarities (>96% recall at e-value <0.0001). Our results support the use of UniRef clusters as a comprehensive and scalable alternative to native sequence databases for similarity searches and reinforces its reliability for use in functional annotation.

Distinct evolutionary rate in the eye field transcription factors found by estimation of ancestral protein structure

Eye-field transcription factors (EFTFs) are a set of genes that compose a regulatory network for eye development in animals, which are highly conserved among various animal phyla. To investigate the processes of conservation and diversification of the transcription factors for eye development, we examined the structural changes in the EFTF proteins by estimating the ancestral sequences with the available genome information. Among the different types of EFTFs, we selected otx2, tbx3, rx1, pax6, six3/6, lhx2 and nr2e1 because they are highly conserved in bilaterian animals. We searched the genome sequences of representative animal phyla for EFTF protein sequences. With deduced ancestral sequences and three-dimensional structures of EFTFs, we traced the evolutionary changes in amino acid residues and found that the DNA-binding domains were always more conserved than other regions, and that the other regions showed distinct evolutionary rates. The EFTF rx1, which resides at the pivotal part of the EFTF network, had a faster evolutionary rate than the others. These results indicated that the evolutionary rates of each protein in the EFTF network, which were expected to be consistent with each other to maintain the interactions in the network, were not constant among or within the factors, but rather, varied to a significant extent.

Common transcriptional mechanisms for visual photoreceptor cell differentiation among Pancrustaceans

A hallmark of visual rhabdomeric photoreceptors is the expression of a rhabdomeric opsin and uniquely associated phototransduction molecules, which are incorporated into a specialized expanded apical membrane, the rhabdomere. Given the extensive utilization of rhabdomeric photoreceptors in the eyes of protostomes, here we address whether a common transcriptional mechanism exists for the differentiation of rhabdomeric photoreceptors. In Drosophila, the transcription factors Pph13 and Orthodenticle (Otd) direct both aspects of differentiation: rhabdomeric opsin transcription and rhabdomere morphogenesis. We demonstrate that the orthologs of both proteins are expressed in the visual systems of the distantly related arthropod species Tribolium castaneum and Daphnia magna and that their functional roles are similar in these species. In particular, we establish that the Pph13 homologs have the ability to bind a subset of Rhodopsin core sequence I sites and that these sites are present in key phototransduction genes of both Tribolium and Daphnia. Furthermore, Pph13 and Otd orthologs are capable of executing deeply conserved functions of photoreceptor differentiation as evidenced by the ability to rescue their respective Drosophila mutant phenotypes. Pph13 homologs are equivalent in their ability to direct both rhabdomere morphogenesis and opsin expression within Drosophila, whereas Otd paralogs demonstrate differential abilities to regulate photoreceptor differentiation. Finally, loss-of-function analyses in Tribolium confirm the conserved requirement of Pph13 and Otd in regulating both rhabdomeric opsin transcription and rhabdomere morphogenesis. Taken together, our data identify components of a regulatory framework for rhabdomeric photoreceptor differentiation in Pancrustaceans, providing a foundation for defining ancestral regulatory modules of rhabdomeric photoreceptor differentiation.

Visual systems are populated by one of two fundamental types of photoreceptors, ciliary and rhabdomeric. Each photoreceptor type is defined by the opsin molecule expressed and the final morphological form adapted to house the phototransduction machinery. Here we address whether a common transcriptional mechanisms exists for the differentiation of rhabdomeric photoreceptors. We demonstrate that orthologs of two Drosophila (fruit fly) transcription factors, Pph13 and Orthodenticle, are expressed in photoreceptors of Pancrustaceans, Tribolium (red flour beetle) and Daphnia (water flea), and are capable of executing conserved functions of rhabdomeric photoreceptor differentiation. In particular, Tribolium and Daphnia orthologs are capable of substituting and rescuing the photoreceptor differentiation defects observed in their corresponding Drosophila mutants. Furthermore, loss of function analysis in Tribolium of both Pph13 and orthodenticle genes demonstrate they regulate opsin transcription and morphogenesis of the photoreceptor apical membrane. Our data illuminate a framework for rhabdomeric photoreceptor differentiation and provide the foundation for defining the ancestral regulatory modules for rhabdomeric differentiation and potential modifications that underlie the functional diversity observed in rhabdomeric photoreceptors.

Cephalopod eye evolution was modulated by the acquisition of Pax-6 splicing variants

Previous studies have reported that the developmental processes of vertebrate eyes are controlled by four Pax-6 splicing variants, each modulating different downstream genes, whereas those of insect eyes are controlled by duplicated Pax-6 genes. Cephalopods belong to the Protostomes but possess a camera-type eye similar to those in vertebrates. We examined Pax-6 variations in the squid and found five types of Pax-6 splicing variants but no duplication of the Pax-6 gene. In the five splicing variants, the splicing patterns were produced by the combination of two additional exons to the ortholog and one jettisoned exon containing most of the Homeobox domain (HD). These five variants show spatio-temporal patterns of gene expression during development in the squid. Our study suggests that cephalopods acquired Pax-6 splicing variants independent of those in vertebrates and that these variants were similarly utilized in the development of the squid eye.

The evolution of complexity in the visual systems of stomatopods: insights from transcriptomics

Stomatopod crustaceans have complex visual systems containing up to 16 different spectral classes of photoreceptors, more than described for any other animal. A previous molecular study of this visual system focusing on the expression of opsin genes found many more transcripts than predicted on the basis of physiology, but was unable to fully document the expressed opsin genes responsible for this diversity. Furthermore, questions remain about how other components of phototransduction cascades are involved. This study continues prior investigations by examining the molecular function of stomatopods' visual systems using new whole eye 454 transcriptome datasets from two species, Hemisquilla californiensis and Pseudosquilla ciliata. These two species represent taxonomic diversity within the order Stomatopoda, as well as variations in the anatomy and physiology of the visual system. Using an evolutionary placement algorithm to annotate the transcriptome, we identified the presence of nine components of the stomatopods' G-protein-coupled receptor (GPCR) phototransduction cascade, including two visual arrestins, subunits of the heterotrimeric G-protein, phospholipase C, transient receptor potential channels, and opsin transcripts. The set of expressed transduction genes suggests that stomatopods utilize a Gq-mediated GPCR-signaling cascade. The most notable difference in expression between the phototransduction cascades of the two species was the number of opsin contigs recovered, with 18 contigs found in retinas of H. californiensis, and 49 contigs in those of P. ciliata. Based on phylogenetic placement and fragment overlap, these contigs were estimated to represent 14 and 33 expressed transcripts, respectively. These data expand the known opsin diversity in stomatopods to clades of arthropod opsins that are sensitive to short wavelengths and ultraviolet wavelengths and confirm the results of previous studies recovering more opsin transcripts than spectrally distinct types of photoreceptors. Many of the recovered transcripts were phylogenetically placed in an evolutionary clade of crustacean opsin sequences that is rapidly expanding as the visual systems from more species are investigated. We discuss these results in relation to the emerging pattern, particularly in crustacean visual systems, of the expression of multiple opsin transcripts in photoreceptors of the same spectral class, and even in single photoreceptor cells.

Blue-light-receptive cryptochrome is expressed in a sponge eye lacking neurons and opsin

Many larval sponges possess pigment ring eyes that apparently mediate phototactic swimming. Yet sponges are not known to possess nervous systems or opsin genes, so the unknown molecular components of sponge phototaxis must differ fundamentally from those in other animals, inspiring questions about how this sensory system functions. Here we present molecular and biochemical data on cryptochrome, a candidate gene for functional involvement in sponge pigment ring eyes. We report that Amphimedon queenslandica, a demosponge, possesses two cryptochrome/photolyase genes, Aq-Cry1 and Aq-Cry2. The mRNA of one gene (Aq-Cry2) is expressed in situ at the pigment ring eye. Additionally, we report that Aq-Cry2 lacks photolyase activity and contains a flavin-based co-factor that is responsive to wavelengths of light that also mediate larval photic behavior. These results suggest that Aq-Cry2 may act in the aneural, opsin-less phototaxic behavior of a sponge.

General trends in the utilization of structural factors contributing to biological complexity

During evolution, proteins containing newly emerged domains and the increasing proportion of multidomain proteins in the full genome-encoded proteome (GEP) have substantially contributed to increasing biological complexity. However, it is not known how these two potential structural factors are preferentially utilized at given physiological states. Here, we classified proteins according to domain number and domain age and explored the general trends across species for the utilization of proteins from GEP to various certain-state proteomes (CSPs, i.e., all the proteins expressed at certain physiological states). We found that multidomain proteins or only older domain-containing proteins are significantly overrepresented in CSPs compared with GEP, which is a trend that is stronger in multicellular organisms than in unicellular organisms. Interestingly, the strengths of overrepresentation decreased during evolution of multicellular eukaryotes. When comparing across CSPs, we found that multidomain proteins are more overrepresented in complex tissues than in simpler ones, whereas no difference among proteins with domains of different ages is evident between complex and simple tissues. Thus, biological complexity under certain conditions is more significantly realized by diverse domain organization than by the emergence of new types of domain. In addition, we found that multidomain or only older domain-containing proteins tend to evolve slowly and generally are under stronger purifying selection, which may partly result from their general overrepresentation trends in CSPs.