Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.)

Genome-wide association studies identifies genetic loci related to fatty acid and branched-chain amino acid metabolism and histone modifications under varying nitrogen treatments in safflower (Carthamus tinctorius)

Effective identification and usage of genetic variation are prerequisites for developing nutrient-efficient cultivars. A collection of 94 safflower (Carthamus tinctorius ) genotypes (G) was investigated for important morphological and photosynthetic traits at four nitrogen (N) treatments. We found significant variation for all the studied traits except chlorophyll b (chl b ) among safflower genotypes, nitrogen treatments and G×N interaction. The examined traits showed a 2.82-50.00% increase in response to N application. Biological yield (BY) reflected a significantly positive correlation with fresh shoot weight (FSW), root length (RL), fresh root weight (FRW) and number of leaves (NOL), while a significantly positive correlation was also observed among carotenoids (C), chlorophyll a (chl a ), chl b and total chlorophyll content (CT) under all treatments. Superior genotypes with respect to plant height (PH), FSW, NOL, RL, FRW and BY were clustered into Group 3, while genotypes with better mean performance regarding chl a , chl b C and CT were clustered into Group 2 as observed in principal component analysis. The identified eight best-performing genotypes could be useful to develop improved nitrogen efficient cultivars. Genome-wide association analysis resulted in 32 marker-trait associations (MTAs) under four treatments. Markers namely DArT-45481731 , DArT-17812864 , DArT-15670279 and DArT-45482737 were found consistent. Protein-protein interaction networks of loci associated with MTAs were related to fatty acid and branched-chain amino acid metabolism and histone modifications.

Adaptation Strategies to Improve the Resistance of Oilseed Crops to Heat Stress Under a Changing Climate: An Overview

Temperature is one of the decisive environmental factors that is projected to increase by 1. 5°C over the next two decades due to climate change that may affect various agronomic characteristics, such as biomass production, phenology and physiology, and yield-contributing traits in oilseed crops. Oilseed crops such as soybean, sunflower, canola, peanut, cottonseed, coconut, palm oil, sesame, safflower, olive etc., are widely grown. Specific importance is the vulnerability of oil synthesis in these crops against the rise in climatic temperature, threatening the stability of yield and quality. The natural defense system in these crops cannot withstand the harmful impacts of heat stress, thus causing a considerable loss in seed and oil yield. Therefore, a proper understanding of underlying mechanisms of genotype-environment interactions that could affect oil synthesis pathways is a prime requirement in developing stable cultivars. Heat stress tolerance is a complex quantitative trait controlled by many genes and is challenging to study and characterize. However, heat tolerance studies to date have pointed to several sophisticated mechanisms to deal with the stress of high temperatures, including hormonal signaling pathways for sensing heat stimuli and acquiring tolerance to heat stress, maintaining membrane integrity, production of heat shock proteins (HSPs), removal of reactive oxygen species (ROS), assembly of antioxidants, accumulation of compatible solutes, modified gene expression to enable changes, intelligent agricultural technologies, and several other agronomic techniques for thriving and surviving. Manipulation of multiple genes responsible for thermo-tolerance and exploring their high expressions greatly impacts their potential application using CRISPR/Cas genome editing and OMICS technology. This review highlights the latest outcomes on the response and tolerance to heat stress at the cellular, organelle, and whole plant levels describing numerous approaches applied to enhance thermos-tolerance in oilseed crops. We are attempting to critically analyze the scattered existing approaches to temperature tolerance used in oilseeds as a whole, work toward extending studies into the field, and provide researchers and related parties with useful information to streamline their breeding programs so that they can seek new avenues and develop guidelines that will greatly enhance ongoing efforts to establish heat stress tolerance in oilseeds.

Unraveling the impact on agronomic traits of the genetic architecture underlying plant-density responses in canola

Plant density defines vegetative architecture and the competition for light between individuals. Brassica napus (canola, rapeseed) presents a radically different plant architecture compared to traditional crops commonly cultivated at high density, and can act as a model system of indeterminate growth. Using a panel of 152 spring-type accessions and a double-haploid population of 99 lines from a cross between the cultivars Lynx and Monty, we performed genome-wide association studies (GWAS) and quantitative trait locus (QTL) mapping for 12 growth and yield traits at two contrasting plant densities of 15 and 60 plants m-2. The most significant associations were found for time to flowering, biomass at harvest, plant height, silique and seed numbers, and seed yield. These were generally independent of plant density, but some density-dependent associations were found in low-density populations. RNA-seq transcriptomic analysis revealed distinctive latent gene-regulatory responses to simulated shade between Lynx and Monty. Having identified candidate genes within the canola QTLs, we further examined their influence on density responses in Arabidopsis lines mutated in certain homologous genes. The results suggested that TCP1 might promote growth independently of plant density, while HY5 could increase biomass and seed yield specifically at high plant density. For flowering time, the results suggested that PIN genes might accelerate flowering in plant a density-dependent manner whilst FT, HY5, and TCP1 might accelerate it in a density-independent. This work highlights the advantages of using agronomic field experiments together with genetic and transcriptomic approaches to decipher quantitative complex traits that potentially mediate improved crop productivity.

Dissecting the Genetic Basis of Lateral and Central Spikelet Development and Grain Traits in Intermedium-Spike Barley (Hordeum vulgare Convar. Intermedium)

