家蚕地方品种遗传多样性及其分子系统学研究

为检测家蚕(Bombyx mori)地方品种间的遗传差异,分析其遗传多样性和系统分化模式,为家蚕的品种鉴定和杂交育种亲本选择提供分子依据,本文对72个家蚕地方品种的淀粉酶基因(Bmamy2)片段进行测序,采用DNA序列分析技术探究了家蚕地方品种在核苷酸序列水平上的遗传差异及其系统分化特征.结果显示,家蚕地方品种在Bmamy2基因变异丰富,碱基突变率为5.6-8.2％,单倍型多样度为0.8294,核苷酸多样度为0.0236±0.00122,表明Bmamy2是一个具有丰富序列多样性的标记基因,能够较好地区分家蚕地方品种间的遗传多样性.系统发育分析、分子方差分析和Network分析等多种分析方法均检测到家蚕地方品种不按地理系统或生态类型聚类,各聚类种群均由来自不同地理系统或生态类型的品种构成,而来自同一地理系统或生态类型的品种则分别归属于不同的聚类种群.推测在家蚕的进化过程中并没有形成按地理系统或生态类型归类的分化格局,家蚕可能具有多种化性混合驯化起源的遗传背景和系统分化特征.     

中国野桑蚕遗传多样性的AFLP分析

采用AFLP技术对我国具有代表性的7个省市的10个野桑蚕(Bombyx mandarina)种群和2个家蚕（Bombyx mori)品种的遗传多样性进行了研究。结果表明：杭州、陕西和重庆3个地区的野桑蚕种群间及种群内个体间的遗传距离都比家蚕品种大。野桑蚕存在广泛的变异,7个省市的10个野桑蚕种群之间的遗传距离为0．164－0．444（平均值为0．3826）,平均杂合度为0．7061,表明野桑蚕自然群体的遗传多样性十分丰富。     

云南洞密蛛的种群遗传结构初探

云南洞密蛛Trogloneta yunnanense（Song&Zhu, 1994）是生活在云贵高原的一种洞穴蜘蛛。本文基于6个洞穴种群159只个体样本的线粒体COⅠ基因,初步探讨了该物种的种群结构和遗传多样性特征。研究结果表明：1)159个样本共检测到16个单倍型,种群间无共享单倍型,遗传多样性呈现出总体高、种群内低、单倍型多样性高、核苷酸多样性低的模式;2)种群间遗传分化极显著,种群内变异小,种群间F_ST值均在0.9以上（P<0.01）,种群间基因流小（Nm=0.01）,种群间变异占总变异的95.75%,种群内变异仅占4.25%;3)中性检验和核苷酸错配分布检验显示,云南洞密蛛未经历种群扩张,群体大小保持稳定。推测洞穴隔离和地下极端环境条件可能维持着云南洞密蛛的种群规模稳定,同时促进其种群间的遗传分化。     

云南中缅树鼩线粒体细胞色素b和D-loop区遗传多样性研究

对3个地方种群的49只样本的线粒体Cyt b基因全序列（1140 bp）及33只样本的控制区D-loop基因区段（745 bp）进行了序列测定。结果表明：Cyt b基因多态性位点有47个,其中单变异位点位点23个,简约信息位点24个。共定义了24个单倍型,其中种群间的共享单倍型有2个（8.33%）,其余均为某个种群所特有,单倍型多样度范围为0.80952（剑川种群）~0.91532（禄劝种群）,核苷酸多样度指数介于0.00326（禄劝种群）~0.00635（剑川种群）之间;D-loop基因多态性位点有18个,其中单变异位点8个,简约信息位点10个。共定义了16个单倍型,无种群间的共享单倍型,单倍型多样度范围为0.76615（禄劝种群）~0.93333（丽江种群）,核苷酸多样度指数介于0.00269（禄劝种群）~0.00583（丽江种群）之间。从各单倍型的TCS网络进化图显示横断山种群位于分支的末端,表现出中缅树鼩由南向北的扩散模式,支持＂岛屿起源＂假说。     

