Evolutionary change of the numbers of homeobox genes in bilateral animals

Wide-scale identification of novel/eliminated genes responsible for evolutionary transformations

BACKGROUND

It is generally accepted that most evolutionary transformations at the phenotype level are associated either with rearrangements of genomic regulatory elements, which control the activity of gene networks, or with changes in the amino acid contents of proteins. Recently, evidence has accumulated that significant evolutionary transformations could also be associated with the loss/emergence of whole genes. The targeted identification of such genes is a challenging problem for both bioinformatics and evo-devo research.

RESULTS

To solve this problem we propose the WINEGRET method, named after the first letters of the title. Its main idea is to search for genes that satisfy two requirements: first, the desired genes were lost/emerged at the same evolutionary stage at which the phenotypic trait of interest was lost/emerged, and second, the expression of these genes changes significantly during the development of the trait of interest in the model organism. To verify the first requirement, we do not use existing databases of orthologs, but rely purely on gene homology and local synteny by using some novel quickly computable conditions. Genes satisfying the second requirement are found by deep RNA sequencing. As a proof of principle, we used our method to find genes absent in extant amniotes (reptiles, birds, mammals) but present in anamniotes (fish and amphibians), in which these genes are involved in the regeneration of large body appendages. As a result, 57 genes were identified. For three of them, c-c motif chemokine 4, eotaxin-like, and a previously unknown gene called here sod4, essential roles for tail regeneration were demonstrated. Noteworthy, we established that the latter gene belongs to a novel family of Cu/Zn-superoxide dismutases lost by amniotes, SOD4.

CONCLUSIONS

We present a method for targeted identification of genes whose loss/emergence in evolution could be associated with the loss/emergence of a phenotypic trait of interest. In a proof-of-principle study, we identified genes absent in amniotes that participate in body appendage regeneration in anamniotes. Our method provides a wide range of opportunities for studying the relationship between the loss/emergence of phenotypic traits and the loss/emergence of specific genes in evolution.

Deciphering Evolutionary Dynamics of WRKY I Genes in Rosaceae Species

WRKY transcription factors participate in various regulation processes at different developmental stages in higher plants. Here, 98 WRKY I genes were identified in seven Rosaceae species. The WRKY I genes are highly enriched in some subgroups and are selectively expanded in Chinese pear [Pyrus bretschneideri (P. bretschneideri)] and apple [Malus domestica (M. domestica)]. By searching for intra-species gene microsynteny, we found the majority of chromosomal segments for WRKY I-containing segments in both P. bretschneideri and M. domestica genomes, while paired segments were hardly identified in the other five genomes. Furthermore, we analyzed the environmental selection pressure of duplicated WRKY I gene pairs, which indicated that the strong purifying selection for WRKY domains may contribute to the stability of its structure and function. The expression patterns of duplication PbWRKY genes revealed that functional redundancy for some of these genes was derived from common ancestry and neo-functionalization or sub-functionalization for some of them. This study traces the evolution of WRKY I genes in Rosaceae genomes and lays the foundation for functional studies of these genes in the future. Our results also show that the rates of gene loss and gain in different Rosaceae genomes are far from equilibrium.

Evaluating the Expression of Candidate Homeobox Genes and Their Role in Local-Site Inflammation in Mucosal Tissue Obtained from Children with Non-Syndromic Cleft Lip and Palate

Craniofacial development including palatogenesis is a complex process which requires an orchestrated and spatiotemporal expression of various genes and factors for proper embryogenesis and organogenesis. One such group of genes essential for craniofacial development is the homeobox genes, transcriptional factors that are commonly associated with congenital abnormalities. Amongst these genes, DLX4, HOXB3, and MSX2 have been recently shown to be involved in the etiology of non-syndromic cleft lip and palate. Hence, we investigated the gene and protein expression of these genes in normal and cleft affected mucosal tissue obtained from 22 children, along with analyzing their role in promoting local-site inflammation using NF-κB. Additionally, we investigated the role of PTX3, which plays a critical role in tissue remodeling and wound repair. We found a residual gene and protein expression of DLX4 in cleft mucosa, although no differences in gene expression levels of HOXB3 and MSX2 were noted. However, a significant increase in protein expression for these genes was noted in the cleft mucosa (p < 0.05), indicating increased cellular proliferation. This was coupled with a significant increase in NF-κB protein expression in cleft mucosa (p < 0.05), highlighting the role of these genes in promotion of pro-inflammatory environment. Finally, no differences in gene expression of PTX3 were noted.

Origins and Evolution of Cuticle Biosynthetic Machinery in Land Plants

The aerial epidermis of land plants is covered with a hydrophobic cuticle that protects the plant against environmental stresses. Although the mechanisms of cuticle biosynthesis have been extensively studied in model plants, particularly in seed plants, the origins and evolution of cuticle biosynthesis are not well understood. In this study, we performed a comparative genomic analysis of core components that mediate cuticle biosynthesis and characterized the chemical compositions and physiological parameters of cuticles from a broad set of embryophytes. Phylogenomic analysis revealed that the cuticle biosynthetic machinery originated in the last common ancestor of embryophytes. Coexpansion and coordinated expression are evident in core genes involved in the biosynthesis of two major cuticle components: the polymer cutin and cuticular waxes. Multispecies analyses of cuticle chemistry and physiology further revealed higher loads of both cutin and cuticular waxes in seed plants than in bryophytes as well as greater proportions of dihydroxy and trihydroxy acids, dicarboxylic acids, very-long-chain alkanes, and >C28 lipophilic compounds. This can be associated with land colonization and the formation of cuticles with enhanced hydrophobicity and moisture retention capacity. These findings provide insights into the evolution of plant cuticle biosynthetic mechanisms.

Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates

Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such "big data" sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new "principled" codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.

Genome-Wide DNA Methylation Analysis of Soybean Curled-Cotyledons Mutant and Functional Evaluation of a Homeodomain-Leucine Zipper (HD-Zip) I Gene GmHDZ20

DNA methylation is a major, conserved epigenetic modification that influences many biological processes. Cotyledons are specialized tissues that provide nutrition for seedlings at the early developmental stage. To investigate the patterns of genomic DNA methylation of germinated cotyledons in soybean (Glycine max) and its effect on cotyledon development, we performed a genome-wide comparative analysis of DNA methylation between the soybean curled-cotyledons (cco) mutant, which has abnormal cotyledons, and its corresponding wild type (WT) by whole-genome bisulfite sequencing. The cco mutant was methylated at more sites but at a slightly lower level overall than the WT on the whole-genome level. A total of 46 CG-, 92 CHG-, and 9723 CHH- (H = A, C, or T) differentially methylated genes (DMGs) were identified in cotyledons. Notably, hypomethylated CHH-DMGs were enriched in the gene ontology term "sequence-specific DNA binding transcription factor activity." We selected a DMG encoding a homeodomain-leucine zipper (HD-Zip) I subgroup transcription factor (GmHDZ20) for further functional characterization. GmHDZ20 localized to the nucleus and was highly expressed in leaf and cotyledon tissues. Constitutive expression of GmHDZ20 in Arabidopsis thaliana led to serrated rosette leaves, shorter siliques, and reduced seed number per silique. A yeast two-hybrid assay revealed that GmHDZ20 physically interacted with three proteins associated with multiple aspects of plant growth. Collectively, our results provide a comprehensive study of soybean DNA methylation in normal and aberrant cotyledons, which will be useful for the identification of specific DMGs that participate in cotyledon development, and also provide a foundation for future in-depth functional study of GmHDZ20 in soybean.

Extensive intraspecies cryptic variation in an ancient embryonic gene regulatory network