Barley (Hordeum vulgare L.) is one of the major grain crops worldwide and considered as a model plant for temperate cereals. One of the barley row-type groups, named intermedium-barley, was used in our previous study where we reported that other genetic loci rather than vrs1 and Int-c could play a role in lateral spikelet development and even in setting grains. To continue this work, we used phenotypic and genotypic data of 254 intermedium-spike barley accessions aimed at dissecting the genetic basis of development and grain traits of lateral and central spikelet using genome wide association (GWAS) analysis. After genotypic data filtering, 8,653 single-nucleotide polymorphism (SNPs) were used for GWAS analysis. A total of 169 significant associations were identified and we focused only on the subset of associations that exceeded the p < 10⁻⁴ threshold. Thirty-three highly significant marker-trait-associations (MTAs), represented in 28 different SNPs on all seven chromosomes for the central and/or lateral spikelet traits; such as kernel length, width, area, weight, unfilled spikelet and 1000-kernel weight, were detected. Highly significant associated markers were anchored physically using barley genome sequencing to identify candidate genes to either contain the SNPs or the closest gene to the SNP position. The results showed that 12 MTAs were specific for lateral spikelet traits, nine MTAs were specific for central spikelet traits and seven MTAs for both central and lateral traits. All together, the GWAS and candidate gene results support our hypothesis that lateral spikelet development could be regulated by loci different from those regulating central spikelet development.

Identification of Regulatory SNPs Associated with Vicine and Convicine Content of Vicia faba Based on Genotyping by Sequencing Data Using Deep Learning

Faba bean (Vicia faba) is a grain legume, which is globally grown for both human consumption as well as feed for livestock. Despite its agro-ecological importance the usage of Vicia faba is severely hampered by its anti-nutritive seed-compounds vicine and convicine (V+C). The genes responsible for a low V+C content have not yet been identified. In this study, we aim to computationally identify regulatory SNPs (rSNPs), i.e., SNPs in promoter regions of genes that are deemed to govern the V+C content of Vicia faba. For this purpose we first trained a deep learning model with the gene annotations of seven related species of the Leguminosae family. Applying our model, we predicted putative promoters in a partial genome of Vicia faba that we assembled from genotyping-by-sequencing (GBS) data. Exploiting the synteny between Medicago truncatula and Vicia faba, we identified two rSNPs which are statistically significantly associated with V+C content. In particular, the allele substitutions regarding these rSNPs result in dramatic changes of the binding sites of the transcription factors (TFs) MYB4, MYB61, and SQUA. The knowledge about TFs and their rSNPs may enhance our understanding of the regulatory programs controlling V+C content of Vicia faba and could provide new hypotheses for future breeding programs.

Evolution in crop improvement approaches and future prospects of molecular markers to CRISPR/Cas9 system

The advent of genetic selection and genome modification method assure about a real novel reformation in biotechnology and genetic engineering. With the extensive capabilities of molecular markers of them being stable, cost-effective and easy to use, they ultimately become a potent tool for variety of applications such a gene targeting, selection, editing, functional genomics; mainly for the improvisation of commercially important crops. Three main benefits of molecular marker in the field of agriculture and crop improvement programmes first, reduction of the duration of breeding programmes, second, they allow creation of new genetic variation and genetic diversity of plants and third most promising benefit is help in production of engineered plant for disease resistance, or resistance from pathogen and herbicides. This review is anticipated to present an outline how the techniques have been evolved from the simple conventional applications of DNA based molecular markers to highly throughput CRISPR technology and geared the crop yield. Techniques like using Zinc Finger Nucleases (ZFNs), Transcription Activator-Like Effector Nucleases (TALENs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR/Cas9) systems have revolutionised in the field of genome editing. These have been promptly accepted in both the research and commercial industry. On the whole, the widespread use of molecular markers with their types, their appliance in plant breeding along with the advances in genetic selection and genome editing together being a novel strategy to boost crop yield has been reviewed.

Genetic Diversity, Population Structure and Marker-Trait Association for 100-Seed Weight in International Safflower Panel Using SilicoDArT Marker Information

Safflower is an important oilseed crop mainly grown in the arid and semi-arid regions of the world. The aim of this study was to explore phenotypic and genetic diversity, population structure, and marker-trait association for 100-seed weight in 94 safflower accessions originating from 26 countries using silicoDArT markers. Analysis of variance revealed statistically significant genotypic effects (p < 0.01), while Turkey samples resulted in higher 100-seed weight compared to Pakistan samples. A Constellation plot divided the studied germplasm into two populations on the basis of their 100-seed weight. Various mean genetic diversity parameters including observed number of alleles (1.99), effective number of alleles (1.54), Shannon's information index (0.48), expected heterozygosity (0.32), and unbiased expected heterozygosity (0.32) for the entire population exhibited sufficient genetic diversity using 12232 silicoDArT markers. Analysis of molecular variance (AMOVA) revealed that most of the variations (91%) in world safflower panel are due to differences within country groups. A model-based structure grouped the 94 safflower accessions into populations A, B, C and an admixture population upon membership coefficient. Neighbor joining analysis grouped the safflower accessions into two populations (A and B). Principal coordinate analysis (PCoA) also clustered the safflower accessions on the basis of geographical origin. Three accessions; Egypt-5, Egypt-2, and India-2 revealed the highest genetic distance and hence might be recommended as candidate parental lines for safflower breeding programs. The mixed linear model i.e., the Q + K model, demonstrated that two DArTseq markers (DArT-45483051 and DArT-15672391) had significant association (p < 0.01) for 100-seed weight. We envisage that identified DArTseq markers associated with 100-seed weight will be helpful to develop high-yielding cultivars of safflower through marker-assisted breeding in the near future.

