BcSNPdb: bovine coding region single nucleotide polymorphisms located proximal to quantitative trait loci

Bovine coding region single nucleotide polymorphisms located proximal to quantitative trait loci were identified to facilitate bovine QTL fine mapping research. A total of 692,763 bovine SNPs was extracted from 39,432 UniGene clusters, and 53,446 candidate SNPs were found to be a depth >3. In order to validate the in silico SNPs experimentally, 186 animals representing 14 breeds and 100 mixed breeds were analyzed. Genotyping of 40 randomly selected candidate SNPs revealed that 43% of these SNPs ranged in frequency from 0.009 to 0.498. To identify non-synonymous SNPs and to correct for possible frameshift errors in the ESTs at the predicted SNP positions, we designed a program that determines coding regions by protein-sequence referencing, and identified 17,735 nsSNPs. The SNPs and bovine quantitative traits loci informations were integrated into a bovine SNP data: BcSNPdb (http://snugenome.snu.ac.kr/BtcSNP/). Currently there are 43 different kinds of quantitative traits available. Thus, these SNPs would serve as valuable resources for exploiting genomic variation that influence economically and agriculturally important traits in cows.     

Mining single nucleotide polymorphisms from EST data of silkworm, Bombyx mori, inbred strain Dazao

We made use of 81,635 expressed sequence tags (ESTs) derived from 12 different cDNA libraries of the silkworm, Bombyx mori, inbred strain Dazao (P50), to identify high-quality candidate single nucleotide polymorphisms (SNPs). By PHRAP assembling, 12,980 contigs containing 11,537 contigs assembled by more than one read were obtained, and 101 candidate SNPs and 27 single base insertions/deletions were identified from 117 contigs assembled from 1576 high-quality reads base-called with PHRED and screened on the basis of the neighborhood quality standard (NQS). Simultaneously, we also predicted 40 SNPs in coding regions (cSNPs), of which 26 were predicted to lead to amino acid non-synonymous variations and 14 synonymous substitutions. Also, the 1.66:1 ratio of transition/transversion is different from that of other insects. As the first SNP analysis of a Lepidoptera, B. mori, the single nucleotide polymorphic density is estimated to be 1.3 x 10(-3) by sequence diversity. This analysis shows that expressed sequences from multiple libraries may provide an abundant source of comparative reads to mine for cSNPs from the silkworm genome.     

The nucleotide sequence of the cDNA coding for the human dihydrofolic acid reductase

The nucleotide sequence of the human dihydrofolic acid reductase (DHFR) reading frame has been derived from the analysis of human DHFR cDNA. This sequence and the corresponding amino acid sequence have been compared with those available for the enzyme and its coding segment from other organisms. There is an 89% nucleotide sequence homology between the human DHFR reading frame and the mouse coding sequence. Furthermore, amino acid-sequence homologies of 74%, 81% and 89% has been found between human DHFR and chicken, bovine and mouse DHFR, respectively.     

Identification and Validation of Single Nucleotide Polymorphisms in Poplar Using Publicly Expressed Sequence Tags

By using assembled expressed sequence tags （ESTs） from 14 different eDNA libraries that contain 84 132 sequences reads, 556 Populus candidate single nucleotide polymorphisms （SNPs） were identified. Because traces were not available from dbEST （http：//www.ncbi.nlm.nih.gov/dbEST/index.html）, stringent filters were used to identify reliable candidate SNPs. Sequences analysis indicated that the main types of substitutions among candidate SNPs were A/G and T/C transitions, which accounted for 22.0% and 30.8%, respectively. One hundred and ten candidate SNPs were tested. As a result, 38 candidate SNPs were confirmed by directed sequencing of PCR products amplified from six different individuals. Thirteen new SNPs in intron regions were found and multiple SNPs were found to be located in both intron and exon regions of four contigs. Heterozygosis was found in all 47 candidate sites and five SNP sites were heterozygous in all six samples. This is the first report of SNP identification in a tree species which reveals that assembled ESTs from multiple libraries of the public database may provide a rich source of comparative sequences for an SNP search in the poplar genome.     

苹果酸度相关基因MdPH1的克隆、表达及亚细胞定位分析

以苹果属（Malus）植物沧江海棠（M.ombrophila Hand.-Mazz）的果实为材料,对其发育过程中苹果酸的含量进行测定,并结合转录组测序的方法筛选控制果实酸度的候选基因。结果显示：MdPH1候选基因的编码区包含2829 bp,编码942个氨基酸;基因组序列全长为4269 bp,包含8个外显子和7个内含子。对10份苹果种质资源中PH1基因序列的分析结果表明,该基因序列中存在22个单核苷酸多态性（SNP）,其中13个位于内含子区,9个位于外显子区;位于最后一个外显子上SNP（G/A）的变异导致了编码氨基酸从缬氨酸变为异亮氨酸。MdPH1蛋白包含8个跨膜结构域,其中蛋白N端包含3个跨膜结构域,C端包含5个跨膜结构域。系统进化分析结果显示,苹果中的PH家族成员与梨（Pyrus communis L.）中的PH家族成员聚集成一簇。组织特异性表达结果发现,MdPH1基因在苹果果实中的表达量最高,其次是叶、花和根,茎中表达量最低。亚细胞定位分析表明MdPH1蛋白定位于液泡膜上。     

The primary structure of chicken liver cytochrome b5 deduced from the DNA sequence of a cDNA clone

A cDNA clone encoding the chicken liver cytochrome b5 was isolated by probing a library with synthetic oligonucleotides based on a partial amino acid sequence of the protein. Determination of the DNA sequence indicated a 414-nucleotide open reading frame which encodes a 138-amino acid residue polypeptide. The open reading frame contains 6 amino acids at the amino terminus which were not present on any of the cytochrome b5 polypeptides for which the amino acid sequence has been determined directly, suggesting that the protein is proteolytically processed to the mature form. The results of genomic Southern analysis were consistent with the presence of two structurally different genes in the chicken genome, raising the possibility that the soluble and membrane-bound forms of the protein are the products of separate genes.     

Confirming single nucleotide polymorphisms from expressed sequence tag datasets derived from three cattle cDNA libraries

Using the Phred/Phrap/Polyphred/Consed pipeline established in the National Livestock Research Institute of Korea, we predicted candidate coding single nucleotide polymorphisms (cSNPs) from 7,600 expressed sequence tags (ESTs) derived from three cDNA libraries (liver, M. longissimus dorsi, and intermuscular fat) of Hanwoo (Korean native cattle) steers. From the 7,600 ESTs, 829 contigs comprising more than two EST reads were assembled using the Phrap assembler. Based on the contig analysis, 201 candidate cSNPs were identified in 129 contigs, in which transitions (69%) outnumbered transversions (31%). To verify whether the predicted cSNPs are real, 17 SNPs involved in lipid and energy metabolism were selected from the ESTs. Twelve of these were confirmed to be real while five were identified as artifacts, possibly due to expressed sequence tag sequence error. Further analysis of the 12 verified cSNPs was performed using the program BLASTX. Five were identified as nonsynonymous cSNPs, five were synonymous cSNPs, and two SNPs were located in 3'-UTRs. Our data indicated that a relatively high SNP prediction rate (71%) from a large EST database could produce abundant cSNPs rapidly, which can be used as valuable genetic markers in cattle.     

高海拔生境下怀玉山高山马铃薯和本土农家薯的全基因组重测序分析

为了探究怀玉山高山马铃薯只限于高海拔生境下的适应机制,本研究以高海拔生境下怀玉山高山马铃薯和怀玉山本土农家薯采后块茎萌发的芽为材料,进行了全基因组重测序分析。结果表明:HAP的总SNP数量为2 835 211,SNP的密度为3.925 212,编码区nsSNP总数为104 456;LAP的总SNP数量为3 968 618,SNP的密度为5.427 855,编码区nsSNP总数为141 060。经CV分析后,LAP和HAP共有319个nsSNP关联了315个基因,其中76个基因具有超过1个的nsSNP。nsSNP的防御反应、蛋白氨基酸磷酸化、蛋白激酶等GOterms显著富集。HAP的总Indel数量为159 797,编码区的Indel数量为1 850;LAP的总Indel数量为154 291,编码区的Indel数量为1 714;两组的总Indel数量远远小于两组的总SNP数。经差异分析后,共有1 629个Indel关联到了726个基因。Indel的防御反应、胁迫应答、细胞凋亡等GOterms显著富集。nsSNP和Indel的GOterms大多与高海拔生境和光合作用相关。nsSNP和Indel分布在端粒附近相对于中心粒具有更高密度的变异。LAP和HAP共获得了1 490个CNV,其中缺失(Deletion)为610个,中性(Neutral)为548个,扩增(Amplification)为332个。本研究结果可为高海拔生境下怀玉山高山马铃薯的SNP和Indel相关标记的开发、优异基因的挖掘提供重要理论依据。     

Cloning and sequence analysis of the human liver rhodanese: comparison with the bovine and chicken enzymes

The cDNA for the human rhodanese (thiosulfate: cyanide sulfurtransferase, EC 2.8.1.1), a nuclearly encoded protein of the mitochondrial matrix, was isolated from a human fetal liver cDNA library. Nucleotide sequence revealed an open reading frame coding for a polypeptide of 295 amino acids, which presented a 57% and 58% identity with the bovine and avian rhodanese, respectively. The analysis of the 5'-ends of the coding region gave no evidence for the presence of a cleavable signal sequence as found in other mitochondrial proteins. A comparison with two available amino acid sequences (cow and chicken) showed that sequence similarity is not restricted to the alpha-helices and beta-structures motifs which are remarkably superimposable in the two halves of bovine rhodanese, but extends to adjacent regions.     

Sequence analysis of the equine SLC26A2 gene locus on chromosome 14q15-->q21

The solute carrier family 26, member 2 (SLC26A2) gene belongs to a family of multifunctional anion exchangers. Mutations in the human SLC26A2 gene are associated with autosomal recessively inherited chondrodysplasias. Hence, we postulate that the equine SLC26A2 could be a candidate gene for conformational traits in horses. An equine BAC clone harboring the SLC26A2 gene was isolated. The complete 142,625 bp insert sequence of this clone was determined by transposon sequencing. Together with the SLC26A2 gene the BAC clone contains four genes, i.e. the macrophage colony stimulating factor 1 receptor precursor (CSF1R), KIAA0194 protein gene similar to the SMF protein (KIAA0194), a tigger transposable element derived 14 (TIGD14), the 3'-5'-cyclic GMP phosphodiesterase alpha-chain (EC 3.1.4.35) and one unidentified open reading frame. The equine SLC26A2 gene encompassing 6,152 bp consists of two exons. The complete open reading frame of 2,211 bp encodes a protein of 736 amino acids. A comparison of the amino acid sequence with other mammalian orthologs revealed homologies with identity in a range between 80% and 88%. By contrast, the equine SLC26A2 protein lacks five C-terminal amino acids. Four single nucleotide polymorphisms (SNP) were identified (three synonymous and one non-synonymous variant Ser210Leu) in the coding region by comparative sequencing of 50 DNA samples representing the German Riding horse. Allele frequencies and distribution were further evaluated in a variety of different breeds: Arabians (for all four SNPs), Old Kladrub Horses, Draught Horses (including Westphalian Draught Horses, Rheinish Westphalian Draught Horses, Saxon-Thuringia Coldbloods, Altmarker Coldbloods), American Saddlebreds, Miniature Horses, Australian Riding Ponies, Appaloosa, Morgan Horses, and Lipizzaner for C629T (Ser210Leu) alone. No animal carrying the homozygous genotype TT has been detected. The overall frequency of the newly described variant T is low (between 2% and 6%). Simulation studies on the protein conformation predict structural protein changes mediated by the SNP.     

A novel cDNA extension procedure. Isolation of chicken fatty acid synthase cDNA clones

We have developed a simple and versatile cDNA extension method using lambda-exonuclease-generated single-stranded DNA as a primer. This plasmid-based cDNA extension method can be used to synthesize unidirectional extensions of the existing cDNA clones or subcloned fragments of the untranslated and exon regions of genomic DNA clones. The method is simple to use and involves no addition of linkers or tailing. We have successfully used this method to isolate 4.6 kilobase pairs of chicken fatty acid synthase cDNA clones, starting from the fragment of a genomic clone coding for the untranslated region of the fatty acid synthase mRNA. About 2.8 kilobase pairs of the cDNA coding for the chicken fatty acid synthase has been sequenced. The sequence has an open reading frame coding for 945 amino acids of the fatty acid synthase. In the sequence, we have identified the enoyl reductase, NADPH binding region, a putative beta-ketoacyl reductase region, and the entire sequences of acyl carrier protein and the thioesterase domains. The arrangement of these partial activities in this sequence confirms the arrangement of these activities as determined through partial proteolytic mapping studies. The amino acid sequence of chicken fatty acid synthase deduced from cDNA sequences shows a high degree of homology with the rat fatty acid synthase sequence, suggesting that these multifunctional proteins are conserved evolutionarily.     

SNPs occur in regions with less genomic sequence conservation

Rates of SNPs (single nucleotide polymorphisms) and cross-species genomic sequence conservation reflect intra- and inter-species variation, respectively. Here, I report SNP rates and genomic sequence conservation adjacent to mRNA processing regions and show that, as expected, more SNPs occur in less conserved regions and that functional regions have fewer SNPs. Results are confirmed using both mouse and human data. Regions include protein start codons, 3' splice sites, 5' splice sites, protein stop codons, predicted miRNA binding sites, and polyadenylation sites. Throughout, SNP rates are lower and conservation is higher at regulatory sites. Within coding regions, SNP rates are highest and conservation is lowest at codon position three and the fewest SNPs are found at codon position two, reflecting codon degeneracy for amino acid encoding. Exon splice sites show high conservation and very low SNP rates, reflecting both splicing signals and protein coding. Relaxed constraint on the codon third position is dramatically seen when separating exonic SNP rates based on intron phase. At polyadenylation sites, a peak of conservation and low SNP rate occurs from 30 to 17 nt preceding the site. This region is highly enriched for the sequence AAUAAA, reflecting the location of the conserved polyA signal. miRNA 3' UTR target sites are predicted incorporating interspecies genomic sequence conservation; SNP rates are low in these sites, again showing fewer SNPs in conserved regions. Together, these results confirm that SNPs, reflecting recent genetic variation, occur more frequently in regions with less evolutionarily conservation.     

Molecular cloning of chicken metallothionein. Deduction of the complete amino acid sequence and analysis of expression using cloned cDNA.

A cDNA library was constructed using RNA isolated from the livers of chickens which had been treated with zinc. This library was screened with a RNA probe complementary to mouse metallothionein-I (MT), and eight chicken MT cDNA clones were obtained. All of the cDNA clones contained nucleotide sequences homologous to regions of the longest (376 bp) cDNA clone. The latter contained an open reading frame of 189 bp, and the deduced amino acid sequence indicates a protein of 63 amino acids of which 20 are cysteine residues. Amino acid composition and partial amino acid sequence analyses of purified chicken MT protein agreed with the amino acid composition and sequence deduced from the cloned cDNA. Amino acid sequence comparisons establish that chicken MT shares extensive homology with mammalian MTs, but is more closely related to the MT-II than to the MT-I isoforms from various mammals. The nucleotide sequence of the coding region of chicken MT shares approximately 70% homology with the consensus sequence for the mammalian MTs. Southern blot analysis of chicken DNA indicates that the chicken MT gene is not a part of a large family of related sequences, but rather is likely to be a unique gene sequence. In the chicken liver, levels of chicken MT mRNA were rapidly induced by metals (Cd2+, Zn2+, Cu2+), glucocorticoids and lipopolysaccharide. MT mRNA was present in low levels in embryonic liver and increased to high levels during the first week after hatching before decreasing again to the basal levels found in adult liver. The results of this study establish that MT is highly conserved between birds and mammals and is regulated in the chicken by agents which also regulate expression of mammalian MT genes. However, in contrast to the mammals, the results suggest the existence of a single isoform of MT in the chicken.     

Snipping polymorphisms from large EST collections in barley (<Emphasis Type="Italic">Hordeum vulgare </Emphasis>L.)

The public EST (expressed sequence tag) databases represent an enormous but heterogeneous repository of sequences, including many from a broad selection of plant species and a wide range of distinct varieties. The significant redundancy within large EST collections makes them an attractive resource for rapid pre-selection of candidate sequence polymorphisms. Here we present a strategy that allows rapid identification of candidate SNPs in barley (Hordeum vulgare L.) using publicly available EST databases. Analysis of 271,630 EST sequences from different cDNA libraries, representing 23 different barley varieties, resulted in the generation of 56,302 tentative consensus sequences. In all, 8171 of these unigene sequences are members of clusters with six or more ESTs. By applying a novel SNP detection algorithm (SNiPpER) to these sequences, we identified 3069 candidate inter-varietal SNPs. In order to verify these candidate SNPs, we selected a small subset of 63 present in 36 ESTs. Of the 63 SNPs selected, we were able to validate 54 (86%) using a direct sequencing approach. For further verification, 28 ESTs were mapped to distinct loci within the barley genome. The polymorphism information content (PIC) and nucleotide diversity () values of the SNPs identified by the SNiPpER algorithm are significantly higher than those that were obtained by random sequencing. This demonstrates the efficiency of our strategy for SNP identification and the cost-efficient development of EST-based SNP-markers.The first two authors contributed equally to this work     

Cloning and characterization of full-length coding sequence (CDS) of the ovine 37/67-kDa laminin receptor (<Emphasis Type="Italic">RPSA</Emphasis>)

The 37-kDa Laminin Receptor Precursor (LRP)/67-kDa Laminin Receptor (LR), also known as ribosomal protein SA (RPSA), had been identified as a putative cell surface receptor for prions. Herein, we isolated the full-length coding sequence (CDS) of the ovine 37/67-kDa LRP/LR gene and submitted it to the GenBank under accession number EF649775. The open reading frame (ORF) of the 37/67-kDa LRP/LR CDS is 885 bp in length, containing six exons encoding a protein of 295 amino acids. The nucleotide sequence presented here is well coincided with the whole ovine genome of the 37/67-kDa LRP/LR previously published. Moreover, we identified four novel single nucleotide polymorphism sites (SNPs) at position 324 in exon 4, positions at 809, 875, and 881 in exon 7, respectively. Further, based on the deduced amino acid sequence alignment of the 37/67-kDa LRP/LR from human, cattle, mice, pig, chicken, and sheep, we also identified three polymorphic amino acid sites (PAAs) at residues 241, 272, and a novel site at residue 270 in the putative indirect prion protein (PrP) interaction region (180–285) on 37/67-kDa LRP/LR. Prediction of protein secondary structure further indicated that PAAs at residues 241, 270 and 272 may cause protein conformation changes as predicted, which may affect on the binding with prion protein. In addition, multiple-tissues RT-PCR results revealed that 37/67-kDa LRP/LR mRNA is expressed in all the 11 selected ovine tissues.