中国东北地区远东鼩鼱遗传多样性与遗传结构

对中国东北地区3个种群(大兴安岭、小兴安岭和长白山山脉)的远东鼩鼱(77个样本)Cyt b基因全序列进行分析,共获得64个单倍型。整体单倍型多态性为0.9920,核苷酸多态性为0.0105,表明该地区远东鼩鼱具有较高遗传多样性,且长白山山脉远东鼩鼱种群遗传多样性明显高于大兴安岭和小兴安岭种群。F-统计量、遗传相似系数和遗传距离分析结果均显示,种群间和采样地间的遗传距离与地理距离基本相符。方差分析显示,种群间的变异占总变异的33.4%,种群内的采样地间变异占总变异的10.2%,采样点内部变异占总变异的56.4%。种群历史分析显示,东北地区远东鼩鼱未经历过数量扩张。从GenBank下载了欧亚其他地区远东鼩鼱序列进行遗传结构研究。远东鼩鼱系统发生树分化为2大分支:一大支主要由大兴安岭和小兴安岭种群构成,两种群具有一定分化;另一大支又分为两个分支。中介网络图显示,远东鼩鼱具有3个谱系:一个谱系主要由大兴安岭和小兴安岭的单倍型样本构成,还包括长白山山脉的4个单倍型样本;另一谱系包括来自于中国东北地区3个种群的个别单倍型,还包括俄罗斯贝加尔湖单倍型和芬兰单倍型;最后一个谱系完全是由长白山山脉单倍型构成。遗传多样性、系统发生树和中介网络图结果均表明,长白山山脉为远东鼩鼱末次冰期避难所。     

江西井冈山地区灰胸竹鸡的遗传多样性研究

灰胸竹鸡Bambusicola thoracica是我国特有鸟类.本文采用聚合链式反应和直接测序的方法测定井冈山地区灰胸竹鸡3个种群线粒体DNA(mtDNA)控制区1142 bp的序列,分析其序列变异和种群遗传多样性.30个样本共发现16个变异位点和10种单倍型,其中hapl广泛分布,占所分析样本的23.33%,是其祖先单倍型.3个种群的平均单倍型多样性和核苷酸多样性分别为0.815和0.00243.青原区种群与其它两个种群遗传分化显著,基因交流受限制.受隔离影响,青原区种群遗传多样性最低.在系统发生树上,10种单倍型形成两支,井冈山种群和永新县种群聚在一起,与其地理位置相一致.     

大石鸡边缘种群的遗传结构

大石鸡（Alectoris magna）分布于青海西部、东部,甘肃中部和宁夏西部的干部和半干旱区,由于环境变化,其种群向甘肃南部森林被砍伐的地区扩散,形成涟缘种群。本研究采用聚合酶链式反应（PCR）和直接测序的方法获得了采自甘肃的大石鸡一个边缘种群和两个中心地理种群共39个个体的线粒体DNA（mtDNA）控制区基因（D－loop）456－457个核苷酸的基因序列,16个变异位点（占整个序列的3.5%）有15个单倍型。边缘扩散种群有3种单倍型,M6与另两个种群共有,单倍型频率为0.108,而M4和M5为其特有,单倍型频率分别为0.081和0.027。边缘种群的单倍型比率和遗传多样性分别为37.5%和0.549。中心地理种群1和种群2各有7个单倍型,除M6炳种群共有,其余为自所特有,单倍型比率分别为46.7%和50.0%,遗传多样性分别是0.729和0.786。边缘种群的单倍型比率和遗传多样性均低于中心地理种群。     

基于线粒体Cyt b基因的雉鸡甘肃亚种的种群遗传结构

