Molecular cytogenetic analysis of the Vigna species distributed in Korea

Vigna plants distributed in Korea were analyzed by molecular cytogenetic fluorescence in situ hybridization (FISH), genomic in situ hybridization (GISH) and rDNA ITS/NTS sequences. FISH revealed that variable 45S rRNA gene loci (one to four) were localized on the terminal regions of chromosomes, while two conserved 5S rRNA gene loci from all species examined, except for rice bean (single locus), were detected. FISH and GISH showed the characteristic organization of rRNA gene loci and genomic homology on the chromosomes, indicating their cytogenetic relatationships. ITS sequence revealed that there was considerable variation in length (190–207 bp in ITS1, 205–221 bp in ITS2) and nucleotide composition (7–67 bp). The 5S rRNA gene unit comprised coding region (118 bp) and extensive sequence heterogeneity (97–221 bp). Phylogenetic analysis of the ITS and NTS sequences demonstrated that the Vigna species are divided into two groups: angularis (V. angularis, V. umbellata, V. nakashimae and V. nipponensis) and unguiculata (V. unguiculata, V. sesquipedalis and V. vexillata). Sequence data also showed that mung bean was closer to the angularis group.     

Genome-wide SNP discovery in mungbean by Illumina HiSeq

Mungbean [Vigna radiata (L.) Wilczek], a self-pollinated diploid plant with 2n = 22 chromosomes, is an important legume crop with a high-quality amino acid profile. Sequence variation at the whole-genome level was examined by comparing two mungbean cultivars, Sunhwanokdu and Gyeonggijaerae 5, using Illumina HiSeq sequencing data. More than 40 billion bp from both mungbean cultivars were sequenced to a depth of 72×. After de novo assembly of Sunhwanokdu contigs by ABySS 1.3.2 (N50 = 9,958 bp), those longer than 10 kb were aligned with Gyeonggijaerae 5 reads using the Burrows–Wheeler Aligner. SAMTools was used for retrieving single nucleotide polymorphisms (SNPs) between Sunhwanokdu and Gyeonggijaerae 5, defining the lowest and highest depths as 5 and 100, respectively, and the sequence quality as 100. Of the 305,504 single-base changes identified, 40,503 SNPs were considered heterozygous in Gyeonggijaerae 5. Among the remaining 265,001 SNPs, 65.9 % (174,579 cases) were transitions and 34.1 % (90,422 cases) were transversions. For SNP validation, a total of 42 SNPs were chosen among Sunhwanokdu contigs longer than 10 kb and sharing at least 80 % sequence identity with common bean expressed sequence tags as determined with est2genome. Using seven mungbean cultivars from various origins in addition to Sunhwanokdu and Gyeonggijaerae 5, most of the SNPs identified by bioinformatics tools were confirmed by Sanger sequencing. These genome-wide SNP markers could enrich the current molecular resources and might be of value for the construction of a mungbean genetic map and the investigation of genetic diversity.     

A genetic linkage map of azuki bean constructed with molecular and morphological markers using an interspecific population (Vigna angularis x V. nakashimae)

A genetic linkage map of azuki bean (Vigna angularis) was constructed with molecular and morphological markers using an F₂ population of an interspecific cross between azuki bean and its wild relative, V. nakashimae. In total, 132 markers (108 RAPD, 19 RFLP and five morphological markers) were mapped in 14 linkage groups covering 1250 cM; ten remained unlinked. The clusters of markers showing distorted segregation were found in linkage groups 2, 8 and 12. By comparing the azuki linkage map with those of mungbean and cowpea, using 20 RFLP common markers, some sets of the markers were found to belong to the same linkage groups of the respective maps, indicating that these linkage blocks are conserved among the three Vigna species. This map provides a tool for markerassisted selection and for studies of genome organization in Vigna species.     

Comparative molecular mapping in Ceratotropis species using an interspecific cross between azuki bean (Vigna angularis) and rice bean (V. umbellata)

