Evolution of the trnF(GAA) gene in Arabidopsis relatives and the brassicaceae family: monophyletic origin and subsequent diversification of a plastidic pseudogene

Recently, we used the 5'-trnL(UAA)-trnF(GAA) region of the chloroplast DNA for phylogeographic reconstructions and phylogenetic analysis among the genera Arabidopsis, Boechera, Rorippa, Nasturtium, and Cardamine. Despite the fact that extensive gene duplications are rare among the chloroplast genome of higher plants, within these taxa the anticodon domain of the trnF(GAA) gene exhibit extensive gene duplications with one to eight tandemly repeated copies in close 5' proximity of the functional gene. Interestingly, even in Arabidopsis thaliana we found six putative pseudogenic copies of the functional trnF gene within the 5'-intergenic trnL-trnF spacer. A reexamination of trnL(UAA)-trnF(GAA) regions from numerous published phylogenetic studies among halimolobine, cardaminoid, and other cruciferous taxa revealed not only extensive trnF gene duplications but also favor the hypothesis about a single origin of trnF pseudogene formation during evolution of the Brassicaceae family 16-21 MYA. Conserved sequence motifs from this tandemly repeated region are codistributed nonrandomly throughout the plastome, and we found some similarities with a DNA sequence duplication in the rps7 gene and its adjacent spacer. Our results demonstrate the potential evolutionary dynamics of a plastidic region generally regarded as highly conserved and probably cotranscribed and, as shown here for several genera among cruciferous plants, greatly characterized by parallel gains and losses of duplicated trnF copies.     

Recombination diversifies chloroplast trnF pseudogenes in Arabidopsis lyrata

Extensive intraspecific variation in the chloroplast trnL(UAA)-trnF(GAA) spacer of model plant Arabidopsis lyrata is caused by multiple copies of a tandemly repeated trnF pseudogene undergoing parallel independent changes in copy number. Linkage disequilibrium and secondary structure analyses indicate that the diversification of pseudogene copies is driven by complex processes of structurally mediated illegitimate recombination. Disperse repeats sharing similar secondary structures interact, facilitating reciprocal exchange of structural motifs between copies via intramolecular and intermolecular recombinations, forming chimeric sequences and iterative expansion and contraction in pseudogene copy numbers. Widely held assumptions that chloroplast sequence evolution is simple and structural changes are informative are violated. Our findings have important implications for the use of this highly variable region in Brassicaceae studies. The reticulate evolution and nonindependent nucleotide substitution render the pseudogene inappropriate for standard phylogenetic reconstruction, but over short evolutionary timescales they may be useful for assessing gene flow, hybridization and introgression.     

Supernetwork identifies multiple events of plastid trnF(GAA) pseudogene evolution in the Brassicaceae

The occurrence of nonfunctional trnF pseudogenes has been rarely described in flowering plants. However, we describe the first large-scale supernetwork for the Brassiccaeae built from gene trees for 5 loci (adh, chs, matK, trnL-F, and ITS) and report multiple independent origins for trnF pseudogenes in crucifers. The duplicated regions of the original trnF gene are comprised of its anticodon domain and several other highly structured motifs not related to the original gene. Length variation of the trnL-F intergenic spacer region in different taxa ranges from 219 to 900 bp as a result of differences in pseudocopy number (1-14). It is speculated that functional constraints favor 2-3 or 5-6 copies, as found in Arabidopsis and Boechera. The phylogenetic distribution of microstructural changes for the trnL-F region supports ancient patterns of divergence in crucifer evolution for some but not all gene loci.     

Plastidic trnFUUC pseudogenes in North American genus Boechera (Brassicaceae): mechanistic aspects of evolution