Genotyping by sequencing provides new insights into the diversity of Napier grass (Cenchrus purpureus) and reveals variation in genome-wide LD patterns between collections

Napier grass is an important tropical forage-grass and of growing potential as an energy crop. One-hundred-five Napier grass accessions, encompassing two independent collections, were subjected to genotyping by sequencing which generated a set of high-density genome-wide markers together with short sequence reads. The reads, averaging 54 nucleotides, were mapped to the pearl millet genome and the closest genes and annotation information were used to select candidate genes linked to key forage traits. 980 highly polymorphic SNP markers, distributed across the genome, were used to assess population structure and diversity with seven-subgroups identified. A few representative accessions were selected with the objective of distributing subsets of a manageable size for further evaluation. Genome-wide linkage disequilibrium (LD) analyses revealed a fast LD-decay, on average 2.54 kbp, in the combined population with a slower LD-decay in the ILRI collection compared with the EMBRAPA collection, the significance of which is discussed. This initiative generated high-density markers with a good distribution across the genome. The diversity analysis revealed the existence of a substantial amount of variation in the ILRI collection and identified some unique materials from the EMBRAPA collection, demonstrating the potential of the overall population for further genetic and marker-trait-association studies.

Genetic dissection of the mechanism of flowering time based on an environmentally stable and specific QTL in Brassica napus

Flowering time is an important agronomic trait that is highly influenced by the environment. To elucidate the genetic mechanism of flowering time in rapeseed (Brassica napus L.), a genome-wide QTL analysis was performed in a doubled haploid population grown in winter, semi-winter and spring ecological conditions. Fifty-five consensus QTLs were identified after combining phenotype and genomic data, including 12 environment-stable QTLs and 43 environment-specific QTLs. Importantly, six major QTLs for flowering time were identified, of which two were considered environment-specific QTLs in spring ecological condition and four were considered environment-stable QTLs in winter and semi-winter ecological conditions. Through QTL comparison, 18 QTLs were colocalized with QTLs from six other published studies. Combining the candidate genes with their functional annotation, in 49 of 55 consensus QTLs, 151 candidate genes in B. napus corresponding to 95 homologous genes in Arabidopsis thaliana related to flowering were identified, including BnaC03g32910D (CO), BnaA02g12130D (FT) and BnaA03g13630D (FLC). Most of the candidate genes were involved in different flowering regulatory pathways. Based on re-sequencing and differences in sequence annotation between the two parents, we found that regions containing some candidate genes have numerous non-frameshift InDels and many non- synonymous mutations, which might directly lead to gene functional variation. Flowering time was negativly correlated with seed yield and thousand seed weight based on a QTL comparison of flowering time and seed yield traits, which has implications in breeding new early-maturing varieties of B. napus. Moreover, a putative flowering regulatory network was constructed, including the photoperiod, circadian clock, vernalization, autonomous and gibberellin pathways. Multiple copies of genes led to functional difference among the different copies of homologous genes, which also increased the complexity of the flowering regulatory networks. Taken together, the present results not only provide new insights into the genetic regulatory network underlying the control of flowering time but also improve our understanding of flowering time regulatory pathways in rapeseed.

Characterizing Palestinian snake melon (Cucumis melo var. flexuosus) germplasm diversity and structure using SNP and DArTseq markers

BACKGROUND

Crop landraces embody a source of beneficial genes potentially providing endurance to environmental stress and other agronomic qualities including yield. Our study included 88 snake melon accessions (Cucumis melo var. flexuosus) collected from 9 districts in the Palestinian West-Bank. These accessions represent four landraces of Palestinian snake melon: Green, and White Baladi, and Green, and White Sahouri.

RESULTS

This is the first report on successful application of genotyping by sequencing in snake melon. Nine thousand seven hundred fifty single-nucleotide polymorphism (SNP) and 7400 DArTseq genetic markers were employed to evaluate genetic biodiversity and population structure of Palestinian snake melon germplasm collection. Clustering based on neighbor-joining-analysis, principle coordinate and Bayesian model implemented in Structure showed that patterns of genetic diversity of snake melon landraces depends on their geographical source and unraveled the presence of two major local landraces (Sahouri, and Baladi) with accessions from each group clustering together. A significant correlation was observed between both types of markers in Mantel correlation test. A significant association between genetic and geographic matrices (P < 0.0001) was also detected. AMOVA indicated that majority of variation (90%) was due to the difference within accessions.

CONCLUSION

The Palestinian landraces seem to have unique genes that may allow the enhancement of the global snake melon gene pool and developments of the plant production worldwide. Our subsequent objective is to detect genotypes with promising qualities and to conduct association mapping studies concentrating on Fusarium-wilt resistance, yield, and environmental stresses.

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.

Genome-wide association analysis of agronomic traits in wheat under drought-stressed and non-stressed conditions

