Concerted and birth-and-death evolution of multigene families

Remarkable conservation of distinct nonclassical MHC class I lineages in divergent amphibian species

Nonclassical MHC class Ib (class Ib) genes are heterogeneous genes encoding molecules that are structurally similar to classical MHC class Ia molecules but with limited tissue distribution and polymorphism. Mammalian class Ib genes have diverse and often uncharacterized functions, and because of their rapid rate of evolution, class Ib phylogeny is difficult to establish. We have conducted an extensive genomic, molecular, and phylogenetic characterization of class Ib genes in two Xenopodinae amphibian species of different genera that diverged from a common ancestor as long ago as primates and rodents (∼65 million years). In contrast with the unsteadiness of mammalian class Ib genes, our results reveal an unusual degree of conservation of most Xenopodinae class Ib gene lineages, including a novel monogenic lineage represented by the divergent Xenopus laevis XNC10 gene and its unequivocal Silurana (Xenopus) tropicalis orthologue, SNC10. The preferential expression of this gene lineage by thymocytes themselves from the onset of thymic organogenesis is consistent with a specialized role of class Ib in early T cell development and suggests such a function is conserved in all tetrapods.

Genetic variation, Fcγ receptors, KIRs and infection: the evolution of autoimmunity

Recent work has emphasised the marked genetic variability that exists in the Fc receptor locus. This variation can contribute to the risk of autoimmune disease in both mice and humans, but can also have a profound impact on defence against infection. Using FcγRIIB and FcγRIIIB as examples, we demonstrate that variations associated with increased susceptibility to autoimmunity may be maintained in populations for their beneficial effect against infection. We examine the KIR locus from the same perspective and highlight similarities between the two loci. Intense selection pressure by pathogens presumably accounts for the marked variability within both regions and leads to susceptibility to autoimmunity for some alleles.

Nematogalectin, a nematocyst protein with GlyXY and galectin domains, demonstrates nematocyte-specific alternative splicing in Hydra

Taxonomically restricted genes or lineage-specific genes contribute to morphological diversification in metazoans and provide unique functions for particular taxa in adapting to specific environments. To understand how such genes arise and participate in morphological evolution, we have investigated a gene called nematogalectin in Hydra, which has a structural role in the formation of nematocysts, stinging organelles that are unique to the phylum Cnidaria. Nematogalectin is a 28-kDa protein with an N-terminal GlyXY domain (glycine followed by two hydrophobic amino acids), which can form a collagen triple helix, followed by a galactose-binding lectin domain. Alternative splicing of the nematogalectin transcript allows the gene to encode two proteins, nematogalectin A and nematogalectin B. We demonstrate that expression of nematogalectin A and B is mutually exclusive in different nematocyst types: Desmonemes express nematogalectin B, whereas stenoteles and isorhizas express nematogalectin B early in differentiation, followed by nematogalectin A. Like Hydra, the marine hydrozoan Clytia also has two nematogalectin transcripts, which are expressed in different nematocyte types. By comparison, anthozoans have only one nematogalectin gene. Gene phylogeny indicates that tandem duplication of nematogalectin B exons gave rise to nematogalectin A before the divergence of Anthozoa and Medusozoa and that nematogalectin A was subsequently lost in Anthozoa. The emergence of nematogalectin A may have played a role in the morphological diversification of nematocysts in the medusozoan lineage.

Complex phylogenetic distribution of a non-canonical genetic code in green algae

BACKGROUND

A non-canonical nuclear genetic code, in which TAG and TAA have been reassigned from stop codons to glutamine, has evolved independently in several eukaryotic lineages, including the ulvophycean green algal orders Dasycladales and Cladophorales. To study the phylogenetic distribution of the standard and non-canonical genetic codes, we generated sequence data of a representative set of ulvophycean green algae and used a robust green algal phylogeny to evaluate different evolutionary scenarios that may account for the origin of the non-canonical code.

RESULTS

This study demonstrates that the Dasycladales and Cladophorales share this alternative genetic code with the related order Trentepohliales and the genus Blastophysa, but not with the Bryopsidales, which is sister to the Dasycladales. This complex phylogenetic distribution whereby all but one representative of a single natural lineage possesses an identical deviant genetic code is unique.

CONCLUSIONS

We compare different evolutionary scenarios for the complex phylogenetic distribution of this non-canonical genetic code. A single transition to the non-canonical code followed by a reversal to the canonical code in the Bryopsidales is highly improbable due to the profound genetic changes that coincide with codon reassignment. Multiple independent gains of the non-canonical code, as hypothesized for ciliates, are also unlikely because the same deviant code has evolved in all lineages. Instead we favor a stepwise acquisition model, congruent with the ambiguous intermediate model, whereby the non-canonical code observed in these green algal orders has a single origin. We suggest that the final steps from an ambiguous intermediate situation to a non-canonical code have been completed in the Trentepohliales, Dasycladales, Cladophorales and Blastophysa but not in the Bryopsidales. We hypothesize that in the latter lineage an initial stage characterized by translational ambiguity was not followed by final reassignment of both stop codons to glutamine. Instead the standard code was retained by the disappearance of the ambiguously decoding tRNAs from the genome. We correlate the emergence of a non-canonical genetic code in the Ulvophyceae to their multinucleate nature.

The neutral theory of molecular evolution in the genomic era

The neutral theory of molecular evolution has been widely accepted and is the guiding principle for studying evolutionary genomics and the molecular basis of phenotypic evolution. Recent data on genomic evolution are generally consistent with the neutral theory. However, many recently published papers claim the detection of positive Darwinian selection via the use of new statistical methods. Examination of these methods has shown that their theoretical bases are not well established and often result in high rates of false-positive and false-negative results. When the deficiencies of these statistical methods are rectified, the results become largely consistent with the neutral theory. At present, genome-wide analyses of natural selection consist of collections of single-locus analyses. However, because phenotypic evolution is controlled by the interaction of many genes, the study of natural selection ought to take such interactions into account. Experimental studies of evolution will also be crucial.