Innovations in metazoan development arise from evolutionary modification of gene regulatory networks (GRNs). We report widespread cryptic variation in the requirement for two key regulatory inputs, SKN-1/Nrf2 and MOM-2/Wnt, into the C. elegans endoderm GRN. While some natural isolates show a nearly absolute requirement for these two regulators, in others, most embryos differentiate endoderm in their absence. GWAS and analysis of recombinant inbred lines reveal multiple genetic regions underlying this broad phenotypic variation. We observe a reciprocal trend, in which genomic variants, or knockdown of endoderm regulatory genes, that result in a high SKN-1 requirement often show low MOM-2/Wnt requirement and vice-versa, suggesting that cryptic variation in the endoderm GRN may be tuned by opposing requirements for these two key regulatory inputs. These findings reveal that while the downstream components in the endoderm GRN are common across metazoan phylogeny, initiating regulatory inputs are remarkably plastic even within a single species.

Origin and Evolution of Core Components Responsible for Monitoring Light Environment Changes during Plant Terrestrialization

Light serves as the source of energy as well as an information signal for photosynthetic plants. During evolution, plants have acquired the ability to monitor environmental light radiation and adjust their developmental patterns to optimally utilize light energy for photosynthesis. The mechanisms of light perception and signal transduction have been comprehensively studied in past decades, mostly in a few model plants, including Arabidopsis thaliana. However, systematic analyses of the origin and evolution of core components involved in light perception and signaling are still lacking. In this study, we took advantage of the recently sequenced genomes and transcriptomes covering all the main Archaeplastida clades in the public domain to identify orthologous genes of core components involved in light perception and signaling and to reconstruct their evolutionary history. Our analyses suggested that acclimation to different distribution of light quality in new environments led to the origination of specific light signaling pathways in plants. The UVR8 (UV Resistance Locus 8) signaling pathway originated during the movement of plants from the deeper sea to shallow water and enabled plants to deal with ultraviolet B light (UV-B). After acquisition of UV-B adaptation, origination of the phytochrome signaling pathway helped plants to colonize water surface where red light became the prominent light energy source. The seedling emergence pathway, which is mediated by a combination of light and phytohormone signals that orchestrate plant growth pattern transitions, originated before the emergence of seed plants. Although cryptochromes and some key components of E3 ubiquitin ligase systems already existed before the divergence of the plant and animal kingdoms, the coevolution and optimization of light perception and downstream signal transduction components, including key transcription factors and E3 ubiquitin ligase systems, are evident during plant terrestrialization.

Homeobox Genes and Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the sixth most common type of cancer, and is the third leading cause of cancer-related deaths each year. It involves a multi-step progression and is strongly associated with chronic inflammation induced by the intake of environmental toxins and/or viral infections (i.e., hepatitis B and C viruses). Although several genetic dysregulations are considered to be involved in disease progression, the detailed regulatory mechanisms are not well defined. Homeobox genes that encode transcription factors with homeodomains control cell growth, differentiation, and morphogenesis in embryonic development. Recently, more aberrant expressions of Homeobox genes were found in a wide variety of human cancer, including HCC. In this review, we summarize the currently available evidence related to the role of Homeobox genes in the development of HCC. The objective is to determine the roles of this conserved transcription factor family and its potential use as a therapeutic target in future investigations.

Bioinformatical analysis and prediction of Nicotiana benthamiana bHLH transcription factors in Phytophthora parasitica resistance

The basic helix-loop-helix (bHLH) family, one of the largest transcription factor groups in plants, regulates many critical developmental processes. However, their functions in plant defense have not been extensively studied in Nicotiana benthamiana, an important model plant species for phytopathology. Here, we identified N. benthamiana bHLH genes (NbbHLHs) using a whole-genome searching approach, and found that the NbbHLHs are highly enriched and some subfamilies are selectively expanded in N. benthamiana. The results showed that gene duplication may be responsible for bHLH family expansion in this plant. Furthermore, we analyzed their expression profiles upon infection with Phytophthora parasitica. Finally, 28 candidate NbbHLHs may play important roles in Phytophthora pathogen resistance using cis-element analysis and protein-interaction network prediction. Taken together, our results established a platform for future studies of the gene family and provide molecular insights into plant immune responses against P. parasitica.

Knockdown of aristaless-like homeobox1 inhibits epithelial-mesenchymal transition through Wnt/β-catenin signaling pathway in melanoma cells

Aristaless-like homeobox1 (ALX1), a member of the ALX family, is capable of mediating survival and development of mesenchyme-derived elements in vertebrates and its mutation will prevent the fusion of frontonasal and maxillary elements. Recently, ALX1 has been reported to be associated with cancer progression. However, the specific roles of ALX1 in melanoma remain unclear. In this study, we investigated the expression pattern and biological functions of ALX1 in melanoma. We found that ALX1 was highly expressed in melanoma tissues and cell lines. Knockdown of ALX1 suppressed the proliferation and invasion of melanoma cells. Furthermore, we showed that ALX1 knockdown reversed the epithelial-mesenchymal transition (EMT) process in melanoma cells, which might be attributed to inactivation of the Wnt/β-catenin pathway. Taken together, this study provided a new insight into the role of ALX1 as a therapeutic target for melanoma treatment.

Evolutionary Analysis of the F-Box Gene Family in Saccharomycetaceae

F-box proteins are a core component of Skp1-Cul1-F-box (SCF) ubiquitin/ligase complexes and are involved in a lot of cellular processes in yeasts. However, the current knowledge of the molecular evolution of the F-box gene family in yeasts remains unclear. In this study, 136 F-box genes were identified in 10 yeast species of the Saccharomycetaceae. In addition to the F-box domain, the other six domains were identified in these F-box proteins. The evolutionary history of F-box gene numbers in 10 Saccharomycetaceae yeasts was reconstructed. Whole-genome duplication, interspersed repeats, and gene loss events were inferred. These events contributed to F-box gene number variation in the 10 yeast species. Eighty-seven and 33 positively selected sites were detected in program Selecton and Datamonkey web-server, respectively. Three of them were considered the significant positively selected sites, and 23 of them had changed radically in amino acid properties by using TreeSAAP. We investigated F-box gene number variation and underlying mechanisms, and selection patterns, all of which were beneficial to deeply understand genome evolution and figure out the function of the F-box proteins.

Zinc Finger-Homeodomain Transcriptional Factors (ZHDs) in Upland Cotton (Gossypium hirsutum): Genome-Wide Identification and Expression Analysis in Fiber Development

Zinc finger-homeodomain (ZHD) genes encode a family of plant-specific transcription factors that not only participate in the regulation of plant growth and development but also play an important role in the response to abiotic stress. The ZHD gene family has been studied in several model plants, including Solanum lycopersicum, Zea mays, Oryza sativa, and Arabidopsis thaliana. However, a comprehensive study of the genes of the ZHD family and their roles in fiber development and pigmentation in upland cotton has not been completed. To address this gap, we selected a brown fiber cultivar for our study; brown color in cotton is one of the most desired colors in the textile industry. The natural colored fibers require less processing and little dying, thereby eliminating dye costs and chemical residues. Using bioinformatics approaches, we identified 37 GhZHD genes from Gossypium hirsutum and then divided these genes into seven groups based on their phylogeny. The GhZHD genes were mostly conserved in each subfamily with minor variations in motif distribution and gene structure. These genes were largely distributed on 19 of the 26 upland cotton chromosomes. Among the Gossypium genomes, the paralogs and orthologs of the GhZHD genes were identified and further characterized. Furthermore, among the paralogs, we observed that the ZHD family duplications in Gossypium genomes (G. hirsutum, G. arboreum, and G. raimondii) were probably derived from segmental duplication or genome-wide duplication (GWD) events. Through a combination of qRT-PCR and proanthocyanidins (PA) accumulation analyses in brown cotton fibers, we concluded that the candidate genes involved in early fiber development and fiber pigment synthesis include the following: GhZHD29, GhZHD35, GhZHD30, GhZHD31, GhZHD11, GhZHD27, GhZHD18, GhZHD15, GhZHD16, GhZHD22, GhZHD6, GhZHD33, GhZHD13, GhZHD5, and GhZHD23. This study delivers insights into the evolution of the GhZHD genes in brown cotton, serves as a valuable resource for further studies, and identifies the conditions necessary for improving the quality of brown cotton fiber.

