Selection for Antimicrobial Peptide Diversity in Frogs Leads to Gene Duplication and Low Allelic Variation

Antimicrobial peptides are highly diverse pathogen-killing molecules. In many taxa, their evolution is characterized by positive selection and frequent gene duplication. It has been proposed that genes encoding antimicrobial peptides might be subject to balancing selection and/or an enhanced mutation rate, but these hypotheses have not been well evaluated because allelic variation has rarely been studied at antimicrobial peptide loci. We present an evolutionary analysis of novel antimicrobial peptide genes from leopard frogs, Rana. Our results demonstrate that a single genome contains multiple homologous copies, among which there is an excess of nonsynonymous nucleotide site divergence relative to that expected from synonymous site divergence. Thus, we confirm the trends of recurrent duplication and positive selection. Allelic variation is quite low relative to interspecies divergence, indicating a recent positive selective sweep with no evidence of balancing selection. Repeated gene duplication, rather than a balanced maintenance of divergent allelic variants at individual loci, appears to be how frogs have responded to selection for a diverse suite of antimicrobial peptides. Our data also support a pattern of enhanced synonymous site substitution in the mature peptide region of the gene, but we cannot conclude that this is due to an elevated mutation rate.     

Characterization of genes with tissue-specific differential expression patterns in Populus

Like many plants, Populus has an evolutionary history in which several, both recent and more ancient, genome duplication events have occurred and, therefore, constitutes an excellent model system for studying the functional evolution of genes. In the present study, we have focused on the properties of genes with tissue-specific differential expression patterns in poplar. We identified the genes by analyzing digital expression profiles derived by mapping 90,000+ expressed sequence tags (ESTs) from 18 sources to the predicted genes of Populus. Our sequence analysis suggests that tissue-specific differentially expressed genes have less diverged paralogs than average, indicating that gene duplication events is an important event in the pathway leading to this type of expression pattern. The functional analysis showed that genes coding for proteins involved in processes of functional importance for the specific tissue(s) in which they are expressed and genes coding for regulatory or responsive proteins are most common among the differentially expressed genes, demonstrating that the expression differentiation process is under strong selective pressure. Thus, our data supports a model where gene duplication followed by gene specialization or expansion of the regulatory and responsive networks leads to tissue-specific differential expression patterns. We have also searched for clustering of genes with similar expression pattern into gene-expression neighborhoods within the Populus genome. However, we could not detect any major clustering among the analyzed genes with highly specific expression patterns. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

The home stretch,a first analysis of the nearly completed genome of Rhodobacter sphaeroides 2.4.1

Rhodobacter sphaeroides 2.4.1 is an α-3 purple nonsulfur eubacterium with an extensive metabolic repertoire. Under anaerobic conditions, it is able to grow by photosynthesis, respiration and fermentation. Photosynthesis may be photoheterotrophic using organic compounds as both a carbon and a reducing source, or photoautotrophic using carbon dioxide as the sole carbon source and hydrogen as the source of reducing power. In addition, R. sphaeroides can grow both chemoheterotrophically and chemoautotrophically. The structural components of this metabolically diverse organism and their modes of integrated regulation are encoded by a genome of ∼4.5 Mb in size. The genome comprises two chromosomes CI and CII (2.9 and 0.9 Mb, respectively) and five other replicons. Sequencing of the genome has been carried out by two groups, the Joint Genome Institute, which carried out shotgun-sequencing of the entire genome and The University of Texas-Houston Medical School, which carried out a targeted sequencing strategy of CII. Here we describe our current understanding of the genome when data from both of these groups are combined. Previous work had suggested that the two chromosomes are equal partners sharing responsibilities for fundamental cellular processes. This view has been reinforced by our preliminary analysis of the virtually completed genome sequence. We also have some evidence to suggest that two of the plasmids, pRS241a and pRS241b encode chromosomal type functions and their role may be more than that of accessory elements, perhaps representing replicons in a transition state. This revised version was published online in June 2006 with corrections to the Cover Date.     

Establishing the genetic relationships between the wild and cultivated olives using a nuclear intron from nitrate reductase (nia-i3)