The origin and maintenance of a plastidic tandem repeat next to the TRNF (UUC) gene were analyzed in the genus BOECHERA in a phylogenetic context and were compared to published analogous examples that emerged in parallel in the Asteraceae and Juncaceae, respectively. Although we identified some features common to these taxonomic groups with respect to structure and origin of the region, obvious differences were encountered, which argue against a specific mechanism or evolutionary principle underlying the parallel origin and maintenance of the TRNF-tandem repeats in those families. In contrast to the situation in the Asteraceae, no reciprocal recombinant repeat types have been observed in the Brassicaceae. Forty copy types, classified into three groups, were isolated from 103 chloroplast haplotypes of BOECHERA and it was demonstrated that they are composed of four subregions of various origins. We discuss various mutation mechanisms such as DNA replication slippage, and inter- and intrachromosomal recombination which were reported to mediate variation in copy numbers and other types of observed sequence length polymorphism. It is shown that the observed molecular structure of the tandem repeat region did not fully fit the particular patterns expected under a scenario of evolution including any of the known mechanisms. Nevertheless, it appeared that intermolecular unequal crossing-over is most likely the driving force in the evolution of this tandem repeat. However, it remains to be explained, why no reciprocal recombinant copy types have been observed. The reconstructed phylogenetic relationships among copies reflected different evolutionary scenarios as follows: (1) A single and ancient origin of copies pre-dates the radiation of BOECHERA. (2) Parallel expansion and shortening of the tandem repeat within different BOECHERA lineages. (3) Conservation of the first copy, as it was the only one present in all chloroplast haplotypes.     

Evolution of trnF(GAA) pseudogenes in cruciferous plants

In several studies we used the 5′-trnL(UAA)–trnF(GAA) region of the chloroplast DNA for phylogeographic reconstructions, gene diversity calculations and phylogenetic analyses among the genera Arabidopsis and Boechera. Despite the fact that extensive gene duplications are rare within the chloroplast genome of higher plants, within several genera of the Brassicaceae the anticodon domain of the trnF(GAA) gene exhibit extensive gene duplications with 1–12 tandemly repeated copies in close 5′-proximity of the functional gene. A recent re-examination and additional analysis of trnL(UAA)–trnF(GAA) regions from numerous cruciferous taxa not only reveal extensive trnF gene duplications, but also favour the hypothesis that in cruciferous taxa at least four independent phylogenetic lineages are characterized by these pseudogenes. Among these lineages there is one major clade of taxa carrying pseudogenes indicating an ancient split in crucifer evolution. In two case studies, Boechera and Arabidopsis, intra- and inter-molecular recombinations have been shown to be the reason for the reciprocal exchange of several similar motifs. However, functional constraints might favour two to three or five to six copies as shown for Arabidopsis and Boechera. Herein, we compare the occurrence and distribution of pseudogene copy number in the framework of a comprehensive survey of cpDNA haplotype variation in Boechera, the former genus Cardaminopsis and Arabidopsis thaliana and comment on the value of such kind of mutations in phylogenetic and evolutionary reconstructions.     

Low variability of internal transcribed spacer rDNA and trnL (UAA) intron sequences of several taxa in the Festuca ovina aggregate (POACEAE)

Identification and classification of numerous Festuca species is still a difficult problem due to the close morphological resemblance. The most difficult fine fescues to identify belong to the Festuca ovina aggregate, which is the largest group in the genus Festuca. Many taxons are considered to be separate species based on quantitative taxonomic characters, differences in ploidy level or the structure of sclerenchyma cells. In order to evaluate the taxonomic value of DNA-based markers, sequence analysis of the internal transcribed spacer (ITS1-5.8S-ITS2) region and the chloroplast trnL (UAA) intron was performed in the ten most problematic fine fescues belonging to the Festuca ovina aggregate. Intraspecific ITS variants were found in a single case while in other cases only intragenomic ITS polymorphisms were detected with 1-2 ambiguous positions. Among the sequences of the trnL (UAA) intron even intragenomic polymorphisms were not detected in any of the Festuca species studied. Thus, the results do not support the species status of these ten taxa.     

Transformation of Chloroplast Ribosomal RNA Genes in Chlamydomonas: Molecular and Genetic Characterization of Integration Events

