Duplication events and the evolution of segmental identity

Genome-wide identification and expression profiling of growth‑regulating factor (GRF) and GRF‑interacting factor (GIF) gene families in chickpea and pigeonpea

The growth-regulating factor (GRF) and GRF-interacting factor (GIF) families encode plant-specific transcription factors and play vital roles in plant development and stress response processes. Although GRF and GIF genes have been identified in various plant species, there have been no reports of the analysis and identification of the GRF and GIF transcription factor families in chickpea (Cicer arietinum) and pigeonpea (Cajanus cajan). The present study identified seven CaGRFs, eleven CcGRFs, four CaGIFs, and four CcGIFs. The identified proteins were grouped into eight and three clades for GRFs and GIFs, respectively based on their phylogenetic relationships. A comprehensive in-silico analysis was performed to determine chromosomal location, sub-cellular localization, and types of regulatory elements present in the putative promoter region. Synteny analysis revealed that GRF and GIF genes showed diploid-polyploid topology in pigeonpea, but not in chickpea. Tissue-specific expression data at the vegetative and reproductive stages of the plant showed that GRFs and GIFs were strongly expressed in tissues like embryos, pods, and seeds, indicating that GRFs and GIFs play vital roles in plant growth and development. This research characterized GRF and GIF families and hints at their primary roles in the chickpea and pigeonpea growth and developmental process. Our findings provide potential gene resources and vital information on GRF and GIF gene families in chickpea and pigeonpea, which will help further understand the regulatory role of these gene families in plant growth and development.

Genome-Wide Investigation of the CRF Gene Family in Maize and Functional Analysis of ZmCRF9 in Response to Multiple Abiotic Stresses

The cytokinin response factors (CRFs) are pivotal players in regulating plant growth, development, and responses to diverse stresses. Despite their significance, comprehensive information on CRF genes in the primary food crop, maize, remains scarce. In this study, a genome-wide analysis of CRF genes in maize was conducted, resulting in the identification of 12 members. Subsequently, we assessed the chromosomal locations, gene duplication events, evolutionary relationships, conserved motifs, and gene structures of all ZmCRF members. Analysis of ZmCRF promoter regions indicated the presence of cis-regulatory elements associated with plant growth regulation, hormone response, and various abiotic stress responses. The expression patterns of maize CRF genes, presented in heatmaps, exhibited distinctive patterns of tissue specificity and responsiveness to multiple abiotic stresses. qRT-PCR experiments were conducted on six selected genes and confirmed the involvement of ZmCRF genes in the plant's adaptive responses to diverse environmental challenges. In addition, ZmCRF9 was demonstrated to positively regulate cold and salt tolerance. Ultimately, we explored the putative interaction partners of ZmCRF proteins. In summary, this systematic overview and deep investigation of ZmCRF9 provides a solid foundation for further exploration into how these genes contribute to the complex interplay of plant growth, development, and responses to stress.

Genome-wide identification, phylogeny, and expression analysis of GRF transcription factors in pineapple (Ananas comosus)

Background

Pineapple is the only commercially grown fruit crop in the Bromeliaceae family and has significant agricultural, industrial, economic, and ornamental value. GRF (growth-regulating factor) proteins are important transcription factors that have evolved in seed plants (embryophytes). They contain two conserved domains, QLQ (Gln, Leu, Gln) and WRC (Trp, Arg, Cys), and regulate multiple aspects of plant growth and stress response, including floral organ development, leaf growth, and hormone responses. The GRF family has been characterized in a number of plant species, but little is known about this family in pineapple and other bromeliads.

Main discoveries

We identified eight GRF transcription factor genes in pineapple, and phylogenetic analysis placed them into five subfamilies (I, III, IV, V, VI). Segmental duplication appeared to be the major contributor to expansion of the AcGRF family, and the family has undergone strong purifying selection during evolution. Relative to that of other gene families, the gene structure of the GRF family showed less conservation. Analysis of promoter cis-elements suggested that AcGRF genes are widely involved in plant growth and development. Transcriptome data and qRT-PCR results showed that, with the exception of AcGRF5, the AcGRFs were preferentially expressed in the early stage of floral organ development and AcGRF2 was strongly expressed in ovules. Gibberellin treatment significantly induced AcGRF7/8 expression, suggesting that these two genes may be involved in the molecular regulatory pathway by which gibberellin promotes pineapple fruit expansion.

Conclusion

AcGRF proteins appear to play a role in the regulation of floral organ development and the response to gibberellin. The information reported here provides a foundation for further study of the functions of AcGRF genes and the traits they regulate.

Evolution and functional analysis of the GRAS family genes in six Rosaceae species

BACKGROUND

GRAS genes formed one of the important transcription factor gene families in plants, had been identified in several plant species. The family genes were involved in plant growth, development, and stress resistance. However, the comparative analysis of GRAS genes in Rosaceae species was insufficient.

RESULTS

In this study, a total of 333 GRAS genes were identified in six Rosaceae species, including 51 in strawberry (Fragaria vesca), 78 in apple (Malus domestica), 41 in black raspberry (Rubus occidentalis), 59 in European pear (Pyrus communis), 56 in Chinese rose (Rosa chinensis), and 48 in peach (Prunus persica). Motif analysis showed the VHIID domain, SAW motif, LR I region, and PFYRE motif were considerably conserved in the six Rosaceae species. All GRAS genes were divided into 10 subgroups according to phylogenetic analysis. A total of 15 species-specific duplicated clades and 3 lineage-specific duplicated clades were identified in six Rosaceae species. Chromosomal localization presented the uneven distribution of GRAS genes in six Rosaceae species. Duplication events contributed to the expression of the GRAS genes, and Ka/Ks analysis suggested the purification selection as a major force during the evolution process in six Rosaceae species. Cis-acting elements and GO analysis revealed that most of the GRAS genes were associated with various environmental stress in six Rosaceae species. Coexpression network analysis showed the mutual regulatory relationship between GRAS and bZIP genes, suggesting the ability of the GRAS gene to regulate abiotic stress in woodland strawberry. The expression pattern elucidated the transcriptional levels of FvGRAS genes in various tissues and the drought and salt stress in woodland strawberry, which were verified by RT-qPCR analysis.

CONCLUSIONS

The evolution and functional analysis of GRAS genes provided insights into the further understanding of GRAS genes on the abiotic stress of Rosaceae species.

Complexity and diversity of motion amplification and control strategies in motile carnivorous plant traps

Similar to animals, plants have evolved mechanisms for elastic energy storage and release to power and control rapid motion, yet both groups have been largely studied in isolation. This is exacerbated by the lack of consistent terminology and conceptual frameworks describing elastically powered motion in both groups. Iconic examples of fast movements can be found in carnivorous plants, which have become important models to study biomechanics, developmental processes, evolution and ecology. Trapping structures and processes vary considerably between different carnivorous plant groups. Using snap traps, suction traps and springboard-pitfall traps as examples, we illustrate how traps mix and match various mechanisms to power, trigger and actuate motions that contribute to prey capture, retention and digestion. We highlight a fundamental trade-off between energetic investment and movement control and discuss it in a functional-ecological context.

Cloning and expression of Hoxc6 gene from Pampus argenteus and its relationship with pelvic fin absence

