Gene Conversion Drives Within Genic Sequences: Concerted Evolution of Ribosomal RNA Genes in Bacteria and Archaea

Multiple copies of a given ribosomal RNA gene family undergo concerted evolution such that sequences of all gene copies are virtually identical within a species although they diverge normally between species. In eukaryotes, gene conversion and unequal crossing over are the proposed mechanisms for concerted evolution of tandemly repeated sequences, whereas dispersed genes are homogenized by gene conversion. However, the homogenization mechanisms for multiple-copy, normally dispersed, prokaryotic rRNA genes are not well understood. Here we compared the sequences of multiple paralogous rRNA genes within a genome in 12 prokaryotic organisms that have multiple copies of the rRNA genes. Within a genome, putative sequence conversion tracts were found throughout the entire length of each individual rRNA genes and their immediate flanks. Individual conversion events convert only a short sequence tract, and the conversion partners can be any paralogous genes within the genome. Interestingly, the genic sequences undergo much slower divergence than their flanking sequences. Moreover, genomic context and operon organization do not affect rRNA gene homogenization. Thus, gene conversion underlies concerted evolution of bacterial rRNA genes, which normally occurs within genic sequences, and homogenization of flanking regions may result from co-conversion with the genic sequence. Received: 31 March 2000 / Accepted: 15 June 2000     

Genomic expansion across the albumin gene family on human chromosome 4q is directional

The albumin gene family arose in a series of duplication events which gave rise to symmetry in its structure. The four genes are tandemly linked on human chromosome 4q in the order: 5'ALB-5'AFP-5'ALF-5'DBP-centromere, and their introns display a symmetrical and repetitive pattern that is shared by members of the gene family. These repetitive motifs provide an internal reference, allowing observations of evolutionary changes within a single line (human) of evolutionary descent. The four genes and three intergenic regions between them increase in size as they get closer to the centromere. An invasion by multiple repetitive DNA elements may account, in part, for this expansion.     

Family business: the multidrug-resistance related protein (MRP) ABC transporter genes in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Despite the completion of the sequencing of the entire genome of Arabidopsis thaliana (L.) Heynh., the exact determination of each single gene and its function remains an open question. This is especially true for multigene families. An approach that combines analysis of genomic structure, expression data and functional genomics to ascertain the role of the members of the multidrug-resistance-related protein ( MRP) gene family, a subfamily of the ATP-binding cassette (ABC) transporters from Arabidopsis is presented. We used cDNA sequencing and alignment-based re-annotation of genomic sequences to define the exact genic structure of all known AtMRP genes. Analysis of promoter regions suggested different induction conditions even for closely related genes. Expression analysis for the entire gene family confirmed these assumptions. Phylogenetic analysis and determination of segmental duplication in the regions of AtMRP genes revealed that the evolution of the extraordinarily high number of ABC transporter genes in plants cannot solely be explained by polyploidisation during the evolution of the Arabidopsis genome. Interestingly MRP genes from Oryza sativa L. (rice; OsMRP) show very similar genomic structures to those from Arabidopsis. Screening of large populations of T-DNA-mutagenised lines of A. thaliana resulted in the isolation of AtMRP insertion mutants. This work opens the way for the defined analysis of a multigene family of important membrane transporters whose broad variety of functions expands their traditional role as cellular detoxifiers.     

Gene Sequence Phylogenies of the Family <Emphasis Type="Italic">Microbacteriaceae</Emphasis>

The type strains of 32 species of 13 genera of the family Microbacteriaceae were analysed with respect to gene-coding phylogeny for DNA gyrase subunit B (gyrB), RNA-polymerase subunit B (rpoB), recombinase A (recA), and polyphosphate kinase (ppk). The resulting gene trees were compared with the 16S rRNA gene phylogeny of the same strains. The topology of neighbour-joining and maximum parsimony phylogenetic trees, based on nucleic-acid sequences and protein sequences of housekeeping genes, differed from one another, and no gene tree was identical to that of the 16S rRNA gene tree. Most genera analysed containing >1 strain formed phylogenetically coherent taxa. The three pathovars of Curtobacterium flaccumfaciens clustered together to the exclusion of the type strains of other Curtobacterium species in all DNA - and protein-based analyses. In no tree did the distribution of a major taxonomic marker, i.e., diaminobutyric acid versus lysine and/or ornithine in the peptidoglycan, or acyl type of peptidoglycan, correlate with the phylogenetic position of the organisms. The changing phylogenetic position of Agrococcus jenensis was unexpected: This strain defined individual lineages in the trees based on 16S rRNA and gyrB and showed identity with Microbacterium saperdae in the other three gene trees.     

Localization of a new prostate-specific antigen-related serine protease gene, KLK4, is evidence for an expanded human kallikrein gene family cluster on chromosome 19q13.3-13.4.