Transformation of chloroplast ribosomal RNA (rRNA) genes in Chlamydomonas has been achieved by the biolistic process using cloned chloroplast DNA fragments carrying mutations that confer antibiotic resistance. The sites of exchange employed during the integration of the donor DNA into the recipient genome have been localized using a combination of antibiotic resistance mutations in the 16S and 23S rRNA genes and restriction fragment length polymorphisms that flank these genes. Complete or nearly complete replacement of a region of the chloroplast genome in the recipient cell by the corresponding sequence from the donor plasmid was the most common integration event. Exchange events between the homologous donor and recipient sequences occurred preferentially near the vector:insert junctions. Insertion of the donor rRNA genes and flanking sequences into one inverted repeat of the recipient genome was followed by intramolecular copy correction so that both copies of the inverted repeat acquired identical sequences. Increased frequencies of rRNA gene transformants were achieved by reducing the copy number of the chloroplast genome in the recipient cells and by decreasing the heterology between donor and recipient DNA sequences flanking the selectable markers. In addition to producing bona fide chloroplast rRNA transformants, the biolistic process induced mutants resistant to low levels of streptomycin, typical of nuclear mutations in Chlamydomonas.     

Phylogenetic Relationships Among Olea Species, Based on Nucleotide Variation at a Non-Coding Chloroplast DNA Region

Abstract: The purpose of this study was to assess nucleotide variation at a non-coding chloroplast DNA region in Olea species, to evaluate their phylogenetic relationships within the Olea genus and, more particularly, to clarify the relationships between cultivated olive (O. europaea) and the other taxa of section Olea. The analysis was made on an intergenic region between the trnT (UGU) and trnL (UAA) 5' exon, within a large single copy region of the chloroplast genome. Site-specific primers were used to amplify the region by PCR. This sequence analysis was applied to the same array of Olea species as assayed by Lumaret et al. (2000^[16]) using cpDNA RFLPs, thus making it possible to compare phylogenetic relationships analysed at two complementary levels of cpDNA variation. On the 666 bp aligned sequence, 8 different haplotypes were defined, with 9 single nucleotide mutations, a different length of a poly-T region and an indel for O. paniculata. Haplotypes were shared by the species pairs O. europaea-O.laperrinei, O. maroccana-O. cerasiformis, O. capensis-O. lancea and O. africana-O.indica. Phylogenetic analyses of these data distinguished four groups: the species Olea capensis and O. lancea, which both belong to subgenus Ligustroides, the Olea forms from southeast Africa, those from Asia and the taxa of northwest Africa and the Mediterranean Basin, which include olive crop. The results are consistent with those previously found using cpDNA RFLPs, with some minor differences observed within each group. They constitute further evidence to clarify the phylogeny of Olea.     

Length Variation in Mitochondrial DNA of the Minnow Cyprinella Spiloptera

Length differences in animal mitochondrial DNA (mtDNA) are common, frequently due to variation in copy number of direct tandem duplications. While such duplications appear to form without great difficulty in some taxonomic groups, they appear to be relatively short-lived, as typical duplication products are geographically restricted within species and infrequently shared among species. To better understand such length variation, we have studied a tandem and direct duplication of approximately 260 bp in the control region of the cyprinid fish, Cyprinella spiloptera. Restriction site analysis of 38 individuals was used to characterize population structure and the distribution of variation in repeat copy number. This revealed two length variants, including individuals with two or three copies of the repeat, and little geographic structure among populations. No standard length (single copy) genomes were found and heteroplasmy, a common feature of length variation in other taxa, was absent. Nucleotide sequence of tandem duplications and flanking regions localized duplication junctions in the phenylalanine tRNA and near the origin of replication. The locations of these junctions and the stability of folded repeat copies support the hypothesized importance of secondary structures in models of duplication formation.     

Excessive variation in Y chromosomal DNA in Rumex acetosa (Polygonaceae)