Hoxc6 gene can be described as having roles in axial patterning in early embryogenesis, and in at least some species, having a contribution to limb positioning. In this study, we cloned and characterised Pampus argenteus Hoxc6. The highly conserved HOXC6 protein sequence contains a homeodomain and a low-complexity region. Expression of Hoxc6 mRNA was measured at different developmental stages and in different tissues by real-time PCR (p < 0.05), and was high during eye capsule and brain differentiation stages, but low in 7 and 13-day-old larvae. Hoxc6 mRNA was more abundant in fin tissue than brain and eye tissues. Western blotting showed that HOXC6 protein levels were high at embryonic stages, but decreased significantly in 7, 13, 16 and 19-day-old larvae, and levels were essentially consistent with those of mRNA measured by real-time PCR in different tissues. In situ hybridisation showed that the Hoxc6 transcript was strongly expressed in the whole brain and anterior part of the body axis in 1-day-old larvae, but in the hindbrain, pectoral fin, mandible and hypothetical pelvic fin region in 7, 13, 16 and 19-day-old organisms. These results clarify the expression and localisation characteristics of Hoxc6 gene in P. argenteus, and provide a theoretical basis for the molecular mechanism of pelvic fin loss in silver pomfret.

Fate of Duplicated Neural Structures

Statistical physics determines the abundance of different arrangements of matter depending on cost-benefit balances. Its formalism and phenomenology percolate throughout biological processes and set limits to effective computation. Under specific conditions, self-replicating and computationally complex patterns become favored, yielding life, cognition, and Darwinian evolution. Neurons and neural circuits sit at a crossroads between statistical physics, computation, and (through their role in cognition) natural selection. Can we establish a statistical physics of neural circuits? Such theory would tell what kinds of brains to expect under set energetic, evolutionary, and computational conditions. With this big picture in mind, we focus on the fate of duplicated neural circuits. We look at examples from central nervous systems, with stress on computational thresholds that might prompt this redundancy. We also study a naive cost-benefit balance for duplicated circuits implementing complex phenotypes. From this, we derive phase diagrams and (phase-like) transitions between single and duplicated circuits, which constrain evolutionary paths to complex cognition. Back to the big picture, similar phase diagrams and transitions might constrain I/O and internal connectivity patterns of neural circuits at large. The formalism of statistical physics seems to be a natural framework for this worthy line of research.

Multiple independent chromosomal fusions accompanied the radiation of the Antarctic teleost genus Trematomus (Notothenioidei:Nototheniidae)

Background

Chromosomal rearrangements are thought to be an important driving force underlying lineage diversification, but their link to speciation continues to be debated. Antarctic teleost fish of the family Nototheniidae (Notothenioidei) diversified in a changing environmental context, which led to ecological, morphological, and genetic differentiation among populations. In addition, extensive chromosomal repatterning accompanied species divergence in several clades. The most striking karyotypic changes involved the recent species radiation (about 10 My) of the genus Trematomus, with chromosomal pair numbers ranging between 29 and 12. These dramatic reductions in chromosome number resulted mostly from large-scale chromosome fusions. Multiple centric and/or tandem fusions have been hypothesized in at least seven of the twelve recognized Trematomus species. To reconstruct their evolutionary history, we employed comparative cytogenomics (BAC-FISH and chromosome painting) to reveal patterns of interspecific chromosomal orthologies across several notothenioid clades.

Results

We defined orthologous chromosomal segments of reference, termed Structural Units (SUs). SUs were identified in a total of 18 notothenioid species. We demonstrated for the first time that SUs were strongly conserved across every specimen examined, with chromosomal syntenies highlighting a paucity of intrachromosomal macro-rearrangements. Multiple independent fusions of these SUs were inferred in the Trematomus species, in contrast to the shared SU fusions in species of the sister lineage Notothenia.

Conclusions

The SU segments were defined units of chromosomal rearrangement in the entire family Nototheiidae, which diverged from the other notothenioid families 20 My ago. Some of the identified chromosomal syntenies within the SUs were even conserved in their closest relatives, the family Eleginopsidae. Comparing the timing of acquisition of the fusions in the closely related genera Notothenia and Trematomus of the nototheniid species family, we conclude that they exhibit distinct chromosomal evolutionary histories, which may be relevant to different speciation scenarios.

MIPhy: identify and quantify rapidly evolving members of large gene families

After transitioning to a new environment, species often exhibit rapid phenotypic innovation. One of the fastest mechanisms for this is duplication followed by specialization of existing genes. When this happens to a member of a gene family, it tends to leave a detectable phylogenetic signature of lineage-specific expansions and contractions. These can be identified by analyzing the gene family across several species and identifying patterns of gene duplication and loss that do not correlate with the known relationships between those species. This signature, termed phylogenetic instability, has been previously linked to adaptations that change the way an organism samples and responds to its environment; conversely, low phylogenetic instability has been previously linked to proteins with endogenous functions. With the increase in genome-level data, there is a need to identify and quantify phylogenetic instability. Here, we present Minimizing Instability in Phylogenetics (MIPhy), a tool that solves this problem by quantifying the incongruence of a gene's evolutionary history. The motivation behind MIPhy was to produce a tool to aid in interpreting phylogenetic trees. It can predict which members of a gene family are under adaptive evolution, working only from a gene tree and the relationship between the species under consideration. While it does not conduct any estimation of positive selection-which is the typical indication of adaptive evolution-the results tend to agree. We demonstrate the usefulness of MIPhy by accurately predicting which members of the mammalian cytochrome P450 gene superfamily metabolize xenobiotics and which metabolize endogenous compounds. Our predictions correlate very well with known substrate specificities of the human enzymes. We also analyze the Caenorhabditis collagen gene family and use MIPhy to predict genes that produce an observable phenotype when knocked down in C. elegans, and show that our predictions correlate well with existing knowledge. The software can be downloaded and installed from https://github.com/dave-the-scientist/miphy and is also available as an online web tool at http://www.miphy.wasmuthlab.org.

Identification of 74 cytochrome P450 genes and co-localized cytochrome P450 genes of the CYP2K, CYP5A, and CYP46A subfamilies in the mangrove killifish Kryptolebias marmoratus

BACKGROUND

The mangrove killifish Kryptolebias marmoratus is the only vertebrate that reproduces by self-fertilizing and is an important model species in genetics and marine ecotoxicology. Using whole-genome and transcriptome sequences, we identified all members of the cytochrome P450 (CYP) family in this model teleost and compared them with those of other teleosts.

RESULTS

A total of 74 cytochrome P450 genes and one pseudogene were identified in K. marmoratus. Phylogenetic analysis indicated that the CYP genes in clan 2 were most expanded, while synteny analysis with other species showed orthologous relationships of CYP subfamilies among teleosts. In addition to the CYP2K expansions, five tandem duplicated gene copies of CYP5A were observed. These features were unique to K. marmoratus.

CONCLUSIONS

These results shed a light on CYP gene evolution, particularly the co-localized CYP2K, CYP5A, and CYP46A subfamilies in fish. Future studies of CYP expression could identify specific endogenous and exogenous environmental factors that triggered the evolution of tandem CYP duplication in K. marmoratus.

Segmental arithmetic: summing up the Hox gene regulatory network for hindbrain development in chordates

Organization and development of the early vertebrate hindbrain are controlled by a cascade of regulatory interactions that govern the process of segmentation and patterning along the anterior-posterior axis via Hox genes. These interactions can be assembled into a gene regulatory network that provides a framework to interpret experimental data, generate hypotheses, and identify gaps in our understanding of the progressive process of hindbrain segmentation. The network can be broadly separated into a series of interconnected programs that govern early signaling, segmental subdivision, secondary signaling, segmentation, and ultimately specification of segmental identity. Hox genes play crucial roles in multiple programs within this network. Furthermore, the network reveals properties and principles that are likely to be general to other complex developmental systems. Data from vertebrate and invertebrate chordate models are shedding light on the origin and diversification of the network. Comprehensive cis-regulatory analyses of vertebrate Hox gene regulation have enabled powerful cross-species gene regulatory comparisons. Such an approach in the sea lamprey has revealed that the network mediating segmental Hox expression was present in ancestral vertebrates and has been maintained across diverse vertebrate lineages. Invertebrate chordates lack hindbrain segmentation but exhibit conservation of some aspects of the network, such as a role for retinoic acid in establishing nested Hox expression domains. These comparisons lead to a model in which early vertebrates underwent an elaboration of the network between anterior-posterior patterning and Hox gene expression, leading to the gene-regulatory programs for segmental subdivision and rhombomeric segmentation. WIREs Dev Biol 2017, 6:e286. doi: 10.1002/wdev.286 For further resources related to this article, please visit the WIREs website.

