Concerted evolution: molecular mechanism and biological implications

Concerted evolution of duplicated mitochondrial control regions in three related seabird species

Background

Many population genetic and phylogenetic analyses of mitochondrial DNA (mtDNA) assume that mitochondrial genomes do not undergo recombination. Recently, concerted evolution of duplicated mitochondrial control regions has been documented in a range of taxa. Although the molecular mechanism that facilitates concerted evolution is unknown, all proposed mechanisms involve mtDNA recombination.

Results

Here, we document a duplication of a large region (cytochrome b, tRNA^Thr, tRNA^Pro, ND6, tRNA^Glu and the control region) in the mitochondrial genome of three related seabird species. To investigate the evolution of duplicate control regions, we sequenced both control region copies (CR1 and CR2) from 21 brown (Sula leucogaster), 21 red-footed (S. sula) and 21 blue-footed boobies (S. nebouxii). Phylogenetic analysis suggested that the duplicated control regions are predominantly evolving in concert; however, approximately 51 base pairs at the 5' end of CR1 and CR2 exhibited a discordant phylogenetic signal and appeared to be evolving independently.

Conclusions

Both the structure of the duplicated region and the conflicting phylogenetic signals are remarkably similar to a pattern found in Thalassarche albatrosses, which are united with boobies in a large clade that includes all procellariiform and most pelecaniform seabirds. Therefore we suggest that concerted evolution of duplicated control regions either is taxonomically widespread within seabirds, or that it has evolved many times.

An intergenic non-coding rRNA correlated with expression of the rRNA and frequency of an rRNA single nucleotide polymorphism in lung cancer cells

BACKGROUND

Ribosomal RNA (rRNA) is a central regulator of cell growth and may control cancer development. A cis noncoding rRNA (nc-rRNA) upstream from the 45S rRNA transcription start site has recently been implicated in control of rRNA transcription in mouse fibroblasts. We investigated whether a similar nc-rRNA might be expressed in human cancer epithelial cells, and related to any genomic characteristics.

METHODOLOGY/PRINCIPAL FINDINGS

Using quantitative rRNA measurement, we demonstrated that a nc-rRNA is transcribed in human lung epithelial and lung cancer cells, starting from approximately -1000 nucleotides upstream of the rRNA transcription start site (+1) and extending at least to +203. This nc-rRNA was significantly more abundant in the majority of lung cancer cell lines, relative to a nontransformed lung epithelial cell line. Its abundance correlated negatively with total 45S rRNA in 12 of 13 cell lines (P = 0.014). During sequence analysis from -388 to +306, we observed diverse, frequent intercopy single nucleotide polymorphisms (SNPs) in rRNA, with a frequency greater than predicted by chance at 12 sites. A SNP at +139 (U/C) in the 5' leader sequence varied among the cell lines and correlated negatively with level of the nc-rRNA (P = 0.014). Modelling of the secondary structure of the rRNA 5'-leader sequence indicated a small increase in structural stability due to the +139 U/C SNP and a minor shift in local configuration occurrences.

CONCLUSIONS/SIGNIFICANCE

The results demonstrate occurrence of a sense nc-rRNA in human lung epithelial and cancer cells, and imply a role in regulation of the rRNA gene, which may be affected by a +139 SNP in the 5' leader sequence of the primary rRNA transcript.

Tracing misalignment and unequal crossing over during evolution of the repetitive region of the SM50 gene in sea urchins

The location of misalignment and unequal crossover involved in concerted evolution of tandemly repetitive sequences is difficult to document owing to the homogeneity of sequences that are subject to this process. However, the repetitive domain of the SM50 gene in sea urchins contains variation, within the gene itself, between alleles, and between species that has allowed us to determine where misalignment and unequal crossing over occurred during evolution of this gene. We have therefore analyzed the SM50 repeat regions in a variety of species to determine where recent changes in repeat numbers have occurred, and from this have deduced the mechanisms that lead to these changes. We next tried to determine whether recent misalignment and unequal crossover has produced allelic variation in current populations of sea urchins. We found SM50 alleles within three species that have different numbers of repeats. This marks the first reported documentation of allelic variation in the number of repeats in the SM50 gene. We also show how a single unequal crossover event could have produced the allelic variation. We have found that substitutions and small deletions in the sequences within the repeats can substantially affect how misalignment occurs, resulting in different patterns of repeats after concerted evolution.

Populations, hybrids and the systematic concepts of species and subspecies in Chagas disease triatomine vectors inferred from nuclear ribosomal and mitochondrial DNA

In Chagas disease, triatomine vectors are the main target for control measures because of the absence of effective drugs. The broad usefulness of nuclear rDNA and mtDNA sequences explains why triatomine studies using these markers have increased so pronouncedly in recent years. This indicates the appropriateness of an updated review about these molecular markers, concentrating on aspects useful for research on Chagas disease vectors. A comparative analysis is presented on the efficiency, weight of their different characteristics, limitations and problems of each of the different DNA markers in the light of the results obtained in studies on populations, hybrids, subspecies and species of the subfamily Triatominae. The use of a standardized composite haplotype code nomenclature for both nuclear rDNA and mtDNA markers is strongly encouraged to avoid difficulties in comparative studies. Triatomine aspects related to concerted evolution, microsatellites, minisatellites and insertions/deletions in nuclear rDNA and silent/non-silent mutations, pseudogenes and weaknesses of partial sequences in mtDNA are analysed. Introgression and hybrids, nuclear and mitochondrial DNA strengths, and compared evolutionary rates of nuclear rDNA and mtDNA in triatomines are discussed. Many conclusions are obtained thanks to the availability, for the first time in triatomines, of a complete sequence of a protein-coding mtDNA gene as ND1 from very numerous triatomine species covering from different populations of a species up to members belonging to different tribes. The evolutionary rates of each nuclear rDNA marker and mtDNA marker are analysed by comparison at subspecies level (intrapopulational, interpopulational, between morphs, and between subspecies) and species level (close and distant species of the same genus, species of different genera, and species of different tribes). Weaknesses of mtDNA for systematic-taxonomic purposes detected recently and newly in insects and triatomines, respectively, are discussed in detail. Emphasis is given to taxonomic units and biological entities presenting well-known problematics, both from the systematic-taxonomic and/or epidemiological-control points of view, as well as to molecular situations which can give rise to erroneous conclusions. All these aspects constitute the background on which the key question about the systematic concepts of species and subspecies in triatomines is focused. The global purpose is to facilitate future work on triatomines by highlighting present gaps, how better choice the appropriate markers, and marker aspects which should be taken into account. Key characteristics as alpha, CI and transformation rate matrices ought to be obtained and noted to get appropriate results and allow correct interpretations. The main aim is to offer a baseline for future fundamental research on triatomines and applied research on transmission, epidemiology and control measures related to Chagas disease vectors.

