意大利蜜蜂染色体DNA非编码区抗白垩病相关的SNP的筛选与验证

【目的】筛选验证意大利蜜蜂Apis mellifera ligustica染色体DNA非编码区与抗白垩病相关的SNP。【方法】本研究将蜜蜂球囊菌Ascosphaera apis孢子接种于人工饲养的意大利蜜蜂3日龄幼虫,根据是否存在白垩病症状进而筛选出抗病个体和易感个体。基于前期重测序结果中意大利蜜蜂第2和11号染色体DNA非编区与抗白垩病相关的SNP信息,利用PCR测序的方法筛选并验证意大利蜜蜂幼虫第2和11号染色体DNA非编码区与幼虫抗白垩病相关的55个SNP。【结果】发现位于意大利蜜蜂第11号染色体LOC100578413基因5′端的非编码区的SNP(T14570310C)在抗病个体中T等位基因频率高于C等位基因频率,且抗病个体中的T等位基因频率显著高于易感幼虫中的T等位基因频率,表明该SNP位点与抗白垩病相关。该分子标记对抗性个体和易感个体的判断结果与前期筛选的编码区SNP(C2587245T)分子标记的结果一致。【结论】筛选并验证意大利蜜蜂第11号染色体DNA非编码区的SNP(T14570310C)与抗白垩病相关。该位点为抗白垩病分子辅助选育提供新的分子标记,在意大利蜜蜂白垩病早期检测和培育白垩病抗性的蜂种方面具有重要意义。     

Interleukin-1beta isolated from a marine fish reveals up-regulated expression in macrophages following activation with lipopolysaccharide and lymphokines.

The gilthead seabream IL-1beta gene consists of five exons/four introns. The complete coding sequence contains a 102 bp 5' untranslated region (UTR), a single open reading frame of 762 bp which translates into a 253 amino acid molecule, and a 407 bp 3'UTR with a polyadenylation signal 14 nucleotides upstream of the poly(A)tail. The seabream sequence has the highest degree of nucleotide (61.7%) and amino acid (53%) identity with the trout IL-1beta sequences. The IL-1beta message was detected by RT-PCR in head-kidney, blood, spleen, liver, gill and peritoneal exudate of both non-infected and Vibrio anguillarum-challenged fish. More importantly, IL-1beta was highly expressed by purified macrophage monolayers and was up-regulated by lipopolysaccharide and lymphocyte-derived macrophage-activating factor stimulation.     

Development and testing of a combined species SNP array for the European seabass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata)

SNP arrays are powerful tools for high-resolution studies of the genetic basis of complex traits, facilitating both selective breeding and population genomic research. The European seabass (Dicentrarchus labrax) and the gilthead seabream (Sparus aurata) are the two most important fish species for Mediterranean aquaculture. While selective breeding programmes increasingly underpin stock supply for this industry, genomic selection is not yet widespread. Genomic selection has major potential to expedite genetic gain, particularly for traits practically impossible to measure on selection candidates, such as disease resistance and fillet characteristics. The aim of our study was to design a combined-species 60 K SNP array for European seabass and gilthead seabream, and to test its performance on farmed and wild populations from numerous locations throughout the species range. To achieve this, high coverage Illumina whole-genome sequencing of pooled samples was performed for 24 populations of European seabass and 27 populations of gilthead seabream. This resulted in a database of ~20 million SNPs per species, which were then filtered to identify high-quality variants and create the final set for the development of the ‘MedFish’ SNP array. The array was then tested by genotyping a subset of the discovery populations, highlighting a high conversion rate to functioning polymorphic assays on the array (92% in seabass; 89% in seabream) and repeatability (99.4–99.7%). The platform interrogates ~30 K markers in each species, includes features such as SNPs previously shown to be associated with performance traits, and is enriched for SNPs predicted to have high functional effects on proteins. The array was demonstrated to be effective at detecting population structure across a wide range of fish populations from diverse geographical origins, and to examine the extent of haplotype sharing among Mediterranean farmed fish populations. In conclusion, the new MedFish array enables efficient and accurate high-throughput genotyping for genome-wide distributed SNPs for each fish species, and will facilitate stock management, population genomics approaches, and acceleration of selective breeding through genomic selection.     