In Oleaceae the most outstanding biological issue is to clarify the taxonomic relationships of cultivated and wild olives. To establish the genetic relationships between the wild (Olea europaea subsp. europaea var. sylvestris (Mill.) Lehr.), the cultivated olive (Olea europaea subsp. europaea var. europaea), and other taxa of the genus Olea (Olea europaea subsp. cuspidata (Wall. ex G. Don) Cif., Olea europaea subsp. cerasiformis G. Kunkel & Sunding, Olea paniculata R. Br.) and other Oleaceae (represented by Ligustrum vulgaris) we carried out the amplification by polymerase chain reaction (PCR) and the sequencing of the third nuclear intron of the nitrate reductase gene (nia-i3). Sequence analyses showed the presence of two different functional variants of the intron (nia1 and nia2) in the Oleaceae, in addition to a shorter non-functional one. Notably, while the shortest and the nia1 variants were present in all the taxa analysed, the nia2 variant was present only in the wild and the cultivated olive. These data confirm the close phylogenetic relationship between wild and cultivated olives and suggest that this gene could be duplicated in these two taxa after its divergence from the remaining Oleaceae. The presence of a target for AflII enzyme in nia2 and its absence in nia1 variant enables easy distinction by PCR-RFLP between, on the one hand, wild and cultivated olive, and on the other the remaining subspecies of the Olea europaea L. complex (O. e. subsp. cuspidata and O. e. subsp. cerasiformis) as well as other Oleaceae (O. paniculata, L. vulgaris L.). Additionally, nia1 sequences provide useful information about phylogeny of the wild and cultivated olives inside the genus Olea.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Evolutionary Analysis for Functional Divergence of the Toll-Like Receptor Gene Family and Altered Functional Constraints

The Toll-like receptor (TLR) gene family consists of type 1 transmembrane receptors, which play essential roles in both innate immunity and adaptive immune response by ligand recognition and signal transduction. Using all available vertebrate TLR protein sequences, we inferred the phylogenetic tree and then characterized critical amino acid residues for functional divergence by detecting altered functional constraints after gene duplications. We found that the extracellular domain of TLR genes showed higher functional divergence than that of the cytoplasmic domain, particularly in the region between leucine-rich repeat (LRR) 10 and LRR 15 of TLR 4. Our finding supports the concept that sequence evolution in the extracellular domain may be responsible for the broad diversity of TLR ligand-binding affinity, providing a testable hypothesis for potential targets that could be verified by further experimentation.     

The Ghost of Selection Past: Rates of Evolution and Functional Divergence of Anciently Duplicated Genes

The duplication of genes and even complete genomes may be a prerequisite for major evolutionary transitions and the origin of evolutionary novelties. However, the evolutionary mechanisms of gene evolution and the origin of novel gene functions after gene duplication have been a subject of many debates. Recently, we compiled 26 groups of orthologous genes, which included one gene from human, mouse, and chicken, one or two genes from the tetraploid Xenopus and two genes from zebrafish. Comparative analysis and mapping data showed that these pairs of zebrafish genes were probably produced during a fish-specific genome duplication that occurred between 300 and 450 Mya, before the teleost radiation (Taylor et al. 2001). As discussed here, many of these retained duplicated genes code for DNA binding proteins. Different models have been developed to explain the retention of duplicated genes and in particular the subfunctionalization model of Force et al. (1999) could explain why so many developmental control genes have been retained. Other models are harder to reconcile with this particular set of duplicated genes. Most genes seem to have been subjected to strong purifying selection, keeping properties such as charge and polarity the same in both duplicates, although some evidence was found for positive Darwinian selection, in particular for Hox genes. However, since only the cumulative pattern of nucleotide substitutions can be studied, clear indications of positive Darwinian selection or neutrality may be hard to find for such anciently duplicated genes. Nevertheless, an increase in evolutionary rate in about half of the duplicated genes seems to suggest that either positive Darwinian selection has occurred or that functional constraints have been relaxed at one point in time during functional divergence. Received: 4 January 2001 / Accepted: 29 March 2001