Evolution of spliceosomal snRNA genes in metazoan animals

While studies of the evolutionary histories of protein families are commonplace, little is known on noncoding RNAs beyond microRNAs and some snoRNAs. Here we investigate in detail the evolutionary history of the nine spliceosomal snRNA families (U1, U2, U4, U5, U6, U11, U12, U4atac, and U6atac) across the completely or partially sequenced genomes of metazoan animals. Representatives of the five major spliceosomal snRNAs were found in all genomes. None of the minor splicesomal snRNAs were detected in nematodes or in the shotgun traces of Oikopleura dioica, while in all other animal genomes at most one of them is missing. Although snRNAs are present in multiple copies in most genomes, distinguishable paralogue groups are not stable over long evolutionary times, although they appear independently in several clades. In general, animal snRNA secondary structures are highly conserved, albeit, in particular, U11 and U12 in insects exhibit dramatic variations. An analysis of genomic context of snRNAs reveals that they behave like mobile elements, exhibiting very little syntenic conservation.

Geographically localised bursts of ribosomal DNA mobility in the grasshopper Podisma pedestris

We report extraordinary variation in the number and the chromosomal location of ribosomal DNA (rDNA) arrays within populations of the alpine grasshopper Podisma pedestris; even greater differences were found between populations. The sites were detected by in situ hybridisation of labelled rDNA to chromosomal preparations. The total number of rDNA sites in an individual varied from three to thirteen. In the most extreme case, individuals from populations only 10 km apart had no rDNA loci in common. A survey of the geographical distribution of this variation identified clusters of populations with relatively similar chromosomal distribution of rDNA loci. These clusters correspond to those identified earlier by analysis of rDNA sequences. To explain this geographical clustering, we reconstructed the post-glacial colonisation of the region by assuming that the species' distribution has ascended to its current altitudinal range as the climate warmed. The reconstruction suggests that each cluster is descended from a colonisation route up a different alpine valley. That history would imply rapid establishment of rDNA differences, conceivably during the last 10,000 years since the last glaciation. The proposal for rapid change is consistent with the extensive within-population variation, which indicates that the processes responsible for the change in rDNA's chromosomal location continue to occur at a higher rate. We discuss whether our reconstruction of colonisation routes implies movement of the hybrid zone, which would indicate that a neo-XY sex chromosome system has spread through extant populations.

Characterization of PR-10 genes from eight Betula species and detection of Bet v 1 isoforms in birch pollen

BACKGROUND

Bet v 1 is an important cause of hay fever in northern Europe. Bet v 1 isoforms from the European white birch (Betula pendula) have been investigated extensively, but the allergenic potency of other birch species is unknown. The presence of Bet v 1 and closely related PR-10 genes in the genome was established by amplification and sequencing of alleles from eight birch species that represent the four subgenera within the genus Betula. Q-TOF LC-MSE was applied to identify which PR-10/Bet v 1 genes are actually expressed in pollen and to determine the relative abundances of individual isoforms in the pollen proteome.

RESULTS

All examined birch species contained several PR-10 genes. In total, 134 unique sequences were recovered. Sequences were attributed to different genes or pseudogenes that were, in turn, ordered into seven subfamilies. Five subfamilies were common to all birch species. Genes of two subfamilies were expressed in pollen, while each birch species expressed a mixture of isoforms with at least four different isoforms. Isoforms that were similar to isoforms with a high IgE-reactivity (Bet v 1a = PR-10.01A01) were abundant in all species except B. lenta, while the hypoallergenic isoform Bet v 1d (= PR-10.01B01) was only found in B. pendula and its closest relatives.

CONCLUSION

Q-TOF LC-MSE allows efficient screening of Bet v 1 isoforms by determining the presence and relative abundance of these isoforms in pollen. B. pendula contains a Bet v 1-mixture in which isoforms with a high and low IgE-reactivity are both abundant. With the possible exception of B. lenta, isoforms identical or very similar to those with a high IgE-reactivity were found in the pollen proteome of all examined birch species. Consequently, these species are also predicted to be allergenic with regard to Bet v 1 related allergies.

Natural hybrids in the marine diatom Pseudo-nitzschia pungens (Bacillariophyceae): genetic and morphological evidence

Hybridization between genetically distinguishable taxa provides opportunities for investigating speciation. While hybridization is a common phenomenon in various macro-organisms, natural hybridization among micro-eukaryotes is barely studied. Here we used a nuclear and a chloroplast molecular marker and morphology to demonstrate the presence of natural hybrids between two genetically and morphologically distinct varieties of the marine planktonic diatom Pseudo-nitzschia pungens (vars. pungens and cingulata) in a contact zone in the northeast Pacific. Cloning and sequencing of the rDNA internal transcribed spacer region revealed strains containing ribotypes from both varieties, indicating hybridization. Both varieties were found to also have different chloroplast-encoded rbcL sequences. Hybrid strains were either hetero- or homoplastidial, as demonstrated by denaturing gradient gel electrophoresis, which is in accordance with expectations based on the mode of chloroplast inheritance in Pseudo-nitzschia. While most hybrids are probably first generation, there are also indications for further hybridization. Morphologically, the hybrids resembled var. pungens for most characters rather than having an intermediate morphology. Further research should focus on the hybridization frequency, by assessing the spatial and temporal extent of the contact zone, and hybrid fitness, to determine the amount of gene flow between the two varieties and its evolutionary consequences.

Intergenic exchange, geographic isolation, and the evolution of bioluminescent color for Pyrophorus click beetles

