Localization of high level of sequence conservation and divergence regions in cotton

In a previous study, we observed that the variations in chromosome size are due to uneven expansion and contraction by comparing the structures and sizes of a pair of homoeologous high-resolution cytogenetic maps of chromosomes 12A and 12D in tetraploid cotton. To reveal the variation at the sequence level, in the present paper, we sequenced two pairs of homoeologous bacterial artificial chromosomes derived from high- to low-variable genomic regions. Comparisons of their sequence variations confirmed that the highly conserved and divergent sequences existed in the distal and pericentric regions, e.g., high- and low-variable genome size regions in these two pairs of cotton homoeologous chromosomes. Sequence analysis also confirmed that the differential accumulation of Gossypium retrotransposable gypsy-like element (Gorge3) accounted for the main contributions for the size difference between the pericentric regions. By fluorescence in situ hybridization analysis, we found that Gorge3 has a bias distribution in the A_T/A proximal regions and is associated with the heterochromatin along the chromosomes in the entire Gossypium genome. These results indicate that, between A_T/A and D_T/D genomes, the distal and pericentric regions usually possess high level of sequence conservation and divergence, respectively, in cotton.     

A High Density Consensus Genetic Map of Tetraploid Cotton That Integrates Multiple Component Maps through Molecular Marker Redundancy Check

A consensus genetic map of tetraploid cotton was constructed using six high-density maps and after the integration of a sequence-based marker redundancy check. Public cotton SSR libraries (17,343 markers) were curated for sequence redundancy using 90% as a similarity cutoff. As a result, 20% of the markers (3,410) could be considered as redundant with some other markers. The marker redundancy information had been a crucial part of the map integration process, in which the six most informative interspecific Gossypium hirsutum×G. barbadense genetic maps were used for assembling a high density consensus (HDC) map for tetraploid cotton. With redundant markers being removed, the HDC map could be constructed thanks to the sufficient number of collinear non-redundant markers in common between the component maps. The HDC map consists of 8,254 loci, originating from 6,669 markers, and spans 4,070 cM, with an average of 2 loci per cM. The HDC map presents a high rate of locus duplications, as 1,292 markers among the 6,669 were mapped in more than one locus. Two thirds of the duplications are bridging homoeologous A_T and D_T chromosomes constitutive of allopolyploid cotton genome, with an average of 64 duplications per A_T/D_T chromosome pair. Sequences of 4,744 mapped markers were used for a mutual blast alignment (BBMH) with the 13 major scaffolds of the recently released Gossypium raimondii genome indicating high level of homology between the diploid D genome and the tetraploid cotton genetic map, with only a few minor possible structural rearrangements. Overall, the HDC map will serve as a valuable resource for trait QTL comparative mapping, map-based cloning of important genes, and better understanding of the genome structure and evolution of tetraploid cotton.     

Chromosomal Locations of 5S and 45S rDNA in Gossypium Genus and Its Phylogenetic Implications Revealed by FISH

We investigated the locations of 5S and 45S rDNA in Gossypium diploid A, B, D, E, F, G genomes and tetraploid genome (AD) using multi-probe fluorescent in situ hybridization (FISH) for evolution analysis in Gossypium genus. The rDNA numbers and sizes, and synteny relationships between 5S and 45S were revealed using 5S and 45S as double-probe for all species, and the rDNA-bearing chromosomes were identified for A, D and AD genomes with one more probe that is single-chromosome-specific BAC clone from G. hirsutum (A₁D₁). Two to four 45S and one 5S loci were found in diploid-species except two 5S loci in G . incanum (E₄), the same as that in tetraploid species. The 45S on the 7th and 9th chromosomes and the 5S on the 9th chromosomes seemed to be conserved in A, D and AD genomes. In the species of B, E, F and G genomes, the rDNA numbers, sizes, and synteny relationships were first reported in this paper. The rDNA pattern agrees with previously reported phylogenetic history with some disagreements. Combined with the whole-genome sequencing data from G . raimondii (D₅) and the conserved cotton karyotype, it is suggested that the expansion, decrease and transposition of rDNA other than chromosome rearrangements might occur during the Gossypium evolution.     

Complete Chloroplast Genome of Tanaecium tetragonolobum: The First Bignoniaceae Plastome

