Chloroplast DNA codon use: Evidence for selection at the psb A locus based on tRNA availability

Codon use in the three sequenced chloroplast genomes (Marchantia, Oryza, and Nicotiana) is examined. The chloroplast has a bias in that codons NNA and NNT are favored over synonymous NNC and NNG codons. This appears to be a consequence of an overall high A + T content of the genome. This pattern of codon use is not followed by the psb A gene of all three genomes and other psb A sequences examined. In this gene, the codon use favors NNC over NNT for twofold degenerate amino acids. In each case the only tRNA coded by the genome is complementary to the NNC codon. This codon use is similar to the codon use by chloroplast genes examined from Chlamydomonas reinhardtii. Since psb A is the major translation product of the chloroplast, this suggests that selection is acting on the codon use of this gene to adapt codons to tRNA availability, as previously suggested for unicellular organisms.     

Comparative studies on codon usage pattern of chloroplasts and their host nuclear genes in four plant species

A detailed comparison was made of codon usage of chloroplast genes with their host (nuclear) genes in the four angiosperm speciesOryza sativa, Zea mays, Triticum aestivum andArabidopsis thaliana. The average GC content of the entire genes, and at the three codon positions individually, was higher in nuclear than in chloroplast genes, suggesting different genomic organization and mutation pressures in nuclear and chloroplast genes. The results of Nc-plots and neutrality plots suggested that nucleotide compositional constraint had a large contribution to codon usage bias of nuclear genes inO. sativa, Z. mays, andT. aestivum, whereas natural selection was likely to be playing a large role in codon usage bias in chloroplast genomes. Correspondence analysis and chi-test showed that regardless of the genomic environment (species) of the host, the codon usage pattern of chloroplast genes differed from nuclear genes of their host species by their AU-richness. All the chloroplast genomes have predominantly A- and/or U-ending codons, whereas nuclear genomes have G-, C- or U-ending codons as their optimal codons. These findings suggest that the chloroplast genome might display particular characteristics of codon usage that are different from its host nuclear genome. However, one feature common to both chloroplast and nuclear genomes in this study was that pyrimidines were found more frequently than purines at the synonymous codon position of optimal codons.     

巨枝叶绿体基因组密码子偏好性分析

该文针对巨桉叶绿体基因组序列,选取其中长于300 nt且以AUG为起始密码子的43个非重复基因作为研究对象,采用CodonW1.4.2软件分析巨桉叶绿体基因组的密码子使用偏好性。结果表明:第3位密码子的平均GC含量为27.97%; ENC的变化范围为39.49～61.00,平均为47.04; RSCU1的密码子有31个,其中29个以A/U结尾;中性分析显示,GC12与GC3无显著相关;回归分析未达到显著性水平; ENCplot分析发现,大部分基因落在曲线上或附近;对应分析表明第1轴的贡献率为17.68%,第2轴的贡献率为11.49%,第3轴、第4轴的贡献率分别为8.00%和5.76%,前4轴累计贡献率达42.93%,第1轴与GC、ENC、CAI达到极显著相关。上述分析结果表明,巨桉叶绿体基因组的密码子偏好较弱,密码子第3位偏好以A或U结尾,选择和突变在巨桉叶绿体基因组密码子偏好中起相对均衡的作用,最终确定UUG、CUU、GUU、UCC、UCA、ACA、UAU、UAA、CAU、AAU、AGA和GGA 12个高频高表达密码子为最优密码子。这为转化叶绿体基因密码子优化,提高表达效率和改良巨桉目标性状奠定了坚实基础。     

越南金花茶叶绿体基因组分析

本研究通过高通量测序技术对越南金花茶(Camellia insularis Orel & Curry)的叶绿体(chloroplast, cp)基因组进行了测序。结果显示,越南金花茶叶绿体基因组长度为156 882 bp,包含132个基因,其中,88个编码蛋白质,36个编码tRNA,8个编码rRNA。密码子偏好性分析显示,编码亮氨酸的密码子数量最多,且第三位密码子显示出了较高的A/U偏好。此外,共鉴定出了67个简单重复序列(simple sequence repeats, SSR),发现越南金花茶SSR偏好使用A和T碱基。除了少数可变区域之外,越南金花茶与其近缘物种的叶绿体基因组反向重复(inverted repeat, IR)边界区域相当保守。系统发育分析表明,越南金花茶与云南金花茶(Camellia fascicularis)亲缘关系最近。金花茶是基因工程培育黄色山茶花的重要材料,本研究为叶绿体工程提供了基本的遗传信息,为深入研究金花茶组植物的进化、物种鉴定和基因组育种研究提供了宝贵的资源。     