Origin and diversification of the L-amino oxidase family in innate immune defenses of animals

L-amino acid oxidases (LAOs), because they produce hydrogen peroxide as a by-product, function in innate immune defenses of both vertebrates and mollusks. Phylogenetic analysis revealed two major subfamilies of LAOs: (1) a subfamily including LAOs from vertebrates and mainly from Terrabacteria and (2) a subfamily including LAOs from mollusks and Hydrobacteria. These subfamilies thus originated early in the history of life, implying that their innate immune functions in vertebrates and mollusks have evolved separately. Mammalian LAOs were found to belong to three separate clades: (1) LAO1, (2) LAO2, and (3) IL4I1. Phylogenetic analysis supported the hypothesis that LAO1 and LAO2 arose by a gene duplication prior to the divergence of marsupials from placental mammals, while IL4I1 duplicated from the ancestor of the LAO1 and LAO2 prior to the divergence of tetrapods from bony fishes. Mammalian IL4I1 clustered with LAOs from bony fishes, and these molecules shared a number of unique sequence features, including both amino acid replacements and a unique two-codon deletion. It is certain such unique features may be functionally important, especially three unique amino acid replacements in close proximity to the putative active site.

Towards the simplification of MHC typing protocols: targeting classical MHC class II genes in a passerine, the pied flycatcher Ficedula hypoleuca

Background

Major Histocompatibility Complex (MHC) has drawn the attention of evolutionary biologists due to its importance in crucial biological processes, such as sexual selection and immune response in jawed vertebrates. However, the characterization of classical MHC genes subjected to the effects of natural selection still remains elusive in many vertebrate groups. Here, we have tested the suitability of flanking intron sequences to guide the selective exploration of classical MHC genes driving the co-evolutionary dynamics between pathogens and their passerine (Aves, Order Passeriformes) hosts.

Findings

Intronic sequences flanking the usually polymorphic exon 2 were isolated from different species using primers sitting on conserved coding regions of MHC class II genes (β chain). Taking the pied flycatcher Ficedula hypoleuca as an example, we demonstrate that careful primer design can evade non-classical MHC gene and pseudogene amplification. At least four polymorphic and expressed loci were co-replicated using a single pair of primers in five non-related individuals (N = 28 alleles). The cross-amplification and preliminary inspection of similar MHC fragments in eight unrelated songbird taxa suggests that similar approaches can also be applied to other species.

Conclusions

Intron sequences flanking the usually polymorphic exon 2 may assist the specific investigation of classical MHC class II B genes in species characterized by extensive gene duplication and pseudogenization. Importantly, the evasion of non-classical MHC genes with a more specific function and non-functional pseudogenes may accelerate data collection and diminish lab costs. Comprehensive knowledge of gene structure, polymorphism and expression profiles may be useful not only for the selective examination of evolutionarily relevant genes but also to restrict chimera formation by minimizing the number of co-amplifying loci.

Rice structural variation: a comparative analysis of structural variation between rice and three of its closest relatives in the genus Oryza

Rapid progress in comparative genomics among the grasses has revealed similar gene content and order despite exceptional differences in chromosome size and number. Large- and small-scale genomic variations are of particular interest, especially among cultivated and wild species, as they encode rapidly evolving features that may be important in adaptation to particular environments. We present a genome-wide study of intermediate-sized structural variation (SV) among rice (Oryza sativa) and three of its closest relatives in the genus Oryza (Oryza nivara, Oryza rufipogon and Oryza glaberrima). We computationally identified regional expansions, contractions and inversions in the Oryza species genomes relative to O. sativa by combining data from paired-end clone alignments to the O. sativa reference genome and physical maps. A subset of the computational predictions was validated using a new approach for BAC size determination. The result was a confirmed catalog of 674 expansions (25-38 Mb) and 611 (4-19 Mb) contractions, and 140 putative inversions (14-19 Mb) between the three Oryza species and O. sativa. In the expanded regions unique to O. sativa we found enrichment in transposable elements (TEs): long terminal repeats (LTRs) were randomly located across the chromosomes, and their insertion times corresponded to the date of the A genome radiation. Also, rice-expanded regions contained an over-representation of single-copy genes related to defense factors in the environment. This catalog of confirmed SV in reference to O. sativa provides an entry point for future research in genome evolution, speciation, domestication and novel gene discovery.

Mosaic gene conversion after a tandem duplication of mtDNA sequence in Diomedeidae (albatrosses)

Although the tandem duplication of mitochondrial (mt) sequences, especially those of the control region (CR), has been detected in metazoan species, few studies have focused on the features of the duplicated sequence itself, such as the gene conversion rate, distribution patterns of the variation, and relative rates of evolution between the copies. To investigate the features of duplicated mt sequences, we partially sequenced the mt genome of 16 Phoebastria albatrosses belonging to three species (P. albatrus, P. nigripes, and P. immutabilis). More than 2,300 base pairs of tandemly-duplicated sequence were shared by all three species. The observed gene arrangement was shared in the three Phoebastria albatrosses and suggests that the duplication event occurred in the common ancestor of the three species. Most of the copies in each individual were identical or nearly identical, and were maintained through frequent gene conversions. By contrast, portions of CR domains I and III had different phylogenetic signals, suggesting that gene conversion had not occurred in those sections after the speciation of the three species. Several lines of data, including the heterogeneity of the rate of molecular evolution, nucleotide differences, and putative secondary structures, suggests that the two sequences in CR domain I are maintained through selection; however, additional studies into the mechanisms of gene conversion and mtDNA synthesis are required to confirm this hypothesis.

Ancient protostome origin of chemosensory ionotropic glutamate receptors and the evolution of insect taste and olfaction