The effect of SNP discovery method and sample size on estimation of population genetic data for Chinese and Indian rhesus macaques (<Emphasis Type="Italic">Macaca mulatta</Emphasis>)

This study was designed to address issues regarding sample size and marker location that have arisen from the discovery of SNPs in the genomes of poorly characterized primate species and the application of these markers to the study of primate population genetics. We predict the effect of discovery sample size on the probability of discovering both rare and common SNPs and then compare this prediction with the proportion of common and rare SNPs discovered when different numbers of individuals are sequenced. Second, we examine the effect of genomic region on estimates of common population genetic data, comparing markers from both coding and non-coding regions of the rhesus macaque genome and the population genetic data calculated from these markers, to measure the degree and direction of bias introduced by SNPs located in coding versus non-coding regions of the genome. We found that both discovery sample size and genomic region surveyed affect SNP marker attributes and population genetic estimates, even when these are calculated from an expanded data set containing more individuals than the original discovery data set. Although none of the SNP detection methods or genomic regions tested in this study was completely uninformative, these results show that each has a different kind of genetic variation that is suitable for different purposes, and each introduces specific types of bias. Given that each SNP marker has an individual evolutionary history, we calculated that the most complete and unbiased representation of the genetic diversity present in the individual can be obtained by incorporating at least 10 individuals into the discovery sample set, to ensure the discovery of both common and rare polymorphisms.     

Characteristics of X- and Y-chromosome specific regions of the amelogenin gene and a PCR-based method for sex identification in red deer (<Emphasis Type="Italic">Cervus elaphus</Emphasis>)

The present study attempts to analyse sequences of the X- and Y-chromosome specific regions of the amelogenin (AMEL) gene in red deer. To this end, primers specific for each form of the gene (AMELX and AMELY) were designed based on bovine genomic sequences and the homologous regions of the genes were sequenced. The obtained sequence of AMELX gene showed high similarity with the corresponding region in cattle (91%) and humans (77%), but this similarity was slightly lower among AMELY genes and showed 87 and 73% of identical nucleotides, respectively. In addition, three single nucleotide polymorphisms (SNPs) were found in the AMELX gene of the female red deer investigated. Comparative analysis of the homologous fragments of the red deer AMELX and AMELY genes confirmed the deletion of an AMELY gene fragment in relation to AMELX. Homology of both sequences was 82% of identical nucleotides in the coding region and 74% in 3′ non-coding sequence. The sequences studied showed considerable similarity to homologous fragments of the human and bovine gene, but the structural differences observed lead us to design PCR-based method for sex identification in red deer, based on the presented sequences.     

Fibroblast growth factor receptor 2 (FGFR2): genomic sequence and variations

Fibroblast growth factor receptors (FGFRs) play an important role in development and tumorigenesis. Mutations in FGFR2 cause more than five craniosynostosis syndromes. The FGFR2 genomic structure is the largest of the FGFR family. We have refined and extended the genomic organization of the FGFR2 gene by sequencing more than 119 kb of PACs, cosmids, and PCR products and assembling a region of approximately 175 kb. Although the gene structure has been reported to include only 20 exons, we have verified the presence of at least 22 exons, some of which are alternatively spliced. The sizes of six exons differed from those reported previously. Comparison of our sequence and those in the NCBI database detected more than 300 potential single nucleotide polymorphisms (SNPs). However, sequencing regions containing 52 of these potential SNPs verified only 14 in PCR products generated from 16 CEPH alleles. In contrast, direct sequencing of the CEPH DNAs revealed 21 other polymorphisms. Only one SNP was found in the 2,926 bp of coding sequence. Twenty-seven SNPs, two insertion polymorphisms and five microsatellite polymorphisms are contained in approximately 16.6 kb of non-coding sequence. These data yield an average of one polymorphism for approximately 488 bp of non-coding sequence examined. This collection of SNP, insertion, and repeat polymorphisms will aid future association studies between the FGFR2 gene and human disease and will enhance mutation detection.     