Bignoniaceae is a Pantropical plant family that is especially abundant in the Neotropics. Members of the Bignoniaceae are diverse in many ecosystems and represent key components of the Tropical flora. Despite the ecological importance of the Bignoniaceae and all the efforts to reconstruct the phylogeny of this group, whole chloroplast genome information has not yet been reported for any members of the family. Here, we report the complete chloroplast genome sequence of Tanaecium tetragonolobum (Jacq.) L.G. Lohmann, which was reconstructed using de novo and referenced-based assembly of single-end reads generated by shotgun sequencing of total genomic DNA in an Illumina platform. The gene order and organization of the chloroplast genome of T. tetragonolobum exhibits the general structure of flowering plants, and is similar to other Lamiales chloroplast genomes. The chloroplast genome of T. tetragonolobum is a circular molecule of 153,776 base pairs (bp) with a quadripartite structure containing two single copy regions, a large single copy region (LSC, 84,612 bp) and a small single copy region (SSC, 17,586 bp) separated by inverted repeat regions (IRs, 25,789 bp). In addition, the chloroplast genome of T. tetragonolobum has 38.3% GC content and includes 121 genes, of which 86 are protein-coding, 31 are transfer RNA, and four are ribosomal RNA. The chloroplast genome of T. tetragonolobum presents a total of 47 tandem repeats and 347 simple sequence repeats (SSRs) with mononucleotides being the most common and di-, tri-, tetra-, and hexanucleotides occurring with less frequency. The results obtained here were compared to other chloroplast genomes of Lamiales available to date, providing new insight into the evolution of chloroplast genomes within Lamiales. Overall, the evolutionary rates of genes in Lamiales are lineage-, locus-, and region-specific, indicating that the evolutionary pattern of nucleotide substitution in chloroplast genomes of flowering plants is complex. The discovery of tandem repeats within T. tetragonolobum and the presence of divergent regions between chloroplast genomes of Lamiales provides the basis for the development of markers at various taxonomic levels. The newly developed markers have the potential to greatly improve the resolution of molecular phylogenies.     

Genome organization in dicots. II. Arabidopsis as a ’bridging species’ to resolve genome evolution events among legumes

Analysis of molecular linkage groups within the soybean (Glycine max L. Merr.) genome reveals many homologous regions, reflecting the ancient polyploidy of soybean. The fragmented arrangement of the duplicated regions suggests that extensive rearrangements, as well as additional duplications, have occurred since the initial polyploidization event. In this study we used comparisons between homoeologous regions in soybean, and the homologous regions in the related diploids Phaseolus vulgaris and Vigna radiata, to elucidate the evolutionary history of the three legume genomes. Our results show that there is not only conservation of large regions of the genomes but that these conserved linkage blocks are also represented twice in the soybean genome. To gain a better understanding of the process of genome evolution in dicots, molecular comparisons have been extended to another well-studied species, Arabidopsis thaliana. Interestingly, the conserved regions we identified in the legume species are also relatively conserved in Arabidopsis. Our results suggest that there is conservation of blocks of DNA between species as distantly related as legumes and brassicas, representing 90 million years of divergence. We also present evidence for an additional, presumably earlier, genome duplication in soybean. These duplicated regions were only recognized by using Arabidopsis as a ’bridging’ species in the genome comparisons. Received: 10 October 2000 / Accepted: 13 January 2001     

Intergenomic single nucleotide polymorphisms as a tool for bacterial artificial chromosome contig building of homoeologous Brassica napus regions

Background

Homoeologous sequences pose a particular challenge if bacterial artificial chromosome (BAC) contigs shall be established for specific regions of an allopolyploid genome. Single nucleotide polymorphisms (SNPs) differentiating between homoeologous genomes (intergenomic SNPs) may represent a suitable screening tool for such purposes, since they do not only identify homoeologous sequences but also differentiate between them.

Results

Sequence alignments between Brassica rapa (AA) and Brassica oleracea (CC) sequences mapping to corresponding regions on chromosomes A1 and C1, respectively were used to identify single nucleotide polymorphisms between the A and C genomes. A large fraction of these polymorphisms was also present in Brassica napus (AACC), an allopolyploid species that originated from hybridisation of A and C genome species. Intergenomic SNPs mapping throughout homoeologous chromosome segments spanning approximately one Mbp each were included in Illumina’s GoldenGate® Genotyping Assay and used to screen multidimensional pools of a Brassica napus bacterial artificial chromosome library with tenfold genome coverage. Based on the results of 50 SNP assays, a BAC contig for the Brassica napus A subgenome was established that spanned the entire region of interest. The C subgenome region was represented in three BAC contigs.