This study determined the population structure and genome-wide marker-trait association of agronomic traits of wheat for drought-tolerance breeding. Ninety-three diverse bread wheat genotypes were genotyped using the Diversity Arrays Technology sequencing (DArTseq) protocol. The number of days-to-heading (DTH), number of days-to-maturity (DTM), plant height (PHT), spike length (SPL), number of kernels per spike (KPS), thousand kernel weight (TKW) and grain yield (GYLD), assessed under drought-stressed and non-stressed conditions, were considered for the study. Population structure analysis and genome-wide association mapping were undertaken based on 16,383 silico DArTs loci with < 10% missing data. The population evaluated was grouped into nine distinct genetic structures. Inter-chromosomal linkage disequilibrium showed the existence of linkage decay as physical distance increased. A total of 62 significant (P < 0.001) marker-trait associations (MTAs) were detected explaining more than 20% of the phenotypic variation observed under both drought-stressed and non-stressed conditions. Significant (P < 0.001) MTA event(s) were observed for DTH, PHT, SPL, SPS, and KPS; under both stressed and non-stressed conditions, while additional significant (P < 0.05) associations were observed for TKW, DTM and GYLD under non-stressed condition. The MTAs reported in this population could be useful to initiate marker-assisted selection (MAS) and targeted trait introgression of wheat under drought-stressed and non-stressed conditions, and for fine mapping and cloning of the underlying genes and QTL.

Genome-wide Association Study Identifies New Loci for Resistance to Leptosphaeria maculans in Canola

Key message "We identified both quantitative and quantitative resistance loci to Leptosphaeria maculans, a fungal pathogen, causing blackleg disease in canola. Several genome-wide significant associations were detected at known and new loci for blackleg resistance. We further validated statistically significant associations in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance in canola." Blackleg, caused by Leptosphaeria maculans, is a significant disease which affects the sustainable production of canola (Brassica napus). This study reports a genome-wide association study based on 18,804 polymorphic SNPs to identify loci associated with qualitative and quantitative resistance to L. maculans. Genomic regions delimited with 694 significant SNP markers, that are associated with resistance evaluated using 12 single spore isolates and pathotypes from four canola stubble were identified. Several significant associations were detected at known disease resistance loci including in the vicinity of recently cloned Rlm2/LepR3 genes, and at new loci on chromosomes A01/C01, A02/C02, A03/C03, A05/C05, A06, A08, and A09. In addition, we validated statistically significant associations on A01, A07, and A10 in four genetic mapping populations, demonstrating that GWAS marker loci are indeed associated with resistance to L. maculans. One of the novel loci identified for the first time, Rlm12, conveys adult plant resistance and mapped within 13.2 kb from Arabidopsis R gene of TIR-NBS class. We showed that resistance loci are located in the vicinity of R genes of Arabidopsis thaliana and Brassica napus on the sequenced genome of B. napus cv. Darmor-bzh. Significantly associated SNP markers provide a valuable tool to enrich germplasm for favorable alleles in order to improve the level of resistance to L. maculans in canola.

Multi-environment QTL studies suggest a role for cysteine-rich protein kinase genes in quantitative resistance to blackleg disease in Brassica napus

BACKGROUND

Resistance to the blackleg disease of Brassica napus (canola/oilseed rape), caused by the hemibiotrophic fungal pathogen Leptosphaeria maculans, is determined by both race-specific resistance (R) genes and quantitative resistance loci (QTL), or adult-plant resistance (APR). While the introgression of R genes into breeding material is relatively simple, QTL are often detected sporadically, making them harder to capture in breeding programs. For the effective deployment of APR in crop varieties, resistance QTL need to have a reliable influence on phenotype in multiple environments and be well defined genetically to enable marker-assisted selection (MAS).

RESULTS

Doubled-haploid populations produced from the susceptible B. napus variety Topas and APR varieties AG-Castle and AV-Sapphire were analysed for resistance to blackleg in two locations over 3 and 4 years, respectively. Three stable QTL were detected in each population, with two loci appearing to be common to both APR varieties. Physical delineation of three QTL regions was sufficient to identify candidate defense-related genes, including a cluster of cysteine-rich receptor-like kinases contained within a 49 gene QTL interval on chromosome A01. Individual L. maculans isolates were used to define the physical intervals for the race-specific R genes Rlm3 and Rlm4 and to identify QTL common to both field studies and the cotyledon resistance response.

CONCLUSION

Through multi-environment QTL analysis we have identified and delineated four significant and stable QTL suitable for MAS of quantitative blackleg resistance in B. napus, and identified candidate genes which potentially play a role in quantitative defense responses to L. maculans.

Genome-wide association analyses reveal complex genetic architecture underlying natural variation for flowering time in canola

Optimum flowering time is the key to maximize canola production in order to meet global demand of vegetable oil, biodiesel and canola-meal. We reveal extensive variation in flowering time across diverse genotypes of canola under field, glasshouse and controlled environmental conditions. We conduct a genome-wide association study and identify 69 single nucleotide polymorphism (SNP) markers associated with flowering time, which are repeatedly detected across experiments. Several associated SNPs occur in clusters across the canola genome; seven of them were detected within 20 Kb regions of a priori candidate genes; FLOWERING LOCUS T, FRUITFUL, FLOWERING LOCUS C, CONSTANS, FRIGIDA, PHYTOCHROME B and an additional five SNPs were localized within 14 Kb of a previously identified quantitative trait loci for flowering time. Expression analyses showed that among FLC paralogs, BnFLC.A2 accounts for ~23% of natural variation in diverse accessions. Genome-wide association analysis for FLC expression levels mapped not only BnFLC.C2 but also other loci that contribute to variation in FLC expression. In addition to revealing the complex genetic architecture of flowering time variation, we demonstrate that the identified SNPs can be modelled to predict flowering time in diverse canola germplasm accurately and hence are suitable for genomic selection of adaptative traits in canola improvement programmes.

