An EST-SSR linkage map of Raphanus sativus and comparative genomics of the Brassicaceae

HASCH - A high-throughput amplicon-based SNP-platform for medicinal cannabis and industrial hemp genotyping applications

BACKGROUND

Cannabis sativa is seeing a global resurgence as a food, fiber and medicinal crop for industrial hemp and medicinal Cannabis industries respectively. However, a widespread moratorium on the use and research of C. sativa throughout most of the 20th century has seen the development of improved cultivars for specific end uses lag behind that of conventional crops. While C. sativa research and development has seen significant investments in the recent past, resulting in a suite of publicly available genomic resources and tools, a versatile and cost-effective mid-density genotyping platform for applied purposes in breeding and pre-breeding is lacking. Here we report on a first mid-density fixed-target SNP platform for C. sativa.

RESULTS

The High-throughput Amplicon-based SNP-platform for medicinal Cannabis and industrial Hemp (HASCH) was designed using a combination of filtering and Integer Linear Programming on publicly available whole-genome sequencing and RNA sequencing data, supplemented with in-house generated genotyping-by-sequencing (GBS) data. HASCH contains 1,504 genome-wide targets of high call rate (97% mean) and even distribution across the genome, designed to be highly informative (> 0.3 minor allele frequency) across both medicinal cannabis and industrial hemp gene pools. Average numbers of mismatch SNP between any two accessions were 251 for medicinal cannabis (N = 116) and 272 for industrial hemp (N = 87). Comparing HASCH data with corresponding GBS data on a collection of diverse C. sativa accessions demonstrated high concordance and resulted in comparable phylogenies and genetic distance matrices. Using HASCH on a segregating F2 population derived from a cross between a tetrahydrocannabinol (THC)-dominant and a cannabidiol (CBD)-dominant accession resulted in a genetic map consisting of 310 markers, comprising 10 linkage groups and a total size of 582.7 cM. Quantitative Trait Locus (QTL) mapping identified a major QTL for CBD content on chromosome 7, consistent with previous findings.

CONCLUSION

HASCH constitutes a versatile, easy to use and cost-effective genotyping solution for the rapidly growing Cannabis research community. It provides consistent genetic fingerprints of 1504 SNPs with wide applicability genetic resource management, quantitative genetics and breeding.

Development and utility of SSR markers based on Brassica sp. whole-genome in triangle of U

Introduction

Simple sequence repeats (SSR), also known as microsatellites, are crucial molecular markers in both animals and plants. Despite extensive previous research on SSRs, the development of microsatellite markers in Brassica crops remains limited and inefficient.

Methods

Krait software was used to identify microsatellites by genome-wide and marker development based on three recently sequenced basic species of Brassica crops in the triangle of U (Brassica rapa, B. nigra and B. oleracea), as well as three allotetraploids (B. juncea, B. napus and B. carinata) using public databases. Subsequently, the primers and the characteristics of microsatellites for most of them were accordingly designed on each chromosome of each of the six Brassica species, and their physical locations were identified,and the cross-transferability of primers have been carried out. In addition, a B-genome specific SSR marker was screened out.

Results

A total of 79341, 92089, 125443, 173964, 173604, and 222160 SSR loci have been identified from the whole genome sequences of Brassica crops within the triangle of U crops, B. rapa (AA), B. nigra (BB), B. oleracea (CC), B. napus (AACC), B. juncea (AABB) and B. carinata (BBCC), respectively. Comparing the number distribution of the three allotetraploid SSR loci in the three subgenomes AA, BB and CC, results indicate that the allotetraploid species have significant reduction in the number of SSR loci in the genome compared with their basic diploid counterparts. Moreover, we compared the basic species with their corresponding varieties, and found that the microsatellite characters between the allotetraploids and their corresponding basic species were very similar or almost identical. Subsequently, each of the 40 SSR primers was employed to investigate the polymorphism potential of B. rapa (85.27%), B. nigra (81.33%) and B. oleracea (73.45%), and B. rapa was found to have a higher cross-transfer rate among the basic species in the triangle of U. Meanwhile, a B-genome specific SSR marker, BniSSR23228 possessing the (AAGGA)3 sequence characteristics was obtained, and it located in chromosome B3 with a total length of 97 bp.

Discussion

In this study, results suggest that the pattern of distribution may be highly conserved during the differentiation of basic Brassica species and their allotetraploid counterparts. Our data indicated that the allotetraploidization process resulted in a significant reduction in SSR loci in the three subgenomes AA, BB and CC. The reasons may be partial gene dominated chromosomal homologous recombination and rearrangement during the evolution of basic diploid species into allotetraploids. This study provides a basis for future genomics and genetic research on the relatedness of Brassica species.

Quantitative trait locus (QTL) analysis and fine-mapping for Fusarium oxysporum disease resistance in Raphanus sativus using GRAS-Di technology

Fusarium wilt is a significant disease in radish, but the genetic mechanisms controlling yellows resistance (YR) are not well understood. This study aimed to identify YR-QTLs and to fine-map one of them using F2:3 populations developed from resistant and susceptible radish parents. In this study, two high-density genetic maps each containing shared co-dominant markers and either female or male dominant markers that spanned 988.6 and 1127.5 cM with average marker densities of 1.40 and 1.53 cM, respectively, were generated using Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di) technology. We identified two YR-QTLs on chromosome R2 and R7, and designated the latter as ForRs1 as the major QTL. Fine mapping narrowed down the ForRs1 locus to a 195 kb region. Among the 16 predicted genes in the delimited region, 4 genes including two receptor-like protein and -kinase genes (RLP/RLK) were identified as prime candidates for ForRs1 based on the nucleotide sequence comparisons between the parents and their predicted functions. This study is the first to use a GRAS-Di for genetic map construction of cruciferous crops and fine map the YR-QTL on the R7 chromosome of radish. These findings will provide groundbreaking insights into radish YR breeding and understanding the genetics of YR mechanism.

Genetic and Transcriptome Analysis of Leaf Trichome Development in Chinese Cabbage (Brassica rapa L. subsp. pekinensis) and Molecular Marker Development

Chinese cabbage (Brassica rapa L. subsp. pekinensis) is one of the vegetables with the largest cultivated area in China and has been a great addition to the daily diet of Chinese people. A genetic map has been constructed in our previous study using the F₂ population of two inbred lines of Chinese cabbage, namely "G291" (a hairy line) and "ZHB" (a hairless line), based on which a candidate gene related to trichome traits was identified on chromosome A06 with a phenotypic variance of 47%. A molecular marker was found to co-segregate with the trichome traits of the F₂ population, which is in the 5'-flanking region of BrGL1, and a corresponding patent has been granted (NO. CN 108545775 B). Transcriptome analysis was carried out on the cotyledon, the first true leaf and the leaf closest to each inflorescence of F₂ individuals of "G291 × ZHB" with or without trichomes, respectively. Ten pathways, including 189 DEGs, were identified to be involved in the development of trichomes in Chinese cabbage, which may be specifically related to the development of leaf trichomes. Most of the pathways were related to the biosynthesis of the secondary metabolites, which may help plants to adapt to the ever-changing external environment. DEGs also enriched the "plant-pathogen interaction" pathway, which is consistent with the conclusion that trichomes are related to the disease resistance of plants. Our study provides a basis for future research on the occurrence and development of trichomes in Chinese cabbage.

