Conservation and divergence in multigene families: alternatives to selection and drift

Identification and characterization of a new family of long satellite DNA, specific of true toads (Anura, Amphibia, Bufonidae)

Amphibians have some of the most variable genome sizes among vertebrates. Genome size variation has been attributed to repetitive and noncoding DNA, including satellite repeats, transposable elements, introns, and nuclear insertions of viral and organelle DNA. In vertebrates, satellite DNAs have been widely described in mammals, but few molecular studies have been carried out in amphibians. Here, we provide a detailed characterization of a new family of satellite DNA, present in all 15 examined species of the family Bufonidae. Southern-blot analysis and PCR reveal that this satellite is formed by monomers of 807 bp, is organized in tandem arrays, and has an AT-content of 57.4%. Phylogenetic analyses show that most clades exhibit species-specific variances, indicating that this satellite DNA has evolved by concerted evolution. The homogenization/fixation process is heterogeneous in Bufonidae, where the genera Bufo and Bufotes do not show species-specific differences, while populations from Rhinella marina exhibit population-specific changes. Additionally, variants of this satellite DNA have been identified in Duttaphrynus melanostictus and R. marina, supporting the 'library hypothesis' (a set, 'library', of satellite DNAs is shared by a species group). Physical mapping in Bufo bufo, Bufo spinosus, Epidalea calamita and Bufotes viridis provides evidence that this repetitive DNA is not dispersed in the karyotype, but accumulated in pericentromeric regions of some chromosomal pairs. This location, together with its presence in the transcriptomes of bufonids, could indicate a role in centromere function or heterochromatin formation and maintenance.

Molecular Dynamics and Evolution of Centromeres in the Genus Equus

The centromere is the chromosomal locus essential for proper chromosome segregation. While the centromeric function is well conserved and epigenetically specified, centromeric DNA sequences are typically composed of satellite DNA and represent the most rapidly evolving sequences in eukaryotic genomes. The presence of satellite sequences at centromeres hampered the comprehensive molecular analysis of these enigmatic loci. The discovery of functional centromeres completely devoid of satellite repetitions and fixed in some animal and plant species represented a turning point in centromere biology, definitively proving the epigenetic nature of the centromere. The first satellite-free centromere, fixed in a vertebrate species, was discovered in the horse. Later, an extraordinary number of satellite-free neocentromeres had been discovered in other species of the genus Equus, which remains the only mammalian genus with numerous satellite-free centromeres described thus far. These neocentromeres arose recently during evolution and are caught in a stage of incomplete maturation. Their presence made the equids a unique model for investigating, at molecular level, the minimal requirements for centromere seeding and evolution. This model system provided new insights on how centromeres are established and transmitted to the progeny and on the role of satellite DNA in different aspects of centromere biology.

Dissecting the Satellite DNA Landscape in Three Cactophilic Drosophila Sequenced Genomes

Eukaryote genomes are replete with repetitive DNAs. This class includes tandemly repeated satellite DNAs (satDNA) which are among the most abundant, fast evolving (yet poorly studied) genomic components. Here, we used high-throughput sequencing data from three cactophilic Drosophila species, D. buzzatii, D. seriema, and D. mojavensis, to access and study their whole satDNA landscape. In total, the RepeatExplorer software identified five satDNAs, three previously described (pBuM, DBC-150 and CDSTR198) and two novel ones (CDSTR138 and CDSTR130). Only pBuM is shared among all three species. The satDNA repeat length falls within only two classes, between 130 and 200 bp or between 340 and 390 bp. FISH on metaphase and polytene chromosomes revealed the presence of satDNA arrays in at least one of the following genomic compartments: centromeric, telomeric, subtelomeric, or dispersed along euchromatin. The chromosomal distribution ranges from a single chromosome to almost all chromosomes of the complement. Fiber-FISH and sequence analysis of contigs revealed interspersion between pBuM and CDSTR130 in the microchromosomes of D. mojavensis Phylogenetic analyses showed that the pBuM satDNA underwent concerted evolution at both interspecific and intraspecific levels. Based on RNA-seq data, we found transcription activity for pBuM (in D. mojavensis) and CDSTR198 (in D. buzzatii) in all five analyzed developmental stages, most notably in pupae and adult males. Our data revealed that cactophilic Drosophila present the lowest amount of satDNAs (1.9-2.9%) within the Drosophila genus reported so far. We discuss how our findings on the satDNA location, abundance, organization, and transcription activity may be related to functional aspects.

Positive and strongly relaxed purifying selection drive the evolution of repeats in proteins

Protein repeats are considered hotspots of protein evolution, associated with acquisition of new functions and novel phenotypic traits, including disease. Paradoxically, however, repeats are often strongly conserved through long spans of evolution. To resolve this conundrum, it is necessary to directly compare paralogous (horizontal) evolution of repeats within proteins with their orthologous (vertical) evolution through speciation. Here we develop a rigorous methodology to identify highly periodic repeats with significant sequence similarity, for which evolutionary rates and selection (dN/dS) can be estimated, and systematically characterize their evolution. We show that horizontal evolution of repeats is markedly accelerated compared with their divergence from orthologues in closely related species. This observation is universal across the diversity of life forms and implies a biphasic evolutionary regime whereby new copies experience rapid functional divergence under combined effects of strongly relaxed purifying selection and positive selection, followed by fixation and conservation of each individual repeat.

