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.     

Development and Evaluation of SoySNP50K,a High-Density Genotyping Array for Soybean

The objective of this research was to identify single nucleotide polymorphisms (SNPs) and to develop an Illumina Infinium BeadChip that contained over 50,000 SNPs from soybean (Glycine max L. Merr.). A total of 498,921,777 reads 35–45bp in length were obtained from DNA sequence analysis of reduced representation libraries from several soybean accessions which included six cultivated and two wild soybean (G. soja Sieb. et Zucc.) genotypes. These reads were mapped to the soybean whole genome sequence and 209,903 SNPs were identified. After applying several filters, a total of 146,161 of the 209,903 SNPs were determined to be ideal candidates for Illumina Infinium II BeadChip design. To equalize the distance between selected SNPs, increase assay success rate, and minimize the number of SNPs with low minor allele frequency, an iteration algorithm based on a selection index was developed and used to select 60,800 SNPs for Infinium BeadChip design. Of the 60,800 SNPs, 50,701 were targeted to euchromatic regions and 10,000 to heterochromatic regions of the 20 soybean chromosomes. In addition, 99 SNPs were targeted to unanchored sequence scaffolds. Of the 60,800 SNPs, a total of 52,041 passed Illumina’s manufacturing phase to produce the SoySNP50K iSelect BeadChip. Validation of the SoySNP50K chip with 96 landrace genotypes, 96 elite cultivars and 96 wild soybean accessions showed that 47,337 SNPs were polymorphic and generated successful SNP allele calls. In addition, 40,841 of the 47,337 SNPs (86%) had minor allele frequencies ≥10% among the landraces, elite cultivars and the wild soybean accessions. A total of 620 and 42 candidate regions which may be associated with domestication and recent selection were identified, respectively. The SoySNP50K iSelect SNP beadchip will be a powerful tool for characterizing soybean genetic diversity and linkage disequilibrium, and for constructing high resolution linkage maps to improve the soybean whole genome sequence assembly.     

Search for and analysis of single nucleotide polymorphisms (SNPs) in rice (Oryza sativa, Oryza rufipogon) and establishment of SNP markers.

We searched for SNPs in 417 regions distributed throughout the genome of three Oryza sativa ssp. japonica cultivars, two indica cultivars, and a wild rice (O. rufipogon). We found 2800 SNPs in approximately 250,000 aligned bases for an average of one SNP every 89 bp, or one SNP every 232 bp between two randomly selected strains. Graphic representation of the frequency of SNPs along each chromosome showed uneven distribution of polymorphism-rich and -poor regions, but little obvious association with the centromere or telomere. The 94 SNPs that we found between the closely related cultivars 'Nipponbare' and 'Koshihikari' can be converted into molecular markers. Our establishment of 213 co-dominant SNP markers distributed throughout the genome illustrates the immense potential of SNPs as molecular markers not only for genome research, but also for molecular breeding of rice.     

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.     

SNP Discovery and Linkage Map Construction in Cultivated Tomato

Few intraspecific genetic linkage maps have been reported for cultivated tomato, mainly because genetic diversity within Solanum lycopersicum is much less than that between tomato species. Single nucleotide polymorphisms (SNPs), the most abundant source of genomic variation, are the most promising source of polymorphisms for the construction of linkage maps for closely related intraspecific lines. In this study, we developed SNP markers based on expressed sequence tags for the construction of intraspecific linkage maps in tomato. Out of the 5607 SNP positions detected through in silico analysis, 1536 were selected for high-throughput genotyping of two mapping populations derived from crosses between ‘Micro-Tom’ and either ‘Ailsa Craig’ or ‘M82’. A total of 1137 markers, including 793 out of the 1338 successfully genotyped SNPs, along with 344 simple sequence repeat and intronic polymorphism markers, were mapped onto two linkage maps, which covered 1467.8 and 1422.7 cM, respectively. The SNP markers developed were then screened against cultivated tomato lines in order to estimate the transferability of these SNPs to other breeding materials. The molecular markers and linkage maps represent a milestone in the genomics and genetics, and are the first step toward molecular breeding of cultivated tomato. Information on the DNA markers, linkage maps, and SNP genotypes for these tomato lines is available at http://www.kazusa.or.jp/tomato/.     