The human tissue kallikrein (KLK) family of serine proteases, which is important in post-translational processing events, currently consists of just three genes-tissue kallikrein (KLK1), KLK2, and prostate-specific antigen (PSA) (KLK3)-clustered at chromosome 19q13. 3-13.4. We identified an expressed sequence tag from an endometrial carcinoma cDNA library with 50% identity to the three known KLK genes. Primers designed to putative exon 2 and exon 3 regions from this novel kallikrein-related sequence were used to polymerase chain reaction-screen five cosmids spanning 130 kb around the KLK locus on chromosome 19. This new gene, which we have named KLK4, is 25 kb downstream of the KLK2 gene and follows a region that includes two other putative KLK-like gene fragments. KLK4 spans 5.2 kb, has an identical genomic structure-five exons and four introns-to the other KLK genes and is transcribed on the reverse strand, in the same direction as KLK1 but opposite to that of KLK2 and KLK3. It encodes a 254-amino acid prepro-serine protease that is most similar (78% identical) to pig enamel matrix serine protease but is also 37% identical to PSA. These data suggest that the human kallikrein gene family locus on chromosome 19 is larger than previously thought and also indicate a greater sequence divergence within this family compared with the highly conserved rodent kallikrein genes.     

Spatial mating networks in insect-pollinated plants

Gene flow in plant populations is largely determined by landscape heterogeneity. Both the shape of the pollination kernel and the spatial distribution of trees affect the distribution of pollen grains and the genotypes they harbour, but little is known about the relative contribution of each of these two factors. Using genetic markers we build a spatial network of pollination events between any two trees in a population of Prunus mahaleb , an insect-pollinated plant. Then, we apply tools from the science of complex networks to characterize the structure of such a mating network. Although the distribution of the number of pollen donors per tree is quite homogeneous, the identity of donors is distributed heterogeneously across the population. This results in a population structured in well-defined modules or compartments, formed by a group of mother trees and their shared pollen donors. Long-distance pollination events decrease the modular structure by favouring mating among all available mates. This increases gene flow across the entire population, reducing its genetic structure, and potentially decreasing the role of genetic drift.     

Identification of Horizontal Gene Transfer from Phylogenetic Gene Trees

We suggest a new procedure to search for the genes with horizontal transfer events in their evolutionary history. The search is based on analysis of topology difference between the phylogenetic trees of gene (protein) groups and the corresponding phylogenetic species trees. Numeric values are introduced to measure the discrepancy between the trees. This approach was applied to analyze 40 prokaryotic genomes classified into 132 classes of orthologs. This resulted in a list of the candidate genes for which the hypothesis of horizontal transfer in evolution looks true.     

Identification of the horizontal gene transfer based on phylogenetic data

We suggest a new procedure to search for the genes with horizontal transfer events in their evolutionary history. The search is based on analysis of topology difference between the phylogenetic trees of gene (protein) groups and the corresponding phylogenetic species trees. Numeric values are introduced to measure the discrepancy between the trees. This approach was applied to analyze 40 prokaryotic genomes classified into 132 classes of orthologs. This resulted in a list of the candidate genes for which the hypothesis of horizontal transfer in evolution looks true.     

The plant glutathione transferase gene family: genomic structure, functions, expression and evolution

Glutathione transferases (GSTs) are ubiquitous, multifunctional proteins encoded by large gene families. In different plant species this gene family is comprised of 25–60 members, that can be grouped into six classes on the basis of sequence identity, gene organization and active site residues in the protein. The Phi and Tau classes are the most represented and are plant specific, while Zeta and Theta GSTs are found also in animals. Despite pronounced sequence and functional diversification, GSTs have maintained a highly conserved three-dimensional structure through evolution. Most GSTs are cytosolic and active as dimers, performing diverse catalytic as well as non-catalytic roles in detoxification of xenobiotics, prevention of oxidative damage and endogenous metabolism. Among their catalytic activities are the conjugation of electrophilic substrates to glutathione, glutathione-dependent isomerizations and reductions of toxic organic hydroperoxides. Their main non-catalytic role is as hormone and flavonoid ligandins. GST genes are predominantly organized in clusters non-randomly distributed in the genome. Phylogenetic studies indicate that plant GSTs have mainly evolved after the divergence of plants, the two prevalent Phi and Tau classes being the result of recent, multiple duplication events.     

Recent duplication and positive selection of the <Emphasis Type="Italic">GAGE</Emphasis> gene family