Significance of satellite DNA revealed by conservation of a widespread repeat DNA sequence among angiosperms

The analysis of plant genome structure and evolution requires comprehensive characterization of repetitive sequences that make up the majority of plant nuclear DNA. In the present study, we analyzed the nature of pCtKpnI-I and pCtKpnI-II tandem repeated sequences, reported earlier in Carthamus tinctorius. Interestingly, homolog of pCtKpnI-I repeat sequence was also found to be present in widely divergent families of angiosperms. pCtKpnI-I showed high sequence similarity but low copy number among various taxa of different families of angiosperms analyzed. In comparison, pCtKpnI-II was specific to the genus Carthamus and was not present in any other taxa analyzed. The molecular structure of pCtKpnI-I was analyzed in various unrelated taxa of angiosperms to decipher the evolutionary conserved nature of the sequence and its possible functional role.

Molecular organization of the 25S-18S rDNA IGS of Fagus sylvatica and Quercus suber: a comparative analysis

The 35S ribosomal DNA (rDNA) units, repeated in tandem at one or more chromosomal loci, are separated by an intergenic spacer (IGS) containing functional elements involved in the regulation of transcription of downstream rRNA genes. In the present work, we have compared the IGS molecular organizations in two divergent species of Fagaceae, Fagus sylvatica and Quercus suber, aiming to comprehend the evolution of the IGS sequences within the family. Self- and cross-hybridization FISH was done on representative species of the Fagaceae. The IGS length variability and the methylation level of 18 and 25S rRNA genes were assessed in representatives of three genera of this family: Fagus, Quercus and Castanea. The intergenic spacers in Beech and Cork Oak showed similar overall organizations comprising putative functional elements needed for rRNA gene activity and containing a non-transcribed spacer (NTS), a promoter region, and a 5′-external transcribed spacer. In the NTS: the sub-repeats structure in Beech is more organized than in Cork Oak, sharing some short motifs which results in the lowest sequence similarity of the entire IGS; the AT-rich region differed in both spacers by a GC-rich block inserted in Cork Oak. The 5′-ETS is the region with the higher similarity, having nonetheless different lengths. FISH with the NTS-5′-ETS revealed fainter signals in cross-hybridization in agreement with the divergence between genera. The diversity of IGS lengths revealed variants from ∼2 kb in Fagus, and Quercus up to 5.3 kb in Castanea, and a lack of correlation between the number of variants and the number of rDNA loci in several species. Methylation of 25S Bam HI site was confirmed in all species and detected for the first time in the 18S of Q. suber and Q. faginea. These results provide important clues for the evolutionary trends of the rDNA 25S-18S IGS in the Fagaceae family.

The cytogenetic architecture of the aphid genome

In recent years aphids, with their well-defined polyphenism, have become favoured as model organisms for the study of epigenetic processes. The availability of the pea aphid (Acyrthosiphon pisum) genome sequence has engendered much research aimed at elucidating the mechanisms by which the phenotypic plasticity of aphids is inherited and controlled. Yet so far this research effort has paid little attention to the cytogenetic processes that play a vital part in the organisation, expression and inheritance of the aphid genome. Aphids have holocentric chromosomes, which have very different properties from the chromosomes with localised centromeres that are found in most other organisms. Here we review the diverse forms of aphid chromosome behaviour that occur during sex determination and male and female meiosis, often in response to environmental changes and mediated by endocrine factors. Remarkable differences occur, even between related species, that could have significant effects on the inheritance of all or parts of the genome. In relation to this, we review the particular features of the distribution of heterochromatin, rDNA genes and other repetitive DNA in aphid chromosomes, and discuss the part that these may play in the epigenetic modification of chromatin structure and function.

Animal plant warfare and secondary metabolite evolution

Abstract

The long-lasting discussion, why plants produce secondary metabolites, which are pharmacologically and toxicologically active towards mammals traces back to the eminent role of medicinal plants in the millennia-old history of manhood. In recent years, the concept of an animal plant warfare emerged, which focused on the co-evolution between plants and herbivores. As a reaction to herbivory, plants developed mechanical defenses such as thorns and hard shells, which paved the way for adapted animal physiques. Plants evolved further defense systems by producing chemicals that exert toxic effects on the animals that ingest them. As a result of this selective pressure, animals developed special enzymes, e.g. cytochrome P450 monooxigenases (CYP450) that metabolize xenobiotic phytochemicals. As a next step in the evolutionary competition between plants and animals, plants evolved to produce non-toxic pro-drugs, which become toxic only after ingestion by animals through metabolization by enzymes such as CYP450. Because these sequestered evolutionary developments call to mind an arms race, the term animal plant warfare has been coined. The evolutionary competition between plants and animals may help to better understand the modes of action of medicinal plants and to foster the efficient and safe use of phytotherapy nowadays.

Graphical abstract

