Ancestral Networks Research Articles

Tracing the genetic footprints of vertebrate landing in non-teleost ray-finned fishes

Rich fossil evidence suggests that many traits and functions related to terrestrial evolution were present long before the ancestor of lobe- and ray-finned fishes. Here, we present genome sequences of the bichir, paddlefish, bowfin, and alligator gar, covering all major early divergent lineages of ray-finned fishes. Our analyses show that these species exhibit many mosaic genomic features of lobe- and ray-finned fishes. In particular, many regulatory elements for limb development are present in these fishes, supporting the hypothesis that the relevant ancestral regulation networks emerged before the origin of tetrapods. Transcriptome analyses confirm the homology between the lung and swim bladder and reveal the presence of functional lung-related genes in early ray-finned fishes. Furthermore, we functionally validate the essential role of a jawed vertebrate highly conserved element for cardiovascular development. Our results imply the ancestors of jawed vertebrates already had the potential gene networks for cardio-respiratory systems supporting air breathing.

Maximum likelihood reconstruction of ancestral networks by integer linear programming.

The study of the evolutionary history of biological networks enables deep functional understanding of various bio-molecular processes. Network growth models, such as the Duplication-Mutation with Complementarity (DMC) model, provide a principled approach to characterizing the evolution of protein-protein interactions (PPIs) based on duplication and divergence. Current methods for model-based ancestral network reconstruction primarily use greedy heuristics and yield sub-optimal solutions. We present a new Integer Linear Programming (ILP) solution for maximum likelihood reconstruction of ancestral PPI networks using the DMC model. We prove the correctness of our solution that is designed to find the optimal solution. It can also use efficient heuristics from general-purpose ILP solvers to obtain multiple optimal and near-optimal solutions that may be useful in many applications. Experiments on synthetic data show that our ILP obtains solutions with higher likelihood than those from previous methods, and is robust to noise and model mismatch. We evaluate our algorithm on two real PPI networks, with proteins from the families of bZIP transcription factors and the Commander complex. On both the networks, solutions from our ILP have higher likelihood and are in better agreement with independent biological evidence from other studies. A Python implementation is available at https://bitbucket.org/cdal/network-reconstruction. Supplementary data are available at Bioinformatics online.

Screening and evaluation of Terminalia catappa and its compound application

Plant determined medications remains an important source, particularly in creating countries, to look at genuine sicknesses roughly 62-80% of the total populace although everything depends on conventional medication for the therapy of specific disease. Indeed, plants produce a various scope of bioactive atoms creation them an ironic wellspring of an alternate kind of drugs. There are hardly any reports and utilization of plants in conventional mending by either ancestral individuals or indigenous network. Rejuvenating plants are the wellspring of extraordinary monetary estimation of everywhere on over the world. Nature has the best word on us a rich plant riches, and an enormous number of assorted kinds of plants develop in various pieces of the nation. Homegrown medication is as yet a pillar of around 75 to 85% the entire populace and the significant aspect of the conventional therapy the utilization of the plant extricate and the dynamic constituents. Among the 7000 types of rejuvenating plants perceived everywhere on over the world in excess of 9000, valuable medicinal plants are found in India. Unfortunately, just not many of them are utilized for their therapeutic worth. Around 1500 plants systematically utilize the conventional arrangement of Indian medication. Notwithstanding, the ethanopharmacologist, microbiologist, botanist and common item physicist world over today, is continually still looking for therapeutic adequacy of the plants on the phytochemicals. Along these lines, the quest for the new phytochemical is the foremost significant important to research the primer phytochemical examination to Terminalia catappa and Syzygium jumbolanam.

Enhancer evolution in chordates: Lessons from functional analyses of cephalochordate cis-regulatory modules.