两种梧桐叶绿体基因组密码子使用偏性分析

为了分析美丽梧桐、云南梧桐叶绿体基因组密码子的使用偏性,该研究通过筛选美丽梧桐、云南梧桐叶绿体基因组中各52条蛋白编码序列,并利用CodonW、CUSP和SPSS软件对其密码子使用模式及偏性进行了分析。结果表明:(1)美丽梧桐、云南梧桐的GC含量分别为38.12%、38.05%,表明叶绿体基因组内富含A/T碱基。(2)有效密码子数(ENC)范围为36.91～56.46、36.55～58.04,表明多数密码子偏性较弱。(3)相对同义密码子(RSCU)分析显示,RSCU1的密码子各有29个,其中28个以A、U结尾。(4)中性绘图显示,GC_3与GC_(12)的相关性不显著,回归曲线斜率分别为0.195和0.304,说明密码子偏好性主要受到自然选择的影响。(5) ENC-plot分析中大部分基因分布于曲线的周围和下方,ENC比值多分布于-0.04～0.10之间,表明突变会影响密码子偏性的形成。此外,17、18个密码子分别被鉴定为美丽梧桐、云南梧桐的最优密码子。以上结果说明美丽梧桐、云南梧桐叶绿体基因组的密码子使用偏性可能受选择和突变共同作用,且使用模式较为相似,但具有一定的差异,可能与适应环境的进化机制有关。     

红锥叶绿体基因组特征及系统发育分析

为确定红锥(Castanopsis hystrix)叶绿体基因组的结构组成情况,判定其在锥属中的进化位置及与同锥属叶绿体基因组的区别,为锥属物种鉴定、遗传多样性分析和资源保护提供相关依据。使用Illumina HiSeq 2500测序平台对红锥叶绿体基因组进行测序,通过生物信息学分析方法进行序列组装、注释和特征分析,并利用R、Python、MISA、CodonW和MEGA 6等生物信息学软件对其基因组结构和数目、密码子偏好性、序列重复、简单重复序列(simple sequence repeat,SSR)位点和系统发育进行分析。结果表明红锥叶绿体基因组大小为153754 bp,呈现四分体结构;共拥有130个基因,包含85个编码基因、37个tRNA基因和8个rRNA基因;通过密码子偏好性分析,平均有效密码子数为55.5,说明其密码子随机性强、偏好性低;通过SSR及长重复片段分析,检测到45个重复序列及111个SSR位点;与近缘种比较,发现其叶绿体基因组序列高度保守,尤其蛋白质编码序列相似度极高;此外,系统发育分析发现红锥与海南锥聚为一支,关系密切。本研究得到了红锥的叶绿体基因组基本情况与系统发育位置,为红锥的物种辨别、天然种群遗传多样性与功能基因组学提供前期研究铺垫。     

Adjustment of the tRNA population to the codon usage in chloroplasts.

In chloroplasts there is a correlation between the amounts of tRNAs specific for a given amino acid and the codons specifying this amino acid. Furthermore, for the amino acids coded for by more than one codon, the population of isoaccepting tRNAs is adjusted to the frequency of synonymous codons used in chloroplast protein genes. A comparison by two-dimensional gel electrophoresis of the tRNA populations extracted from chloroplasts and from chloroplast polysomes shows that all chloroplast tRNAs are involved in protein biosynthesis.     

紫花苜蓿叶绿体基因组密码子偏好性分析

为分析紫花苜蓿叶绿体基因组密码子偏好性的使用模式,该文以紫花苜蓿叶绿体基因组中筛选到的49条蛋白质编码序列为研究对象,利用CodonW、CUSP、CHIPS、SPSS等软件对其密码子的使用模式和偏好性进行研究。结果表明:(1)紫花苜蓿叶绿体基因的第3位密码子的平均GC含量为26.44%,有效密码子数(ENC)在40.6～51.41之间,多数密码子的偏好性较弱。(2)相对同义密码子使用度(RSCU)分析发现,RSCU>1 的密码子数目有30个,以A、U结尾的有29个,说明了紫花苜蓿叶绿体基因组A或U出现的频率较高。(3)中性分析发现,GC3与 GC12的相关性不显著,表明密码子偏性主要受自然选择的影响; ENC-plot 分析发现一部分基因落在曲线的下方及周围,表明突变也影响了部分密码子偏性的形成。此外,有17个密码子被鉴定为紫花苜蓿叶绿体基因组的最优密码子。紫花苜蓿叶绿体基因组的密码子偏好性可能受自然选择和突变的共同作用。该研究将为紫花苜蓿叶绿体基因工程的开展和目标性状的遗传改良奠定基础。     