Ionotropic glutamate receptors (iGluRs) are a highly conserved family of ligand-gated ion channels present in animals, plants, and bacteria, which are best characterized for their roles in synaptic communication in vertebrate nervous systems. A variant subfamily of iGluRs, the Ionotropic Receptors (IRs), was recently identified as a new class of olfactory receptors in the fruit fly, Drosophila melanogaster, hinting at a broader function of this ion channel family in detection of environmental, as well as intercellular, chemical signals. Here, we investigate the origin and evolution of IRs by comprehensive evolutionary genomics and in situ expression analysis. In marked contrast to the insect-specific Odorant Receptor family, we show that IRs are expressed in olfactory organs across Protostomia--a major branch of the animal kingdom that encompasses arthropods, nematodes, and molluscs--indicating that they represent an ancestral protostome chemosensory receptor family. Two subfamilies of IRs are distinguished: conserved "antennal IRs," which likely define the first olfactory receptor family of insects, and species-specific "divergent IRs," which are expressed in peripheral and internal gustatory neurons, implicating this family in taste and food assessment. Comparative analysis of drosophilid IRs reveals the selective forces that have shaped the repertoires in flies with distinct chemosensory preferences. Examination of IR gene structure and genomic distribution suggests both non-allelic homologous recombination and retroposition contributed to the expansion of this multigene family. Together, these findings lay a foundation for functional analysis of these receptors in both neurobiological and evolutionary studies. Furthermore, this work identifies novel targets for manipulating chemosensory-driven behaviours of agricultural pests and disease vectors.

Evolution of the Cinnamyl/Sinapyl Alcohol Dehydrogenase (CAD/SAD) gene family: the emergence of real lignin is associated with the origin of Bona Fide CAD

Lignin plays a vital role in plant adaptation to terrestrial environments. The cinnamyl alcohol dehydrogenase (CAD) catalyzes the last step in monolignol biosynthesis and might have contributed to the lignin diversity in plants. To investigate the evolutionary history and functional differentiation of the CAD gene family, we made a comprehensive evolutionary analysis of this gene family from 52 species, including bacteria, early eukaryotes and green plants. The phylogenetic analysis, together with gene structure and function, indicates that all members of land plants, except two of moss, could be divided into three classes. Members of Class I (bona fide CAD), generally accepted as the primary genes involved in the monolignol biosynthesis, are all from vascular plants, and form a robustly supported monophyletic group with the lycophyte CADs at the basal position. This class is also conserved in the predicted three-dimensional structure and the residues constituting the substrate-binding pocket of the proteins. Given that Selaginella has real lignin, the above evidence strongly suggests that the earliest occurrence of the bona fide CAD in the lycophyte could be directly correlated with the origin of lignin. Class II comprises members more similar to the aspen sinapyl alcohol dehydrogenase gene, and includes three groups corresponding to lycophyte, gymnosperm, and angiosperm. Class III is conserved in land plants. The three classes differ in patterns of evolution and expression, implying that functional divergence has occurred among them. Our study also supports the hypothesis of convergent evolution of lignin biosynthesis between red algae and vascular plants.

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well-known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.

Characterisation of MHC class II DRB genes in the northern tree shrew (Tupaia belangeri)

Genes of the major histocompatibility complex (MHC) mainly code for proteins of the immune system of jawed vertebrates. In particular, MHC class I and II cell surface proteins are crucial for the self/non-self discrimination of the adaptive immune system and are the most polymorphic genes in vertebrates. Positive selection, gene duplications and pseudogenes shape the face of the MHC and reflect a highly dynamic evolution. Here, we present for the first time data of the highly polymorphic MHC class II DRB exon 2 of a representative of the mammalian order scandentia, the northern tree shrew Tupaia belangeri. We found up to eight different alleles per individual and determined haplotype constitution by intensively studying their inheritance. The alleles were assigned to four putative loci, all of which were polymorphic. Only the most polymorphic locus was subject to positive selection within the antigen binding sites and only alleles of this locus were transcribed.

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.

Evolutionary conservation of amino acid composition in paralogous insect vitellogenins

Comparison of paralogous vitellogenins in 10 insect species representing six orders showed a remarkable degree of conservation of amino acid composition in spite of sequence differences. For example, the correlation between the percentages of the 20 amino acids in two vitellogenins from the beetle Tribolium castaneum was 0.975, even though the two amino acid sequences differed from each other at 49.4% of sites. There was a positive correlation between the frequency of occurrence of reciprocal pairs of amino acids in more distantly related paralogs, and this correlation was generally strongest when both of the amino acids in the pair were nutritionally essential. These results imply that conservation of amino acid composition occurs through amino acid replacements that result in a balanced loss and gain of each amino acid residue. Thus insect vitellogenins seem to be subject to an unusual kind of purifying selection, where the amino acid content is conserved rather than the sequence per se, selection apparently arising from the nutritional needs of the developing embryo appears to be responsible for maintaining the balance of amino acids.

R2 dynamics in Triops cancriformis (Bosc, 1801) (Crustacea, Branchiopoda, Notostraca): turnover rate and 28S concerted evolution

The R2 retrotransposon is here characterized in bisexual populations of the European crustacean Triops cancriformis. The isolated element matches well with the general aspects of the R2 family and it is highly differentiated from that of the congeneric North American Triops longicaudatus. The analysis of 5' truncations indicates that R2 dynamics in T. cancriformis populations show a high turnover rate as observed in Drosophila simulans. For the first time in the literature, though, individuals harboring truncation variants, but lacking the complete element, are found. Present results suggest that transposition-mediated deletion mechanisms, possibly involving genomic turnover processes acting on rDNAs, can dramatically decrease the copy number or even delete R2 from the ribosomal locus. The presence of R2 does not seem to impact on the nucleotide variation of inserted 28S rDNA with respect to the uninserted genes. On the other hand, a low level of polymorphism characterizes rDNA units because new 28S variants continuously spread across the ribosomal array. Again, the interplay between transposition-mediated deletion and molecular drive may explain this pattern.