Rumex acetosa is one of the few angiosperms that possesses sex chromosomes. The same types of abundant repetitive sequences cover both heterochromatic Y chromosomes present in males. The aim of this study was to investigate genetic variation in paternally inherited Y chromosomal DNA and in maternally inherited cpDNA, and to find out whether the examined genomic regions are suited to a phylogeographic study in R. acetosa. DNA sequence polymorphisms present in the 850-bp heterochromatic segment on the Y chromosomes were compared to variation in the 409-bp long chloroplast section (trnL- trnF spacer) in R. acetosa originating from several European locations and from the Altai mountains in Russia. A great amount of genetic variation was detected within the Y chromosomal region while only four chloroplast genotypes were detected. Although the chloroplast haplotypes possessed some geographic pattern, no clear phylogeographic pattern was detected based on the variable Y chromosomes. The mean Y chromosomal nucleotide diversity among all samples equaled 6.6 %, and the mean proportion of polymorphic sites per individual equaled 8.2 % among SNP sites and 1.7 % among all sites investigated. The high number of substitutions detected in the Y chromosomal DNA shows that this heterochromatic sequence has a high mutation rate. The diversity pattern indicates that gene flow via pollen is extensive and it blurs any geographical pattern in the Y chromosomal variation. The high number of repeats and uncertainty concerning the extent of recombination between the two Y chromosomes impair the usability of the Y chromosomal segment for phylogeographic or population genetic studies.     

Length Variation, Heteroplasmy and Sequence Divergence in the Mitochondrial DNA of Four Species of Sturgeon (Acipenser)

The extent of mtDNA length variation and heteroplasmy as well as DNA sequences of the control region and two tRNA genes were determined for four North American sturgeon species: Acipenser transmontanus, A. medirostris, A. fulvescens and A. oxyrhnychus. Across the Continental Divide, a division in the occurrence of length variation and heteroplasmy was observed that was concordant with species biogeography as well as with phylogenies inferred from restriction fragment length polymorphisms (RFLP) of whole mtDNA and pairwise comparisons of unique sequences of the control region. In all species, mtDNA length variation was due to repeated arrays of 78-82-bp sequences each containing a D-loop strand synthesis termination associated sequence (TAS). Individual repeats showed greater sequence conservation within individuals and species rather than between species, which is suggestive of concerted evolution. Differences in the frequencies of multiple copy genomes and heteroplasmy among the four species may be ascribed to differences in the rates of recurrent mutation. A mechanism that may offset the high rate of mutation for increased copy number is suggested on the basis that an increase in the number of functional TAS motifs might reduce the frequency of successfully initiated H-strand replications.     

Deep genomic analysis of Coelastrella saipanensis (Scenedesmaceae,Chlorophyta): comparative chloroplast genomics of Scenedesmaceae

Many species belonging to the coccoid green algae genus Coelastrella are considered potential candidates for the large-scale production of natural pigments and biofuels. However, little is known about the structural, functional and molecular aspects of the chloroplast genomes (cpDNAs) of this genus. In the present study, the complete sequence of the cpDNA of strain FACHB-2138, which was further identified as Coelastrella saipanensis Hanagata based on morphological and molecular analyses, was elucidated. The 196 140 bp cpDNA sequence that was assembled as a circular map was found to possess the typical quadripartite structure. The two identical copies of 11 897 bp inverted repeat (IR) sequences were separated from one another by single copy regions. The large single copy region (LSC) was 104 949 bp, whereas the small single copy region (SSC) was 67 397 bp. The cpDNA encoded a total of 96 unique genes, which included 67 protein-coding genes, three rRNA genes and 26 tRNA genes. A total of 19 group I introns were annotated in this genome. Comparative analyses with three species from the family Scenedesmaceae showed C. saipanensis had a slightly expanded genome, higher GC content and less skewed distribution of its genes between the two DNA strands than that of the other three species. The cpDNA data deduced from the present study helps to expand our present understanding of plant systematics and phylogenetic reconstruction, and identify the possible biotechnological applications of the species belonging to the studied taxa.     

A variable minisatellite sequence in the chloroplast genome of Sorbus L. (Rosaceae: Maloideae).