Genome-wide identification and expression analyses of the homeobox transcription factor family during ovule development in seedless and seeded grapes

Seedless grapes are of considerable importance for the raisin and table grape industries. Previous transcriptome analyses of seed development in grape revealed that genes encoding homeobox transcription factors were differentially regulated in seedless compared with seeded grape during seed development. In the present study, we identified a total of 73 homeobox-like genes in the grapevine genome and analyzed the genomic content and expression profiles of these genes. Based on domain architecture and phylogenetic analyses grape homeobox genes can be classified into eleven subfamilies. An analysis of the exon-intron structures and conserved motifs provided further insight into the evolutionary relationships between these genes. Evaluation of synteny indicated that segmental and tandem duplications have contributed greatly to the expansion of the grape homeobox gene superfamily. Synteny analysis between the grape and Arabidopsis genomes provided a potential functional relevance for these genes. The tissue-specific expression patterns of homeobox genes suggested roles in both vegetative and reproductive tissues. Expression profiling of these genes during the course of ovule development in seeded and seedless cultivars suggested a potential role in ovule abortion associated with seedlessness. This study will facilitate the functional analysis of these genes and provide new resources for molecular breeding of seedless grapes.

Nitric oxide regulates homeoprotein OTX1 and OTX2 expression in the rat myenteric plexus after intestinal ischemia-reperfusion injury

Neuronal and inducible nitric oxide synthase (nNOS and iNOS) play a protective and damaging role, respectively, on the intestinal neuromuscular function after ischemia-reperfusion (I/R) injury. To uncover the molecular pathways underlying this dichotomy we investigated their possible correlation with the orthodenticle homeobox proteins OTX1 and OTX2 in the rat small intestine myenteric plexus after in vivo I/R. Homeobox genes are fundamental for the regulation of the gut wall homeostasis both during development and in pathological conditions (inflammation, cancer). I/R injury was induced by temporary clamping the superior mesenteric artery under anesthesia, followed by 24 and 48 h of reperfusion. At 48 h after I/R intestinal transit decreased and was further reduced by N^ω-propyl-l-arginine hydrochloride (NPLA), a nNOS-selective inhibitor. By contrast this parameter was restored to control values by 1400W, an iNOS-selective inhibitor. In longitudinal muscle myenteric plexus (LMMP) preparations, iNOS, OTX1, and OTX2 mRNA and protein levels increased at 24 and 48 h after I/R. At both time periods, the number of iNOS- and OTX-immunopositive myenteric neurons increased. nNOS mRNA, protein levels, and neurons were unchanged. In LMMPs, OTX1 and OTX2 mRNA and protein upregulation was reduced by 1400W and NPLA, respectively. In myenteric ganglia, OTX1 and OTX2 staining was superimposed with that of iNOS and nNOS, respectively. Thus in myenteric ganglia iNOS- and nNOS-derived NO may promote OTX1 and OTX2 upregulation, respectively. We hypothesize that the neurodamaging and neuroprotective roles of iNOS and nNOS during I/R injury in the gut may involve corresponding activation of molecular pathways downstream of OTX1 and OTX2.NEW & NOTEWORTHY Intestinal ischemia-reperfusion (I/R) injury induces relevant alterations in myenteric neurons leading to dismotility. Nitrergic neurons seem to be selectively involved. In the present study the inference that both neuronal and inducible nitric oxide synthase (nNOS and iNOS) expressing myenteric neurons may undergo important changes sustaining derangements of motor function is reinforced. In addition, we provide data to suggest that NO produced by iNOS and nNOS regulates the expression of the vital transcription factors orthodenticle homeobox protein 1 and 2 during an I/R damage.

The Genomes of Two Bat Species with Long Constant Frequency Echolocation Calls

Bats can perceive the world by using a wide range of sensory systems, and some of the systems have become highly specialized, such as auditory sensory perception. Among bat species, the Old World leaf-nosed bats and horseshoe bats (rhinolophoid bats) possess the most sophisticated echolocation systems. Here, we reported the whole-genome sequencing and de novo assembles of two rhinolophoid bats – the great leaf-nosed bat (Hipposideros armiger) and the Chinese rufous horseshoe bat (Rhinolophus sinicus). Comparative genomic analyses revealed the adaptation of auditory sensory perception in the rhinolophoid bat lineages, probably resulting from the extreme selectivity used in the auditory processing by these bats. Pseudogenization of some vision-related genes in rhinolophoid bats was observed, suggesting that these genes have undergone relaxed natural selection. An extensive contraction of olfactory receptor gene repertoires was observed in the lineage leading to the common ancestor of bats. Further extensive gene contractions can be observed in the branch leading to the rhinolophoid bats. Such concordance suggested that molecular changes at one sensory gene might have direct consequences for genes controlling for other sensory modalities. To characterize the population genetic structure and patterns of evolution, we re-sequenced the genome of 20 great leaf-nosed bats from four different geographical locations of China. The result showed similar sequence diversity values and little differentiation among populations. Moreover, evidence of genetic adaptations to high altitudes in the great leaf-nosed bats was observed. Taken together, our work provided a useful resource for future research on the evolution of bats.

The Stepwise Increase in the Number of Transcription Factor Families in the Precambrian Predated the Diversification of Plants On Land

The colonization of the land by streptophytes and their subsequent radiation is a major event in Earth history. We report a stepwise increase in the number of transcription factor (TF) families and subfamilies in Archaeplastida before the colonization of the land. The subsequent increase in TF number on land was through duplication within existing TF families and subfamilies. Almost all subfamilies of the Homeodomain (HD) and basic Helix-Loop-Helix (bHLH) had evolved before the radiation of extant land plant lineages from a common ancestor. We demonstrate that the evolution of these TF families independently followed similar trends in both plants and metazoans; almost all extant HD and bHLH subfamilies were present in the first land plants and in the last common ancestor of bilaterians. These findings reveal that the majority of innovation in plant and metazoan TF families occurred in the Precambrian before the Phanerozoic radiation of land plants and metazoans.

Characterization of Rice Homeobox Genes, OsHOX22 and OsHOX24, and Over-expression of OsHOX24 in Transgenic Arabidopsis Suggest Their Role in Abiotic Stress Response

Homeobox transcription factors are well known regulators of plant growth and development. In this study, we carried out functional analysis of two candidate stress-responsive HD-ZIP I class homeobox genes from rice, OsHOX22, and OsHOX24. These genes were highly up-regulated under various abiotic stress conditions at different stages of rice development, including seedling, mature and reproductive stages. The transcript levels of these genes were enhanced significantly in the presence of plant hormones, including abscisic acid (ABA), auxin, salicylic acid, and gibberellic acid. The recombinant full-length and truncated homeobox proteins were found to be localized in the nucleus. Electrophoretic mobility shift assay established the binding of these homeobox proteins with specific DNA sequences, AH1 (CAAT(A/T)ATTG) and AH2 (CAAT(C/G)ATTG). Transactivation assays in yeast revealed the transcriptional activation potential of full-length OsHOX22 and OsHOX24 proteins. Homo- and hetero-dimerization capabilities of these proteins have also been demonstrated. Further, we identified putative novel interacting proteins of OsHOX22 and OsHOX24 via yeast-two hybrid analysis. Over-expression of OsHOX24 imparted higher sensitivity to stress hormone, ABA, and abiotic stresses in the transgenic Arabidopsis plants as revealed by various physiological and phenotypic assays. Microarray analysis revealed differential expression of several stress-responsive genes in transgenic lines as compared to wild-type. Many of these genes were found to be involved in transcriptional regulation and various metabolic pathways. Altogether, our results suggest the possible role of OsHOX22/OsHOX24 homeobox proteins as negative regulators in abiotic stress responses.

Genome-wide analysis of homeobox gene family in legumes: identification, gene duplication and expression profiling

