A post-labeling method for multiplexed and multicolored genotyping analysis of SSR, indel and SNP markers in single tube with bar-coded split tag (BStag)

Allelic haplotype combinations at the MS-P1 region, including P-class pentatricopeptide repeat family genes, influence wide phenotypic variation in pollen grain number through a cytoplasmic male sterility model in citrus

In citrus breeding programs, male sterility is an important trait for developing seedless varieties. Sterility associated with the male sterile cytoplasm of Kishu mandarin (Kishu-cytoplasm) has been proposed to fit the cytoplasmic male sterility (CMS) model. However, it remains undetermined whether CMS in citrus is controlled by interactions between sterile cytoplasm and nuclear restorer-of-fertility (Rf) genes. Accordingly, mechanisms underlying the control of the wide phenotypic variation in pollen number for breeding germplasm should be elucidated. This study aimed to identify complete linkage DNA markers responsible for male sterility at the MS-P1 region based on fine mapping. Two P-class pentatricopeptide repeat (PPR) family genes were identified as candidates for Rf based on predicted mitochondrial localization and higher expression in a male fertile variety/selected strain than in a male sterile variety. Eleven haplotypes (HT1-HT11) at the MS-P1 region were defined based on genotyping of DNA markers. Association analysis of diplotypes at the MS-P1 region and the number of pollen grains per anther (NPG) in breeding germplasms harboring Kishu-cytoplasm revealed that the diplotypes in this region influenced NPG. Among these haplotypes, HT1 is a non-functional restorer-of-fertility (rf) haplotype; HT2, a less-functional Rf; HT3-HT5 are semi-functional Rfs; and HT6 and HT7 are functional Rfs. However, the rare haplotypes HT8-HT11 could not be characterized. Therefore, P-class PPR family genes in the MS-P1 region may constitute the nuclear Rf genes within the CMS model, and a combination of the seven haplotypes could contribute to phenotypic variation in the NPG of breeding germplasms. These findings reveal the genomic mechanisms of CMS in citrus and will contribute to seedless citrus breeding programs by selecting candidate seedless seedlings using the DNA markers at the MS-P1 region.

Parentage analysis of tea cultivars in Japan based on simple sequence repeat markers

Tea cultivars have been bred by individual selection of landraces and by crossbreeding, but the validation of the parentage is limited. In this study, we performed parentage analysis of 79 tea cultivars in Japan based on SSR markers to confirm or identify the parent-offspring relationships among them. The effectiveness of nine SSR markers for parentage analysis was validated by comparing them to the existing cleaved amplified polymorphic sequence markers. The former markers were detectable more alleles than the latter. Simulation of parentage analysis of the tea cultivars predicted biparental origins for 12 cultivars ('Houshun', 'Mie ryokuhou no. 1', 'Surugawase', 'Tenmyo', 'Yamanoibuki', 'Harumidori', 'Koushun', 'Minekaori', 'Okumusashi', 'Saemidori', 'Sofu', and 'Toyoka'), in the first five of which candidate parents of yet-to-be-defined pedigree were newly identified. Comparisons of a total of 41 SSR genotypes confirmed the newly-identified parentages of 'Asahi' for 'Tenmyo', 'Rokurou' for 'Houshun', 'Surugawase', and 'Yamanoibuki', and 'Yamatomidori' for 'Mie ryokuhou no. 1'. The maternity of seven cultivars out of the 12 was also confirmed with chloroplast DNA sequences. Uniparental origins were confirmed for 25 cultivars. This parentage analysis has improved our knowledge of tea pedigrees and will aid in the development of new cultivars.

Do aspermic (parthenogenetic) Fasciola forms have the ability to reproduce their progeny via parthenogenesis?