Genetic Relationship of Brassicaceae Hybrids with Various Resistance to Blackleg Is Disclosed by the Use of Molecular Markers

Brassica napus is an important oil source. Its narrow gene pool can be widened by interspecific hybridization with the Brassicaceae species. One of the agronomically important traits, that can be transferred through the hybridization, is the resistance to blackleg, a dangerous disease mainly caused by Leptosphaeria maculans. Hybrid individuals can be analyzed with various molecular markers, including Simple Sequence Repeats (SSR). We investigated the genetic similarity of 32 Brassicaceae hybrids and 19 parental components using SSR markers to reveal their genetic relationship. Furthermore, we compared the field resistance to blackleg of the interspecific progenies. The tested set of 15 SSR markers proved to be useful in revealing the genetic distances in the Brassicaceae hybrids and species. However, genetic similarity of the studied hybrids could not be correlated with the level of field resistance to L. maculans. Moreover, our studies confirmed the usefulness of the Brassicaceae hybrids in terms of blackleg management.

Construction of an Intragenic SSR-Based Linkage Map and QTL Mapping for Agronomic Traits in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Chinese cabbage (Brassica rapa L. ssp. pekinensis) is one of the most widely cultivated and economically important vegetables in China. Constructing an effective genetic linkage map and mapping quantitative trait loci (QTLs) related to yield and leafy head morphology is of great importance for molecular breeding of Chinese cabbage. Using two diverse Chinese cabbage inbred lines, ZHB and G291, as parents, an F2 segregating population consisting of 240 individuals was prepared for genetic map construction and phenotype investigation in this study. The two parents are significantly different in both shape and size. Sixteen important agronomic traits of F2 individuals were investigated. A genetic map of 105 intragenic simple sequence repeat (SSR) markers distributed across 10 linkage groups (LGs) was constructed, which was 2034.1 cM in length and had an average inter-locus distance of 21.75 cM. We identified 48 QTLs for the tested important agronomic traits on the studied LGs, with LOD scores of 2.51–12.49, which explained the phenotypic variance of 3.41–26.66%. The QTLs identified in this study will facilitate further genetic analysis and marker-assisted genetic improvement of Chinese cabbage.

QTL mapping for Fusarium wilt resistance based on the whole-genome resequencing and their association with functional genes in Raphanus sativus

Two major QTL associated with resistance to Fusarium wilt (FW) were identified using whole-genome resequencing. Sequence variations and gene expression level differences suggest that TIR-NBS and LRR-RLK are candidate genes associated with FW-resistance. Fusarium wilt (FW) caused by Fusarium oxysporum f. sp. raphani is an important disease in radish, leading to severe decrease in yield and quality. YR4 as a novel genetic source to resistant to FW was confirmed through screening with five pathogen isolates. We have generated F₂ and F_2:3 populations segregated with FW resistance using YR4 and YR18 inbred lines. The disease symptom was evaluated in F_2:3 population (n = 180) in three independent studies over two years. We identified 4 QTL including the two major QTL (FoRsR7.1^59A and FoRsR9.3^59A). FoRsR7.1^59A and FoRsR9.3^59A were detected in three replicated experiments. FoRsR7.1^59A was delimited to the 2.18-Mb physical interval on chromosome R07, with a high LOD value (5.17-12.84) and explained phenotypic variation (9.34%-27.97%). The FoRsR9.3^59A represented relatively low LOD value (3.38-4.52) and explained phenotypic variation (6.24%-8.82%). On the basis of the re-sequencing data for the parental lines, we identified five putative resistance-related genes and 13 unknown genes with sequence variations at the gene and protein levels. A semi-quantitative RT-PCR analysis revealed that Rs382940 (TIR-NBS) and Rs382200 (RLK) were expressed only in 'YR4' from 0 to 6 days after the inoculation. Moreover, Rs382950 (TIR-NBS-LRR) was more highly expressed in 'YR4' from 3 to 6 days after the inoculation. These three genes might be important for FW-resistance in radish. We identified several markers based on these potential candidate genes. The marker set should be useful for breeding system to introduce the FW resistance loci from 'YR4' to improve tolerance to FW.

Genotyping-by-Sequencing-Based Genome-Wide Association Studies of Fusarium Wilt Resistance in Radishes (Raphanus sativus L.)

Fusarium wilt (FW) is a fungal disease that causes severe yield losses in radish production. The most effective method to control the FW is the development and use of resistant varieties in cultivation. The identification of marker loci linked to FW resistance are expected to facilitate the breeding of disease-resistant radishes. In the present study, we applied an integrated framework of genome-wide association studies (GWAS) using genotyping-by-sequencing (GBS) to identify FW resistance loci among a panel of 225 radish accessions, including 58 elite breeding lines. Phenotyping was conducted by manual inoculation of seedlings with the FW pathogen, and scoring for the disease index was conducted three weeks after inoculation during two constitutive years. The GWAS analysis identified 44 single nucleotide polymorphisms (SNPs) and twenty putative candidate genes that were significantly associated with FW resistance. In addition, a total of four QTLs were identified from F2 population derived from a FW resistant line and a susceptible line, one of which was co-located with the SNPs on chromosome 7, detected in GWAS study. These markers will be valuable for molecular breeding programs and marker-assisted selection to develop FW resistant varieties of R. sativus.

Fine-mapping of a major QTL (Fwr1) for fusarium wilt resistance in radish

A major radish QTL (Fwr1) for fusarium wilt resistance was fine-mapped. Sequence and expression analyses suggest that a gene encoding a serine/arginine-rich protein kinase is a candidate gene for Fwr1. Fusarium wilt resistance locus 1 (Fwr1) is a major quantitative trait locus (QTL) mediating the resistance of radish inbred line 'B2' to Fusarium oxysporum, which is responsible for fusarium wilt. We previously detected Fwr1 on radish linkage group 3 (i.e., chromosome 5). In this study, a high-resolution genetic map of the Fwr1 locus was constructed by analyzing 354 recombinant F₂ plants derived from a cross between 'B2' and '835', the latter of which is susceptible to fusarium wilt. The Fwr1 QTL was fine-mapped to a 139.8-kb region between markers FM82 and FM87 in the middle part of chromosome 5. Fifteen candidate genes were predicted in this region based on a sequence comparison with the 'WK10039' radish reference genome. Additionally, we examined the time-course expression patterns of these predicted genes following an infection by the fusarium wilt pathogen. The ORF4 expression level was significantly higher in the resistant 'B2' plants than in the susceptible '835' plants. The ORF4 sequence was predicted to encode a serine/arginine-rich protein kinase and includes SNPs that result in nonsynonymous mutations, which may have important functional consequences. This study reveals a novel gene responsible for fusarium wilt resistance in radish. Further analyses of this gene may elucidate the molecular mechanisms underlying the fusarium wilt resistance of plants.

Construction of genetic linkage map and identification of a novel major locus for resistance to pine wood nematode in Japanese black pine (Pinus thunbergii)

BACKGROUND

Pine wilt disease (PWD), which is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, is currently the greatest threat to pine forests in Europe and East Asian countries including Japan. Constructing a detailed linkage map of DNA markers and identifying PWD resistance genes/loci lead to improved resistance in Pinus thunbergii, as well as other Pinus species that are also susceptible to PWD.