Conclusions

This proof-of-concept study shows that sequence resources of diploid progenitor genomes can be used to deduce intergenomic SNPs suitable for multiplex polymerase chain reaction (PCR)-based screening of multidimensional BAC pools of a polyploid organism. Owing to their high abundance and ease of identification, intergenomic SNPs represent a versatile tool to establish BAC contigs for homoeologous regions of a polyploid genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-560) contains supplementary material, which is available to authorized users.     

Sequencing and Utilization of the Gossypium Genomes

Revealing the genetic underpinnings of cotton productivity will require understanding both the prehistoric evolution of spinnable fibers, and the results of independent domestication processes in both the Old and New Worlds. Progress toward a reference sequence for the smallest Gossypium genome is a logical stepping-stone toward revealing diversity in the remaining seven genomes (A, B, C, E, F, G, K) that permitted Gossypium species to adapt to a wide range of ecosystems in warmer arid regions of the world, and toward identifying the emergent properties that account for the superior productivity and quality of tetraploid cottons. The greatest challenge facing the cotton community is not genome sequencing per se but the conversion of sequence to knowledge.     

Molecular evolution and phylogeny of the RPB2 gene in the genus Hordeum

Background and Aims

It is known that the miniature inverted-repeat terminal element (MITE) preferentially inserts into low-copy-number sequences or genic regions. Characterization of the second largest subunit of low-copy nuclear RNA polymerase II (RPB2) has indicated that MITE and indels have shaped the homoeologous RPB2 loci in the St and H genome of Eymus species in Triticeae. The aims of this study was to determine if there is MITE in the RPB2 gene in Hordeum genomes, and to compare the gene evolution of RPB2 with other diploid Triticeae species. The sequences were used to reconstruct the phylogeny of the genus Hordeum.

Methods

RPB2 regions from all diploid species of Hordeum, one tetraploid species (H. brevisubulatum) and ten accessions of diploid Triticeae species were amplified and sequenced. Parsimony analysis of the DNA dataset was performed in order to reveal the phylogeny of Hordeum species.

Key Results

MITE was detected in the Xu genome. A 27–36 bp indel sequence was found in the I and Xu genome, but deleted in the Xa and some H genome species. Interestingly, the indel length in H genomes corresponds well to their geographical distribution. Phylogenetic analysis of the RPB2 sequences positioned the H and Xa genome in one monophyletic group. The I and Xu genomes are distinctly separated from the H and Xa ones. The RPB2 data also separated all New World H genome species except H. patagonicum ssp. patagonicum from the Old World H genome species.

Conclusions

MITE and large indels have shaped the RPB2 loci between the Xu and H, I and Xa genomes. The phylogenetic analysis of the RPB2 sequences confirmed the monophyly of Hordeum. The maximum-parsimony analysis demonstrated the four genomes to be subdivided into two groups.Key words: Molecular evolution, RPB2, Hordeum, transposable element, phylogeny     

Comparative RFLP mapping of meadow and tall fescue

 Molecular markers based on restriction fragment length polymorphism (RFLP) were used to construct a genetic linkage map in diploid meadow fescue, Festuca pratensis Huds. (2n=2x=14, genomic designation PP), and to compare its genomic relationship with a related species, hexaploid tall fescue (Festuca arundinacea Schreb.; 2n=6x=42, PPG₁G₁G₂G₂). Using a collection of 66 tall-fescue (heterologous) markers, an RFLP linkage map was constructed in F. pratensis. This map, which has a total length of 280.1 cM, includes seven linkage groups. A comparison of 33 markers that were mapped in both F. pratensis and F. arundinacea detected highly conserved linkage groups between these two species. Our data are consistent with the proposal that one of the genomes of F. arundinacea was derived from F. pratensis. However, since significant changes in marker sequences, map distances, and homoeologous linkage groups were also detected between the two species, it appears that the P genome diverged substantially during evolution from the diploid to the hexaploid Festuca. Received: 23 May 1997 / Accepted: 15 January 1998     

Isolation and identification of Triticum aestivum L. em. Thell. cv Chinese Spring-T. peregrinum Hackel disomic chromosome addition lines