Across Far East Asia, aspermic Fasciola forms are found endemically. They have abnormal spermatogenesis and oogenesis, and it is presumed that their progeny are produced parthenogenetically and clonally. Because of this, they are also termed parthenogenic Fasciola forms. Currently, there is no evidence that they do indeed reproduce parthenogenetically and clonally. In this study, the multilocus genetic type (MLG) in 12 microsatellite markers of adult flukes and their subsequent progeny larvae were analysed using two laboratory aspermic Fasciola triploid strains. The MLGs of adults and their larvae were identical for all markers evaluated, suggesting that these strains reproduce their progeny clonally. The deviation between theoretical and actual frequency within the larvae genotype of the Fh_6 locus resulted in the inability for self-fertilization within individual adult flukes. These findings strongly suggested that aspermic Fasciola forms reproduce their progeny by means of parthenogenesis, possibly gynogenesis.

Identification of quantitative trait loci for powdery mildew resistance in highly resistant cucumber (Cucumis sativus L.) using ddRAD-seq analysis

Powdery mildew, caused by Podosphaera xanthii (syn. Sphaerotheca fuliginea ex Fr. Poll.), is one of the most economically important foliar diseases in cucumber (Cucumis sativus L.). Cucumber parental line 'Kyuri Chukanbohon Nou 5 Go', developed from weedy cucumber line CS-PMR1, is highly resistant to powdery mildew and is promising breeding material. We performed quantitative trait locus (QTL) analysis using double-digest restriction-site-associated DNA sequencing (ddRAD-Seq) in a population from a cross between 'Kyuri Chukanbohon Nou 5 Go' and the Japanese native cultivar 'Kaga-aonaga-fushinari', which is susceptible to powdery mildew. The resistance of the population and its parents was evaluated using leaf disc assays and image analysis. We detected one major QTL on Chr. 5 that was effective at both 20°C and 25°C and one minor QTL on Chr. 1 effective at 20°C. We detected two additional QTLs in subpopulation: one on Chr. 3 effective at 20°C and one on Chr. 5 effective at both 20°C and 25°C in a position different from the major QTL. The resistance alleles at all four QTLs were contributed by 'Kyuri Chukanbohon Nou 5 Go'. The results of this study can be used to develop practical DNA markers tightly linked to genes for powdery mildew resistance.

Genetic diversity and population structure analyses based on microsatellite DNA of parthenogenetic Fasciola flukes obtained from cattle and sika deer in Japan

Understanding the population structure of Fasciola flukes in domestic and wild animals is important for determining the extent of cross-infection between them. Although the parthenogenetic Fasciola flukes in Japan have been shown to comprise five genetic types based on the ribosomal internal transcribed spacer 1 (ITS1) and mitochondrial NADH dehydrogenase subunit 1 (nad1) regions, these genetic regions are not suitable for analyzing their population structure. In the present study, the genetic diversity and population structure of the parthenogenetic Fasciola flukes in Japan were studied using microsatellite DNA, ITS1, and nad1 regions. A total of 144 parthenogenetic Fasciola flukes, obtained from cattle and sika deer in 16 localities, were individually analyzed using PCR-RFLP for ITS1, PCR-direct sequence analysis for nad1, and post-labeling PCR and capillary electrophoresis for microsatellite DNA regions. The flukes showed higher genetic diversity in the microsatellite DNA regions than ITS1 and nad1. The population structures of parthenogenetic Fasciola flukes were unclear, however, it was suggested that the flukes are more diverse populations. We hypothesized that their distribution throughout Japan is closely related to livestock movement dependent on human activity. Moreover, it is considered that cross-infection of the flukes between cattle and sika deer possibly has occurred in the past.

Genetic diversity and population structure of feral rapeseed (Brassica napus L.) in Japan