The chloroplast genome is now known to be more variable than was once thought. Reports of RFLP (restriction fragment length polymorphism) and sequence variation, as well as variation in chloroplast microsatellites, are common. Here, data are presented on the variability of a minisatellite sequence in the chloroplast genome of Sorbus species. RFLP analysis of a PCR product comprising the region between the trnM and rbcL genes of nine Sorbus species identified seven size variants. Sequencing revealed the observed size polymorphism to be due to differences in the number of copies of an imperfect 9-bp motif. A more intensive survey of the variability of the minisatellite was undertaken in populations of Sorbus aucuparia. The potential uses of such regions in chloroplast DNA are discussed and a possible mechanism for the evolution of the minisatellite is presented.     

Complete mitochondrial genome of the longtooth grouper Epinephelus bruneus (Perciformes, Serranidae)

We determined the complete nucleotide sequence of the mitochondrial (mt) genome for the longtooth grouper, Epinephelus bruneus (Perciformes, Serranidae). This mt genome, consisting of 16,686 base pairs (bp), encoded genes for 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a noncoding control region as those found in other vertebrates, with the gene order identical to that of typical vertebrates. A major noncoding region between the trnP and trnF genes (991 bp) was considered to be the control region (D-loop). Within this sequence, 22 copies of a 17-bp tandem repeat element, 5'-TGATATTACATATATGC-3', were identified in the control region unlike previous reported Epinephelus species.     

Phylogenetic relationships in Monanthes (Crassulaceae) based on morphological,chloroplast and nuclear DNA variation

Phylogenetic relationships of all 10 recognized taxa of the genus Monanthes which is endemic to the Canary Islands and Salvage Islands, were investigated using the four data sets: morphology, sequences of the chloroplast DNA trnL (UAA) — trnF (GAA) intergenic spacer, ITS2 sequences of the nuclear ribosomal region and Random Amplified Polymorphic DNAs (RAPDs). In contrast to the molecular data, the morphological data were internally inconsistent which probably resulted from parallel or convergent evolution of morphological characters. The molecular data sets indicated that the genus is not monophyletic due to inclusion of the annual M. icterica which is the putative sister taxon of Aichryson, that M. polyphylla is the sister taxon of the perennial species of the genus, that M. muralis from Hierro is of allotetraploid origin and that M. lowei, M. minima, M. brachycaulos, M. laxiflora, M. anagensis and M. muralis from La Palma are closely related. Combined ITS2 sequences, trnL — trnF sequences and morphological data indicated that the relationships among three types of perennial growth forms, i.e. tiny rosettes, small, branched shrubs and diffuse branches shrublets, are highly dependent on the outgroup used. After deletion of the most distant outgroup and a taxon of alleged hybrid origin (M. muralis) relationships among the growth forms of Monanthes still could not be consistently resolved.     

喜马红景天叶绿体基因组特征及其系统发育分析

以青藏高原特有药用植物——喜马红景天（Rhodiola himalensis）为试验材料,利用高通量测序技术对喜马红景天进行叶绿体基因组测序、组装和注释,获得完整的叶绿体基因组。结果显示：喜马红景天叶绿体基因组全长为151 074 bp,GC含量为37.8%,具有1个长单拷贝区、1个短单拷贝区和1对反向重复区的典型四分体结构,其序列长度分别为82 309、17 017、25 874 bp;叶绿体基因组共编码130个基因,其中编码蛋白的基因86个、编码tRNA的基因37个、编码rRNA的基因7个;叶绿体基因组共检测出25 513个密码子,其中编码亮氨酸（Leu）的密码子占比最大;喜马红景天IRa和IRb区的rps19和ycf1基因缺失,长单拷贝区的trnH基因收缩;喜马红景天与圣地红景天（R. sacra）亲缘关系最近;短单拷贝区域的单核苷酸多态性（SNP）变异频率最高。本研究报道了喜马红景天的叶绿体基因组,并对其进行了组装、注释和序列分析,为今后开展喜马红景天的遗传多样性研究和合理开发利用提供理论依据。     

Hairpins involving both inverted and direct repeats are associated with homoplasious indels in non-coding chloroplast DNA of Taraxacum (Lactuceae: Asteraceae).