Homeobox genes encode transcription factors that are known to play a major role in different aspects of plant growth and development. In the present study, we identified homeobox genes belonging to 14 different classes in five legume species, including chickpea, soybean, Medicago, Lotus and pigeonpea. The characteristic differences within homeodomain sequences among various classes of homeobox gene family were quite evident. Genome-wide expression analysis using publicly available datasets (RNA-seq and microarray) indicated that homeobox genes are differentially expressed in various tissues/developmental stages and under stress conditions in different legumes. We validated the differential expression of selected chickpea homeobox genes via quantitative reverse transcription polymerase chain reaction. Genome duplication analysis in soybean indicated that segmental duplication has significantly contributed in the expansion of homeobox gene family. The Ka/Ks ratio of duplicated homeobox genes in soybean showed that several members of this family have undergone purifying selection. Moreover, expression profiling indicated that duplicated genes might have been retained due to sub-functionalization. The genome-wide identification and comprehensive gene expression profiling of homeobox gene family members in legumes will provide opportunities for functional analysis to unravel their exact role in plant growth and development.

Massive expansion of Ubiquitination-related gene families within the Chlamydiae

Gene loss, gain, and transfer play an important role in shaping the genomes of all organisms; however, the interplay of these processes in isolated populations, such as in obligate intracellular bacteria, is less understood. Despite a general trend towards genome reduction in these microbes, our phylogenomic analysis of the phylum Chlamydiae revealed that within the family Parachlamydiaceae, gene family expansions have had pronounced effects on gene content. We discovered that the largest gene families within the phylum are the result of rapid gene birth-and-death evolution. These large gene families are comprised of members harboring eukaryotic-like ubiquitination-related domains, such as F-box and BTB-box domains, marking the largest reservoir of these proteins found among bacteria. A heterologous type III secretion system assay suggests that these proteins function as effectors manipulating the host cell. The large disparity in copy number of members in these families between closely related organisms suggests that nonadaptive processes might contribute to the evolution of these gene families. Gene birth-and-death evolution in concert with genomic drift might represent a previously undescribed mechanism by which isolated bacterial populations diversify.

Rapid birth-and-death evolution of imprinted snoRNAs in the Prader-Willi syndrome locus: implications for neural development in Euarchontoglires

Imprinted small nucleolar RNAs (snoRNAs) are only found in eutherian genomes and closely related to brain functions. A complex human neurological disease, Prader-Willi syndrome (PWS), is primarily attributed to the deletion of imprinted snoRNAs in chromosome 15q11-q13. Here we investigated the snoRNA repertoires in the PWS locus of 12 mammalian genomes and their evolution processes. A total of 613 imprinted snoRNAs were identified in the PWS homologous loci and the gene number was highly variable across lineages, with a peak in Euarchontoglires. Lineage-specific gene gain and loss events account for most extant genes of the HBII-52 (SNORD115) and the HBII-85 (SNORD116) gene family, and remarkable high gene-birth rates were observed in the primates and the rodents. Meanwhile, rapid sequence substitution occurred only in imprinted snoRNA genes, rather than their flanking sequences or the protein-coding genes located in the same imprinted locus. Strong selective constraints on the functional elements of these imprinted snoRNAs further suggest that they are subjected to birth-and-death evolution. Our data suggest that the regulatory role of HBII-52 on 5-HT_2CR pre-mRNA might originate in the Euarchontoglires through adaptive process. We propose that the rapid evolution of PWS-related imprinted snoRNAs has contributed to the neural development of Euarchontoglires.

Comparative genomics provide insights into evolution of trichoderma nutrition style

Saprotrophy on plant biomass is a recently developed nutrition strategy for Trichoderma. However, the physiology and evolution of this new nutrition strategy is still elusive. We report the deep sequencing and analysis of the genome of Trichoderma longibrachiatum, an efficient cellulase producer. The 31.7-Mb genome, smallest among the sequenced Trichoderma species, encodes fewer nutrition-related genes than saprotrophic T. reesei (Tr), including glycoside hydrolases and nonribosomal peptide synthetase-polyketide synthase. Homology and phylogenetic analyses suggest that a large number of nutrition-related genes, including GH18 chitinases, β-1,3/1,6-glucanases, cellulolytic enzymes, and hemicellulolytic enzymes, were lost in the common ancestor of T. longibrachiatum (Tl) and Tr. dN/dS (ω) calculation indicates that all the nutrition-related genes analyzed are under purifying selection. Cellulolytic enzymes, the key enzymes for saprotrophy on plant biomass, are under stronger purifying selection pressure in Tl and Tr than in mycoparasitic species, suggesting that development of the nutrition strategy of saprotrophy on plant biomass has increased the selection pressure. In addition, aspartic proteases, serine proteases, and metalloproteases are subject to stronger purifying selection pressure in Tl and Tr, suggesting that these enzymes may also play important roles in the nutrition. This study provides insights into the physiology and evolution of the nutrition strategy of Trichoderma.

Gene duplication and fragment recombination drive functional diversification of a superfamily of cytoplasmic effectors in Phytophthora sojae

Phytophthora and other oomycetes secrete a large number of putative host cytoplasmic effectors with conserved FLAK motifs following signal peptides, termed crinkling and necrosis inducing proteins (CRN), or Crinkler. Here, we first investigated the evolutionary patterns and mechanisms of CRN effectors in Phytophthora sojae and compared them to two other Phytophthora species. The genes encoding CRN effectors could be divided into 45 orthologous gene groups (OGG), and most OGGs unequally distributed in the three species, in which each underwent large number of gene gains or losses, indicating that the CRN genes expanded after species evolution in Phytophthora and evolved through pathoadaptation. The 134 expanded genes in P. sojae encoded family proteins including 82 functional genes and expressed at higher levels while the other 68 genes encoding orphan proteins were less expressed and contained 50 pseudogenes. Furthermore, we demonstrated that most expanded genes underwent gene duplication or/and fragment recombination. Three different mechanisms that drove gene duplication or recombination were identified. Finally, the expanded CRN effectors exhibited varying pathogenic functions, including induction of programmed cell death (PCD) and suppression of PCD through PAMP-triggered immunity or/and effector-triggered immunity. Overall, these results suggest that gene duplication and fragment recombination may be two mechanisms that drive the expansion and neofunctionalization of the CRN family in P. sojae, which aids in understanding the roles of CRN effectors within each oomycete pathogen.

Genome-wide identification and expression profiling of ankyrin-repeat gene family in maize

Members of the ankyrin repeats (ANK) gene family encode ANK domain that are common in diverse organisms and play important roles in cell growth and development, such as cell-cell signal transduction and cell cycle regulation. Recently, genome-wide identification and evolutionary analyses of the ANK gene family have been carried out in Arabidopsis and rice. However, little is known regarding the ANK genes in the entire maize genome. In this study, we described the identification and structural characterization of 71 ANK genes in maize (ZmANK). Then, comprehensive bioinformatics analyses of ZmANK genes family were performed including phylogenetic, domain and motif analysis, chromosomal localization, intron/exon structural patterns, gene duplications and expression profiling. Domain composition analyses showed that ZmANK genes formed ten subfamilies. Five tandem duplications and 14 segmental duplications were identified in ZmANK genes. Furthermore, we took comparative analysis of the total ANK gene family in Arabidopsis, rice and maize, ZmANKs were more closely paired with OsANKs than with AtANKs. At last, expression profile analyses were performed. Forty-one members of ZmANK genes held EST sequences records. Semi-quantitative expression and microarray data analysis of these 41 ZmANK genes demonstrated that ZmANK genes exhibit a various expression pattern, suggesting that functional diversification of ZmANK genes family. The results will present significant insights to explore ANK genes expression and function in future studies in maize.

Phylogenetic and transcriptomic analysis of chemosensory receptors in a pair of divergent ant species reveals sex-specific signatures of odor coding