Unusually high numbers of ribosomal RNA genes in copepods (Arthropoda: Crustacea) and their relationship to genome size

We report on copy numbers of 18S ribosomal RNA genes in three species of copepods (Crustacea: Copepoda), two of which possess an unusual arrangement in which 5S genes are included within the 18S-5.8S-28S repeat unit. Slot blots of genomic and standard DNA were hybridized with an 18S rRNA gene probe constructed from one of the marine species and hybridization was quantified using chemiluminescence. Diploid 18S rRNA gene copy numbers are estimated as ca. 15 300 and 33 500 in the marine species Calanus finmarchicus (13.0 pg DNA in 2C adult nuclei) and C. glacialis (24.2 pg DNA), respectively, and ca. 840 and 730 in two freshwater populations of Mesocyclops edax (both ca. 3 pg DNA) from Virginia and Nova Scotia, respectively. The roughly proportional relationship between 2C somatic nuclear DNA contents and rRNA gene copy number in the sibling species C. finmarchicus and C. glacialis may reflect polytenic replication of entire genomes during abrupt speciation events. Copy numbers may also reflect differential losses during embryonic chromatin diminution.

Organization and molecular evolution of CENP-A--associated satellite DNA families in a basal primate genome

Centromeric regions in many complex eukaryotic species contain highly repetitive satellite DNAs. Despite the diversity of centromeric DNA sequences among species, the functional centromeres in all species studied to date are marked by CENP-A, a centromere-specific histone H3 variant. Although it is well established that families of multimeric higher-order alpha satellite are conserved at the centromeres of human and great ape chromosomes and that diverged monomeric alpha satellite is found in old and new world monkey genomes, little is known about the organization, function, and evolution of centromeric sequences in more distant primates, including lemurs. Aye-Aye (Daubentonia madagascariensis) is a basal primate and is located at a key position in the evolutionary tree to study centromeric satellite transitions in primate genomes. Using the approach of chromatin immunoprecipitation with antibodies directed to CENP-A, we have identified two satellite families, Daubentonia madagascariensis Aye-Aye 1 (DMA1) and Daubentonia madagascariensis Aye-Aye 2 (DMA2), related to each other but unrelated in sequence to alpha satellite or any other previously described primate or mammalian satellite DNA families. Here, we describe the initial genomic and phylogenetic organization of DMA1 and DMA2 and present evidence of higher-order repeats in Aye-Aye centromeric domains, providing an opportunity to study the emergence of chromosome-specific modes of satellite DNA evolution in primate genomes.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

Background

Ribosomal 5S genes are well known for the critical role they play in ribosome folding and functionality. These genes are thought to evolve in a concerted fashion, with high rates of homogenization of gene copies. However, the majority of previous analyses regarding the evolutionary process of rDNA repeats were conducted in invertebrates and plants. Studies have also been conducted on vertebrates, but these analyses were usually restricted to the 18S, 5.8S and 28S rRNA genes. The recent identification of divergent 5S rRNA gene paralogs in the genomes of elasmobranches and teleost fishes indicate that the eukaryotic 5S rRNA gene family has a more complex genomic organization than previously thought. The availability of new sequence data from lower vertebrates such as teleosts and elasmobranches enables an enhanced evolutionary characterization of 5S rDNA among vertebrates.

Results

We identified two variant classes of 5S rDNA sequences in the genomes of Potamotrygonidae stingrays, similar to the genomes of other vertebrates. One class of 5S rRNA genes was shared only by elasmobranches. A broad comparative survey among 100 vertebrate species suggests that the 5S rRNA gene variants in fishes originated from rounds of genome duplication. These variants were then maintained or eliminated by birth-and-death mechanisms, under intense purifying selection. Clustered multiple copies of 5S rDNA variants could have arisen due to unequal crossing over mechanisms. Simultaneously, the distinct genome clusters were independently homogenized, resulting in the maintenance of clusters of highly similar repeats through concerted evolution.

Conclusions

We believe that 5S rDNA molecular evolution in fish genomes is driven by a mixed mechanism that integrates birth-and-death and concerted evolution.

Composition and epigenetic markers of heterochromatin in the aphid Aphis nerii (Hemiptera: Aphididae)

A detailed karyotype analysis of the oleander aphid Aphis nerii focusing on the distribution, molecular composition and epigenetic modifications of heterochromatin was done in order to better understand the structure and evolution of holocentric/holokinetic chromosomes in aphids. The female karyotype (2n = 8) consisted of 3 pairs of autosomes and a pair of X chromosomes that were the longest elements in the karyotype and carried a single, terminally located nucleolar organizer region. Males showed 2n = 7 chromosomes due to the presence of a single X chromosome. Heterochromatin was located in the X chromosomes only and consisted of 4 satellite DNAs that have been identified. A. nerii constitutive heterochromatin was enriched in mono-, di- and tri-methylated H3 histones and HP1 proteins but, interestingly, it lacked DNA methylation that was widespread in euchromatic chromosomal regions. These results suggest that aphid heterochromatin is assembled and condensed without any involvement of DNA methylation.

The library model for satellite DNA evolution: a case study with the rodents of the genus Ctenomys (Octodontidae) from the Iberá marsh, Argentina