滇黄精叶绿体全基因组序列及其密码子使用偏性分析

为探究滇黄精(Polygonatum kingianum)叶绿体全基因组特征和密码子使用偏性,利用第二代测序技术对滇黄精嫩叶进行测序,再经组装与注释后得到其叶绿体基因组全序列,通过MISA、EMBOSS和CodonW等软件对滇黄精叶绿体全基因组的SSR位点、系统发育及密码子偏好性进行分析。结果表明,滇黄精完整叶绿体基因组长度为155 852 bp,基因组平均GC含量为37.7%,其大、小单拷贝区(LSC)长度分别为84 633和185 25 bp,反向重复区长度为26 347 bp,注释了132个基因,包括86个蛋白编码基因、38个tRNA基因和8个核糖rRNA基因。叶绿体基因组中共有69个SSR位点,绝大多数属于单碱基重复的A/T类型。系统发育分析表明滇黄精与格脉黄精(P. tessellatum)亲缘关系近,可能与分布地域有关。密码子偏好性分析表明,滇黄精叶绿体基因组密码子使用模式受到自然选择影响大于突变因素,最终确定9个最优密码子。因此, 滇黄精叶绿体基因组遗传结构和系统发育位置及其密码子偏倚的分析,为叶绿体基因工程研究提供理论依据。     

Synonymous codon usage in chloroplast genome of Coffea arabica

Synonymous codon usage of 53 protein coding genes in chloroplast genome of Coffea arabica was analyzed for the first time to find out the possible factors contributing codon bias. All preferred synonymous codons were found to use A/T ending codons as chloroplast genomes are rich in AT. No difference in preference for preferred codons was observed in any of the two strands, viz., leading and lagging strands. Complex correlations between total base compositions (A, T, G, C, GC) and silent base contents (A₃, T₃, G₃, C₃, GC₃) revealed that compositional constraints played crucial role in shaping the codon usage pattern of C. arabica chloroplast genome. ENC Vs GC₃ plot grouped majority of the analyzed genes on or just below the left side of the expected GC₃ curve indicating the influence of base compositional constraints in regulating codon usage. But some of the genes lie distantly below the continuous curve confirmed the influence of some other factors on the codon usage across those genes. Influence of compositional constraints was further confirmed by correspondence analysis as axis 1 and 3 had significant correlations with silent base contents. Correlation of ENC with axis 1, 4 and CAI with 1, 2 prognosticated the minor influence of selection in nature but exact separation of highly and lowly expressed genes could not be seen. From the present study, we concluded that mutational pressure combined with weak selection influenced the pattern of synonymous codon usage across the genes in the chloroplast genomes of C. arabica.     

沙枣叶绿体全基因组序列及其使用密码子偏性分析

该研究以2株野生沙枣(Elaeagnus angustifolia Linn.)嫩枝经温室水培后的嫩叶为材料,采用CTAB法分别提取总DNA,并利用第二代测序技术进行总DNA从头测序,组装后得到2株沙枣叶绿体基因组全序列,并详细分析了其蛋白质编码基因密码子使用的偏好性及其原因,为沙枣叶绿体基因工程和分子系统进化等研究奠定基础。结果显示:(1)组装得到沙枣叶绿体基因组序列全长150 546 bp,由长度为81 113 bp的长单拷贝(LSC)区域和25 494 bp的短单拷贝(SSC)区域,以及1对分隔开它们的长18 445 bp的反向重复序列(IRS)组成;注释共得到132个基因,包括86个蛋白编码基因、38个tRNA基因和8个rRNA基因。(2)沙枣叶绿体基因组蛋白编码基因密码子的第三位碱基GC含量(GC_3)为28.47%,明显低于整个叶绿体基因组GC含量(37%),也低于第一位(GC_1)和第二位(GC_2)碱基的GC含量,说明密码子对AT碱基结尾有偏好性;其中, UCU、CCU、UGU、GCU、CUU、GAU、UCA和UAA为最优密码子。(3)同义密码子相对使用频率(RSCU)分析发现,影响密码子使用模式的因素并不单一,密码子的偏好性受到突变、选择及其他因素的共同影响,并且自然选择表达引起的序列差异比突变对密码子偏好性的影响要显著;中性绘图分析、有效密码子数(ENC-plot)分析和奇偶偏好性(PR2-plot)分析表明,沙枣叶绿体基因组使用密码子的偏性受选择的影响更大。(4)通过最大似然法、最大简约法和贝叶斯方法对胡颓子科6个物种和1个枣的叶绿体基因序列构建系统发育树,与它们使用密码子偏性聚类的结果一致,表明叶绿体基因组使用密码子偏性与物种的亲缘关系相关。     