Ants are a highly successful family of insects that thrive in a variety of habitats across the world. Perhaps their best-known features are complex social organization and strict division of labor, separating reproduction from the day-to-day maintenance and care of the colony, as well as strict discrimination against foreign individuals. Since these social characteristics in ants are thought to be mediated by semiochemicals, a thorough analysis of these signals, and the receptors that detect them, is critical in revealing mechanisms that lead to stereotypic behaviors. To address these questions, we have defined and characterized the major chemoreceptor families in a pair of behaviorally and evolutionarily distinct ant species, Camponotus floridanus and Harpegnathos saltator. Through comprehensive re-annotation, we show that these ant species harbor some of the largest yet known repertoires of odorant receptors (Ors) among insects, as well as a more modest number of gustatory receptors (Grs) and variant ionotropic glutamate receptors (Irs). Our phylogenetic analyses further demonstrate remarkably rapid gains and losses of ant Ors, while Grs and Irs have also experienced birth-and-death evolution to different degrees. In addition, comparisons of antennal transcriptomes between sexes identify many chemoreceptors that are differentially expressed between males and females and between species. We have also revealed an agonist for a worker-enriched OR from C. floridanus, representing the first case of a heterologously characterized ant tuning Or. Collectively, our analysis reveals a large number of ant chemoreceptors exhibiting patterns of differential expression and evolution consistent with sex/species-specific functions. These differentially expressed genes are likely associated with sex-based differences, as well as the radically different social lifestyles observed between C. floridanus and H. saltator, and thus are targets for further functional characterization. Our findings represent an important advance toward understanding the molecular basis of social interactions and the differential chemical ecologies among ant species.

Diversity, phylogeny and expression patterns of Pou and Six homeodomain transcription factors in hydrozoan jellyfish Craspedacusta sowerbyi

Formation of all metazoan bodies is controlled by a group of selector genes including homeobox genes, highly conserved across the entire animal kingdom. The homeobox genes from Pou and Six classes are key members of the regulation cascades determining development of sensory organs, nervous system, gonads and muscles. Besides using common bilaterian models, more attention has recently been targeted at the identification and characterization of these genes within the basal metazoan phyla. Cnidaria as a diploblastic sister group to bilateria with simple and yet specialized organs are suitable models for studies on the sensory organ origin and the associated role of homeobox genes. In this work, Pou and Six homeobox genes, together with a broad range of other sensory-specific transcription factors, were identified in the transcriptome of hydrozoan jellyfish Craspedacusta sowerbyi. Phylogenetic analyses of Pou and Six proteins revealed cnidarian-specific sequence motifs and contributed to the classification of individual factors. The majority of the Craspedacusta sowerbyi Pou and Six homeobox genes are predominantly expressed in statocysts, manubrium and nerve ring, the tissues with sensory and nervous activities. The described diversity and expression patterns of Pou and Six factors in hydrozoan jellyfish highlight their evolutionarily conserved functions. This study extends the knowledge of the cnidarian genome complexity and shows that the transcriptome of hydrozoan jellyfish is generally rich in homeodomain transcription factors employed in the regulation of sensory and nervous functions.

Origins and evolution of microRNA genes in plant species

MicroRNAs (miRNAs) are among the most important regulatory elements of gene expression in animals and plants. However, their origin and evolutionary dynamics have not been studied systematically. In this paper, we identified putative miRNA genes in 11 plant species using the bioinformatic technique and examined their evolutionary changes. Our homology search indicated that no miRNA gene is currently shared between green algae and land plants. The number of miRNA genes has increased substantially in the land plant lineage, but after the divergence of eudicots and monocots, the number has changed in a lineage-specific manner. We found that miRNA genes have originated mainly by duplication of preexisting miRNA genes or protein-coding genes. Transposable elements also seem to have contributed to the generation of species-specific miRNA genes. The relative importance of these mechanisms in plants is quite different from that in Drosophila species, where the formation of hairpin structures in the genomes seems to be a major source of miRNA genes. This difference in the origin of miRNA genes between plants and Drosophila may be explained by the difference in the binding to target mRNAs between plants and animals. We also found that young miRNA genes are less conserved than old genes in plants as well as in Drosophila species. Yet, nearly half of the gene families in the ancestor of flowering plants have been lost in at least one species examined. This indicates that the repertoires of miRNA genes have changed more dynamically than previously thought during plant evolution.

Evolution of the Alx homeobox gene family: parallel retention and independent loss of the vertebrate Alx3 gene

The Alx gene family is implicated in craniofacial development and comprises two to four homeobox genes in each vertebrate genome analyzed. Using phylogenetics and comparative genomics, we show that the common ancestor of jawed vertebrates had three Alx genes descendent from the two-round genome duplications (Alx1, Alx3, Alx4), compared with a single amphioxus gene. Later in evolution one of the paralogues, Alx3, was lost independently from at least three different vertebrate lineages, whereas Alx1 and Alx4 were consistently retained. Comparison of spatial gene expression patterns reveals that the three mouse genes have equivalent craniofacial expression to the two chick and frog genes, suggesting that redundancy compensated for gene loss. We suggest that multiple independent loss of one Alx gene was predisposed by extensive and persistent overlap in gene expression between Alx paralogues. Even so, it is unclear whether it was coincidence or evolutionary bias that resulted in the same Alx gene being lost on each occasion, rather than different members of the gene family.

The homeodomain complement of the ctenophore Mnemiopsis leidyi suggests that Ctenophora and Porifera diverged prior to the ParaHoxozoa

Background

The much-debated phylogenetic relationships of the five early branching metazoan lineages (Bilateria, Cnidaria, Ctenophora, Placozoa and Porifera) are of fundamental importance in piecing together events that occurred early in animal evolution. Comparisons of gene content between organismal lineages have been identified as a potentially useful methodology for phylogenetic reconstruction. However, these comparisons require complete genomes that, until now, did not exist for the ctenophore lineage. The homeobox superfamily of genes is particularly suited for these kinds of gene content comparisons, since it is large, diverse, and features a highly conserved domain.

Results

We have used a next-generation sequencing approach to generate a high-quality rough draft of the genome of the ctenophore Mnemiopsis leidyi and subsequently identified a set of 76 homeobox-containing genes from this draft. We phylogenetically categorized this set into established gene families and classes and then compared this set to the homeodomain repertoire of species from the other four early branching metazoan lineages. We have identified several important classes and subclasses of homeodomains that appear to be absent from Mnemiopsis and from the poriferan Amphimedon queenslandica. We have also determined that, based on lineage-specific paralog retention and average branch lengths, it is unlikely that these missing classes and subclasses are due to extensive gene loss or unusually high rates of evolution in Mnemiopsis.

Conclusions

This paper provides a first glimpse of the first sequenced ctenophore genome. We have characterized the full complement of Mnemiopsis homeodomains from this species and have compared them to species from other early branching lineages. Our results suggest that Porifera and Ctenophora were the first two extant lineages to diverge from the rest of animals. Based on this analysis, we also propose a new name - ParaHoxozoa - for the remaining group that includes Placozoa, Cnidaria and Bilateria.

Expansion of hexose transporter genes was associated with the evolution of aerobic fermentation in yeasts

The genetic basis of organisms' adaptation to different environments is a central issue of molecular evolution. The budding yeast Saccharomyces cerevisiae and its relatives predominantly ferment glucose into ethanol even in the presence of oxygen. This was suggested to be an adaptation to glucose-rich habitats, but the underlying genetic basis of the evolution of aerobic fermentation remains unclear. In S. cerevisiae, the first step of glucose metabolism is transporting glucose across the plasma membrane, which is carried out by hexose transporter proteins. Although several studies have recognized that the rate of glucose uptake can affect how glucose is metabolized, the role of HXT genes in the evolution of aerobic fermentation has not been fully explored. In this study, we identified all members of the HXT gene family in 23 fully sequenced fungal genomes, reconstructed their evolutionary history to pinpoint gene gain and loss events, and evaluated their adaptive significance in the evolution of aerobic fermentation. We found that the HXT genes have been extensively amplified in the two fungal lineages that have independently evolved aerobic fermentation. In contrast, reduction of the number of HXT genes has occurred in aerobic respiratory species. Our study reveals a strong positive correlation between the copy number of HXT genes and the strength of aerobic fermentation, suggesting that HXT gene expansion has facilitated the evolution of aerobic fermentation.

