Nucleomorph ribosomal DNA and telomere dynamics in chlorarachniophyte algae

PhyloChromoMap, a Tool for Mapping Phylogenomic History along Chromosomes, Reveals the Dynamic Nature of Karyotype Evolution in Plasmodium falciparum

The genome of Plasmodium falciparum, the causative agent of malaria in Africa, has been extensively studied since it was first fully sequenced in 2002. However, many open questions remain, including understanding the chromosomal context of molecular evolutionary changes (e.g., relationship between chromosome map and phylogenetic conservation, patterns of gene duplication, and patterns of selection). Here, we present PhyloChromoMap, a method that generates a phylogenomic map of chromosomes from a custom-built bioinformatics pipeline. Using P. falciparum 3D7 as a model, we analyze 2,116 genes with homologs in up to 941 diverse eukaryotic, bacterial and archaeal lineages. We estimate the level of conservation along chromosomes based on conservation across clades, and identify “young” regions (i.e., those with recent or fast evolving genes) that are enriched in subtelomeric regions as compared with internal regions. We also demonstrate that patterns of molecular evolution for paralogous genes differ significantly depending on their location as younger paralogs tend to be found in subtelomeric regions whereas older paralogs are enriched in internal regions. Combining these observations with analyses of synteny, we demonstrate that subtelomeric regions are actively shuffled among chromosome ends, which is consistent with the hypothesis that these regions are prone to ectopic recombination. We also assess patterns of selection by comparing dN/dS ratios of gene family members in subtelomeric versus internal regions, and we include the important antigenic gene family var. These analyses illustrate the highly dynamic nature of the karyotype of P. falciparum, and provide a method for exploring genome dynamics in other lineages.

Nucleomorph Genome Sequences of Two Chlorarachniophytes, Amorphochlora amoebiformis and Lotharella vacuolata

Many algal groups acquired complex plastids by the uptake of green and red algae through multiple secondary endosymbioses. As a result of gene loss and transfer during the endosymbiotic processes, algal endosymbiont nuclei disappeared in most cases. However, chlorarachniophytes and cryptophytes still possess a relict nucleus, so-called the nucleomorph, of the green and red algal endosymbiont, respectively. Nucleomorph genomes are an interesting and suitable model to study the reductive evolution of endosymbiotically derived genomes. To date, nucleomorph genomes have been sequenced in four cryptophyte species and two chlorarachniophyte species, including Bigelowiella natans (373 kb) and Lotharella oceanica (610 kb). In this study, we report complete nucleomorph genome sequences of two chlorarachniophytes, Amorphochlora amoebiformis and Lotharella vacuolata, to gain insight into the reductive evolution of nucleomorph genomes in the chlorarachniophytes. The nucleomorph genomes consist of three chromosomes totaling 374 and 432 kb in size in A. amoebiformis and L. vacuolata, respectively. Comparative analyses among four chlorarachniophyte nucleomorph genomes revealed that these sequences share 171 function-predicted genes (86% of total 198 function-predicted nucleomorph genes), including the same set of genes encoding 17 plastid-associated proteins, and no evidence of a recent nucleomorph-to-nucleus gene transfer was found. This suggests that chlorarachniophyte nucleomorph genomes underwent most of their reductive evolution prior to the radiation of extent members of the group. However, there are slight variations in genome size, GC content, duplicated gene number, and subtelomeric regions among the four nucleomorph genomes, suggesting that the genomes might be undergoing changes that do not affect the core functions in each species.

Nucleomorph and plastid genome sequences of the chlorarachniophyte Lotharella oceanica: convergent reductive evolution and frequent recombination in nucleomorph-bearing algae

BACKGROUND

Nucleomorphs are residual nuclei derived from eukaryotic endosymbionts in chlorarachniophyte and cryptophyte algae. The endosymbionts that gave rise to nucleomorphs and plastids in these two algal groups were green and red algae, respectively. Despite their independent origin, the chlorarachniophyte and cryptophyte nucleomorph genomes share similar genomic features such as extreme size reduction and a three-chromosome architecture. This suggests that similar reductive evolutionary forces have acted to shape the nucleomorph genomes in the two groups. Thus far, however, only a single chlorarachniophyte nucleomorph and plastid genome has been sequenced, making broad evolutionary inferences within the chlorarachniophytes and between chlorarachniophytes and cryptophytes difficult. We have sequenced the nucleomorph and plastid genomes of the chlorarachniophyte Lotharella oceanica in order to gain insight into nucleomorph and plastid genome diversity and evolution.