Genome-wide identification of SSR and SNP markers from the non-heading Chinese cabbage for comparative genomic analyses

BACKGROUND

Non-heading Chinese cabbage (NHCC), belonging to Brassica, is an important leaf vegetable in Asia. Although genetic analyses have been performed through conventional selection and breeding efforts, the domestication history of NHCC and the genetics underlying its morphological diversity remain unclear. Thus, the reliable molecular markers representative of the whole genome are required for molecular-assisted selection in NHCC.

RESULTS

A total of 20,836 simple sequence repeats (SSRs) were detected in NHCC, containing repeat types from mononucleotide to nonanucleotide. The average density was 62.93 SSRs/Mb. In gene regions, 5,435 SSRs were identified in 4,569 genes. A total of 5,008 primer pairs were designed, and 74 were randomly selected for validation. Among these, 60 (81.08%) were polymorphic in 18 Cruciferae. The number of polymorphic bands ranged from two to five, with an average of 2.70 for each primer. The average values of the polymorphism information content, observed heterozygosity, Hardy-Weinberg equilibrium, and Shannon's information index were 0.2970, 0.4136, 0.5706, and 0.5885, respectively. Four clusters were classified according to the unweighted pair-group method with arithmetic average cluster analysis of 18 genotypes. In addition, a total of 1,228,979 single nucleotide polymorphisms (SNPs) were identified in the NHCC through a comparison with the genome of Chinese cabbage, and the average SNP density in the whole genome was 4.33/Kb. The number of SNPs ranged from 341,939 to 591,586 in the 10 accessions, and the average heterozygous SNPs ratio was ~42.53%. All analyses showed these markers were high quality and reliable. Therefore, they could be used in the construction of a linkage map and for genetic diversity studies for NHCC in future.

CONCLUSIONS

This is the first systematic and comprehensive analysis and identification of SSRs in NHCC and 17 species. The development of a large number of SNP and SSR markers was successfully achieved for NHCC. These novel markers are valuable for constructing genetic linkage maps, comparative genome analysis, quantitative trait locus (QTL) mapping, genome-wide association studies, and marker-assisted selection in NHCC breeding system research.

Development and use of molecular markers: past and present

Molecular markers, due to their stability, cost-effectiveness and ease of use provide an immensely popular tool for a variety of applications including genome mapping, gene tagging, genetic diversity diversity, phylogenetic analysis and forensic investigations. In the last three decades, a number of molecular marker techniques have been developed and exploited worldwide in different systems. However, only a handful of these techniques, namely RFLPs, RAPDs, AFLPs, ISSRs, SSRs and SNPs have received global acceptance. A recent revolution in DNA sequencing techniques has taken the discovery and application of molecular markers to high-throughput and ultrahigh-throughput levels. Although, the choice of marker will obviously depend on the targeted use, microsatellites, SNPs and genotyping by sequencing (GBS) largely fulfill most of the user requirements. Further, modern transcriptomic and functional markers will lead the ventures onto high-density genetic map construction, identification of QTLs, breeding and conservation strategies in times to come in combination with other high throughput techniques. This review presents an overview of different marker technologies and their variants with a comparative account of their characteristic features and applications.

SNP markers-based map construction and genome-wide linkage analysis in Brassica napus

An Illumina Infinium array comprising 5306 single nucleotide polymorphism (SNP) markers was used to genotype 175 individuals of a doubled haploid population derived from a cross between Skipton and Ag-Spectrum, two Australian cultivars of rapeseed (Brassica napus L.). A genetic linkage map based on 613 SNP and 228 non-SNP (DArT, SSR, SRAP and candidate gene markers) covering 2514.8 cM was constructed and further utilized to identify loci associated with flowering time and resistance to blackleg, a disease caused by the fungus Leptosphaeria maculans. Comparison between genetic map positions of SNP markers and the sequenced Brassica rapa (A) and Brassica oleracea (C) genome scaffolds showed several genomic rearrangements in the B. napus genome. A major locus controlling resistance to L. maculans was identified at both seedling and adult plant stages on chromosome A07. QTL analyses revealed that up to 40.2% of genetic variation for flowering time was accounted for by loci having quantitative effects. Comparative mapping showed Arabidopsis and Brassica flowering genes such as Phytochrome A/D, Flowering Locus C and agamous-Like MADS box gene AGL1 map within marker intervals associated with flowering time in a DH population from Skipton/Ag-Spectrum. Genomic regions associated with flowering time and resistance to L. maculans had several SNP markers mapped within 10 cM. Our results suggest that SNP markers will be suitable for various applications such as trait introgression, comparative mapping and high-resolution mapping of loci in B. napus.

Quantitative trait loci for thermal time to flowering and photoperiod responsiveness discovered in summer annual-type Brassica napus L

