湿地松转录组SSR分析及EST-SSR标记开发

【目的】湿地松是南方地区优先推广的优质产脂树种。但湿地松分子遗传基础薄弱,基因组序列信息匮乏,影响了湿地松基因组学的深入研究。目前,湿地松分子研究所用的SSR标记主要来自其他近缘种或利用公共数据库中有限的基因序列资源开发的SSR标记,其多态性和通用性较差。为解决这一问题, 笔者根据湿地松转录组测序数据开发EST-SSR位点,并揭示其在转录组序列中的分布类型及特征,为湿地松分子标记辅助育种奠定基础。【方法】利用MISA软件对转录组序列进行SSRs查找和分布特征分析。查找标准参数设置为: 单核苷酸重复>10次,二核苷酸重复>6次,三、四、五、六核苷酸重复>5次。根据SSR位点两端的保守区域,利用Primer3.0设计并随机挑选120对SSR引物,通过琼脂糖电泳和毛细管电泳对来自美国和吉安的113份家系个体进行遗传多样性分析,确定引物多态性。【结果】79 574条unigenes序列中搜索到3 818个SSR位点,出现频率为4.80%,平均18.27 kb出现1个SSR位点,3 373个unigenes含有SSR位点,SSR发生频率(含SSR位点的序列数/搜索序列总数)为4.24%,其中2 980条序列含1个SSR位点,含1个以上SSR位点的序列有393条。在检测到的3 818个SSR标记中,单核苷酸分布最多,其次是二核苷酸和三核苷酸,SSR数量分别占总数的63.54%、19.15%和16.27%,而四、五、六核苷酸重复类型所占比例较小,分别为0.52%、0.13%和0.31%。SSR重复单元的重复次数分布在5~22次之间,除单核苷酸重复外的1 391个SSR中,重复5次的数量最多,为498个(35.80%);重复6次和7次的次之,分别为417个(29.98%)和198个(14.23%);重复10次以上的仅有38个(2.73%)。在检测到的731个二核苷酸重复SSR中,最常出现的重复单元为AT/AT,数量为491个(12.86%),AG/CT和AC/GT两种类型的重复单元出现的次数次之,分别为156(4.09%)和81个(2.12%)。在检测到的621个三核苷酸重复中,AAT/ATT是出现频率最高的单元,共139个(3.64%),其次是AAG/CTT,共122个(3.20%)。3 818个SSR中有24.59%的位置未知,其余的SSR则分布在非编码区域(untranslated region,UTR)或者编码区(coding sequence,CDS)上,分布数量表现为3'UTR>5'UTR>CDS。参试的120对SSR引物,有92对扩增成功(76.78%),其中24对呈现多态(20%)。24对引物(13个二核苷酸重复、7个三核苷酸重复和4个四核苷酸重复)共检测出81个等位基因,每个位点的等位基因数为2~9,平均为3.38个。多态性信息含量(polymorphism information content,PIC)为0.103~0.726,平均为0.349。【结论】通过对湿地松转录组数据的挖掘,共获得 3 818个SSR位点,主要重复单元为AT/AT和AAT/ATT,可扩增出多态性位点的引物重复单元以二、三核苷酸重复为主。基于湿地松转录组序列的SSR标记开发是可行的。     

青钱柳叶片转录组数据组装及基因功能注释

青钱柳是一种民间常用中药,青钱柳叶中有多种生物活性的次生代谢物,如有机酸、黄酮类、皂苷类、铱类、精油、无机元素等,但对于青钱柳次生代谢产物合成的分子机制仍未见有报道.使用Illumina公司Hiseq 4000平台对青钱柳叶进行转录组测序,对 reads 进行拼接,得到50 126个unigenes平均长度为 1 247 bp,有19 875 (39.65%) unigenes由NCBI非冗余数据库进行注释的,23 716 (47.31%) unigenes被注释到GO数据库,14 950个(29.82%) unigenes匹配到COG 功能组.进一步的分析结果显示,6 012 (11.20%)个 unigenes被富集到 254个 KEGG 代谢通路,其中1 212个unigenes参与了次生代谢产物的生物合成,并且发现126个unigenes参与异戊二烯代谢,包括萜烯类、香茅醛、呋喃酮、苷的代谢途径.此外,总共检测到21 089个简单序列重复(SSRs).通过对老叶与嫩叶的转录组比较分析结果显示,在青钱柳不同生长发育时期的叶片中,基因表达上调占主异作用的过程主要涉及萜类和多肽的代谢、辅酶和维生素的代谢和氨基酸代谢,而基因表达下调占主异作用的过程主要涉及信号传输、类脂物代谢与能量代谢.该转录组数据可为青钱柳重要生物活性成分的生物合成和调控提供参考.     