采用聚合酶链式反应(PCR)和直接测序法分析雉鸡甘肃亚种的种群遗传结构。测序获得了8个种群79个样本的长819bp的线粒体细胞色素b(Cytb)基因序列,发现7个变异位点,界定了6种单倍型(Haplotype),其中Hap-1是共享单倍型,占79.75%。在8个种群中,除ZJC和CHX种群外,其它种群的核苷酸多样性和单倍型多样性都很低,种群间遗传分化小,基因流大。ZJC种群由于亚种间的基因渗透,导致其遗传多样性最高。CHX种群受徽成盆地周边山地的森林隔离作用,单倍型最多。分子方差分析(AMOVA)表明,遗传变异主要发生在种群内,占86.44%。种群间和地理种群组间差异不显著。地理种群组的核苷酸多样性(π)同形态上的变异规律一致。FS值检测表明,雉鸡甘肃亚种曾发生过快速扩张。     

赤水河两种荷马条鳅属鱼类的遗传多样性及谱系生物地理学过程分析

研究以线粒体Cyt b基因为分子标记,对赤水河两种荷马条鳅属（Homatula）鱼类（红尾荷马条鳅和短体荷马条鳅）的遗传多样性及种群结构进行了分析;同时,结合了对在长江上游其他几个水系分布的同种鱼类进行比较,分析其生物地理学过程。遗传多样性分析结果表明,5个水系135尾红尾荷马条鳅Cyt b基因序列共检测出42个单倍型,单倍型多样性和核苷酸多样性分别为0.936和0.00493;其中,赤水河种群的分别为0.891和0.00208。3个水系52尾短体荷马条鳅Cyt b基因序列共检测出12个单倍型,单倍型多样性和核苷酸多样性分别0.821和0.01105;赤水河种群的分别为0.646和0.00390。基于ML和BI法构建的单倍型分子系统发育树结果表明,两物种各自构成单系,且得到较强支持。红尾荷马条鳅各种群地理格局分布明显,赤水河种群为并系类群位于分支基部。在短体荷马条鳅支系中,岷江与沱江两个水系的个体相互聚类在一起,而赤水河群体聚成一个分支。赤水河与其他水系不存在共享单倍型,表现了明显的隔离和差异性的地理分布格局。由于这些水系之间地理位置相距较远,推测这种格局的形成不是地质运动造成的水系隔离,而是历史时期水位的高低变化造成鱼类种群的扩散和隔离。错配分析支持赤水河两物种种群扩张的推断,但中性检验却并非全部支持,显示种群历史相对复杂。     

基于nrDNA GapC基因内含子序列的资源冷杉遗传多样性研究

通过nrDNA GapC基因内含子序列的测序和分析,揭示资源冷杉遗传多样性水平和种群分化的强弱并推断其进化历史。资源冷杉3个种群的34个个体共获得70条GapC基因内含子序列,有8个核苷酸变异位点,鉴别出12种单倍型。银竹老山、大院和舜黄山种群分别有10种、6种和7种单倍型;有7个个体得到了2种以上的单倍型。资源冷杉物种水平的单倍型多样性（h）为0.817 0,种群的在0.683 3～0.883 1之间;物种水平的核苷酸多样性（π）为0.003 90,种群的在0.002 63～0.003 82之间。种群遗传分化研究结果（G_st=0.103,p<0.05）说明种群间存在显著的遗传分化,分子方差分析（AMOVA）结果也证明虽然大多数的核苷酸多态性（88.64％,p<0.001）来源于种群内个体间的变异,但是仍有显著性比例存在种群间（11.36%,p<0.05）。Tajima’s D、Fu &; Li’s D、Fu &; Li’s Fs中性检验结果都表明资源冷杉在物种和种群水平上都不拒绝中性进化。资源冷杉单倍型的谱系没有出现地区特异性谱系分支,核苷酸不配对分析（mismatch analysis）结果表明没有发生近期的群体扩张,资源冷杉的N_st（0.131）与G_st（0.103）没有显著性差异（p>0.05）,表明种群没有明显的地理结构。推测现存资源冷杉是相对较近的时间内片断化的产物。     

Evidence of selection at melanin synthesis pathway loci during silkworm domestication

