Comparative genomics of the vertebrate insulin/TOR signal transduction pathway: A network-level analysis of selective pressures

Complexity of biological function relies on large networks of interacting molecules. However, the evolutionary properties of these networks are not fully understood. It has been shown that selective pressures depend on the position of genes in the network. We have previously shown that in the Drosophila insulin/target of rapamycin (TOR) signal transduction pathway there is a correlation between the pathway position and the strength of purifying selection, with the downstream genes being most constrained. In this study, we investigated the evolutionary dynamics of this well-characterized pathway in vertebrates. More specifically, we determined the impact of natural selection on the evolution of 72 genes of this pathway. We found that in vertebrates there is a similar gradient of selective constraint in the insulin/TOR pathway to that found in Drosophila. This feature is neither the result of a polarity in the impact of positive selection nor of a series of factors affecting selective constraint levels (gene expression level and breadth, codon bias, protein length, and connectivity). We also found that pathway genes encoding physically interacting proteins tend to evolve under similar selective constraints. The results indicate that the architecture of the vertebrate insulin/TOR pathway constrains the molecular evolution of its components. Therefore, the polarity detected in Drosophila is neither specific nor incidental of this genus. Hence, although the underlying biological mechanisms remain unclear, these may be similar in both vertebrates and Drosophila.

Constitutive expression of FGF4 disrupts the development of the eye and the anterior CNS during mouse embryogenesi, but does not influence the expressionof shh in these areas.

The biological consequences of constitutive fibroblast growth factor-4 (fgf4) expression have been analysed during anterior CNS development of mouse chimeric embryos. Severe mutant embryos exhibit exencephaly, absence of eye development and anomalous differentiation of neuropithelium. These embryos also show ectopic limb buds resembling the early phases of limb development. Because our results show that anterior CNS in those chimeric embrios does not express shh ectopically, we suggest that malformations may be due to interference between the ectopic expression of fgf4 in the cephalic area and the receptors for the members of the FGF family that regulate brain and eye development, namely fgf8. If this is correct, the results indirectly suport the crucial role of fgf8 in patterning the anterior CNS.

Maintenance of A/P body regions in planarians by TCEN49, a putative cystine-knot neurotrophin.

In freshwater planarians, the protein TCEN49 has been linked to the regional specification of the central body region, which includes the pharynx.

Comment on 'The Molecular Evolutionary Patterns of the Insulin/FOXO Signaling Pathway'

Letter to the Editor on Wang M, Wang Q, Wang Z, Zhang X, Pan Y. The molecular evolutionary patterns of the insulin/FOXO signaling pathway

Competition between SOCS36E and Drk modulates Sevenless receptor tyrosine kinase activity

Modulation of signalling pathways can trigger different cellular responses, including differences in cell fate. This modulation can be achieved by controlling the pathway activity with great precision to ensure robustness and reproducibility of the specification of cell fate. The development of the photoreceptor R7 in the Drosophila melanogasterretina has become a model in which to investigate the control of cell signalling. During R7 specification, a burst of Ras small GTPase (Ras) and mitogen-activated protein kinase (MAPK) controlled by Sevenless receptor tyrosine kinase (Sev) is required. Several cells in each ommatidium express sev. However, the spatiotemporal expression of the boss ligand and the action of negative regulators of the Sev pathway will restrict the R7 fate to a single cell. The Drosophila suppressor of cytokine signalling 36E (SOCS36E) protein contains an SH2 domain and acts as a Sev signalling attenuator. By contrast, downstream of receptor kinase (Drk), the fly homolog of the mammalian Grb2 adaptor protein, which also contains an SH2 domain, acts as a positive activator of the pathway. Here, we apply the Förster resonance energy transfer (FRET) assay to transfected Drosophila S2 cells and demonstrate that Sev binds directly to either the suppressor protein SOCS36E or the adaptor protein Drk. We propose a mechanistic model in which the competition between these two proteins for binding to the same docking site results in either attenuation of the Sev transduction in cells that should not develop R7 photoreceptors or amplification of the Ras-MAPK signal only in the R7 precursor.

RNA interference is ineffective as a routine method for gene silencing in chick embryos as monitored by fgf8 silencing.

The in vivo accessibility of the chick embryo makes it a favoured model system for experimental developmental biology. Although the range of available techniques now extends to miss-expression of genes through in ovo electroporation, it remains difficult to knock out individual gene expression. Recently, the possibility of silencing gene expression by RNAi in chick embryos has been reported. However, published studies show only discrete quantitative differences in the expression of the endogenous targeted genes and unclear morphological alterations. To elucidate whether the tools currently available are adequate to silence gene expression sufficiently to produce a clear and specific null-like mutant phenotype, we have performed several experiments with different molecules that trigger RNAi: dsRNA, siRNA, and shRNA produced from a plasmid coexpressing green fluorescent protein as an internal marker. Focussing on fgf8 expression in the developing isthmus, we show that no morphological defects are observed, and that fgf8 expression is neither silenced in embryos microinjected with dsRNA nor in embryos microinjected and electroporated with a pool of siRNAs. Moreover, fgf8 expression was not significantly silenced in most isthmic cells transformed with a plasmid producing engineered shRNAs to fgf8. We also show that siRNA molecules do not spread significantly from cell to cell as reported for invertebrates, suggesting the existence of molecular differences between different model systems that may explain the different responses to RNAi. Although our results are basically in agreement with previously reported studies, we suggest, in contrast to them, that with currently available tools and techniques the number of cells in which fgf8 gene expression is decreased, if any, is not sufficient to generate a detectable mutant phenotype, thus making RNAi useless as a routine method for functional gene analysis in chick embryos.

