A sequence-based genetic linkage map as a reference for Brassica rapa pseudochromosome assembly

Genome survey and high-resolution genetic map provide valuable genetic resources for Fenneropenaeus chinensis

Fenneropenaeus chinensis is one of the most important aquaculture species in China. Research on its genomic and genetic structure not only helps us comprehend the genetic basis of complex economic traits, but also offers theoretical guidance in selective breeding. In the present study, a genome survey sequencing was performed to generate a rough reference genome utilized for groping preliminary genome characteristics and facilitate linkage and quantitative trait locus (QTL) mapping. Linkage mapping was conducted using a reduced-representation sequencing method 2b-RAD. In total, 36,762 SNPs were genotyped from 273 progenies in a mapping family, and a high-resolution linkage map was constructed. The consensus map contained 12,884 markers and spanned 5257.81 cM with an average marker interval of 0.41 cM, which was the first high-resolution genetic map in F. chinensis to our knowledge. QTL mapping and association analysis were carried out in 29 characters including body size, sex and disease resistance. 87 significant QTLs were detected in several traits and they were also evaluated by association analysis. Results of this study provide us valuable suggestions in genetic improvement and breeding of new varieties and also lay a basic foundation for further application of cloning of economic genes in selective breeding program and marker-assisted selection.

A high-density genetic map construction and sex-related loci identification in Chinese Giant salamander

BACKGROUND

The Chinese giant salamander Andrias davidianus is an important amphibian species in China because of its increasing economic value, protection status and special evolutionary position from aquatic to terrestrial animal. Its large genome presents challenges to genetic research. Genetic linkage mapping is an important tool for genome assembly and determination of phenotype-related loci.

RESULTS

In this study, we constructed a high-density genetic linkage map using ddRAD sequencing technology to obtain SNP genotyping data of members from an full-sib family which sex had been determined. A total of 10,896 markers were grouped and oriented into 30 linkage groups, representing 30 chromosomes of A. davidianus. The genetic length of LGs ranged from 17.61 cM (LG30) to 280.81 cM (LG1), with a mean inter-locus distance ranging from 0.11(LG3) to 0.48 cM (LG26). The total genetic map length was 2643.10 cM with an average inter-locus distance of 0.24 cM. Three sex-related loci and four sex-related markers were found on LG6 and LG23, respectively.

CONCLUSION

We constructed the first High-density genetic linkage map and identified three sex-related loci in the Chinese giant salamander. Current results are expected to be a useful tool for future genomic studies aiming at the marker-assisted breeding of the species.

Construction of a High-Density Genetic Map and Identification of Quantitative Trait Loci for Nitrite Tolerance in the Pacific White Shrimp (Litopenaeus vannamei)

Nitrite is a major environmental toxin in aquaculture systems that disrupts multiple physiological functions in aquatic animals. Although nitrite tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available. In this study, we constructed a high-density genetic map of Litopenaeus vannamei with 17,242 single nucleotide polymorphism markers spanning 6,828.06 centimorgans (cM), with an average distance of 0.4 cM between adjacent markers on 44 linkage groups (LGs). Using this genetic map, we identified two markers associated with nitrite tolerance. We then sequenced the transcriptomes of the most nitrite-tolerant and nitrite-sensitive individuals from each of four genetically distinct L. vannamei families (LV-I–4). We found 2,002, 1,983, 1,954, and 1,867 differentially expressed genes in families LV-1, LV-2, LV-3, and LV-4, respectively. By integrating QTL and transcriptomics analyses, we identified a candidate gene associated with nitrite tolerance. This gene was annotated as solute carrier family 26 member 6 (SLC26A6). RNA interference (RNAi) analysis demonstrated that SLC26A6 was critical for nitrite tolerance in L. vannamei. The present study increases our understanding of the molecular mechanisms underlying nitrite tolerance in shrimp and provides a basis for molecular-marker-assisted shrimp breeding.

Genetic Analysis and Fine Mapping of a Spontaneously Mutated Male Sterility Gene in Brassica rapa ssp. chinensis