Gene duplication is an evolutionary process in which the emergent property of the whole can become greater and different than the sum of its parts. One potential outcome for gene duplication is for loci to evolve different, yet related functions. In this case, intergenic exchange can shuffle blocks of differentiated nucleotides between paralogues to create new alleles and phenotypes rather than simply homogenize loci. Bioluminescent click beetles in the genus Pyrophorus (Coleoptera: Elateridae) provide an opportunity to explore the creative potential of intergenic exchange for gene family evolution. Pyrophorus beetles bioluminesce different light colors from a pair of dorsal light organs and a ventral light organ. The light organs are under the separate genetic control of dorsal and ventral luciferase loci. Here, we report that intergenic exchange is common between dorsal and ventral loci for beetles from Jamaica (P. plagiophthalamus), the Dominican Republic (P. mellifluous), Belize (P. luscus), and Trinidad (P. noctilucus). We also present evidence that periods of past geographic isolation for beetles on Jamaica, probably acting in concert with selection, built differentiated blocks of substitutions within dorsal and ventral P. plagiophthalamus luciferase loci. Gene flow and intergenic exchange subsequently shuffled these substitutions between dorsal and ventral loci to produce new color phenotypes on Jamaica, including a yellow-green polymorphism. We discuss the possibility of a previously unrecognized emergent evolutionary property of intergenic exchange for luciferase involving cycles of bioluminescent color change related to differences in selective constrains acting on dorsal versus ventral loci. We also explore whether intergenic exchange may commonly create novel variation and the potential for cyclic evolution in other multigene family systems.

Evolutionary origin and genomic organisation of runt-domain containing genes in arthropods

BACKGROUND

Gene clusters, such as the Hox gene cluster, are known to have critical roles in development. In eukaryotes gene clusters arise primarily by tandem gene duplication and divergence. Genes within a cluster are often co-regulated, providing selective pressure to maintain the genome organisation, and this co-regulation can result in temporal or spatial co-linearity of gene expression. It has been previously noted that in Drosophila melanogaster, three of the four runt-domain (RD) containing genes are found in a relatively tight cluster on chromosome 1, raising the possibility of a putative functional RD gene cluster in D. melanogaster.

RESULTS

To investigate the possibility of such a gene cluster, orthologues of the Drosophila melanogaster RD genes were identified in several endopterygotan insects, two exopterygotan insects and two non-insect arthropods. In all insect species four RD genes were identified and orthology was assigned to the Drosophila sequences by phylogenetic analyses. Although four RD genes were found in the crustacean D. pulex, orthology could not be assigned to the insect sequences, indicating independent gene duplications from a single ancestor following the split of the hexapod lineage from the crustacean lineage.In insects, two chromosomal arrangements of these genes was observed; the first a semi-dispersed cluster, such as in Drosophila, where lozenge is separated from the core cluster of three RD genes often by megabases of DNA. The second arrangement was a tight cluster of the four RD genes, such as in Apis mellifera.This genomic organisation, particularly of the three core RD genes, raises the possibility of shared regulatory elements. In situ hybridisation of embryonic expression of the four RD genes in Drosophila melanogaster and the honeybee A. mellifera shows no evidence for either spatial or temporal co-linearity of expression during embryogenesis.

CONCLUSION

All fully sequenced insect genomes contain four RD genes and orthology can be assigned to these genes based on similarity to the D. melanogaster protein sequences. Examination of the genomic organisation of these genes provides evidence for a functional RD gene cluster. RD genes from non-insect arthropods are also clustered, however the lack of orthology between these and insect RD genes suggests this cluster is likely to have resulted from a duplication event independent from that which created the insect RD gene cluster. Analysis of embryonic RD gene expression in two endopterygotan insects, A. mellifera and D. melanogaster, did not show evidence for coordinated gene expression, therefore while the functional significance of this gene cluster remains unknown its maintenance during insect evolution implies some functional significance to the cluster.

Controlled chaos of polymorphic mucins in a metazoan parasite (Schistosoma mansoni) interacting with its invertebrate host (Biomphalaria glabrata)

Invertebrates were long thought to possess only a simple, effective and hence non-adaptive defence system against microbial and parasitic attacks. However, recent studies have shown that invertebrate immunity also relies on immune receptors that diversify (e.g. in echinoderms, insects and mollusks (Biomphalaria glabrata)). Apparently, individual or population-based polymorphism-generating mechanisms exists that permit the survival of invertebrate species exposed to parasites. Consequently, the generally accepted arms race hypothesis predicts that molecular diversity and polymorphism also exist in parasites of invertebrates. We investigated the diversity and polymorphism of parasite molecules (Schistosoma mansoni Polymorphic Mucins, SmPoMucs) that are key factors for the compatibility of schistosomes interacting with their host, the mollusc Biomphalaria glabrata. We have elucidated the complex cascade of mechanisms acting both at the genomic level and during expression that confer polymorphism to SmPoMuc. We show that SmPoMuc is coded by a multi-gene family whose members frequently recombine. We show that these genes are transcribed in an individual-specific manner, and that for each gene, multiple splice variants exist. Finally, we reveal the impact of this polymorphism on the SmPoMuc glycosylation status. Our data support the view that S. mansoni has evolved a complex hierarchical system that efficiently generates a high degree of polymorphism—a “controlled chaos”—based on a relatively low number of genes. This contrasts with protozoan parasites that generate antigenic variation from large sets of genes such as Trypanosoma cruzi, Trypanosoma brucei and Plasmodium falciparum. Our data support the view that the interaction between parasites and their invertebrate hosts are far more complex than previously thought. While most studies in this matter have focused on invertebrate host diversification, we clearly show that diversifying mechanisms also exist on the parasite side of the interaction. Our findings shed new light on how and why invertebrate immunity develops.

Contrary to the traditional view that immunity in invertebrates is limited to non-specific mechanisms, recent studies have shown that they have diverse, specific immune receptors. An example is provided by the FREPs of the mollusk Biomphalaria glabrata, polymorphic members of the immunoglobulin superfamily. This capacity for an individual or population-based polymorphic immune response raises the question of whether a corresponding polymorphism exists in parasites of invertebrates, as would be expected in an “arms race” between host and parasite. We have indeed identified such polymorphic molecules in Schistosoma mansoni, a flatworm parasite of B. glabrata, by comparing two strains of schistosome that are respectively compatible and incompatible with the same mollusk host strain. However, in contrast to antigenic variation in protozoan parasites that is based on an extensive gene repertoire, we show here that a high level of polymorphism in these S. mansoni polymorphic mucins (SmPoMucs) is generated from a low number of genes by a complex cascade of mechanisms, a “controlled chaos”.