Evolution of a bitter taste receptor gene cluster in a New World sparrow

Bitter taste perception likely evolved as a protective mechanism against the ingestion of harmful compounds in food. The evolution of the taste receptor type 2 (TAS2R) gene family, which encodes the chemoreceptors that are directly responsible for the detection of bitter compounds, has therefore been of considerable interest. Though TAS2R repertoires have been characterized for a number of species, to date the complement of TAS2Rs from just one bird, the chicken, which had a notably small number of TAS2Rs, has been established. Here, we used targeted mapping and genomic sequencing in the white-throated sparrow (Zonotrichia albicollis) and sample sequencing in other closely related birds to reconstruct the history of a TAS2R gene cluster physically linked to the break points of an evolutionary chromosomal rearrangement. In the white-throated sparrow, this TAS2R cluster encodes up to 18 functional bitter taste receptors and likely underwent a large expansion that predates and/or coincides with the radiation of the Emberizinae subfamily into the New World. In addition to signatures of gene birth-and-death evolution within this cluster, estimates of Ka/Ks for the songbird TAS2Rs were similar to those previously observed in mammals, including humans. Finally, comparison of the complete genomic sequence of the cluster from two common haplotypes in the white-throated sparrow revealed a number of nonsynonymous variants and differences in functional gene content within this species. These results suggest that interspecies and intraspecies genetic variability does exist in avian TAS2Rs and that these differences could contribute to variation in bitter taste perception in birds.

Origins and evolution of microRNA genes in Drosophila species

MicroRNAs (miRs) regulate gene expression at the posttranscriptional level. To obtain some insights into the origins and evolutionary patterns of miR genes, we have identified miR genes in the genomes of 12 Drosophila species by bioinformatics approaches and examined their evolutionary changes. The results showed that the extant and ancestral Drosophila species had more than 100 miR genes and frequent gains and losses of miR genes have occurred during evolution. Although many miR genes appear to have originated from random hairpin structures in intronic or intergenic regions, duplication of miR genes has also contributed to the generation of new miR genes. Estimating the rate of nucleotide substitution of miR genes, we have found that newly arisen miR genes have a substitution rate similar to that of synonymous nucleotide sites in protein-coding genes and evolve almost neutrally. This suggests that most new miR genes have not acquired any important function and would become inactive. By contrast, old miR genes show a substitution rate much lower than the synonymous rate. Moreover, paired and unpaired nucleotide sites of miR genes tend to remain unchanged during evolution. Therefore, once miR genes acquired their functions, they appear to have evolved very slowly, maintaining essentially the same structures for a long time.

Genome-wide classification and evolutionary analysis of the bHLH family of transcription factors in Arabidopsis, poplar, rice, moss, and algae

Basic helix-loop-helix proteins (bHLHs) are found throughout the three eukaryotic kingdoms and constitute one of the largest families of transcription factors. A growing number of bHLH proteins have been functionally characterized in plants. However, some of these have not been previously classified. We present here an updated and comprehensive classification of the bHLHs encoded by the whole sequenced genomes of Arabidopsis (Arabidopsis thaliana), Populus trichocarpa, Oryza sativa, Physcomitrella patens, and five algae species. We define a plant bHLH consensus motif, which allowed the identification of novel highly diverged atypical bHLHs. Using yeast two-hybrid assays, we confirm that (1) a highly diverged bHLH has retained protein interaction activity and (2) the two most conserved positions in the consensus play an essential role in dimerization. Phylogenetic analysis permitted classification of the 638 bHLH genes identified into 32 subfamilies. Evolutionary and functional relationships within subfamilies are supported by intron patterns, predicted DNA-binding motifs, and the architecture of conserved protein motifs. Our analyses reveal the origin and evolutionary diversification of plant bHLHs through differential expansions, domain shuffling, and extensive sequence divergence. At the functional level, this would translate into different subfamilies evolving specific DNA-binding and protein interaction activities as well as differential transcriptional regulatory roles. Our results suggest a role for bHLH proteins in generating plant phenotypic diversity and provide a solid framework for further investigations into the role carried out in the transcriptional regulation of key growth and developmental processes.

Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

Count: evolutionary analysis of phylogenetic profiles with parsimony and likelihood

SUMMARY

Count is a software package for the analysis of numerical profiles on a phylogeny. It is primarily designed to deal with profiles derived from the phyletic distribution of homologous gene families, but is suited to study any other integer-valued evolutionary characters. Count performs ancestral reconstruction, and infers family- and lineage-specific characteristics along the evolutionary tree. It implements popular methods employed in gene content analysis such as Dollo and Wagner parsimony, propensity for gene loss, as well as probabilistic methods involving a phylogenetic birth-and-death model.

AVAILABILITY

Count is available as a stand-alone Java application, as well as an application bundle for MacOS X, at the web site http://www.iro.umontreal.ca/ approximately csuros/gene_content/count.html. It can also be launched using Java Webstart from the same site. The software is distributed under a BSD-style license. Source code is available upon request from the author.

Evolvability and Speed of Evolutionary Algorithms in Light of Recent Developments in Biology

Biological and artificial evolutionary systems exhibit varying degrees of evolvability and different rates of evolution. Such quantities can be affected by various factors. Here, we review some evolutionary mechanisms and discuss new developments in biology that can potentially improve evolvability or accelerate evolution in artificial systems. Biological notions are discussed to the degree they correspond to notions in Evolutionary Computation. We hope that the findings put forward here can be used to design computational models of evolution that produce significant gains in evolvability and evolutionary speed.

The evolution of a high copy gene array in Arabidopsis