Development and evaluation of single-nucleotide polymorphism markers in allotetraploid rapeseed (<Emphasis Type="Italic">Brassica napus</Emphasis> L.)

Single-nucleotide polymorphisms (SNPs) and insertion–deletions (INDELs) are currently the important classes of genetic markers for major crop species. In this study, methods for developing SNP markers in rapeseed (Brassica napus L.) and their in silico mapping and use for genotyping are demonstrated. For the development of SNP and INDEL markers, 181 fragments from 121 different gene sequences spanning 86 kb were examined. A combination of different selection methods (genome-specific amplification, hetero-duplex analysis and sequence analysis) allowed the detection of 18 singular fragments that showed a total of 87 SNPs and 6 INDELs between 6 different rapeseed varieties. The average frequency of sequence polymorphism was estimated to be one SNP every 247 bp and one INDEL every 3,583 bp. Most SNPs and INDELs were found in non-coding regions. Polymorphism information content values for SNP markers ranged between 0.02 and 0.50 in a set of 86 varieties. Using comparative genetics data for B. napus and Arabidopsis thaliana, an allocation of SNP markers to linkage groups in rapeseed was achieved: a unique location was determined for seven gene sequences; two and three possible locations were found for six and four sequences, respectively. The results demonstrate the usefulness of existing genomic resources for SNP discovery in rapeseed.     

Combined sequence and sequence-structure-based methods for analyzing RAAS gene SNPs: a computational approach

Abstract

The renin–angiotensin–aldosterone system (RAAS) plays a key role in the regulation of blood pressure (BP). Mutations on the genes that encode components of the RAAS have played a significant role in genetic susceptibility to hypertension and have been intensively scrutinized. The identification of such probably causal mutations not only provides insight into the RAAS but may also serve as antihypertensive therapeutic targets and diagnostic markers. The methods for analyzing the SNPs from the huge dataset of SNPs, containing both functional and neutral SNPs is challenging by the experimental approach on every SNPs to determine their biological significance. To explore the functional significance of genetic mutation (SNPs), we adopted combined sequence and sequence-structure-based SNP analysis algorithm. Out of 3864 SNPs reported in dbSNP, we found 108 missense SNPs in the coding region and remaining in the non-coding region. In this study, we are reporting only those SNPs in coding region to be deleterious when three or more tools are predicted to be deleterious and which have high RMSD from the native structure. Based on these analyses, we have identified two SNPs of REN gene, eight SNPs of AGT gene, three SNPs of ACE gene, two SNPs of AT1R gene, three SNPs of CYP11B2 gene and three SNPs of CMA1 gene in the coding region were found to be deleterious. Further this type of study will be helpful in reducing the cost and time for identification of potential SNP and also helpful in selecting potential SNP for experimental study out of SNP pool.     

Identification of 142 single nucleotide polymorphisms in 41 candidate genes for rheumatoid arthritis in the Japanese population