水稻单核苷酸多态性及其应用现状

单核苷酸多态性（single nucleotide polymorphisms, SNPs）在水稻中数量多,分布密度高,遗传稳定性高。水稻SNPs的发现方法主要有对样本DNA的PCR产物直接测序、从SSR区段检测SNPs和从基因组序列直接搜索等。目前已有多种基因分型技术运用到了水稻SNPs检测,SNPs检测的高度自动化使水稻SNPs基因分型非常方便。单核苷酸多态性在水稻遗传图谱的构建、基因克隆和功能基因组学研究、标记辅助选择育种、遗传资源分类及物种进化等方面的应用具有巨大潜力。     

Population genetic structure of Penaeus monodon, in relation to monsoon current patterns in Southwest, East and Andaman coastal waters of India

The black tiger shrimp (Penaeus monodon), a commercially important penaeid species, is widely distributed across the Indo-Pacific region. Genetic diversity in P. monodon collected from eight geographical regions in Southwest, East and Andaman coastal waters of India (N = 418) was investigated using 10 polymorphic microsatellite loci. Average observed heterozygosity at sampled loci were high, ranging from 0.643 (Coromandel Coast) to 0.753 (South Andaman). Pairwise F_ST (ranged from 0.005 to 0.078) and R_ST (ranged from 0.005 to 0.171) estimates revealed surprisingly strong and statistically significant genetic structure among tiger shrimp populations. A synthetic map generated by multidimensional scaling shows an apparent cline in allele frequencies paralleling the roughly circular flow of surface currents in the Bay of Bengal. Significant heterozygote deficiencies were noted in most population samples at most loci. Andaman Island sites showed the highest diversity. Recognition of high genetic diversity and distinct population structuring of P. monodon in Indian seas has important implications for future domestication of this species in India, for two reasons: identification of the best wild founding stocks for aquaculture and, subsequently, the potential impacts of release of domesticates to the wild, either accidentally or deliberately (i.e. for stock enhancement).     

High-throughput genotyping in citrus accessions using an SNP genotyping array

We developed a 384 multiplexed SNP array, named CitSGA-1, for the genotyping of Citrus cultivars, and evaluated the performance and reliability of the genotyping. SNPs were surveyed by direct sequence comparison of the sequence tagged site (STS) fragment amplified from genomic DNA of cultivars representing the genetic diversity of citrus breeding in Japan. Among 1497 SNPs candidates, 384 SNPs for a high-throughput genotyping array were selected based on physical parameters of Illumina’s bead array criteria. The assay using CitSGA-1 was applied to a hybrid population of 88 progeny and 103 citrus accessions for breeding in Japan, which resulted in 73,726 SNP calls. A total of 351 SNPs (91 %) could call different genotypes among the DNA samples, resulting in a success rate for the assay comparable to previously reported rates for other plant species. To confirm the reliability of SNP genotype calls, parentage analysis was applied, and it indicated that the number of reliable SNPs and corresponding STSs were 276 and 213, respectively. The multiplexed SNP genotyping array reported here will be useful for the efficient construction of linkage map, for the detection of markers for marker-assisted breeding, and for the identification of cultivars.     

Limited RNA Editing in Exons of Mouse Liver and Adipose

Several studies have investigated RNA–DNA differences (RDD), presumably due to RNA editing, with conflicting results. We report a rigorous analysis of RDD in exonic regions in mice, taking into account critical biases in RNA-Seq analysis. Using deep-sequenced F1 reciprocal inbred mice, we mapped 40 million RNA-Seq reads per liver sample and 180 million reads per adipose sample. We found 7300 apparent hepatic RDDs using a multiple-site mapping procedure, compared with 293 RDD found using a unique-site mapping procedure. After filtering for repeat sequence, splice junction proximity, undirectional strand, and extremity read bias, 63 RDD remained. In adipose tissue unique-site mapping identified 1667 RDD, and after applying the same four filters, 188 RDDs remained. In both tissues, the filtering procedure increased the proportion of canonical (A-to-I and C-to-U) editing events. The genomic DNA of 12 RDD sites among the potential 63 hepatic RDD was tested by Sanger sequencing, three of which proved to be due to unreferenced SNPs. We validated seven liver RDD with Sequenom technology, including two noncanonical, Gm5424 C-to-I(G) and Pisd I(G)-to-A RDD. Differences in diet, sex, or genetic background had very modest effects on RDD occurrence. Only a small number of apparent RDD sites overlapped between liver and adipose, indicating a high degree of tissue specificity. Our findings underscore the importance of properly filtering for bias in RNA-Seq investigations, including the necessity of confirming the DNA sequence to eliminate unreferenced SNPs. Based on our results, we conclude that RNA editing is likely limited to hundreds of events in exonic RNA in liver and adipose.     