A genetic linkage map was developed with 86 F₂ plants derived from an interspecific cross between azuki bean (Vigna angularis, 2n=2x=22) and rice bean (V. umbellata, 2n=2x=22). In total, 14 linkage groups, each containing more than 4 markers, were constructed with one phenotypic, 114 RFLP and 74 RAPD markers. The total map size was 1702 cM, and the average distance between markers was 9.7 cM. The loci showing significant deviation from the expected ratio clustered in several linkage groups. Most of the skewed loci were due to the predominance of rice bean alleles. The azuki-rice bean linkage map was compared with other available maps of Vigna species in subgenus Ceratotropis. Based on the lineage of the common mapped markers, 7 and 16 conserved linkage blocks were found in the interspecific map of azuki bean ×V. nakashimae and mungbean map, respectively. Although the present map is not fully saturated, it may facilitate gene tagging, QTL mapping and further useful gene transfer for azuki bean breeding. Received: 20 March 1999 / Accepted: 29 April 1999     

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 single nucleotide polymorphism markers in the large and complex rubber tree genome using next-generation sequence data

The development of single nucleotide polymorphism (SNP) markers provides the opportunity to improve many areas of plant breeding and population genetics. Unfortunately, for species such as the rubber tree (Hevea brasiliensis), the use of next-generation sequencing for genomic SNP discovery is very difficult because of the large genome size and the abundance of repeated sequences. Access to a set of validated SNP markers is a significant advantage for rubber researchers who wish to apply SNPs in scientific research. Here, we performed genomic sequencing of H. brasiliensis and generated 10,993,648 short reads, which were assembled into 10,071 contigs (N50 = 3078) by a de novo assembly strategy. A total of 2446 contigs presented no hits in the current H. brasiliensis genome assembly and may therefore be considered novel genomic sequences of rubber tree. A total of 143 putative polymorphic positions were selected, gene annotations were available for 58.7 % of the markers, and all of the sequences could be anchored to the released H. brasiliensis genome. These SNPs were validated in eight genotypes of H. brasiliensis and 15 F1 plants from a mapping population, resulting in 30 (20.9 %) positions correctly classified. The analysis revealed key candidate genes responsible for defence mechanisms and provided markers for further genetic improvement of Hevea in breeding programmes.     

Transcriptome sequencing to produce SNP-based genetic maps of onion

We used the Roche-454 platform to sequence from normalized cDNA libraries from each of two inbred lines of onion (OH1 and 5225). From approximately 1.6 million reads from each inbred, 27,065 and 33,254 cDNA contigs were assembled from OH1 and 5225, respectively. In total, 3,364 well supported single nucleotide polymorphisms (SNPs) on 1,716 cDNA contigs were identified between these two inbreds. One SNP on each of 1,256 contigs was randomly selected for genotyping. OH1 and 5225 were crossed and 182 gynogenic haploids extracted from hybrid plants were used for SNP mapping. A total of 597 SNPs segregated in the OH1 × 5225 haploid family and a genetic map of ten linkage groups (LOD ≥8) was constructed. Three hundred and thirty-nine of the newly identified SNPs were also mapped using a previously developed segregating family from BYG15-23 × AC43, and 223 common SNPs were used to join the two maps. Because these new SNPs are in expressed regions of the genome and commonly occur among onion germplasms, they will be useful for genetic mapping, gene tagging, marker-aided selection, quality control of seed lots, and fingerprinting of cultivars.     

Use of an expressed sequence tag-based method for single nucleotide polymorphism identification and discrimination of Citrus species and cultivars

