Genome-wide single nucleotide polymorphism and Insertion-Deletion discovery through next-generation sequencing of reduced representation libraries in common bean

Single nucleotide polymorphisms (SNPs) and insertions-deletions (InDels) are valuable molecular markers for genomics and genetics studies and molecular breeding. The advent of next-generation sequencing techniques has enabled researchers to approach high-throughput and cost-effective SNP and InDel discovery on a genomic scale. In this report, 36 common bean genotypes grown in Canada were used to construct reduced representation libraries for next-generation sequencing. Using 76 million sequence reads generated by the Illumina HiSeq 2000 Sequencing System, we identified a total of 43,698 putative SNPs and 1,267 putative InDels. Of the SNPs, 43,504 were bi-allelic and 194 were tri-allelic, and the InDels comprised 574 insertions and 693 deletions. The putative bi-allelic SNPs were distributed across all 11 chromosomes with the highest number of SNPs observed in chromosome 2 (4,788), and the lowest in chromosome 10 (2,941). With the aid of the recent release of the first chromosome-scale version of Phaseolus vulgaris, 24,907 bi-allelic SNPs, 79 tri-allelic SNPs, 315 insertions, and 377 deletions were located in 8,758, 77, 273, and 364 genes, respectively. Among these 24,907 bi-allelic SNPs, 7,168 nonsynonymous bi-allelic SNPs were identified within 36 common bean genotypes that were located in 4,303 genes. A total of 113 putative SNPs were randomly chosen for validation using high-resolution melt analysis. Of the 113 candidate SNPs, 105 (92.9 %) contained the predicted SNPs.     

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     

First-generation SNP/InDel markers tagging loci for pathogen resistance in the potato genome

A panel of 17 tetraploid and 11 diploid potato genotypes was screened by comparative sequence analysis of polymerase chain reaction (PCR) products for single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (InDels), in regions of the potato genome where genes for qualitative and/or quantitative resistance to different pathogens have been localized. Most SNP and InDel markers were derived from bacterial artificial chromosome (BAC) insertions that contain sequences similar to the family of plant genes for pathogen resistance having nucleotide-binding-site and leucine-rich-repeat domains (NBS-LRR-type genes). Forty-four such NBS-LRR-type genes containing BAC-insertions were mapped to 14 loci, which tag most known resistance quantitative trait loci (QTL) in potato. Resistance QTL not linked to known resistance-gene-like (RGL) sequences were tagged with other markers. In total, 78 genomic DNA fragments with an overall length of 31 kb were comparatively sequenced in the panel of 28 genotypes. 1498 SNPs and 127 InDels were identified, which corresponded, on average, to one SNP every 21 base pairs and one InDel every 243 base pairs. The nucleotide diversity of the tetraploid genotypes (pi = 0.72 x 10(-3)) was lower when compared with diploid genotypes (pi = 2.31 x 10(-3)). RGL sequences showed higher nucleotide diversity when compared with other sequences, suggesting evolution by divergent selection. Information on sequences, sequence similarities, SNPs and InDels is provided in a database that can be queried via the Internet.     

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.     

A catalog of neutral and deleterious polymorphism in yeast

The abundance and identity of functional variation segregating in natural populations is paramount to dissecting the molecular basis of quantitative traits as well as human genetic diseases. Genome sequencing of multiple organisms of the same species provides an efficient means of cataloging rearrangements, insertion, or deletion polymorphisms (InDels) and single-nucleotide polymorphisms (SNPs). While inbreeding depression and heterosis imply that a substantial amount of polymorphism is deleterious, distinguishing deleterious from neutral polymorphism remains a significant challenge. To identify deleterious and neutral DNA sequence variation within Saccharomyces cerevisiae, we sequenced the genome of a vineyard and oak tree strain and compared them to a reference genome. Among these three strains, 6% of the genome is variable, mostly attributable to variation in genome content that results from large InDels. Out of the 88,000 polymorphisms identified, 93% are SNPs and a small but significant fraction can be attributed to recent interspecific introgression and ectopic gene conversion. In comparison to the reference genome, there is substantial evidence for functional variation in gene content and structure that results from large InDels, frame-shifts, and polymorphic start and stop codons. Comparison of polymorphism to divergence reveals scant evidence for positive selection but an abundance of evidence for deleterious SNPs. We estimate that 12% of coding and 7% of noncoding SNPs are deleterious. Based on divergence among 11 yeast species, we identified 1,666 nonsynonymous SNPs that disrupt conserved amino acids and 1,863 noncoding SNPs that disrupt conserved noncoding motifs. The deleterious coding SNPs include those known to affect quantitative traits, and a subset of the deleterious noncoding SNPs occurs in the promoters of genes that show allele-specific expression, implying that some cis-regulatory SNPs are deleterious. Our results show that the genome sequences of both closely and distantly related species provide a means of identifying deleterious polymorphisms that disrupt functionally conserved coding and noncoding sequences.     

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.     