Rapeseed (Brassica napus L.) is one of the most economically important oilseed crops worldwide. In Japan, it has been cultivated for more than a century and has formed many feral populations. The aim of this study was to elucidate the genetic diversity of feral rapeseeds by genotyping 537 individuals (among which 130 were determined to be genetically modified) sampled from various regions in Japan. Analysis of 30 microsatellite markers amplified 334 alleles and indicated moderate genetic diversity and high inbreeding (expected heterozygosity, 0.50; observed heterozygosity, 0.16; inbreeding coefficient within individuals, 0.68) within the feral populations. The Mantel test showed only an insignificant weak positive correlation between geographic distance and genetic distance. Analysis of molecular variance showed a greater genetic diversity among individuals than between populations. These results are in accordance with population structure assessed by using principal coordinate analysis and the program STRUCTURE, which showed that the 537 individuals could be assigned to 8 genetic clusters with very large genetic differences among individuals within the same geographic population, and that among feral individuals, many are closely related to rapeseed accessions in the NARO Genebank but some have unknown origins. These unique feral rapeseeds are likely to be affected by strong selection pressure. The results for genetically modified individuals also suggest that they have two different sources and have a considerable degree of diversity, which might be explained by hybridization with nearby individuals and separation of hybrid cultivars. The information obtained in this study could help improve the management of feral rapeseed plants in Japan.

Genetic diversity of leafy amaranth (Amaranthus tricolor L.) resources in Vietnam

Leafy amaranths, which are consumed as traditional food in Asia and Africa, are now considered among the most promising vegetables. In Vietnam, leafy amaranths, particularly Amaranthus tricolor L., are important summer vegetables due to their excellent nutritional values and high tolerance to biotic and abiotic stresses. However, this species has not been subjected to systematic breeding. Here we describe species identification and evaluation of the genetic diversity of Vietnamese amaranth collection by using matK and simple sequence repeats (SSR) markers. Our phylogenetic analysis based on the matK marker classified the species of 68% of the accessions, of which 120 belonged to A. tricolor. We developed 21 SSR markers, which amplified a total of 153 alleles in 294 A. tricolor accessions originating from Vietnam and overseas, with a mean allelic richness of 7.29 per marker, observed heterozygosity of 0.14, expected heterozygosity of 0.38, and polymorphic information content of 0.35. The STRUCTURE and F_ST analysis indicated a positive relationship between geographic distance and genetic differentiation among most of the overseas groups and the Vietnamese collection, but not among geographic groups within the Vietnamese collection. Vietnamese amaranths could be divided into two major types, one common in East Asia and the other one unique to Vietnam.

Classification of "nabana" (Brassica rapa) cultivars and landraces based on simple sequence repeat markers

Brassica rapa or B. napus vegetables for eating as young inflorescences and stalks are called "nabana". Japanese nabana includes "flower-bud type" and "stem-and-leaf type". Chinese and European types are also known (cai-xin, zicaitai, and broccoletto). We classified nabana belonging to B. rapa and other B. rapa vegetables. In a simple sequence repeat-based phylogram, 49 ingroup samples were classified into four groups (I-IV). Flower-bud and stem-and-leaf types were separated into groups I and III, respectively, with a slight overlap in group II. Cai-xin and non-heading Chinese cabbages were included in group IV. Broccoletto was placed in group III, close to turnips. Zicaitai cultivars were included in group II. We tested for clubroot resistance (CR) and its marker genotypes in nabana because of their agronomical importance. Ten cultivars were resistant to group 4 pathogen but not to group 2. Most of the CR cultivars had heterozygous resistance alleles in the CRb and Crr1 loci, consistent with inoculation tests. Our results suggest that Japanese nabana lines and foreign types were differentiated according to their consumption parts and cultivar origins, respectively. This study elucidates the relationships and CR properties of nabana and provides valuable information for the breeding of nabana cultivars.

QTL mapping of male sterility and transmission pattern in progeny of Satsuma mandarin

