A universal method for automated gene mapping

Small insertions or deletions (InDels) constitute a ubiquituous class of sequence polymorphisms found in eukaryotic genomes. Here, we present an automated high-throughput genotyping method that relies on the detection of fragment-length polymorphisms (FLPs) caused by InDels. The protocol utilizes standard sequencers and genotyping software. We have established genome-wide FLP maps for both Caenorhabditis elegans and Drosophila melanogaster that facilitate genetic mapping with a minimum of manual input and at comparatively low cost.     

Development of PCR-based SNP markers for rice blast resistance genes at the<Emphasis Type="Italic"> Piz</Emphasis> locus

We assessed the utility of single-nucleotide polymorphisms (SNPs) and small insertion/deletion polymorphisms (InDels) as DNA markers in genetic analysis and breeding of rice. Toward this end, we surveyed SNPs and InDels in the chromosomal region containing the Piz and Piz-t rice blast resistance genes and developed PCR-based markers for typing the SNPs. Analysis of sequences from a blast-susceptible Japanese cultivar and two cultivars each containing one of these genes revealed that SNPs are abundant in the Piz and Piz-t regions (on average, one SNP every 248 bp), but the number of InDels was much lower. The dense distribution of SNPs facilitated the generation of SNP markers in the vicinity of the genes. For typing these SNPs, we used a modified allele-specific PCR method. Of the 49 candidate allele-specific markers, 33 unambiguously and reproducibly discriminated between the two alleles. We used the markers for mapping the Piz and Piz-t genes and evaluating the size of DNA segments introgressed from the Piz donor cultivar in Japanese near-isogenic lines containing Piz. Our findings suggest that, because of its ability to generate numerous markers within a target region and its simplicity in assaying genotypes, SNP genotyping with allele-specific PCR is a valuable tool for gene mapping, map-based cloning, and marker-assisted selection in crops, especially rice.Communicated by D.J. Mackill     

Sequence polymorphisms in porcine homologs of murine coat colour‐related genes

Herein, we report the variability among 57 porcine homologs of murine coat colour‐related genes. We identified single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within 44 expressed gene sequences by aligning eight pig complementary DNA (cDNA) samples. The sequence alignment revealed a total of 485 SNPs and 15 InDels. The polymorphisms were then validated by performing matrix‐assisted laser desorption/ionization time‐of‐flight mass spectrometry (MALDI‐TOF MS) with reference DNA samples obtained from 384 porcine individuals. Of the 384 individuals, three parents of the experimental F₂ family were included to detect polymorphisms between them for linkage mapping. We also genotyped previously reported polymorphisms of 12 genes, and one SNP each in three genes that were detected by performing a BLAST search of the Trace database. A total of 211 SNPs and three InDels were successfully genotyped from our porcine DNA panel. We detected SNPs in 33 of the 44 genes among the parents of an experimental F₂ family and then constructed a linkage map of the 33 genes for this family. The linkage assignment of each gene to the porcine chromosomes was consistent with the location of the BAC clone in the porcine genome and the corresponding gene sequence. We confirmed complete substitutions of EDNRB and MLPH in the Jinhua and Clawn miniature breeds, respectively. Furthermore, we identified polymorphic alleles exclusive to each pig group: 13 for Jinhua, two for Duroc, three for Meishan, four for the Japanese wild boar, one for the Clawn miniature pig and four for the Potbelly pig.     

Genomewide discovery of DNA polymorphisms in rice cultivars with contrasting drought and salinity stress response and their functional relevance

Next‐generation sequencing technologies provide opportunities to understand the genetic basis of phenotypic differences, such as abiotic stress response, even in the closely related cultivars via identification of large number of DNA polymorphisms. We performed whole‐genome resequencing of three rice cultivars with contrasting responses to drought and salinity stress (sensitive IR64, drought‐tolerant Nagina 22 and salinity‐tolerant Pokkali). More than 356 million 90‐bp paired‐end reads were generated, which provided about 85% coverage of the rice genome. Applying stringent parameters, we identified a total of 1 784 583 nonredundant single‐nucleotide polymorphisms (SNPs) and 154 275 InDels between reference (Nipponbare) and the three resequenced cultivars. We detected 401 683 and 662 509 SNPs between IR64 and Pokkali, and IR64 and N22 cultivars, respectively. The distribution of DNA polymorphisms was found to be uneven across and within the rice chromosomes. One‐fourth of the SNPs and InDels were detected in genic regions, and about 3.5% of the total SNPs resulted in nonsynonymous changes. Large‐effect SNPs and InDels, which affect the integrity of the encoded protein, were also identified. Further, we identified DNA polymorphisms present in the differentially expressed genes within the known quantitative trait loci. Among these, a total of 548 SNPs in 232 genes, located in the conserved functional domains, were identified. The data presented in this study provide functional markers and promising target genes for salinity and drought tolerance and present a valuable resource for high‐throughput genotyping and molecular breeding for abiotic stress traits in rice.     