Local gene duplication is a prominent mechanism of gene copy number expansion. Elucidating the mechanisms by which local duplicates arise is necessary in understanding the evolution of genomes and their host organisms. Chromosome one of Arabidopsis thaliana contains an 81-gene array subdivided into 27 triplet units (t-units), with each t-unit containing three pre-transfer RNA genes. We utilized phylogenetic tree reconstructions and comparative genomics to order the events leading to the array's formation, and propose a model using unequal crossing-over as the primary mechanism of array formation. The model is supported by additional phylogenetic information from intergenic spacer sequences separating each t-unit, comparative analysis to an orthologous array of 12 t-units in the sister taxa Arabidopsis lyrata, and additional modeling using a stochastic simulation of orthologous array divergence. Lastly, comparative phylogenetic analysis demonstrates that the two orthologous t-unit arrays undergo concerted evolution within each taxa and are likely fluctuating in copy number under neutral evolutionary drift. These findings hold larger implications for future research concerning gene and genome evolution.

Duplication, concerted evolution and purifying selection drive the evolution of mosquito vitellogenin genes

Background

Mosquito vitellogenin (Vtg) genes belong to a small multiple gene family that encodes the major yolk protein precursors required for egg production. Multiple Vtg genes have been cloned and characterized from several mosquito species, but their origin and molecular evolution are poorly understood.

Results

Here we used in silico and molecular cloning techniques to identify and characterize the evolution of the Vtg gene family from the genera Culex, Aedes/Ochlerotatus, and Anopheles. We identified the probable ancestral Vtg gene among different mosquito species by its conserved association with a novel gene approximately one kilobase upstream of the start codon. Phylogenetic analysis indicated that the Vtg gene family arose by duplication events, but that the pattern of duplication was different in each mosquito genera. Signatures of purifying selection were detected in Culex, Aedes and Anopheles. Gene conversion is a major driver of concerted evolution in Culex, while unequal crossover is likely the major driver of concerted evolution in Anopheles. In Aedes, smaller fragments have undergone gene conversion events.

Conclusions

The study shows concerted evolution and purifying selection shaped the evolution of mosquito Vtg genes following gene duplication. Additionally, similar evolutionary patterns were observed in the Vtg genes from other invertebrate and vertebrate organisms, suggesting that duplication, concerted evolution and purifying selection may be the major evolutionary forces driving Vtg gene evolution across highly divergent taxa.

Circadian input kinases and their homologs in cyanobacteria: evolutionary constraints versus architectural diversification

The circadian input kinase A (cikA) gene encodes a protein relaying environmental signal to the central circadian oscillator in cyanobacteria. The CikA protein has a variable architecture and usually consists of four tandemly arrayed domains: GAF, histidine kinase (HisKA), histidine kinase-like ATPase (HATPase_c), and a pseudo-receiver (REC). Among them, HisKA and HATPase_c are the least polymorphic, and REC is not present in heterocystic filamentous cyanobacteria. CikA contains several conserved motifs that are likely important for circadian function. There are at least three types of circadian systems, each of which possesses a different set of circadian genes. The originally described circadian system (kaiABC system) possesses both cikA and kaiA, while the others lack either only cikA (kaiABC (Delta)) or both (kaiBC). The results we obtained allowed us to approximate the time of the cikA origin to be about 2600-2200 MYA and the time of its loss in the species with the kaiABC (Delta) or kaiBC system between 1100 and 600 MYA. Circadian specialization of CikA, as opposed to its non-circadian homologs, is a result of several factors, including the unique conserved domain architecture and high evolutionary constraints of some domains and regions, which were previously identified as critical for the circadian function of the gene.

Unique structural characteristics and evolution of a cluster of venom phospholipase A2 isozyme genes of Protobothrops flavoviridis snake

Protobothrops flavoviridis (Crotalinae) venom gland phospholipase A(2) (PLA(2)) isozyme genes have evolved in an accelerated manner to acquire diverse physiological activities in their products. For elucidation of the multiplication mechanism of PLA(2) genes, a 25,026 bp genome segment harboring five PLA(2) isozyme genes was obtained from Amami-Oshima P. flavoviridis liver and sequenced. The gene PfPLA 2 encoded [Lys(49)]PLA(2) called BPII, the gene PfPLA 4 neurotoxic [Asp(49)]PLA(2) called PLA-N, the gene PfPLA 5 basic [Asp(49)]PLA(2) called PLA-B, and PfPLA 1(psi) and PfPLA 3(psi) were the inactivated genes. The 5' truncated reverse transcriptase (RT) elements, whose intact forms constitute long interspersed nuclear elements (LINEs), were found in close proximity to the 3' end of PLA(2) genes and named PLA(2) gene-coupled RT fragments (PcRTFs). The facts that PcRTFs have the stem-loop and repetitive sequence in the 3' untranslated region (UTR) which is characteristic of CR1 LINEs suggest that PcRTFs are the debris of P. flavoviridis ancestral CR1 LINEs, denoted as PfCR1s. Since the associated pairs of PLA(2) genes and PcRTFs are arranged in tandem in the 25,026 bp segment, it is thought that an ancestral PLA(2) gene-PfCR1 unit (PfPLA-PfCR1) which was produced by retrotransposition of PfCR1 by itself to the 3' end of PLA(2) gene duplicated several times to form a multimer of PfPLA-PfCR1, a cluster of PLA(2) genes, in the period after Crotalinae and Viperinae snakes branched off. Recombinational hot spot of a 37bp segment, named Scomb, was found in the region 548 bp upstream from the TATA box of PLA(2) genes. Thus, it could be assumed that multiplication of PfPLA-PfCR1 occurred by unequal crossing over of the segment, -Scomb-PfPLA-PfCR1-Scomb-. The PfCR1 moieties were afterward disrupted in the 5' portion to PcRTFs. The detection of two types of PcRTFs different in length which were produced by elimination of two definitive sequences in PfCR1 moiety possibly by gene conversion clearly supports such process but not multiplication of the PLA(2) gene-PcRTF unit.