Male sterility has been widely used in hybrid seed production in Brassica, but not in B. rapa ssp. chinensis, and genetic models of male sterility for this subspecies are unclear. We discovered a spontaneous mutant in B. rapa ssp. chinensis A series of progeny tests indicated that male sterility in B. rapa ssp. chinensis follows a three-allele model with BrMsa , BrMsb , and BrMsc The male sterility locus has been mapped to chromosome A07 in BC1 and F2 populations through genotyping by sequencing. Fine mapping in a total of 1,590 F2 plants narrowed the male sterility gene BrMs to a 400 kb region, with two SNP markers only 0.3 cM from the gene. Comparative gene mapping shows that the Ms gene in B. rapa ssp. pekinensis is different from the BrMs gene of B. rapa ssp. chinensis, despite that both genes are located on chromosome A07. Interestingly, the DNA sequence orthologous to a male sterile gene in Brassica napus, BnRf, is within 400 kb of the BrMs locus. The BnRf orthologs of B. rapa ssp. chinensis were sequenced, and one KASP marker (BrMs_indel) was developed for genotyping based on a 14 bp indel at intron 4. Cosegregation of male sterility and BrMs_indel genotypes in the F2 population indicated that BnRf from B. napus and BrMs from B. rapa are likely to be orthologs. The BrMs_indel marker developed in this study will be useful in marker-assisted selection for the male sterility trait.

Genome-wide characterization of ALDH Superfamily in Brassica rapa and enhancement of stress tolerance in heterologous hosts by BrALDH7B2 expression

Aldehyde dehydrogenase (ALDH) carries out oxidation of toxic aldehydes using NAD⁺/NADP⁺ as cofactors. In the present study, we performed a genome-wide identification and expression analysis of genes in the ALDH gene family in Brassica rapa. A total of 23 ALDH genes in the superfamily have been identified according to the classification of ALDH Gene Nomenclature Committee (AGNC). They were distributed unevenly across all 10 chromosomes. All the 23 Brassica rapa ALDH (BrALDH) genes exhibited varied expression patterns during treatments with abiotic stress inducers and hormonal treatments. The relative expression profiles of ALDH genes in B. rapa showed that they are predominantly expressed in leaves and stem suggesting their function in the vegetative tissues. BrALDH7B2 showed a strong response to abiotic stress and hormonal treatments as compared to other ALDH genes; therefore, it was overexpressed in heterologous hosts, E. coli and yeast to study its possible function under abiotic stress conditions. Over-expression of BrALDH7B2 in heterologous systems, E. coli and yeast cells conferred significant tolerance to abiotic stress treatments. Results from this work demonstrate that BrALDH genes are a promising and untapped genetic resource for crop improvement and could be deployed further in the development of drought and salinity tolerance in B. rapa and other economically important crops.

A Genomic Variation Map Provides Insights into the Genetic Basis of Spring Chinese Cabbage (Brassica rapa ssp. pekinensis) Selection

Chinese cabbage is the most consumed leafy crop in East Asian countries. However, premature bolting induced by continuous low temperatures severely decreases the yield and quality of the Chinese cabbage, and therefore restricts its planting season and geographic distribution. In the past 40 years, spring Chinese cabbage with strong winterness has been selected to meet the market demand. Here, we report a genome variation map of Chinese cabbage generated from the resequencing data of 194 geographically diverse accessions of three ecotypes. In-depth analyses of the selection sweeps and genome-wide patterns revealed that spring Chinese cabbage was selected from a specific population of autumn Chinese cabbage around the area of Shandong peninsula in northern China. We identified 23 genomic loci that underwent intensive selection, and further demonstrated by gene expression and haplotype analyses that the incorporation of elite alleles of VERNALISATION INSENTIVE 3.1 (BrVIN3.1) and FLOWER LOCUS C 1 (BrFLC1) is a determinant genetic source of variation during selection. Moreover, we showed that the quantitative response of BrVIN3.1 to cold due to the sequence variations in the cis elements of the BrVIN3.1 promoter significantly contributes to bolting-time variation in Chinese cabbage. Collectively, our study provides valuable insights into the genetic basis of spring Chinese cabbage selection and will facilitate the breeding of bolting-resistant varieties by molecular-marker-assisted selection, transgenic or gene editing approaches.

High-density genetic map construction and QTLs identification for plant height in white jute (Corchorus capsularis L.) using specific locus amplified fragment (SLAF) sequencing

Background

Genetic mapping and quantitative trait locus (QTL) detection are powerful methodologies in plant improvement and breeding. White jute (Corchorus capsularis L.) is an important industrial raw material fiber crop because of its elite characteristics. However, construction of a high-density genetic map and identification of QTLs has been limited in white jute due to a lack of sufficient molecular markers. The specific locus amplified fragment sequencing (SLAF-seq) strategy combines locus-specific amplification and high-throughput sequencing to carry out de novo single nuclear polymorphism (SNP) discovery and large-scale genotyping. In this study, SLAF-seq was employed to obtain sufficient markers to construct a high-density genetic map for white jute. Moreover, with the development of abundant markers, genetic dissection of fiber yield traits such as plant height was also possible. Here, we present QTLs associated with plant height that were identified using our newly constructed genetic linkage groups.

Results