RESULTS

A total F₁ mapping population of 188 individuals derived from a cross between the PWD-resistant P. thunbergii varieties 'Tanabe 54' (resistant rank 2 to PWD) and 'Tosashimizu 63' (resistant rank 4 to PWD) was inoculated with PWN, and was evaluated for disease symptoms. To perform linkage analysis for PWN resistance, a set of three maps was constructed; two parental maps generated using the integrated two-way pseudo-testcross method, and a consensus map with population-type cross-pollination. The linkage map of 'Tanabe 54' consisted of 167 loci, and covered 14 linkage groups (LGs), with a total genetic distance of 1214.6 cM. The linkage map of 'Tosashimizu 63' consisted of 252 loci, and covered 14 LGs, with a total genetic distance of 1422.1 cM. The integrated consensus map comprised 12 LGs with the basic chromosome number of P. thunbergii, and a total genetic distance of 1403.6 cM. Results from quantitative trait loci (QTL) analysis using phenotype data and linkage maps indicated that PWN resistance is controlled by a single dominant allele, which was derived from the 'Tanabe 54' female parent. This major QTL was located on linkage group 3 and was designated PWD1 for PINE WILT DISEASE 1.

CONCLUSIONS

The PWD1 locus is a major resistance QTL located on the Pinus consensus LG03 that acts in a dominant manner to confer pine wood nematode resistance. Information from the present study will be useful for P. thunbergii breeding programs to improve resistance to PWD, and also to help identify susceptibility genes in Pinus species.

The radish genome database (RadishGD): an integrated information resource for radish genomics

Radish (Raphanus sativus L.) is an important root vegetable crop in the family Brassicaceae, which provides diverse nutrients for human health and is closely related to the Brassica crop species. Recently, we sequenced and assembled the radish genome into nine chromosome pseudomolecules. In addition, we developed diverse genomic resources, including genetic maps, molecular markers, transcriptome, genome-wide methylation and variome data. In this study, we describe the radish genome database (RadishGD), including details of data sets that we generated and the web interface that allows access to these data. RadishGD comprises six major units that enable researchers and general users to search, browse and analyze the radish genomic data in an integrated manner. The Search unit provides gene structures and sequences for gene models through keyword or BLAST searches. The Genome browser displays graphic representations of gene models, mRNAs, repetitive sequences, genome-wide methylation and variomes among various genotypes. The Functional annotation unit offers gene ontology, plant ontology, pathway and gene family information for gene models. The Genetic map unit provides information about markers and their genetic locations using two types of genetic maps. The Expression unit presents transcriptional characteristics and methylation levels for each gene in 18 tissues. All sequence data incorporated into RadishGD can be downloaded from the Data resources unit. RadishGD will be continually updated to serve as a community resource for radish genomics and breeding research.

Transcriptome sequencing assisted discovery and computational analysis of novel SNPs associated with flowering in Raphanus sativus in-bred lines for marker-assisted backcross breeding

The sequencing of radish genome aids in the better understanding and tailoring of traits associated with economic importance. In order to accelerate the genomics assisted breeding and genetic selection, transcriptomes of 33 radish inbred lines with diverse traits were sequenced for the development of single nucleotide polymorphic (SNP) markers. The sequence reads ranged from 2,560,543,741 bp to 20,039,688,139 bp with the GC (%) of 47.80-49.34 and phred quality score (Q30) of 96.47-97.54%. A total of 4951 polymorphic SNPs were identified among the accessions after stringent filtering and 298 SNPs with efficient marker assisted backcross breeding (MAB) markers were generated from the polymorphic SNPs. Further, functional annotations of SNPs revealed the effects and importance of the SNPs identified in the flowering process. The SNPs were predominantly associated with the four major flowering related transcription factors such as MYB, MADS box (AG), AP2/EREB, and bHLH. In addition, SNPs in the vital flowering integrator gene (FT) and floral repressors (EMBRYONIC FLOWER 1, 2, and FRIGIDA) were identified among the radish inbred lines. Further, 50 SNPs were randomly selected from 298 SNPs and validated using Kompetitive Allele Specific PCR genotyping system (KASP) in 102 radish inbred lines. The homozygosity of the inbred lines varied from 56 to 96% and the phylogenetic analysis resulted in the clustering of inbred lines into three subgroups. Taken together, the SNP markers identified in the present study can be utilized for the discrimination, seed purity test, and adjusting parental combinations for breeding in radish.

Gene Introgression in Weeds Depends on Initial Gene Location in the Crop: Brassica napus-Raphanus raphanistrum Model

The effect of gene location within a crop genome on its transfer to a weed genome remains an open question for gene flow assessment. To elucidate this question, we analyzed advanced generations of intergeneric hybrids, derived from an initial pollination of known oilseed rape varieties (Brassica napus, AACC, 2n = 38) by a local population of wild radish (Raphanus raphanistrum, RrRr, 2n = 18). After five generations of recurrent pollination, 307 G5 plants with a chromosome number similar to wild radish were genotyped using 105 B. napus specific markers well distributed along the chromosomes. They revealed that 49.8% of G5 plants carried at least one B. napus genomic region. According to the frequency of B. napus markers (0-28%), four classes were defined: Class 1 (near zero frequency), with 75 markers covering ∼70% of oilseed rape genome; Class 2 (low frequency), with 20 markers located on 11 genomic regions; Class 3 (high frequency), with eight markers on three genomic regions; and Class 4 (higher frequency), with two adjacent markers detected on A10. Therefore, some regions of the oilseed rape genome are more prone than others to be introgressed into wild radish. Inheritance and growth of plant progeny revealed that genomic regions of oilseed rape could be stably introduced into wild radish and variably impact the plant fitness (plant height and seed number). Our results pinpoint that novel technologies enabling the targeted insertion of transgenes should select genomic regions that are less likely to be introgressed into the weed genome, thereby reducing gene flow.

High-throughput development of simple sequence repeat markers for genetic diversity research in Crambe abyssinica

BACKGROUND

The allohexaploid Crambe abyssinica (crambe) is an oilseed crop that has been recognized for its potential value in the chemical industry, particularly in terms of producing high-erucic acid content vegetable oil. However, as an understudied crop, improvement of crambe has been hampered by the lack of genetic and genomic information to enhance its yield, oil quality and resistance against biotic and abiotic stress. Development of molecular markers is therefore of great significance to facilitate genetic improvement of crambe.

RESULTS

In this study, high-throughput sequencing was performed to generate sequences for the transcriptome and genome of a widely planted crambe cultivar, Galactica. A total of 186,778 expressed sequence tag (EST) contigs as 8,130,350 genomic contigs were assembled as well. Altogether, 82,523 pairs of primers were designed in the flanking sequences of the simple sequence repeat (SSR) within these contigs. Virtual PCR analysis showed that a fraction of these primers could be mapped onto the genomes of related species of Brassica, including Brassica rapa, B. oleraceae and B. napus. Genetic diversity analysis using a subset of 166 markers on 30 independent C. abyssinica accessions exhibited that 1) 95 % of the designed SSRs were polymorphic among these accessions; 2) the polymorphism information content (PIC) value of the markers ranged from 0.13 to 0.89; 3) the genetic distances (coefficient NEI72) between accessions varied from 0.06 to 0.36. Cluster analysis subsequent on the accessions demonstrated consistency with crambe breeding history. F-statistics analysis revealed a moderate level of genetic differentiation in C. abyssinica (Gst = 0.3934) and a accordingly low estimated gene flow (Nm = 0.7709).

