Most plastidic genomes have four regions, namely large single copy region (LSC, 80 Kb), small single copy region (SSC, 20 kb) and two inverted repeat regions (IR, 25 kb). The single copy region is separated by two IRs. This structural conservation however breaks in some plants such as Vicia faba (1) and Cryptomeria japonica (2, 3) by loss of an IR, and in Euglena gracilis that has three tandem repeats (4).
Variations among different species provide large information for the phylogenetic studies. Chloroplasts have low mutation rate with great deal of conservation in their genome size and structure, gene content and organization. Few differences have been reported in the same species, but significant differences could be detected between the different species in genome size and gene orientation (5). It has been reported that, chloroplast genes like 16S, 23S, ndhB, psbA, psbD, psaB, pasA, psbC, psbB and rbcL are appropriate to study the relationship among higher plants; ycf1, ycf2, accD, matK, rpoC2 and ndhF are more suitable to study the relationship of the close species (5).
Transplastomics have proved to be a powerful tool to improve the plant genetic architecture with high expression of the foreign protein, low risk of the pollen pollution (6) and no gene silencing. Therefore and In addition to phylogenetic analysis based on plastidic genomes, it is imperative to understand the chloroplast genome in order to logically design our next generation transplastomics. Accordingly, chloroplast genomes of many species have been sequenced (7, 8, 9, 10, 11, 12, 13, 14).
Gossypium includes 52 species that are divided to eight diploid genome A-G and K (2n=26), and one allotetraploid genome (AADD, 2n=52). G. barbadense and G. hirsutum are extensively cultivated in the world and their chloroplast genome sequences have been published. G. thurberi belongs to D genome, which is an important wild bio-source for the cotton breeding and genetic research. Present study was conducted to study and compare the complete chloroplast sequence of G. thurberi, analyses of its genome structure, gene content and organization, repeat sequence and codon usage. Meanwhile the comparison of the three sequenced cotton species was performed.
3. Materials and Methods
3.1. Chloroplast Sequence
Compleate chloroplast genome sequence of Gossypium thurberi with accession number NC_015204.1 downloaded from NCBI (http://www.ncbi.nlm.nih.gov/nuccore/?term=Gossy pium thurberi).
3.2. Genome Assembly and Gene Annotation
The available sequence was assembled by DNAman (Version 22.214.171.1248). Gene annotation was mainly performed by DOGMA (Dual Organellar Geno Me Annotator, http://dogma.ccbb.utexas.edu/; Wyman, 2004). DOGMA uses BLAST against 11 plant chloroplast database (Adiantum capillus-veneris, Arabidopsis thaliana, Chlorella vulgaris, Lotus japonicus, Marchantia polymorpha, Mesostigma viride, Nephroselmis olivacea, Nicotiana tabacum, Oenothera elata, Oryza sativa, Pinus thunbergii, Psilotum nudum, Spinacia oleracea, Triticum aestivum, Zea mays). Identity cutoff for protein coding genes was set at 60%. Identity cutoff for RNAs was set at 80%. The map of genome structure and gene distribution were carried out using OGDRAW V1.1 (OrganellarGenomeDRAW, http://ogdraw.mpimp-golm.mpg.de/), which takes a Genbank file or a special accession number (15).
3.3. Chloroplast Genome Analysis
Relative synonymous codon usage (RSCU) of different codons in each gene sample was calculated by codonW in Mobyle (http://mobyle.pasteur.fr/cgi-bin/portal.py). To determine the repeat sequence and location, an online version of REPuter (http://bibiserv. techfak.uni-bielefeld.de/reputer/) was used (16). Searching condition was followed as Saski (3).
4.1. Overall Structure
Chloroplast genome of G. thurberi (Figure 1) has a conserved quadripartite structure. Total genome is a circular DNA molecule of 160,264 bp, which is shorter than G. barbadnese (17) and G. hirsutum (18). The two single copy regions are separated by the two inverted repeats. The whole genome was analyzed (Table 1). The large single copy region is 88,737 bp, the small single copy is 20, 271 bp and the two inverted repeats are 25,628 bp each. The coding region is 91,485 bp in length, accounting for 57.08% of the whole plastidic genome, which is similar to Gossypium hirsutum by 56.46% (18), Bambusa oldhamii by 53.4% (13) and Dendrocalamus latiflorus by 53.4% (13), genus Megaleranthis 52.4% (5), genus Alsophila 53.2% (8). where as it is smaller than Glycine max (60%). G. thurberi plastidic genome codes for proteins (49.76%), tRNA genes (1.73%) and rRNA (5.60%), similar to Manihot esculenta (19), cucumber (20) and coffee (21). The non-coding region is 70,351 bp in length (43.90% of the genome). The proportions of intergenic spacers and intron are 31.15% and 12.75%, respectively.
4.2. Repeat Sequence
Chloroplast genome structures are similar to prokaryotes, it has been considered uncommon to have large scale of repeat sequences in these genomes. Here, PEPuter was used to detect the repeat sequence of cp genome of G. thurberi. Four types of repeats were detected; forward (direct) match, reverse match, complements match and palindromic (inverted) match. Sixty six repeats having more than 30 bp in length are listed in (Table 1). There are 2 complementary repeats, 28 forward repeats, 23 inverted repeats and 13 reverse repeats. Most of the repeats are located at ycf2 and intergenic spacers (IGS), and few located at trnS and introns. The largest repeat is 72 bp, which is located at psaB and psaA, while the most of the repeats are 30-40 bp. In addition to the four types of repeats, there are few simple sequence repeats (SSRs).
Simple sequence repeats were screened in G. thurberi chloroplast genome and 67 cpSSRs (³10 bp) were obtained. Most of the SSRs are mononucleotide repeats, while 16 dinucleotide repeats and 3 trinucleotide repeats. The longest repeat is the repeat of “CT”, which is 14 bp, but the most of the repeats are C and T having 13 bp (Table 2).
4.3. Gene Content and Codon Usage
Genes coded by the cp of G. thurberi are listed (Table 3). Among the total 79 protein coded genes, there are 4 rRNA genes, 30 tRNA genes and 113 single genes; out of which 20 genes are duplicated, locating at IR. According to the gene function, all genes can be classified as genes of the functional genetic system, the photosynthetic system, the biosynthesis and some with unknown function. In G. thurberi cp genome, five genes with unknown function (ycf gene) were detected and considered as to be essential in plants, which were highly conserved between species (22). Intrestingly, two genes, namely rps12 has an intron (5). rps12 was separated (by an intron) into two fragments with one exon locating at LSC (5¢-end) and the other at 3¢-end at IR. matK is 1.5 Kbp in length, and was found in the intron of trnK-UUU, which is the only gene located in an intron and encodes maturase K. This gene has both conserved and variable fragments (23). Thus, it is frequently used in phylogenetic studies (24, 25, 23, 14)
The codon usage was analyzed (Table 4). ATG and TGG code for methionine as the start codon and tryptophane, respectively with RSCU value=1. RSCU values of the three terminal codons TAA, TGA and TAG are 1.76, 0.64 and 0.6, respectively. Accoding to RSCU value, G. thurberi prefers TAA as its stop codon. The RSCU values grater than 1 indicates greater codon frequency. Most of the codons prefer A or T at the third position. The analysis of the composition for the codons showed that A+T content at the third position was 72.6%, similar to what was reported for Alophila (8) and Panax schinseng Nees (1).
5.1. SSR in cpDNA
Despite the high level of conservation in cp genome SSRs are evident as stated in previous reports (20) cpSSRs are useful in analysis of genetic diversity (26) due to their greater efficiency as opposed to genomic SSRs (26). Furthermore and due to the greater level of conservation, the information of the other species can be used to design specific primers for a species with unknown sequence data (27, 28, 14).
5.2. Gene Loss in Chloroplast
During the course of evolution, loss and gain of genetic material have been noted for cpDNA. For instance ycf15, a non-functional gene in other plants (10, 22, 12) is also present in G. thurberi. The other example is infA, most mobile gene between chloroplast and nuclear genome, that codes for a translation initial factor 1 (29, 30, 31). In our study similar to cassava (19) G. hirsutum (18), the infA was absent. However, some others had the infA as a pseudogene (17, 22), while in others infA appeared as an intact gene (21). Similar to G. hirsutum (18) and G. barbadense (17) and angiosperms, trnP-GGG was absent in G. thurberi. However, trnP-GGG was reported in Cryptomeria japonica (2). Thus it can be suggested that the gene has lost before the divergence of angiosperms. The other gene that worth considering was rpl22 that codes for the large subunit of ribosomal protein 22. rpl22 is present in G. thurberi chloroplast genome similar to G. barbadense (17), but has been reported to be absent in G. hirsutum (18) and 3 legumes, namely Glycine, Lotus and Medicago (3). Therefore, its analysis may shed some light on the evolution of Gossypium.
5.3. Extent of IR
The border of the IR is usually different between species and the IR expansion and contraction are important as far as genome size is concerned. IR expansion often leads to larger sizes of genome. Usually the pseudogenes are residing at the junction of IR and LSC/ SSC. The differences of the junctions among G. thurberi and five other species were analyzed (Figure 2).
The IRb/LSC junction was found within rps19 in G. thurberi, Solanum lycopersicum, Arabidopsis thaliana and Spinacia oleracea, indicating that rps19 was duplicated at the junction at IRa and LSC. This duplication is very common in plants (20, 13). G. thurberi (8 bp) and Spinacia (144 bp) have the shortest and the longest duplications, respectively.
On the border of IRb and SSC, G. thurberi is similar with G. barbadense, having 79 bp of ycf1 fragment on the IRb border. Solanum, Arabidopsis and Spinacia have the same type of overlapping of ycf1 and ndhF at the junction. The longest overlap is in Arabidopsis with 37 bp and the shortest is in Solanum with 17 bp. The overlap is also found in Cucumis (20). The ycf1 is located at the junction of SSC/IRb. So ycf1 was duplicated in IRb at the border of IRb and SSC. In the Spinacia, ycf1 has the longest duplication with 1445 bp. In G. thurberi and G. barbadense, ycf1 has the shortest duplication with 79 bp. Gossypium hirsutum is in opposite direction of SSC as compared to other five speciesand therefore; ycf1 is located at the junction of IRb/SSC with 98 bp.
In the G. thurberi, 18 genes were found containing one or two introns, which is the same as in Panax schinseng Nees (1). In contrast to G. thurberi, introns were absent in rpoC1 and clpP in B. Oldhamii and D. Latiflorus (13). The number and location of the intron in chloroplast seems to be conserved. The comparison (Table 5) of introns among G. thurberi, G. hirsutum and G. barbadense shows that 18 genes have one or two introns in the cp genome; 6 of which are tRNA coding genes and the rest are protein-coding genes. The longest intron is located in trnK-UUU with 2542 bp in G. thurberi, which is the only intron in cotton with another gene, matK, inside. The smallest intron is located in rpl12-3end with 536 bp, which is situated in IR. Genes of ycf3 and clpP are located at LSC and are divided by two introns. Small variations among the introns were noted in three cotton species; intron lengths for ycf3-1, ycf3-2, rps12_3end, ndhB and ndhA are conserved, while the others have small variations. These intron sequences have a high identity, especially ycf3-2 with a 100% sequence identity among the three cotton species.
5.5. GC Content
The GC content of G. thurberi cp genome is 37.22 %, similar to other plants, such as 37.86% in Solanum lycopersicum, 37.85% in Nicotiana tabacum, 37.56% in Atropa belladonna, 37.25% in G. hirsutum and 34% in Glycine max. Both coding and non-coding regions are low in GC content (32) with 40.35% and 33.08%, respectively in G. thurberi. Variation in GC content among four different regions in G. thurberi cpDNA was observed (Table 6) and IR was the richest (42.95%), similar to an earlier report (9). It is supposed that ribosomal genes (rrna4.5, rrna5, rrna16, rrna23) and coding regions (19, 8 ,10, 5) are responsible for high GC content in IR. GC content distribution of the each region is similar with other species (20, 10, 5). According to Gao (8) GC content was uneven across cp genome in Alsophila. In this study, we cut genomes of the three Gossypium species into 1 kb-unit to compare uint to uint GC content across whole chloroplast genome. Our result showd that the distribution of GC is similar across the whole genome (Figure 3). At SSC region G. barbadense and G. thurberi are similar but different from G. hirsutum because of different direction of SSC. Across the whole genomes of the three cotton species, different fragments and even the adjacent fragments share different GC contents.
Gao (8) reported that GC contents in the chloroplast genomes are not the same between genes in different functional groups; rRNA (55.18%)>tRNA (54.55%)>photosynthetic (43.85%)>genetic system (40.80%)>NADH (39.54%). In G. thurberi, similar data was obtained. In the coding region, the rRNA genes have the highest GC content (55.51%) and the protein genes have the lowest (37.69%). In the non-coding region, GC content of IGS and intron is 31.58% and 36.75%, respectively. The non-coding regions experienced a fast evolution, thus the non-coding region is richer in GC than coding regions.
GC content is an important feature of a genome that is correlated to the number of microRNA binding sites (33), functional elements physical location (34), recombination rate and gene distribution (35), organelle RNA editing (36) and gene expression regulation (37). GC content varies in the 5¢UTR and 3¢UTR (34). GC content has a rare relationship with replication timing in human genome (38). GC content also have an effect on RNAi, because it is highly correlated to RNAi target site accessibility and negatively correlated with RNAi activity (39).
Low GC content is a significant feature of plastid genomes, which is possibly formed after endosymbiosis by DNA replication and repair (32). In viruses GC content is not dependent on genes constitution, but it is correlated with its location (40). Whether or not this exists in chloroplast genome needs more efferts on further studies.