An F₈ population consisting of 100 lines was developed. In total, 69,446 high-quality SLAFs were detected of which 5,074 SLAFs were polymorphic; 913 polymorphic markers were used for the construction of a genetic map. The average coverage for each SLAF marker was 43-fold in the parents, and 9.8-fold in each F₈ individual. A linkage map was constructed that contained 913 SLAFs on 11 linkage groups (LGs) covering 1621.4 cM with an average density of 1.61 cM per locus. Among the 11 LGs, LG1 was the largest with 210 markers, a length of 406.34 cM, and an average distance of 1.93 cM between adjacent markers. LG11 was the smallest with only 25 markers, a length of 29.66 cM, and an average distance of 1.19 cM between adjacent markers. ‘SNP_only’ markers accounted for 85.54% and were the predominant markers on the map. QTL mapping based on the F₈ phenotypes detected 11 plant height QTLs including one major effect QTL across two cultivation locations, with each QTL accounting for 4.14–15.63% of the phenotypic variance.

Conclusions

To our knowledge, the linkage map constructed here is the densest one available to date for white jute. This analysis also identified the first QTL in white jute. The results will provide an important platform for gene/QTL mapping, sequence assembly, genome comparisons, and marker-assisted selection breeding for white jute.

A genome-wide SNP-based genetic map and QTL mapping for agronomic traits in Chinese cabbage

The aim of this work was to construct a high-resolution genetic map for the dissection of complex morphological and agronomic traits in Chinese cabbage (Brassica rapa L. syn. B. campestris). Chinese cabbage, an economically important vegetable, is a good model plant for studies on the evolution of morphologic variation. Herein, two high-generation inbred Chinese cabbage lines, 'Huangxiaoza' and 'Bqq094-11', were crossed. Then restriction-site-associated DNA sequencing (RAD-seq) was performed on the parents and 120 F₂ individuals. A genetic map containing 711 bins representing 3985 single nucleotide polymorphism (SNP) markers was constructed. By using WinQTL with composite interval mapping (CIM) and mixed-model based composite interval mapping (MCIM) analysis via QTLNetwork, quantitative trait loci (QTL) linked to 16 genetic traits related to plant size, color and leaf characteristics were mapped to 10 linkage groups. The high density genetic map and QTL identified for morphological and agronomic traits lay the groundwork for functional gene mapping, map-based cloning and marker-assisted selection (MAS) in Chinese cabbage.

Fine mapping and candidate gene analysis of an anthocyanin-rich gene, BnaA.PL1, conferring purple leaves in Brassica napus L

Because of the advantages of anthocyanins, the genetics and breeding of crops rich in anthocyanins has become a hot research topic. However, due to the lack of anthocyanin-related mutants, no regulatory genes have been mapped in Brassica napus. In this study, we first report the characterization of a B. napus line with purple leaves and the fine mapping and candidate screening of the BnaA.PL1 gene. The amount of anthocyanins in the purple leaf line was six times higher than that in a green leaf line. A genetic analysis indicated that the purple character was controlled by an incomplete dominant gene. Through map-based cloning, we localized the BnaA.PL1 gene to a 99-kb region at the end of B. napus chromosome A03. Transcriptional analysis of 11 genes located in the target region revealed that the expression level of only the BnAPR2 gene in seedling leaves decreased from purple to reddish green to green individuals, a finding that was consistent with the measured anthocyanin accumulation levels. Molecular cloning and sequence analysis of BnAPR2 showed that the purple individual-derived allele contained 17 variants. Markers co-segregating with BnaA.PL1 were developed from the sequence of BnAPR2 and were validated in the BC4P2 population. These results suggested that BnAPR2, which encodes adenosine 5'-phosphosulfate reductase, is likely to be a valuable candidate gene. This work may lay the foundation for the marker-assisted selection of B. napus vegetables that are rich in anthocyanins and for an improved understanding of the molecular mechanisms controlling anthocyanin accumulation in Brassica.

Development of cleaved amplified polymorphic sequence markers and a CAPS-based genetic linkage map in watermelon (Citrullus lanatus [Thunb.] Matsum. and Nakai) constructed using whole-genome re-sequencing data

Cleaved amplified polymorphic sequence (CAPS) markers are useful tools for detecting single nucleotide polymorphisms (SNPs). This study detected and converted SNP sites into CAPS markers based on high-throughput re-sequencing data in watermelon, for linkage map construction and quantitative trait locus (QTL) analysis. Two inbred lines, Cream of Saskatchewan (COS) and LSW-177 had been re-sequenced and analyzed by Perl self-compiled script for CAPS marker development. 88.7% and 78.5% of the assembled sequences of the two parental materials could map to the reference watermelon genome, respectively. Comparative assembled genome data analysis provided 225,693 and 19,268 SNPs and indels between the two materials. 532 pairs of CAPS markers were designed with 16 restriction enzymes, among which 271 pairs of primers gave distinct bands of the expected length and polymorphic bands, via PCR and enzyme digestion, with a polymorphic rate of 50.94%. Using the new CAPS markers, an initial CAPS-based genetic linkage map was constructed with the F2 population, spanning 1836.51 cM with 11 linkage groups and 301 markers. 12 QTLs were detected related to fruit flesh color, length, width, shape index, and brix content. These newly CAPS markers will be a valuable resource for breeding programs and genetic studies of watermelon.