New insights into the role of Jmjd3 and Utx in axial skeletal formation in mice

Jmjd3 and Utx are demethylases specific for lysine 27 of histone H3. Previous reports indicate that Jmjd3 is essential for differentiation of various cell types, such as macrophages and epidermal cells in mice, whereas Utx is involved in cancer and developmental diseases in humans and mice, as well as Hox regulation in zebrafish and nematodes. Here, we report that Jmjd3, but not Utx, is involved in axial skeletal formation in mice. A Jmjd3 mutant embryo (Jmjd3^Δ18/Δ18), but not a catalytically inactive Utx truncation mutant (Utx^-/y), showed anterior homeotic transformation. Quantitative RT-PCR and microarray analyses showed reduced Hox expression in both Jmjd3^Δ18/Δ18 embryos and tailbuds, whereas levels of Hox activators, such as Wnt signaling factors and retinoic acid synthases, did not decrease, which suggests that Jmjd3 plays a regulatory role in Hox expression during axial patterning. Chromatin immunoprecipitation analyses of embryo tailbud tissue showed trimethylated lysine 27 on histone H3 to be at higher levels at the Hox loci in Jmjd3^Δ18/Δ18 mutants compared with wild-type tailbuds. In contrast, trimethylated lysine 4 on histone H3 levels were found to be equivalent in wild-type and Jmjd3^Δ18/Δ18 tailbuds. Demethylase-inactive Jmjd3 mutant embryos showed the same phenotype as Jmjd3^Δ18/Δ18 mice. These results suggest that the demethylase activity of Jmjd3, but not that of Utx, affects mouse axial patterning in concert with alterations in Hox gene expression.-Naruse, C., Shibata, S., Tamura, M., Kawaguchi, T., Abe, K., Sugihara, K., Kato, T., Nishiuchi, T., Wakana, S., Ikawa, M., Asano, M. New insights into the role of Jmjd3 and Utx in axial skeletal formation in mice.

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.

Specialization of a polyphenism switch gene following serial duplications in Pristionchus nematodes

Polyphenism is an extreme manifestation of developmental plasticity, requiring distinct developmental programs and the addition of a switch mechanism. Because the genetic basis of polyphenism switches has only begun to be understood, how their mechanisms arise is unclear. In the nematode Pristionchus pacificus, which has a mouthpart polyphenism specialized for alternative diets, a gene (eud-1) executing the polyphenism switch was recently identified as the product of lineage-specific duplications. Here, we infer the role of gene duplications in producing a switch gene. Using reverse genetics and population genetic analyses, we examine evidence for competing scenarios of degeneration and complementation, neutral evolution, and functional specialization. Of the daughter genes, eud-1 alone has assumed switch-like regulation of the mouth polyphenism. Measurements of life-history traits in single, double, and triple sulfatase mutants did not, given a benign environment, identify alternative or complementary roles for eud-1 paralogs. Although possible roles are still unknown, selection analyses of the sister species and 104 natural isolates of P. pacificus detected purifying selection on the genes, suggesting their functionality by their fixation and evolutionary maintenance. Our approach shows the tractability of reverse genetics in a nontraditional model system to study evolution by gene duplication.

Whole Genome Sequencing of the Asian Arowana (Scleropages formosus) Provides Insights into the Evolution of Ray-Finned Fishes

The Asian arowana (Scleropages formosus) is of commercial importance, conservation concern, and is a representative of one of the oldest lineages of ray-finned fish, the Osteoglossomorpha. To add to genomic knowledge of this species and the evolution of teleosts, the genome of a Malaysian specimen of arowana was sequenced. A draft genome is presented consisting of 42,110 scaffolds with a total size of 708 Mb (2.85% gaps) representing 93.95% of core eukaryotic genes. Using a k-mer-based method, a genome size of 900 Mb was also estimated. We present an update on the phylogenomics of fishes based on a total of 27 species (23 fish species and 4 tetrapods) using 177 orthologous proteins (71,360 amino acid sites), which supports established relationships except that arowana is placed as the sister lineage to all teleost clades (Bayesian posterior probability 1.00, bootstrap replicate 93%), that evolved after the teleost genome duplication event rather than the eels (Elopomorpha). Evolutionary rates are highly heterogeneous across the tree with fishes represented by both slowly and rapidly evolving lineages. A total of 94 putative pigment genes were identified, providing the impetus for development of molecular markers associated with the spectacular colored phenotypes found within this species.

Roles of Hox genes in the patterning of the central nervous system of Drosophila

One of the key aspects of functional nervous systems is the restriction of particular neural subtypes to specific regions, which permits the establishment of differential segment-specific neuromuscular networks. Although Hox genes play a major role in shaping the anterior-posterior body axis during animal development, our understanding of how they act in individual cells to determine particular traits at precise developmental stages is rudimentary. We have used the abdominal leucokinergic neurons (ABLKs) to address this issue. These neurons are generated during both embryonic and postembryonic neurogenesis by the same progenitor neuroblast, and are designated embryonic and postembryonic ABLKs, respectively. We report that the genes of the Bithorax-Complex, Ultrabithorax (Ubx) and abdominal-A (abd-A) are redundantly required to specify the embryonic ABLKs. Moreover, the segment-specific pattern of the postembryonic ABLKs, which are restricted to the most anterior abdominal segments, is controlled by the absence of Abdominal-B (Abd-B), which we found was able to repress the expression of the neuropeptide leucokinin. We discuss this and other examples of how Hox genes generate diversity within the central nervous system of Drosophila.

The regulation of Hox gene expression during animal development

Hox genes encode a family of transcriptional regulators that elicit distinct developmental programmes along the head-to-tail axis of animals. The specific regional functions of individual Hox genes largely reflect their restricted expression patterns, the disruption of which can lead to developmental defects and disease. Here, we examine the spectrum of molecular mechanisms controlling Hox gene expression in model vertebrates and invertebrates and find that a diverse range of mechanisms, including nuclear dynamics, RNA processing, microRNA and translational regulation, all concur to control Hox gene outputs. We propose that this complex multi-tiered regulation might contribute to the robustness of Hox expression during development.

Genome-wide identification, evolutionary and expression analysis of the aspartic protease gene superfamily in grape

Background

Aspartic proteases (APs) are a large family of proteolytic enzymes found in almost all organisms. In plants, they are involved in many biological processes, such as senescence, stress responses, programmed cell death, and reproduction. Prior to the present study, no grape AP gene(s) had been reported, and their research on woody species was very limited.

Results

In this study, a total of 50 AP genes (VvAP) were identified in the grape genome, among which 30 contained the complete ASP domain. Synteny analysis within grape indicated that segmental and tandem duplication events contributed to the expansion of the grape AP family. Additional analysis between grape and Arabidopsis demonstrated that several grape AP genes were found in the corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the divergence of grape and Arabidopsis. Phylogenetic relationships of the 30 VvAPs with the complete ASP domain and their Arabidopsis orthologs, as well as their gene and protein features were analyzed and their cellular localization was predicted. Moreover, expression profiles of VvAP genes in six different tissues were determined, and their transcript abundance under various stresses and hormone treatments were measured. Twenty-seven VvAP genes were expressed in at least one of the six tissues examined; nineteen VvAPs responded to at least one abiotic stress, 12 VvAPs responded to powdery mildew infection, and most of the VvAPs responded to SA and ABA treatments. Furthermore, integrated synteny and phylogenetic analysis identified orthologous AP genes between grape and Arabidopsis, providing a unique starting point for investigating the function of grape AP genes.