RESULTS

The L. oceanica nucleomorph genome was found to consist of three linear chromosomes totaling ~610 kilobase pairs (kbp), much larger than the 373 kbp nucleomorph genome of the model chlorarachniophyte Bigelowiella natans. The L. oceanica plastid genome is 71 kbp in size, similar to that of B. natans. Unexpectedly long (~35 kbp) sub-telomeric repeat regions were identified in the L. oceanica nucleomorph genome; internal multi-copy regions were also detected. Gene content analyses revealed that nucleomorph house-keeping genes and spliceosomal intron positions are well conserved between the L. oceanica and B. natans nucleomorph genomes. More broadly, gene retention patterns were found to be similar between nucleomorph genomes in chlorarachniophytes and cryptophytes. Chlorarachniophyte plastid genomes showed near identical protein coding gene complements as well as a high level of synteny.

CONCLUSIONS

We have provided insight into the process of nucleomorph genome evolution by elucidating the fine-scale dynamics of sub-telomeric repeat regions. Homologous recombination at the chromosome ends appears to be frequent, serving to expand and contract nucleomorph genome size. The main factor influencing nucleomorph genome size variation between different chlorarachniophyte species appears to be expansion-contraction of these telomere-associated repeats rather than changes in the number of unique protein coding genes. The dynamic nature of chlorarachniophyte nucleomorph genomes lies in stark contrast to their plastid genomes, which appear to be highly stable in terms of gene content and synteny.

Evolutionary mobility of the ribosomal DNA array in yeasts

The ribosomal DNA (rDNA) of eukaryotes is organized as large tandem arrays. Here, we compare the genomic locations of rDNA among yeast species and show that, despite its huge size (>1 Mb), the rDNA array has moved around the genome several times within the family Saccharomycetaceae. We identify an ancestral, nontelomeric, rDNA site that is conserved across many species including Saccharomyces cerevisiae. Within the genus Lachancea, however, the rDNA apparently transposed from the ancestral site to a new site internal to a different chromosome, becoming inserted into a short intergenic region beside a tRNA gene. In at least four other yeast lineages, the rDNA moved from the ancestral site to telomeric locations. Remarkably, both the ancestral rDNA site and the new site in Lachancea are adjacent to protein-coding genes whose products maintain the specialized chromatin structure of rDNA (HMO1 and CDC14, respectively). In almost every case where the rDNA was lost from the ancestral site, the entire array disappeared without any other rearrangements in the region, leaving just an intergenic spacer of less than 2 kb. The mechanism by which this large and complex locus moves around the genome is unknown, but we speculate that it may involve the formation of double-strand DNA breaks by Fob1 protein or the formation of extrachromosomal rDNA circles.

Emergence and expansion of TFIIB-like factors in the plant kingdom

Many gene families in higher plants have expanded in number, giving rise to diverse protein paralogs with specialized biochemical functions. For instance, plant general transcription factors such as TFIIB have expanded in number and in some cases perform specialized transcriptional functions in the plant cell. To date, no comprehensive genome-wide identification of the TFIIB gene family has been conducted in the plant kingdom. To determine the extent of TFIIB expansion in plants, I used the remote homology program HHPred to search for TFIIB homologs in the plant kingdom and performed a comprehensive analysis of eukaryotic TFIIB gene families. I discovered that higher plants encode more than 10 different TFIIB-like proteins. In particular, Arabidopsis thaliana encodes 14 different TFIIB-like proteins and predicted domain architectures of the newly identified TFIIB-like proteins revealed that they have unique modular domain structures that are divergent in sequence and size. Phylogenetic analysis of selected eukaryotic organisms showed that most life forms encode three major TFIIB subfamilies that include TFIIB, Brf, Rrn7/TAF1B/MEE12 subfamilies, while all plants and some algae species encode one or two additional TFIIB-related protein subfamilies. A subset of A. thaliana GTFs have also expanded in number, indicating that GTF diversification and expansion is a general phenomenon in higher plants. Together, these findings were used to generate a model for the evolutionary history of TFIIB-like proteins in eukaryotes.