Time of flowering is a key adaptive trait in plants and is conditioned by the interaction of genes and environmental cues including length of photoperiod, ambient temperature and vernalisation. Here we investigated the photoperiod responsiveness of summer annual-types of Brassica napus (rapeseed, canola). A population of 131 doubled haploid lines derived from a cross between European and Australian parents was evaluated for days to flowering, thermal time to flowering (measured in degree-days) and the number of leaf nodes at flowering in a compact and efficient glasshouse-based experiment with replicated short and long day treatments. All three traits were under strong genetic control with heritability estimates ranging from 0.85–0.93. There was a very strong photoperiod effect with flowering in the population accelerated by 765 degree-days in the long day versus short day treatments. However, there was a strong genetic correlation of line effects (0.91) between the long and short day treatments and relatively low genotype x treatment interaction indicating that photoperiod had a similar effect across the population. Bivariate analysis of thermal time to flowering in short and long days revealed three main effect quantitative trait loci (QTLs) that accounted for 57.7% of the variation in the population and no significant interaction QTLs. These results provided insight into the contrasting adaptations of Australian and European varieties. Both parents responded to photoperiod and their alleles shifted the population to earlier flowering under long days. In addition, segregation of QTLs in the population caused wide transgressive segregation in thermal time to flowering. Potential candidate flowering time homologues located near QTLs were identified with the aid of the Brassica rapa reference genome sequence. We discuss how these results will help to guide the breeding of summer annual types of B. napus adapted to new and changing environments.

Genome-wide delineation of natural variation for pod shatter resistance in Brassica napus

Resistance to pod shattering (shatter resistance) is a target trait for global rapeseed (canola, Brassica napus L.), improvement programs to minimise grain loss in the mature standing crop, and during windrowing and mechanical harvest. We describe the genetic basis of natural variation for shatter resistance in B. napus and show that several quantitative trait loci (QTL) control this trait. To identify loci underlying shatter resistance, we used a novel genotyping-by-sequencing approach DArT-Seq. QTL analysis detected a total of 12 significant QTL on chromosomes A03, A07, A09, C03, C04, C06, and C08; which jointly account for approximately 57% of the genotypic variation in shatter resistance. Through Genome-Wide Association Studies, we show that a large number of loci, including those that are involved in shattering in Arabidopsis, account for variation in shatter resistance in diverse B. napus germplasm. Our results indicate that genetic diversity for shatter resistance genes in B. napus is limited; many of the genes that might control this trait were not included during the natural creation of this species, or were not retained during the domestication and selection process. We speculate that valuable diversity for this trait was lost during the natural creation of B. napus. To improve shatter resistance, breeders will need to target the introduction of useful alleles especially from genotypes of other related species of Brassica, such as those that we have identified.

Constructing a dense genetic linkage map and mapping QTL for the traits of flower development in Brassica carinata

An integrated dense genetic linkage map was constructed in a B. carinata population and used for comparative genome analysis and QTL identification for flowering time. An integrated dense linkage map of Brassica carinata (BBCC) was constructed in a doubled haploid population based on DArT-Seq(TM) markers. A total of 4,031 markers corresponding to 1,366 unique loci were mapped including 639 bins, covering a genetic distance of 2,048 cM. We identified 136 blocks and islands conserved in Brassicaceae, which showed a feature of hexaploidisation representing the suggested ancestral crucifer karyotype. The B and C genome of B. carinata shared 85 % of commonly conserved blocks with the B genome of B. nigra/B. juncea and 80 % of commonly conserved blocks with the C genome of B. napus, and shown frequent structural rearrangements such as insertions and inversions. Up to 24 quantitative trait loci (QTL) for flowering and budding time were identified in the DH population. Of these QTL, one consistent QTL (qFT.B4-2) for flowering time was identified in all of the environments in the J block of the B4 linkage group, where a group of genes for flowering time were aligned in A. thaliana. Another major QTL for flowering time under a winter-cropped environment was detected in the E block of C6, where the BnFT-C6 gene was previously localised in B. napus. This high-density map would be useful not only to reveal the genetic variation in the species with QTL analysis and genome sequencing, but also for other applications such as marker-assisted selection and genomic selection, for the African mustard improvement.

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.

New evidence of ancestral polyploidy in the Genistoid legume Lupinus angustifolius L. (narrow-leafed lupin)

This is the first clear evidence of duplication and/or triplication of large chromosomal regions in a genome of a Genistoid legume, the most basal clade of Papilionoid legumes. Lupinus angustifolius L. (narrow-leafed lupin) is the most widely cultivated species of Genistoid legume, grown for its high-protein grain. As a member of this most basal clade of Papilionoid legumes, L. angustifolius serves as a useful model for exploring legume genome evolution. Here, we report an improved reference genetic map of L. angustifolius comprising 1207 loci, including 299 newly developed Diversity Arrays Technology markers and 54 new gene-based PCR markers. A comparison between the L. angustifolius and Medicago truncatula genomes was performed using 394 sequence-tagged site markers acting as bridging points between the two genomes. The improved L. angustifolius genetic map, the updated M. truncatula genome assembly and the increased number of bridging points between the genomes together substantially enhanced the resolution of synteny and chromosomal colinearity between these genomes compared to previous reports. While a high degree of syntenic fragmentation was observed that was consistent with the large evolutionary distance between the L. angustifolius and M. truncatula genomes, there were striking examples of conserved colinearity of loci between these genomes. Compelling evidence was found of large-scale duplication and/or triplication in the L. angustifolius genome, consistent with one or more ancestral polyploidy events.

Quantitative trait loci for thermal time to flowering and photoperiod responsiveness discovered in summer annual-type Brassica napus L