Evolution of olfactory receptor genes in primates dominated by birth-and-death process

Olfactory receptor (OR) is a large family of G protein–coupled receptors that can detect odorant in order to generate the sense of smell. They constitute one of the largest multiple gene families in animals including primates. To better understand the variation in odor perception and evolution of OR genes among primates, we computationally identified OR gene repertoires in orangutans, marmosets, and mouse lemurs and investigated the birth-and-death process of OR genes in the primate lineage. The results showed that 1) all the primate species studied have no more than 400 intact OR genes, fewer than rodents and canine; 2) Despite the similar number of OR genes in the genome, the makeup of the OR gene repertoires between different primate species is quite different as they had undergone dramatic birth-and-death evolution with extensive gene losses in the lineages leading to current species; 3) Apes and Old World monkey (OWM) have similar fraction of pseudogenes, whereas New World monkey (NWM) have fewer pseudogenes. To measure the selective pressure that had affected the OR gene repertoires in primates, we compared the ratio of nonsynonymous with synonymous substitution rates by using 70 one-to-one orthologous quintets among five primate species. We found that OR genes showed relaxed selective constraints in apes (humans, chimpanzees, and orangutans) than in OWMs (macaques) and NWMs (marmosets). We concluded that OR gene repertoires in primates have evolved in such a way to adapt to their respective living environments. Differential selective constraints might play important role in the primate OR gene evolution in each primate species.

The phylogeny and expression pattern of APETALA2-like genes in rice

The multigene families undergo birth-and-death evolution and thus contribute to biological innovations. The APETALA2-like genes belong to the euAP2 group of the AP2 gene family. These genes are characterized by several distinct motifs and exist in ferns, gymnosperms, and angiosperms. The phylogenetic analysis indicated that these genes have undergone the birth-and-death evolution. The five APETALA2-like genes in rice (Oryza sativa L.) display redundant but distinct expression patterns as demonstrated by RT-PCR and in situ hybridization. The potential functions of these genes were discussed on the basis of phylogenetic and expression pattern.

The Homeodomain Resource: a comprehensive collection of sequence, structure, interaction, genomic and functional information on the homeodomain protein family

The Homeodomain Resource is a curated collection of sequence, structure, interaction, genomic and functional information on the homeodomain family. The current version builds upon previous versions by the addition of new, complete sets of homeodomain sequences from fully sequenced genomes, the expansion of existing curated homeodomain information and the improvement of data accessibility through better search tools and more complete data integration. This release contains 1534 full-length homeodomain-containing sequences, 93 experimentally derived homeodomain structures, 101 homeodomain protein-protein interactions, 107 homeodomain DNA-binding sites and 206 homeodomain proteins implicated in human genetic disorders.Database URL: The Homeodomain Resource is freely available and can be accessed at http://research.nhgri.nih.gov/homeodomain/

Dynamic evolution of bitter taste receptor genes in vertebrates

Background

Sensing bitter tastes is crucial for many animals because it can prevent them from ingesting harmful foods. This process is mainly mediated by the bitter taste receptors (T2R), which are largely expressed in the taste buds. Previous studies have identified some T2R gene repertoires, and marked variation in repertoire size has been noted among species. However, the mechanisms underlying the evolution of vertebrate T2R genes remain poorly understood.

Results

To better understand the evolutionary pattern of these genes, we identified 16 T2R gene repertoires based on the high coverage genome sequences of vertebrates and studied the evolutionary changes in the number of T2R genes during birth-and-death evolution using the reconciled-tree method. We found that the number of T2R genes and the fraction of pseudogenes vary extensively among species. Based on the results of phylogenetic analysis, we showed that T2R gene families in teleost fishes are more diverse than those in tetrapods. In addition to the independent gene expansions in teleost fishes, frogs and mammals, lineage-specific gene duplications were also detected in lizards. Furthermore, extensive gains and losses of T2R genes were detected in each lineage during their evolution, resulting in widely differing T2R gene repertoires.

Conclusion

These results further support the hypotheses that T2R gene repertoires are closely related to the dietary habits of different species and that birth-and-death evolution is associated with adaptations to dietary changes.

The evolution of animal chemosensory receptor gene repertoires: roles of chance and necessity

Chemosensory receptors are essential for the survival of organisms that range from bacteria to mammals. Recent studies have shown that the numbers of functional chemosensory receptor genes and pseudogenes vary enormously among the genomes of different animal species. Although much of the variation can be explained by the adaptation of organisms to different environments, it has become clear that a substantial portion is generated by genomic drift, a random process of gene duplication and deletion. Genomic drift also generates a substantial amount of copy-number variation in chemosensory receptor genes within species. It seems that mutation by gene duplication and inactivation has important roles in both the adaptive and non-adaptive evolution of chemosensation.

Reconciliation with non-binary species trees

Reconciliation extracts information from the topological incongruence between gene and species trees to infer duplications and losses in the history of a gene family. The inferred duplication-loss histories provide valuable information for a broad range of biological applications, including ortholog identification, estimating gene duplication times, and rooting and correcting gene trees. While reconciliation for binary trees is a tractable and well studied problem, there are no algorithms for reconciliation with non-binary species trees. Yet a striking proportion of species trees are non-binary. For example, 64% of branch points in the NCBI taxonomy have three or more children. When applied to non-binary species trees, current algorithms overestimate the number of duplications because they cannot distinguish between duplication and incomplete lineage sorting. We present the first algorithms for reconciling binary gene trees with non-binary species trees under a duplication-loss parsimony model. Our algorithms utilize an efficient mapping from gene to species trees to infer the minimum number of duplications in O(|V(G) | x (k(S) + h(S))) time, where |V(G)| is the number of nodes in the gene tree, h(S) is the height of the species tree and k(S) is the size of its largest polytomy. We present a dynamic programming algorithm which also minimizes the total number of losses. Although this algorithm is exponential in the size of the largest polytomy, it performs well in practice for polytomies with outdegree of 12 or less. We also present a heuristic which estimates the minimal number of losses in polynomial time. In empirical tests, this algorithm finds an optimal loss history 99% of the time. Our algorithms have been implemented in NOTUNG, a robust, production quality, tree-fitting program, which provides a graphical user interface for exploratory analysis and also supports automated, high-throughput analysis of large data sets.

Differential loss of embryonic globin genes during the radiation of placental mammals

The differential gain and loss of genes from homologous gene families represents an important source of functional variation among the genomes of different species. Differences in gene content between species are primarily attributable to lineage-specific gene gains via duplication and lineage-specific losses via deletion or inactivation. Here, we use a comparative genomic approach to investigate this process of gene turnover in the beta-globin gene family of placental mammals. By analyzing genomic sequence data from representatives of each of the main superordinal clades of placental mammals, we were able to reconstruct pathways of gene family evolution during the basal radiation of this physiologically and morphologically diverse vertebrate group. Our analysis revealed that an initial expansion of the nonadult portion of the beta-globin gene cluster in the ancestor of placental mammals was followed by the differential loss and retention of ancestral gene lineages, thereby generating variation in the complement of embryonic globin genes among contemporary species. The sorting of epsilon-, gamma-, and eta-globin gene lineages among the basal clades of placental mammals has produced species differences in the functional types of hemoglobin isoforms that can be synthesized during the course of embryonic development.

NK homeobox genes with choanocyte-specific expression in homoscleromorph sponges