紫九牛叶绿体基因组密码子偏好性分析

为确定瑶药紫九牛叶绿体基因组密码子的使用模式及其成因,该研究以紫九牛叶绿体基因组50条蛋白质编码序列为研究对象,利用Codon W 1.4.2和在线软件CUSP和Chips分析其密码子偏好性。结果表明:(1)RSCU>1的密码子有29个,其中有28个以A/U结尾,说明叶绿体基因组的同义密码子中偏好以A/U结尾。(2)紫九牛叶绿体基因组密码子的GC含量GC₁(47.38%)>GC₂(39.81%)>GC₃(29.60%),ENC值大于45的有40个,说明紫九牛叶绿体基因组存在较弱的偏性。(3)中性绘图分析和ENC-plot分析说明了紫九牛叶绿体基因组密码子的偏好性既受到选择的作用,又受到突变因素的影响。(4)通过构建的高低基因表达库最终确定了15个最优密码子,分别为UUG、AUU、GUU、GUA、UCU、 CCU、ACU、ACA、GCU、CAA、AAC、GAA、UGU、CGU和GGU。该研究为紫九牛叶绿体基因组的确定以及遗传多样性分析提供了依据。     

多花海棠叶绿体基因组分析

多花海棠（Malus floribunda Siebold.）是世界范围内广泛栽培的苹果属物种,具有较高的观赏价值和育种意义。对其进行叶绿体基因组比较分析,有利于完善苹果属系统进化以及种质利用的研究内容。基于全基因组测序数据,组装获得一个完整的具有四分体结构的多花海棠叶绿体基因组。该基因组包括大单拷贝区（88 142 bp）、反向重复区B （26 353 bp）、小单拷贝区（19 189 bp）与反向重复区A （26 353 bp）,共计160 037 bp。多花海棠叶绿体全基因组共注释到111个基因,包括78个蛋白编码基因、29个tRNA基因和4个rRNA基因。此外,在其基因组中识别到大量的重复序列,与三叶海棠和变叶海棠略有差异。通过计算相对同义密码子使用度,发现其高频密码子共30种,并且密码子具有偏向A/T结尾的使用模式。种间序列比对、边界分析的结果表明,大单拷贝区序列变异较大,8种苹果属植物SC区与IR区扩张收缩情况整体上较为相似。基于叶绿体基因组序列的系统进化分析,将多花海棠、湖北海棠和变叶海棠聚为一类。多花海棠叶绿体基因组的研究可为今后遗传标记开发与种质资源利用等提供数据支持。     

The close proximity ofEscherichia coli genes: Consequences for stop codon and synonymous codon use

It is shown that synonymous codon usage is less biased in favor of those codons preferred by highly expressed genes at the end ofEscherichia coli genes than in the middle. This appears to be due to the close proximity of manyE. coli genes. It is shown that a substantial number of genes overlap either the Shine-Dalgarno sequence or the coding sequence of the next gene on the chromosome and that the codons that overlap have lower synonymous codon bias than those which do not. It is also shown that there is an increase in the frequency of A-ending codons, and a decrease in the frequency of G-ending codons at the end ofE. coli genes that lie close to another gene. It is suggested that these trends in composition could be associated with selection against the formation of mRNA secondary structure near the start of the next gene on the chromosome. Stop codon use is also affected by the close proximity of genes; many genes are forced to use TGA and TAG stop codons because they terminate either within the Shine-Dalgarno or coding sequence of the next gene on the chromosome. The implications these results have for the evolution of synonymous codon use are discussed.     

樟树叶绿体基因组密码子偏好性分析