De novo assembly and characterization of the complete chloroplast genome of radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is an edible root vegetable crop that is cultivated worldwide and whose genome has been sequenced. Here we report the complete nucleotide sequence of the radish cultivar WK10039 chloroplast (cp) genome, along with a de novo assembly strategy using whole genome shotgun sequence reads obtained by next generation sequencing. The radish cp genome is 153,368 bp in length and has a typical quadripartite structure, composed of a pair of inverted repeat regions (26,217 bp each), a large single copy region (83,170 bp), and a small single copy region (17,764 bp). The radish cp genome contains 87 predicted protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Sequence analysis revealed the presence of 91 simple sequence repeats (SSRs) in the radish cp genome.     

Single nucleotide polymorphisms in the mitochondrial control region are associated with metabolic phenotypes and oxidative stress

Objective

To identify the mitochondrial DNA (mtDNA) single nucleotide polymorphisms (SNPs) in the control region and elucidate their role in metabolic phenotypes and oxidative stress.

Methods

A total of 861 nondiabetic subjects were enrolled, including 250 impaired fasting glucose (IFG) and 370 obese subjects (body mass index [BMI] > 25 kg/m²). Antioxidant status presented as total free thiol level was determined from serum samples. DNA was extracted from peripheral blood leucocytes, and the sequences were analyzed using the DNASTAR software. SNPs were identified by comparison with the Cambridge Reference Sequence.

Results

After adjusting odds ratios for age, sex, and BMI, the selected independently significant SNPs indicated 4 susceptible SNPs: SNP-16126C and SNP-16261T, which were related to abdominal obesity (P = 0.009; 0.06); SNP-16390A, related to hypertension (HTN) (P = 0.007); and SNP-16092C, related to decreased antioxidant capacity (P = 0.015). In the obese subgroup, 3 susceptible SNPs included SNP-16189C and SNP-16260T, which showed significantly higher IFG prevalence (P = 0.016 and 0.024, respectively), and SNP-16519C, which was significantly higher in the HTN group (P = 0.036). As to protective SNPs, 5 protective SNPs were identified in all subjects but only one SNP-16093C is consistent in obese group, which showed a significantly lower prevalence in patients with abdominal obesity and was associated with a higher antioxidant status (P < 0.001).

Conclusion

SNPs in the mtDNA control region are associated with metabolic phenotypes and oxidative stress markers. Some SNPs are relating to the interaction between obesity and genetic factors. The beneficial effects of these protective SNPs were insignificant and some susceptible SNPs became dominant within the obese subgroup. Subjects harboring these SNPs should avoid excessive weight gain.     

Molecular characterization of two y-type high molecular weight glutenin subunit alleles 1Ay12 and 1Ay8 from cultivated einkorn wheat (Triticum monococcum ssp. monococcum)

Two y-type high molecular weight glutenin subunits (HMW-GSs) 1Ay12^? and 1Ay8^? from the two accessions PI560720 and PI345186 of cultivated einkorn wheat (Triticum monococcum ssp. monococcum, AA, 2n = 2x = 14), were identified by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). The mobility of 1Ay12^? and 1Ay8^? was similar to that of 1Dy12 and 1By8 from common wheat Chinese Spring, respectively. Their ORFs respectively consisted of 1812 bp and 1935 bp, encoding 602 and 643 amino acid residues with the four typical structural domains of HMW-GS including signal peptide, conserved N-, and C-terminal and central repetitive domains. Compared with the most similar active 1Ay alleles previous published, there were a total of 15 SNPs and 2 InDels in them. Their encoding functions were confirmed by successful heterogeneous expression. The two novel 1Ay alleles were named as 1Ay12^? and 1Ay8^? with the accession No. JQ318694 and JQ318695 in GenBank, respectively. The two alleles were classed into the two distinct groups, Phe-type and Cys-type, which might be relevant to the differentiation of Glu-A1-2 alleles. Of which, 1Ay8^? belonged to Cys-type group, and its protein possessed an additional conserved cysteine residue in central repetitive region besides the six common ones in N- and C-terminal regions of Phe-type group, and was the second longest in all the known active 1Ay alleles. These results suggested that the subunit 1Ay8^? of cultivated einkorn wheat accession PI345186 might have a potential ability to strengthen the gluten polymer interactions and be a valuable genetic resource for wheat quality improvement.     