Discovery of genome-wide DNA polymorphisms in a landrace cultivar of Japonica rice by whole-genome sequencing

Molecular breeding approaches are of growing importance to crop improvement. However, closely related cultivars generally used for crossing material lack sufficient known DNA polymorphisms due to their genetic relatedness. Next-generation sequencing allows the identification of a massive number of DNA polymorphisms such as single nucleotide polymorphisms (SNPs) and insertions-deletions (InDels) between highly homologous genomes. Using this technology, we performed whole-genome sequencing of a landrace of japonica rice, Omachi, which is used for sake brewing and is an important source for modern cultivars. A total of 229 million reads, each comprising 75 nucleotides of the Omachi genome, was generated with 45-fold coverage and uniquely mapped to 89.7% of the Nipponbare genome, a closely related cultivar. We identified 132,462 SNPs, 16,448 insertions and 19,318 deletions between the Omachi and Nipponbare genomes. An SNP array was designed to validate 731 selected SNPs, resulting in validation rates of 95 and 88% for the Omachi and Nipponbare genomes, respectively. Among the 577 SNPs validated in both genomes, 532 are entirely new SNP markers not previously reported between related rice cultivars. We also validated InDels on a part of chromosome 2 as DNA markers and successfully genotyped five japonica rice cultivars. Our results present the methodology and extensive data on SNPs and InDels available for whole-genome genotyping and marker-assisted breeding. The polymorphism information between Omachi and Nipponbare is available at NGRC_Rice_Omachi (http://www.nodai-genome.org/oryza_sativa_en.html).     

Development of InDel markers for Brassica rapa based on whole-genome re-sequencing

Genome-wide detection of short insertion/deletion length polymorphisms (InDels, <5 bp) in Brassica rapa (named the A genome) was performed by comparing whole-genome re-sequencing data from two B. rapa accessions, L144 and Z16, to the reference genome sequence of Chiifu-401-42. In total, we identified 108,558 InDel polymorphisms between Chiifu-401-42 and L144, 26,795 InDels between Z16 and Chiifu-401-42, and 26,693 InDels between L144 and Z16. From these, 639 InDel polymorphisms of 3–5 bp in length between L144 and Z16 were selected for experimental validation; 491 (77 %) yielded single PCR fragments and showed polymorphisms, 7 (1 %) did not amplify a product, and 141 (22 %) showed no polymorphism. For further validation of these intra-specific InDel polymorphisms, 503 candidates, randomly selected from the 639 InDels, were screened across seven accessions representing different B. rapa cultivar groups. Of these assayed markers, 387 (77 %) were polymorphic, 111 (22 %) were not polymorphic and 5 (1 %) did not amplify a PCR product. Furthermore, we randomly selected 518 InDel markers to validate their polymorphism in B. napus (the AC genome) and B. juncea (the AB genome), of which more than 90 % amplified a PCR product; 132 (25 %) showed polymorphism between the two B. napus accessions and 41 (8 %) between the two B. juncea accessions. This set of novel PCR-based InDel markers will be a valuable resource for genetic studies and breeding programs in B. rapa.     

Selection of 29 highly informative InDel markers for human identification and paternity analysis in Chinese Han population by the SNPlex genotyping system

The interest of forensic researchers in single nucleotide polymorphism (SNP) has been attracted because of its potential advantages, such as low mutation rates, amenable to high-throughput automated platform and the improved application in the analysis of degraded samples. In this paper, 29 highly informative insertion/deletion (InDel, a special kind of SNP) markers were selected from the dbSNP () according to the given criteria. 109 unrelated Chinese Han subjects were genotyped for the 29 InDels with SNPlex genotyping system. The allele frequency data revealed that the combined power of discrimination for the 29 InDel markers was 0.999999999990867 and the combined probability of paternity exclusion (PE) was 0.9930. Sensitivity studies were performed to evaluate the flexibility of the SNPlex genotyping system on the set of 29 InDels. Highly reproducible results could be obtained with 40–100 ng genomic DNA and the proportion of total allele drop-in was significantly increased when the amount of DNA added to PCR was lower than 35 ng. These results suggested that the set of 29 InDels was useful in paternity analysis or human identification in the future.     

Identification and characterization of nucleotide variations in the genome of Ziziphus jujuba (Rhamnaceae) by next generation sequencing

In this study, single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) in the genome of Ziziphus jujuba were identified using sequences generated by the Roche 454 GS-FLX sequencer. A total of, 573,141 reads were produced with an average read length of 360 bp. After quality control, 258,754 of the filtered reads were assembled into 23,864 contigs, and 293,458 remained as singletons. Using the contig assemblies as a reference, 17,160 SNPs and 478 InDels were identified. Among the SNPs, transitions occurred three times more frequently than transversions. In transitions, the number of C/T and G/A transitions was similar. Among the transversions, A/T was the most abundant, and C/G was much rarer than any of the other types of transversions, accounting for only about half the numbers of A/C, A/T and G/T transversions. For the InDels, mononucleotide changes amounted to 64.4 % of the total number of InDels. In general, the frequency of detected InDels decreased as the length of the InDels increased. This study provides valuable marker resources for future genetic studies of Ziziphus spp.     