Chordates comprise three major groups, cephalochordates (amphioxus), tunicates (urochordates), and vertebrates. Since cephalochordates were the early branching group, comparisons between amphioxus and other chordates help us to speculate about ancestral chordates. Here, I summarize accumulating data from functional studies analyzing amphioxus cis-regulatory modules (CRMs) in model systems of other chordate groups, such as mice, chickens, clawed frogs, fish, and ascidians. Conservatism and variability of CRM functions illustrate how gene regulatory networks have evolved in chordates. Amphioxus CRMs, which correspond to CRMs deeply conserved among animal phyla, govern reporter gene expression in conserved expression domains of the putative target gene in host animals. In addition, some CRMs located in similar genomic regions (intron, upstream, or downstream) also possess conserved activity, even though their sequences are divergent. These conservative CRM functions imply ancestral genomic structures and gene regulatory networks in chordates. However, interestingly, if expression patterns of amphioxus genes do not correspond to those of orthologs of experimental models, some amphioxus CRMs recapitulate expression patterns of amphioxus genes, but not those of endogenous genes, suggesting that these amphioxus CRMs are close to the ancestral states of chordate CRMs, while vertebrates/tunicates innovated new CRMs to reconstruct gene regulatory networks subsequent to the divergence of the cephalochordates. Alternatively, amphioxus CRMs may have secondarily lost ancestral CRM activity and evolved independently. These data help to solve fundamental questions of chordate evolution, such as neural crest cells, placodes, a forebrain/midbrain, and genome duplication. Experimental validation is crucial to verify CRM functions and evolution.

Dissecting Adaptive Traits with Nested Association Mapping: Genetic Architecture of Inflorescence Morphology in Sorghum.

In the cereal crop sorghum (Sorghum bicolor) inflorescence morphology variation underlies yield variation and confers adaptation across precipitation gradients, but its genetic basis is poorly understood. We characterized the genetic architecture of sorghum inflorescence morphology using a global nested association mapping (NAM) population (2200 recombinant inbred lines) and 198,000 phenotypic observations from multi-environment trials for four inflorescence morphology traits (upper branch length, lower branch length, rachis length, and rachis diameter). Trait correlations suggest that lower and upper branch length are under somewhat independent control, while lower branch length and rachis diameter are highly pleiotropic. Joint linkage and genome-wide association mapping revealed an oligogenic architecture with 1–22 QTL per trait, each explaining 0.1–5.0% of variation across the entire NAM population. There is a significant enrichment (2.twofold) of QTL colocalizing with grass inflorescence gene homologs, notably with orthologs of maize Ramosa2 and rice Aberrant Panicle Organization1 and TAWAWA1. Still, many QTL do not colocalize with inflorescence gene homologs. In global georeferenced germplasm, allelic variation at the major inflorescence QTL is geographically patterned but only weakly associated with the gradient of annual precipitation. Comparison of NAM with diversity panel association suggests that naive association models may capture some true associations not identified by mixed linear models. Overall, the findings suggest that global inflorescence diversity in sorghum is largely controlled by oligogenic, epistatic, and pleiotropic variation in ancestral regulatory networks. The findings also provide a basis for genomics-enabled breeding of locally-adapted inflorescence morphology.

MAPPS: A Web-Based Tool for Metabolic Pathway Prediction and Network Analysis in the Postgenomic Era.