Deploying QTL-seq for rapid delineation of a potential candidate gene underlying major trait-associated QTL in chickpea

A rapid high-resolution genome-wide strategy for molecular mapping of major QTL(s)/gene(s) regulating important agronomic traits is vital for in-depth dissection of complex quantitative traits and genetic enhancement in chickpea. The present study for the first time employed a NGS-based whole-genome QTL-seq strategy to identify one major genomic region harbouring a robust 100-seed weight QTL using an intra-specific 221 chickpea mapping population (desi cv. ICC 7184 × desi cv. ICC 15061). The QTL-seq-derived major SW QTL (CaqSW1.1) was further validated by single-nucleotide polymorphism (SNP) and simple sequence repeat (SSR) marker-based traditional QTL mapping (47.6% R² at higher LOD >19). This reflects the reliability and efficacy of QTL-seq as a strategy for rapid genome-wide scanning and fine mapping of major trait regulatory QTLs in chickpea. The use of QTL-seq and classical QTL mapping in combination narrowed down the 1.37 Mb (comprising 177 genes) major SW QTL (CaqSW1.1) region into a 35 kb genomic interval on desi chickpea chromosome 1 containing six genes. One coding SNP (G/A)-carrying constitutive photomorphogenic9 (COP9) signalosome complex subunit 8 (CSN8) gene of these exhibited seed-specific expression, including pronounced differential up-/down-regulation in low and high seed weight mapping parents and homozygous individuals during seed development. The coding SNP mined in this potential seed weight-governing candidate CSN8 gene was found to be present exclusively in all cultivated species/genotypes, but not in any wild species/genotypes of primary, secondary and tertiary gene pools. This indicates the effect of strong artificial and/or natural selection pressure on target SW locus during chickpea domestication. The proposed QTL-seq-driven integrated genome-wide strategy has potential to delineate major candidate gene(s) harbouring a robust trait regulatory QTL rapidly with optimal use of resources. This will further assist us to extrapolate the molecular mechanism underlying complex quantitative traits at a genome-wide scale leading to fast-paced marker-assisted genetic improvement in diverse crop plants, including chickpea.     

Molecular Mapping and Characterization of Genes Governing Time to Flowering, Seed Weight, and Plant Height in an Intraspecific Genetic Linkage Map of Chickpea (Cicer arietinum)

Drought is the major constraint to chickpea productivity worldwide. Utilizing early flowering genotypes and larger seed size have been suggested as strategies for breeding in drought zones. Therefore, this study aimed to identify potential markers linked to days-to-flowering, 100-seed weight, and plant height in a chickpea intraspecific F_2:3 population derived from the cross ILC3279 × ICCV2. A closely linked marker (TA117) on linkage group LG3 was identified for the days-to-flowering trait, explaining 33% of the variation. In relation to plant height, a quantitative trait loci (QTL) was located in LG3, close to the Ts5 marker, that explained 29% of phenotypic variation. A QTL for 100-seed weight located in LG4, close to TA176, explained 51% of variation. The identification of a locus linked both to high 100-seed weight and days-to-flowering may account for the correlation observed between these traits in this and other breeding attempts.     

Genome-Wide Analysis of Genetic Variations and the Detection of Rich Variants of NBS-LRR Encoding Genes in Common Wild Rice Lines

Common wild rice (Oryza rufipogon Griff.) is invaluable genetic resource for rice resistance breeding. Whole-genome re-sequencing was conducted to systematically analyze the variations in two new inbred lines (Huaye 3 and Huaye 4) developed from a common wild rice. A total of 4,841,127 SNPs, 1,170,479 InDels, 24,080 structural variations (SVs), and 298 copy number variations (CNVs) were identified in three materials. Approximately 16.24 and 5.64% of the total SNPs and InDels of Huaye 3 and Huaye 4 were located in genic regions, respectively. Together, 12,486 and 15,925 large-effect SNPs, and 12,417 and 14,513 large-effect InDels, which affect the integrity of the encoded protein, were identified in Huaye 3 and Huaye 4, respectively. The distribution map of 194 and 245 NBS-LRR encoding homologs was constructed across 12 rice chromosomes. Further, GO enrichment analysis of the homologs with identical genotype variations in Huaye 3 and Huaye 4 revealed 67, 82, and 58 homologs involved in cell death, response to stress, and both terms, respectively. Comparative analysis displayed that 550 out of 652 SNPs and 129 out of 147 InDels were present in a widely used blast-susceptible rice variety (LTH). Protein-protein interaction analysis revealed a strong interaction between NBS-LRR candidates and several known R genes. One homolog of disease resistance protein (RPM1) was involved in the plant-pathogen interaction pathway. Artificial inoculation of disease/insect displayed resistance phenotypes against rice blast and brown planthopper in two lines. The results will provide allele-specific markers for rice molecular breeding.     