CONCLUSION

Application of high-throughput sequencing technology has facilitated SSR marker development, which was successfully employed in evaluating genetic diversity of C. abyssinica as demonstrated in our study. Results showed these molecular markers were robust and provided powerful tools for assessing genetic diversity and estimating crambe breeding history. Moreover, the SSR primers and sequence information developed in the study are freely available to the research community.

Characterizing the transcriptome and molecular markers information for roach, Rutilus rutilus

Rutilus rutilus (roach) is native to most of Europe and western Asia, and the Irtysh River basin in Sinkiang, northwest China is the marginal area of their natural distribution. The wide distribution and unique characteristic of this species makes it an ideal model for analysing ecological and comparative genomics. However, the limited genome sequences available for this species have hindered these investigations. Transcriptomes from the brains and livers of five individuals collected from the Irtysh River basin were sequenced using Illumina paired-end sequencing technology. A collection of 132,289 unigenes for this species were obtained using a de novo assembly method based on nearly 120 million clean reads encompassing more than 14 Gb data. Approximately 37.5% (49,656), 27.1% (35,867) and 21.2% (27,987) of the transcriptome had homologues deposited in Nt, Nr and Swiss-Prot, respectively; 12.3% (16,328) were assigned to eukaryotic orthologous groups of proteins classifications, and 21.5% (28,429) harboured Interpro domains. On the basis of the assembled transcriptome, we detected 177,493 single-nuc leotide variation resident in 39.3% (52,029) of the sequences and 20.8% (27,497) of the sequences harbouring 36,639 simple sequence repeats. The identified molecular markers are a basis for further ecological analysis, and the transcriptome reported here allows for more extensive evolutionary analyses of the Cyprinidae, the most species-rich family of freshwater fishes.

Identification of Flowering-Related Genes Responsible for Differences in Bolting Time between Two Radish Inbred Lines

Late bolting after cold exposure is an economically important characteristic of radish (Raphanus sativus L.), an important Brassicaceae root vegetable crop. However, little information is available regarding the genes and pathways that govern flowering time in this species. We performed high-throughput RNA sequencing analysis to elucidate the molecular mechanisms that determine the differences in flowering times between two radish lines, NH-JS1 (late bolting) and NH-JS2 (early bolting). In total, 71,188 unigenes were identified by reference-guided assembly, of which 309, 788, and 980 genes were differentially expressed between the two inbred lines after 0, 15, and 35 days of vernalization, respectively. Among these genes, 218 homologs of Arabidopsis flowering-time (Ft) genes were identified in the radish, and 49 of these genes were differentially expressed between the two radish lines in the presence or absence of vernalization treatment. Most of the Ft genes up-regulated in NH-JS1 vs. NH-JS2 were repressors of flowering, such as RsFLC, consistent with the late-bolting phenotype of NH-JS1. Although, the functions of genes down-regulated in NH-JS1 were less consistent with late-bolting characteristics than the up-regulated Ft genes, several Ft enhancer genes, including RsSOC1, a key floral integrator, showed an appropriate expression to the late-bolting phenotype. In addition, the patterns of gene expression related to the vernalization pathway closely corresponded with the different bolting times of the two inbred lines. These results suggest that the vernalization pathway is conserved between radish and Arabidopsis.

Quantitative Trait Loci for Morphological Traits and their Association with Functional Genes in Raphanus sativus

Identification of quantitative trait loci (QTLs) governing morphologically important traits enables to comprehend their potential genetic mechanisms in the genetic breeding program. In this study, we used 210 F2 populations derived from a cross between two radish inbred lines (Raphanus sativus) "835" and "B2," including 258 SSR markers were used to detect QTLs for 11 morphological traits that related to whole plant, leaf, and root yield in 3 years of replicated field test. Total 55 QTLs were detected which were distributed on each linkage group of the Raphanus genome. Individual QTLs accounted for 2.69-12.6 of the LOD value, and 0.82-16.25% of phenotypic variation. Several genomic regions have multiple traits that clustered together, suggested the existence of pleiotropy linkage. Synteny analysis of the QTL regions with A. thaliana genome selected orthologous genes in radish. InDels and SNPs in the parental lines were detected in those regions by Illumina genome sequence. Five identified candidate gene-based markers were validated by co-mapping with underlying QTLs affecting different traits. Semi-quantitative reverse transcriptase PCR analysis showed the different expression levels of these five genes in parental lines. In addition, comparative QTL analysis with B. rapa revealed six common QTL regions and four key major evolutionarily conserved crucifer blocks (J, U, R, and W) harboring QTL for morphological traits. The QTL positions identified in this study will provide a valuable resource for identifying more functional genes when whole radish genome sequence is released. Candidate genes identified in this study that co-localized in QTL regions are expected to facilitate in radish breeding programs.

The radish genome and comprehensive gene expression profile of tuberous root formation and development

Understanding the processes that regulate plant sink formation and development at the molecular level will contribute to the areas of crop breeding, food production and plant evolutionary studies. We report the annotation and analysis of the draft genome sequence of the radish Raphanus sativus var. hortensis (long and thick root radish) and transcriptome analysis during root development. Based on the hybrid assembly approach of next-generation sequencing, a total of 383 Mb (N50 scaffold: 138.17 kb) of sequences of the radish genome was constructed containing 54,357 genes. Syntenic and phylogenetic analyses indicated that divergence between Raphanus and Brassica coincide with the time of whole genome triplication (WGT), suggesting that WGT triggered diversification of Brassiceae crop plants. Further transcriptome analysis showed that the gene functions and pathways related to carbohydrate metabolism were prominently activated in thickening roots, particularly in cell proliferating tissues. Notably, the expression levels of sucrose synthase 1 (SUS1) were correlated with root thickening rates. We also identified the genes involved in pungency synthesis and their transcription factors.

Construction of a reference genetic map of Raphanus sativus based on genotyping by whole-genome resequencing

This manuscript provides a genetic map of Raphanus sativus that has been used as a reference genetic map for an ongoing genome sequencing project. The map was constructed based on genotyping by whole-genome resequencing of mapping parents and F 2 population. Raphanus sativus is an annual vegetable crop species of the Brassicaceae family and is one of the key plants in the seed industry, especially in East Asia. Assessment of the R. sativus genome provides fundamental resources for crop improvement as well as the study of crop genome structure and evolution. With the goal of anchoring genome sequence assemblies of R. sativus cv. WK10039 whose genome has been sequenced onto the chromosomes, we developed a reference genetic map based on genotyping of two parents (maternal WK10039 and paternal WK10024) and 93 individuals of the F2 mapping population by whole-genome resequencing. To develop high-confidence genetic markers, ~83 Gb of parental lines and ~591 Gb of mapping population data were generated as Illumina 100 bp paired-end reads. High stringent sequence analysis of the reads mapped to the 344 Mb of genome sequence scaffolds identified a total of 16,282 SNPs and 150 PCR-based markers. Using a subset of the markers, a high-density genetic map was constructed from the analysis of 2,637 markers spanning 1,538 cM with 1,000 unique framework loci. The genetic markers integrated 295 Mb of genome sequences to the cytogenetically defined chromosome arms. Comparative analysis of the chromosome-anchored sequences with Arabidopsis thaliana and Brassica rapa revealed that the R. sativus genome has evident triplicated sub-genome blocks and the structure of gene space is highly similar to that of B. rapa. The genetic map developed in this study will serve as fundamental genomic resources for the study of R. sativus.