User Guidelines for the Brassica Database: BRAD

The genome sequence of Brassica rapa was first released in 2011. Since then, further Brassica genomes have been sequenced or are undergoing sequencing. It is therefore necessary to develop tools that help users to mine information from genomic data efficiently. This will greatly aid scientific exploration and breeding application, especially for those with low levels of bioinformatic training. Therefore, the Brassica database (BRAD) was built to collect, integrate, illustrate, and visualize Brassica genomic datasets. BRAD provides useful searching and data mining tools, and facilitates the search of gene annotation datasets, syntenic or non-syntenic orthologs, and flanking regions of functional genomic elements. It also includes genome-analysis tools such as BLAST and GBrowse. One of the important aims of BRAD is to build a bridge between Brassica crop genomes with the genome of the model species Arabidopsis thaliana, thus transferring the bulk of A. thaliana gene study information for use with newly sequenced Brassica crops.

A High-Density Genetic Linkage Map for Cucumber (Cucumis sativus L.): Based on Specific Length Amplified Fragment (SLAF) Sequencing and QTL Analysis of Fruit Traits in Cucumber

High-density genetic linkage map plays an important role in genome assembly and quantitative trait loci (QTL) fine mapping. Since the coming of next-generation sequencing, makes the structure of high-density linkage maps much more convenient and practical, which simplifies SNP discovery and high-throughput genotyping. In this research, a high-density linkage map of cucumber was structured using specific length amplified fragment sequencing, using 153 F2 populations of S1000 × S1002. The high-density genetic map composed 3,057 SLAFs, including 4,475 SNP markers on seven chromosomes, and spanned 1061.19 cM. The average genetic distance is 0.35 cM. Based on this high-density genome map, QTL analysis was performed on two cucumber fruit traits, fruit length and fruit diameter. There are 15 QTLs for the two fruit traits were detected.

A transposon insertion in FLOWERING LOCUS T is associated with delayed flowering in Brassica rapa

Long days and vernalization accelerate the transition from vegetative growth to reproductive growth in Brassica rapa. Bolting before plants reach the harvesting stage is a serious problem in B. rapa vegetable crop cultivation. The genetic dissection of flowering time is important for breeding of premature bolting-resistant B. rapa crops. Using a recombinant inbred line (RIL) population, we twice detected two major quantitative trait loci (QTLs) for flowering time in two different growing seasons that were located on chromosomes A02 and A07, respectively. We hypothesized that an orthologue of the Arabidopsis thaliana FLOWERING LOCUS T (FT) gene, named as BrFT2, was the candidate gene underlying the QTL localized to A07. A transposon insertion in the second intron of BrFT2 was detected in one of the parental lines, which was predicted to generate a loss-of-function allele. Transcription analysis revealed that the BrFT2 transcript was not present in the parental line that harbored the mutated allele. RILs carrying only the mutated BrFT2 allele showed delayed flowering regardless of growing seasons when compared to RILs carrying the wild-type BrFT2 allele. These data suggest that BrFT2 is involved in flowering time regulation in controlling flowering time in B. rapa.

Genome survey and high-density genetic map construction provide genomic and genetic resources for the Pacific White Shrimp Litopenaeus vannamei

The Pacific white shrimp Litopenaeus vannamei is the dominant crustacean species in global seafood mariculture. Understanding the genome and genetic architecture is useful for deciphering complex traits and accelerating the breeding program in shrimp. In this study, a genome survey was conducted and a high-density linkage map was constructed using a next-generation sequencing approach. The genome survey was used to identify preliminary genome characteristics and to generate a rough reference for linkage map construction. De novo SNP discovery resulted in 25,140 polymorphic markers. A total of 6,359 high-quality markers were selected for linkage map construction based on marker coverage among individuals and read depths. For the linkage map, a total of 6,146 markers spanning 4,271.43 cM were mapped to 44 sex-averaged linkage groups, with an average marker distance of 0.7 cM. An integration analysis linked 5,885 genome scaffolds and 1,504 BAC clones to the linkage map. Based on the high-density linkage map, several QTLs for body weight and body length were detected. This high-density genetic linkage map reveals basic genomic architecture and will be useful for comparative genomics research, genome assembly and genetic improvement of L. vannamei and other penaeid shrimp species.