The domesticated silkworm (Bombyx mori) was domesticated from wild silkworm (Bombyx mandarina) more than 5,000 years ago. During domestication, body color between B. mandarina and B. mori changed dramatically. However, the molecular mechanism of the silkworm body color transition is not known. In the present study, we examined within- and between-species nucleotide diversity for eight silkworm melanin synthesis pathway genes, which play a key role in cuticular pigmentation of insects. Our results showed that the genetic diversity of B. mori was significantly lower than that of B. mandarina and 40.7% of the genetic diversity of wild silkworm was lost in domesticated silkworm. We also examined whether position effect exists among melanin synthesis pathway genes in B. mandarina and B. mori. We found that the upstream genes have significantly lower levels of genetic diversity than the downstream genes, supporting a functional constraint hypothesis (FCH) of metabolic pathway, that is, upstream enzymes are under greater selective constraint than downstream enzymes because upstream enzymes participate in biosynthesis of a number of metabolites. We also investigated whether some of the melanin synthesis pathway genes experienced selection during domestication. Neutrality test, coalescent simulation, as well as network and phylogenetic analyses showed that tyrosine hydroxylase (TH) gene was a domestication locus. Sequence analysis further suggested that a putative expression enhancer (Abd-B-binding site) in the intron of TH gene might be disrupted during domestication. TH is the rate-limiting enzyme of melanin synthesis pathway in insects. Real-time polymerase chain reaction assay did show that the relative expression levels of TH gene in B. mori were significantly lower than that in B. mandarina at three different developmental stages, which is consistent with light body color of domesticated silkworm relative to wild silkworm. Therefore, we speculated that expression change of TH gene may contribute to the body color transition from B. mandarina to B. mori. Our results emphasize the exceptional role of gene expression regulation in morphological transition of domesticated animals.     

M chromosome of the wild silkworm, Bombyx mandarina (n = 27), corresponds to two chromosomes in the domesticated silkworm, Bombyx mori (n = 28).

Chromosomes of Bombyx mori (n = 28) and of Bombyx mandarina (n = 27) were studied cytogenetically to resolve the origin of the large M chromosome in the Japaneses type of B. mandarina. In the F1 progeny from the reciprocal cross between B. mandarina and B. mori, the mitotic chromosome number was 2n = 55, and a chromosome configuration of 26 bivalents plus 1 trivalent was observed at metaphase I of germ cells. The trivalent chromosome consisted of the M chromosome from B. mandarina and two chromosomes from B. mori. When males of B. mori were mated to the F1 females, nuclei with two types of chromosome number (2n = 55 and 2n = 56) and two sets of chromosome pairs (26 bivalents plus 1 trivalent versus 28 bivalents) were observed in the metaphase I stage. Linkage analysis showed that the 14th chromosome of B. mori was involved in these two types of chromosome segregation. This result indicates that the M chromosome in B. mandarina arose from a fusion between a chromosome corresponding to the 14th linkage group and another, yet unidentified linkage group.     

Nucleotide sequence variation in mitochondrial COI gene among 147 silkworm (Bombyx mori) strains from Japanese, Chinese, European and moltinism classes

We characterized the nucleotide sequences of PCR-amplified mitochondrial COI fragments of 147 silkworm (Bombyx mori) strains that have been maintained in the National Institute of Agrobiological Sciences. Coding sequences (714 bp) of the 147 COI fragments were classified into eight haplotypes based on nucleotide differences at eight segregating sites. No length variation was identified in this region. The 5'-noncoding region showed different features, wherein changes in the number of Ts in the T-stretch, together with two base substitutions, were observed. As a result, the 147 COI noncoding sequences were classified into six haplotypes. Combining the coding and noncoding regions, we identified 14 haplotypes. One of the 14 haplotypes, Hap1A was exclusively abundant in the Japanese native strain class, while this haplotype was less frequent in the other three native strain classes. This finding suggests that the Japanese strain class underwent significant genetic differentiation from the Chinese, European, and moltinism classes, when the each class is regarded as a population. Comparison of the nucleotide sequences to those of B. mandarina (which inhabits Japan) revealed changes that are significantly larger than those within either B. mori or B. mandarina. Furthermore, we detected no common haplotypes between them, which suggests the concept of suppressed gene flow between the two species.     