In vivo functional analysis of a nuclear restorer PPR protein

BACKGROUND

Nuclear restorers of cytoplasmic male fertility (CMS) act to suppress the male sterile phenotype by down-regulating the expression of novel CMS-specifying mitochondrial genes. One such restorer gene is Rfo, which restores fertility to the radish Ogura or ogu CMS. Rfo, like most characterized restorers, encodes a pentatricopeptide repeat (PPR) protein, a family of eukaryotic proteins characterized by tandem repeats of a 35 amino acid motif. While over 400 PPR genes are found in characterized plant genomes and the importance of this gene family in organelle gene expression is widely recognized, few detailed in vivo assessments of primary structure-function relationships in this protein family have been conducted.

RESULTS

In contrast to earlier studies, which identified 16 or 17 PPR domains in the Rfo protein, we now find, using a more recently developed predictive tool, that Rfo has 18 repeat domains with the additional domain N-terminal to the others. Comparison of transcript sequences from pooled rfo/rfo plants with pooled Rfo/Rfo plants of a mapping population led to the identification of a non-restoring rfo allele with a 12 bp deletion in the fourth domain. Introduction into ogu CMS plants of a genetic construct in which this deletion had been introduced into Rfo led to a partial loss in the capacity to produce viable pollen, as assessed by vital staining, pollen germination and the capacity for seed production following pollination of CMS plants. The degree of viable pollen production among different transgenic plants roughly correlated with the copy number of the introduced gene and with the reduction of the levels of the ORF138 CMS-associated protein. All other constructs tested, including one in which only the C-terminal PPR repeat was deleted and another in which this repeat was replaced by the corresponding domain of the related, non-restoring gene, PPR-A, failed to result in any measure of fertility restoration.

CONCLUSIONS

The identification of the additional PPR domain in Rfo indicates that the protein, apart from its N-terminal mitochondrial targeting presequence, consists almost entirely of PPR repeats. The newly identified rfo allele carries the same 4 amino acid deletion as that found in the neighboring, related, non-restoring PPR gene, PPR-A. Introduction of this four amino acid deletion into a central domain the Rfo protein, however, only partially reduces its restoration capacity, even though this alteration might be expected to alter the spacing between the adjoining repeats. All other tested alterations, generated by deleting specific PPR repeats or exchanging repeats with corresponding domains of PPR-A, led to a complete loss of restorer function. Overall we demonstrate that introduction of targeted alterations of Rfo into ogu CMS plants provides a sensitive in vivo readout for analysis of the relationship between primary structure and biological function in this important family of plant proteins.

A complex recombination pattern in the genome of allotetraploid Brassica napus as revealed by a high-density genetic map

Polyploidy plays a crucial role in plant evolution. Brassica napus (2n = 38, AACC), the most important oil crop in the Brassica genus, is an allotetraploid that originated through natural doubling of chromosomes after the hybridization of its progenitor species, B. rapa (2n = 20, AA) and B. oleracea (2n = 18, CC). A better understanding of the evolutionary relationship between B. napus and B. rapa, B. oleracea, as well as Arabidopsis, which has a common ancestor with these three species, will provide valuable information about the generation and evolution of allopolyploidy. Based on a high-density genetic map with single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers, we performed a comparative genomic analysis of B. napus with Arabidopsis and its progenitor species B. rapa and B. oleracea. Based on the collinear relationship of B. rapa and B. oleracea in the B. napus genetic map, the B. napus genome was found to consist of 70.1% of the skeleton components of the chromosomes of B. rapa and B. oleracea, with 17.7% of sequences derived from reciprocal translocation between homoeologous chromosomes between the A- and C-genome and 3.6% of sequences derived from reciprocal translocation between non-homologous chromosomes at both intra- and inter-genomic levels. The current study thus provides insights into the formation and evolution of the allotetraploid B. napus genome, which will allow for more accurate transfer of genomic information from B. rapa, B. oleracea and Arabidopsis to B. napus.

Kazusa Marker DataBase: a database for genomics, genetics, and molecular breeding in plants