Evolution of a new function by degenerative mutation in cephalochordate steroid receptors

Gene duplication is the predominant mechanism for the evolution of new genes. Major existing models of this process assume that duplicate genes are redundant; degenerative mutations in one copy can therefore accumulate close to neutrally, usually leading to loss from the genome. When gene products dimerize or interact with other molecules for their functions, however, degenerative mutations in one copy may produce repressor alleles that inhibit the function of the other and are therefore exposed to selection. Here, we describe the evolution of a duplicate repressor by simple degenerative mutations in the steroid hormone receptors (SRs), a biologically crucial vertebrate gene family. We isolated and characterized the SRs of the cephalochordate Branchiostoma floridae, which diverged from other chordates just after duplication of the ancestral SR. The B. floridae genome contains two SRs: BfER, an ortholog of the vertebrate estrogen receptors, and BfSR, an ortholog of the vertebrate receptors for androgens, progestins, and corticosteroids. BfSR is specifically activated by estrogens and recognizes estrogen response elements (EREs) in DNA; BfER does not activate transcription in response to steroid hormones but binds EREs, where it competitively represses BfSR. The two genes are partially coexpressed, particularly in ovary and testis, suggesting an ancient role in germ cell development. These results corroborate previous findings that the ancestral steroid receptor was estrogen-sensitive and indicate that, after duplication, BfSR retained the ancestral function, while BfER evolved the capacity to negatively regulate BfSR. Either of two historical mutations that occurred during BfER evolution is sufficient to generate a competitive repressor. Our findings suggest that after duplication of genes whose functions depend on specific molecular interactions, high-probability degenerative mutations can yield novel functions, which are then exposed to positive or negative selection; in either case, the probability of neofunctionalization relative to gene loss is increased compared to existing models.

Evolution of genes and repeats in the Nimrod superfamily

The recently identified Nimrod superfamily is characterized by the presence of a special type of EGF repeat, the NIM repeat, located right after a typical CCXGY/W amino acid motif. On the basis of structural features, nimrod genes can be divided into three types. The proteins encoded by Draper-type genes have an EMI domain at the N-terminal part and only one copy of the NIM motif, followed by a variable number of EGF-like repeats. The products of Nimrod B-type and Nimrod C-type genes (including the eater gene) have different kinds of N-terminal domains, and lack EGF-like repeats but contain a variable number of NIM repeats. Draper and Nimrod C-type (but not Nimrod B-type) proteins carry a transmembrane domain. Several members of the superfamily were claimed to function as receptors in phagocytosis and/or binding of bacteria, which indicates an important role in the cellular immunity and the elimination of apoptotic cells. In this paper, the evolution of the Nimrod superfamily is studied with various methods on the level of genes and repeats. A hypothesis is presented in which the NIM repeat, along with the EMI domain, emerged by structural reorganizations at the end of an EGF-like repeat chain, suggesting a mechanism for the formation of novel types of repeats. The analyses revealed diverse evolutionary patterns in the sequences containing multiple NIM repeats. Although in the Nimrod B and Nimrod C proteins show characteristics of independent evolution, many internal NIM repeats in Eater sequences seem to have undergone concerted evolution. An analysis of the nimrod genes has been performed using phylogenetic and other methods and an evolutionary scenario of the origin and diversification of the Nimrod superfamily is proposed. Our study presents an intriguing example how the evolution of multigene families may contribute to the complexity of the innate immune response.

Culture-independent methods for identifying microbial communities in cheese

This review focuses on the culture-independent methods available for the description of both bacterial and fungal communities in cheese. Important steps of the culture-independent strategy, which relies on bulk DNA extraction from cheese and polymerase chain reaction (PCR) amplification of selected sequences, are discussed. We critically evaluate the identification techniques already used for monitoring microbial communities in cheese, including PCR-denaturing gradient gel electrophoresis (PCR-DGGE), PCR-temporal temperature gradient gel electrophoresis (PCR-TTGE) or single-strand conformation polymorphism-PCR (SSCP-PCR) as well as some other techniques that remain to be adapted to the study of cheese communities. Further, our analysis draws attention to the lack of data available on suitable DNA sequences for identifying fungal communities in cheese and proposes some potential DNA targets.

Human U2 snRNA genes exhibit a persistently open transcriptional state and promoter disassembly at metaphase

In mammals, small multigene families generate spliceosomal U snRNAs that are nearly as abundant as rRNA. Using the tandemly repeated human U2 genes as a model, we show by footprinting with DNase I and permanganate that nearly all sequences between the enhancer-like distal sequence element and the initiation site are protected during interphase whereas the upstream half of the U2 snRNA coding region is exposed. We also show by chromatin immunoprecipitation that the SNAPc complex, which binds the TATA-like proximal sequence element, is removed at metaphase but remains bound under conditions that induce locus-specific metaphase fragility of the U2 genes, such as loss of CSB, BRCA1, or BRCA2 function, treatment with actinomycin D, or overexpression of the tetrameric p53 C terminus. We propose that the U2 snRNA promoter establishes a persistently open state to facilitate rapid reinitiation and perhaps also to bypass TFIIH-dependent promoter melting; this open state would then be disassembled to allow metaphase chromatin condensation.

Patterns of genetic structure among Hawaiian corals of the genus Pocillopora yield clusters of individuals that are compatible with morphology

Six variable sequence markers are developed and analyzed to find out species boundaries in Hawaiian corals of the genus Pocillopora: the putative mitochondrial control region; a recently discovered, hypervariable mitochondrial open reading frame; the internal transcribed spacer 2 (ITS2), located in the nuclear ribosomal DNA; three nuclear introns of calmodulin, elongation factor-1alpha and the ATP synthase beta subunit. Using the first two markers, we identify five distinct mitochondrial lineages and these lineages are compatible with morphology. The situation is more complex with nuclear markers since more than two haplotypes are observed in some individuals. To detect clusters of individuals, haplotype networks are constructed with additional connections drawn between co-occurring haplotypes to delineate potential fields for recombination: few clusters of nuclear haplotypes are found to correspond to clusters of individuals, but those that are detected (mostly in the ITS2 dataset) are also compatible with morphology.