Time of flowering is a key adaptive trait in plants and is conditioned by the interaction of genes and environmental cues including length of photoperiod, ambient temperature and vernalisation. Here we investigated the photoperiod responsiveness of summer annual-types of Brassica napus (rapeseed, canola). A population of 131 doubled haploid lines derived from a cross between European and Australian parents was evaluated for days to flowering, thermal time to flowering (measured in degree-days) and the number of leaf nodes at flowering in a compact and efficient glasshouse-based experiment with replicated short and long day treatments. All three traits were under strong genetic control with heritability estimates ranging from 0.85-0.93. There was a very strong photoperiod effect with flowering in the population accelerated by 765 degree-days in the long day versus short day treatments. However, there was a strong genetic correlation of line effects (0.91) between the long and short day treatments and relatively low genotype x treatment interaction indicating that photoperiod had a similar effect across the population. Bivariate analysis of thermal time to flowering in short and long days revealed three main effect quantitative trait loci (QTLs) that accounted for 57.7% of the variation in the population and no significant interaction QTLs. These results provided insight into the contrasting adaptations of Australian and European varieties. Both parents responded to photoperiod and their alleles shifted the population to earlier flowering under long days. In addition, segregation of QTLs in the population caused wide transgressive segregation in thermal time to flowering. Potential candidate flowering time homologues located near QTLs were identified with the aid of the Brassica rapa reference genome sequence. We discuss how these results will help to guide the breeding of summer annual types of B. napus adapted to new and changing environments.

A high-density SNP map for accurate mapping of seed fibre QTL in Brassica napus L

A high density genetic linkage map for the complex allotetraploid crop species Brassica napus (oilseed rape) was constructed in a late-generation recombinant inbred line (RIL) population, using genome-wide single nucleotide polymorphism (SNP) markers assayed by the Brassica 60 K Infinium BeadChip Array. The linkage map contains 9164 SNP markers covering 1832.9 cM. 1232 bins account for 7648 of the markers. A subset of 2795 SNP markers, with an average distance of 0.66 cM between adjacent markers, was applied for QTL mapping of seed colour and the cell wall fiber components acid detergent lignin (ADL), cellulose and hemicellulose. After phenotypic analyses across four different environments a total of 11 QTL were detected for seed colour and fiber traits. The high-density map considerably improved QTL resolution compared to the previous low-density maps. A previously identified major QTL with very high effects on seed colour and ADL was pinpointed to a narrow genome interval on chromosome A09, while a minor QTL explaining 8.1% to 14.1% of variation for ADL was detected on chromosome C05. Five and three QTL accounting for 4.7% to 21.9% and 7.3% to 16.9% of the phenotypic variation for cellulose and hemicellulose, respectively, were also detected. To our knowledge this is the first description of QTL for seed cellulose and hemicellulose in B. napus, representing interesting new targets for improving oil content. The high density SNP genetic map enables navigation from interesting B. napus QTL to Brassica genome sequences, giving useful new information for understanding the genetics of key seed quality traits in rapeseed.

Construction of integrated linkage map of a recombinant inbred line population of white lupin (Lupinus albus L.)

We report the development of a Diversity Arrays Technology (DArT) marker panel and its utilisation in the development of an integrated genetic linkage map of white lupin (Lupinus albus L.) using an F8 recombinant inbred line population derived from Kiev Mutant/P27174. One hundred and thirty-six DArT markers were merged into the first genetic linkage map composed of 220 amplified fragment length polymorphisms (AFLPs) and 105 genic markers. The integrated map consists of 38 linkage groups of 441 markers and spans a total length of 2,169 cM, with an average interval size of 4.6 cM. The DArT markers exhibited good genome coverage and were associated with previously identified genic and AFLP markers linked with quantitative trait loci for anthracnose resistance, flowering time and alkaloid content. The improved genetic linkage map of white lupin will aid in the identification of markers for traits of interest and future syntenic studies.

A consensus map of rapeseed (Brassica napus L.) based on diversity array technology markers: applications in genetic dissection of qualitative and quantitative traits

Background

Dense consensus genetic maps based on high-throughput genotyping platforms are valuable for making genetic gains in Brassica napus through quantitative trait locus identification, efficient predictive molecular breeding, and map-based gene cloning. This report describes the construction of the first B. napus consensus map consisting of a 1,359 anchored array based genotyping platform; Diversity Arrays Technology (DArT), and non-DArT markers from six populations originating from Australia, Canada, China and Europe. We aligned the B. napus DArT sequences with genomic scaffolds from Brassica rapa and Brassica oleracea, and identified DArT loci that showed linkage with qualitative and quantitative loci associated with agronomic traits.

Results

The integrated consensus map covered a total of 1,987.2 cM and represented all 19 chromosomes of the A and C genomes, with an average map density of one marker per 1.46 cM, corresponding to approximately 0.88 Mbp of the haploid genome. Through in silico physical mapping 2,457 out of 3,072 (80%) DArT clones were assigned to the genomic scaffolds of B. rapa (A genome) and B. oleracea (C genome). These were used to orientate the genetic consensus map with the chromosomal sequences. The DArT markers showed linkage with previously identified non-DArT markers associated with qualitative and quantitative trait loci for plant architecture, phenological components, seed and oil quality attributes, boron efficiency, sucrose transport, male sterility, and race-specific resistance to blackleg disease.

Conclusions

The DArT markers provide increased marker density across the B. napus genome. Most of the DArT markers represented on the current array were sequenced and aligned with the B. rapa and B. oleracea genomes, providing insight into the Brassica A and C genomes. This information can be utilised for comparative genomics and genomic evolution studies. In summary, this consensus map can be used to (i) integrate new generation markers such as SNP arrays and next generation sequencing data; (ii) anchor physical maps to facilitate assembly of B. napus genome sequences; and (iii) identify candidate genes underlying natural genetic variation for traits of interest.