We report that the GAGE gene family of human Cancer/testis antigen (CTA) genes is likely to be in an early stage of its evolution. Members of this gene family are tandemly arranged on the X chromosome only in human, chimpanzee and macaque genomes and share a very high similarity. Phylogenetic trees show that the GAGE gene family began to duplicate after the split of human and chimpanzee. The estimated ages of the duplication events range from 4 million years ago to the present. The Ka/Ks values between the duplicates are significantly greater than 1, indicating that the mutation rate is higher in coding regions than non-coding regions of the genes, which suggests that the GAGE gene family is under positive selection. These findings indicate that the GAGE gene family may be a newly formed gene family undergoing rapid functional evolution. Yang Liu and Qiyun Zhu contributed equally to this work.     

Regulation of cellular chromatin state

The identity and functionality of eukaryotic cells is defined not just by their genomic sequence which remains constant between cell types, but by their gene expression profiles governed by epigenetic mechanisms. Epigenetic controls maintain and change the chromatin state throughout development, as exemplified by the setting up of cellular memory for the regulation and maintenance of homeotic genes in proliferating progenitors during embryonic development. Higher order chromatin structure in reversibly arrested adult stem cells also involves epigenetic regulation and in this review we highlight common trends governing chromatin states, focusing on quiescence and differentiation during myogenesis. Together, these diverse developmental modules reveal the dynamic nature of chromatin regulation providing fresh insights into the role of epigenetic mechanisms in potentiating development and differentiation.     

Phylogenetic congruence and discordance among one morphological and three molecular data sets from Pontederiaceae

A morphological data set and three sources of data from the chloroplast genome (two genes and a restriction site survey) were used to reconstruct the phylogenetic history of the pickerelweed family Pontederiaceae. The chloroplast data converged towards a single tree, presumably the true chloroplast phylogeny of the family. Unrooted trees estimated from each of the three chloroplast data sets were identical or extremely similar in shape to each other and mostly robustly supported. There was no evidence of significant heterogeneity among the data sets, and the few topological differences seen among unrooted trees from each chloroplast data set are probably artifacts of sampling error on short branches. Despite well-documented differences in rates of evolution for different characters in individual data sets, equally weighted parsimony permits accurate reconstructions of chloroplast relationships in Pontederiaceae. A separate morphology-based data set yielded trees that were very different from the chloroplast trees. Although there was substantial support from the morphological evidence for several major clades supported by chloroplast trees, most of the conflicting phylogenetic structure on the morphology trees was not robust. Nonetheless, several statistical tests of incongruence indicate significant heterogeneity between molecules and morphology. The source of this apparent incongruence appears to be a low ratio of phylogenetic signal to noise in the morphological data.     

Mammalian alcohol dehydrogenase of higher classes: analyses of human ADH5 and rat ADH6

Alcohol dehydrogenases (ADH) of classes V and VI, ADH5 and ADH6, have been defined in man and rodents, respectively. Sequence data have been obtained at cDNA and genomic levels, but limited data are available for functionality and substrate repertoire. The low positional identity (65%) between the two ADHs, place them into separate classes. We have shown that the ADH5 gene yields two differently processed mRNAs and harbors a gene organization identical to other mammalian ADHs. This is probably due to an alternative splicing in the eighth intron that results in a shorter message missing the ninth exon or a normal message with the expected number of codons. The isolated rat ADH6 cDNA was found to be fused to ADH2 at the 5'-end. The resulting main open reading frame translates into an N-terminally extended polypeptide. In vitro translation results in a polypeptide of about 42 kDa and further, protein was possible to express in COS cells as a fusion product with Green Fluorescent Protein. Both ADH5 and ADH6 show genes and gene products that are processed comparably to other mammalian ADHs and the deduced amino acid sequences indicate a lack of ethanol dehydrogenase activity that probably explains why no corresponding proteins have been isolated. The functionality of these ADHs is therefore still an enigma.     

Reconciling gene and genome duplication events: using multiple nuclear gene families to infer the phylogeny of the aquatic plant family Pontederiaceae

Most plant phylogenetic inference has used DNA sequence data from the plastid genome. This genome represents a single genealogical sample with no recombination among genes, potentially limiting the resolution of evolutionary relationships in some contexts. In contrast, nuclear DNA is inherently more difficult to employ for phylogeny reconstruction because major mutational events in the genome, including polyploidization, gene duplication, and gene extinction can result in homologous gene copies that are difficult to identify as orthologs or paralogs. Gene tree parsimony (GTP) can be used to infer the rooted species tree by fitting gene genealogies to species trees while simultaneously minimizing the estimated number of duplications needed to reconcile conflicts among them. Here, we use GTP for five nuclear gene families and a previously published plastid data set to reconstruct the phylogenetic backbone of the aquatic plant family Pontederiaceae. Plastid-based phylogenetic studies strongly supported extensive paraphyly of Eichhornia (one of the four major genera) but also depicted considerable ambiguity concerning the true root placement for the family. Our results indicate that species trees inferred from the nuclear genes (alone and in combination with the plastid data) are highly congruent with gene trees inferred from plastid data alone. Consideration of optimal and suboptimal gene tree reconciliations place the root of the family at (or near) a branch leading to the rare and locally restricted E. meyeri. We also explore methods to incorporate uncertainty in individual gene trees during reconciliation by considering their individual bootstrap profiles and relate inferred excesses of gene duplication events on individual branches to whole-genome duplication events inferred for the same branches. Our study improves understanding of the phylogenetic history of Pontederiaceae and also demonstrates the utility of GTP for phylogenetic analysis.     