Seedlessness is one of the important traits in citrus breeding. Male sterility derived from Satsuma mandarin (Citrus unshiu) has been used in Japanese citrus breeding programs to obtain seedless cultivars. The efficiency of seedless cultivar breeding would be improved by developing a selection marker linked to seedlessness. In this study, we performed QTL mapping in ‘Okitsu No. 46’ × ‘Okitsu No. 56’ (O46-O56) crosses for the number of pollen grains per anther (NPG) and apparent pollen fertility (APF), two traits used as an index of male sterility, and detected a candidate QTL for NPG (MS-P1) on linkage group 8 with a significant LOD score (7.31) and 47% of variance explained. The QTL for APF (MS-F1) was detected on linkage group 6 with a significant LOD score (5.71) and 63.6% of variance explained. The role of both MS-P1 in reducing NPG and MS-F1 in decreasing APF were confirmed with the ‘Okitsu No.46’ × ‘Kara’ (O46-K) cross. Pedigree analysis inferred that both MS-P1 and MS-F1 in ‘Okitsu No. 46’ were derived from kunenbo (Citrus nobilis) through hassaku (C. hassaku) and ‘Sweet Spring’. Cytoplasm analysis revealed that both male-sterile ‘Sweet Spring’ and ‘Okitsu No. 46’ have cytoplasm derived from Kishu (C. kinokuni hort. ex Tanaka), but the cytoplasm of male-sterile kunenbo and hassaku were derived from other varieties rather than Kishu. These results suggest that MS-P1 and MS-F1 primarily reduce the NPG and decrease APF, but their expression requires a cytoplasm derived from Kishu. These findings will improve our understanding of the molecular mechanism of male sterility in citrus and help to develop a DNA marker for seedless breeding in citrus.

Analysis of genetic diversity of rapeseed genetic resources in Japan and core collection construction

Diversity analysis of rapeseed accessions preserved in the Japanese Genebank can provide valuable information for breeding programs. In this study, 582 accessions were genotyped with 30 SSR markers covering all 19 rapeseed chromosomes. These markers amplified 311 alleles (10.37 alleles per marker; range, 3-39). The genetic diversity of Japanese accessions was lower than that of overseas accessions. Analysis of molecular variance indicated significant genetic differentiation between Japanese and overseas accessions. Small but significant differences were found among geographical groups in Japan, and genetic differentiation tended to increase with geographical distance. STRUCTURE analysis indicated the presence of two main genetic clusters in the NARO rapeseed collection. With the membership probabilities threshold, 227 accessions mostly originating from overseas were assigned to one subgroup, and 276 accessions mostly originating from Japan were assigned to the other subgroup. The remaining 79 accessions are assigned to admixed group. The core collection constructed comprises 96 accessions of diverse origin. It represents the whole collection well and thus it may be useful for rapeseed genetic research and breeding programs. The core collection improves the efficiency of management, evaluation, and utilization of genetic resources.

Microsatellite markers for Nuphar japonica (Nymphaeaceae), an aquatic plant in the agricultural ecosystem of Japan

PREMISE OF THE STUDY

Nuphar species (Nymphaeaceae) are representative aquatic plants in irrigation ponds in Japanese agricultural ecosystems. We developed 15 polymorphic microsatellite markers for N. japonica and confirmed their utility for its close relatives N. oguraensis var. akiensis and N. ×saijoensis, which originated from natural hybridization between N. japonica and N. oguraensis.

METHODS AND RESULTS

Genetic variation was characterized in 15 polymorphic loci in three populations of N. japonica. The average number of alleles per locus was 3.47 (range = 2-9; n = 32), and the average expected heterozygosity per locus was 0.84 (range = 0.5-1.0); 11 loci were amplified in N. oguraensis var. akiensis and 15 in N. ×saijoensis.

CONCLUSIONS

The polymorphic microsatellite markers developed in this study will be useful for investigating the levels of genetic diversity within remnant populations of Nuphar taxa and could provide a valuable tool for conservation genetics of these taxa.

Hybrid Origins of Citrus Varieties Inferred from DNA Marker Analysis of Nuclear and Organelle Genomes