T2prhd: a tool to study the patterns of repeat evolution

Background

The models developed to characterize the evolution of multigene families (such as the birth-and-death and the concerted models) have also been applied on the level of sequence repeats inside a gene/protein. Phylogenetic reconstruction is the method of choice to study the evolution of gene families and also sequence repeats in the light of these models. The characterization of the gene family evolution in view of the evolutionary models is done by the evaluation of the clustering of the sequences with the originating loci in mind. As the locus represents positional information, it is straightforward that in the case of the repeats the exact position in the sequence should be used, as the simple numbering according to repeat order can be misleading.

Results

We have developed a novel rapid visual approach to study repeat evolution, that takes into account the exact repeat position in a sequence. The "pairwise repeat homology diagram" visualizes sequence repeats detected by a profile HMM in a pair of sequences and highlights their homology relations inferred by a phylogenetic tree. The method is implemented in a Perl script (t2prhd) available for downloading at and is also accessible as an online tool at . The power of the method is demonstrated on the EGF-like and fibronectin-III-like (Fn-III) domain repeats of three selected mammalian Tenascin sequences.

Conclusion

Although pairwise repeat homology diagrams do not carry all the information provided by the phylogenetic tree, they allow a rapid and intuitive assessment of repeat evolution. We believe, that t2prhd is a helpful tool with which to study the pattern of repeat evolution. This method can be particularly useful in cases of large datasets (such as large gene families), as the command line interface makes it possible to automate the generation of pairwise repeat homology diagrams with the aid of scripts.

5S ribosomal RNA genes in six species of Mediterranean grey mullets: genomic organization and phylogenetic inference

This paper describes a study of the 5S ribosomal RNA genes (5S rDNA) in a group of 6 species belonging to 4 genera of Mugilidae. In these 6 species, the relatively short 5S rDNA repeat units, generated by PCR and ranging in size from 219 to 257 bp, show a high level of intragenomic homogeneity of both coding and spacer regions (NTS-I). Phylogenetic reconstructions based on this data set highlight the greater phylogenetic and genetic diversity of Mugil cephalus and Oedalechilus labeo compared with the genera Liza and Chelon. Comparative sequence analysis revealed significant conservation of the short 5S rDNA repeat units across Chelon and Liza. Moreover, a second size class of 5S rDNA repeat units, ranging from roughly 800 to 1100 bp, was produced in the Liza and Chelon samples. Only short 5S rDNA repeat units were found in M. cephalus and O. labeo. The sequences of the long 5S rDNA repeat units, obtained in Chelon labrosus and Liza ramada, differ owing to the presence of 2 large insertion/deletions (indels) in the spacers (NTS-II) and show considerable sequence identity with NTS-I spacers. Interspecific sequence variation of NTS-II spacers, excluding the indels, is low. Southern-blot hybridization patterns suggest an intermixed arrangement of short and long repeat units within a single chromosome locus.

Intraspecific concerted evolution of the rDNA ITS1 in Anopheles farauti sensu stricto (Diptera: Culicidae) reveals recent patterns of population structure

We examined the intraindividual variation present in the first ribosomal internal transcribed spacer (ITS1) of Anopheles farauti to determine the level of divergence among populations for this important malarial vector. We isolated 187 clones from 70 individuals and found regional variation among four internal tandem repeats. The data were partitioned prior to analysis given the presence of a paralogous ITS2 sequence, called the 5'-subrepeat, inserted in the ITS1 of most clones. A high level of homogenization and population differentiation was observed for this repeat, which indicates a higher rate of turnover relative to the adjacent 'core' region. Bayesian analysis was performed using several substitutional models on both a combined and a partitioned data set. On the whole, the ITS1 phylogeny and geographic origin of the samples appear to be congruent. Some interesting exceptions indicate the spread of variant repeats between populations and the retention of ancestral polymorphism. Our data clearly demonstrate concerted evolution at the intraspecific level despite intraindividual variation and a complex internal repeat structure from a species that occupies a continuous coastal distribution. A high rate of genomic turnover in combination with a high level of sequence divergence appears to be a major factor leading to its concerted evolution within these populations.

Intra-specific and intra-sporocarp ITS variation of ectomycorrhizal fungi as assessed by rDNA sequencing of sporocarps and pooled ectomycorrhizal roots from a Quercus woodland

The Internal Transcribed Spacer (ITS) regions of ribosomal DNA are widely used as markers for phylogenetic analyses and environmental sampling from a variety of organisms including fungi, plants, and animals. In theory, concerted evolution homogenizes multicopy genes so that little or no variation exists within populations or individuals. However, contrary to theory, ITS variation has been confirmed in populations and individuals from a diverse range of eukaryotes. The presence of intraspecific and intra-individual variation in multicopy genes has important implications for ecological and phylogenetic studies, yet relatively little is known about natural variation of these genes, particularly at the community level. In this study, we examined intraspecific and intra-sporocarp ITS variation by DNA sequencing from sporocarps and pooled roots from 68 species of ectomycorrhizal fungi collected at a single site in a Quercus woodland. We detected ITS variation in 27 species, roughly 40% of the taxa examined. Although intraspecific ITS variation was generally low (0.16-2.85%, mean = 0.74%), it was widespread within this fungal community. We detected ITS variation in both sporocarps and ectomycorrhizal roots, and variation was present within species of Ascomycota and Basidiomycota, two distantly related lineages within the Fungi. We discuss the implications of such widespread ITS variability with special reference to DNA-based environmental sampling from diverse fungal communities.

Euglena light-harvesting complexes are encoded by multifarious polyprotein mRNAs that evolve in concert