Development of Species-Specific InDel Markers in Citrus

Citrus taxonomy is complex owing to the existence of a wide range of species: Poncirus is used mainly for rootstock; Fortunella produces small fruit and edible pericarp; and Citrus comprises the most widespread fruit crop species worldwide. Rapidly increasing genome resources from different citrus species facilitate the development of convenient and genome-wide molecular markers that can be applied to both inter- and intra-species analyses. In this study, by comparing the genome sequences of four citrus species, a set of 1958 InDels were identified and 453 candidate InDels were converted into PCR-based markers. Among these candidate InDels, 268 (65%) exhibited length polymorphisms from 30 bp to 200 bp when applied to seven species from the genera Poncirus, Fortunella and Citrus. Seven InDel markers exhibited high intraspecific polymorphisms in a natural pummelo population. The results showed that the InDel markers are effective for both inter- and intra-specific variation and identification analyses. These InDel markers are expected to be applied to germplasm identification, phylogenetic analysis, genetic diversity evaluation and marker-assisted breeding in citrus.     

珍稀濒危飘带兜兰叶绿体全基因组种内变异研究(附录)

飘带兜兰(Paphiopedilum parishii)分布范围狭窄,仅在中国、缅甸、泰国以及老挝有少量分布。近年来,因生境破坏和人为滥采而导致飘带兜兰野生种群极度缩减。为开发种内多态性的分子标记用于保护生物学研究,该研究对飘带兜兰4个野生个体经测序、组装、注释获得的叶绿体基因组序列,与已公布的飘带兜兰2个个体的叶绿体全基因组序列进行比对,分析飘带兜兰叶绿体基因组的种内差异。结果表明:(1)飘带兜兰叶绿体基因组具有典型被子植物叶绿体基因组环状四分体结构,基因组长度为154 403～154 809 bp,共编码129个基因,包括78个蛋白质编码基因、39个tRNA基因、8个rRNA基因,以及4个假基因。(2)在飘带兜兰6个个体叶绿体基因组中检测到103～107个SSRs(simple sequence repeats)位点,其中21个SSR位点具有多态性。此外,在6个个体叶绿体基因组中还检测到60个长序列重复,包括17～21个正向重复、18～29个反向重复、9～16个回文重复、4～9个互补重复。(3)通过比较6个个体叶绿体基因组序列的核苷酸多样性,共发现70处变异,包括10个SNPs(single nucleotide polymorphism)、60个插入缺失(InDels)。其中,有3个SNP位点发生了非同义替换,导致编码功能基因的氨基酸发生改变; 19个插入缺失多态性较高,具有开发为分子标记的潜力。(4)通过计算核苷酸多样性值(P_i)共发现8个有变异的区域,P_i值为0～0.006 32,其中变异度较大的是rps3-rpl22、trnL-UAC-rpl32、rpoB-trnC-GCA以及ycf4,这些高变区可开发为分子标记用于评估飘带兜兰遗传多样性。(5)系统发生分析结果表明,飘带兜兰6个个体叶绿体基因组序列聚在一起,与长瓣兜兰互为姐妹群。综上表明,飘带兜兰叶绿体基因组的SSRs、长序列重复、SNPs、InDels以及核苷酸序列呈现了足够的种内多样性,可开发成分子标记用于该种的系统演化及保护生物学研究。     