Evolutionary dynamics of the 5S rDNA gene family in the mussel Mytilus: mixed effects of birth-and-death and concerted evolution

In higher eukaryotes, the gene family encoding the 5S ribosomal RNA (5S rRNA) has been used (together with histones) to showcase the archetypal example of a gene family subject to concerted evolution. However, recent studies have revealed conspicuous features challenging the predictions of this model, including heterogeneity of repeat units, the presence of functional 5S gene variants as well as the existence of 5S rDNA divergent pseudogenes lacking traces of homogenization. In the present work, we have broadened the scope in the evolutionary study of ribosomal gene families by studying the 5S rRNA family in mussels, a model organism which stands out among other animals due to the heterogeneity it displays regarding sequence and organization. To this end, 48 previously unknown 5S rDNA units (coding and spacer regions) were sequenced in five mussel species, leading to the characterization of two new types of units (referred to here as small-beta 5S rDNA and gamma-5S rDNA) coexisting in the genome with alpha and beta rDNA units. The intense genetic dynamics of this family is further supported by the first description of an association between gamma-5S rDNA units and tRNA genes. Molecular evolutionary and phylogenetic analyses revealed an extensive lack of homology among spacer sequences belonging to different rDNA types, suggesting the presence of independent evolutionary pathways leading to their differentiation. Overall, our results suggest that the long-term evolution of the 5S rRNA gene family in mussels is most likely mediated by a mixed mechanism involving the generation of genetic diversity through birth-and-death, followed by a process of local homogenization resulting from concerted evolution in order to maintain the genetic identities of the different 5S units, probably after their transposition to independent chromosomal locations.

Phylogenetic detection of numerous gene duplications shared by animals, fungi and plants

BACKGROUND

Gene duplication is considered a major driving force for evolution of genetic novelty, thereby facilitating functional divergence and organismal diversity, including the process of speciation. Animals, fungi and plants are major eukaryotic kingdoms and the divergences between them are some of the most significant evolutionary events. Although gene duplications in each lineage have been studied extensively in various contexts, the extent of gene duplication prior to the split of plants and animals/fungi is not clear.

RESULTS

Here, we have studied gene duplications in early eukaryotes by phylogenetic relative dating. We have reconstructed gene families (with one or more orthogroups) with members from both animals/fungi and plants by using two different clustering strategies. Extensive phylogenetic analyses of the gene families show that, among nearly 2,600 orthogroups identified, at least 300 of them still retain duplication that occurred before the divergence of the three kingdoms. We further found evidence that such duplications were also detected in some highly divergent protists, suggesting that these duplication events occurred in the ancestors of most major extant eukaryotic groups.

CONCLUSIONS

Our phylogenetic analyses show that numerous gene duplications happened at the early stage of eukaryotic evolution, probably before the separation of known major eukaryotic lineages. We discuss the implication of our results in the contexts of different models of eukaryotic phylogeny. One possible explanation for the large number of gene duplication events is one or more large-scale duplications, possibly whole genome or segmental duplication(s), which provides a genomic basis for the successful radiation of early eukaryotes.

Tandem duplication of the FLC locus and the origin of a new gene in Arabidopsis related species and their functional implications in allopolyploids

Flowering time is an important adaptive trait and varies among Arabidopsis thaliana and its related species, including allopolyploids that are formed between A. thaliana and Arabidopsis arenosa. FLOWERING LOCUS C (FLC) inhibits early flowering in A. thaliana. A previous study showed that late-flowering A. arenosa contained two or more FLC alleles that were differentially expressed in Arabidopsis allotetraploids, but the genomic organization and evolution of FLC locus were unknown. Comparative sequence and evolutionary analyses were performed in FLC-containing genomic regions in A. thaliana, A. arenosa and Arabidopsis lyrata, and expression of FLC loci and alleles was examined in Arabidopsis allopolyploids. The FLC locus was tandemly duplicated in A. lyrata and triplicated in A. arenosa, and the tandem duplication event occurred after divergence from A. thaliana. Although FLC duplicates were highly conserved, their upstream sequences rapidly diverged. The third FLC copy in A. arenosa acquired a new splicing site through a point mutation in the intron and generated the new exon followed by an early stop codon, resulting in a novel MADS box gene. Flowering time variation in Arabidopsis allopolyploids is probably related to the expression diversity and/or copy number of multiple FLC loci. Moreover, exonization of intronic sequence is a mechanism for the origin of new genes.

Species recognition and clinical relevance of the zygomycetous genus Lichtheimia (syn. Absidia pro parte, Mycocladus)

The zygomycete genus Lichtheimia (syn. Absidia pro parte, Mycocladus) consists of saprotrophic fungi inhabiting soil or dead plant material. Lichtheimia corymbifera (syn. Absidia corymbifera, Mycocladus corymbifer) and Lichtheimia ramosa (syn. Absidia ramosa, Mycocladus ramosus) may cause fulminant infections in patients with impaired immunity. The present study investigated the species boundaries in Lichtheimia using genealogical concordance phylogenetic species recognition (by comparison of the genealogies of the internal transcribed spacer [ITS] sequence, the D1/D2 region of the large subunit [LSU], and actin), biological species recognition by mating tests, as well as morphological and physiological characteristics. The three molecular markers used were selected by evaluating the polymorphisms and paralogies of several loci, including those for beta-tubulin, translation elongation factor 1alpha, the two largest subunits of the RNA polymerase II (RPB1 and RPB2), the mitochondrial cytochrome c oxidase subunit I (COI), and the mitochondrial small-subunit (mtSSU) rDNA, among four strains belonging to different putative species. Comparing the genealogies of the ITS, LSU, and actin genes, we recognized seven phylogenetic species. However, mating tests did not show intrinsic reproductive barriers for two pairs of the phylogenetic species. Therefore, we regard five species in Lichtheima to be confirmed: Lichtheimia corymbifera, L. ornata comb. nov., L. ramosa, L. hyalospora, and L. sphaerocystis sp. nov. Only the first three species seem to have clinical relevance. Lichtheimia blakesleeana is reduced to a synonym of Lichtheimia hyalospora. We provide a detailed description of Lichtheimia sphaerocystis sp. nov. and a key for the identification of all accepted species identified in the present study on the basis of their morphological traits and growth at different temperatures.