In order to provide useful genomic information for agronomical plants, we have established a database, the Kazusa Marker DataBase (http://marker.kazusa.or.jp). This database includes information on DNA markers, e.g., SSR and SNP markers, genetic linkage maps, and physical maps, that were developed at the Kazusa DNA Research Institute. Keyword searches for the markers, sequence data used for marker development, and experimental conditions are also available through this database. Currently, 10 plant species have been targeted: tomato (Solanum lycopersicum), pepper (Capsicum annuum), strawberry (Fragaria × ananassa), radish (Raphanus sativus), Lotus japonicus, soybean (Glycine max), peanut (Arachis hypogaea), red clover (Trifolium pratense), white clover (Trifolium repens), and eucalyptus (Eucalyptus camaldulensis). In addition, the number of plant species registered in this database will be increased as our research progresses. The Kazusa Marker DataBase will be a useful tool for both basic and applied sciences, such as genomics, genetics, and molecular breeding in crops.

Comparative analysis of the radish genome based on a conserved ortholog set (COS) of Brassica

This manuscript provides a Brassica conserved ortholog set (COS) that can be used as diagnostic cross-species markers as well as tools for genetic mapping and genome comparison of the Brassicaceae. A conserved ortholog set (COS) is a collection of genes that are conserved in both sequence and copy number between closely related genomes. COS is a useful resource for developing gene-based markers and is suitable for comparative genome mapping. We developed a COS for Brassica based on proteome comparisons of Arabidopsis thaliana, B. rapa, and B. oleracea to establish a basis for comparative genome analysis of crop species in the Brassicaceae. A total of 1,194 conserved orthologous single-copy genes were identified from the genomes based on whole-genome BLASTP analysis. Gene ontology analysis showed that most of them encoded proteins with unknown function and chloroplast-related genes were enriched. In addition, 152 Brassica COS primer sets were applied to 16 crop and wild species of the Brassicaceae and 57.9-92.8 % of them were successfully amplified across the species representing that a Brassica COS can be used as diagnostic cross-species markers of diverse Brassica species. We constructed a genetic map of Raphanus sativus by analyzing the segregation of 322 COS genes in an F2 population (93 individuals) of Korean cultivars (WK10039 × WK10024). Comparative genome analysis based on the COS genes showed conserved genome structures between R. sativus and B. rapa with lineage-specific rearrangement and fractionation of triplicated subgenome blocks indicating close evolutionary relationship and differentiation of the genomes. The Brassica COS developed in this study will play an important role in genetic, genomic, and breeding studies of crop Brassicaceae species.

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.

Comparative genomics of Brassicaceae crops

The family Brassicaceae is one of the major groups of the plant kingdom and comprises diverse species of great economic, agronomic and scientific importance, including the model plant Arabidopsis. The sequencing of the Arabidopsis genome has revolutionized our knowledge in the field of plant biology and provides a foundation in genomics and comparative biology. Genomic resources have been utilized in Brassica for diversity analyses, construction of genetic maps and identification of agronomic traits. In Brassicaceae, comparative sequence analysis across the species has been utilized to understand genome structure, evolution and the detection of conserved genomic segments. In this review, we focus on the progress made in genetic resource development, genome sequencing and comparative mapping in Brassica and related species. The utilization of genomic resources and next-generation sequencing approaches in improvement of Brassica crops is also discussed.

Transcriptome sequencing of two parental lines of cabbage (Brassica oleracea L. var. capitata L.) and construction of an EST-based genetic map

Background

Expressed sequence tag (EST)-based markers are preferred because they reflect transcribed portions of the genome. We report the development of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers derived from transcriptome sequences in cabbage, and their utility for map construction.

Results

Transcriptome sequences were obtained from two cabbage parental lines, C1184 and C1234, which are susceptible and resistant to black rot disease, respectively, using the 454 platform. A total of 92,255 and 127,522 reads were generated and clustered into 34,688 and 40,947 unigenes, respectively. We identified 2,405 SSR motifs from the unigenes of the black rot-resistant parent C1234. Trinucleotide motifs were the most abundant (66.15%) among the repeat motifs. In addition, 1,167 SNPs were detected between the two parental lines. A total of 937 EST-based SSR and 97 SNP-based dCAPS markers were designed and used for detection of polymorphism between parents. Using an F₂ population, we built a genetic map comprising 265 loci, and consisting of 98 EST-based SSRs, 21 SNP-based dCAPS, 55 IBP markers derived from B. rapa genome sequence and 91 public SSRs, distributed on nine linkage groups spanning a total of 1,331.88 cM with an average distance of 5.03 cM between adjacent loci. The parental lines used in this study are elite breeding lines with little genetic diversity; therefore, the markers that mapped in our genetic map will have broad spectrum utility.

Conclusions

This genetic map provides additional genetic information to the existing B. oleracea map. Moreover, the new set of EST-based SSR and dCAPS markers developed herein is a valuable resource for genetic studies and will facilitate cabbage breeding. Additionally, this study demonstrates the usefulness of NGS transcriptomes for the development of genetic maps even with little genetic diversity in the mapping population.

Plant Genome DataBase Japan (PGDBj): a portal website for the integration of plant genome-related databases

The Plant Genome DataBase Japan (PGDBj, http://pgdbj.jp/?ln=en) is a portal website that aims to integrate plant genome-related information from databases (DBs) and the literature. The PGDBj is comprised of three component DBs and a cross-search engine, which provides a seamless search over the contents of the DBs. The three DBs are as follows. (i) The Ortholog DB, providing gene cluster information based on the amino acid sequence similarity. Over 500,000 amino acid sequences of 20 Viridiplantae species were subjected to reciprocal BLAST searches and clustered. Sequences from plant genome DBs (e.g. TAIR10 and RAP-DB) were also included in the cluster with a direct link to the original DB. (ii) The Plant Resource DB, integrating the SABRE DB, which provides cDNA and genome sequence resources accumulated and maintained in the RIKEN BioResource Center and National BioResource Projects. (iii) The DNA Marker DB, providing manually or automatically curated information of DNA markers, quantitative trait loci and related linkage maps, from the literature and external DBs. As the PGDBj targets various plant species, including model plants, algae, and crops important as food, fodder and biofuel, researchers in the field of basic biology as well as a wide range of agronomic fields are encouraged to perform searches using DNA sequences, gene names, traits and phenotypes of interest. The PGDBj will return the search results from the component DBs and various types of linked external DBs.

Comprehensive analysis of expressed sequence tags from cultivated and wild radish (Raphanus spp.)

BACKGROUND

Radish (Raphanus sativus L., 2n = 2× = 18) is an economically important vegetable crop worldwide. A large collection of radish expressed sequence tags (ESTs) has been generated but remains largely uncharacterized.

RESULTS

In this study, approximately 315,000 ESTs derived from 22 Raphanus cDNA libraries from 18 different genotypes were analyzed, for the purpose of gene and marker discovery and to evaluate large-scale genome duplication and phylogenetic relationships among Raphanus spp. The ESTs were assembled into 85,083 unigenes, of which 90%, 65%, 89% and 89% had homologous sequences in the GenBank nr, SwissProt, TrEMBL and Arabidopsis protein databases, respectively. A total of 66,194 (78%) could be assigned at least one gene ontology (GO) term. Comparative analysis identified 5,595 gene families unique to radish that were significantly enriched with genes related to small molecule metabolism, as well as 12,899 specific to the Brassicaceae that were enriched with genes related to seed oil body biogenesis and responses to phytohormones. The analysis further indicated that the divergence of radish and Brassica rapa occurred approximately 8.9-14.9 million years ago (MYA), following a whole-genome duplication event (12.8-21.4 MYA) in their common ancestor. An additional whole-genome duplication event in radish occurred at 5.1-8.4 MYA, after its divergence from B. rapa. A total of 13,570 simple sequence repeats (SSRs) and 28,758 high-quality single nucleotide polymorphisms (SNPs) were also identified. Using a subset of SNPs, the phylogenetic relationships of eight different accessions of Raphanus was inferred.

CONCLUSION

Comprehensive analysis of radish ESTs provided new insights into radish genome evolution and the phylogenetic relationships of different radish accessions. Moreover, the radish EST sequences and the associated SSR and SNP markers described in this study represent a valuable resource for radish functional genomics studies and breeding.

Comparative mapping of Raphanus sativus genome using Brassica markers and quantitative trait loci analysis for the Fusarium wilt resistance trait

Fusarium wilt (FW), caused by the soil-borne fungal pathogen Fusarium oxysporum is a serious disease in cruciferous plants, including the radish (Raphanus sativus). To identify quantitative trait loci (QTL) or gene(s) conferring resistance to FW, we constructed a genetic map of R. sativus using an F2 mapping population derived by crossing the inbred lines '835' (susceptible) and 'B2' (resistant). A total of 220 markers distributed in 9 linkage groups (LGs) were mapped in the Raphanus genome, covering a distance of 1,041.5 cM with an average distance between adjacent markers of 4.7 cM. Comparative analysis of the R. sativus genome with that of Arabidopsis thaliana and Brassica rapa revealed 21 and 22 conserved syntenic regions, respectively. QTL mapping detected a total of 8 loci conferring FW resistance that were distributed on 4 LGs, namely, 2, 3, 6, and 7 of the Raphanus genome. Of the detected QTL, 3 QTLs (2 on LG 3 and 1 on LG 7) were constitutively detected throughout the 2-year experiment. QTL analysis of LG 3, flanked by ACMP0609 and cnu_mBRPGM0085, showed a comparatively higher logarithm of the odds (LOD) value and percentage of phenotypic variation. Synteny analysis using the linked markers to this QTL showed homology to A. thaliana chromosome 3, which contains disease-resistance gene clusters, suggesting conservation of resistance genes between them.

Construction of a chromosome-assigned, sequence-tagged linkage map for the radish, Raphanus sativus L. and QTL analysis of morphological traits

The radish displays great morphological variation but the genetic factors underlying this variability are mostly unknown. To identify quantitative trait loci (QTLs) controlling radish morphological traits, we cultivated 94 F4 and F5 recombinant inbred lines derived from a cross between the rat-tail radish and the Japanese radish cultivar 'Harufuku' inbred lines. Eight morphological traits (ovule and seed numbers per silique, plant shape, pubescence and root formation) were measured for investigation. We constructed a map composed of 322 markers with a total length of 673.6 cM. The linkage groups were assigned to the radish chromosomes using disomic rape-radish chromosome-addition lines. On the map, eight and 10 QTLs were identified in 2008 and 2009, respectively. The chromosome-linkage group correspondence, the sequence-specific markers and the QTLs detected here will provide useful information for further genetic studies and for selection during radish breeding programs.

Deciphering the diploid ancestral genome of the Mesohexaploid Brassica rapa

The genus Brassica includes several important agricultural and horticultural crops. Their current genome structures were shaped by whole-genome triplication followed by extensive diploidization. The availability of several crucifer genome sequences, especially that of Chinese cabbage (Brassica rapa), enables study of the evolution of the mesohexaploid Brassica genomes from their diploid progenitors. We reconstructed three ancestral subgenomes of B. rapa (n = 10) by comparing its whole-genome sequence to ancestral and extant Brassicaceae genomes. All three B. rapa paleogenomes apparently consisted of seven chromosomes, similar to the ancestral translocation Proto-Calepineae Karyotype (tPCK; n = 7), which is the evolutionarily younger variant of the Proto-Calepineae Karyotype (n = 7). Based on comparative analysis of genome sequences or linkage maps of Brassica oleracea, Brassica nigra, radish (Raphanus sativus), and other closely related species, we propose a two-step merging of three tPCK-like genomes to form the hexaploid ancestor of the tribe Brassiceae with 42 chromosomes. Subsequent diversification of the Brassiceae was marked by extensive genome reshuffling and chromosome number reduction mediated by translocation events and followed by loss and/or inactivation of centromeres. Furthermore, via interspecies genome comparison, we refined intervals for seven of the genomic blocks of the Ancestral Crucifer Karyotype (n = 8), thus revising the key reference genome for evolutionary genomics of crucifers.

RadishBase: a database for genomics and genetics of radish

Radish is an economically important vegetable crop. During the past several years, large-scale genomics and genetics resources have been accumulated for this species. To store, query, analyze and integrate these radish resources efficiently, we have developed RadishBase (http://bioinfo.bti.cornell.edu/radish), a genomics and genetics database of radish. Currently the database contains radish mitochondrial genome sequences, expressed sequence tag (EST) and unigene sequences and annotations, biochemical pathways, EST-derived single nucleotide polymorphism (SNP) and simple sequence repeat (SSR) markers, and genetic maps. RadishBase is designed to enable users easily to retrieve and visualize biologically important information through a set of efficient query interfaces and analysis tools, including the BLAST search and unigene annotation query interfaces, and tools to classify unigenes functionally, to identify enriched gene ontology (GO) terms and to visualize genetic maps. A database containing radish pathways predicted from unigene sequences is also included in RadishBase. The tools and interfaces in RadishBase allow efficient mining of recently released and continually expanding large-scale radish genomics and genetics data sets, including the radish genome sequences and RNA-seq data sets.

QTL analysis using SNP markers developed by next-generation sequencing for identification of candidate genes controlling 4-methylthio-3-butenyl glucosinolate contents in roots of radish, Raphanus sativus L

SNP markers for QTL analysis of 4-MTB-GSL contents in radish roots were developed by determining nucleotide sequences of bulked PCR products using a next-generation sequencer. DNA fragments were amplified from two radish lines by multiplex PCR with six primer pairs, and those amplified by 2,880 primer pairs were mixed and sequenced. By assembling sequence data, 1,953 SNPs in 750 DNA fragments, 437 of which have been previously mapped in a linkage map, were identified. A linkage map of nine linkage groups was constructed with 188 markers, and five QTLs were detected in two F(2) populations, three of them accounting for more than 50% of the total phenotypic variance being repeatedly detected. In the identified QTL regions, nine SNP markers were newly produced. By synteny analysis of the QTLs regions with Arabidopsis thaliana and Brassica rapa genome sequences, three candidate genes were selected, i.e., RsMAM3 for production of aliphatic glucosinolates linked to GSL-QTL-4, RsIPMDH1 for leucine biosynthesis showing strong co-expression with glucosinolate biosynthesis genes linked to GSL-QTL-2, and RsBCAT4 for branched-chain amino acid aminotransferase linked to GSL-QTL-1. Nucleotide sequences and expression of these genes suggested their possible function in 4MTB-GSL biosynthesis in radish roots.

Development of Capsicum EST-SSR markers for species identification and in silico mapping onto the tomato genome sequence

Capsicum spp. are widely cultivated for use as vegetables and spices. The Kihara Institute for Biological Research, Yokohama City University, Japan, has stocks of approximately 800 lines of Capsicum spp. collected from various regions of Central and South America, the regions of origin for Capsicum spp. In this study, 5,751 primer pairs for simple sequence repeat markers, based on 118,060 publicly available sequences of expressed sequence tags of Capsicum annuum, were designed and subjected to a similarity search against the genomic sequence of tomato, a model Solanaceae species. Nucleotide sequences spanning 2,245 C. annuum markers were successfully mapped onto the tomato genome, and 96 of these, which spanned the entire tomato genome, were selected for further analysis. In genotyping analysis, 60 out of the 77 markers that produced specific DNA amplicons showed polymorphism among the Capsicum lines examined. On the basis of the resulting data, the 192 tested lines were grouped into five main clusters. The additional sequencing analysis of the plastid genes, matK and rbcL, divided the resources into three groups. As a result, 19 marker loci exhibited genotypes specific to species and cluster, suggesting that the DNA markers are useful for species identification. Information on the DNA markers will contribute to Capsicum genetics, genomics, and breeding. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-012-9774-z) contains supplementary material, which is available to authorized users.

Transcriptome analysis of the roots at early and late seedling stages using Illumina paired-end sequencing and development of EST-SSR markers in radish

The tuberous root of radish is an important vegetable, but insufficient transcriptomic and genomic data are currently available to understand the molecular mechanisms underlying tuberous root formation and development. High-throughput transcriptomic sequencing is essential to generate a large transcript sequence data set for gene discovery and molecular marker development. In this study, a total of 107.3 million clean reads were generated using Illumina paired-end sequencing technology. De novo assembly generated 61,554 unigenes with an average length of 820 bp. Based on a sequence similarity search with known proteins or nucleotides, 85.51 % (52,634), 90.18 % (55,507) and 54 % (33,242) consensus sequences showed homology with sequences in the Nr, Nt and Swiss-Prot databases, respectively. Of these annotated unigenes, 21,109 and 17,343 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. A total of 27,809 unigenes were assigned to 123 pathways in the Kyoto Encyclopedia of Genes and Genomes database. Analysis of transcript differences between libraries from the early and late seedling developmental stages demonstrated that starch and sucrose metabolism and phenylpropanoid biosynthesis may be the dominant metabolic events during tuberous root formation and plant hormones probably play critical roles in regulation of this developmental process. In total, 14,641 potential EST-SSRs were identified among the unigenes, and 12,733 primer pairs for 2,511 SSR were obtained. Summarily, this study gave us a clue to understand the radish tuberous root formation and development, and also provided us with a valuable sequence resource for novel gene discovery and marker-assisted selective breeding in radish.

KEY MESSAGE

De novo assembled and characterized the radish tuberous root transcriptome; explored the mechanism of radish tuberous root formation; development of EST-SSR markers in radish.

Construction of a high-resolution linkage map of Rfd1, a restorer-of-fertility locus for cytoplasmic male sterility conferred by DCGMS cytoplasm in radish (Raphanus sativus L.) using synteny between radish and Arabidopsis genomes

Cytoplasmic male sterility caused by Dongbu cytoplasmic and genic male-sterility (DCGMS) cytoplasm and its nuclear restorer-of-fertility locus (Rfd1) with a linked molecular marker (A137) have been reported in radish (Raphanus sativus L.). To construct a linkage map of the Rfd1 locus, linked amplified fragment length polymorphism (AFLP) markers were screened using bulked segregant analysis. A 220-bp linked AFLP fragment sequence from radish showed homology with an Arabidopsis coding sequence. Using this Arabidopsis gene sequence, a simple PCR marker (A220) was developed. The A137 and A220 markers flanked the Rfd1 locus. Two homologous Arabidopsis genes with both marker sequences were positioned on Arabidopsis chromosome-3 with an interval of 2.4 Mb. To integrate the Rfd1 locus into a previously reported expressed sequence tag (EST)-simple sequence repeat (SSR) linkage map, the radish EST sequences located in three syntenic blocks within the 2.4-Mb interval were used to develop single nucleotide polymorphism (SNP) markers for tagging each block. The SNP marker in linkage group-2 co-segregated with male fertility in an F(2) population. Using radish ESTs positioned in linkage group-2, five intron length polymorphism (ILP) markers and one cleaved amplified polymorphic sequence (CAPS) marker were developed and used to construct a linkage map of the Rfd1 locus. Two closely linked markers delimited the Rfd1 locus within a 985-kb interval of Arabidopsis chromosome-3. Synteny between the radish and Arabidopsis genomes in the 985-kb interval were used to develop three ILP and three CAPS markers. Two ILP markers further delimited the Rfd1 locus to a 220-kb interval of Arabidopsis chromosome-3.

Genetic linkage map construction and QTL mapping of cadmium accumulation in radish (Raphanus sativus L.)

Cadmium (Cd) is a widespread soil pollutant and poses a significant threat to human health via the food chain. Large phenotypic variations in Cd concentration of radish roots and shoots have been observed. However, the genetic and molecular mechanisms of Cd accumulation in radish remain to be elucidated. In this study, a genetic linkage map was constructed using an F(2) mapping population derived from a cross between a high Cd-accumulating cultivar NAU-Dysx and a low Cd-accumulating cultivar NAU-Yh. The linkage map consisted of 523 SRAP, RAPD, SSR, ISSR, RAMP, and RGA markers and had a total length of 1,678.2 cM with a mean distance of 3.4 cM between two markers. All mapped markers distributed on nine linkage groups (LGs) having sizes between 134.7 and 236.8 cM. Four quantitative trait loci (QTLs) for root Cd accumulation were mapped on LGs 1, 4, 6, and 9, which accounted for 9.86 to 48.64 % of all phenotypic variance. Two QTLs associated with shoot Cd accumulation were detected on LG1 and 3, which accounted for 17.08 and 29.53 % of phenotypic variance, respectively. A major-effect QTL, qRCd9 (QTL for root Cd accumulation on LG9), was identified on LG 9 flanked by NAUrp011_754 and EM5me6_286 markers with a high LOD value of 23.6, which accounted for 48.64 % of the total phenotypic variance in Cd accumulation of F(2) lines. The results indicated that qRCd9 is a novel QTL responsible for controlling root Cd accumulation in radish, and the identification of specific molecular markers tightly linked to the major QTL could be further applied for marker-assisted selection (MAS) in low-Cd content radish breeding program.

EST-derived genic molecular markers: development and utilization for generating an advanced transcript map of chickpea

Well-saturated linkage maps especially those based on expressed sequence tag (EST)-derived genic molecular markers (GMMs) are a pre-requisite for molecular breeding. This is especially true in important legumes such as chickpea where few simple sequence repeats (SSR) and even fewer GMM-based maps have been developed. Therefore, in this study, 2,496 ESTs were generated from chickpea seeds and utilized for the development of 487 novel EST-derived functional markers which included 125 EST-SSRs, 151 intron targeted primers (ITPs), 109 expressed sequence tag polymorphisms (ESTPs), and 102 single nucleotide polymorphisms (SNPs). Whereas ESTSSRs, ITPs, and ESTPs were developed by in silico analysis of the developed EST sequences, SNPs were identified by allele resequencing and their genotyping was performedusing the Illumina GoldenGate Assay. Parental polymorphism was analyzed between C. arietinum ICC4958 and C. reticulatum PI489777, parents of the reference chickpea mapping population, using a total of 872 markers: 487 new gene-based markers developed in this study along with 385 previously published markers, of which 318 (36.5%) were found to be polymorphic and were used for genotyping. The genotypic data were integrated with the previously published data of 108 markers and an advanced linkage map was generated that contained 406 loci distributed on eight linkage groups that spanned 1,497.7 cM. The average marker density was 3.68 cM and the average number of markers per LG was 50.8. Among the mapped markers, 303 new genomic locations were defined that included 177 gene-based and 126 gSSRs (genomic SSRs) thereby producing the most advanced gene-rich map of chickpea solely based on co-dominant markers.

Size and location of radish chromosome regions carrying the fertility restorer Rfk1 gene in spring turnip rape

In spring turnip rape (Brassica rapa L. spp. oleifera), the most promising F1 hybrid system would be the Ogu-INRA CMS/Rf system. A Kosena fertility restorer gene Rfk1, homolog of the Ogura restorer gene Rfo, was successfully transferred from oilseed rape into turnip rape and that restored the fertility in female lines carrying Ogura cms. The trait was, however, unstable in subsequent generations. The physical localization of the radish chromosomal region carrying the Rfk1 gene was investigated using genomic in situ hybridization (GISH) and bacterial artificial chromosome-fluorescence in situ hybridization (BAC-FISH) methods. The metaphase chromosomes were hybridized using radish DNA as the genomic probe and BAC64 probe, which is linked with Rfo gene. Both probes showed a signal in the chromosome spreads of the restorer line 4021-2 Rfk of turnip rape but not in the negative control line 4021B. The GISH analyses clearly showed that the turnip rape restorer plants were either monosomic (2n=2x=20+1R) or disomic (2n=2x=20+2R) addition lines with one or two copies of a single alien chromosome region originating from radish. In the BAC-FISH analysis, double dot signals were detected in subterminal parts of the radish chromosome arms showing that the fertility restorer gene Rfk1 was located in this additional radish chromosome. Detected disomic addition lines were found to be unstable for turnip rape hybrid production. Using the BAC-FISH analysis, weak signals were sometimes visible in two chromosomes of turnip rape and a homologous region of Rfk1 in chromosome 9 of the B. rapa A genome was verified with BLAST analysis. In the future, this homologous area in A genome could be substituted with radish chromosome area carrying the Rfk1 gene.