Sweet orange [Citrus sinensis (L.) Osbeck] represents the most important Citrus species, followed by clementine (C. clementina Hort. ex Tan.). Citrus species and genotypes are difficult to recognize as they have a moderate level of diversity due to nucellar selection, vegetative propagation and origin by single spontaneous mutation. Despite the large number of available sequences and the existence of a draft assembly of sweet orange and clementine, there are currently no single nucleotide polymorphism (SNP) databases for Citrus species. For this purpose, the QualitySNP software was used to discover SNPs in 19 Citrus species starting from 540,000 expressed sequence tags (ESTs) assembled in 52,000 contigs. The vast majority of ESTs, contigs and SNPs were found in C. clementina and C. sinensis: 4,400 out of 16,000 contigs (27 %) of C. clementina and 4,100 out of 17,000 contigs (24 %) of C. sinensis contained putative SNPs. A total of 3,634 sequences were associated with enzymes belonging to 121 metabolic KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, among which the secondary metabolite pathway was the most represented. A total of 163 SNPs from 52 contigs and genes of specific functional categories were validated and 81 polymorphic sites were found. Thirty-seven selected SNPs, validated by Sanger sequencing, confirmed that polymorphisms were mainly between species, while poor within-species variability was discovered. This work provides a collection of 15,879 putative SNP markers that could be exploited by the Citrus community. Furthermore, the validated SNPs associated with specific genes could be used for functional genetic studies in germplasm diversity analysis, mapping and breeding.     

Development of genomic simple sequence repeat markers in faba bean by next-generation sequencing

Faba bean (Vicia faba L.) is an important food legume crop with a huge genome. Development of genetic markers for faba bean is important to study diversity and for molecular breeding. In this study, we used Next Generation Sequencing (NGS) technology for the development of genomic simple sequence repeat (SSR) markers. A total of 14,027,500 sequence reads were obtained comprising 4,208 Mb. From these reads, 56,063 contigs were assembled (16,367 Mb) and 2138 SSRs were identified. Mono and dinucleotides were the most abundant, accounting for 57.5 % and 20.9 % of all SSR repeats, respectively. A total of 430 primer pairs were designed from contigs larger than 350 nucleotides and 50 primers pairs were tested for validation of SSR locus amplification. Nearly all (96 %) of the markers were found to produce clear amplicons and to be reproducible. Thirty-nine SSR markers were then applied to 46 faba bean accessions from worldwide origins, resulting in 161 alleles with 87.5 % polymorphism, and an average of 4.1 alleles per marker. Gene diversity (GD) of the markers ranged from 0 to 0.48 with an average of 0.27. Testing of the markers showed that they were useful in determining genetic relationships and population structure in faba bean accessions.     

The genetic contribution of CIDEA polymorphisms, haplotypes and loci interaction to obesity in a Han Chinese population

To investigate the association of tag-SNPs and haplotype structures of the CIDEA gene with obesity in a Han Chinese population. Five single nucleotide polymorphisms (SNPs) (rs1154588/V115F, rs4796955/SNP1, rs8092502/SNP2, rs12962340/SNP3 and rs7230480/SNP4) in the CIDEA gene were genotyped in a case–control study. Genotyping was performed using the sequenom matrix-assisted laser desorption/ionization time-of-flight mass spectrometry iPLEX platform. There were significant differences between the obese and control groups in genotype distributions of V115F (P < 0.001), SNP1 (P = 0.006) and SNP2 (P = 0.005). Carriers of V115F-TT, SNP1-GG and SNP2-CC genotypes had a 2.84-fold (95 % CI 1.73–4.66), 2.19-fold (95 % CI 1.09–4.38) and 4.37-fold (95 % CI 1.21–15.08) increased risk for obesity, respectively. Haplotype analysis showed that GTTC (SNP1/SNP2/V115F/SNP4) had 1.41-fold (95 % CI 1.02–1.95) increased risk for obesity; whereas, haplotype TTGC had 0.48-fold (95 % CI 0.24–0.96) decreased risk for obesity. Using the multifactor dimensionality reduction method, the best model including SNP1, SNP2, V115F and SNP4 polymorphisms was identified with a maximum testing accuracy to 59.32 % and a perfect cross-validation consistency of 10/10 (P = 0.011). Logistic analysis indicated that there was a significant interaction between SNP1 and V115F associated with obesity. Subjects having both genotypes of SNP1/GG and V115F/TT were more susceptible to obesity in the Han Chinese population (OR 2.66, 95 %: 1.22–5.80). Genotypes of V115F/TT, SNP1/GG and SNP2/CC and haplotype GTTC of CIDEA gene were identified as risk factors for obesity in the Han Chinese population. The interaction between SNP1 and V115F could play a joint role in the development of obesity.     