Sequence variation in 3′UTR region of crtRB1 gene and its effect on β-carotene accumulation in maize kernel

β-carotene fortification of maize has emerged as a potential, long-term and sustainable approach to alleviate vitamin A deficiency in humans. Among the several genes involved in the carotenoid biosynthetic pathway, the 543 bp allele at crtRB1 3′TE (Transposable Element) gene (allele 1, without insertion) is associated with higher β-carotene accumulation. Estimation of β-carotene through high performance liquid chromatography showed that the CIMMYT genotypes with allele 1 had high kernel β-carotene content whereas the Indian inbreds with the same allele had low β-carotene content. To know the reason for this variation, allele 1 of crtRB1 3′TE gene was sequenced from a set of 11 diverse maize inbreds collected from CIMMYT and Indian germplasm. The sequence data of the allele 1 revealed the presence of 13 single nucleotide polymorphisms (SNPs) and 7 insertions and deletions (InDels). Exonic region had two SNPs, intronic region had one SNP and one InDel, whereas 3′-untranslated region (UTR) region of the gene showed 10 SNPs and 6 InDels. Among the several SNPs and InDels, SNP4, SNP13, InDel6 and InDel7 identified in the 3′-UTR region clearly differentiated the high and the low β-carotene genotypes. These 3′-UTR polymorphisms in allele 1 of the crtRB1 3′TE gene could be associated with the variation in kernel β-carotene accumulation by regulating the translation and stability of the mRNA. The SNPs and the InDels associated with higher level of β-carotene will be used as a gene-based marker(s) in selection of genotypes and to develop biofortified maize hybrids to alleviate vitamin A deficiency in humans.     

Whole genome resequencing in tomato reveals variation associated with introgression and breeding events

Background

One of the goals of genomics is to identify the genetic loci responsible for variation in phenotypic traits. The completion of the tomato genome sequence and recent advances in DNA sequencing technology allow for in-depth characterization of genetic variation present in the tomato genome. Like many self-pollinated crops, cultivated tomato accessions show a low molecular but high phenotypic diversity. Here we describe the whole-genome resequencing of eight accessions (four cherry-type and four large fruited lines) chosen to represent a large range of intra-specific variability and the identification and annotation of novel polymorphisms.

Results

The eight genomes were sequenced using the GAII Illumina platform. Comparison of the sequences with the reference genome yielded more than 4 million single nucleotide polymorphisms (SNPs). This number varied from 80,000 to 1.5 million according to the accessions. Almost 128,000 InDels were detected. The distribution of SNPs and InDels across and within chromosomes was highly heterogeneous revealing introgressions from wild species and the mosaic structure of the genomes of the cherry tomato accessions. In-depth annotation of the polymorphisms identified more than 16,000 unique non-synonymous SNPs. In addition 1,686 putative copy-number variations (CNVs) were identified.

Conclusions

This study represents the first whole genome resequencing experiment in cultivated tomato. Substantial genetic differences exist between the sequenced tomato accessions and the reference sequence. The heterogeneous distribution of the polymorphisms may be related to introgressions that occurred during domestication or breeding. The annotated SNPs, InDels and CNVs identified in this resequencing study will serve as useful genetic tools, and as candidate polymorphisms in the search for phenotype-altering DNA variations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-791) contains supplementary material, which is available to authorized users.     

Genome-wide discovery of DNA polymorphism in Brassica rapa

Single nucleotide polymorphisms (SNPs) and/or insertion/deletions (InDels) are frequent sequence variations in the plant genome, which can be developed as molecular markers for genetic studies on crop improvement. The ongoing Brassica rapa genome sequencing project has generated vast amounts of sequence data useful in genetic research. Here, we report a genome-wide survey of DNA polymorphisms in the B. rapa genome based on the 557 bacterial artificial clone sequences of B. rapa ssp. pekinensis cv. Chiifu. We identified and characterized 21,311 SNPs and 6,753 InDels in the gene space of the B. rapa genome by re-sequencing 1,398 sequence-tagged sites (STSs) in eight genotypes. Comparison of our findings with a B. rapa genetic linkage map confirmed that STS loci were distributed randomly over the B. rapa whole genome. In the 1.4 Mb of aligned sequences, mean nucleotide polymorphism and diversity were θ = 0.00890 and π = 0.00917, respectively. Additionally, the nucleotide diversity in introns was almost three times greater than that in exons, and the frequency of observed InDel was almost 17 times higher in introns than in exons. Information regarding SNPs/InDels obtained here will provide an important resource for genetic studies and breeding programs of B. rapa.     