Light-harvesting complexes (LHCs) are a superfamily of chlorophyll- and carotenoid-binding proteins that are responsible for the capture of light energy and its transfer to the photosynthetic reaction centers. Unlike those of most eukaryotes, the LHCs of Euglena gracilis are translated from large mRNAs, producing polyprotein precursors consisting of multiple concatenated LHC subunits that are separated by conserved decapeptide linkers. These precursors are posttranslationally targeted to the chloroplast and cleaved into individual proteins. We analyzed expressed sequence tags from Euglena to further characterize the structural features of the LHC polyprotein-coding genes and to examine the evolution of this multigene family. Of the 19 different LHC transcriptional units we detected, 17 encoded polyproteins composed of both tandem and nontandem repeats of LHC subunits; organizations that likely occurred through unequal crossing-over. Of the 2 nonpolyprotein-encoding LHC transcripts detected, 1 evolved from the truncation of a polyprotein-coding gene. Duplication of LHC polyprotein-coding genes was particularly important in the LHCI gene family where multiple paralogous sequences were detected. Intriguingly, several of the individual LHC-coding subunits both within and between transcriptional units appeared to be evolving in concert, suggesting that gene conversion has been a significant mechanism for LHC evolution in Euglena.

The second internal transcribed spacer of nuclear ribosomal DNA as a tool for Latin American anopheline taxonomy - a critical review

Among the molecular markers commonly used for mosquito taxonomy, the internal transcribed spacer 2 (ITS2) of the ribosomal DNA is useful for distinguishing among closely-related species. Here we review 178 GenBank accession numbers matching ITS2 sequences of Latin American anophelines. Among those, we found 105 unique sequences corresponding to 35 species. Overall the ITS2 sequences distinguish anopheline species, however, information on intraspecific and geographic variations is scarce. Intraspecific variations ranged from 0.2% to 19% and our analysis indicates that misidentification and/or sequencing errors could be responsible for some of the high values of divergence. Research in Latin American malaria vector taxonomy profited from molecular data provided by single or few field capture mosquitoes. However we propose that caution should be taken and minimum requirements considered in the design of additional studies. Future studies in this field should consider that: (1) voucher specimens, assigned to the DNA sequences, need to be deposited in collections, (2) intraspecific variations should be thoroughly evaluated, (3) ITS2 and other molecular markers, considered as a group, will provide more reliable information, (4) biological data about vector populations are missing and should be prioritized, (5) the molecular markers are most powerful when coupled with traditional taxonomic tools.

Intra-isolate heterogeneity of the ITS region of rDNA in Pythium helicoides

Heterogeneity of the rDNA ITS region in Pythium helicoides and the phylogenetic relationship between P. helicoides and closely related species were investigated. In PCR-RFLP analysis of the rDNA ITS region of six P. helicoides isolates investigated, including the type culture, intraspecific variation was found at the HhaI site. The total length of fragments was longer than before cutting, indicating sequence heterogeneity within isolates. Digestion of the cloned rDNA ITS region derived from seven isolates with HhaI revealed polymorphisms among and within single zoospore isolates, and variability of the region was also present among the clones derived from the same isolate. To test whether the rDNA ITS region of closely related species and other regions in the genome of P. helicoides are also variable, the rDNA ITS region of P. ultimum and the cytochrome oxydase II (cox II) gene encoded in mitochondria were sequenced. P. ultimum had little variation in the rDNA ITS region. The cox II gene sequences of both species revealed only a low intraspecific variability and no intra-isolate variation. In the phylogenic tree based on the rDNA ITS sequences, all clones of P. helicoides formed one large clade that was distinct from the clades comprising morphologically similar species, such as P. oedochilum and P. ostracodes, and was closely related to P. chamaehyphon rather than the other species.

An evaluation of LSU rDNA D1-D2 sequences for their use in species identification

Background

Identification of species via DNA sequences is the basis for DNA taxonomy and DNA barcoding. Currently there is a strong focus on using a mitochondrial marker for this purpose, in particular a fragment from the cytochrome oxidase I gene (COI). While there is ample evidence that this marker is indeed suitable across a broad taxonomic range to delineate species, it has also become clear that a complementation by a nuclear marker system could be advantageous. Ribosomal RNA genes could be suitable for this purpose, because of their global occurrence and the possibility to design universal primers. However, it has so far been assumed that these genes are too highly conserved to allow resolution at, or even beyond the species level. On the other hand, it is known that ribosomal gene regions harbour also highly divergent parts. We explore here the information content of two adjacent divergence regions of the large subunit ribosomal gene, the D1-D2 region.

Results

Universal primers were designed to amplify the D1-D2 region from all metazoa. We show that amplification products in the size between 800–1300 bp can be obtained across a broad range of animal taxa, provided some optimizations of the PCR procedure are implemented. Although the ribosomal genes occur in multiple copies in the genomes, we find generally very little intra-individual polymorphism (<< 0.1% on average) indicating that concerted evolution is very effective in most cases. Studies in two fish taxa (genus Cottus and genus Aphyosemion) show that the D1-D2 LSU sequence can resolve even very closely related species with the same fidelity as COI sequences. In one case we can even show that a mitochondrial transfer must have occurred, since the nuclear sequence confirms the taxonomic assignment, while the mitochondrial sequence would have led to the wrong classification. We have further explored whether hybrids between species can be detected with the nuclear sequence and we show for a test case of natural hybrids among cyprinid fish species (Alburnus alburnus and Rutilus rutilus) that this is indeed possible.

Conclusion

The D1-D2 LSU region is a suitable marker region for applications in DNA based species identification and should be considered to be routinely used as a marker complementing broad scale studies based on mitochondrial markers.

Comparative genomic analysis of the mouse and rat amylase multigene family

The rat and mouse amylase gene families were characterized using sequence data from the UCSC genome assembly. We found that the rat genome contains one amylase-1 and two amylase-2 genes, lying close to one another on the same chromosome. Detailed analysis revealed at least six additional amylase pseudogenes in the rat genome in the region adjacent to the amylase-2 genes. In contrast, the mouse has one amylase-1 gene and five amylase-2 genes; the latter are tandemly and systematically arranged on the same chromosome and were generated by segmental duplication. Detailed analysis revealed that the mouse has two amylase pseudogenes, located 5' to the five amylase-2 segments. Thus, the amylase genes of mouse and rat tend to be amplified; the sequences of some of them are fixed while others have become pseudogenes during evolution. This is the second report of amylase genomic organization in mammals and the first in the rodents.

On the origin of the sweet-smelling Parma violet cultivars (Violaceae): wide intraspecific hybridization, sterility, and sexual reproduction