Data on nonbilaterian animals (sponges, cnidarians, and ctenophores) have suggested that Antennapedia (ANTP) class homeobox genes played a crucial role in the early diversification of animal body plans. Estimates of ancestral gene diversity within this important class of developmental regulators have been mostly based on recent analyses of the complete genome of a demosponge species, leading to the proposal that all ANTP families found in nonsponges animals (eumetazoans) derived from an ancestral "proto-NK" six-gene cluster. However, a single sponge species cannot reveal ancestral metazoan traits, in particular because lineage-specific gene duplications or losses are likely to have occurred during the long history of the Porifera. We thus looked for ANTP genes by degenerate polymerase chain reaction search in five species belonging to the Homoscleromorpha, a sponge lineage recently phylogenetically classified outside demosponges and characterized by unique histological features. We identified new genes of the ANTP class called HomoNK. Our phylogenetic analyses placed HomoNK (without significant support) close to the NK6 and NK7 families of cnidarian and bilaterian ANTP genes and did not recover the monophyly of the proposed "proto-NK" cluster. Our expression analyses of the HomoNK gene OlobNK in adult Oscarella lobularis showed that this gene is a strict marker of choanocytes, the most typical sponge cell type characterized by an apical flagellum surrounded by a collar of microvilli. These results are discussed in the light of the predominant neurosensory expression of NK6 and NK7 genes in bilaterians and of the recent proposal that choanocytes could be the sponge homologs of sensory cells.

Early sexual origins of homeoprotein heterodimerization and evolution of the plant KNOX/BELL family

Developmental mechanisms that yield multicellular diversity are proving to be well conserved within lineages, generating interest in their origins in unicellular ancestors. We report that molecular regulation of the haploid-diploid transition in Chlamydomonas, a unicellular green soil alga, shares common ancestry with differentiation pathways in land plants. Two homeoproteins, Gsp1 and Gsm1, contributed by gametes of plus and minus mating types respectively, physically interact and translocate from the cytosol to the nucleus upon gametic fusion, initiating zygote development. Their ectopic expression activates zygote development in vegetative cells and, in a diploid background, the resulting zygotes undergo a normal meiosis. Gsm1/Gsp1 dyads share sequence homology with and are functionally related to KNOX/BELL dyads regulating stem-cell (meristem) specification in land plants. We propose that combinatorial homeoprotein-based transcriptional control, a core feature of the fungal/animal radiation, may have originated in a sexual context and enabled the evolution of land-plant body plans.

Behavioral and synaptic defects in C. elegans lacking the NK-2 homeobox gene ceh-28

C. elegans pharyngeal behavior consists of two distinct types of muscle contractions, termed pumping and peristalsis. Pumping ingests and concentrates bacteria in the anterior pharyngeal lumen, and it is occasionally followed by a transient peristaltic contraction that carries ingested bacteria through the posterior pharyngeal isthmus. These behaviors are controlled by a small pharyngeal nervous system consisting of 20 neurons that is almost completely independent of the extra-pharyngeal nervous system. The cholinergic motor neuron M4 controls peristalsis via synapses with the posterior isthmus muscles. Here we show that the NK-2 family homeobox gene ceh-28 is expressed in M4, where it regulates synapse assembly and peristalsis. ceh-28 mutants exhibit frequent and prolonged peristalses, and treatment with agonists or antagonists of muscarinic acetylcholine receptors can phenocopy or suppress ceh-28 mutant defects, respectively. Synapses in ceh-28 mutant M4 cells are irregularly spaced and sized, and they are abnormally located along the full length of the isthmus. We suggest that CEH-28 inhibits synaptogenesis, and that ceh-28 mutant behavioral defects result from excessive or ectopic stimulation of muscarinic acetylcholine receptors in the isthmus muscles.

Gene family evolution across 12 Drosophila genomes

Comparison of whole genomes has revealed large and frequent changes in the size of gene families. These changes occur because of high rates of both gene gain (via duplication) and loss (via deletion or pseudogenization), as well as the evolution of entirely new genes. Here we use the genomes of 12 fully sequenced Drosophila species to study the gain and loss of genes at unprecedented resolution. We find large numbers of both gains and losses, with over 40% of all gene families differing in size among the Drosophila. Approximately 17 genes are estimated to be duplicated and fixed in a genome every million years, a rate on par with that previously found in both yeast and mammals. We find many instances of extreme expansions or contractions in the size of gene families, including the expansion of several sex- and spermatogenesis-related families in D. melanogaster that also evolve under positive selection at the nucleotide level. Newly evolved gene families in our dataset are associated with a class of testes-expressed genes known to have evolved de novo in a number of cases. Gene family comparisons also allow us to identify a number of annotated D. melanogaster genes that are unlikely to encode functional proteins, as well as to identify dozens of previously unannotated D. melanogaster genes with conserved homologs in the other Drosophila. Taken together, our results demonstrate that the apparent stasis in total gene number among species has masked rapid turnover in individual gene gain and loss. It is likely that this genomic revolving door has played a large role in shaping the morphological, physiological, and metabolic differences among species.

Though comparative genome sequencing has revealed vast similarities in the total number of genes contained within closely related species, this similarity hides enormous complexities in the identity and number of constituent proteins. Species can differ in their complement of genes through both gene duplication and loss. Here we investigated the gain and loss of genes from the genomes of 12 fully sequenced Drosophila (fruit flies). We find high rates of gain and loss in all species and estimate that approximately one new gene is gained or lost every 60,000 years. We also find several hundred cases of extremely rapid gene turnover, with dozens of genes gained or lost in only a few million years. The highest turnover in gene number occurs in genes involved in sex and reproduction. Taken together, our results demonstrate that the apparent stasis in total gene number among species has masked rapid turnover in individual gene gain and loss. It is likely that this evolutionary revolving door has played a large role in shaping the morphological, physiological, and metabolic differences among species.

Bias in phylogenetic tree reconciliation methods: implications for vertebrate genome evolution

Background

Comparative genomic studies are revealing frequent gains and losses of whole genes via duplication and pseudogenization. One commonly used method for inferring the number and timing of gene gains and losses reconciles the gene tree for each gene family with the species tree of the taxa considered. Recent studies using this approach have found a large number of ancient duplications and recent losses among vertebrate genomes.

Results

I show that tree reconciliation methods are biased when the inferred gene tree is not correct. This bias places duplicates towards the root of the tree and losses towards the tips of the tree. I demonstrate that this bias is present when tree reconciliation is conducted on both multiple mammal and Drosophila genomes, and that lower bootstrap cut-off values on gene trees lead to more extreme bias. I also suggest a method for dealing with reconciliation bias, although this method only corrects for the number of gene gains on some branches of the species tree.

Conclusion

Based on the results presented, it is likely that most tree reconciliation analyses show biases, unless the gene trees used are exceptionally well-resolved and well-supported. These results cast doubt upon previous conclusions that vertebrate genome history has been marked by many ancient duplications and many recent gene losses.

Molecular evolution of glutathione S-transferases in the genus Drosophila

As classical phase II detoxification enzymes, glutathione S-transferases (GSTs) have been implicated in insecticide resistance and may have evolved in response to toxins in the niche-defining feeding substrates of Drosophila species. We have annotated the GST genes of the 12 Drosophila species with recently sequenced genomes and analyzed their molecular evolution. Gene copy number variation is attributable mainly to unequal crossing-over events in the large delta and epsilon clusters. Within these gene clusters there are also GST genes with slowly diverging orthologs. This implies that they have their own unique functions or have spatial/temporal expression patterns that impose significant selective constraints. Searches for positively selected sites within the GSTs identified G171K in GSTD1, a protein that has previously been shown to be capable of metabolizing the insecticide DDT. We find that the same radical substitution (G171K) in the substrate-binding domain has occurred at least three times in the Drosophila radiation. Homology-modeling places site 171 distant from the active site but adjacent to an alternative DDT-binding site. We propose that the parallel evolution observed at this site is an adaptive response to an environmental toxin and that sequencing of historical alleles suggests that this toxin was not a synthetic insecticide.

Genomic drift and copy number variation of sensory receptor genes in humans