Allelic immunoglobulin VH genes in two mouse strains: possible germline gene recombination.

The nucleotide sequence of two germline immunoglobulin heavy chain variable region (VH) genes of mouse BALB/c origin was determined. These two genes are highly homologous to each other. They both have the unusual codon CCT for proline at position 7, which so far has been found only in a specific set of VH genes, called the NPb family. We show that the two VH genes belong to this set. One of our BALB/c genes, VH124, is more homologous to a C57BL/6 NPb VH gene than to any BALB/c VH gene, and we propose that these two genes are alleles. A comparison of the substitutions between these two genes with published sequences of all other BALB/c and C57BL/6 NPb VH genes reveals evidence for past homologous recombination events between related germline VH genes Homologous recombination may play an important role in the diversification of germline immunoglobulin VH genes.     

High-multiplicity of chitinase genes in Streptomyces coelicolor A3(2).

Six different genes for chitinase from ordered cosmids of the chromosome of Streptomyces coelicolor A3(2) were identified by hybridization, using the chitinase genes from other Streptomyces spp. as probes, and cloned. The genes were sequenced and analyzed. The genes, together with an additional chitinase gene obtained from the data bank, can be classified into either family 18 or family 19 of the glycosyl hydrolase classification. The five chitinases that fall into family 18 show diversity in their multiple domain structures as well as in the amino acid sequences of their catalytic domains. The remaining two chitinases are members of family 19 chitinases, since their C-terminus shares more than 70% identity with the catalytic domain of ChiC of Streptomyces griseus, the sole gene for family 19 chitinase so far found in an organism other than higher plants.     

Nucleotide sequence analysis of the human salivary protein genes HIS1 and HIS2, and evolution of the STATH/HIS gene family

Human histatins are a family of low-M(r), neutral to very basic, histidine-rich salivary polypeptides. They probably function as part of the nonimmune host defense system in the oral cavity. A 39-kb region of DNA containing the HIS1 and HIS2 genes was isolated from two human genomic phage libraries as a series of overlapping clones. The nucleotide sequences of the HIS1 gene and part of the HIS2(1) gene were determined. The transcribed region of HIS1 spans 8.5 kb and contains six exons and five introns. The HIS1 and HIS2(1) genes exhibit 89% overall sequence identity, with exon sequences exhibiting 95% identity. The two loci probably arose by a gene duplication event approximately 15-30 Mya. The HIS1 sequence data were also compared with that of STATH. Human statherin is a low-M(r) acidic phosphoprotein that acts as an inhibitor of precipitation of calcium phosphate salts in the oral cavity. The HIS1 and STATH genes show nearly identical overall gene structures. The HIS1 and STATH loci exhibit 77%-81% sequence identity in intron DNA and 80%-88% sequence identity in noncoding exons but only 38%-43% sequence identity in the protein-coding regions of exons 4 and 5. These unusual data suggest that HIS1, HIS2, and STATH belong to a single gene family exhibiting accelerated evolution between the HIS and STATH coding sequences.      

Whole chloroplast genome comparison of rice,maize, and wheat: implications for chloroplast gene diversification and phylogeny of cereals

The fully sequenced chloroplast genomes of maize (subfamily Panicoideae), rice (subfamily Bambusoideae), and wheat (subfamily Pooideae) provide the unique opportunity to investigate the evolution of chloroplast genes and genomes in the grass family (Poaceae) by whole-genome comparison. Analyses of nucleotide sequence variations in 106 cereal chloroplast genes with tobacco sequences as the outgroup suggested that (1) most of the genic regions of the chloroplast genomes of maize, rice, and wheat have evolved at similar rates; (2) RNA genes have highly conservative evolutionary rates relative to the other genes; (3) photosynthetic genes have been under strong purifying selection; (4) between the three cereals, 14 genes which account for about 28% of the genic region have evolved with heterogeneous nucleotide substitution rates; and (5) rice genes tend to have evolved more slowly than the others at loci where rate heterogeneity exists. Although the mechanism that underlies chloroplast gene diversification is complex, our analyses identified variation in nonsynonymous substitution rates as a genetic force that generates heterogeneity, which is evidence of selection in chloroplast gene diversification at the intrafamilial level. Phylogenetic trees constructed with the variable nucleotide sites of the chloroplast genes place maize basal to the rice-wheat clade, revealing a close relationship between the Bambusoideae and Pooideae.