Analyses of RFLPs, isozymes, morphological markers and chromosome pairing were used to isolate 12 Triticum aestivum cv Chinese Spring (genomes A, B, and D)-T. peregrinum (genomes S^v and U^v) disomic chromosome addition lines. The evidence obtained indicates that each of the 12 lines contains an intact pair of T. peregrinum chromosomes. One monosomic addition line, believed to contain an intact 6S^v chromosome, was also isolated. A CS-7U^v chromosome addition line was not obtained. Syntenic relationships in common with the standard Triticeae arrangement were found for five of the seven S^v genome chromosomes. The exceptions were 4S^v and 7S^v. A reciprocal translocation exists between 4S¹ and 7S¹ in T. longissimum and evidence was obtained that the same translocation exists in T. peregrinum. In contrast, evidence for syntenic relationships in common with the standard Triticeae arrangements were found for only one U^v chromosome of T. peregrinum.; namely, chromosome 2U^v. All other U^v genome chromosomes are involved in at least one translocation, and the same translocations were found in the U genome of T. umbellulatum. Evidence was also obtained indicating that the centromeric regions of 4U and 4U^v are homoeologous to the centromeric regions of Triticeae homoeologous group-6 chromosomes, that the centromeric regions of 6U and 6U^v are homoeologous to the centromeric regions of group-4 chromosomes, and that 4U and 4U^v are more closely related overall to Triticeae homoeologous group-6 chromosomes than they are to group-4 chromosomes.     

Genome mapping of polyploid tall fescue (Festuca arundinacea Schreb.) with RFLP markers

Genetic mapping using molecular markers such as restriction fragment length polymorphisms (RFLPs) has become a powerful tool for plant geneticists and breeders. Like many economically important polyploid plant species, detailed genetic studies of hexaploid tall fescue (Festuca arundinacea Schreb.) are complicated, and no genetic map has been established. We report here the first tall fescue genetic map. This map was generated from an F₂ population of HD28-56 by Kentucky-31 and contains 108 RFLP markers. Although the two parental plants were heterozygous, the perennial and tillering growth habit, high degree of RFLP, and disomic inheritance of tall fescue enabled us to identify the segregating homologous alleles. The map covers 1274 cM on 19 linkage groups with an average of 5 loci per linkage group (LG) and 17.9 cM between loci. Mapping the homoeologous loci detected by the same probe allowed us to identify five homoeologous groups within which the gene orders were found to be generally conserved among homoeologous chromosomes. An exception was homoeologous group 5, in which only 2 of the 3 homoeologous chromosomes were identified. Using 12 genome-specific probes, we were able to assign several linkage groups to one of the three genomes (PG₁G₂) in tall fescue. All the loci detected by the 11 probes specific to the G₁ and/or G₂ genomes, with one exception, identified loci located on 4 chromosomes of two homoeologous groups (LG2a, LG2c, LG3a, and LG3c). A P-genome-specific probe was used to map a locus on LG5c. Comparative genome mapping with maize probes indicated that homoeologous group 3 and 2 chromosomes in tall fescue corresponded to maize chromosome 1. Difficulties and advantages of applying RFLP technology in polyploids with high levels of heterozygosity are discussed.Journal Series No. 12, 190     

Identification and Characterization of Microsatellite Markers Derived from the Whole Genome Analysis of Taenia solium

Background

Infections with Taenia solium are the most common cause of adult acquired seizures worldwide, and are the leading cause of epilepsy in developing countries. A better understanding of the genetic diversity of T. solium will improve parasite diagnostics and transmission pathways in endemic areas thereby facilitating the design of future control measures and interventions. Microsatellite markers are useful genome features, which enable strain typing and identification in complex pathogen genomes. Here we describe microsatellite identification and characterization in T. solium, providing information that will assist in global efforts to control this important pathogen.

Methods

For genome sequencing, T. solium cysts and proglottids were collected from Huancayo and Puno in Peru, respectively. Using next generation sequencing (NGS) and de novo assembly, we assembled two draft genomes and one hybrid genome. Microsatellite sequences were identified and 36 of them were selected for further analysis. Twenty T. solium isolates were collected from Tumbes in the northern region, and twenty from Puno in the southern region of Peru. The size-polymorphism of the selected microsatellites was determined with multi-capillary electrophoresis. We analyzed the association between microsatellite polymorphism and the geographic origin of the samples.

Results