QTL mapping of leafy heads by genome resequencing in the RIL population of Brassica rapa

Leaf heads of cabbage (Brassica oleracea), Chinese cabbage (B. rapa), and lettuce (Lactuca sativa) are important vegetables that supply mineral nutrients, crude fiber and vitamins in the human diet. Head size, head shape, head weight, and heading time contribute to yield and quality. In an attempt to investigate genetic basis of leafy head in Chinese cabbage (B. rapa), we took advantage of recent technical advances of genome resequencing to perform quantitative trait locus (QTL) mapping using 150 recombinant inbred lines (RILs) derived from the cross between heading and non-heading Chinese cabbage. The resequenced genomes of the parents uncovered more than 1 million SNPs. Genotyping of RILs using the high-quality SNPs assisted by Hidden Markov Model (HMM) generated a recombination map. The raw genetic map revealed some physical assembly error and missing fragments in the reference genome that reduced the quality of SNP genotyping. By deletion of the genetic markers in which recombination rates higher than 20%, we have obtained a high-quality genetic map with 2209 markers and detected 18 QTLs for 6 head traits, from which 3 candidate genes were selected. These QTLs provide the foundation for study of genetic basis of leafy heads and the other complex traits.

Genome-wide DNA polymorphisms in two cultivars of mei (Prunus mume sieb. et zucc.)

Background

Mei (Prunus mume Sieb. et Zucc.) is a famous ornamental plant and fruit crop grown in East Asian countries. Limited genetic resources, especially molecular markers, have hindered the progress of mei breeding projects. Here, we performed low-depth whole-genome sequencing of Prunus mume ‘Fenban’ and Prunus mume ‘Kouzi Yudie’ to identify high-quality polymorphic markers between the two cultivars on a large scale.

Results

A total of 1464.1 Mb and 1422.1 Mb of ‘Fenban’ and ‘Kouzi Yudie’ sequencing data were uniquely mapped to the mei reference genome with about 6-fold coverage, respectively. We detected a large number of putative polymorphic markers from the 196.9 Mb of sequencing data shared by the two cultivars, which together contained 200,627 SNPs, 4,900 InDels, and 7,063 SSRs. Among these markers, 38,773 SNPs, 174 InDels, and 418 SSRs were distributed in the 22.4 Mb CDS region, and 63.0% of these marker-containing CDS sequences were assigned to GO terms. Subsequently, 670 selected SNPs were validated using an Agilent’s SureSelect solution phase hybridization assay. A subset of 599 SNPs was used to assess the genetic similarity of a panel of mei germplasm samples and a plum (P. salicina) cultivar, producing a set of informative diversity data. We also analyzed the frequency and distribution of detected InDels and SSRs in mei genome and validated their usefulness as DNA markers. These markers were successfully amplified in the cultivars and in their segregating progeny.

Conclusions

A large set of high-quality polymorphic SNPs, InDels, and SSRs were identified in parallel between ‘Fenban’ and ‘Kouzi Yudie’ using low-depth whole-genome sequencing. The study presents extensive data on these polymorphic markers, which can be useful for constructing high-resolution genetic maps, performing genome-wide association studies, and designing genomic selection strategies in mei.

Shifts in the evolutionary rate and intensity of purifying selection between two Brassica genomes revealed by analyses of orthologous transposons and relics of a whole genome triplication

Recent sequencing of the Brassica rapa and Brassica oleracea genomes revealed extremely contrasting genomic features such as the abundance and distribution of transposable elements between the two genomes. However, whether and how these structural differentiations may have influenced the evolutionary rates of the two genomes since their split from a common ancestor are unknown. Here, we investigated and compared the rates of nucleotide substitution between two long terminal repeats (LTRs) of individual orthologous LTR-retrotransposons, the rates of synonymous and non-synonymous substitution among triplicated genes retained in both genomes from a shared whole genome triplication event, and the rates of genetic recombination estimated/deduced by the comparison of physical and genetic distances along chromosomes and ratios of solo LTRs to intact elements. Overall, LTR sequences and genic sequences showed more rapid nucleotide substitution in B. rapa than in B. oleracea. Synonymous substitution of triplicated genes retained from a shared whole genome triplication was detected at higher rates in B. rapa than in B. oleracea. Interestingly, non-synonymous substitution was observed at lower rates in the former than in the latter, indicating shifted densities of purifying selection between the two genomes. In addition to evolutionary asymmetry, orthologous genes differentially regulated and/or disrupted by transposable elements between the two genomes were also characterized. Our analyses suggest that local genomic and epigenomic features, such as recombination rates and chromatin dynamics reshaped by independent proliferation of transposable elements and elimination between the two genomes, are perhaps partially the causes and partially the outcomes of the observed inter-specific asymmetric evolution.