Single-nucleotide polymorphisms (SNPs) can make an important contribution to our understanding of genetic backgrounds that may influence medical conditions and ethnic diversity. We undertook a systematic survey of genomic DNA for SNPs located not only in coding sequences but also in non-coding regions (e.g., introns and 5' flanking regions) of selected genes. Using DNA samples from 48 Japanese patients with rheumatoid arthritis (RA) as templates, we surveyed 41 genes that represent candidates for RA, screening a total of 104 kb of DNA (30 kb of coding sequences and 74 kb of non-coding DNA). Within this 104 kb of genomic sequences we identified 163 polymorphisms (1 per 638 bases on average), of which 142 were single-nucleotide substitutions and the remainder, insertions or deletions. Of the coding SNPs, 52% were non-synonymous substitutions, and non-conservative amino acid changes were observed in a quarter of those. Sixty-nine polymorphisms showed high frequencies for minor alleles (more than 15%) and 20 revealed low frequencies (<5%). Our results indicated a greater average distance between SNPs than others have reported, but this disparity may reflect the type of genes surveyed and/or the relative ethnic homogeneity of our test population.     

A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

We describe a new computer program, SnpEff, for rapidly categorizing the effects of variants in genome sequences. Once a genome is sequenced, SnpEff annotates variants based on their genomic locations and predicts coding effects. Annotated genomic locations include intronic, untranslated region, upstream, downstream, splice site, or intergenic regions. Coding effects such as synonymous or non-synonymous amino acid replacement, start codon gains or losses, stop codon gains or losses, or frame shifts can be predicted. Here the use of SnpEff is illustrated by annotating ~356,660 candidate SNPs in ~117 Mb unique sequences, representing a substitution rate of ~1/305 nucleotides, between the Drosophila melanogaster w(1118); iso-2; iso-3 strain and the reference y(1); cn(1) bw(1) sp(1) strain. We show that ~15,842 SNPs are synonymous and ~4,467 SNPs are non-synonymous (N/S ~0.28). The remaining SNPs are in other categories, such as stop codon gains (38 SNPs), stop codon losses (8 SNPs), and start codon gains (297 SNPs) in the 5'UTR. We found, as expected, that the SNP frequency is proportional to the recombination frequency (i.e., highest in the middle of chromosome arms). We also found that start-gain or stop-lost SNPs in Drosophila melanogaster often result in additions of N-terminal or C-terminal amino acids that are conserved in other Drosophila species. It appears that the 5' and 3' UTRs are reservoirs for genetic variations that changes the termini of proteins during evolution of the Drosophila genus. As genome sequencing is becoming inexpensive and routine, SnpEff enables rapid analyses of whole-genome sequencing data to be performed by an individual laboratory.     

Novel single nucleotide polymorphism identification in interleukin-6 gene of Pakistani sheep

The present study aimed to identify single-nucleotide polymorphism (SNP) in coding and non-coding regions of interleukin-6 (IL-6) gene of Pakistani sheep. The IL-6 gene of 205 animals from nine sheep breeds were sequenced for screening of SNP. Characterizing the IL-6 gene revealed thirteen SNP sites within the intronic region of IL-6 gene. The novel SNPs found in the present study can serve as genetic marker for association studies with susceptibility/resistance to parasite infection in sheep. This is first report of SNP polymorphism of IL-6 gene of Pakistani sheep.     

Mammalian keratin gene families: organisation of genes coding for the B2 high-sulphur proteins of sheep wool.

We have isolated two genomic clones containing three B2 high-sulphur keratin genes from a sheep genomic library constructed in Charon 4A. These genes do not contain intervening sequences. Two genes, encoding the B2A and B2D proteins are closely linked in the genome, being separated by 1.9 kb, and are transcribed in the same direction. Although there is extensive sequence conservation in the 5' non-coding and coding regions, the 3' non-coding regions diverge both in length and sequence. Within the 5' non-coding region adjacent to the initiating AUG there is a highly conserved 18 bp sequence which is also present in another gene coding for a member of a different, unrelated high-sulphur keratin family. In the B2A-B2D intergene region, tightly linked to the B2D gene, there is a putative, divergently transcribed gene.     

Development of a set of SNP markers present in expressed genes of the apple

Molecular markers associated with gene coding regions are useful tools for bridging functional and structural genomics. Due to their high abundance in plant genomes, single nucleotide polymorphisms (SNPs) are present within virtually all genomic regions, including most coding sequences. The objective of this study was to develop a set of SNPs for the apple by taking advantage of the wealth of genomics resources available for the apple, including a large collection of expressed sequenced tags (ESTs). Using bioinformatics tools, a search for SNPs within an EST database of approximately 350,000 sequences developed from a variety of apple accessions was conducted. This resulted in the identification of a total of 71,482 putative SNPs. As the apple genome is reported to be an ancient polyploid, attempts were made to verify whether those SNPs detected in silico were attributable either to allelic polymorphisms or to gene duplication or paralogous or homeologous sequence variations. To this end, a set of 464 PCR primer pairs was designed, PCR was amplified using two subsets of plants, and the PCR products were sequenced. The SNPs retrieved from these sequences were then mapped onto apple genetic maps, including a newly constructed map of a Royal Gala x A689-24 cross and a Malling 9 x Robusta 5, map using a bin mapping strategy. The SNP genotyping was performed using the high-resolution melting (HRM) technique. A total of 93 new markers containing 210 coding SNPs were successfully mapped. This new set of SNP markers for the apple offers new opportunities for understanding the genetic control of important horticultural traits using quantitative trait loci (QTL) or linkage disequilibrium analysis. These also serve as useful markers for aligning physical and genetic maps, and as potential transferable markers across the Rosaceae family.     

The structure of a single unit of ribosomal RNA gene (rDNA) including intergenic subrepeats in the Australian bulldog ant Myrmecia croslandi (Hymenoptera: Formicidae)

A complete single unit of a ribosomal RNA gene (rDNA) of M. croslandi was sequenced. The ends of the 18S, 5.8S and 28S rRNA genes were determined by using the sequences of D. melanogaster rDNAs as references. Each of the tandemly repeated rDNA units consists of coding and non-coding regions whose arrangement is the same as that of D. melanogaster rDNA. The intergenic spacer (IGS) contains, as in other species, a region with subrepeats, of which the sequences are different from those previously reported in other insect species. The length of IGSs was estimated to be 7-12 kb by genomic Southern hybridization, showing that an rDNA repeating unit of M. croslandi is 14-19 kb-long. The sequences of the coding regions are highly conserved, whereas IGS and ITS (internal transcribed spacer) sequences are not. We obtained clones with insertions of various sizes of R2 elements, the target sequence of which was found in the 28S rRNA coding region. A short segment in the IGS that follows the 3' end of the 28S rRNA gene was predicted to form a secondary structure with long stems.     

Identification of 187 single nucleotide polymorphisms (SNPs) among 41 candidate genes for ischemic heart disease in the Japanese population

To investigate whether common variants in the human genetic background are associated with pathogenesis of ischemic heart diseases, we systematically surveyed 41 possible candidate genes for single-nucleotide polymorphisms (SNPs) by directly sequencing 96 independent alleles at each locus, derived from 48 unrelated Japanese patients with myocardial infarction, including 25.8 kb 5' flanking regions, 56.8 kb exonic and 35.4 kb intronic sequences, and 1.8 kb 3' flanking regions. In this genomic DNA of nearly 120 kb, we identified 187 SNPs: 55 in 5' flanking regions, seven in 5' untranslated regions (UTRs), 52 in coding elements, 64 in introns, eight in 3' UTRs, and one in a 3' flanking region. Among the 52 coding SNPs, 26 were non-synonymous changes. Allelic frequencies of some of the polymorphisms were significantly different from those reported in European populations. For example, the Q506R substitution in the coagulation factor V gene, the so-called "Leiden mutation", has a reported frequency of 2.3% in Europeans, but we detected the Leiden mutation in none of the Japanese genomes that we investigated. The allelic frequencies of the -33A>G SNP in the thrombomodulin gene were also very different; this allele occurred at a 12% frequency in the Japanese patients that we examined, although it had been detected in none of 82 Caucasians reported previously. These data support the hypothesis that some SNPs are specific to particular ethnic groups.     

Structure of the rat prolactin gene

The organization and sequence of the rat preprolactin gene has been investigated. Analysis of two different plasmids containing pituitary cDNA inserts has provided the complete 681-nucleotide coding sequence of preprolactin as well as 17 nucleotides preceding the initiation codon and 90 nucleotides following the termination codon. Digestion of rat chromosomal DNA with the restriction endonuclease Eco RI followed by size fractionation and hybridization to a labeled prolactin cDNA probe has demonstrated that prolactin genomic sequences are located on 6.0-, 3.9-, and 2.9-kilobase fragments. The 6.0- and 3.9-kilobase fragments were isolated from a library of cloned rat DNA fragments. The sequence of more than 1800 nucleotides of the cloned DNA has been determined. The sequenced region contains coding regions of 180 and 189 nucleotides which specify the COOH-terminal 123 amino acids of the 227-amino-acid sequence of rat preprolactin. These coding regions are separated by an intervening sequence of 597 nucleotides. At least one other large intervening sequence separates this region from the region coding for the NH2-terminal portion of preprolactin. Hybridization experiments suggested that the intervening sequences of the rat prolactin gene contain DNA sequences which are repeated elsewhere in the rat genome.     

Genome-wide single-nucleotide polymorphism map for Candida albicans

Single-nucleotide polymorphisms (SNPs) are essential tools for studying a variety of organismal properties and processes, such as recombination, chromosomal dynamics, and genome rearrangement. This paper describes the development of a genome-wide SNP map for Candida albicans to study mitotic recombination and chromosome loss. C. albicans is a diploid yeast which propagates primarily by clonal mitotic division. It is the leading fungal pathogen that causes infections in humans, ranging from mild superficial lesions in healthy individuals to severe, life-threatening diseases in patients with suppressed immune systems. The SNP map contains 150 marker sequences comprising 561 SNPs and 9 insertions-deletions. Of the 561 SNPs, 437 were transition events while 126 were transversion events, yielding a transition-to-transversion ratio of 3:1, as expected for a neutral accumulation of mutations. The average SNP frequency for our data set was 1 SNP per 83 bp. The map has one marker placed every 111 kb, on average, across the 16-Mb genome. For marker sequences located partially or completely within coding regions, most contained one or more nonsynonymous substitutions. Using the SNP markers, we identified a loss of heterozygosity over large chromosomal fragments in strains of C. albicans that are frequently used for gene manipulation experiments. The SNP map will be useful for understanding the role of heterozygosity and genome rearrangement in the response of C. albicans to host environments.     

cDNA of cattle alpha S1-casein: cloning and nucleotide sequence

The nucleotide sequence of recombinant plasmids representing a full-size cDNA of cow alpha s1-casein was investigated. The corresponding mRNA consists of 1133 nucleotides except for poly(A) and includes 642 nucleotides of the coding region, 63 nucleotides of 5'- and 428 nucleotides of the 3'-noncoding regions. A comparative analysis of nucleotide sequences of cow alpha s1-casein and guinea pig B-casein showed that the homology in the 5'-nontranslatable region is 90.5%, that of a precasein single peptide is 82.22%, while that of the major polypeptide in the coding region is 64% without taking into account the blank spaces. The homology is higher in the 3'-noncoding region than in the coding region and makes up to 72%. The data obtained testify to the high degree of conservatism of sequences in casein mRNA noncoding regions as well as to functional and regulatory role of these sequences in gene expression of caseins.     

The human thyroglobulin gene contains two 15-17 kb introns near its 3''-end.

We have cloned overlapping segments of the human thyroglobulin gene from a genomic cosmid library. Restriction mapping and electron microscopy show that a region of 38 kb at or near the 3'-end of this gene encodes only 850 nucleotides or 10% of the messenger RNA (mRNA) sequence. The region contains five exons of 130-210 nucleotides, split by introns of 1 to 15-17 kb. This represents the lowest ratio of coding to non-coding DNA (2.2%) found thus far in any eukaryotic gene. Blot hybridization under non-stringent conditions shows the presence of only one copy of this gene in the human genome and the absence of other closely related sequences.     

Polymorphism attribution of cSNPs in cancer-related genes located in loss regions with a high frequency of HCC between HBV and health groups

Cancer-related genes harbored in the loss regions containing a high frequency of hepatocellular carcinoma (HCC) were selected.Related information was gathered and the coding single nucleotide polymorphism (cSNP) sequences were obtained from the single nucleotide polymorphism (SNP) database.The appropriate primers and oligonucleotide probes were then designed in accordance with the SNP sites,and subsequently,the gene chips for detecting SNPs were constructed.Genomic DNA was extracted from blood samples of healthy controls and from patients with HBV infection.The sequences,including the SNPs,were amplified via polymerase chain reaction (PCR) and labeled using digoxigenin deoxyuridine tri-phosphate (Dig-dUTP).The labeled products were then hybridized with the SNP chips.Results confirmed that the differences in allele frequencies of three SNPs EGFL3 (rs947345),Caspase9 (rs2308950),and E2F2 (rs3218171) were distinct between HBV-infected patients and controls,suggesting that these SNPs ocuring in high frequency in HBV-infected individuals may be associated with susceptibility to HCC.