Characterization, development and exploitation of EST-derived microsatellites in Gossypium raimondii Ulbrich

COTTON IS AN IMPORTANT GLOBAL CASH CROP. IN THE RECENTYEARS, MOLECULAR MARKER TECHNOLOGY HAS BEEN WIDELY APPLIED TO SUCH STUDIES ON COTTON AS GENETIC MAP- PING[1―4], VARIETY PURITY DETECTION[5], GENETIC DIVERSITY ANALYSIS[6], MOLECULAR MARKER-ASSISTED BR…     

松属编码区微卫星特征和相应基因功能分析

从GenBank下载453 892条松属EST序列,序列组装后得到20 886条contig。用Sputnik软件从这些contig中查找了2 678个微卫星,其中3碱基重复微卫星占的比例最高,为59.2%,而其他重复长度的微卫星都相对较低,比例分别为:2碱基重复微卫星占12.0%,4碱基重复微卫星占13.3%,5碱基重复微卫星占15.5%。3碱基重复微卫星变化引起的基因读码框改变最小,松属树种基因区3碱基重复微卫星的富集显示了强烈的密码子选择效应。此次研究还对查找到的微卫星进行了引物设计和扩增分析。实验结果显示,设计的微卫星引物在云南松中的扩增成功率是72.9%。从扩增成功的引物中进一步选取了155对引物,对14个松属树种和1个黄杉属树种进行了引物通用性实验分析,结果显示155对引物在14个松属树种间的通用性在71.0%以上,而在黄杉属树种中的通用性只有25.2%。对松属树种中含有微卫星的基因进行了功能分类研究,结果显示基因在是否保留微卫星序列方面有显著分化,微卫星参与了如细胞成分分类的共质体组成、病毒颗粒及病毒颗粒组成、生物节律调控, 以及生长素转运蛋白等生物学过程。     

油茶基因组微卫星特征分析

对油茶基因组约10%覆盖度的DNA序列进行微卫星查找,共获得11 344个重复单元长度为1～6碱基的微卫星。在此基础上,通过对这些微卫星序列分析发现:在油茶基因组中长度为二核苷酸的微卫星重复单元最为丰富,占27.1%;在单碱基重复和二碱基重复这两种类型中,最主要的优势重复单元分别是A/T以及AT/TA、AG/TC。三碱基、四碱基、五碱基重复类型中,(AAN)n、(AAAN)n和(AAAAN)n为对应的优势重复单元,这些优势重复单元中富含碱基A和T。油茶基因组中变异程度高的微卫星(长度≥20 bp)约占11.7%。分析还发现,除单核苷酸重复微卫星外,油茶基因组微卫星长度的变异速率与重复单元长度呈负相关,即油茶基因组中长度较短的微卫星变异速率较快,而较长的重复单元变异速度较慢,相对较为稳定。     

Comparative expression analysis of three genes from the Arabidopsis vacuolar Na^＋/H^＋ antiporter （AtNHX） family in relation to abiotic stresses