Parma violets are reputed for their double, fragrant flowers and have been cultivated for centuries in Europe. However, due to a rather atypical morphology their taxonomic affinity has not been clarified. Authors have proposed an origin from three possible species, Viola alba, V. odorata, or V. suavis, or a hybrid origin. Using both ITS sequence variation and allozyme variation in 14 putative loci, we showed that the Parma violet cultivars have their origin within Viola alba and that they are best included in the Mediterranean subsp. dehnhardtii. There is no trace of interspecific hybridization. However, the cultivars appear to have a single origin in a wide hybrid within V. alba, involving parental plants from the eastern and western Mediterranean region; historical literature sources seem to indicate Turkey and Italy, respectively. The Parma violet cultivars possess high levels of allozyme heterozygosity and to some extent also within-individual ITS sequence variation. Losses of heterozygosity and within-individual ITS sequence variation in some of the cultivars indicate subsequent rare events of sexual reproduction, presumably through cleistogamous seed set. We unambiguously identify the closest wild relative of this group of cultivars, allowing growers to develop new selection procedures, and show a peculiar molecular process associated with human selection.

Multiple rRNA Variants in a Single Spore of the Microsporidian Nosema bombi

To understand the source of the multiple DNA sequence variants of Nosema bombi ribosomal RNA (rRNA) found in a single bumble bee host, we PCR amplified, cloned, and sequenced the partial rRNA gene from 125 clones, which were derived from four out of 46 spores individually isolated from a single host by laser microdissection. At least two rRNA variants, characterized by either (GTTT)(2) or (GTTT)(3) repeat units within the internal transcribed spacer (ITS) region, were found per spore in approximately equal proportions, variants which were also found in approximately equal proportions in 55 clones of the two DNA extracts of multiple spores from the same host. Firstly, we demonstrate for the first time that DNA sequences can be obtained from single-binucleate microsporidia. Secondly, it appears that concerted evolution has not homogenized the sequences of all rRNA copies within a single N. bombi spore or even within a single nucleus. We thereby demonstrate unequivocally that two or more rRNA sequence variants exist per N. bombi spore, and urge caution in the use of multicopy rRNA genes for population genetic and phylogenetic analysis of this and other Microsporidia unless homologous copies can be reliably typed.

An unusual mode of concerted evolution of the EGF-TM7 receptor chimera EMR2

The epidermal growth factor (EGF)-TM7 receptors CD97, EMR1, EMR2, EMR3, and EMR4 form a group of adhesion class heptahelical molecules predominantly expressed by cells of the immune system. These receptors bind cellular ligands through EGF-like domains, localized N-terminal to a large extracellular region. Remarkably, EMR2 possesses a chimeric structure with a seven-span transmembrane (TM7) region most related to EMR3 and an EGF domain region nearly identical to CD97. By comparing EGF-TM7 receptors in primates and dogs, we identified an intriguing pattern of concerted evolution, apparently mediated by gene conversion, among EMR2 and the oppositely orientated and physically adjacent genes CD97 and EMR3. This concerted evolution has continuously maintained the chimeric structure of EMR2 since early mammal radiation. Most highly conserved between EMR2 and CD97 is the fourth EGF domain, which mediates binding to chondroitin sulfate, a ligand specificity shared by both receptors. Another ligand, CD55, is bound effectively only by CD97. We show that different molecular mechanisms (mutations vs. alternative splicing) prevent CD55 binding by EMR2 in hominoids. Our findings illustrate how various and partially opposing evolutionary events have shaped the structure and ligand specificity of a modern mammalian gene family.

The evolution of hexapod engrailed-family genes: evidence for conservation and concerted evolution

Phylogenetic analyses imply that multiple engrailed-family gene duplications occurred during hexapod evolution, a view supported by previous reports of only a single engrailed-family gene in members of the grasshopper genus Schistocerca and in the beetle Tribolium castaneum. Here, we report the cloning of a second engrailed-family gene from Schistocerca gregaria and present evidence for two engrailed-family genes from four additional hexapod species. We also report the existence of a second engrailed-family gene in the Tribolium genome. We suggest that the engrailed and invected genes of Drosophila melanogaster have existed as a conserved gene cassette throughout holometabolous insect evolution. In total 11 phylogenetically diverse hexapod orders are now known to contain species that possess two engrailed-family paralogues, with in each case only one paralogue encoding the RS-motif, a characteristic feature of holometabolous insect invected proteins. We propose that the homeoboxes of hexapod engrailed-family paralogues are evolving in a concerted fashion, resulting in gene trees that overestimate the frequency of gene duplication. We present new phylogenetic analyses using non-homeodomain amino acid sequence that support this view. The S. gregaria engrailed-family paralogues provide strong evidence that concerted evolution might in part be explained by recurrent gene conversion. Finally, we hypothesize that the RS-motif is part of a serine-rich domain targeted for phosphorylation.

Sequence comparison of distal and proximal ribosomal DNA arrays in rice (Oryza sativa L.) chromosome 9S and analysis of their flanking regions

Rice (Oryza sativa ssp. japonica cv. Nipponbare) harbors a ribosomal RNA gene (rDNA) cluster in the nucleolar-organizing region at the telomeric end of the short arm of chromosome 9. We isolated and sequenced two genomic clones carrying rice rDNA fragments from this region. The rice rDNA repeat units could be classified into three types based on length, which ranged from 7,928 to 8,934 bp. This variation was due to polymorphism in the number of 254-bp subrepeats in the intergenic spacer (IGS). Polymerase chain reaction (PCR) analysis suggested that the rDNA units in rice vary widely in length and that the copy number of the subrepeats in the IGS ranges from 1 to 12 in the rice genome. PCR and Southern blot analyses showed that most rDNA units have three intact and one truncated copies of the subrepeats in the IGS, and distal (telomere-side) rDNA units have more subrepeats than do proximal (centromere-side) ones. Both genomic clones we studied contained rDNA-flanking DNA sequences of either telomeric repeats (5'-TTTAGGG-3') or a chromosome-specific region, suggesting that they were derived from the distal or proximal end, respectively, of the rDNA cluster. A similarity search indicated that retrotransposons appeared more frequently in a 500-kb portion of the proximal rDNA-flanking region than in other subtelomeric regions or sequenced regions of the genome. This study reveals the repetitive nature of the telomeric end of the short arm of chromosome 9, which consists of telomeric repeats, an rDNA array, and a retrotransposon-rich chromosomal region.