On the basement of the library model of satellite DNA evolution is the differential amplification of subfamilies through lineages diversification. However, this idea has rarely been explored from an experimental point of view. In the present work, we analyzed copy number and sequence variability of RPCS (repetitive PvuII Ctenomys sequence), the major satellite DNA present in the genomes of the rodents of the genus Ctenomys, in a closely related group of species and forms inhabiting the Iberá marsh in Argentina. We studied the dependence of these two parameters at the intrapopulation level because in the case of interbreeding genomes, differences in RPCS copy number are due to recent amplification/contraction events. We found an inverse relationship among RPCS copy number and sequence variability: amplifications lead to a decrease in sequence variability, by means of biased homogenization of the overall satellite DNA, prevailing few variants. On the contrary, the contraction events that involve tandems of homogeneous monomers contribute-by default-minor variants to become "evident", which otherwise were undetectable. On the other hand, all the RPCS sequence variants are totally or partially shared by all the studied populations. As a whole, these results are comprehensible if these RPCS variants preexisted in the common ancestor of this Ctenomys group.

Elimination of 5S DNA unit classes in newly formed allopolyploids of the genera Aegilops and Triticum

Two classes of 5S DNA units, namely the short (containing units of 410 bp) and the long (containing units of 500 bp), are recognized in species of the wheat (the genera Aegilops and Triticum) group. While every diploid species of this group contains 2 unit classes, the short and the long, every allopolyploid species contains a smaller number of unit classes than the sum of the unit classes of its parental species. The aim of this study was to determine whether the reduction in these unit classes is due to the process of allopolyploidization, that is, interspecific or intergeneric hybridization followed by chromosome doubling, and whether it occurs during or soon after the formation of the allopolyploids. To study this, the number and types of unit classes were determined in several newly formed allotetraploids, allohexaploids, and an allooctoploid of Aegilops and Triticum. It was found that elimination of unit classes of 5S DNA occurred soon (in the first 3 generations) after the formation of the allopolyploids. This elimination was reproducible, that is, the same unit classes were eliminated in natural and synthetic allopolyploids having the same genomic combinations. No further elimination occurred in the unit classes of the 5S DNA during the life of the allopolyploid. The genetic and evolutionary significance of this elimination as well as the difference in response to allopolyploidization of 5S DNA and rDNA are discussed.

The evolution of the hexaploid grass Zingeriakochii (Mez) Tzvel. (2n=12) was accompanied by complex hybridization and uniparental loss of ribosomal DNA

In the grass tribe Poeae a small group of taxa occur with an exceptionally low chromosome number of 2n=2x=4 belonging to the closely related genera Colpodium and Zingeria. To understand the formation of polyploids in this group we analyzed the evolution of allohexaploid Zingeriakochii (2n=12) and its presumable ancestral species. Genomic insitu hybridization demonstrated that Z.kochii evolved from an interspecific hybrid involving species closely related to contemporary Z.biebersteiniana (2n=4) and Colpodiumversicolor (2n=4) and a third unknown species. Following allopolyploidization of Z.kochii the biebersteiniana-like parental chromosomes underwent loss of ribosomal DNA. No interlocus homogenization of 45S rDNA took place in Z.kochii and phylogenetic analysis showed that C.versicolor contributed its genome to Z.kochii relatively recently. Insitu hybridization was particularly effective in understanding the allopolyploid evolution in Zingeria while the analysis of ITS sequences alone would have resulted in a wrong interpretation of the allopolyploid history of the genus.

Low rates of homogenization of the DBC-150 satellite DNA family restricted to a single pair of microchromosomes in species from the Drosophila buzzatii cluster

A satellite DNA family, termed DBC-150, comprises slightly GC-rich repeat units of approximately 150 bp that were isolated (by DNA digestions or PCR) from the genome of all seven Drosophila species from the buzzatii cluster (repleta group). The presence of subrepeats suggests that part of the extant DBC-150 monomer originated by the duplication of small sequence motifs. The DBC-150 family is compared to the previously described pBuM satDNA family, an abundant component of the genome of five species of the cluster. The two families are different in several aspects, including primary structure, A + T content, intraspecific and interspecific variability and rates of homogenization (or nucleotide spread). The data indicate a lower rate of homogenization (and absence of complete concerted evolution) of the DBC-150 compared to the pBuM family. FISH on metaphase chromosomes revealed that the DBC-150 family is located exclusively in the microchromosomes. To our knowledge this is the first record of a complex Drosophila satDNA restricted to a single pair of microchromosomes. The observed low rates of homogenization of the DBC-150 family might be related to a presumed reduction or suppression of meiotic recombination in the microchromosomes.

Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.

Wax plants disentangled: a phylogeny of Hoya (Marsdenieae, Apocynaceae) inferred from nuclear and chloroplast DNA sequences