Na~ /H~ antiporters (NHX) are ubiquitous transmembrane proteins that play a key role in salt tolerance of plants. In this study, the sequence of 3 Arabidopsis NHX gene (AtNHX2―4) were compared with other AtNHX members. Putative cis-elements analysis identified elements that have been associated with stress responses. The activities of the promoters AtNHX2―4 were studied in transgenic plants carrying corresponding promoter-β-glucuronidase (GUS) fusions. The AtNHX2 promoter-GUS analysis indicated that AtNHX2 was expressed in constitutive pattern with high GUS activity in roots and leaves. AtNHX2 promoter activity was not up-regulated by NaCl or abscisic acid (ABA), in contrast to the AtNHX1 promoter which was previously studied. The AtNHX3 and AtNHX4 promoters showed tissue-specific activities. Strong GUS activity was detected in roots and vascular bundles of the stele in plants carry-ing an AtNHX4 promoter-GUS fusion, and GUS activity increased under salt stress suggesting a func-tion related to salt tolerance. Transgenic plants carrying the AtNHX3 promoter-GUS fusion showed strong GUS activity in petals, stamens and tops of siliques, suggesting a possible role of AtNHX3 in flower and seed development. Results of histochemical analysis suggested that AtNHX2―4 are involved in divergent functions and are differentially regulated under abiotic stress. The structure of AtNHX4 was predicted to include 12 transmembrane regions and a NHX domain. Overexpression of AtNHX4 in Arabidopsis transgenic lines confers greater salt tolerance than in wild type plants. These results suggest that AtNHX4 may encode a putative vacuolar NHX that plays an important role in salt tolerance.     

Molecular identification and phylosenetic position of Cuora yunnanensis

The Yunnan box turtle （Cuora yunnanensis, Boulenger, 1906）, which has drawn much attention in conservation biology, was regarded as extinct since it was previously known only from 12 specimens collected in Yunnan, China, before 1908. Recently, live specimens have been discovered which are suggested to be C. yunnanensis. To determine whether the newly discovered specimens are really C. yunnanensis, we have established a molecular phylogeny, with a 1725-bp fragment of mitochondrial DNA, using samples from three live individuals of C. yunnanensis, together with sequence data from a museum specimen of C. yunnanensis （MNHN 1907.10） and other members of the genus Cuora. We found that the three newly discovered individuals and the old museum specimen of C. yunnanensis are very similar both in morphology and in mitochondrial DNA sequence, suggesting that the three new individuals are the very C. yunnanensis, and thus the species is not extinct. Our phylogenetic analysis also demonstrates that C. yunnanensis is not of recent hybrid origin, but rather represents a distinct evolutionary lineage.     

Dissecting anxiety-related QTLs in mice by univariate and multivariate mapping

To detect genes underlying anxiety-related traits in mice,we performed univariate and multivariate QTL mapping analyses of phenotypes obtained from 71 mice of the BXD recombinant inbred (RI) strains (n=528 mice) and their parental strains (C57BL/6J and DBA/2J).Separate and joint mapping analyses were carried out using a linkage map composed of 506 simple sequence repeats (SSRs).The main QTL effects,interactions between pairs of QTLs (epistasis),and their environmental interactions were estimated.The results showed that anxiety-related traits were influenced by multiple QTLs (five main effect QTLs and three epistatic QTLs).Ten potential anxiety-related candidate genes within the QTL intervals on chromosomes 5,13 and 15 were identified.Some of these genes have been reported previously to be associated with the anxiety response.Based on our results,it is suggested that the multivariate QTL mapping approach improves the statistical power for detecting QTL and the precision of parameter estimation.Moreover,multivariate mapping can also detect pleiotropic QTL effects.     

孝顺竹笋箨全长转录组测序分析

【目的】揭示竹子笋箨的生物学功能与其衰老的分子基础。【方法】利用PacBio Sequel 3代全长转录组测序技术结合生物信息学方法,对孝顺竹不同衰老阶段笋箨全长转录本进行分析。【结果】共获得106 148 条平均长度为3 615 bp的全长转录本。多数转录本长度分布在1.4~7.0 kb。NCBI等注释显示,97.34%的转录本具有注释结果。Mercator注释分析显示笋箨转录本覆盖了其全部34 个类别,其中与蛋白类别相关的序列最多,达到了18 224 条;与microRNA类别相关的最少,仅有9 条。在重要功能基因方面,共获得了2 489 条激素代谢及信号转导相关序列; 8 804 条转录因子序列中393 条与光合作用相关的转录本。其中,注释到的189 条共30 类NAC(NAM/ATAF/CUC)转录因子基因,93.33%的NAC基因,共计28类,其表达量与笋箨衰老呈相关性, 包括已被证实在叶片衰老中具有重要调控作用的NAC002、NAC016、NAC017、NAC029(NAP)、NAC042、NAC055与NAC083等7个NAC转录因子与NAC014等14个被报道在叶片衰老中具有潜在作用的NAC基因。NAC025、NAC028、NAC045、NAC061、NAC086、NAC103与NAC1L (NAC 1 Like) 7 个NAC转录因子表达与孝顺竹笋箨衰老呈正相关,为新发现的正调控笋箨衰老的潜在转录因子。此外,利用MISA程序在76 499 个序列中检测到简单序列重复(simple sequence repeats,SSR)位点,主要以单核苷酸重复(SSRs)为主,占据了全部(SSRs)的55.2%。利用PLEK等软件共获得2 769个长链非编码RNA(long non-coding RNA, lncRNA)。【结论】竹子笋箨基因表达多样化,并具有完整的光合系统基因,显示其具有潜在光合功能;NAC转录因子在孝顺竹笋箨衰老中具有潜在重要调控作用。本研究首次解析了竹子笋箨的全长转录本特征,为今后深入分析竹子笋箨功能及其衰老的分子机制奠定基础。     