Conclusions

The genome-wide identification, evolutionary and expression analyses of grape AP genes provide a framework for future analysis of AP genes in defining their roles during stress response. Integrated synteny and phylogenetic analyses provide novel insight into the functions of less well-studied genes using information from their better understood orthologs.

Lost in the map

Organismal development and evolution are complex, multifaceted processes that depend intimately on context. They are subject to environmental influences, chance appearance and fixation of mutations, and numerous other idiosyncrasies. Genomics is detailing the molecular signature of effects of these mechanisms on phenotypes, but because numerous distinct evolutionary explanations can produce a given genomic pattern, the molecular details, rather than elucidating process, typically distract from explanatory insight and contribute little to predictive capability. While genomic research has burgeoned, direct study of evolutionary and developmental processes has lagged. We advocate for reinvigoration of direct study of process, along with refocusing of attention on questions of broad biological import, as more productive of urgently needed insights, which genomic approaches are not providing.

Genome-wide identification and analysis of the TIFY gene family in grape

BACKGROUND

The TIFY gene family constitutes a plant-specific group of genes with a broad range of functions. This family encodes four subfamilies of proteins, including ZML, TIFY, PPD and JASMONATE ZIM-Domain (JAZ) proteins. JAZ proteins are targets of the SCF(COI1) complex, and function as negative regulators in the JA signaling pathway. Recently, it has been reported in both Arabidopsis and rice that TIFY genes, and especially JAZ genes, may be involved in plant defense against insect feeding, wounding, pathogens and abiotic stresses. Nonetheless, knowledge concerning the specific expression patterns and evolutionary history of plant TIFY family members is limited, especially in a woody species such as grape.

METHODOLOGY/PRINCIPAL FINDINGS

A total of two TIFY, four ZML, two PPD and 11 JAZ genes were identified in the Vitis vinifera genome. Phylogenetic analysis of TIFY protein sequences from grape, Arabidopsis and rice indicated that the grape TIFY proteins are more closely related to those of Arabidopsis than those of rice. Both segmental and tandem duplication events have been major contributors to the expansion of the grape TIFY family. In addition, synteny analysis between grape and Arabidopsis demonstrated that homologues of several grape TIFY genes were found in the corresponding syntenic blocks of Arabidopsis, suggesting that these genes arose before the divergence of lineages that led to grape and Arabidopsis. Analyses of microarray and quantitative real-time RT-PCR expression data revealed that grape TIFY genes are not a major player in the defense against biotrophic pathogens or viruses. However, many of these genes were responsive to JA and ABA, but not SA or ET.

CONCLUSION

The genome-wide identification, evolutionary and expression analyses of grape TIFY genes should facilitate further research of this gene family and provide new insights regarding their evolutionary history and regulatory control.

Integration of the genetic map and genome assembly of fugu facilitates insights into distinct features of genome evolution in teleosts and mammals