鸡6个功能基因microRNA靶标区域SNP的生物信息学预测

GDF-8、IGF-I、IGF-III、IGF2R、IGFBP2和GHR是鸡的重要经济性状候选基因。利用miRanda和Targetscan软件预测6个基因3′UTR潜在的microRNA靶标, 并发掘靶标区域SNP位点。结果表明: 在6个基因的26个microRNA靶标区域, 共检测到125个SNP位点, 在靶标及其5′和3′邻接等长侧翼区分别检测到47个、44个和35个SNPs位点, 其中12个SNP定位于靶标种子序列互补区。种子序列互补区及其3′侧翼区的SNP位点可能会影响microRNA的调控, 导致家禽的表型变异     

DNA sequence polymorphism of the Rhg4 candidate gene conferring resistance to soybean cyst nematode in Chinese domesticated and wild soybeans

Rhg4 is one of the major resistant genes conferring resistance to soybean cyst nematode races 1, 3 and 4. In order to better understand its sequence diversity among different Chinese soybean populations and the impact of human activities on it, we designed 5 primer sets based on its sequence deposited in Genbank (Genbank accession No. AF506518) to obtain the Rhg4 sequence from 104 Chinese cultivated and wild soybean genotypes, and then analyzed the DNA sequence polymorphism in different Chinese soybean populations. The alignment of Rhg4 sequence included 5,216 nucleotide base pairs. A total of 67 single nucleotide polymorphisms (SNPs) including 59 single base changes and 8 DNA insertion-deletions (InDels) were identified with a SNP frequency of 1/78. Except for a 14-base InDel, there were 29 SNPs in coding regions, and among them, 13 were non-synonymous (9 in functional domains with 1 in a leucine-rich repeats region, 2 in a transmembrane region and 6 in a Ser/Thr kinase domain). The probability of substitution at each site was not the same, there were two hot spots, one was in the 5'-untranslated region between positions 124 and 804, and the other was in the region between positions 2520 and 3733. Sequence diversity analysis among 104 soybean genotypes showed π?=?0.00102 and θ?=?0.00218 for Rhg4. A domestication bottleneck was found because of lower sequence diversity and 58% unique SNPs loss in landraces compared with Glycine soja. Intensive selection increased the sequence diversity of cultivars, which had higher diversity and more unique SNPs than landraces. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-012-9703-1) contains supplementary material, which is available to authorized users.     

Construction of EMS-Induced Peanut Mutant Libraries and Identification of Pod-Related Traits Mutant Lines

Peanut (Arachis hypogaea L.) is an oil and economic crop of vital importance, and peanut pod is the key organ influencing the yield and processing quality. Hence, the Pod-related traits (PRTs) are considered as important agronomic traits in peanut breeding. To broaden the variability of PRTs in current peanut germplasms, three elite peanut cultivars were used to construct Ethyl methane sulfonate (EMS)-induced mutant libraries in this study. The optimal EMS treatment conditions for the three peanut varieties were determined. It was found that the median lethal dose (LD50) of EMS treatment varied greatly among different genotypes. Finally, the EMS-induced peanut mutant libraries were constructed and a total of 124 mutant lines for PRTs were identified and evaluated. Furthermore, “M-8070”, one of the mutant lines for pod constriction, was re-sequenced via high-throughput sequencing technology. The genome-wide variations between “M-8070” and its wild parent “Fuhua 8” (FH 8) were detected. 2994 EMS-induced single nucleotide polymorphisms (SNPs) and 1188 insertion-deletions (InDels) between “M-8070” and its wild parent were identified. The predominant SNP mutation type was C/G to T/A transitions, while the predominant InDel mutation type was “1-bp”. We analyzed the distribution of identified mutations and annotated their functions. Most of the mutations (91.68% of the SNPs and 77.69% of the InDels) were located in the intergenic region. 72 SNPs were identified in the exonic region, leading to 27 synonymous, 43 non-synonymous and 2 stop-gain variation for gene structure. 13 Indels were identified in the exonic region, leading to 4 frame-shift, 8 non-frame-shift and 1 stop-gain variations of genes. These mutations may lead to the phenotypic variation of “M-8070”. Our study provided valuable resources for peanut improvement and functional genomic research.