Studies of new EST-SSRs derived from Gossypium barbadense

Existing cotton EST-SSR markers are mostly derived from Gossypium arboreum and Gossypium hir-sutum, but EST-SSR markers from Gossypium barbadense are scarce. One hundred and nineteen EST-SSRs were developed based on 98 unique ESTs from a cDNA library constructed in our laboratory using developing fibers from G. barbadense cv. Pima3-79. Among the SSRs, trinucleotide AAG appeared at a high frequency of 11.76%. 36 accessions (consisting of 13 diploids of the A genome, 11 diploids of the D genome and 12 allotetraploids of the AD genome) were employed to test new EST-SSRs. 76 EST-SSRs were successfully amplified, and 313 polymorphic fragments were yielded, with an average of 4.11 fragments per primer pair. The PIC ranged from 0.17 to 0.95 with an average of 0.53. Based on Jaccard’s genetic similarity coefficient, these 36 accessions were clustered into three groups. 21 EST-SSRs exhibited polymorphisms in BC1 population ((Emian22 × Pima3-79) × Emian22), 24 polymor- phic loci were generated, while 22 of the 24 polymorphic loci were integrated with our interspecific BC1 backbone genetic linkage map, and anchored in 12 chromosomes. This study effectively proved that EST-SSRs from G. barbadense are valuable for genetic diversity analysis and genetic mapping.     

菠菜转录因子MYB基因家族的鉴定及表达分析

利用生物信息学手段鉴定了76个具有典型R结构的菠菜转录因子（Spinacia oleracea） MYB,其中包括72条R2R3-MYB基因（2R-MYB）和4条R1R2R3-MYB基因（3R-MYB）。通过生物信息学对菠菜MYB转录因子家族成员的理化性质、染色体定位、结构域序列保守性和系统进化关系进行了分析,结果表明：菠菜MYB家族有32个基因位于染色体正链,另外44个基因位于染色体反链;MYB的DNA结合域中的保守域主要位于两个R重复序列的第二和第三螺旋之间,结合域中每个R重复的第一和第二色氨酸之间的氨基酸序列相对不保守;根据菠菜、拟南芥及甜菜的MYB家族系统进化关系可以推测,菠菜MYB家族中56个成员可以按功能划分为4类,在菠菜的生长发育过程中可能起着重要的调节作用,其余成员中有7个MYB基因可能参与菠菜响应氮素浓度的氮素利用及生长发育进程。     

中国紫菜养殖

中国养殖紫菜,在目前,主要是群众性廉式养殖坛紫菜P.haitanensis和条斑紫菜P.yezoensis。其养殖可分为两个阶段:首先,在室內培养、促熟壳斑藻conchocelis和全人工采壳孢子conchospore苗(坛紫菜通常田喷壳孢子水法);然后,在潮间带进行半浮动筏(簾)式养成。     

Cloning and expression profiles of 15 genes encoding WRKY transcription factor in wheat （Triticum aestivem L.）