Patterns of nucleotide diversity at the regions encompassing the Drosophila insulin-like peptide (dilp) genes: demography vs positive selection in Drosophila melanogaster.

In Drosophila, the insulin-signaling pathway controls some life history traits, such as fertility and lifespan, and it is considered to be the main metabolic pathway involved in establishing adult body size. Several observations concerning variation in body size in the Drosophila genus are suggestive of its adaptive character. Genes encoding proteins in this pathway are, therefore, good candidates to have experienced adaptive changes and to reveal the footprint of positive selection. The Drosophila insulin-like peptides (DILPs) are the ligands that trigger the insulin-signaling cascade. In Drosophila melanogaster, there are several peptides that are structurally similar to the single mammalian insulin peptide. The footprint of recent adaptive changes on nucleotide variation can be unveiled through the analysis of polymorphism and divergence. With this aim, we have surveyed nucleotide sequence variation at the dilp1-7 genes in a natural population of D. melanogaster. The comparison of polymorphism in D. melanogaster and divergence from D. simulans at different functional classes of the dilp genes provided no evidence of adaptive protein evolution after the split of the D. melanogaster and D. simulans lineages. However, our survey of polymorphism at the dilp gene regions of D. melanogaster has provided some evidence for the action of positive selection at or near these genes. The regions encompassing the dilp1-4 genes and the dilp6 gene stand out as likely affected by recent adaptive events.

Gradients and regeneration: the case of TNEX59 in planarians (Platyhelminthes, Tricladida)

During development, multicellular organisms determine and then differentiate regions that constitute the body's architecture. These regions are established and controlled by a number of molecules, including nuclear factors, that drive the organism from the egg to its final shape. We studied the molecules involved in the regionalisation of the freshwater planarian body (Platyhelminthes, Turbelleria, Tricladida).

Comparative genomics of Hedgehog loci in chordates and the origin of Shh regulatory novelties

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications.

Broken colinearity of the amphioxus Hox cluster.

Background In most eumetazoans studied so far, Hox genes determine the identity of structures along the main body axis. They are usually linked in genomic clusters and, in the case of the vertebrate embryo, are expressed with spatial and temporal colinearity. Outside vertebrates, temporal colinearity has been reported in the cephalochordate amphioxus (the least derived living relative of the chordate ancestor) but only for anterior and central genes, namely Hox1 to Hox4 and Hox6. However, most of the Hox gene expression patterns in amphioxus have not been reported. To gain global insights into the evolution of Hox clusters in chordates, we investigated a more extended expression profile of amphioxus Hox genes. Results Here we report an extended expression profile of the European amphioxus Branchiostoma lanceolatum Hox genes and describe that all Hox genes, except Hox13, are expressed during development. Interestingly, we report the breaking of both spatial and temporal colinearity for at least Hox6 and Hox14, which thus have escaped from the classical Hox code concept. We show a previously unidentified Hox6 expression pattern and a faint expression for posterior Hox genes in structures such as the posterior mesoderm, notochord, and hindgut. Unexpectedly, we found that amphioxus Hox14 had the most divergent expression pattern. This gene is expressed in the anterior cerebral vesicle and pharyngeal endoderm. Amphioxus Hox14 expression represents the first report of Hox gene expression in the most anterior part of the central nervous system. Nevertheless, despite these divergent expression patterns, amphioxus Hox6 and Hox14 seem to be still regulated by retinoic acid. Conclusions Escape from colinearity by Hox genes is not unusual in either vertebrates or amphioxus and we suggest that those genes escaping from it are probably associated with the patterning of lineage-specific morphological traits, requiring the loss of those developmental constraints that kept them colinear.

Comment on 'The Molecular Evolutionary Patterns of the Insulin/FOXO Signaling Pathway'

Letter to the Editor on Wang M, Wang Q, Wang Z, Zhang X, Pan Y. The molecular evolutionary patterns of the insulin/FOXO signaling pathway

Uncovering the functional constraints underlying the genomic organisation of the Odorant-Binding Protein genes