Containment of extended length polymorphisms in silk proteins

The spider silk gene family to the current date has been developed by gene duplication and homogenization events as well as conservation of crucial sequence parts. These evolutionary processes have created an amazing diversity of silk types each associated with specific properties and functions. In addition, they have led to allelic and gene variants within a species as exemplified by the major ampullate spidroin 1 gene of Nephila clavipes. Due to limited numbers of individuals screened to date little is known about the extent of these heterogeneities and how they are finally manifested in the proteins. Using expanded sample sizes, we show that sequence variations expressed as deletions or insertions of tri-nucleotides lead to different sized and structured repetitive units throughout a silk protein. Moreover, major ampullate spidroins 1 can quite dramatically differ in their overall lengths; however, extreme variants do not spread widely in a spider population. This suggests that a certain size range stabilized by purifying selection is important for spidroin 1 gene integrity and protein function. More than one locus for spidroin 1 genes possibly exist within one individual genome, which are homogenized in size, are differentially expressed and give a spider a certain degree of adaptation on silk's composition and properties. Such mechanisms are shared to a lesser extent by the second major ampullate spidroin gene.

Nonallelic gene conversion in the genus Drosophila

Nonallelic gene conversion has been proposed as a major force in homogenizing the sequences of paralogous genes. In this work, we investigate the extent and characteristics of gene conversion among gene families in nine species of the genus Drosophila. We carried out a genome-wide study of 2855 gene families (including 17,742 genes) and determined that conversion events involved 2628 genes. The proportion of converted genes ranged across species from 1 to 9% when paralogs of all ages were included. Although higher levels of gene conversion were found among young gene duplicates, at most 1-2% of the coding sequences of these duplicates were affected by conversion. Using a second approach relying on gene family size changes and gene-tree/species-tree reconciliation methods, we estimate that only 1-15% of gene trees are misled by gene conversion, depending on the lineage considered. Several features of paralogous genes correlate with gene conversion, such as intra-/interchromosomal location, level of nucleotide divergence, and GC content, although we found no definitive evidence for biased substitution patterns. After considering species-specific differences in the age and distance between paralogs, we found a highly significant difference in the amount of gene conversion among species. In particular, members of the melanogaster group showed the lowest proportion of converted genes. Our data therefore suggest underlying differences in the mechanistic basis of gene conversion among species.

Analysis of copy number variations among diverse cattle breeds

Genomic structural variation is an important and abundant source of genetic and phenotypic variation. Here, we describe the first systematic and genome-wide analysis of copy number variations (CNVs) in modern domesticated cattle using array comparative genomic hybridization (array CGH), quantitative PCR (qPCR), and fluorescent in situ hybridization (FISH). The array CGH panel included 90 animals from 11 Bos taurus, three Bos indicus, and three composite breeds for beef, dairy, or dual purpose. We identified over 200 candidate CNV regions (CNVRs) in total and 177 within known chromosomes, which harbor or are adjacent to gains or losses. These 177 high-confidence CNVRs cover 28.1 megabases or approximately 1.07% of the genome. Over 50% of the CNVRs (89/177) were found in multiple animals or breeds and analysis revealed breed-specific frequency differences and reflected aspects of the known ancestry of these cattle breeds. Selected CNVs were further validated by independent methods using qPCR and FISH. Approximately 67% of the CNVRs (119/177) completely or partially span cattle genes and 61% of the CNVRs (108/177) directly overlap with segmental duplications. The CNVRs span about 400 annotated cattle genes that are significantly enriched for specific biological functions, such as immunity, lactation, reproduction, and rumination. Multiple gene families, including ULBP, have gone through ruminant lineage-specific gene amplification. We detected and confirmed marked differences in their CNV frequencies across diverse breeds, indicating that some cattle CNVs are likely to arise independently in breeds and contribute to breed differences. Our results provide a valuable resource beyond microsatellites and single nucleotide polymorphisms to explore the full dimension of genetic variability for future cattle genomic research.

The rate of unequal crossing over in the dumpy gene from Drosophila melanogaster

The PIGSFEAST (PF) exon of the Drosophila dumpy gene is undergoing concerted evolution by the process of unequal crossing over. We have developed a long-range PCR-based assay to amplify the approximately 12 kb long exon which contains variable numbers of 303 or 306 nt long repeats in a tandem array. We applied this procedure to mutation accumulation lines of Drosophila melanogaster established by M. Wayne and L. Higgins. Nine new repeat length variants were found in these lines allowing us to measure the rate of unequal crossing over in the PF exon. The rate, which for several reasons is an underestimate, is 7.05 x 10(-4) exchanges per generation.

Genome annotation and comparative analyses of the odorant-binding proteins and chemosensory proteins in the pea aphid Acyrthosiphon pisum

Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) are two families of small water-soluble proteins, abundant in the aqueous fluid surrounding olfactory receptor neurons in insect antennae. OBPs are involved in the first step of olfactory signal transduction, carrying airborne semiochemicals to the odorant receptors and can be classified into three groups: Classic OBPs, Plus-C OBPs and Atypical OBPs. Here, we identified and annotated genes encoding putative OBPs and CSPs in the pea aphid Acyrthosiphon pisum using bioinformatics. This identified genes encoding 13 Classic and two Plus-C OBPs and 13 CSPs. Homologous OBP sequences were also identified in nine other aphid species, allowing us to compare OBPs across several aphid and non-aphid species. We show that, although OBP sequences are divergent within a species and between different orders, there is a high similarity between orthologs within a range of aphid species. Furthermore, the phylogenetic relationships between OBP orthologs reflect the divergence of aphid evolution lineages. Our results support the 'birth-and-death' model as the major mechanism explaining aphid OBP sequence evolution, with the main force acting on the evolution being purifying selection.