WRKY proteins are involved in various physiological processes, including biotic and abiotic stress responses, hormone responses and development. However, no systematic identification, expression and function analysis of WRKY genes in wheat were reported. In this study, we isolated 15 wheat cDNAs with complete open reading frame （ORF） encoding putative WRKY proteins using in silico cloning. Phylogenetic analysis indicated that the 15 wheat WRKY genes belonged to three major WRKY groups. Expression analysis revealed that most genes expressed drastically in leaf, except Ta WRKYIO which expressed in crown intensively. Four genes were strongly up-regulated with the senescence of leaves. Eight genes were responsive to low temperature, high temperature, NaCl or PEG treatment. Moreover, differential expression patterns were also observed between wheat hybrid and its parents, and some genes were more responsive to PEG treatment in the hybrid. These results demonstrated that wheat WRKY genes are involved in leaf senescing and abiotic stresses. And the changed expression of these WRKY genes in hybrid might contribute to the heterosis by improving the stress tolerance in hybrids.     

Submerged culture of Magnetospirillum gryphiswaldense under N2-fixing condition and regulation of activity of nitrogen fixation

A submerged culture technique for Magnetospirillum gryphiswaldense under the nitrogen-fixing condition (microaerobic and N-limited) was set up. In N-limited medium with Na-lactate as a sole carbon source, the optical density (A_600nm) and activity of nitrogen fixation of cells were 1.3 and 217 nmol of ethylene produced per hour per A_600nm respectively within 21 h by three times of feeds. The pH and temperature were controlled at 7.2 and 30℃ respectively, and the oxygen concentration was controlled by sparging with N₂ containing 0.4%—0.8% of O₂. The activity of nitrogen fixation of cells was obviously inhibited by oxygen and ammonium. It indicated that the posttranslational regulation of nitrogenase existed in M. gryphiswaldense.     

红树林土壤中假单胞菌群的重金属耐受活性

本研究从广西钦州茅尾海红树林保护区土壤中分离出的假单胞菌群着手，对其进行重金属耐受性能筛选，旨在获得可对重金属污染环境进行高效修复的菌株。选用传统稀释涂布法和16S rRNA基因序列的系统发育分析红树林土壤中的可培养假单胞菌；通过肉汤稀释法对菌株展开6种重金属耐受特性研究，对其进行耐受基因检测，并进一步测试重金属耐受活性显著菌株的重金属吸附性能。结果表明，从红树林土壤中共获得23株假单胞菌，18株菌对至少一种重金属表现出耐受活性，总阳性率为78.26%；10株菌的基因组DNA中检测出重金属耐受基因。其中，对铅和/或锌耐受活性显著的11株菌中，Y234、Y123、Y264、Y225和Y346具有显著的吸附性能，铅和/或锌的吸附率均大于50%。广西钦州茅尾海红树林保护区土壤中可培养假单胞菌群丰富，其重金属耐受活性显著，具有较大的挖掘潜力。     

不同轮伐期对巨尾桉人工林碳固存的影响

【目的】深入探讨不同轮伐期对巨尾桉人工林碳固存的影响,为从应对全球气候变化的视角确定合理轮伐期提供理论依据。【方法】以轮伐期为短(7a)、中(13a)和长周期(21a)的巨尾桉人工林为研究对象,通过对不同轮伐期桉树林分生物量、碳固存、年平均固碳量的分析,揭示不同轮伐期对桉树林分碳固存的影响。【结果】巨尾桉人工林的生物量碳随着轮伐期的延长而逐渐增加,由7a轮伐期的(75.81±5.12)t·C/hm~2增至13a轮伐期的(180.11±19.97)t·C/hm~2以及21a轮伐期的(192.08±16.50)t·C/hm~2,方差分析表明,13a和21a轮伐期的总生物量碳显著高于7a轮伐期,而13a和21a轮伐期之间的差异不显著。巨尾桉人工林土壤有机碳随轮伐期延长而显著降低,由7a轮伐期的(89.99±0.35)t·C/hm~2、13a轮伐期的(85.42±0.76)t·C/hm~2下降到21a轮伐期的(74.64±0.24)t·C/hm~2。7~13a仍是巨尾桉人工林固碳能力迅速增长期,年平均总生物量碳由7a时的10.78t·C/(hm~2·a)迅速提高到13a的19.54t·C/(hm~2·a),增长81%;21a时巨尾桉人工林进入固碳能力下降期,年平均总生物量碳降至3.78t·C/(hm~2·a),固碳能力只是13a的19.34%。【结论】在南亚热带,巨尾桉人工林的最佳轮伐期确定在13a左右较为适宜,这与经济效益的最大化一致。