Most indigenous citrus varieties are assumed to be natural hybrids, but their parentage has so far been determined in only a few cases because of their wide genetic diversity and the low transferability of DNA markers. Here we infer the parentage of indigenous citrus varieties using simple sequence repeat and indel markers developed from various citrus genome sequence resources. Parentage tests with 122 known hybrids using the selected DNA markers certify their transferability among those hybrids. Identity tests confirm that most variant strains are selected mutants, but we find four types of kunenbo (Citrus nobilis) and three types of tachibana (Citrus tachibana) for which we suggest different origins. Structure analysis with DNA markers that are in Hardy-Weinberg equilibrium deduce three basic taxa coinciding with the current understanding of citrus ancestors. Genotyping analysis of 101 indigenous citrus varieties with 123 selected DNA markers infers the parentages of 22 indigenous citrus varieties including Satsuma, Temple, and iyo, and single parents of 45 indigenous citrus varieties, including kunenbo, C. ichangensis, and Ichang lemon by allele-sharing and parentage tests. Genotyping analysis of chloroplast and mitochondrial genomes using 11 DNA markers classifies their cytoplasmic genotypes into 18 categories and deduces the combination of seed and pollen parents. Likelihood ratio analysis verifies the inferred parentages with significant scores. The reconstructed genealogy identifies 12 types of varieties consisting of Kishu, kunenbo, yuzu, koji, sour orange, dancy, kobeni mikan, sweet orange, tachibana, Cleopatra, willowleaf mandarin, and pummelo, which have played pivotal roles in the occurrence of these indigenous varieties. The inferred parentage of the indigenous varieties confirms their hybrid origins, as found by recent studies.

Isolation and characterization of microsatellite loci for the endangered wetland plant Adenophora palustris (Campanulaceae)

PREMISE OF THE STUDY

Adenophora palustris (Campanulaceae) is an endangered wetland plant species in Japan. Although it is widely distributed in East Asia, only six extant populations are known in Japan, with fewer than 1000 individuals in total. We developed 15 microsatellite markers for this species and confirmed their utility for the closely related species A. triphylla var. japonica.

METHODS AND RESULTS

Ten polymorphic loci were characterized for genetic variation within three populations of A. palustris. The number of alleles per locus ranged from four to 15, with an average of 9.3; the expected heterozygosity ranged from 0.48 to 0.89, with an average of 0.74. Nine loci were successfully amplified in A. triphylla var. japonica, and three of these loci showed polymorphism.

CONCLUSIONS

These markers are useful for investigating genetic diversity and gene flow within and among remnant populations of A. palustris in Japan, and the results will provide crucial information for conservation.

Multicolor-based discrimination of 21 short tandem repeats and amelogenin using four fluorescent universal primers

The aim of this study was to develop a cost-effective genotyping method using high-quality DNA for human identification. A total of 21 short tandem repeats (STRs) and amelogenin were selected, and fluorescent fragments at 22 loci were simultaneously amplified in a single-tube reaction using locus-specific primers with 24-base universal tails and four fluorescent universal primers. Several nucleotide substitutions in universal tails and fluorescent universal primers enabled the detection of specific fluorescent fragments from the 22 loci. Multiplex polymerase chain reaction (PCR) produced intense FAM-, VIC-, NED-, and PET-labeled fragments ranging from 90 to 400 bp, and these fragments were discriminated using standard capillary electrophoretic analysis. The selected 22 loci were also analyzed using two commercial kits (the AmpFLSTR Identifiler Kit and the PowerPlex ESX 17 System), and results for two loci (D19S433 and D16S539) were discordant between these kits due to mutations at the primer binding sites. All genotypes from the 100 samples were determined using 2.5 ng of DNA by our method, and the expected alleles were completely recovered. Multiplex 22-locus genotyping using four fluorescent universal primers effectively reduces the costs to less than 20% of genotyping using commercial kits, and our method would be useful to detect silent alleles from commercial kit analysis.

Microsatellite marker development by multiplex ion torrent PGM sequencing: a case study of the endangered Odorrana narina complex of frogs