为分析樟树(Cinnamomum camphora)叶绿体基因组密码子偏好性使用模式,该研究利用CodonW、EMBOSS、R语言等软件和程序,对53条樟树叶绿体基因组密码子使用模式及偏好性进行了系统分析。结果表明:樟树叶绿体基因的有效密码子数(ENC)在36.82～59.30之间,表明密码子的偏好性较弱。相对同义密码子使用度(RSCU)分析发现RSCU>1的密码子有32个,其中28个以A、U结尾,表明第3位密码子偏好使用A和U碱基。中性绘图分析发现GC₃与GC₁₂的相关性不显著,回归曲线斜率为0.049,说明密码子偏好性主要受到自然选择的影响。ENC-plot分析发现大部分基因落在曲线的下方,同样表明选择是影响密码子偏好性的主要因素。该研究发现共有9个密码子(UUU、CUU、UCA、ACA、UAU、AAU、GAU、UGA、GGA)被鉴定为樟树叶绿体基因组的最优密码子。     

槲蕨属叶绿体基因组密码子偏好性分析

叶绿体基因组密码子偏好性使用模式往往影响基因表达效率,为促进药用植物叶绿体基因工程的发展,提高槲蕨药用品质,该研究以川滇槲蕨、栎叶槲蕨和槲蕨三个近缘药用植物为材料,使用CodonW、CUSP和SPSS等软件分析其叶绿体基因组编码基因密码子使用偏好性,筛选出三个物种的最优密码子。结果表明:川滇槲蕨、栎叶槲蕨和槲蕨叶绿体基因组的有效密码子数(ENC)范围分别为40.10~61、40.33~61和40.15~61,其密码子偏好性较弱;ndhE、rpl22、rpl14、rpl20、ccsA、rps4和rpl16编码基因的ENC值差异较大,表明近缘物种中,部分基因的密码子偏好性存在一定差异;三个物种编码基因的ENC频数集中于-0.1~0.1之间,说明槲蕨属基因密码子偏好性主要受到突变的影响。川滇槲蕨12个最优密码子有6个和栎叶槲蕨相同,分别是UCU、ACU、GCU、CAA、AAA和GAU;栎叶槲蕨10个最优密码子有2个和槲蕨相同,分别是UUA和AUU,而川滇槲蕨与槲蕨无相同的最优密码子。该研究结果可为槲蕨属药用植物基因工程中外源基因的改良及其表达奠定基础。     

Evidence from simulation studies for selective constraints on the codon usage of the Angiosperm psbA gene

The codon usage of the Angiosperm psbA gene is atypical for flowering plant chloroplast genes but similar to the codon usage observed in highly expressed plastid genes from some other Plantae, particularly Chlorobionta, lineages. The pattern of codon bias in these genes is suggestive of selection for a set of translationally optimal codons but the degree of bias towards these optimal codons is much weaker in the flowering plant psbA gene than in high expression plastid genes from lineages such as certain green algal groups. Two scenarios have been proposed to explain these observations. One is that the flowering plant psbA gene is currently under weak selective constraints for translation efficiency, the other is that there are no current selective constraints and we are observing the remnants of an ancestral codon adaptation that is decaying under mutational pressure. We test these two models using simulations studies that incorporate the context-dependent mutational properties of plant chloroplast DNA. We first reconstruct ancestral sequences and then simulate their evolution in the absence of selection on codon usage by using mutation dynamics estimated from intergenic regions. The results show that psbA has a significantly higher level of codon adaptation than expected while other chloroplast genes are within the range predicted by the simulations. These results suggest that there have been selective constraints on the codon usage of the flowering plant psbA gene during Angiosperm evolution.     

蔓赤车叶绿体基因组结构及系统进化分析

蔓赤车(Pellionia scabra)属于荨麻科赤车属,是一种具有高营养价值的优质野菜。本研究以蔓赤车为材料,基于高通量技术方法完成叶绿体基因组测序、组装注释、结构解析及构建系统发育树,以此深入研究蔓赤车叶绿体基因组特征。结果表明,蔓赤车叶绿体基因组大小为153220 bp,GC含量为36.4%,属于典型的四分体结构。共注释到130个基因,其中85个蛋白编码基因,37个转运RNA基因,8个核糖体RNA基因;其中,有15个基因包含1个内含子,2个基因包含2个内含子,rps12存在反式剪接情况。蔓赤车叶绿体基因组可分为光合作用(43个)、自我复制(64个)、其他编码蛋白(7个)以及未知功能(4个)4大类基因。蔓赤车叶绿体基因组共检测出51073个密码子,其中编码亮氨酸(Leu)的密码子占比最大,密码子偏向使用A和U两种碱基。检测到72个简单重复序列位点,其中有58个单核苷酸、12个二核苷酸、1个三核苷酸和1个四核苷酸4种不同类型的简单重复序列(simple sequence repeats,SSRs)。蔓赤车IRb/SSC边界存在ycf1基因扩张现象。系统进化树显示,蔓赤车(Pellionia scabra OL800583)与庐山楼梯草(Elatostema stewardii MZ292972)、盘托楼梯草(Elatostema dissectum MK227819)、光叶楼梯草(Elatostema laevissimum var.laevissimum MN189961)亲缘关系最密切。基于蔓赤车叶绿体基因组的研究,旨在为蔓赤车物种鉴定、遗传进化及基因组学研究奠定理论基础。     