Hoya (Marsdenieae, Apocynaceae) includes at least 200 species distributed from India to the Pacific Islands. We here infer major species groups in the genus based on combined sequences from the chloroplast atpB-rbcL spacer, the trnL region, and nuclear ribosomal DNA ITS region for 42 taxa of Hoya and close relatives. To assess levels of ITS polymorphism, ITS sequences for a third of the accessions were obtained by cloning. Most ITS clones grouped by species, indicating that speciation in Hoya usually predates ITS duplication. One ITS sequence of H. carnosa, however, grouped with a sequence of the morphologically similar H. pubicalyx, pointing to recent hybridization or the persistence of paralogous copies through a speciation event. The topology resulting from the combined chloroplast and nuclear data recovers some morphology-based sections, such as Acanthostemma and Eriostemma, as well as a well-supported Australian/New Guinean clade. The combined data also suggest that morphological adaptations for ant-symbiosis evolved at least three times within Hoya.

Satellite DNA from the Y chromosome of the malaria vector Anopheles gambiae

Satellite DNA is an enigmatic component of genomic DNA with unclear function that has been regarded as "junk." Yet, persistence of these tandem highly repetitive sequences in heterochromatic regions of most eukaryotic chromosomes attests to their importance in the genome. We explored the Anopheles gambiae genome for the presence of satellite repeats and identified 12 novel satellite DNA families. Certain families were found in close juxtaposition within the genome. Six satellites, falling into two evolutionarily linked groups, were investigated in detail. Four of them were experimentally confirmed to be linked to the Y chromosome, whereas their relatives occupy centromeric regions of either the X chromosome or the autosomes. A complex evolutionary pattern was revealed among the AgY477-like satellites, suggesting their rapid turnover in the A. gambiae complex and, potentially, recombination between sex chromosomes. The substitution pattern suggested rolling circle replication as an array expansion mechanism in the Y-linked 53-bp satellite families. Despite residing in different portions of the genome, the 53-bp satellites share the same monomer lengths, apparently maintained by molecular drive or structural constraints. Potential functional centromeric DNA structures, consisting of twofold dyad symmetries flanked by a common sequence motif, have been identified in both satellite groups.

Molecular evolution of rDNA external transcribed spacer and phylogeny of sect. Petota (genus Solanum)