The endangered Ryukyu tip-nosed frog Odorrana narina and its related species, Odorrana amamiensis, Odorrana supranarina, and Odorrana utsunomiyaorum, belong to the family Ranidae and are endemically distributed in Okinawa (O. narina), Amami and Tokunoshima (O. amamiensis), and Ishigaki and Iriomote (O. supranarina and O. utsunomiyaorum) Islands. Because of varying distribution patterns, this species complex is an intrinsic model for speciation and adaptation. For effective conservation and molecular ecological studies, further genetic information is needed. For rapid, cost-effective development of several microsatellite markers for these and 2 other species, we used next-generation sequencing technology of Ion Torrent PGM™. Distribution patterns of repeat motifs of microsatellite loci in these modern frog species (Neobatrachia) were similarly skewed. We isolated and characterized 20 new microsatellite loci of O. narina and validated cross-amplification in the three-related species. Seventeen, 16, and 13 loci were cross-amplified in O. amamiensis, O. supranarina, and O. utsunomiyaorum, respectively, reflecting close genetic relationships between them. Mean number of alleles and expected heterozygosity of newly isolated loci varied depending on the size of each inhabited island. Our findings suggested the suitability of Ion Torrent PGM™ for microsatellite marker development. The new markers developed for the O. narina complex will be applicable in conservation genetics and molecular ecological studies.

An efficient technique for primer development and application that integrates fluorescent labeling and multiplex PCR

PREMISE OF THE STUDY

Development of genetic markers can be costly and time-consuming, especially when multiple primer pairs are fluorescently labeled. This step was streamlined by combining two techniques in the same PCR reaction: (1) custom-labeling of primers by the investigator and (2) multiplexing multiple primers together in the same reaction. •

METHODS AND RESULTS

This technique was successfully used to develop microsatellite markers in several plant species. Microsatellites amplified with this multiplexing process were identical to those generated from PCR using individual primer pairs and with traditional methods using a priori labeled fluorescent primers. Tests of PCR cycling programs revealed that conditions recommended for the commercial kit generated stronger fragment peaks than the previously recommended cycling protocol. •

CONCLUSIONS

This technique is an efficient and economical way to fluorescently label multiple microsatellite primers in the same reaction. It is also applicable to other markers used in PCR amplification of genetic material.

Genotyping of 38 insertion/deletion polymorphisms for human identification using universal fluorescent PCR

Short insertion/deletion (Indel) polymorphisms of approximately 2-6 bp are useful as biallelic markers for forensic analysis, and the application of Indel genotyping as a supplementary tool would improve human identification accuracy. We examined the allele frequencies of 37 autosomal Indels in the Japanese population and developed a novel dual-color genotyping method for human identification on the basis of universal fluorescent PCR, including the sex-typing amelogenin locus. Target genomic fragment sizes for 38 Indels were 49-143 bp. We analyzed these Indels in 100 Japanese individuals using the M13(-47) sequence as a universal primer. For dual-color genotyping, we designed a novel universal primer with high amplification efficiency and specificity. Using FAM-labeled M13(-47) and HEX-labeled modified M13(-47) primers, fluorescent signals at all loci were clearly distinguished in two independent multiplex PCRs. Average minor allele frequency was 0.39, and accumulated matching probability was 2.12 × 10(-15). Complete profiles were successfully amplified with as little as 0.25 ng of DNA. This method provides robust, sensitive, and cost-effective genotyping for human identification.

Universal fluorescent labeling of amplification products using locked nucleic acids

Amplification/hybridization-based genetic analyses using primers containing locked nucleic acids (LNAs) present many benefits. Here, we developed a novel design for universal fluorescent PCR using LNAs. Universal fluorescent PCR generates intermediate nonlabeled fragments and final fluorescent fragments in a two-step amplification process that uses locus-specific primers with universal tails and universal fluorescent primers. In this study, a few standard nucleotides were replaced with LNAs only in the fluorescent universal primers. The sequence of the fluorescent universal primer significantly affected the amplification efficiency. For primers with three LNAs, the fluorescent primers with stable M13(-47) sequences provided the most efficient signal (approximately tenfold higher than the primers with M13(-21) sequences at lower Tm values). Moreover, AT-rich LNA substitutions in the fluorescent primers produced much lower amplification efficiencies than GC-rich substitutions. GC-rich LNAs produced greater differences in Tm values among primers, and resulted in the preferential production of fluorescently labeled amplicons. The specificity and sensitivity of LNA-containing fluorescent primers were assessed by genotyping eight STRs in Japanese individuals, and full STR profiles could be generated using as little as 0.25 ng of genomic DNA. The method permitted clear discrimination of alleles and represents sensitive STR genotyping at a reduced cost.