The compact genome of fugu (Takifugu rubripes) has been used widely as a reference genome for understanding the evolution of vertebrate genomes. However, the fragmented nature of the fugu genome assembly has restricted its use for comparisons of genome architecture in vertebrates. To extend the contiguity of the assembly to the chromosomal level, we have generated a comprehensive genetic map of fugu and anchored the scaffolds of the assembly to the 22 chromosomes of fugu. The map consists of 1,220 microsatellite markers that provide anchor points to 697 scaffolds covering 86% of the genome assembly (http://www.fugu-sg.org/). The integrated genome map revealed a higher recombination rate in fugu compared with other vertebrates and a wide variation in the recombination rate between sexes and across chromosomes of fugu. We used the extended assembly to explore recent rearrangement events in the lineages of fugu, Tetraodon, and medaka and compared them with rearrangements in three mammalian (human, mouse, and opossum) lineages. Between the two pufferfishes, fugu has experienced fewer chromosomal rearrangements than Tetraodon. The gene order is more highly conserved in the three teleosts than in mammals largely due to a lower rate of interchromosomal rearrangements in the teleosts. These results provide new insights into the distinct patterns of genome evolution between teleosts and mammals. The consolidated genome map and the genetic map of fugu are valuable resources for comparative genomics of vertebrates and for elucidating the genetic basis of the phenotypic diversity of ~25 species of Takifugu that evolved within the last 5 My.

Prenatal development of central rhythm generation

Foetal breathing in mice results from prenatal activity of the two coupled hindbrain oscillators considered to be responsible for respiratory rhythm generation after birth: the pre-Bötzinger complex (preBötC) is active shortly before the onset of foetal breathing; the parafacial respiratory group (e-pF in embryo) starts activity one day earlier. Transcription factors have been identified that are essential to specify neural progenitors and lineages forming each of these oscillators during early development of the neural tube: Hoxa1, Egr2 (Krox20), Phox2b, Lbx1 and Atoh1 for the e-pF; Dbx1 and Evx1 for the preBötC which eventually grow contralateral axons requiring expression of Robo3. Inactivation of the genes encoding these factors leads to mis-specification of these neurons and distinct breathing abnormalities: apneic patterns and loss of central chemosensitivity for the e-pF (central congenital hypoventilation syndrome, CCHS, in humans), complete loss of breathing for the preBötC, right-left desynchronized breathing in Robo3 mutants. Mutations affecting development in more rostral (pontine) respiratory territories change the shape of the inspiratory drive without affecting the rhythm. Other (primordial) embryonic oscillators start in the mouse three days before the e-pF, to generate low frequency (LF) rhythms that are probably required for activity-dependent development of neurones at embryonic stages; in the foetus, however, they are actively silenced to avoid detrimental interaction with the on-going respiratory rhythm. Altogether, these observations provide a strong support to the previously proposed hypothesis that the functional organization of the respiratory generator is specified at early stages of development and is dual in nature, comprising two serially non-homologous oscillators.

Differential selection within the Drosophila retinal determination network and evidence for functional divergence between paralog pairs

The retinal determination (RD) network in Drosophila comprises 14 known nuclear proteins that include DNA-binding proteins, transcriptional coactivators, kinases, and phosphatases. The composition of the network varies considerably throughout the animal kingdom, with the network in several basal insects having fewer members and with vertebrates having potentially significantly higher numbers of RD genes. One important contributing factor for the variation in gene number within the network is gene duplication. For example, 10 members of the RD network in Drosophila are derived from duplication events. Here we present an analysis of the coding regions of the five pairs of duplicate genes from within the RD network of several different Drosophila species. We demonstrate that there is differential selection across the coding regions of all RD genes. Additionally, some of the most significant differences in ratios of non-silent-to-silent site substitutions (d(N)/d(S)) between paralog pairs are found within regions that have no ascribed function. Previous structure/function analyses of several duplicate genes have identified areas within one gene that contain novel activities when compared with its paralog. The evolutionary analysis presented here identifies these same areas in the paralogs as being under high levels of relaxed selection. We suggest that sequence divergence between paralogs and selection signatures can be used as a reasonable predictor of functional changes in rapidly evolving motifs.

Genetic factors determining the functional organization of neural circuits controlling rhythmic movements the murine embryonic parafacial rhythm generator

In mammals, fetal movements governed by central pattern generators are essential for the development of adaptive behaviors such as breathing, walking, and chewing, which are vital after birth. Combining targeted mutations and genetic fate mapping can help to define the molecular determinants that control the development of these central pattern generators. In this chapter, recent results are presented on the embryonic parafacial (e-pF) rhythm generator, one of the two oscillators involved in controlling the breathing behavior and chemosensitive responsiveness.

Anion exchanger 1b, but not sodium-bicarbonate cotransporter 1b, plays a role in transport functions of zebrafish H+-ATPase-rich cells

Similar to mammalian proximal tubular cells, H(+)-ATPase rich (HR) cells in zebrafish skin and gills are also responsible for Na(+) uptake and acid secretion functions. However, the basolateral transport pathways in HR cells are still unclear. In the present study, we tested the hypothesis if there are specific slc4 members involved in basolateral ion transport pathways in HR cells. Fourteen isoforms were identified in the zebrafish(z) slc4 family, and the full-length cDNAs of two novel isoforms, zslc4a1b (anion exchanger, zAE1b) and zslc4a4b (Na(+)/HCO(3)(-) cotransporter, zNBCe1b), were sequenced. mRNA signals of zslc4a1b and zslc4a4b were mainly detected in certain groups of ionocytes in zebrafish skin/gills. Further double immunocytochemistry or in situ hybridization demonstrated that zAE1b, but not zNBCe1b, was localized to basolateral membranes of HR cells. Acclimation to low-Na(+) or acidic environments stimulated the mRNA expression of zslc4a1b in zebrafish gills, and loss-of-function of zslc4a1b with specific morpholinos caused significant decreases in both the whole body Na(+) content and the skin H(+) activity in the morphants. On the basis of these results, it was concluded that zAE1b, but not zNBCe1b, is involved in the basolateral transport pathways in Na(+) uptake/acid secretion mechanisms in zebrafish HR cells.

Analyses of nuclearly encoded mitochondrial genes suggest gene duplication as a mechanism for resolving intralocus sexually antagonistic conflict in Drosophila

Gene duplication is probably the most important mechanism for generating new gene functions. However, gene duplication has been overlooked as a potentially effective way to resolve genetic conflicts. Here, we analyze the entire set of Drosophila melanogaster nuclearly encoded mitochondrial duplicate genes and show that both RNA- and DNA-mediated mitochondrial gene duplications exhibit an unexpectedly high rate of relocation (change in location between parental and duplicated gene) as well as an extreme tendency to avoid the X chromosome. These trends are likely related to our observation that relocated genes tend to have testis-specific expression. We also infer that these trends hold across the entire Drosophila genus. Importantly, analyses of gene ontology and functional interaction networks show that there is an overrepresentation of energy production-related functions in these mitochondrial duplicates. We discuss different hypotheses to explain our results and conclude that our findings substantiate the hypothesis that gene duplication for male germline function is likely a mechanism to resolve intralocus sexually antagonistic conflicts that we propose are common in testis. In the case of nuclearly encoded mitochondrial duplicates, our hypothesis is that past sexually antagonistic conflict related to mitochondrial energy function in Drosophila was resolved by gene duplication.

Cloning and characterization of a zebrafish homologue of human AQP1: a bifunctional water and gas channel

The mammalian aquaporins AQP1, AQP4, and AQP5 have been shown to function not only as water channels but also as gas channels. Zebrafish have two genes encoding an AQP1 homologue, aqp1a and aqp1b. In the present study, we cloned the cDNA that encodes the zebrafish protein Aqp1a from the 72-h postfertilization (hpf) embryo of Danio rerio, as well as from the swim bladder of the adult. The deduced amino-acid sequence of aqp1a consists of 260 amino acids and is 59% identical to human AQP1. By analyzing the genomic DNA sequence, we identified four exons in the aqp1a gene. By in situ hybridization, aqp1a is expressed transiently in the developing vasculature and in erythrocytes from 16 to 48 h of development. Later, at 72 hpf, aqp1a is expressed in dermal ionocytes and in the swim bladder. Western blot analysis of adult tissues reveals that Aqp1a is most highly expressed in the eye and swim bladder. Xenopus oocytes expressing aqp1a have a channel-dependent (*) osmotic water permeability (P(f)(*)) that is indistinguishable from that of human AQP1. On the basis of the magnitude of the transient change in surface pH (ΔpH(S)) that were recorded as the oocytes were exposed to either CO(2) or NH(3), we conclude that zebrafish Aqp1a is permeable to both CO(2) and NH(3). The ratio (ΔpH(S)(*))((CO)2)/P(f)(*) is about half that of human AQP1, and the ratio (ΔpH(S)(*))(NH3)/P(f)(*) is about one-quarter that of human AQP1. Thus, compared with human AQP1, zebrafish Aqp1a has about twice the selectivity for CO(2) over NH(3).

The birth of new genes by RNA- and DNA-mediated duplication during mammalian evolution

Gene duplication has long been recognized as a major force in genome evolution and has recently been recognized as an important source of individual variation. For many years, the origin of functional gene duplicates was assumed to be whole or partial genome duplication events, but recently retrotransposition has also been shown to contribute new functional protein coding genes and siRNA's. In this study, we utilize pseudogenes to recreate more complete gene family histories, and compare the rates of RNA and DNA-mediated duplication and new functional gene formation in five mammalian genomes. We find that RNA-mediated duplication occurs at a much higher and more variable rate than DNA-mediated duplication, and gives rise to many more duplicated sequences over time. We show that, while the chance of RNA-mediated duplicates becoming functional is much lower than that of their DNA-mediated counterparts, the higher rate of retrotransposition leads to nearly equal contributions of new genes by each mechanism. We also find that functional RNA-mediated duplicates are closer to neighboring genes than non-functional RNA-mediated copies, consistent with co-option of regulatory elements at the site of insertion. Overall, new genes derived from DNA and RNA-mediated duplication mechanisms are under similar levels of purifying selective pressure, but have broadly different functions. RNA-mediated duplication gives rise to a diversity of genes but is dominated by the highly expressed genes of RNA metabolic pathways. DNA-mediated duplication can copy regulatory material along with the protein coding region of the gene and often gives rise to classes of genes whose function are dependent on complex regulatory information. This mechanistic difference may in part explain why we find that mammalian protein families tend to evolve by either one mechanism or the other, but rarely by both. Supplementary Material has been provided (see online Supplementary Material at www.liebertonline.com ).

Gene duplication, transfer, and evolution in the chloroplast genome

In addition to the nuclear genome, organisms have organelle genomes. Most of the DNA present in eukaryotic organisms is located in the cell nucleus. Chloroplasts have independent genomes which are inherited from the mother. Duplicated genes are common in the genomes of all organisms. It is believed that gene duplication is the most important step for the origin of genetic variation, leading to the creation of new genes and new gene functions. Despite the fact that extensive gene duplications are rare among the chloroplast genome, gene duplication in the chloroplast genome is an essential source of new genetic functions and a mechanism of neo-evolution. The events of gene transfer between the chloroplast genome and nuclear genome via duplication and subsequent recombination are important processes in evolution. The duplicated gene or genome in the nucleus has been the subject of several recent reviews. In this review, we will briefly summarize gene duplication and evolution in the chloroplast genome. Also, we will provide an overview of gene transfer events between chloroplast and nuclear genomes.

Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Background

The 12-oxo-phytodienoic acid reductases (OPRs) are enzymes that catalyze the reduction of double-bonds in α, β-unsaturated aldehydes or ketones and are part of the octadecanoid pathway that converts linolenic acid to jasmonic acid. In plants, OPRs belong to the old yellow enzyme family and form multigene families. Although discoveries about this family in Arabidopsis and other species have been reported in some studies, the evolution and function of multiple OPRs in plants are not clearly understood.

Results

A comparative genomic analysis was performed to investigate the phylogenetic relationship, structural evolution and functional divergence among OPR paralogues in plants. In total, 74 OPR genes were identified from 11 species representing the 6 major green plant lineages: green algae, mosses, lycophytes, gymnosperms, monocots and dicots. Phylogenetic analysis showed that seven well-conserved subfamilies exist in plants. All OPR genes from green algae were clustered into a single subfamily, while those from land plants fell into six other subfamilies, suggesting that the events leading to the expansion of the OPR family occurred in land plants. Further analysis revealed that lineage-specific expansion, especially by tandem duplication, contributed to the current OPR subfamilies in land plants after divergence from aquatic plants. Interestingly, exon/intron structure analysis showed that the gene structures of OPR paralogues exhibits diversity in intron number and length, while the intron positions and phase were highly conserved across different lineage species. These observations together with the phylogenetic tree revealed that successive single intron loss, as well as indels within introns, occurred during the process of structural evolution of OPR paralogues. Functional divergence analysis revealed that altered functional constraints have occurred at specific amino acid positions after diversification of the paralogues. Most notably, significant functional divergence was also found in all pairs, except for the II/IV, II/V and V/VI pairs. Strikingly, analysis of the site-specific profiles established by posterior probability revealed that the positive-selection sites and/or critical amino acid residues for functional divergence are mainly distributed in α-helices and substrate binding loop (SBL), indicating the functional importance of these regions for this protein family.

Conclusion

This study highlights the molecular evolution of the OPR gene family in all plant lineages and indicates critical amino acid residues likely relevant for the distinct functional properties of the paralogues. Further experimental verification of these findings may provide valuable information on the OPRs' biochemical and physiological functions.

Heterologous expression of human mPRalpha, mPRbeta and mPRgamma in yeast confirms their ability to function as membrane progesterone receptors

The nuclear progesterone receptor (nPR) mediates many of the physiological effects of progesterone by regulating the expression of genes, however, progesterone also exerts non-transcriptional (non-genomic) effects that have been proposed to rely on a receptor that is distinct from nPR. Several members of the progestin and AdipoQ-Receptor (PAQR) family were recently identified as potential mediators of these non-genomic effects. Membranes from cells expressing these proteins, called mPRalpha, mPRbeta and mPRgamma, were shown to specifically bind progesterone and have G-protein coupled receptor (GPCR) characteristics, although other studies dispute these findings. To clarify the role of these mPRs in non-genomic progesterone signaling, we established an assay for PAQR functional evaluation using heterologous expression in Saccharomyces cerevisiae. Using this assay, we demonstrate unequivocally that mPRalpha, mPRbeta and mPRgamma can sense and respond to progesterone with EC(50) values that are physiologically relevant. Agonist profiles also show that mPRalpha, mPRbeta and mPRgamma are activated by ligands, such as 17alpha-hydroxyprogesterone, that are known to activate non-genomic pathways but not nPR. These results strongly suggest that these receptors may indeed function as the long-sought-after membrane progesterone receptors. Additionally, we show that two uncharacterized PAQRs, PAQR6 and PAQR9, are also capable of responding to progesterone. These mPR-like PAQRs have been renamed as mPRdelta (PAQR6) and mPRvarepsilon (PAQR9). Additional characterization of mPRgamma and mPRalpha indicates that their progesterone-dependent signaling in yeast does not require heterotrimeric G-proteins, thus calling into question the characterization of the mPRs as a novel class of G-protein coupled receptor.

The loci of evolution: how predictable is genetic evolution?

Is genetic evolution predictable? Evolutionary developmental biologists have argued that, at least for morphological traits, the answer is a resounding yes. Most mutations causing morphological variation are expected to reside in the cis-regulatory, rather than the coding, regions of developmental genes. This "cis-regulatory hypothesis" has recently come under attack. In this review, we first describe and critique the arguments that have been proposed in support of the cis-regulatory hypothesis. We then test the empirical support for the cis-regulatory hypothesis with a comprehensive survey of mutations responsible for phenotypic evolution in multicellular organisms. Cis-regulatory mutations currently represent approximately 22% of 331 identified genetic changes although the number of cis-regulatory changes published annually is rapidly increasing. Above the species level, cis-regulatory mutations altering morphology are more common than coding changes. Also, above the species level cis-regulatory mutations predominate for genes not involved in terminal differentiation. These patterns imply that the simple question "Do coding or cis-regulatory mutations cause more phenotypic evolution?" hides more interesting phenomena. Evolution in different kinds of populations and over different durations may result in selection of different kinds of mutations. Predicting the genetic basis of evolution requires a comprehensive synthesis of molecular developmental biology and population genetics.

Gene duplication and transfer events in plant mitochondria genome

Gene or genome duplication events increase the amount of genetic material available to increase the genomic, and thereby phenotypic, complexity of organisms during evolution. Gene duplication and transfer events have been important to molecular evolution in all three domains of life, and may be the first step in the emergence of new gene functions. Gene transfer events have been proposed as another accelerator of evolution. The duplicated gene or genome, mainly nuclear, has been the subject of several recent reviews. In addition to the nuclear genome, organisms have organelle genomes, including mitochondrial genome. In this review, we briefly summarize gene duplication and transfer events in the plant mitochondrial genome.

Subfunctionalization of duplicated zebrafish pax6 genes by cis-regulatory divergence

Gene duplication is a major driver of evolutionary divergence. In most vertebrates a single PAX6 gene encodes a transcription factor required for eye, brain, olfactory system, and pancreas development. In zebrafish, following a postulated whole-genome duplication event in an ancestral teleost, duplicates pax6a and pax6b jointly fulfill these roles. Mapping of the homozygously viable eye mutant sunrise identified a homeodomain missense change in pax6b, leading to loss of target binding. The mild phenotype emphasizes role-sharing between the co-orthologues. Meticulous mapping of isolated BACs identified perturbed synteny relationships around the duplicates. This highlights the functional conservation of pax6 downstream (3') control sequences, which in most vertebrates reside within the introns of a ubiquitously expressed neighbour gene, ELP4, whose pax6a-linked exons have been lost in zebrafish. Reporter transgenic studies in both mouse and zebrafish, combined with analysis of vertebrate sequence conservation, reveal loss and retention of specific cis-regulatory elements, correlating strongly with the diverged expression of co-orthologues, and providing clear evidence for evolution by subfunctionalization.

Consequences of hoxb1 duplication in teleost fish

Vertebrate evolution is characterized by gene and genome duplication events. There is strong evidence that a whole-genome duplication occurred in the lineage leading to the teleost fishes. We have focused on the teleost hoxb1 duplicate genes as a paradigm to investigate the consequences of gene duplication. Previous analysis of the duplicated zebrafish hoxb1 genes suggested they have subfunctionalized. The combined expression pattern of the two zebrafish hoxb1 genes recapitulates the expression pattern of the single Hoxb1 gene of tetrapods, possibly due to degenerative changes in complementary cis-regulatory elements of the duplicates. Here we have tested the hypothesis that all teleost duplicates had a similar fate post duplication, by examining hoxb1 genes in medaka and striped bass. Consistent with this theory, we found that the ancestral Hoxb1 expression pattern is subdivided between duplicate genes in a largely similar fashion in zebrafish, medaka, and striped bass. Further, our analysis of hoxb1 genes reveals that sequence changes in cis-regulatory regions may underlie subfunctionalization in all teleosts, although the specific changes vary between species. It was previously shown that zebrafish hoxb1 duplicates have also evolved different functional capacities. We used misexpression to compare the functions of hoxb1 duplicates from zebrafish, medaka and striped bass. Unexpectedly, we found that some biochemical properties, which were paralog specific in zebrafish, are conserved in both duplicates of other species. This work suggests that the fate of duplicate genes varies across the teleost group.

Ancestral and recently recruited global control of the Hox genes in development

Genes from the Hox family are involved in the common task of providing nascent embryonic tissues with their positional identity. They are organised in clusters in most species. Mouse Hox genes are regulated in part by gene-proximal regulatory elements, but owe several of their essential properties to the use of global regulatory elements located outside the complexes. The clustered Hox genes in that sense behave as a single large locus. Genomic and sequence data from different animal species suggest that a concerted regulation of the Hox clusters, inherently coupled to their patterning properties, originated early during evolution and pre-figured the temporal colinearity of expression of vertebrate Hox genes. In addition, vertebrates have recruited novel global mechanisms to control the expression of linear subsets of Hox genes in specific embryonic structures. Several of such novel global regulatory circuits have recently been characterised at the molecular genetic level in the mouse.

PACAP-related peptide (PRP)--molecular evolution and potential functions

PACAP-related peptide (PRP) and PACAP are structurally related peptides that are encoded in the same transcripts. In the past, it was believed that the mammalian PRPs are evolved from GHRHs in non-mammals. With the recent discovery of authentic GHRH and receptor genes in frog and fish, this review aims to (1) coin the name of all GHRH-like peptides in previous literature as PRPs and (2) provide the background for new research direction for PRP in vertebrates. As a goldfish receptor highly specific for PRP with distinct tissue distribution has previously been characterized, it is highly possible that PRP plays a physiological role in non-mammalian vertebrates and the function of PRP has somehow been lost in mammals as a consequence of the loss of its receptor in the genome. This information may provide clues to elucidate functions of PRP in the future.

Phylogenomic analyses of KCNA gene clusters in vertebrates: why do gene clusters stay intact?

Background

Gene clusters are of interest for the understanding of genome evolution since they provide insight in large-scale duplications events as well as patterns of individual gene losses. Vertebrates tend to have multiple copies of gene clusters that typically are only single clusters or are not present at all in genomes of invertebrates. We investigated the genomic architecture and conserved non-coding sequences of vertebrate KCNA gene clusters. KCNA genes encode shaker-related voltage-gated potassium channels and are arranged in two three-gene clusters in tetrapods. Teleost fish are found to possess four clusters. The two tetrapod KNCA clusters are of approximately the same age as the Hox gene clusters that arose through duplications early in vertebrate evolution. For some genes, their conserved retention and arrangement in clusters are thought to be related to regulatory elements in the intergenic regions, which might prevent rearrangements and gene loss. Interestingly, this hypothesis does not appear to apply to the KCNA clusters, as too few conserved putative regulatory elements are retained.

Results

We obtained KCNA coding sequences from basal ray-finned fishes (sturgeon, gar, bowfin) and confirmed that the duplication of these genes is specific to teleosts and therefore consistent with the fish-specific genome duplication (FSGD). Phylogenetic analyses of the genes suggest a basal position of the only intron containing KCNA gene in vertebrates (KCNA7). Sistergroup relationships of KCNA1/2 and KCNA3/6 support that a large-scale duplication gave rise to the two clusters found in the genome of tetrapods. We analyzed the intergenic regions of KCNA clusters in vertebrates and found that there are only a few conserved sequences shared between tetrapods and teleosts or between paralogous clusters. The orthologous teleost clusters, however, show sequence conservation in these regions.

Conclusion

The lack of overall conserved sequences in intergenic regions suggests that there are either other processes than regulatory evolution leading to cluster conservation or that the ancestral regulatory relationships among genes in KCNA clusters have been changed together with their regulatory sites.

Polyglutamine variation in a flowering time protein correlates with island age in a Hawaiian plant radiation

Background

A controversial topic in evolutionary developmental biology is whether morphological diversification in natural populations can be driven by expansions and contractions of amino acid repeats in proteins. To promote adaptation, selection on protein length variation must overcome deleterious effects of multiple correlated traits (pleiotropy). Thus far, systems that demonstrate this capacity include only ancient or artificial morphological diversifications. The Hawaiian Islands, with their linear geological sequence, present a unique environment to study recent, natural radiations. We have focused our research on the Hawaiian endemic mints (Lamiaceae), a large and diverse lineage with paradoxically low genetic variation, in order to test whether a direct relationship between coding-sequence repeat diversity and morphological change can be observed in an actively evolving system.

Results

Here we show that in the Hawaiian mints, extensive polyglutamine (CAG codon repeat) polymorphism within a homolog of the pleiotropic flowering time protein and abscisic acid receptor FCA tracks the natural environmental cline of the island chain, consequent with island age, across a period of 5 million years. CAG expansions, perhaps following their natural tendency to elongate, are more frequent in colonists of recently-formed, nutrient-rich islands than in their forebears on older, nutrient-poor islands. Values for several quantitative morphological variables related to reproductive investment, known from Arabidopsis fca mutant studies, weakly though positively correlate with increasing glutamine tract length. Together with protein modeling of FCA, which indicates that longer polyglutamine tracts could induce suboptimally mobile functional domains, we suggest that CAG expansions may form slightly deleterious alleles (with respect to protein function) that become fixed in founder populations.

Conclusion

In the Hawaiian mint FCA system, we infer that contraction of slightly deleterious CAG repeats occurred because of competition for resources along the natural environmental cline of the island chain. The observed geographical structure of FCA variation and its correlation with morphologies expected from Arabidopsis mutant studies may indicate that developmental pleiotropy played a role in the diversification of the mints. This discovery is important in that it concurs with other suggestions that repetitive amino acid motifs might provide a mechanism for driving morphological evolution, and that variation at such motifs might permit rapid tuning to environmental change.

A new time-scale for ray-finned fish evolution

The Actinopterygii (ray-finned fishes) is the largest and most diverse vertebrate group, but little is agreed about the timing of its early evolution. Estimates using mitochondrial genomic data suggest that the major actinopterygian clades are much older than divergence dates implied by fossils. Here, the timing of the evolutionary origins of these clades is reinvestigated using morphological, and nuclear and mitochondrial genetic data. Results indicate that existing fossil-based estimates of the age of the crown-group Neopterygii, including the teleosts, Lepisosteus (gar) and Amia (bowfin), are at least 40 Myr too young. We present new palaeontological evidence that the neopterygian crown radiation is a Palaeozoic event, and demonstrate that conflicts between molecular and morphological data for the age of the Neopterygii result, in part, from missing fossil data. Although our molecular data also provide an older age estimate for the teleost crown, this range extension remains unsupported by the fossil evidence. Nuclear data from all relevant clades are used to demonstrate that the actinopterygian whole-genome duplication event is teleost-specific. While the date estimate of this event overlaps the probable range of the teleost stem group, a correlation between the genome duplication and the large-scale pattern of actinopterygian phylogeny remains elusive.

The rise and fall of Hox gene clusters

Although all bilaterian animals have a related set of Hox genes, the genomic organization of this gene complement comes in different flavors. In some unrelated species, Hox genes are clustered; in others, they are not. This indicates that the bilaterian ancestor had a clustered Hox gene family and that, subsequently, this genomic organization was either maintained or lost. Remarkably, the tightest organization is found in vertebrates, raising the embarrassingly finalistic possibility that vertebrates have maintained best this ancestral configuration. Alternatively, could they have co-evolved with an increased ;organization' of the Hox clusters, possibly linked to their genomic amplification, which would be at odds with our current perception of evolutionary mechanisms? When discussing the why's and how's of Hox gene clustering, we need to account for three points: the mechanisms of cluster evolution; the underlying biological constraints; and the developmental modes of the animals under consideration. By integrating these parameters, general conclusions emerge that can help solve the aforementioned dilemma.

Viable mice with compound mutations in the Wnt/Dvl pathway antagonists nkd1 and nkd2

Gradients of Wnt/beta-catenin signaling coordinate development and physiological homeostasis in metazoan animals. Proper embryonic development of the fruit fly Drosophila melanogaster requires the Naked cuticle (Nkd) protein to attenuate a gradient of Wnt/beta-catenin signaling across each segmental anlage. Nkd inhibits Wnt signaling by binding the intracellular protein Dishevelled (Dsh). Mice and humans have two nkd homologs, nkd1 and nkd2, whose encoded proteins can bind Dsh homologs (the Dvl proteins) and inhibit Wnt signaling. To determine whether nkd genes are necessary for murine development, we replaced nkd exons that encode Dvl-binding sequences with IRES-lacZ/neomycin cassettes. Mutants homozygous for each nkd(lacZ) allele are viable with slightly reduced mean litter sizes. Surprisingly, double-knockout mice are viable, with subtle alterations in cranial bone morphology that are reminiscent of mutation in another Wnt/beta-catenin antagonist, axin2. Our data show that nkd function in the mouse is dispensable for embryonic development.

Rapid birth-death evolution specific to xenobiotic cytochrome P450 genes in vertebrates

Genes vary greatly in their long-term phylogenetic stability and there exists no general explanation for these differences. The cytochrome P450 (CYP450) gene superfamily is well suited to investigating this problem because it is large and well studied, and it includes both stable and unstable genes. CYP450 genes encode oxidase enzymes that function in metabolism of endogenous small molecules and in detoxification of xenobiotic compounds. Both types of enzymes have been intensively studied. My analysis of ten nearly complete vertebrate genomes indicates that each genome contains 50-80 CYP450 genes, which are about evenly divided between phylogenetically stable and unstable genes. The stable genes are characterized by few or no gene duplications or losses in species ranging from bony fish to mammals, whereas unstable genes are characterized by frequent gene duplications and losses (birth-death evolution) even among closely related species. All of the CYP450 genes that encode enzymes with known endogenous substrates are phylogenetically stable. In contrast, most of the unstable genes encode enzymes that function as xenobiotic detoxifiers. Nearly all unstable CYP450 genes in the mouse and human genomes reside in a few dense gene clusters, forming unstable gene islands that arose by recurrent local gene duplication. Evidence for positive selection in amino acid sequence is restricted to these unstable CYP450 genes, and sites of selection are associated with substrate-binding regions in the protein structure. These results can be explained by a general model in which phylogenetically stable genes have core functions in development and physiology, whereas unstable genes have accessory functions associated with unstable environmental interactions such as toxin and pathogen exposure. Unstable gene islands in vertebrates share some functional properties with bacterial genomic islands, though they arise by local gene duplication rather than horizontal gene transfer.

Vertebrate vitellogenin gene duplication in relation to the "3R hypothesis": correlation to the pelagic egg and the oceanic radiation of teleosts

The spiny ray-finned teleost fishes (Acanthomorpha) are the most successful group of vertebrates in terms of species diversity. Their meteoric radiation and speciation in the oceans during the late Cretaceous and Eocene epoch is unprecedented in vertebrate history, occurring in one third of the time for similar diversity to appear in the birds and mammals. The success of marine teleosts is even more remarkable considering their long freshwater ancestry, since it implies solving major physiological challenges when freely broadcasting their eggs in the hyper-osmotic conditions of seawater. Most extant marine teleosts spawn highly hydrated pelagic eggs, due to differential proteolysis of vitellogenin (Vtg)-derived yolk proteins. The maturational degradation of Vtg involves depolymerization of mainly the lipovitellin heavy chain (LvH) of one form of Vtg to generate a large pool of free amino acids (FAA 150-200 mM). This organic osmolyte pool drives hydration of the ooctye while still protected within the maternal ovary. In the present contribution, we have used Bayesian analysis to examine the evolution of vertebrate Vtg genes in relation to the "3R hypothesis" of whole genome duplication (WGD) and the functional end points of LvH degradation during oocyte maturation. We find that teleost Vtgs have experienced a post-R3 lineage-specific gene duplication to form paralogous clusters that correlate to the pelagic and benthic character of the eggs. Neo-functionalization allowed one paralogue to be proteolyzed to FAA driving hydration of the maturing oocytes, which pre-adapts them to the marine environment and causes them to float. The timing of these events matches the appearance of the Acanthomorpha in the fossil record. We discuss the significance of these adaptations in relation to ancestral physiological features, and propose that the neo-functionalization of duplicated Vtg genes was a key event in the evolution and success of the teleosts in the oceanic environment.

Phylogenetic relationships and gene expression pattern of three different cathepsin L (Ctsl) isoforms in zebrafish: Ctsla is the putative yolk processing enzyme

Certain cysteine proteases, such as cathepsin L (Ctsl), have been involved in yolk processing mechanisms in oocytes and embryos of lower vertebrates. In zebrafish (Danio rerio), three different ctsl genes, ctsla, ctslb and ctslc, have been found in the genome, but their pattern of expression, as well as information on which the encoded enzymes are potentially involved in yolk absorption during embryogenesis, is unknown. Here, phylogenetic and gene structure analysis revealed that zebrafish ctsla and ctslb genes are similar, showing a highly conserved structure in comparison with human ctsl, while ctslc presents different exon organization together with an earlier evolution. Thus, ctslc appears to be evolved from a common ancestral ctsl-like gene, possibly through an early duplication event, whereas ctsla and ctslb may be originated from a second duplication mechanism. Zebrafish ctsla, ctslb and ctslc also showed different patterns of mRNA expression during embryogenesis and in adult tissues. While Ctsla transcripts were accumulated in embryos throughout development and in the adult ovary, those encoding Ctslb were detected only in embryos around the time of hatching as previously reported, and those for Ctslc appeared only in larvae and in some adult tissues, but not in the ovary. In zebrafish and killifish (Fundulus heteroclitus) embryos, Ctsla mRNA was first detected in blastomers, and later in development it was localized in cells of the yolk syncytial layer, an embryonic structure involved in yolk absorption. These data therefore suggested that Ctsla is most likely the putative protease involved in yolk processing in fish embryos, while Ctslc seems not to be required during early embryogenesis in zebrafish.

The role of aquaporin 3 in teleost fish

The aquaporin isoform, AQP3 has now been identified in a number of different teleost fish species, with additional DNA sequence information on AQP3 genes in further fish species available in genome databases. In zebrafish (Danio rerio), the AQP3 gene is present as two duplicate isoforms resulting from a teleostean fish genome-wide duplication. A further splicoform/isoform has also been identified in rainbow trout (Oncorhynchus mykiss). The identification of these AQP3 isoforms in other fish species is consequently explored. The role of AQP3 in physiological/osmoregulatory processes, in various teleost organs is then described. In teleost gill, AQP3 is expressed in 'chloride' cells, and in some species, in other epithelial cell types, where it may have a number of different functions including the prevention of dehydration. In eel esophagus, immunohistochemistry shows that AQP3 is expressed in surface epithelial cells in the anterior esophagus, but in mucus cells within the epithelium of the posterior esophagus. In eel intestine, AQP3 is found in macrophage-like cells and probably plays no part in osmoregulatory processes. In the rectum, as in the posterior esophagus AQP3 is expressed in mucus cells. In eel kidney, AQP3 is expressed in a subset of renal tubules, and localizes to the apical pole of tubule cells. There is no apparent change in the location or protein abundance of renal AQP3 following the acclimation of eels from freshwater to seawater.