A high-throughput apple SNP genotyping platform using the GoldenGate™ assay

EST data generated from 14 apple genotypes were downloaded from NCBI and mapped against a reference EST assembly to identify Single Nucleotide Polymorphisms (SNPs). Mapping of these SNPs was undertaken using 90% of sequence similarity and minimum coverage of four reads at each SNP position. In total, 37,807 SNPs were identified with an average of one SNP every 187 bp from a total of 6888 unique EST contigs. Identified SNPs were checked for flanking sequences of ≥ 60 bp along both sides of SNP alleles for reliable design of a custom high-throughput genotyping assay. A total of 12,299 SNPs, representing 6525 contigs, fit the selected criterion of ≥ 60 bp sequences flanking a SNP position. Of these, 1411 SNPs were validated using four apple genotypes. Based on genotyping assays, it was estimated that 60% of SNPs were valid SNPs, while 26% of SNPs might be derived from paralogous regions.     

Draft Sequences of the Radish (Raphanus sativus L.) Genome

Radish (Raphanus sativus L., n = 9) is one of the major vegetables in Asia. Since the genomes of Brassica and related species including radish underwent genome rearrangement, it is quite difficult to perform functional analysis based on the reported genomic sequence of Brassica rapa. Therefore, we performed genome sequencing of radish. Short reads of genomic sequences of 191.1 Gb were obtained by next-generation sequencing (NGS) for a radish inbred line, and 76,592 scaffolds of ≥300 bp were constructed along with the bacterial artificial chromosome-end sequences. Finally, the whole draft genomic sequence of 402 Mb spanning 75.9% of the estimated genomic size and containing 61,572 predicted genes was obtained. Subsequently, 221 single nucleotide polymorphism markers and 768 PCR-RFLP markers were used together with the 746 markers produced in our previous study for the construction of a linkage map. The map was combined further with another radish linkage map constructed mainly with expressed sequence tag-simple sequence repeat markers into a high-density integrated map of 1,166 cM with 2,553 DNA markers. A total of 1,345 scaffolds were assigned to the linkage map, spanning 116.0 Mb. Bulked PCR products amplified by 2,880 primer pairs were sequenced by NGS, and SNPs in eight inbred lines were identified.     

Development of single nucleotide polymorphism markers specific to Apis mellifera (Hymenoptera: Apidae) line displaying high hygienic behavior against Varroa destructor,an ectoparasitic mite

To control Varroa destructor, an ectoparasitic mite, a honey bee line possessing high hygienic behavior (HHB) against this mite has been bred in South Korea. However, a method that can diagnose and assess the HHB line from control (the low hygienic behavior, LHB) line has not been reported yet. Thus, the objective of this study was to develop single nucleotide polymorphism (SNP) markers through whole-genome sequencing of worker bees from HHB line of A. mellifera caucasica and LHB line of A. m. carnica (Hymenoptera: Apidae). A total of 319,445,977 sequence reads were mapped to the known A. mellifera reference genome (average coverage of 87.46%). In 2,316,128 and 3,266,756 SNPs from HHB and LHB line, respectively, 20 SNPs that showed homozygosity in each line were selected and eight SNPs were used to diagnose the HHB line either by typical PCR-restriction fragment length polymorphism or allele-specific PCR. Six of remaining SNPs were of different sizes, enabling relatively easy differentiation of these two honey bee lines on typical agarose gel. Another remaining six SNPs had different sequences, including SNP sites. These SNP markers can be used to diagnose and assess V. destructor-specific HHB line of honey bees.