Direct Chloroplast Sequencing: Comparison of Sequencing Platforms and Analysis Tools for Whole Chloroplast Barcoding

Direct sequencing of total plant DNA using next generation sequencing technologies generates a whole chloroplast genome sequence that has the potential to provide a barcode for use in plant and food identification. Advances in DNA sequencing platforms may make this an attractive approach for routine plant identification. The HiSeq (Illumina) and Ion Torrent (Life Technology) sequencing platforms were used to sequence total DNA from rice to identify polymorphisms in the whole chloroplast genome sequence of a wild rice plant relative to cultivated rice (cv. Nipponbare). Consensus chloroplast sequences were produced by mapping sequence reads to the reference rice chloroplast genome or by de novo assembly and mapping of the resulting contigs to the reference sequence. A total of 122 polymorphisms (SNPs and indels) between the wild and cultivated rice chloroplasts were predicted by these different sequencing and analysis methods. Of these, a total of 102 polymorphisms including 90 SNPs were predicted by both platforms. Indels were more variable with different sequencing methods, with almost all discrepancies found in homopolymers. The Ion Torrent platform gave no apparent false SNP but was less reliable for indels. The methods should be suitable for routine barcoding using appropriate combinations of sequencing platform and data analysis.     

Adzuki bean: A new resource of biocatalyst for asymmetric reduction of aromatic ketones with high stereoselectivity and substrate tolerance

A new resource of biocatalyst for asymmetric reduction of aromatic ketones has been discovered for the first time from a common plant seed, adzuki bean, i.e. Phaseolus angularis (Willd.) W.F. Wight. The study investigated the best methods to prepare the biocatalyst and its ability to reduce ketones. Our results indicated that the biocatalyst from adzuki bean could reduce various aromatic ketones at relatively high concentrations (e.g. 100 mM), exhibiting excellent stereoselectivity (>98% e.e.). In addition, it was found that NADPH acts as the reducing cofactor, which can be regenerated by the crude enzyme system itself using glucose as an auxiliary substrate.     

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.     

Rapid Development of Microsatellite Markers for Callosobruchus chinensis Using Illumina Paired-End Sequencing

Background

The adzuki bean weevil, Callosobruchus chinensis L., is one of the most destructive pests of stored legume seeds such as mungbean, cowpea, and adzuki bean, which usually cause considerable loss in the quantity and quality of stored seeds during transportation and storage. However, a lack of genetic information of this pest results in a series of genetic questions remain largely unknown, including population genetic structure, kinship, biotype abundance, and so on. Co-dominant microsatellite markers offer a great resolving power to determine these events. Here, we report rapid microsatellite isolation from C. chinensis via high-throughput sequencing.

Principal Findings

In this study, 94,560,852 quality-filtered and trimmed reads were obtained for the assembly of genome using Illumina paired-end sequencing technology. In total, the genome with total length of 497,124,785 bp, comprising 403,113 high quality contigs was generated with de novo assembly. More than 6800 SSR loci were detected and a suit of 6303 primer pair sequences were designed and 500 of them were randomly selected for validation. Of these, 196 pair of primers, i.e. 39.2%, produced reproducible amplicons that were polymorphic among 8 C. chinensis genotypes collected from different geographical regions. Twenty out of 196 polymorphic SSR markers were used to analyze the genetic diversity of 18 C. chinensis populations. The results showed the twenty SSR loci were highly polymorphic among these populations.

Conclusions

This study presents a first report of genome sequencing and de novo assembly for C. chinensis and demonstrates the feasibility of generating a large scale of sequence information and SSR loci isolation by Illumina paired-end sequencing. Our results provide a valuable resource for C. chinensis research. These novel markers are valuable for future genetic mapping, trait association, genetic structure and kinship among C. chinensis.