Genetic and physical mapping of flowering time loci in canola (Brassica napus L.)

We identified quantitative trait loci (QTL) underlying variation for flowering time in a doubled haploid (DH) population of vernalisation-responsive canola (Brassica napus L.) cultivars Skipton and Ag-Spectrum and aligned them with physical map positions of predicted flowering genes from the Brassica rapa genome. Significant genetic variation in flowering time and response to vernalisation were observed among the DH lines from Skipton/Ag-Spectrum. A molecular linkage map was generated comprising 674 simple sequence repeat, sequence-related amplified polymorphism, sequence characterised amplified region, Diversity Array Technology, and candidate gene based markers loci. QTL analysis indicated that flowering time is a complex trait and is controlled by at least 20 loci, localised on ten different chromosomes. These loci each accounted for between 2.4 and 28.6% of the total genotypic variation for first flowering and response to vernalisation. However, identification of consistent QTL was found to be dependant upon growing environments. We compared the locations of QTL with the physical positions of predicted flowering time genes located on the sequenced genome of B. rapa. Some QTL associated with flowering time on A02, A03, A07, and C06 may represent homologues of known flowering time genes in Arabidopsis; VERNALISATION INSENSITIVE 3, APETALA1, CAULIFLOWER, FLOWERING LOCUS C, FLOWERING LOCUS T, CURLY LEAF, SHORT VEGETATIVE PHASE, GA3 OXIDASE, and LEAFY. Identification of the chromosomal location and effect of the genes influencing flowering time may hasten the development of canola varieties having an optimal time for flowering in target environments such as for low rainfall areas, via marker-assisted selection.

Genetic and physical mapping of flowering time loci in canola (Brassica napus L.)

We identified quantitative trait loci (QTL) underlying variation for flowering time in a doubled haploid (DH) population of vernalisation-responsive canola (Brassica napus L.) cultivars Skipton and Ag-Spectrum and aligned them with physical map positions of predicted flowering genes from the Brassica rapa genome. Significant genetic variation in flowering time and response to vernalisation were observed among the DH lines from Skipton/Ag-Spectrum. A molecular linkage map was generated comprising 674 simple sequence repeat, sequence-related amplified polymorphism, sequence characterised amplified region, Diversity Array Technology, and candidate gene based markers loci. QTL analysis indicated that flowering time is a complex trait and is controlled by at least 20 loci, localised on ten different chromosomes. These loci each accounted for between 2.4 and 28.6% of the total genotypic variation for first flowering and response to vernalisation. However, identification of consistent QTL was found to be dependant upon growing environments. We compared the locations of QTL with the physical positions of predicted flowering time genes located on the sequenced genome of B. rapa. Some QTL associated with flowering time on A02, A03, A07, and C06 may represent homologues of known flowering time genes in Arabidopsis; VERNALISATION INSENSITIVE 3, APETALA1, CAULIFLOWER, FLOWERING LOCUS C, FLOWERING LOCUS T, CURLY LEAF, SHORT VEGETATIVE PHASE, GA3 OXIDASE, and LEAFY. Identification of the chromosomal location and effect of the genes influencing flowering time may hasten the development of canola varieties having an optimal time for flowering in target environments such as for low rainfall areas, via marker-assisted selection.

A genetic linkage map of Brassica carinata constructed with a doubled haploid population

Brassica carinata is an important oilseed crop with unique favourable traits that are desirable for other Brassica crops. However, given the limited research into genetic resources in B. carinata, knowledge of the genetic structure of this species is relatively poor. Nine homozygous, genetically distinct accessions of B. carinata were obtained via microspore culture, from which two divergent doubled haploid (DH) lines were used to develop a DH mapping population that consisted of 183 lines. The mapping population showed segregation of multiple traits of interest. A genetic map was constructed with PCR-based markers, and a total of 212 loci, which covered 1,703 cM, were assigned to eight linkage groups in the B genome and nine linkage groups in the C genome, which allowed comparison with genetic maps of other important Brassica species that contain the B/C genome(s). Loci for two Mendelian-inherited traits related to pigmentation (petal and anther tip colour) and one quantitative trait (seed coat colour) were identified using the linkage map. The significance of the mapping population in the context of genetic improvement of Brassica crops is discussed.

Comparative analysis of FLC homologues in Brassicaceae provides insight into their role in the evolution of oilseed rape

We identified nine FLOWERING LOCUS C homologues (BnFLC) in Brassica napus and found that the coding sequences of all BnFLCs were relatively conserved but the intronic and promoter regions were more divergent. The BnFLC homologues were mapped to six of 19 chromosomes. All of the BnFLC homologues were located in the collinear region of FLC in the Arabidopsis genome except BnFLC.A3b and BnFLC.C3b, which were mapped to noncollinear regions of chromosome A3 and C3, respectively. Four of the homologues were associated significantly with quantitative trait loci for flowering time in two mapping populations. The BnFLC homologues showed distinct expression patterns in vegetative and reproductive organs, and at different developmental stages. BnFLC.A3b was differentially expressed between the winter-type and semi-winter-type cultivars. Microsynteny analysis indicated that BnFLC.A3b might have been translocated to the present segment in a cluster with other flowering-time regulators, such as a homologue of FRIGIDA in Arabidopsis. This cluster of flowering-time genes might have conferred a selective advantage to Brassica species in terms of increased adaptability to diverse environments during their evolution and domestication process.