The 5(') external transcribed spacer (ETS) region of ribosomal DNA of 30 species of Solanum sect. Petota and the European Solanum dulcamara were compared. Two structural elements can be distinguished in the ETS: (i). a variable region (VR), demonstrating significant structural rearrangements and (ii). a conservative region (CR), evolving mainly by base substitutions. In VR, a conservative element (CE) with similarity to the ETS of distantly related Nicotiana is present. The ancestral organization of ETS (variant A) was found for non-tuber-bearing species of ser. Etuberosa, tuber-bearing wild potatoes of Central American ser. Bulbocastana, Pinnatisecta, and Polyadenia and S. dulcamara. Duplication of CE took place in the ETS of species from ser. Commersoniana and Circaeifolia (variant B). South American diploids and Mexican polyploids from superser. Rotata also possess two CE, and additionally two duplications around CE1 are present in VR (variant C). Three major lineages could be distinguished: non-tuber-bearing species of ser. Etuberosa, tuber-bearing Central American diploids and all South American species radiated from a common ancestor at early stages of evolution, indicating a South American origin of the tuber-bearing species. Later, Central and South American diploids evolved further as independent lineages. South American species form a monophyletic group composed of series with both stellata and rotata flower morphology. Solanum commersonii represents a sister taxon for all rotata species, whereas ser. Circaeifolia diverged earlier. Two main groups, C1 and C2, may be distinguished for species possessing ETS variant C. C1 contains ser. Megistacroloba, Conicibaccata, Maglia, and Acaulia, whereas all diploids of ser. Tuberosa are combined into C2. A closer relationship of Solanum chacoense (ser. Yungasensa) to the C2 group was found. The origin of polyploid species Solanum maglia, Solanum acaule, Solanum tuberosum, Solanum iopetalum, and Solanum demissum is discussed.

Regulatory motifs are present in the ITS1 of some flatworm species

Particular sequence motifs can act as transcription regulators. Because the total regulatory effects of such motifs can be related to their abundance, their presence might be expected at locations within the genome where sequences are repeated. Multiple repeats that vary in number among individuals occur within the ribosomal first internal transcribed spacer (ITS1) in some species in three trematode genera: Paragonimus, Schistosoma and Dolichosaccus. In all of these genera we found in ITS1, sequences identical to known enhancer motifs. We also searched for, and identified, known regulatory motifs in published ITS1 sequences of other parasitic flatworms including Echinostoma spp. (Trematoda) and Echinococcus spp. (Cestoda) which lack multiple repeats in ITS1. We present three lines of evidence that this widespread occurrence of such motifs within the ITS1 of parasitic flatworms may indicate a functional role in regulating tissue- or stage-specific transcription of ribosomal genes. First, these motifs are identical to ones whose functional roles have been established using in vitro assays of transcriptional rates. Second, in all 18 species investigated here, between one and three different regulatory motifs were identified. In 14 of these 18 species, the probability that at least one of these motifs occurred because of the random assortment of bases within the regions investigated was 10% or less. In 12 of these 14 species, the probability was 5% or less. Third, the evolutionary divergence of flatworm species investigated is quite ancient. Therefore, the interspecific distribution of motifs observed here, in a rapidly evolving region such as ITS1, is unlikely to be attributable solely to shared evolutionary histories. These results, therefore, suggest a broader functional role for the ITS1 than previously thought.

Uniparental loss of ribosomal DNA in the allotetraploid grass Zingeria trichopoda (2n = 8)

Analysis of the grass Zingeria trichopoda (2n = 8, 2C = 5.3 pg) revealed a dynamic evolution with the following characteristics. (i) Genomic in situ hybridization (GISH) demonstrates that Z. trichopoda evolved from an interspecific hybrid involving a species like contemporary Zingeria biebersteiniana (2n = 4) and a second species with a similar low number of chromosomes. The nucleus of Z. trichopoda is spatially organized at the genome level and the two parental genomes occupy distinct and separate domains of lateral arrangements. (ii) The copy number of the Z. biebersteiniana specific pericentromeric tandem repeat family Zbcen1 is drastically reduced in Z. trichopoda. (iii) GISH in combination with labeled rDNA sequences simultaneously discriminated the two parental genomes and the corresponding 5S and 45S rDNA sites. Hence, following allopolyploidization of Z. trichopoda the Z. biebersteiniana like parental chromosomes probably underwent drastic loss of 45S rDNA. This could have arisen either through the loss of Z. biebersteiniana derived 45S rDNA or through Z. trichopoda genome-wide homogenization of Z. biebersteiniana type 45S rDNA and subsequent elimination of 45S rDNA loci from Z. biebersteiniana derived chromosomes. Finally, 5S rDNA loci are present in both subgenomes of Z. trichopoda and the chromosomal position of these loci is similar for both Z. biebersteiniana and the Z. biebersteiniana like parental genome of Z. trichopoda.

Highly repetitive DNA families restricted to germ cells in a Japanese hagfish (Eptatretus burgeri): a hierarchical and mosaic structure in eliminated chromosomes

It is known that in eight hagfishes chromosome elimination occurs during early embryogenesis. The eliminated chromosomes are mostly C-band positive, so that none of the somatic cells have any C-band-positive chromatin. Recently, some highly repetitive DNA sequences have been reported as eliminated elements in these hagfishes based on molecular biological methods. However, no germline-restricted repetitive DNA have been directly isolated from the Japanese hagfish Eptatretus burgeri, from which approximately 21% of the total DNA is eliminated from presumptive somatic cells. Through electrophoretic investigation after digestion with restriction endonucleases, two DNA families that are restricted to germline DNA were isolated. Molecular cloning and sequence analysis revealed that these families are composed of closely related sequences of 64 and 57bp in length, respectively. Southern blot hybridization revealed that the two DNA families are restricted to germline DNA and were thus named EEEb1 and EEEb2, respectively. Moreover, these eliminated elements were highly and tandemly repeated, and it is predicted that they might amplify by saltatory replication and have evolved in a concerted manner. By densitometric scanning, EEEb1 and EEEb2 were found to amount to make up approximately 18.5 and 0.024% of the total germline genomic DNA, accounting for 88.6% of the total eliminated DNA. A fluorescence in situ hybridization experiment demonstrated that EEEb1 is located on all C-band-positive chromosomes that are limited to germ cells, suggesting that EEEb1 is the primary component of eliminated DNA of E. burgeri.

The phylogeography and connectivity of the latitudinally widespread scleractinian coral Plesiastrea versipora in the Western Pacific

Whereas terrestrial animal populations might show genetic connectivity within a continent, marine species, such as hermatypic corals, may have connectivity stretching to all corners of the planet. We quantified the genetic variability within and among populations of the widespread scleractinian coral, Plesiastrea versipora along the eastern Australian seaboard (4145 km) and the Ryukyu Archipelago (Japan, 681 km) using sequences of internal transcribed spacers (ITS1-2) from ribosomal DNA. Geographic patterns in genetic variability were deduced from a nested clade analysis (NCA) performed on a parsimony network haplotype. This analysis allowed the establishment of geographical associations in the distribution of haplotypes within the network cladogram, therefore allowing us to deduce phylogeographical patterns based under models of restricted gene flow, fragmentation and range expansion. No significant structure was found among Ryukyu Archipelago populations. The lack of an association between the positions of haplotypes in the cladogram with geographical location of these populations may be accounted for by a high level of gene flow of P. versipora within this region, probably due to the strong Kuroshio Current. In contrast, strong geographical associations were apparent among populations of P. versipora along the south-east coast of Australia. This pattern of restricted genetic connectivity among populations of P. versipora on the eastern seaboard of Australia seems to be associated with the present surface ocean current (the East Australian Current) on this side of the south-western Pacific Ocean.

Cytochromes P450 of insects: the tip of the iceberg

The cytochrome P450-dependent monooxygenases are an extremely important metabolic system involved in the metabolism of endogenous compounds and xenobiotics. Collectively, P450 monooxygenases can metabolize numerous substrates and carry out multiple oxidative reactions. The large number of substrates metabolized is due to the plethora of P450 isoforms and to the broad substrate specificity of some isoforms. Monooxygenases of insects have several functional roles, including growth, development, feeding and protection against xenobiotics, including resistance to pesticides and tolerance to plant toxins. This review begins with background information about P450s and their evolution, followed by a discussion of the extraordinary diversity of insect P450s. Given the enormous interest in studying individual P450s, we then provide a synopsis of the different methods that have been used in their isolation and the substrates that are known to be metabolized. We conclude by summarizing the lessons we have learned from the study of individual insect P450s, including their roles in insecticide resistance, plant-insect interactions and insect physiology. However, these studies are just the 'tip of the iceberg'. Our knowledge continues to expand at a rapid pace, suggesting that the next decade will outpace the last in terms of improving our understanding of the cytochromes P450 of insects.

The evolution of controlled multitasked gene networks: the role of introns and other noncoding RNAs in the development of complex organisms

Eukaryotic phenotypic diversity arises from multitasking of a core proteome of limited size. Multitasking is routine in computers, as well as in other sophisticated information systems, and requires multiple inputs and outputs to control and integrate network activity. Higher eukaryotes have a mosaic gene structure with a dual output, mRNA (protein-coding) sequences and introns, which are released from the pre-mRNA by posttranscriptional processing. Introns have been enormously successful as a class of sequences and comprise up to 95% of the primary transcripts of protein-coding genes in mammals. In addition, many other transcripts (perhaps more than half) do not encode proteins at all, but appear both to be developmentally regulated and to have genetic function. We suggest that these RNAs (eRNAs) have evolved to function as endogenous network control molecules which enable direct gene-gene communication and multitasking of eukaryotic genomes. Analysis of a range of complex genetic phenomena in which RNA is involved or implicated, including co-suppression, transgene silencing, RNA interference, imprinting, methylation, and transvection, suggests that a higher-order regulatory system based on RNA signals operates in the higher eukaryotes and involves chromatin remodeling as well as other RNA-DNA, RNA-RNA, and RNA-protein interactions. The evolution of densely connected gene networks would be expected to result in a relatively stable core proteome due to the multiple reuse of components, implying that cellular differentiation and phenotypic variation in the higher eukaryotes results primarily from variation in the control architecture. Thus, network integration and multitasking using trans-acting RNA molecules produced in parallel with protein-coding sequences may underpin both the evolution of developmentally sophisticated multicellular organisms and the rapid expansion of phenotypic complexity into uncontested environments such as those initiated in the Cambrian radiation and those seen after major extinction events.

CYP9E2, CYP4C21 and related pseudogenes from German cockroaches, Blattella germanica: implications for molecular evolution, expression studies and nomenclature of P450s

The cDNAs of two novel P450s (CYP9E2 and CYP4C21) were isolated from German cockroaches, Blattella germanica. Both CYP9E2 and CYP4C21 are typical microsomal P450s and their deduced amino acid sequences share a number of common characteristics with other members of the P450 superfamily. Northern blot analyses using a CYP9E2 or CYP4C21 probe showed that 'CYP9E2' and 'CYP4C21' were expressed at all life stages. Two pseudogenes related to CYP9E2 and three pseudogenes related to CYP4C21 were also isolated. These represent the first P450 pseudogenes from an insect other than Drosophila melanogaster. The relative number of P450 pseudogenes in B. germanica is apparently higher than in D. melanogaster. The implications of these results for the molecular evolution, expression studies and nomenclature of P450s are discussed.

Diversity of the immunoglobulin heavy chain in the Atlantic salmon (Salmo salar L.) is contributed by genes from two parallel IgH isoloci

Immunoglobulin heavy chain (IgH) variable (V) region cDNAs from the Atlantic salmon, Salmo salar L., have been isolated and analysed with respect to diversity and transcription of the two parallel IgH isoloci in this species. A total of nine V(H) families were defined according to the 80% identity criterion, of which seven were highly related (>80% identity) to the V(H) families defined in rainbow trout and arctic charr. The variability of the CDR1 and 2 was low, although mutational hot-spot consensus sequences were accumulated in these regions. The CDR3 showed largest variability, expressing at least eight different groups of D motifs diversified by fusion of the D motifs, possible N and P nucleotide insertions and exonuclease activity. Presumably functional transcripts expressing D motifs in all three reading frames were identified for two of the motifs. The cDNAs were mapped to either of the two parallel loci, and sequence analysis revealed that the repertoire of V(H) segments was contributed by transcription of genes from both of the IgH isoloci. Transcription of genes from both isoloci generated no obvious effects on variability in the CDR3 of the Atlantic salmon IgH chains, although one additional J(H)-segment with altered N-terminal was generated by the process of duplication and divergence. Thus, the issue of biological significance of the two IgH isoloci remains unclear.

Evolutionary dynamics of the SGM transposon family in the Drosophila obscura species group

SGM (Drosophila subobscura, Drosophila guanche, and Drosophila madeirensis) transposons are a family of transposable elements (TEs) in Drosophila with some functional and structural similarities to miniature inverted-repeat transposable elements (MITEs). These elements were recently active in D. subobscura and D. madeirensis (1-2 MYA), but in D. guanche (3-4 MYA), they gave rise to a species-specifically amplified satellite DNA making up approximately 10% of its genome. SGM elements were already active in the common ancestor of all three species, giving rise to the A-type specific promoter section of the P:-related neogene cluster. SGM sequences are similar to elements found in other obscura group species, such as the ISY elements in D. miranda and the ISamb elements in Drosophila ambigua. SGM elements are composed of different sequence modules, and some of them, i.e., LS and LS-core, are found throughout the Drosophila and Sophophora radiation with similarity to more distantly related TEs. The LS-core module is highly enriched in the noncoding sections of the Drosophila melanogaster genome, suggesting potential regulatory host gene functions. The SGM elements can be considered as a model system elucidating the evolutionary dynamics of mobile elements in their arms race with host-directed silencing mechanisms and their evolutionary impact on the structure and composition of their respective host genomes.

Molecular evolution of the opioid/orphanin gene family

Gene duplication is a recurring theme in the evolution of vertebrate polypeptide hormones and neuropeptides. These duplication events can lead to the formation of gene families in which divergence of function is the usual outcome. In the case of the opioid/orphanin family of genes, duplication events have proceeded along two paths: (a) an apparent duplication of function as seen in the analgesic activity of Proenkephalin and Prodynorphin end-products; and (b) divergence of function as seen in the nociceptic activity of Proorphanin end-products or the melanocortin (color change and chronic stress regulation) activity of Proopiomelanocortin end-products. Although genes coding for Proopiomelanocortin, Proenkephalin, Prodynorphin, and Proorphanin have been extensively studied in mammals, the distribution and radiation of these genes in nonmammalian vertebrates is less well understood. This review will present the hypothesis that the radiation of the opioid/orphanin gene family is the result of the duplication and divergence of the Proenkephalin gene during the radiation of the chordates. To evaluate the Proenkephalin gene duplication hypothesis, a 3'RACE procedure was used to screen for the presence of Prodynorphin-related, Proenkephalin-related, and Proorphanin-related cDNAs expressed in the brains of nonmammalian vertebrates.

Convergent evolution of antifreeze glycoproteins in Antarctic notothenioid fish and Arctic cod

Antarctic notothenioid fishes and several northern cods are phylogenetically distant (in different orders and superorders), yet produce near-identical antifreeze glycoproteins (AFGPs) to survive in their respective freezing environments. AFGPs in both fishes are made as a family of discretely sized polymers composed of a simple glycotripeptide monomeric repeat. Characterizations of the AFGP genes from notothenioids and the Arctic cod show that their AFGPs are both encoded by a family of polyprotein genes, with each gene encoding multiple AFGP molecules linked in tandem by small cleavable spacers. Despite these apparent similarities, detailed analyses of the AFGP gene sequences and substructures provide strong evidence that AFGPs in these two polar fishes in fact evolved independently. First, although Antarctic notothenioid AFGP genes have been shown to originate from a pancreatic trypsinogen, Arctic cod AFGP genes share no sequence identity with the trypsinogen gene, indicating trypsinogen is not the progenitor. Second, the AFGP genes of the two fish have different intron-exon organizations and different spacer sequences and, thus, different processing of the polyprotein precursors, consistent with separate genomic origins. Third, the repetitive AFGP tripeptide (Thr-Ala/Pro-Ala) coding sequences are drastically different in the two groups of genes, suggesting that they arose from duplications of two distinct, short ancestral sequences with a different permutation of three codons for the same tripeptide. The molecular evidence for separate ancestry is supported by morphological, paleontological, and paleoclimatic evidence, which collectively indicate that these two polar fishes evolved their respective AFGPs separately and thus arrived at the same AFGPs through convergent evolution.

Evolution of antifreeze glycoprotein gene from a trypsinogen gene in Antarctic notothenioid fish

Freezing avoidance conferred by different types of antifreeze proteins in various polar and subpolar fishes represents a remarkable example of cold adaptation, but how these unique proteins arose is unknown. We have found that the antifreeze glycoproteins (AFGPs) of the predominant Antarctic fish taxon, the notothenioids, evolved from a pancreatic trypsinogen. We have determined the likely evolutionary process by which this occurred through characterization and analyses of notothenioid AFGP and trypsinogen genes. The primordial AFGP gene apparently arose through recruitment of the 5' and 3' ends of an ancestral trypsinogen gene, which provided the secretory signal and the 3' untranslated region, respectively, plus de novo amplification of a 9-nt Thr-Ala-Ala coding element from the trypsinogen progenitor to create a new protein coding region for the repetitive tripeptide backbone of the antifreeze protein. The small sequence divergence (4-7%) between notothenioid AFGP and trypsinogen genes indicates that the transformation of the proteinase gene into the novel ice-binding protein gene occurred quite recently, about 5-14 million years ago (mya), which is highly consistent with the estimated times of the freezing of the Antarctic Ocean at 10-14 mya, and of the main phyletic divergence of the AFGP-bearing notothenioid families at 7-15 mya. The notothenioid trypsinogen to AFGP conversion is the first clear example of how an old protein gene spawned a new gene for an entirely new protein with a new function. It also represents a rare instance in which protein evolution, organismal adaptation, and environmental conditions can be linked directly.