Ribosome-binding sites on chloroplastrbcL andpsbA mRNAs and light-induced initiation of D1 translation

Chloroplast ribosome-binding sites were identified on the plastidrbcL andpsbA mRNAs using toeprint analysis. TherbcL translation initiation domain is highly conserved and contains a prokaryotic Shine-Dalgarno (SD) sequence (GGAGG) located 4 to 12 nucleotides upstream of the initiator AUG. Toeprint analysis ofrbcL mRNA associated with plastid polysomes revealed strong toeprint signals 15 nucleotides downstream from the AUG indicating ribosome binding at the translation initiation site.Escherichia coli 30S ribosomes generated similar toeprint signals when mixed withrbcL mRNA in the presence of initiator tRNA. These results indicate that plastid SD sequences are functional in chloroplast translation initiation. ThepsbA initiator region lacks a SD sequence within 12 nucleotides of the initiator AUG. However, toeprint analysis of soluble and membrane polysome-associatedpsbA mRNA revealed ribosomes bound to the initiator region.E. coli 30S ribosomes did not associate with thepsbA translation initiation region.E. coli and chloroplast ribosomes bind to an upstream region which contains a conserved SD-like sequence. Therefore, translation initiation onpsbA mRNA may involve the transient binding of chloroplast ribosomes to this upstream SD-like sequence followed by scanning to localize the initiator AUG. Illumination 8-day-old dark-grown barley seedlings caused an increase in polysome-associatedpsbA mRNA and the abundance of initiation complexes bound topsbA mRNA. These results demonstrate that light modulates D1 translation initiation in plastids of older dark-grown barley seedlings.     

极小种群植物川柿叶绿体基因组特征与系统发育分析（附录）

川柿(Diospyros sutchuensis)为极小种群和国家重点保护野生植物,分布范围狭窄,种群数量极少。目前,川柿基因组信息缺乏,在柿属(Diospyros)中的系统亲缘关系不明确。该研究通过Illumina平台对川柿叶绿体基因组进行测序,应用Getorganellev1.7.3.4和PGA软件对基因组进行组装和注释,使用DnaSP6.12.03软件进行多序列对比分析,并使用REPuter、Tandem Reapeats Finder和MISA软件进行重复序列分析,使用CodonW1.4和EasyCodemL软件分别进行密码子偏好性和选择压力分析。同时,基于4个不同的叶绿体基因组序列数据集,使用IQtree软件分析川柿与11个柿属物种的系统发育关系。结果表明:(1)川柿叶绿体基因组全长157 917 bp,包含1对26 111 bp的反向重复区、大单拷贝区(87 303 bp)和小单拷贝区(18 392 bp),GC碱基含量为37.4%。(2)川柿叶绿体基因组共注释到113个基因,包括79个蛋白编码基因、30个tRNA基因和4个rRNA基因; 共检测到49个长重复序列、27个串联重复序列和34个简单重复序列; 蛋白编码基因中高频密码子31个,多数密码子末位碱基为A或U,编码亮氨酸的密码子使用最多; 基因组编码区比非编码区更为保守,10个高变热点区域可作为潜在的分子标记; 蛋白编码基因中有8个基因(ndhB、ndhG、ndhI、rbcL、rpoB、petB、petD、rps12)受到正选择压力。(3)系统发育分析显示,川柿与老鸦柿(D. rhombifolia)和乌柿(D. cathayensis)亲缘关系最为密切,它们与海南柿(D. hainanensis)共同形成一个单系分支。该研究结果既为川柿及柿属种质资源鉴定、遗传多样性保护以及种群恢复等提供了叶绿体基因组资源,也为阐明川柿的系统进化提供了重要的分子信息。