Sequence variation in 2.2 kb of non-coding regions of the chloroplast genome of eight dandelions (Taraxacum: Lactuceae) from Asia and Europe is interpreted in the light of the phylogenetic signal of base substitutions vs. indels (insertions-deletions). The four non-coding regions displayed a total of approximately 30 structural mutations of which 9 are potentially phylogenetically informative. Insertions, deletions, and an inversion were found that involved consecutive stretches of up to 172 bases. When compared to phylogenetic relationships of the chloroplast genomes based on nucleotide substitutions only, many homoplasious indels (33%) were detected that differed considerably in length and did not comprise simple sequence repeats typically associated with replication slippage. Though many indels in the intergenic spacers were associated with direct repeats, frequently, the variable stretches participated in inverted repeat stabilized hairpins. In each intergenic spacer or intron examined, nucleotide stretches ranging from 30 to 60 bp were able to fold into stabilized secondary structures. When these indels were homoplasious, they always ranked among the most stabilized hairpins in the non-coding regions. The association of higher order structures that involve both classes of repeats and parallel structural mutations in hot spot regions of the chloroplast genome can be used to differentiate among mutations that differ in phylogenetic reliability.     

Population genetics of the Y chromosome of Drosophila melanogaster: rDNA variation and phenotypic correlates

One means of examining the evolutionary significance of molecular variation on the Y chromosome is to identify phenotypes specifically affected by Y-linked genes, and to quantify the phenotypic variation and its correlation to the molecular variation. The functional importance of the Y-linked array of rRNA genes is demonstrated by the ability of Y chromosome to rescue X-linked bobbed lethal alleles, whose lethality is seen in homozygous females. Because low numbers of X-linked rDNA gene copies result in increased developmental time and shortened bristles, and because there is considerable natural variation in Y-linked copy number, a careful examination of Y-linked variation in these two traits may uncover a mode of selection acting on the multigene family. In this study, 36 Y-chromosome replacement lines were tested to detect subtle variation in bristle phenotypes and developmental rates. Correlations among these traits, rDNA gene copy number, and intergenic sequence length were quantified. The absence of significant correlations between phenotypic characters and rDNA copy number of intergenic sequence length suggests that the extant molecular variation in Y-linked rDNA can have at most very small selective effects.     

Gene copy number variation and its significance in cyanobacterial phylogeny

ABSTRACT: BACKGROUND: In eukaryotes, variation in gene copy numbers is often associated with deleterious effects, but may also have positive effects. For prokaryotes, studies on gene copy number variation are rare. Previous studies have suggested that high numbers of rRNA gene copies can be advantageous in environments with changing resource availability, but further association of gene copies and phenotypic traits are not documented. We used one of the morphologically most diverse prokaryotic phyla to test whether numbers of gene copies are associated with levels of cell differentiation. RESULTS: We implemented a search algorithm that identified 44 genes with highly conserved copies across 22 fully sequenced cyanobacterial taxa. For two very basal cyanobacterial species, Gloeobacter violaceus and a thermophilic Synechococcus species, distinct phylogenetic positions previously found were supported by identical protein coding gene copy numbers. Furthermore, we found that increased ribosomal gene copy numbers showed a strong correlation to cyanobacteria capable of terminal cell differentiation. Additionally, we detected extremely low variation of 16S rRNA sequence copies within the cyanobacteria. We compared our results for 16S rRNA to three other eubacterial phyla (Chroroflexi, Spirochaetes and Bacteroidetes). Based on Bayesian phylogenetic inference and the comparisons of genetic istances, we could confirm that cyanobacterial 16S rRNA paralogs and orthologs show significantly stronger conservation than found in other eubacterial phyla. Conclusions: A higher number of ribosomal operons could potentially provide an advantage to terminally differentiated cyanobacteria. Furthermore, we suggest that 16S rRNA gene copies in cyanobacteria are homogenized by both concerted evolution and purifying selection. In addition, the small ribosomal subunit in cyanobacteria appears to evolve at extraordinary slow evolutionary rates, an observation that has been made previously for morphological characteristics of cyanobacteria.