Animal olfactory systems have a critical role for the survival and reproduction of individuals. In insects, the odorant-binding proteins (OBPs) are encoded by a moderately sized gene family, and mediate the first steps of the olfactory processing. Most OBPs are organized in clusters of a few paralogs, which are conserved over time. Currently, the biological mechanism explaining the close physical proximity among OBPs is not yet established. Here, we conducted a comprehensive study aiming to gain insights into the mechanisms underlying the OBP genomic organization. We found that the OBP clusters are embedded within large conserved arrangements. These organizations also include other non-OBP genes, which often encode proteins integral to plasma membrane. Moreover, the conservation degree of such large clusters is related to the following: 1) the promoter architecture of the confined genes, 2) a characteristic transcriptional environment, and 3) the chromatin conformation of the chromosomal region. Our results suggest that chromatin domains may restrict the location of OBP genes to regions having the appropriate transcriptional environment, leading to the OBP cluster structure. However, the appropriate transcriptional environment for OBP and the other neighbor genes is not dominated by reduced levels of expression noise. Indeed, the stochastic fluctuations in the OBP transcript abundance may have a critical role in the combinatorial nature of the olfactory coding process.

Acoel flatworms are not Platyhelminthes: evidence from phylogenomics

Acoel flatworms are small marine worms traditionally considered to belong to the phylum Platyhelminthes. However, molecular phylogenetic analyses suggest that acoels are not members of Platyhelminthes, but are rather extant members of the earliest diverging Bilateria. This result has been called into question, under suspicions of a long branch attraction (LBA) artefact. Here we re-examine this problem through a phylogenomic approach using 68 different protein-coding genes from the acoel Convoluta pulchra and 51 metazoan species belonging to 15 different phyla. We employ a mixture model, named CAT, previously found to overcome LBA artefacts where classical models fail. Our results unequivocally show that acoels are not part of the classically defined Platyhelminthes, making the latter polyphyletic. Moreover, they indicate a deuterostome affinity for acoels, potentially as a sister group to all deuterostomes, to Xenoturbellida, to Ambulacraria, or even to chordates. However, the weak support found for most deuterostome nodes, together with the very fast evolutionary rate of the acoel Convoluta pulchra, call for more data from slowly evolving acoels (or from its sister-group, the Nemertodermatida) to solve this challenging phylogenetic problem.

Genome-wide analysis of adaptive molecular evolution in the carnivorous plant Utricularia gibba

The genome of the bladderwort Utricularia gibba provides an unparalleled opportunity to uncover the adaptive landscape of an aquatic carnivorous plant with unique phenotypic features such as absence of roots, development of water-filled suction bladders, and a highly ramified branching pattern. Despite its tiny size, the U. gibba genome accommodates approximately as many genes as other plant genomes. To examine the relationship between the compactness of its genome and gene turnover, we compared the U. gibba genome with that of four other eudicot species, defining a total of 17,324 gene families (orthogroups). These families were further classified as either 1) lineage-specific expanded/contracted or 2) stable in size. The U. gibba-expanded families are generically related to three main phenotypic features: 1) trap physiology, 2) key plant morphogenetic/developmental pathways, and 3) response to environmental stimuli, including adaptations to life in aquatic environments. Further scans for signatures of protein functional specialization permitted identification of seven candidate genes with amino acid changes putatively fixed by positive Darwinian selection in the U. gibba lineage. The Arabidopsis orthologs of these genes (AXR, UMAMIT41, IGS, TAR2, SOL1, DEG9, and DEG10) are involved in diverse plant biological functions potentially relevant for U. gibba phenotypic diversification, including 1) auxin metabolism and signal transduction, 2) flowering induction and floral meristem transition, 3) root development, and 4) peptidases. Taken together, our results suggest numerous candidate genes and gene families as interesting targets for further experimental confirmation of their functional and adaptive roles in the U. gibba's unique lifestyle and highly specialized body plan.

Developmental and tissue-specific involvement of peroxisome proliferator-activated receptor-alpha in the control of mouse uncoupling protein-3 gene expression.

Uncoupling protein-3 (UCP3) is a member of the mitochondrial carrier family expressed preferentially in skeletal muscle and heart. It appears to be involved in metabolic handling of fatty acids in a way that minimizes excessive production of reactive oxygen species. Fatty acids are powerful regulators of UCP3 gene transcription. We have found that the role of peroxisome proliferator-activated receptor-α (PPARα) on the control of UCP3 gene expression depends on the tissue and developmental stage. In adults, UCP3 mRNA expression is unaltered in skeletal muscle from PPARα-null mice both in basal conditions and under the stimulus of starvation. In contrast, UCP3 mRNA is down-regulated in adult heart both in fed and fasted PPARα-null mice. This occurs despite the increased levels of free fatty acids caused by fasting in PPARα-null mice. In neonates, PPARα-null mice show impaired UCP3 mRNA expression in skeletal muscle in response to milk intake, and this is not a result of reduced free fatty acid levels. The murine UCP3 promoter is activated by fatty acids through either PPARα or PPARδ but not by PPARγ or retinoid X receptor alone. PPARδ-dependent activation could be a potential compensatory mechanism to ensure appropriate expression of UCP3 gene in adult skeletal muscle in the absence of PPARα. However, among transcripts from other PPARα and PPARδ target genes, only those acutely induced by milk intake in wild-type neonates were altered in muscle or heart from PPARα-null neonates. Thus, PPARα-dependent regulation is required for appropriate gene regulation of UCP3 as part of the subset of fatty-acid-responsive genes in neonatal muscle and heart.