Gene family size conservation is a good indicator of evolutionary rates

The evolution of duplicate genes has been a topic of broad interest. Here, we propose that the conservation of gene family size is a good indicator of the rate of sequence evolution and some other biological properties. By comparing the human–chimpanzee–macaque orthologous gene families with and without family size conservation, we demonstrate that genes with family size conservation evolve more slowly than those without family size conservation. Our results further demonstrate that both family expansion and contraction events may accelerate gene evolution, resulting in elevated evolutionary rates in the genes without family size conservation. In addition, we show that the duplicate genes with family size conservation evolve significantly more slowly than those without family size conservation. Interestingly, the median evolutionary rate of singletons falls in between those of the above two types of duplicate gene families. Our results thus suggest that the controversy on whether duplicate genes evolve more slowly than singletons can be resolved when family size conservation is taken into consideration. Furthermore, we also observe that duplicate genes with family size conservation have the highest level of gene expression/expression breadth, the highest proportion of essential genes, and the lowest gene compactness, followed by singletons and then by duplicate genes without family size conservation. Such a trend accords well with our observations of evolutionary rates. Our results thus point to the importance of family size conservation in the evolution of duplicate genes.

Massive gene losses in Asian cultivated rice unveiled by comparative genome analysis

Background

Rice is one of the most important food crops in the world. With increasing world demand for food crops, there is an urgent need to develop new cultivars that have enhanced performance with regard to yield, disease resistance, and so on. Wild rice is expected to provide useful genetic resources that could improve the present cultivated species. However, the quantity and quality of these unexplored resources remain unclear. Recent accumulation of the genomic information of both cultivated and wild rice species allows for their comparison at the molecular level. Here, we compared the genome sequence of Oryza sativa ssp. japonica with sets of bacterial artificial chromosome end sequences (BESs) from two wild rice species, O. rufipogon and O. nivara, and an African rice species, O. glaberrima.

Results

We found that about four to five percent of the BESs of the two wild rice species and about seven percent of the African rice could not be mapped to the japonica genome, suggesting that a substantial number of genes have been lost in the japonica rice lineage; however, their close relatives still possess their counterpart genes. We estimated that during evolution, O. sativa has lost at least one thousand genes that are still preserved in the genomes of the other species. In addition, our BLASTX searches against the non-redundant protein sequence database showed that disease resistance-related proteins were significantly overrepresented in the close relative-specific genomic portions. In total, 235 unmapped BESs of the three relatives matched 83 non-redundant proteins that contained a disease resistance protein domain, most of which corresponded to an NBS-LRR domain.

Conclusion

We found that the O. sativa lineage appears to have recently experienced massive gene losses following divergence from its wild ancestor. Our results imply that the domestication process accelerated large-scale genomic deletions in the lineage of Asian cultivated rice and that the close relatives of cultivated rice have the potential to restore the lost traits.

Differential expansion of the merozoite surface protein (msp)-7 gene family in Plasmodium species under a birth-and-death model of evolution

MSP-7 is a surface protein expressed by the Plasmodium merozoite as part of a protein-complex involved in initial interaction between merozoite and erythrocyte. Contigs of seven Plasmodium species were analyzed in order to identify all msp-7 family genes. The search identified annotated and unannotated open reading frames (ORFs) and showed an uneven number of msp-7 genes among the different species of the Plasmodium lineage. A phylogenetic analysis established the presence of at least two ancestral genes and identified various lineage- and species-specific duplication events. An estimation of synonymous (d(S)) and non-synonymous substitutions (d(N)) showed higher d(S) values compared to d(N) values, suggesting the action of purifying selection on these genes, moreover no changes in omega (evolutive rates) were found in codon models test. These data together with the data obtained from the Gu's type-I functional divergence test and comparisons between evolutionary rates among orthologous and paralogous genes suggest functional redundancy. Finally, an analysis of recombination events suggests that several sequences are undergoing such process and that this mechanism could therefore be playing an important role in the emergence of new sequences. We conclude that evolution of the msp-7 family is in agreement with a birth-and-death model of evolution, as msp-7 genes have expanded until reaching an optimal gene copy number in each Plasmodium species in order to adapt to different niches.

Evolutionary dynamics of olfactory receptor genes in chordates: interaction between environments and genomic contents

Olfaction is essential for the survival of animals. Versatile odour molecules in the environment are received by olfactory receptors (ORs), which form the largest multigene family in vertebrates. Identification of the entire repertories of OR genes using bioinformatics methods from the whole-genome sequences of diverse organisms revealed that the numbers of OR genes vary enormously, ranging from ~1,200 in rats and ~400 in humans to ~150 in zebrafish and ~15 in pufferfish. Most species have a considerable fraction of pseudogenes. Extensive phylogenetic analyses have suggested that the numbers of gene gains and losses are extremely large in the OR gene family, which is a striking example of the birth-and-death evolution. It appears that OR gene repertoires change dynamically, depending on each organism's living environment. For example, higher primates equipped with a well-developed vision system have lost a large number of OR genes. Moreover, two groups of OR genes for detecting airborne odorants greatly expanded after the time of terrestrial adaption in the tetrapod lineage, whereas fishes retain diverse repertoires of genes that were present in aquatic ancestral species. The origin of vertebrate OR genes can be traced back to the common ancestor of all chordate species, but insects, nematodes and echinoderms utilise distinctive families of chemoreceptors, suggesting that chemoreceptor genes have evolved many times independently in animal evolution.