Construction of a high-density genetic map for sesame based on large scale marker development by specific length amplified fragment (SLAF) sequencing

BACKGROUND

The genetics and molecular biology of sesame has only recently begun to be studied even though sesame is an important oil seed crop. A high-density genetic map for sesame has not been published yet due to a lack of sufficient molecular markers. Specific length amplified fragment sequencing (SLAF-seq) is a recently developed high-resolution strategy for large-scale de novo SNP discovery and genotyping. SLAF-seq was employed in this study to obtain sufficient markers to construct a high-density genetic map for sesame.

RESULTS

In total, 28.21 Gb of data containing 201,488,285 pair-end reads was obtained after sequencing. The average coverage for each SLAF marker was 23.48-fold in the male parent, 23.38-fold in the female parent, and 14.46-fold average in each F2 individual. In total, 71,793 high-quality SLAFs were detected of which 3,673 SLAFs were polymorphic and 1,272 of the polymorphic markers met the requirements for use in the construction of a genetic map. The final map included 1,233 markers on the 15 linkage groups (LGs) and was 1,474.87 cM in length with an average distance of 1.20 cM between adjacent markers. To our knowledge, this map is the densest genetic linkage map to date for sesame. 'SNP_only' markers accounted for 87.51% of the markers on the map. A total of 205 markers on the map showed significant (P < 0.05) segregation distortion.

CONCLUSIONS

We report here the first high-density genetic map for sesame. The map was constructed using an F2 population and the SLAF-seq approach, which allowed the efficient development of a large number of polymorphic markers in a short time. Results of this study will not only provide a platform for gene/QTL fine mapping, map-based gene isolation, and molecular breeding for sesame, but will also serve as a reference for positioning sequence scaffolds on a physical map, to assist in the process of assembling the sesame genome sequence.

RNA-seq based SNPs in some agronomically important oleiferous lines of Brassica rapa and their use for genome-wide linkage mapping and specific-region fine mapping

Background

Brassica rapa (AA) contains very diverse forms which include oleiferous types and many vegetable types. Genome sequence of B. rapa line Chiifu (ssp. pekinensis), a leafy vegetable type, was published in 2011. Using this knowledge, it is important to develop genomic resources for the oleiferous types of B. rapa. This will allow more involved molecular mapping, in-depth study of molecular mechanisms underlying important agronomic traits and introgression of traits from B. rapa to major oilseed crops - B. juncea (AABB) and B. napus (AACC). The study explores the availability of SNPs in RNA-seq generated contigs of three oleiferous lines of B. rapa - Candle (ssp. oleifera, turnip rape), YSPB-24 and Tetra (ssp. trilocularis, Yellow sarson) and their use in genome-wide linkage mapping and specific-region fine mapping using a RIL population between Chiifu and Tetra.

Results

RNA-seq was carried out on the RNA isolated from young inflorescences containing unopened floral buds, floral axis and small leaves, using Illumina paired-end sequencing technology. Sequence assembly was carried out using the Velvet de-novo programme and the assembled contigs were organised against Chiifu gene models, available in the BRAD-CDS database. RNA-seq confirmed the presence of more than 17,000 single-copy gene models described in the BRAD database. The assembled contigs and the BRAD gene models were analyzed for the presence of SSRs and SNPs. While the number of SSRs was limited, more than 0.2 million SNPs were observed between Chiifu and the three oleiferous lines. Assays for SNPs were designed using KASPar technology and tested on a F₇-RIL population derived from a Chiifu x Tetra cross. The design of the SNP assays were based on three considerations - the 50 bp flanking region of the SNPs should be strictly similar, the SNP should have a read-depth of ≥7 and no exon/intron junction should be present within the 101 bp target region. Using these criteria, a total of 640 markers (580 for genome-wide mapping and 60 for specific-region mapping) marking as many genes were tested for mapping. Out of 640 markers that were tested, 594 markers could be mapped unambiguously which included 542 markers for genome-wide mapping and 42 markers for fine mapping of the tet-o locus that is involved with the trait tetralocular ovary in the line Tetra.

Conclusion

A large number of SNPs and PSVs are present in the transcriptome of B. rapa lines for genome-wide linkage mapping and specific-region fine mapping. Criteria used for SNP identification delivered markers, more than 93% of which could be successfully mapped to the F₇–RIL population of Chiifu x Tetra cross.

Development of InDel markers for Brassica rapa based on whole-genome re-sequencing