野桑蚕性信息素结合蛋白基因pbp1、气味受体基因or1和or3的克隆及其与家蚕同源基因的进化分析

昆虫能够特异性识别同类异性。雄蚕蛾对雌蚕蛾感知定位过程中, 性信息素结合蛋白PBP1、气味受体OR1和OR3起重要作用。为研究家蚕Bombyx mori和野桑蚕Bombyx mandarina杂交困难的分子机制, 了解性识别相关基因的进化, 本研究克隆得到了野桑蚕的性信息素结合蛋白基因pbp1（GenBank注册号：GQ246497）和气味受体基因or1（Genank注册号：GQ246496）和or3（GenBank注册号：GQ246498）。序列分析发现, 家蚕与野桑蚕相比, pbp1基因存在4个SNP位点, 分别为C10A, A40T, T270C和A333G, 其中2个SNP位点引起氨基酸的改变, 分别为Q→K和N→Y; or1基因存在5个SNP位点, 分别为T910C, A1147C, A1192T, T1276C和G1282A, 其中1个SNP位点引起氨基酸F→L的改变; or3基因存在4个SNP位点, 分别为A507G, A513G, T605C和G672A, 其中1个SNP位点引起氨基酸I→T的改变。3个基因的遗传距离很近, 进化速率也很慢。氨基酸的分子量和等电点有细微差异或无差异。PHD预测的二级结构表明, 变异位点对附近区域的结构没有任何影响, 功能位点也没有变化。推测家蚕与野桑蚕之间, 这些基因功能可能没有差异, 即二者的雌雄性个体间可以相互感知、识别, 这与实验观察结果一致。     

Diversity in copy number and structure of a silkworm morphogenetic gene as a result of domestication

The carotenoid-binding protein (CBP) of the domesticated silkworm, Bombyx mori, a major determinant of cocoon color, is likely to have been substantially influenced by domestication of this species. We analyzed the structure of the CBP gene in multiple strains of B. mori, in multiple individuals of the wild silkworm, B. mandarina (the putative wild ancestor of B. mori), and in a number of other lepidopterans. We found the CBP gene copy number in genomic DNA to vary widely among B. mori strains, ranging from 1 to 20. The copies of CBP are of several types, based on the presence of a retrotransposon or partial deletion of the coding sequence. In contrast to B. mori, B. mandarina was found to possess a single copy of CBP without the retrotransposon insertion, regardless of habitat. Several other lepidopterans were found to contain sequences homologous to CBP, revealing that this gene is evolutionarily conserved in the lepidopteran lineage. Thus, domestication can generate significant diversity of gene copy number and structure over a relatively short evolutionary time.     

Survey of long terminal repeat retrotransposons of domesticated silkworm (Bombyx mori)

Long terminal retrotransposons are major components of eukaryotic transposable elements. We have surveyed the long terminal repeats (LTR) retrotransposons of domesticated silkworm (Bombyx mori) by mining the data produced by Bombyx mori Genome Sequencing Project. At least 29 separate families of LTR retrotransposons are identified in this survey, comprising of 11.8% of the complete sequence. Families of domesticated silkworm LTR retrotransposons can be mainly classified into three groups: gypsy-like, copia-like, Pao-Bel. Fourteen families identified consist of gypsy-like elements, four families consist of copia-like elements and seven families consist of Pao-Bel elements. In addition to the three groups of LTR retrotransposons, two families of unusual non-coding elements are identified in the genome of this species. Further phylogenetic analysis of RT domain indicates that the elements of B.mori show high diversity and can form different clades in each group. An analysis of sequence variation from different families reveals distinct patterns of variation for the elements belonging to three groups. The analysis of the domesticated silkworm LTR retrotransposons should assist in our understanding of the roles of retroelement in lepidopteron insect genome evolution.     