The predicted size of the hybrid (proglottid genome combined with cyst genome) T. solium genome was 111 MB with a GC content of 42.54%. A total of 7,979 contigs (>1,000 nt) were obtained. We identified 9,129 microsatellites in the Puno-proglottid genome and 9,936 in the Huancayo-cyst genome, with 5 or more repeats, ranging from mono- to hexa-nucleotide. Seven microsatellites were polymorphic and 29 were monomorphic within the analyzed isolates. T. solium tapeworms were classified into two genetic groups that correlated with the North/South geographic origin of the parasites.

Conclusions/Significance

The availability of draft genomes for T. solium represents a significant step towards the understanding the biology of the parasite. We report here a set of T. solium polymorphic microsatellite markers that appear promising for genetic epidemiology studies.     

Nucleotide diversity in the two co-resident genomes of allopolyploid cotton

Genetic diversity within and among populations lies at the heart of evolution. Unraveling the extent to which each intrinsic or extrinsic factor determines levels of diversity among genes, populations, and species is challenging, given the difficulty of isolating any single potentially important variable from all others. Allopolyploid species provide an opportunity to disentangle external and intrinsic factors, as the two (or more) homoeologous genomes co-occur in the same nucleus, often exhibiting high collinearity along homoeologous chromosomes. Here we evaluate the pace of molecular evolution and intraspecific, intragenomic diversity in two species of allopolyploid Gossypium, G. hirsutum and G. barbadense, using several hundred genes sequenced from multiple accessions of each species. Genic diversity in both species is low, having been influenced both by the polyploid bottleneck and a domestication bottleneck (for cultivated accessions), but with a directional bias in homoeolog diversity favoring the same genome in both allopolyploids. Total diversity is remarkably similar for the two homoeologous genomes overall, but the two copies of many gene pairs have accumulated statistically different diversity levels, and in a biased fashion with respect to genome. Domesticated accessions show reduced diversity in both genomes, as expected, but with a much greater reduction in one of the two homoeologous genomes. Furthermore, this biased reduction affects opposite homoeologous genomes in the two species. Interspecific introgression has played a role in shaping diversity within each species. Introgression was only detected for certain accessions, and only from G. barbadense into G. hirsutum in one of the two co-resident genomes.     

Analysis of [Gossypium capitis-viridis × (G.hirsutum × G.australe)2] Trispecific Hybrid and Selected Characteristics

Speciation is always a contentious and challenging issue following with the presence of gene flow. In Gossypium, there are many valuable resources and wild diploid cotton especially C and B genome species possess some excellent traits which cultivated cotton always lacks. In order to explore character transferring rule from wild cotton to upland tetraploid cotton, the [G. capitis-viridis × (G. hirsutum × G. australe)²] triple hybrid was synthesized by interspecies hybridization and chromosome doubling. Morphology comparisons were measured among this hybrid and its parents. It showed that trispecific hybrid F₁ had some intermediate morphological characters like leaf style between its parents and some different characters from its parents, like crawl growth characteristics and two kind flower color. It is highly resistant to insects comparing with other cotton species by four year field investigation. By cytogenetic analysis, triple hybrid was further confirmed by meiosis behavior of pollen mother cells. Comparing with regular meiosis of its three parents, it was distinguished by the occurrence of polyads with various numbers of unbalanced microspores and finally generating various abnormal pollen grains. All this phenomenon results in the sterility of this hybrid. This hybrid was further identified by SSR marker from DNA molecular level. It showed that 98 selected polymorphism primers amplified effective bands in this hybrids and its parents. The genetic proportion of three parents in this hybrid is 47.8% from G. hirsutum, 14.3% from G. australe, 7.0% from G. capitis-viridis, and 30.9% recombination bands respectively. It was testified that wild genetic material has been transferred into cultivated cotton and this new germplasm can be incorporated into cotton breeding program.     

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.     

Spatio‐temporal patterns of genome evolution in allotetraploid species of the genus Oryza

Despite knowledge that polyploidy is widespread and a major evolutionary force in flowering plant diversification, detailed comparative molecular studies on polyploidy have been confined to only a few species and families. The genus Oryza is composed of 23 species that are classified into ten distinct ‘genome types’ (six diploid and four polyploid), and is emerging as a powerful new model system to study polyploidy. Here we report the identification, sequence and comprehensive comparative annotation of eight homoeologous genomes from a single orthologous region (Adh1–Adh2) from four allopolyploid species representing each of the known Oryza genome types (BC, CD, HJ and KL). Detailed comparative phylogenomic analyses of these regions within and across species and ploidy levels provided several insights into the spatio‐temporal dynamics of genome organization and evolution of this region in ‘natural’ polyploids of Oryza. The major findings of this study are that: (i) homoeologous genomic regions within the same nucleus experience both independent and parallel evolution, (ii) differential lineage‐specific selection pressures do not occur between polyploids and their diploid progenitors, (iii) there have been no dramatic structural changes relative to the diploid ancestors, (iv) a variation in the molecular evolutionary rate exists between the two genomes in the BC complex species even though the BC and CD polyploid species appear to have arisen <2 million years ago, and (v) there are no clear distinctions in the patterns of genome evolution in the diploid versus polyploid species.     

Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes

BackgroundSNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes.ResultsWe construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis.ConclusionsWe report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0678-1) contains supplementary material, which is available to authorized users.     

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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-374) contains supplementary material, which is available to authorized users.     

The Genome Sequences of Cellulomonas fimi and “Cellvibrio gilvus” Reveal the Cellulolytic Strategies of Two Facultative Anaerobes,Transfer of “Cellvibrio gilvus” to the Genus Cellulomonas,and Proposal of Cellulomonas gilvus sp. nov

Actinobacteria in the genus Cellulomonas are the only known and reported cellulolytic facultative anaerobes. To better understand the cellulolytic strategy employed by these bacteria, we sequenced the genome of the Cellulomonas fimi ATCC 484^T. For comparative purposes, we also sequenced the genome of the aerobic cellulolytic “Cellvibrio gilvus” ATCC 13127^T. An initial analysis of these genomes using phylogenetic and whole-genome comparison revealed that “Cellvibrio gilvus” belongs to the genus Cellulomonas. We thus propose to assign “Cellvibrio gilvus” to the genus Cellulomonas. A comparative genomics analysis between these two Cellulomonas genome sequences and the recently completed genome for Cellulomonas flavigena ATCC 482^T showed that these cellulomonads do not encode cellulosomes but appear to degrade cellulose by secreting multi-domain glycoside hydrolases. Despite the minimal number of carbohydrate-active enzymes encoded by these genomes, as compared to other known cellulolytic organisms, these bacteria were found to be proficient at degrading and utilizing a diverse set of carbohydrates, including crystalline cellulose. Moreover, they also encode for proteins required for the fermentation of hexose and xylose sugars into products such as ethanol. Finally, we found relatively few significant differences between the predicted carbohydrate-active enzymes encoded by these Cellulomonas genomes, in contrast to previous studies reporting differences in physiological approaches for carbohydrate degradation. Our sequencing and analysis of these genomes sheds light onto the mechanism through which these facultative anaerobes degrade cellulose, suggesting that the sequenced cellulomonads use secreted, multidomain enzymes to degrade cellulose in a way that is distinct from known anaerobic cellulolytic strategies.     

Interspecies Insertion Polymorphism Analysis Reveals Recent Activity of Transposable Elements in Extant Coelacanths

Coelacanths are lobe-finned fish represented by two extant species, Latimeria chalumnae in South Africa and Comoros and L. menadoensis in Indonesia. Due to their intermediate phylogenetic position between ray-finned fish and tetrapods in the vertebrate lineage, they are of great interest from an evolutionary point of view. In addition, extant specimens look similar to 300 million-year-old fossils; because of their apparent slowly evolving morphology, coelacanths have been often described as « living fossils ». As an underlying cause of such a morphological stasis, several authors have proposed a slow evolution of the coelacanth genome. Accordingly, sequencing of the L. chalumnae genome has revealed a globally low substitution rate for protein-coding regions compared to other vertebrates. However, genome and gene evolution can also be influenced by transposable elements, which form a major and dynamic part of vertebrate genomes through their ability to move, duplicate and recombine. In this work, we have searched for evidence of transposition activity in coelacanth genomes through the comparative analysis of orthologous genomic regions from both Latimeria species. Comparison of 5.7 Mb (0.2%) of the L. chalumnae genome with orthologous Bacterial Artificial Chromosome clones from L. menadoensis allowed the identification of 27 species-specific transposable element insertions, with a strong relative contribution of CR1 non-LTR retrotransposons. Species-specific homologous recombination between the long terminal repeats of a new coelacanth endogenous retrovirus was also detected. Our analysis suggests that transposon activity is responsible for at least 0.6% of genome divergence between both Latimeria species. Taken together, this study demonstrates that coelacanth genomes are not evolutionary inert: they contain recently active transposable elements, which have significantly contributed to post-speciation genome divergence in Latimeria.