Genome-wide detection of short insertion/deletion length polymorphisms (InDels, <5 bp) in Brassica rapa (named the A genome) was performed by comparing whole-genome re-sequencing data from two B. rapa accessions, L144 and Z16, to the reference genome sequence of Chiifu-401-42. In total, we identified 108,558 InDel polymorphisms between Chiifu-401-42 and L144, 26,795 InDels between Z16 and Chiifu-401-42, and 26,693 InDels between L144 and Z16. From these, 639 InDel polymorphisms of 3-5 bp in length between L144 and Z16 were selected for experimental validation; 491 (77%) yielded single PCR fragments and showed polymorphisms, 7 (1%) did not amplify a product, and 141 (22%) showed no polymorphism. For further validation of these intra-specific InDel polymorphisms, 503 candidates, randomly selected from the 639 InDels, were screened across seven accessions representing different B. rapa cultivar groups. Of these assayed markers, 387 (77%) were polymorphic, 111 (22%) were not polymorphic and 5 (1%) did not amplify a PCR product. Furthermore, we randomly selected 518 InDel markers to validate their polymorphism in B. napus (the AC genome) and B. juncea (the AB genome), of which more than 90% amplified a PCR product; 132 (25%) showed polymorphism between the two B. napus accessions and 41 (8%) between the two B. juncea accessions. This set of novel PCR-based InDel markers will be a valuable resource for genetic studies and breeding programs in B. rapa.

Mapping and candidate-gene screening of the novel Turnip mosaic virus resistance gene retr02 in Chinese cabbage (Brassica rapa L.)

The extreme resistance to Turnip mosaic virus observed in the Chinese cabbage (Brassica rapa) line, BP8407, is monogenic and recessive. Bulked segregant analysis was carried out to identify simple sequence repeat and Indel markers linked to this recessive resistance gene, termed recessive Turnip mosaic virus resistance 02 (retr02). Mapping of PCR-specific Indel markers on 239 individuals of a BP8407 × Ji Zao Chun F(2) population, located this resistance gene to a 0.9-cM interval between two Indel markers (BrID10694 and BrID101309) and in scaffold000060 or scaffold000104 on chromosome A04 of the B. rapa genome. Eleven eukaryotic initiation factor 4E (eIF4E) and 14 eukaryotic initiation factor 4G (eIF4G) genes are predicted in the B. rapa genome. A candidate gene, Bra035393 on scaffold000104, was predicted within the mapped resistance locus. The gene encodes the eIF(iso)4E protein. Bra035393 was sequenced in BP8407 and Ji Zao Chun. A polymorphism (A/G) was found in exon 3 between BP8407 and Ji Zao Chun. This gene was analysed in four resistant and three susceptible lines. A correlation was observed between the amino acid substitution (Gly/Asp) in the eIF(iso)4E protein and resistance/susceptibility. eIF(iso)4E has been shown previously to interact with the TuMV genome-linked protein, VPg.

A reference linkage map for Eucalyptus

Background

Genetic linkage maps are invaluable resources in plant research. They provide a key tool for many genetic applications including: mapping quantitative trait loci (QTL); comparative mapping; identifying unlinked (i.e. independent) DNA markers for fingerprinting, population genetics and phylogenetics; assisting genome sequence assembly; relating physical and recombination distances along the genome and map-based cloning of genes. Eucalypts are the dominant tree species in most Australian ecosystems and of economic importance globally as plantation trees. The genome sequence of E. grandis has recently been released providing unprecedented opportunities for genetic and genomic research in the genus. A robust reference linkage map containing sequence-based molecular markers is needed to capitalise on this resource. Several high density linkage maps have recently been constructed for the main commercial forestry species in the genus (E. grandis, E. urophylla and E. globulus) using sequenced Diversity Arrays Technology (DArT) and microsatellite markers. To provide a single reference linkage map for eucalypts a composite map was produced through the integration of data from seven independent mapping experiments (1950 individuals) using a marker-merging method.

Results

The composite map totalled 1107 cM and contained 4101 markers; comprising 3880 DArT, 213 microsatellite and eight candidate genes. Eighty-one DArT markers were mapped to two or more linkage groups, resulting in the 4101 markers being mapped to 4191 map positions. Approximately 13% of DArT markers mapped to identical map positions, thus the composite map contained 3634 unique loci at an average interval of 0.31 cM.

Conclusion

The composite map represents the most saturated linkage map yet produced in Eucalyptus. As the majority of DArT markers contained on the map have been sequenced, the map provides a direct link to the E. grandis genome sequence and will serve as an important reference for progressing eucalypt research.

Biased gene fractionation and dominant gene expression among the subgenomes of Brassica rapa