Unique hybrids between the fungal pathogens Cryptococcus neoformans and Cryptococcus gattii

Cryptococcus neoformans and Cryptococcus gattii are yeasts that cause meningoencephalitis, but that differ in host range and geographical distribution. Cryptococcus neoformans occurs world-wide and mostly infects immunocompromised patients, whereas C. gattii occurs mainly in (sub)tropical regions and infects healthy individuals. Anomalous C. neoformans strains were isolated from patients. These strains were found to be monokaryotic, and diploid or aneuploid. Amplified Fragment Length Polymorphism (AFLP) and sequence analyses indicated that AFLP genotypes 2 (C. neoformans) and 4 (C. gattii) were present. The strains were serologically BD. Mating- and serotype-specific PCR reactions showed that the strains were MATa-serotype D/MATalpha-serotype B. This study is the first to describe naturally occurring hybrids between C. neoformans and C. gattii.

Cytogenetics of a new cytotype of African Mus (subgenus Nannomys) minutoides (Rodentia, Muridae) from Kenya: C- and G- banding and distribution of (TTAGGG)n telomeric sequences

We present the results of a cytogenetic study on Mus (Nannomys) minutoides from Kenya by means of C- and G- banding and in-situ fluorescence hybridization (FISH) to localize the telomeric sequences. The karyotype is characterized by the occurrence of several Rb chromosomes Rb(1.X), Rb(1.Y). Rb(2.17), Rb(3.13), Rb(4.10), Rb(5.11), Rb(6.7), Rb(8.12), not previously described for this species. This finding suggests a high level of chromosomal diversification, which means it is possible to consider this cytotype as a new, well-differentiated, chromosomal lineage within the subgenus. The C-banding of the metaphases illustrated conspicuous blocks of centromeric heterochromatin at the paracentromeric regions of all telocentric chromosomes. Centromeric heterochromatin is not visible on all biarmed chromosomes. Following hybridization with telomeric probes, bright interstitial telomeric sequence (ITS) fluorescence signals are evident at the pericentromeric area of all Rb chromosomes, with the exception of Rb(2.17). Considering the localization of the C-positive heterochromatin and of the telomeric sequences, the events leading to the Kenyan cytotype from an all-telocentric condition probably included two steps: first, fusion without loss of heterochromatin and pericentromeric telomeric sequences; second, the reduction of the C-positive satellite DNA followed by the amplification of telomeric sequences in the C-negative paracentromeric region of Rb chromosomes. The presence of a single Rb(2.17) without ITS indicates possible variations of this mechanism.

Genetic map-based analysis of genome structure in the homosporous fern Ceratopteris richardii

Homosporous ferns have extremely high chromosome numbers relative to flowering plants, but the species with the lowest chromosome numbers show gene expression patterns typical of diploid organisms, suggesting that they may be diploidized ancient polyploids. To investigate the role of polyploidy in fern genome evolution, and to provide permanent genetic resources for this neglected group, we constructed a high-resolution genetic linkage map of the homosporous fern model species, Ceratopteris richardii (n = 39). Linkage map construction employed 488 doubled haploid lines (DHLs) that were genotyped for 368 RFLP, 358 AFLP, and 3 isozyme markers. Forty-one linkage groups were recovered, with average spacing between markers of 3.18 cM. Most loci (approximately 76%) are duplicated and most duplicates occur on different linkage groups, indicating that as in other eukaryotic genomes, gene duplication plays a prominent role in shaping the architecture of fern genomes. Although past polyploidization is a potential mechanism for the observed abundance of gene duplicates, a wide range in the number of gene duplicates as well as the absence of large syntenic regions consisting of duplicated gene copies implies that small-scale duplications may be the primary mode of gene duplication in C. richardii. Alternatively, evidence of past polyploidization(s) may be masked by extensive chromosomal rearrangements as well as smaller-scale duplications and deletions following polyploidization(s).

Repetitive DNA in tropical tasar silkworm Antheraea mylitta

Antheraea mylitta is an endemic insect species producing the world famous tasar silk. Its populations occupying different ecological and geographical regions show certain degree of phenotypic variability for which they are known as 'eco-races'. In order to understand the genetic variability and phylogenetic relationship among the different eco-races we characterized a repetitive TaqI genomic DNA fragment as a genetic marker. The sequence analysis and Southern hybridization show the repetitive nature of TaqI DNA fragment, designated as A. mylitta TaqI family repeat, AmTFR. The PCR amplification of AmTFR reveals its presence in all the tested eco-races of tasar silkworm and some other silk producing insects. The AmTFR is evenly distributed in all 31 meiotic metaphase I chromosomes as observed by fluorescent in situ hybridization. The AmTFR based phylogenetic analysis of the eco-races is not congruent with the morphological variations and their geographical distribution.

The orthology of HLA-E and H2-Qa1 is hidden by their concerted evolution with other MHC class I molecules

BACKGROUND

Whether MHC molecules undergo concerted evolution or not has been the subject of a long-standing debate.

RESULTS

By comparing sequences of eight functional homologues of HLA-E from primates and rodents with those of MHC class Ia molecules from the same eight species, we find that different portions of MHC class I molecules undergo different patterns of evolution. By focusing our analyses sequentially on these various portions, we have obtained clear evidence for concerted evolution of MHC class I molecules, suggesting the occurrence of extensive interallelic and intergenic exchanges. Intra-species homogenisation of sequences is particularly noticeable at the level of exon 4, which codes for the alpha3 domain, but our results suggest that homogenisation also concerns certain residues of the alpha1-alpha2 codomain that lie outside the antigen recognition site.

CONCLUSION

A model is presented in which Darwinian selective pressures due to pathogens could, at the same time, favour diversification of MHC class Ia molecules and promote concerted evolution of separate loci by spreading advantageous motifs arising by mutations in individual MHC molecules to other alleles and to other loci of the MHC region. This would also allow MHC molecules to co-evolve with the proteins with which they interact to fulfil their functions of antigen presentation and regulation of NK cell activity. One of the raisons d'être of the MHC may therefore be to favour at the same time both diversification of MHC class Ia molecules and homogenisation of the whole pool of MHC class I molecules (Ia and Ib) involved in antigen presentation.

REVIEWERS

This article was reviewed by Stephan Beck, Lutz Walter and Pierre Pontarotti.