The number of sensory receptor genes varies extensively among different mammalian species. This variation is believed to be caused partly by physiological requirements of animals and partly by genomic drift due to random duplication and deletion of genes. If the contribution of genomic drift is substantial, each species should contain a significant amount of copy number variation (CNV). We therefore investigated CNVs in sensory receptor genes among 270 healthy humans by using published CNV data. The results indicated that olfactory receptor (OR), taste receptor type 2, and vomeronasal receptor type 1 genes show a high level of intraspecific CNVs. In particular, >30% of the approximately 800 OR gene loci in humans were polymorphic with respect to copy number, and two randomly chosen individuals showed a copy number difference of approximately 11 in functional OR genes on average. There was no significant difference in the amount of CNVs between functional and nonfunctional OR genes. Because pseudogenes are expected to evolve in a neutral fashion, this observation suggests that functional OR genes also have evolved in a similar manner with respect to copy number change. In addition, we found that the evolutionary change of copy number of OR genes approximately follows the Gaussian process in probability theory, and the copy number divergence between populations has increased with evolutionary time. We therefore conclude that genomic drift plays an important role for generating intra- and interspecific CNVs of sensory receptor genes. Similar results were obtained when all annotated genes were analyzed.

Hox, Wnt, and the evolution of the primary body axis: insights from the early-divergent phyla

The subkingdom Bilateria encompasses the overwhelming majority of animals, including all but four early-branching phyla: Porifera, Ctenophora, Placozoa, and Cnidaria. On average, these early-branching phyla have fewer cell types, tissues, and organs, and are considered to be significantly less specialized along their primary body axis. As such, they present an attractive outgroup from which to investigate how evolutionary changes in the genetic toolkit may have contributed to the emergence of the complex animal body plans of the Bilateria. This review offers an up-to-date glimpse of genome-scale comparisons between bilaterians and these early-diverging taxa. Specifically, we examine these data in the context of how they may explain the evolutionary development of primary body axes and axial symmetry across the Metazoa. Next, we re-evaluate the validity and evolutionary genomic relevance of the zootype hypothesis, which defines an animal by a specific spatial pattern of gene expression. Finally, we extend the hypothesis that Wnt genes may be the earliest primary body axis patterning mechanism by suggesting that Hox genes were co-opted into this patterning network prior to the last common ancestor of cnidarians and bilaterians.

Polymorphisms in the homeobox gene OTX2 may be a risk factor for bipolar disorder

We investigated the possible involvement of OTX2, a homeobox gene crucial for forebrain development, in the pathogenesis of schizophrenia and bipolar disorder. The disruption of this gene results in cortical malformations and causes serotonergic and dopaminergic cells in the midbrain to be expressed in aberrant locations. Resequencing of DNA from OTX2 exons and surrounding introns from 60 individuals (15 schizophrenia, 15 bipolar disorder, 15 depression, and 15 control) revealed two intronic polymorphisms, rs2277499 (C/T) and rs28757218 (G/T), but no other variations. The minor allele of rs2277499 (T) did not associate with clinical diagnosis. However, using a Taqman genotyping assay, we found the rs28757218 minor allele (T) in 30 out of 720 (4.2%) individuals with bipolar disorder but only in 6 out of 526 (1.1%) control individuals (odds ratio 3.5, 95% confidence interval 1.4-10.4, P = 0.003). On the other hand, the rs28757218 minor allele was only found in 6 out of 458 (1.3%) individuals with schizophrenia. All individuals with the rs28757218 polymorphism were heterozygous for the allele. Based on this positive case-control association finding, we conclude that variations in OTX2 might confer risk for the development of bipolar disorder.

Classification and nomenclature of all human homeobox genes

Background

The homeobox genes are a large and diverse group of genes, many of which play important roles in the embryonic development of animals. Increasingly, homeobox genes are being compared between genomes in an attempt to understand the evolution of animal development. Despite their importance, the full diversity of human homeobox genes has not previously been described.

Results

We have identified all homeobox genes and pseudogenes in the euchromatic regions of the human genome, finding many unannotated, incorrectly annotated, unnamed, misnamed or misclassified genes and pseudogenes. We describe 300 human homeobox loci, which we divide into 235 probable functional genes and 65 probable pseudogenes. These totals include 3 genes with partial homeoboxes and 13 pseudogenes that lack homeoboxes but are clearly derived from homeobox genes. These figures exclude the repetitive DUX1 to DUX5 homeobox sequences of which we identified 35 probable pseudogenes, with many more expected in heterochromatic regions. Nomenclature is established for approximately 40 formerly unnamed loci, reflecting their evolutionary relationships to other loci in human and other species, and nomenclature revisions are proposed for around 30 other loci. We use a classification that recognizes 11 homeobox gene 'classes' subdivided into 102 homeobox gene 'families'.

Conclusion

We have conducted a comprehensive survey of homeobox genes and pseudogenes in the human genome, described many new loci, and revised the classification and nomenclature of homeobox genes. The classification scheme may be widely applicable to homeobox genes in other animal genomes and will facilitate comparative genomics of this important gene superclass.

Visualizing differences in phylogenetic information content of alignments and distinction of three classes of long-branch effects

BACKGROUND

Published molecular phylogenies are usually based on data whose quality has not been explored prior to tree inference. This leads to errors because trees obtained with conventional methods suppress conflicting evidence, and because support values may be high even if there is no distinct phylogenetic signal. Tools that allow an a priori examination of data quality are rarely applied.

RESULTS

Using data from published molecular analyses on the phylogeny of crustaceans it is shown that tree topologies and popular support values do not show existing differences in data quality. To visualize variations in signal distinctness, we use network analyses based on split decomposition and split support spectra. Both methods show the same differences in data quality and the same clade-supporting patterns. Both methods are useful to discover long-branch effects. We discern three classes of long branch effects. Class I effects consist of attraction of terminal taxa caused by symplesiomorphies, which results in a false monophyly of paraphyletic groups. Addition of carefully selected taxa can fix this effect. Class II effects are caused by drastic signal erosion. Long branches affected by this phenomenon usually slip down the tree to form false clades that in reality are polyphyletic. To recover the correct phylogeny, more conservative genes must be used. Class III effects consist of attraction due to accumulated chance similarities or convergent character states. This sort of noise can be reduced by selecting less variable portions of the data set, avoiding biases, and adding slower genes.

CONCLUSION

To increase confidence in molecular phylogenies an exploratory analysis of the signal to noise ratio can be conducted with split decomposition methods. If long-branch effects are detected, it is necessary to discern between three classes of effects to find the best approach for an improvement of the raw data.

Extensive gains and losses of olfactory receptor genes in mammalian evolution

Odor perception in mammals is mediated by a large multigene family of olfactory receptor (OR) genes. The number of OR genes varies extensively among different species of mammals, and most species have a substantial number of pseudogenes. To gain some insight into the evolutionary dynamics of mammalian OR genes, we identified the entire set of OR genes in platypuses, opossums, cows, dogs, rats, and macaques and studied the evolutionary change of the genes together with those of humans and mice. We found that platypuses and primates have <400 functional OR genes while the other species have 800–1,200 functional OR genes. We then estimated the numbers of gains and losses of OR genes for each branch of the phylogenetic tree of mammals. This analysis showed that (i) gene expansion occurred in the placental lineage each time after it diverged from monotremes and from marsupials and (ii) hundreds of gains and losses of OR genes have occurred in an order-specific manner, making the gene repertoires highly variable among different orders. It appears that the number of OR genes is determined primarily by the functional requirement for each species, but once the number reaches the required level, it fluctuates by random duplication and deletion of genes. This fluctuation seems to have been aided by the stochastic nature of OR gene expression.

The new mutation theory of phenotypic evolution

Recent studies of developmental biology have shown that the genes controlling phenotypic characters expressed in the early stage of development are highly conserved and that recent evolutionary changes have occurred primarily in the characters expressed in later stages of development. Even the genes controlling the latter characters are generally conserved, but there is a large component of neutral or nearly neutral genetic variation within and between closely related species. Phenotypic evolution occurs primarily by mutation of genes that interact with one another in the developmental process. The enormous amount of phenotypic diversity among different phyla or classes of organisms is a product of accumulation of novel mutations and their conservation that have facilitated adaptation to different environments. Novel mutations may be incorporated into the genome by natural selection (elimination of preexisting genotypes) or by random processes such as genetic and genomic drift. However, once the mutations are incorporated into the genome, they may generate developmental constraints that will affect the future direction of phenotypic evolution. It appears that the driving force of phenotypic evolution is mutation, and natural selection is of secondary importance.