Design and validation of a high-throughput assay to detect codon 146 polymorphisms in the caprine prion protein gene

In sheep, scrapie susceptibility is so strongly associated with single nucleotide polymorphisms (SNPs) in the gene encoding the prion protein (PrP) that this linkage constitutes the basis for selective breeding strategies directed toward controlling the disease. For goats, in contrast, the association between scrapie susceptibility/resistance and variations in the PrP gene is far weaker, with only a few identified SNPs showing an influence on scrapie susceptibility. A recent survey of PrP genotypes in Cypriot goats, however, revealed the existence of a robust association between polymorphisms at codon 146 of the caprine PrP gene and resistance/susceptibility to natural scrapie. Here we describe here a high-throughput assay, based on homogeneous MassExtend technology coupled with mass spectrometry, for genotyping codon 146 of the caprine PrP gene. Our results demonstrate that this assay exhibits high accuracy, reproducibility, and repeatability, thereby making it suitable for large-scale SNP genotyping, as required for scrapie surveillance programs.     

Performance evaluation of an in-house panel containing 59 autosomal InDels for forensic identification in Chinese Hui and Mongolian groups

In recent years, a novel multiplex system containing two mini-short tandem repeats, 59 autosomal InDels, two Y-chromosomal InDels, and the Amelogenin gene with all amplicons less than 200 bp has been constructed and validated by ourselves for forensic degration sample, and its forensic application efficiency has been studied in Chinese some populations. Herein, the population genetic polymorphisms of these loci were investigated in Chinese Hui (n = 249) and Mongolian (n = 222) ethnic groups using direct multiplex amplification and capillary electrophoresis platform. The forensic identification efficiencies of this self-developed system were further evaluated in these two groups. And the results showed that the values of the combined power of discrimination were 0.9999999999999999999999999999006 (Hui) and 0.999999999999999999999999999738 (Mongolian), respectively. Moreover, the combined power of exclusion values were 0.99999817 (Hui) and 0.99999779 (Mongolian). The 59 autosomal InDels used in this study exhibited high forensic identification efficiencies in 10 East Asian populations, which was also expected to be a new powerful tool for identifying degraded biological materials in East Asian populations.     

Development of Genetic Markers in Eucalyptus Species by Target Enrichment and Exome Sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.     

Fluorescence resonance energy transfer (FRET) using ssDNA binding fluorescent dye

There is a need for simple and inexpensive methods for genotyping single nucleotide polymorphisms (SNPs) and short insertion/deletion variations (InDels). In this work, I demonstrate that a single-stranded DNA (ssDNA) binding dye can be used as a donor fluorophore for fluorescence resonance energy transfer (FRET). The method presented is a homogenous assay in which detection is based on the FRET from the fluorescence of the ssDNA dye bound to the unmodified detection primer to the fluorescent nucleotide analog incorporated into this detection primer during cyclic template directed primer extension reaction. Collection of the FRET emission spectrum with a scanning fluorescence spectrophotometer allows powerful data analysis. The fluorescence emission signal is modified by the optical properties of the assay vessel. This seems to be a completely neglected parameter. By proper selection of the optical properties of the assay plate one can improve the detection of the fluorescence emission signal.     

Genome-wide DNA polymorphism and transcriptome analysis of an early-maturing rice mutant

In order to develop a rice population with improved important traits such as flowering time, we developed 2,911 M₂ targeting-induced local lesions in genomes (TILLING) lines by irradiating rice seeds with γ-rays. In all, 15 M₃ lines were obtained from 3 different M₂ lines that exhibited an early-maturing phenotype: these plants matured approximately 25 days faster than wild-type (WT) plants. To identify genome-wide DNA polymorphisms, we performed whole-genome resequencing of both the plant types, i.e., WT and early-maturing TILLING 1 (EMT1), and obtained mapped reads of 118,488,245 bp (99.53 %) and 128,489,860 bp (99.72 %), respectively; Nipponbare was used as the reference genome. We obtained 63,648 and 147,728 single nucleotide polymorphisms (SNPs) and 33,474 and 31,082 insertions and deletions (InDels) for the WT and EMT1, respectively. Interestingly, there was a higher number of SNPs (2.6-fold) and slightly lower number of InDels (0.9-fold) in EMT1 than in WT. The expression of at least 202 structurally altered genes was changed in EMT1, and functional enrichment analysis of these genes revealed that their molecular functions were related to flower development. These results might provide a critical insight into the regulatory pathways of rice flowering.