The genome sequence of silkworm, Bombyx mori.

We performed threefold shotgun sequencing of the silkworm (Bombyx mori) genome to obtain a draft sequence and establish a basic resource for comprehensive genome analysis. By using the newly developed RAMEN assembler, the sequence data derived from whole-genome shotgun (WGS) sequencing were assembled into 49,345 scaffolds that span a total length of 514 Mb including gaps and 387 Mb without gaps. Because the genome size of the silkworm is estimated to be 530 Mb, almost 97% of the genome has been organized in scaffolds, of which 75% has been sequenced. By carrying out a BLAST search for 50 characteristic Bombyx genes and 11,202 non-redundant expressed sequence tags (ESTs) in a Bombyx EST database against the WGS sequence data, we evaluated the validity of the sequence for elucidating the majority of silkworm genes. Analysis of the WGS data revealed that the silkworm genome contains many repetitive sequences with an average length of <500 bp. These repetitive sequences appear to have been derived from truncated transposons, which are interspersed at 2.5- to 3-kb intervals throughout the genome. This pattern suggests that silkworm may have an active mechanism that promotes removal of transposons from the genome. We also found evidence for insertions of mitochondrial DNA fragments at 9 sites. A search for Bombyx orthologs to Drosophila genes controlling sex determination in the WGS data revealed 11 Bombyx genes and suggested that the sex-determining systems differ profoundly between the two species.     

Molecular phylogeny of the domesticated silkworm,Bombyx mori, based on the sequences of mitochondrialcytochrome b genes

Pupae from the Chinese wild mulberry silkworm,Bombyx mandarina, and 11 representative strains of the domesticated silkworm,Bombyx mori were selected for preparation of mitochondrial DNA. The 5′-end fragments ofcytochrome b genes (Cytb) were generated by polymerase chain reaction products and sequenced directly. The homologous sequences of the JapaneseB. mandarina and three strains ofB. mori were from the GenBank database. The sequences of the 16 silkworm strains were analysed with DNASTAR software and a phylogenic tree was constructed using PHYLIP software. The result showed that: (i) The sequence divergence between the strains ofB. mori and the JapaneseB. mandarina was larger (5.4% ≈ 5.8%) compared with that between strains ofB. mori and the ChineseB. mandarina (0.8% ≈ 1.9%). Analysis of clustering also showed that the sequences ofB. mori strains and ChineseB. mandarina clustered into group (B group), while that of JapaneseB. mandarina (A group) was outside this cluster. This may be evidence for the hypothesis thatB. mori originated from ChineseB. mandarina. (ii) Among 14 strains ofB. mori, sequence divergence was small and the most divergence was seen between strains Yanhe-1 and Chuxiong, whose sequences branched off from those of the otherB. mori strains on the phylogenetic tree. From this and from historical records, we infer that the strains Yanhe-1 and Chuxiong originated independently from southwest China.     

Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina

Bombyx mori and Bombyx mandarina are morphologically and physiologically similar. In this study, we compared the nucleotide variations in the complete mitochondrial (mt) genomes between the domesticated silkmoth, B. mori, and its wild ancestors, Chinese B. mandarina (ChBm) and Japanese B. mandarina (JaBm). The sequence divergence and transition mutation ratio between B. mori and ChBm are significantly smaller than those observed between B. mori and JaBm. The preference of transition by DNA strands between B. mori and ChBm is consistent with that between B. mori and JaBm, however, the regional variation in nucleotide substitution rate shows a different feature. These results suggest that the ChBm mt genome is not undergoing the same evolutionary process as JaBm, providing evidence for selection on mtDNA. Moreover, investigation of the nucleotide sequence divergence in the A+T-rich region of Bombyx mt genomes also provides evidence for the assumption that the A+T-rich region might not be the fastest evolving region of the mtDNA of insects.