Comparative and evolutionary analyses of metabolic networks have a wide range of applications, ranging from research into metabolic evolution through to practical applications in drug development, synthetic biology, and biodegradation. We present MAPPS: Metabolic network Analysis and Pathway Prediction Server (https://mapps.lums.edu.pk), a web-based tool to study functions and evolution of metabolic networks using traditional and 'omics data sets. MAPPS provides diverse functionalities including an interactive interface, graphical visualization of results, pathway prediction and network comparison, identification of potential drug targets, in silico metabolic engineering, host-microbe interactions, and ancestral network building. Importantly, MAPPS also allows users to upload custom data, thus enabling metabolic analyses on draft and custom genomes, and has an 'omics pipeline to filter pathway results, making it relevant in today's postgenomic era.

The evolution of the gene regulatory networks patterning the Drosophila Blastoderm.

The Drosophila blastoderm gene regulatory network is one of the best studied networks in biology. It is composed of a series of tiered sub-networks that act sequentially to generate a primary segmental pattern. Many of these sub-networks have been studied in other arthropods, allowing us to reconstruct how each of them evolved over the transition from the arthropod ancestor to the situation seen in Drosophila today. I trace the evolution of each of these networks, showing how some of them have been modified significantly in Drosophila relative to the ancestral state while others are largely conserved across evolutionary timescales. I compare the putative ancestral arthropod segmentation network with that found in Drosophila and discuss how and why it has been modified throughout evolution, and to what extent this modification is unusual.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade inCaenorhabditisEndoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

A genome-wide assessment of the ancestral neural crest gene regulatory network

The neural crest (NC) is an embryonic cell population that contributes to key vertebrate-specific features including the craniofacial skeleton and peripheral nervous system. Here we examine the transcriptional and epigenomic profiles of NC cells in the sea lamprey, in order to gain insight into the ancestral state of the NC gene regulatory network (GRN). Transcriptome analyses identify clusters of co-regulated genes during NC specification and migration that show high conservation across vertebrates but also identify transcription factors (TFs) and cell-adhesion molecules not previously implicated in NC migration. ATAC-seq analysis uncovers an ensemble of cis-regulatory elements, including enhancers of Tfap2B, SoxE1 and Hox-α2 validated in the embryo. Cross-species deployment of lamprey elements identifies the deep conservation of lamprey SoxE1 enhancer activity, mediating homologous expression in jawed vertebrates. Our data provide insight into the core GRN elements conserved to the base of the vertebrates and expose others that are unique to lampreys.

A novel evolutionary model for constructing gene coexpression networks with comprehensive features

BackgroundUncovering the evolutionary principles of gene coexpression network is important for our understanding of the network topological property of new genes. However, most existing evolutionary models only considered the evolution of duplication genes and only based on the degree of genes, ignoring the other key topological properties. The evolutionary mechanism by which how are new genes integrated into the ancestral networks are not yet to be comprehensively characterized. Herein, based on the human ribonucleic acid-sequencing (RNA-seq) data, we develop a new evolutionary model of gene coexpression network which considers the evolutionary process of both duplication genes and de novo genes.ResultsBased on the human RNA-seq data, we construct a gene coexpression network consisting of 8061 genes and 638624 links. We find that there are 1394 duplication genes and 126 de novo genes in the network. Then based on human gene age data, we reproduce the evolutionary process of this gene coexpression network and develop a new evolutionary model. We find that the generation rates of duplication genes and de novo genes are approximately 3.58/Myr (Myr=Million year) and 0.31/Myr, respectively. Based on the average degree and coreness of parent genes, we find that the gene duplication is a random process. Eventually duplication genes only inherit 12.89% connections from their parent genes and the retained connections have a smaller edge betweenness. Moreover, we find that both duplication genes and de novo genes prefer to develop new interactions with genes which have a large degree and a large coreness. Our proposed model can generate an evolutionary network when the number of newly added genes or the length of evolutionary time is known.ConclusionsGene duplication and de novo genes are two dominant evolutionary forces in shaping the coexpression network. Both duplication genes and de novo genes develop new interactions through a “rich-gets-richer" mechanism in terms of degree and coreness. This mechanism leads to the scale-free property and hierarchical architecture of biomolecular network. The proposed model is able to construct a gene coexpression network with comprehensive biological characteristics.

Topological evolution of coexpression networks by new gene integration maintains the hierarchical and modular structures in human ancestors.

We analyze the global structure and evolution of human gene coexpression networks driven by new gene integration. When the Pearson correlation coefficient is greater than or equal to 0.5, we find that the coexpression network consists of 334 small components and one "giant" connected subnet comprising of 6317 interacting genes. This network shows the properties of power-law degree distribution and small-world. The average clustering coefficient of younger genes is larger than that of the elderly genes (0.6685 vs. 0.5762). Particularly, we find that the younger genes with a larger degree also show a property of hierarchical architecture. The younger genes play an important role in the overall pivotability of the network and this network contains few redundant duplicate genes. Moreover, we find that gene duplication and orphan genes are two dominant evolutionary forces in shaping this network. Both the duplicate genes and orphan genes develop new links through a "rich-gets-richer" mechanism. With the gradual integration of new genes into the ancestral network, most of the topological structure features of the network would gradually increase. However, the exponent of degree distribution and modularity coefficient of the whole network do not change significantly, which implies that the evolution of coexpression networks maintains the hierarchical and modular structures in human ancestors.

Evolution of hierarchy in bacterial metabolic networks

Flow hierarchy is a useful way to characterize the movement of information and matter throughout a network. Hierarchical network organizations are shown to arise when there is a cost of maintaining links in the network. A similar constraint exists in metabolic networks, where costs come from reduced efficiency of nonspecific enzymes or from producing unnecessary enzymes. Previous analyses of bacterial metabolic networks have been used to predict the minimal nutrients that a bacterium needs to grow, its mutualistic relationships with other bacteria, and its major ecological niche.We use metabolic network inference to obtain metabolite flow graphs of 2935 bacterial metabolic networks and find that flow hierarchy evolves independently of modularity and other network properties. By inferring the ancestral metabolic networks and estimating the hierarchical character of the inferred network, we show that hierarchical structure first increased and later decreased over evolutionary history. Furthermore, hierarchical structure in the network is associated with slower growth rates; bacteria with hierarchy scores above the median grow on average 2.25 times faster than those with hierarchy scores below the median.

Bone-like features in skate suggest a novel elasmobranch synapomorphy and deep homology of trabecular mineralization patterns

Bone is a defining characteristic of the vertebrate skeleton, and while chondrichthyans (sharks, skates, and other cartilaginous fishes) are vertebrates, they are hypothesized to have lost the ability to make bone during their evolution. Multiple descriptions of a bone-like tissue in neural arches of vertebrae in various shark species (selachians), however, challenge this hypothesis. Here, we extend this argument by analyzing vertebrae of two members of the batoids (the little skate Leucoraja erinacea and Eaton’s skate Bathyraja eatonii), the sister group to selachians within elasmobranchs. Micro-CT images showed a bone-like mineralization pattern in neural arches of each skate species, and histological analyses confirmed that this bone-like tissue surrounded a cartilage core, exactly as described in sharks. Another mineralization pattern identified in skate vertebrae was distinct from the polygonal tesseral and areolar patterns that classically are associated with the chondrichthyan endoskeleton. Many regions of the vertebrae, including the neural spine and transverse processes, showed this perichondral mineralization pattern, termed here trabecular tesseral. Other than the cartilage core of the neural arch, all mineralized tissues in skate vertebrae had flattened cells surrounded by matrix with bone-like histology. Analyses of quantitative microstructural parameters revealed that, compared to rat vertebrae, the bone-like mineralization pattern in the neural arches of skate vertebrae was more similar to compact bone than trabecular bone. In contrast, the thickness of the trabecular tesseral pattern was more similar to trabecular bone than compact bone of rat vertebrae. In conclusion, a bone-like tissue in neural arches of skate vertebrae appears to be a novel elasmobranch synapomorphy. We propose that the trabecular tesseral mineralization pattern in the skate might have deep homology to the mineralization pattern utilized in trabecular bone. Statement of SignificanceMineralization patterns of skeletal tissues have not been investigated thoroughly in all vertebrate clades. Despite their designation as ‘cartilaginous fish’, chondrichthyans clearly evolved from ancestral vertebrates that made bone. The consensus that chondrichthyans lost the ability to make bone during their evolution, however, is challenged by reports of bone and bone-like tissues in the neural arches of vertebrae in extant sharks (selachians). Here, we provide evidence from micro-CT imaging and histological analyses to support our hypothesis that a bone-like tissue is present in the neural arches of batoids (the sister group to selachians within elasmobranchs). These results argue strongly that the neural arch bone-like tissue is a previously unknown synapomorphy of elasmobranchs. In addition to the bone-like mineralization pattern identified in the neural arches, micro-CT images also showed a novel mineralization pattern which we described as trabecular tesseral. Quantitative microstructural features shared between trabecular tesseral pattern and trabecular bone (from homologous rat vertebrae) suggest that both patterns might derive from an ancestral gene network driving trabecular mineralization (i.e., deep homology).

Towards reconstructing the ancestral brain gene-network regulating caste differentiation in ants.

Specialized queens and life-time unmated workers evolved once in the common ancestor of all ants, but whether caste development across ants continues to be at least partly regulated by a single core set of genes remains obscure. We analysed brain transcriptomes from five ant species (three subfamilies) and reconstructed the origins of genes with caste-biased expression. Ancient genes predating the Neoptera were more likely to regulate gyne (virgin queen) phenotypes, while caste differentiation roles of younger, ant-lineage-specific genes varied. Transcriptome profiling showed that the ancestral network for caste-specific gene-regulation has been maintained, but that signatures of common ancestry are obscured by later modifications. Adjusting for such differences, we identified a core gene-set that: 1. consistently displayed similar directions and degrees of caste-differentiated expression, and 2. have mostly not been reported as being involved in caste differentiation. These core regulatory genes exist in the genomes of ant species that secondarily lost the queen caste, but expression differences for reproductive and sterile workers are minor and similar to social paper wasps that lack differentiated castes. Many caste-biased ant genes have caste-differentiated expression in honeybees, but directions of caste bias were uncorrelated, as expected when permanent castes evolved independently in both lineages.

Resurrected Protein Interaction Networks Reveal the Innovation Potential of Ancient Whole-Genome Duplication.

The evolution of plants is characterized by whole-genome duplications, sometimes closely associated with the origin of large groups of species. The gamma (γ) genome triplication occurred at the origin of the core eudicots, which comprise ∼75% of flowering plants. To better understand the impact of whole-genome duplication, we studied the protein interaction network of MADS domain transcription factors, which are key regulators of reproductive development. We reconstructed, synthesized, and tested the interactions of ancestral proteins immediately before and closely after the triplication and directly compared these ancestral networks to the extant networks of Arabidopsis thaliana and tomato (Solanum lycopersicum). We found that gamma expanded the MADS domain interaction network more strongly than subsequent genomic events. This event strongly rewired MADS domain interactions and allowed for the evolution of new functions and installed robustness through new redundancy. Despite extensive rewiring, the organization of the network was maintained through gamma. New interactions and protein retention compensated for its potentially destructive impact on network organization. Post gamma, the network evolved from an organization around the single hub SEP3 to a network organized around multiple hubs and well-connected proteins lost, rather than gained, interactions. The data provide a resource for comparative developmental biology in flowering plants.

Species-Specific Changes in a Primate Transcription Factor Network Provide Insights into the Molecular Evolution of the Primate Prefrontal Cortex.

The human prefrontal cortex (PFC) differs from that of other primates with respect to size, histology, and functional abilities. Here, we analyzed genome-wide expression data of humans, chimpanzees, and rhesus macaques to discover evolutionary changes in transcription factor (TF) networks that may underlie these phenotypic differences. We determined the co-expression networks of all TFs with species-specific expression including their potential target genes and interaction partners in the PFC of all three species. Integrating these networks allowed us inferring an ancestral network for all three species. This ancestral network as well as the networks for each species is enriched for genes involved in forebrain development, axonogenesis, and synaptic transmission. Our analysis allows us to directly compare the networks of each species to determine which links have been gained or lost during evolution. Interestingly, we detected that most links were gained on the human lineage, indicating increase TF cooperativity in humans. By comparing network changes between different tissues, we discovered that in brain tissues, but not in the other tissues, the human networks always had the highest connectivity. To pinpoint molecular changes underlying species-specific phenotypes, we analyzed the sub-networks of TFs derived only from genes with species-specific expression changes in the PFC. These sub-networks differed significantly in structure and function between the human and chimpanzee. For example, the human-specific sub-network is enriched for TFs implicated in cognitive disorders and for genes involved in synaptic plasticity and cognitive functions. Our results suggest evolutionary changes in TF networks that might have shaped morphological and functional differences between primate brains, in particular in the human PFC.

The regulatory ancestral network of surgical meshes.

BackgroundAll surgical meshes entering the U.S. market have been cleared for clinical use by the 510(k) process of the Food and Drug Administration (FDA), in which devices simply require proof of “substantial equivalence” to predicate devices, without the need for clinical trials. However, recalled meshes associated with adverse effects may, indirectly, continue to serve as predicates for new devices raising concerns over the safety of the 510(k) route.MethodologyHere we assess the potential magnitude of this problem by determining the ancestral network of equivalence claims linking recently cleared surgical meshes. Using the FDA website we identified all surgical meshes cleared by the 510(k) route between January 2013 and December 2015 along with all listed predicates for these devices. Using a network approach, we trace the ancestry of predicates across multiple generations of equivalence claims and identify those meshes connected to devices that have since recalled from the market along with the reason for their recall.ConclusionsWe find that the 77 surgical meshes cleared between 2013 and 2015 are based on 771 interconnected predicate claims of equivalence from 400 other devices. The vast majority of these devices (97%) are descended from only six surgical meshes that were present on the market prior to 1976. One of these ancestral meshes alone, provided the basis of 183 subsequent devices. Furthermore, we show that 16% of recently cleared devices are connected through equivalence claims to the 3 predicate meshes that have been recalled for design and material related flaws causing serious adverse events. Taken together, our results show that surgical meshes are connected through a tangled web of equivalency claims and many meshes recently cleared by the FDA have connections through chains of equivalency to devices which have been recalled from the market due to concerns over clinical safety. These findings raise concerns over the efficacy of the 510(k) route in ensuring patient safety.

Structure versus time in the evolutionary diversification of avian carotenoid metabolic networks.

Historical associations of genes and proteins are thought to delineate pathways available to subsequent evolution; however, the effects of past functional involvements on contemporary evolution are rarely quantified. Here, we examined the extent to which the structure of a carotenoid enzymatic network persists in avian evolution. Specifically, we tested whether the evolution of carotenoid networks was most concordant with phylogenetically structured expansion from core reactions of common ancestors or with subsampling of biochemical pathway modules from an ancestral network. We compared structural and historical associations in 467 carotenoid networks of extant and ancestral species and uncovered the overwhelming effect of pre-existing metabolic network structure on carotenoid diversification over the last 50million years of avian evolution. Over evolutionary time, birds repeatedly subsampled and recombined conserved biochemical modules, which likely maintained the overall structure of the carotenoid metabolic network during avian evolution. These findings explain the recurrent convergence of evolutionary distant species in carotenoid metabolism and weak phylogenetic signal in avian carotenoid evolution. Remarkable retention of an ancient metabolic structure throughout extensive and prolonged ecological diversification in avian carotenoid metabolism illustrates a fundamental requirement of organismal evolution - historical continuity of a deterministic network that links past and present functional associations of its components.

Origin of the Animal Circadian Clock: Diurnal and Light-Entrained Gene Expression in the Sponge Amphimedon queenslandica

The circadian clock is a molecular network that coordinates organismal behavior and physiology with daily environmental changes in the day-night cycle. In eumetazoans (bilaterians + cnidarians), this network appears to be largely conserved, yet different from other known eukaryotic circadian networks. To determine if the eumetazoan circadian network has an older origin, we ask here whether orthologues comprising this network are expressed in a manner consistent with a role in regulating circadian patterns in a representative of an earlier-branching animal lineage, the sponge Amphimedon queenslandica. The A. queenslandica genome encodes orthologues of many eumetazoan circadian genes, including two cryptochrome genes that encode flavoproteins, three Timeout genes, and two PAR-bZIP and seven bHLH-PAS transcription factor genes. There is no apparent Cycle orthologue, although we can identify three closely related ARNT genes. Of the putative circadian genes, only AqPARa and AqCry2 have a consistent oscillating diurnal expression profile, and the rhythmic expression of both these genes is partially lost when the animals are exposed to constant light or darkness. Expression of the other putative circadian genes, in particular AqClock, is neither diurnally-oscillating nor light-dependent. AqPARa and AqCry2 are also temporally and spatially co-expressed throughout embryonic and larval development. Transcripts of these genes are enriched first in cells comprising the larval posterior pigment ring, which is a simple photosensory organ that is responsible for the negative phototactic behaviour displayed by larvae, and subsequently in the larval epithelial and subepithelial layers. The combined findings of no clear Cycle orthologue and of PAR-bZIP and cryptochrome being the only orthologues expressed in a pattern consistent with a circadian role suggests that either (i) the ancestral metazoan circadian network was simpler than the eumetazoan network, or (ii) that this sponge has lost some components, as has occurred in some other animals such as Hydra.

Parallel reorganization of protein function in the spindle checkpoint pathway through evolutionary paths in the fitness landscape that appear neutral in laboratory experiments.

Regulatory networks often increase in complexity during evolution through gene duplication and divergence of component proteins. Two models that explain this increase in complexity are: 1) adaptive changes after gene duplication, such as resolution of adaptive conflicts, and 2) non-adaptive processes such as duplication, degeneration and complementation. Both of these models predict complementary changes in the retained duplicates, but they can be distinguished by direct fitness measurements in organisms with short generation times. Previously, it has been observed that repeated duplication of an essential protein in the spindle checkpoint pathway has occurred multiple times over the eukaryotic tree of life, leading to convergent protein domain organization in its duplicates. Here, we replace the paralog pair in S. cerevisiae with a single-copy protein from a species that did not undergo gene duplication. Surprisingly, using quantitative fitness measurements in laboratory conditions stressful for the spindle-checkpoint pathway, we find no evidence that reorganization of protein function after gene duplication is beneficial. We then reconstruct several evolutionary intermediates from the inferred ancestral network to the extant one, and find that, at the resolution of our assay, there exist stepwise mutational paths from the single protein to the divergent pair of extant proteins with no apparent fitness defects. Parallel evolution has been taken as strong evidence for natural selection, but our results suggest that even in these cases, reorganization of protein function after gene duplication may be explained by neutral processes.