Polyploidization, both ancient and recent, is frequent among plants. A “two-step theory" was proposed to explain the meso-triplication of the Brassica “A" genome: Brassica rapa. By accurately partitioning of this genome, we observed that genes in the less fractioned subgenome (LF) were dominantly expressed over the genes in more fractioned subgenomes (MFs: MF1 and MF2), while the genes in MF1 were slightly dominantly expressed over the genes in MF2. The results indicated that the dominantly expressed genes tended to be resistant against gene fractionation. By re-sequencing two B. rapa accessions: a vegetable turnip (VT117) and a Rapid Cycling line (L144), we found that genes in LF had less non-synonymous or frameshift mutations than genes in MFs; however mutation rates were not significantly different between MF1 and MF2. The differences in gene expression patterns and on-going gene death among the three subgenomes suggest that “two-step" genome triplication and differential subgenome methylation played important roles in the genome evolution of B. rapa.

Genetic analysis of morphological traits in a new, versatile, rapid-cycling Brassica rapa recombinant inbred line population

A recombinant inbred line (RIL) population was produced based on a wide cross between the rapid-cycling and self-compatible genotypes L58, a Caixin vegetable type, and R-o-18, a yellow sarson oil type. A linkage map based on 160 F7 lines was constructed using 100 Single nucleotide polymorphisms (SNPs), 130 AFLP®, 27 InDel, and 13 publicly available SSR markers. The map covers a total length of 1150 centiMorgan (cM) with an average resolution of 4.3 cM/marker. To demonstrate the versatility of this new population, 17 traits, related to plant architecture and seed characteristics, were subjected to quantitative trait loci (QTL) analysis. A total of 47 QTLs were detected, each explaining between 6 and 54% of the total phenotypic variance for the concerned trait. The genetic analysis shows that this population is a useful new tool for analyzing genetic variation for interesting traits in B. rapa, and for further exploitation of the recent availability of the B. rapa whole genome sequence for gene cloning and gene function analysis.

BRAD, the genetics and genomics database for Brassica plants

BACKGROUND

Brassica species include both vegetable and oilseed crops, which are very important to the daily life of common human beings. Meanwhile, the Brassica species represent an excellent system for studying numerous aspects of plant biology, specifically for the analysis of genome evolution following polyploidy, so it is also very important for scientific research. Now, the genome of Brassica rapa has already been assembled, it is the time to do deep mining of the genome data.

DESCRIPTION

BRAD, the Brassica database, is a web-based resource focusing on genome scale genetic and genomic data for important Brassica crops. BRAD was built based on the first whole genome sequence and on further data analysis of the Brassica A genome species, Brassica rapa (Chiifu-401-42). It provides datasets, such as the complete genome sequence of B. rapa, which was de novo assembled from Illumina GA II short reads and from BAC clone sequences, predicted genes and associated annotations, non coding RNAs, transposable elements (TE), B. rapa genes' orthologous to those in A. thaliana, as well as genetic markers and linkage maps. BRAD offers useful searching and data mining tools, including search across annotation datasets, search for syntenic or non-syntenic orthologs, and to search the flanking regions of a certain target, as well as the tools of BLAST and Gbrowse. BRAD allows users to enter almost any kind of information, such as a B. rapa or A. thaliana gene ID, physical position or genetic marker.

CONCLUSION

BRAD, a new database which focuses on the genetics and genomics of the Brassica plants has been developed, it aims at helping scientists and breeders to fully and efficiently use the information of genome data of Brassica plants. BRAD will be continuously updated and can be accessed through http://brassicadb.org.

Extensive chromosome homoeology among Brassiceae species were revealed by comparative genetic mapping with high-density EST-based SNP markers in radish (Raphanus sativus L.)

A linkage map of expressed sequence tag (EST)-based markers in radish (Raphanus sativus L.) was constructed using a low-cost and high-efficiency single-nucleotide polymorphism (SNP) genotyping method named multiplex polymerase chain reaction–mixed probe dot-blot analysis developed in this study. Seven hundred and forty-six SNP markers derived from EST sequences of R. sativus were assigned to nine linkage groups with a total length of 806.7 cM. By BLASTN, 726 markers were found to have homologous genes in Arabidopsis thaliana, and 72 syntenic regions, which have great potential for utilizing genomic information of the model species A. thaliana in basic and applied genetics of R. sativus, were identified. By construction and analysis of the genome structures of R. sativus based on the 24 genomic blocks within the Brassicaceae ancestral karyotype, 23 of the 24 genomic blocks were detected in the genome of R. sativus, and half of them were found to be triplicated. Comparison of the genome structure of R. sativus with those of the A, B, and C genomes of Brassica species and that of Sinapis alba L. revealed extensive chromosome homoeology among Brassiceae species, which would facilitate transfer of the genomic information from one Brassiceae species to another.