Genome-wide SNP identification and QTL mapping for black rot resistance in cabbage

Study of cabbage antioxidant system response on early infection stage of Xanthomonas campestris pv. campestris

Black rot, caused by Xanthomonas campestris pv. campestris (Xcc) significantly affects the production of cabbage and other cruciferous vegetables. Plant antioxidant system plays an important role in pathogen invasion and is one of the main mechanisms underlying resistance to biological stress. Therefore, it is important to study the resistance mechanisms of the cabbage antioxidant system during the early stages of Xcc. In this study, 108 CFU/mL (OD600 = 0.1) Xcc race1 was inoculated on "zhonggan 11" cabbage using the spraying method. The effects of Xcc infection on the antioxidant system before and after Xcc inoculation (0, 1, 3, and 5 d) were studied by physiological indexes determination, transcriptome and metabolome analyses. We concluded that early Xcc infection can destroy the balance of the active oxygen metabolism system, increase the generation of free radicals, and decrease the scavenging ability, leading to membrane lipid peroxidation, resulting in the destruction of the biofilm system and metabolic disorders. In response to Xcc infection, cabbage clears a series of reactive oxygen species (ROS) produced during Xcc infection via various antioxidant pathways. The activities of antioxidant enzymes such as superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) increased after Xcc infection, and the ROS scavenging rate increased. The biosynthesis of non-obligate antioxidants, such as ascorbic acid (AsA) and glutathione (GSH), is also enhanced after Xcc infection. Moreover, the alkaloid and vitamin contents increased significantly after Xcc infection. We concluded that cabbage could resist Xcc invasion by maintaining the stability of the cell membrane system and improving the biosynthesis of antioxidant substances and enzymes after infection by Xcc. Our results provide theoretical basis and data support for subsequent research on the cruciferous vegetables resistance mechanism and breeding to Xcc.

Characterisation of New Foxunavirus Phage Murka with the Potential of Xanthomonas campestris pv. campestris Control

Phages of phytopathogenic bacteria are considered to be promising agents for the biological control of bacterial diseases in plants. This paper reports on the isolation and characterisation of a new Xanthomonas campestris pv. campestris phage, Murka. Phage morphology and basic kinetic characteristics of the infection were determined, and a phylogenomic analysis was performed. The phage was able to lyse a reasonably broad range (64%, 9 of the 14 of the Xanthomonas campestris pv. campestris strains used in the study) of circulating strains of the cabbage black rot pathogen. This lytic myovirus has a DNA genome of 44,044 bp and contains 83 predicted genes. Taxonomically, it belongs to the genus Foxunavirus. This bacteriophage is promising for use as a possible means of biological control of cabbage black rot.

Molecular Marker Development for the Rapid Differentiation of Black Rot Causing Xanthomonas campestris pv. campestris Race 7

Xanthomonas campestris pv. campestris (Xcc) is a plant pathogen of Brassica crops that causes black rot disease throughout the world. At present, 11 physiological races of Xcc (races 1-11) have been reported. The conventional method of using differential cultivars for Xcc race detection is not accurate and it is laborious and time-consuming. Therefore, the development of specific molecular markers has been used as a substitute tool because it offers an accurate and reliable result, particularly a quick diagnosis of Xcc races. Previously, our laboratory has successfully developed race-specific molecular markers for Xcc races 1-6. In this study, specific molecular markers to identify Xcc race 7 have been developed. In the course of study, whole genome sequences of several Xcc races, X. campestris pv. incanae, X. campestris pv. raphani, and X. campestris pv. vesicatoria were aligned to identify variable regions like sequence-characterized amplified regions and insertions and deletions specific to race 7. Primer pairs were designed targeting these regions and validated against 22 samples. The polymerase chain reaction analysis revealed that three primer pairs specifically amplified the DNA fragment corresponding to race 7. The obtained finding clearly demonstrates the efficiency of the newly developed markers in accurately detecting Xcc race 7 among the other races. These results indicated that the newly developed marker can successfully and rapidly detect Xcc race 7 from other races. This study represents the first report on the successful development of specific molecular markers for Xcc race 7.

Climate and Soil Microsite Conditions Determine Local Adaptation in Declining Silver Fir Forests

Ongoing climatic change is threatening the survival of drought-sensitive tree species, such as silver fir (Abies alba). Drought-induced dieback had been previously explored in this conifer, although the role played by tree-level genetic diversity and its relationship with growth patterns and soil microsite conditions remained elusive. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to describe different genetic characteristics of five silver fir forests in the Spanish Pyrenees, including declining and non-declining trees. Single nucleotide polymorphisms (SNPs) were used to investigate the relationships between genetics, dieback, intraspecific trait variation (functional dendrophenotypic traits and leaf traits), local bioclimatic conditions, and rhizosphere soil properties. While there were no noticeable genetic differences between declining and non-declining trees, genome-environment associations with selection signatures were abundant, suggesting a strong influence of climate, soil physicochemical properties, and soil microbial diversity on local adaptation. These results provide novel insights into how genetics and diverse environmental factors are interrelated and highlight the need to incorporate genetic data into silver fir forest dieback studies to gain a better understanding of local adaptation.

A GBS-based genetic linkage map and quantitative trait loci (QTL) associated with resistance to Xanthomonas campestris pv. campestris race 1 identified in Brassica oleracea

The production of Brassica oleracea, an important vegetable crop, is severely affected by black rot disease caused by the bacterial pathogen Xanthomonas campestris pv. campestris. Resistance to race 1, the most virulent and widespread race in B. oleracea, is under quantitative control; therefore, identifying the genes and genetic markers associated with resistance is crucial for developing resistant cultivars. Quantitative trait locus (QTL) analysis of resistance in the F₂ population developed by crossing the resistant parent BR155 with the susceptible parent SC31 was performed. Sequence GBS approach was used to develop a genetic linkage map. The map contained 7,940 single nucleotide polymorphism markers consisting of nine linkage groups spanning 675.64 cM with an average marker distance of 0.66 cM. The F_2:3 population (N = 126) was evaluated for resistance to black rot disease in summer (2020), fall (2020), and spring (2021). QTL analysis, using a genetic map and phenotyping data, identified seven QTLs with LOD values between 2.10 and 4.27. The major QTL, qCaBR1, was an area of overlap between the two QTLs identified in the 2^nd and 3^rd trials located at C06. Among the genes located in the major QTL interval, 96 genes had annotation results, and eight were found to respond to biotic stimuli. We compared the expression patterns of eight candidate genes in susceptible (SC31) and resistant (BR155) lines using qRT-PCR and observed their early and transient increases or suppression in response to Xanthomonas campestris pv. campestris inoculation. These results support the involvement of the eight candidate genes in black rot resistance. The findings of this study will contribute towards marker-assisted selection, additionally the functional analysis of candidate genes may elucidate the molecular mechanisms underlying black rot resistance in B. oleracea.

Whole-genome resequencing of Coffea arabica L. (Rubiaceae) genotypes identify SNP and unravels distinct groups showing a strong geographical pattern

BACKGROUND

Coffea arabica L. is an economically important agricultural crop and the most popular beverage worldwide. As a perennial crop with recalcitrant seed, conservation of the genetic resources of coffee can be achieved through the complementary approach of in-situ and ex-situ field genebank. In Ethiopia, a large collection of C. arabica L. germplasm is preserved in field gene banks. Here, we report the whole-genome resequencing of 90 accessions from Choche germplasm bank representing garden and forest-based coffee production systems using Illumina sequencing technology.

RESULTS

The genome sequencing generated 6.41 billion paired-end reads, with a mean of 71.19 million reads per sample. More than 93% of the clean reads were mapped onto the C. arabica L. reference genome. A total of 11.08 million variants were identified, among which 9.74 million (87.9%) were SNPs (Single nucleotide polymorphisms) and 1.34 million (12.1%) were InDels. In all accessions, genomic variants were unevenly distributed across the coffee genome. The phylogenetic analysis using the SNP markers displayed distinct groups.

CONCLUSIONS

Resequencing of the coffee accessions has allowed identification of genetic markers, such as SNPs and InDels. The SNPs discovered in this study might contribute to the variation in important pathways of genes for important agronomic traits such as caffeine content, yield, disease, and pest in coffee. Moreover, the genome resequencing data and the genetic markers identified from 90 accessions provide insight into the genetic variation of the coffee germplasm and facilitate a broad range of genetic studies.

Advances in Multi-Omics Approaches for Molecular Breeding of Black Rot Resistance in Brassica oleracea L

Brassica oleracea is one of the most important species of the Brassicaceae family encompassing several economically important vegetables produced and consumed worldwide. But its sustainability is challenged by a range of pathogens, among which black rot, caused by Xanthomonas campestris pv. campestris (Xcc), is the most serious and destructive seed borne bacterial disease, causing huge yield losses. Host-plant resistance could act as the most effective and efficient solution to curb black rot disease for sustainable production of B. oleracea. Recently, 'omics' technologies have emerged as promising tools to understand the host-pathogen interactions, thereby gaining a deeper insight into the resistance mechanisms. In this review, we have summarized the recent achievements made in the emerging omics technologies to tackle the black rot challenge in B. oleracea. With an integrated approach of the omics technologies such as genomics, proteomics, transcriptomics, and metabolomics, it would allow better understanding of the complex molecular mechanisms underlying black rot resistance. Due to the availability of sequencing data, genomics and transcriptomics have progressed as expected for black rot resistance, however, other omics approaches like proteomics and metabolomics are lagging behind, necessitating a holistic and targeted approach to address the complex questions of Xcc-Brassica interactions. Genomic studies revealed that the black rot resistance is a complex trait and is mostly controlled by quantitative trait locus (QTL) with minor effects. Transcriptomic analysis divulged the genes related to photosynthesis, glucosinolate biosynthesis and catabolism, phenylpropanoid biosynthesis pathway, ROS scavenging, calcium signalling, hormonal synthesis and signalling pathway are being differentially expressed upon Xcc infection. Comparative proteomic analysis in relation to susceptible and/or resistance interactions with Xcc identified the involvement of proteins related to photosynthesis, protein biosynthesis, processing and degradation, energy metabolism, innate immunity, redox homeostasis, and defence response and signalling pathways in Xcc-Brassica interaction. Specifically, most of the studies focused on the regulation of the photosynthesis-related proteins as a resistance response in both early and later stages of infection. Metabolomic studies suggested that glucosinolates (GSLs), especially aliphatic and indolic GSLs, its subsequent hydrolysis products, and defensive metabolites synthesized by jasmonic acid (JA)-mediated phenylpropanoid biosynthesis pathway are involved in disease resistance mechanisms against Xcc in Brassica species. Multi-omics analysis showed that JA signalling pathway is regulating resistance against hemibiotrophic pathogen like Xcc. So, the bonhomie between omics technologies and plant breeding is going to trigger major breakthroughs in the field of crop improvement by developing superior cultivars with broad-spectrum resistance. If multi-omics tools are implemented at the right scale, we may be able to achieve the maximum benefits from the minimum. In this review, we have also discussed the challenges, future prospects, and the way forward in the application of omics technologies to accelerate the breeding of B. oleracea for disease resistance. A deeper insight about the current knowledge on omics can offer promising results in the breeding of high-quality disease-resistant crops.

Inheritance of Black Rot Resistance and Development of Molecular Marker Linked to Xcc Races 6 and 7 Resistance in Cabbage

Black rot, caused by Xanthomonas campestris pv. campestris (Xcc), produces V-shaped chlorotic lesions on the leaves of cabbage (Brassica oleracea var. capitata L.), causing darkened veins and drastically reducing yield and quality. Of the 11 Xcc races identified, races 1, 4, and 6 are predominant globally. In the present study, we aimed to develop a molecular marker linked to black rot resistance against Xcc races 6 and 7. Crossed between black rot-resistant ('SCNU-C-3470') and -susceptible ('SCNU-C-3328') lines obtained 186 F₂ plants. Resistance to Xcc race 6 segregated in a 3:1 (susceptible:resistant) ratio in the F₂ population, which is consistent with a monogenic recessive trait. Nucleotide-binding site (NBS) leucine rich repeat (LRR)-encoding resistance (R) genes play a crucial role in plant defenses to various pathogens. The candidate R gene (Bol031422) located on chromosome C08, previously reported by our research group, was cloned and sequenced in resistant and susceptible cabbage lines. The R gene Bol031422 consisted of a single exon with a 3 bp insertion/deletions (InDels), a 292 bp polymorphism (an insertion in the exon of the resistant line relative to the susceptible line) and several single nucleotide polymorphisms (SNPs). Here, we developed the InDel marker BR6-InDel to assess linkage between variation at Bol031422 and resistance to Xcc races 6 and 7. This marker will help cabbage breeders develop cabbage cultivars resistant to Xcc races 6 and 7.

Marker-Assisted Pyramiding of Downy Mildew-Resistant Gene Ppa3 and Black Rot-Resistant Gene Xca1bo in Popular Early Cauliflower Variety Pusa Meghna

Cauliflower is an important extensively grown cool season vegetable in India. Black rot and downy mildew are major devastating diseases reducing yield and quality of the crop. To tackle these through host plant resistance, a marker-assisted backcross breeding method was followed to pyramid a black rot-resistant gene (Xca1bo) and a downy mildew-resistant gene (Ppa3) from donors BR-161 and BR-2, respectively, into the background of Pusa Meghna cauliflower cultivar. Marker-assisted backcross breeding was followed up to BC₂ generation using SCAR marker ScOPO-04₈₃₃ and SSR marker BoGMS0624 for black rot and downy mildew resistance genes in foreground selection, respectively. In background selection, at each stage of backcrossing, 47 parental polymorphic SSR markers were used. The graphical genotyping of the five two-gene (Xca1boXca1boPpa3Ppa3) homozygous BC₂F₂ plants showed an average recovery of 85.44% of the Pusa Meghna genome with highest genome recovery of 91.7%. The genome contribution of donor parents (BR-161 and BR-2) was 8.26 with 6.34% of residual heterozygousity. The backcross derived pyramided lines BC₂F_2:3-7-16 and BC₂F_2:3-7-33 showed high resistance to both the diseases and exhibited higher yield and vitamin C content as compared with recipient parent Pusa Meghna. It is, therefore, evident from this study that resistant genes can be introgressed successfully into a Pusa Meghna cultivar without any yield penalty, benefitting farmers with reduced input cost and consumers with chemical residue free produce. Besides, the pyramided lines carrying dominant resistant genes can be exploited in a hybridization programme to develop hybrid(s) in cauliflower.

Genetic Diversity and Genome-Wide Association Study of Seed Aspect Ratio Using a High-Density SNP Array in Peanut (Arachis hypogaea L.)

Peanut (Arachis hypogaea L.) is one of the important oil crops of the world. In this study, we aimed to evaluate the genetic diversity of 384 peanut germplasms including 100 Korean germplasms and 284 core collections from the United States Department of Agriculture (USDA) using an Axiom_Arachis array with 58K single-nucleotide polymorphisms (SNPs). We evaluated the evolutionary relationships among 384 peanut germplasms using a genome-wide association study (GWAS) of seed aspect ratio data processed by ImageJ software. In total, 14,030 filtered polymorphic SNPs were identified from the peanut 58K SNP array. We identified five SNPs with significant associations to seed aspect ratio on chromosomes Aradu.A09, Aradu.A10, Araip.B08, and Araip.B09. AX-177640219 on chromosome Araip.B08 was the most significantly associated marker in GAPIT and Regularization method. Phosphoenolpyruvate carboxylase (PEPC) was found among the eleven genes within a linkage disequilibrium (LD) of the significant SNPs on Araip.B08 and could have a strong causal effect in determining seed aspect ratio. The results of the present study provide information and methods that are useful for further genetic and genomic studies as well as molecular breeding programs in peanuts.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

An extended KASP-SNP resource for molecular breeding in Chinese cabbage(Brassica rapa L. ssp. pekinensis)

Kompetitive allele-specific PCR (KASP) is a cost-effective single-step SNP genotyping technology, With an objective to enhance the marker repertoire and develop high efficient KASP-SNP markers in Chinese cabbage, we re-sequenced four Chinese cabbage doubled haploid (DH) lines, Y177-47, Y635-10, Y510-1 and Y510-9, and generated a total of more than 38.5 billion clean base pairs. A total of 827,720 SNP loci were identified with an estimated density of 3,217 SNPs/Mb. Further, a total of 387,354 SNPs with at least 30 bp to the next most adjacent SNPs on either side were selected as resource for KASP markers. From this resource, 258 (96.27%) of 268 SNP loci were successfully transformed into KASP-SNP markers using a Roche LightCycler 480-II instrument. Among these markers, 221 (85.66%) were co-dominant markers, 220 (85.27%) were non-synonymous SNPs, and 257 (99.6%) were newly developed markers. In addition, 53 markers were applied for genotyping of 34 Brassica rapa accessions. Cluster analysis separated these 34 accessions into three clusters based on heading types. The millions of SNP loci, a large set of resource for KASP markers, as well as the newly developed KASP markers in this study may facilitate further genetic and molecular breeding studies in Brassica rapa.

Transcriptomic analysis of resistant and susceptible cabbage lines reveals differential expressions and candidate genes involved in cabbage early responses to black rot

Cabbage (Brassica oleracea var. capitata) is one of the most important cruciferous leafy vegetable crops and is widely cultivated all over the world. Its yield and quality are often affected by diseases such as cabbage black rot. 2R is a cabbage line that is newly resistant to black rot, which was created by radiation mutagenesis and backcross transfer. However, the underlying molecular bases and mechanisms of early-phase response of different resistant cabbage lines against black rot infections remain unknown. Here, we completed a comprehensive transcriptome profile analysis between resistant (2R) and susceptible (2T) cabbage lines after black rot inoculations. The results showed that the typical V-shaped lesions were found in inoculated plants after 15 days, and the symptoms in the susceptible cabbage lines (2T) were significant severe than that of the resistant 2R line. A total of 10,030 differentially expressed genes (DEGs) were identified, of which 384 DEGs were found to overlap in resistant and susceptible cabbage lines after black rot infections, suggesting those DEGs may play more important roles in cabbage early responses to black rot infections. We ranked the expression levels of DEGs among the four comparison sets of resistant and susceptible cabbage lines and, interestingly, found the top ten differential expression genes contained NBS-LRR genes, protein kinase genes and expansin genes. These findings provide a comprehensive differential transcriptome profile between resistant and susceptible cabbage lines and indicate some genes play key roles in the regulation of early response to black rot infections, which will help to understand the molecular resistance of cabbage against these infections.

An update on the arsenal: mining resistance genes for disease management of Brassica crops in the genomic era

Brassica species include many economically important crops that provide nutrition and health-promoting substances to humans worldwide. However, as with all crops, their production is constantly threatened by emerging viral, bacterial, and fungal diseases, whose incidence has increased in recent years. Traditional methods of control are often costly, present limited effectiveness, and cause environmental damage; instead, the ideal approach is to mine and utilize the resistance genes of the Brassica crop hosts themselves. Fortunately, the development of genomics, molecular genetics, and biological techniques enables us to rapidly discover and apply resistance (R) genes. Herein, the R genes identified in Brassica crops are summarized, including their mapping and cloning, possible molecular mechanisms, and application in resistance breeding. Future perspectives concerning how to accurately discover additional R gene resources and efficiently utilize these genes in the genomic era are also discussed.

Heterozygous Single Nucleotide Polymorphic Loci in Haploid Gametophytes of Gracilariopsis lemaneiformis (Rhodophyta)

Gracilariopsis lemaneiformis is an important commercial macroalga. Whole-genome resequencing was conducted to identify single nucleotide polymorphisms (SNPs) in parental gametophytes and 60 F₁ gametophytes of Gp. lemaneiformis in this study, and 9,989 SNPs located in nonrepetitive sequences were obtained. Among these SNPs, 92.02% of loci were identified as having a heterozygous genotype in at least one gametophyte, and 48.07% of loci had identical heterozygous genotypes in the 62 gametophytes. For each gametophyte, the proportions of homozygous and heterozygous loci ranged between 13.74 and 21.61% (mean of 17.04%) and between 66.36 and 83.59% (mean of 77.12%), respectively. The remainder were missing loci, representing an average 5.84%. Sources of heterozygous SNPs were free of exogenous DNA contamination, cross contamination among individuals, plastid and mitochondrial sequences, chimeras of different thallus parts or different cells, and repetitive sequences. Genotypes of heterozygous SNPs were verified by Sanger sequencing of PCR products and monoclones. Duplications of chromosomal rearrangements in the genome of Gp. lemaneiformis might explain the presence of heterozygous SNPs in haploid gametophytes.

PCR identification of genes of resistance to black rot in white cabbage using SSR-markers

The present study revealed the polymorphism of SSR loci by the resistance of white cabbage to black rot in contrasting isogenic samples of white cabbage. 2 informative SSR markers were selected: Ol10-C01 and Ol11-H06 for ranking breeding samples based on resistance to Xanthomonas campestris pv. campestris Dows. The microsatellite marker Ol10-C01 was tested on breeding samples and reveals polymorphism between them; therefore, it can be recommended for practical breeding for programs to develop black rot-resistant hybrids of cabbage, which will solve the problem of import substitution and healthy food (environmentally friendly products, grown without the use of means of chemical protection).

Comparative transcriptomics analysis uncovers alternative splicing events and molecular markers in cabbage (Brassica oleracea L.)

Alternative splicing (AS) events were identified and verified in cabbage by comparative transcriptome analysis. The corresponding markers were developed and the germplasm resources were identified. Alternative splicing (AS) is a central regulatory mechanism that greatly contributes to plant gene expression and transcriptome diversity. A large body of evidence has shown that AS complexity is relevant for plant development, evolution, complexity, and adaptation. Both insertion/deletion (InDel) and single nucleotide polymorphism (SNP) are typically co-dominant inheritance markers and have abundant polymorphisms. These have been widely used for marker-assisted selection, genetic mapping, and germplasm identification in plants. However, little is known about the molecular mechanisms underlying AS events and the development of markers including SNP and InDel from the cabbage transcriptome. In this study, three cabbage transcriptome datasets were collected and aligned to the cabbage reference genome to analyze AS events and marker development. 31,524 AS events were identified from three cabbage genotypes, accounting for 20.8% of the total cabbage genes. Alternative 3' splice site donor (A3SS) was the most frequent type of the four main AS events in cabbage. 70,475 InDels and 706,269 SNPs were identified with average frequencies of 1 InDel/6.9 kb and 1 SNP/0.7 kb, respectively. 71,942 potential SSRs were identified in 53,129 assembled unigenes with a density of 1 SSR/6.8 kb. The ratio of SNPs with synonymous/non-synonymous mutations was 1:0.65. 142 InDels and 36 SNPs were randomly selected and validated via Sanger sequencing and polymorphism was found among 66.2% of the InDels and 78.6% of the SNPs. Furthermore, 35 informative InDel markers were successfully used for genetic diversity analysis on 36 cabbage accessions. These results facilitate understanding of the molecular regulation mechanism underlying AS events in cabbage. They also provide molecular marker resource data for genetic mapping construction and germplasm identification, and facilitate the genetic improvement of cabbage via breeding.

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.

Dissecting quantitative resistance to Xanthomonas campestris pv. campestris in leaves of Brassica oleracea by QTL analysis

Black rot, caused by the bacterium Xanthomonas campestris pv. campestris (Xcc), produces important economic losses in crops of Brassica oleracea worldwide. Resistance to race 1, the most virulent and widespread in B. oleracea, is under quantitative control. Knowledge about the genetics of this resistance would help in designing strategies to control initial stages of invasion and development of the disease. QTL analysis of the resistance in the BolTBDH mapping population was performed. Resistance was measured with five traits related to initial stages of the invasion, success of infection and spread of the pathogen. Four single-trait QTLs of resistance were found, from which one represent novel variation. After performing multi-trait QTL, we concluded that spread of Xcc is related to the size of the leaf. Individuals from the mapping population follow two different strategies to cope with the spread of the disease: reducing lesion size or maintain more area of the leaf photosynthetically active, being more tolerant to Xcc invasion. Mechanisms underlying variation for resistance may be related to different aspects of plant immunity, including the synthesis of glucosinolates and phenolics.

Identification of NBS-encoding genes linked to black rot resistance in cabbage (Brassica oleracea var. capitata)

Heading cabbage is a nutritionally rich and economically important cruciferous vegetable. Black rot disease, caused by the bacterium Xanthomonas campestris pv. campestris, reduces both the yield and quality of the cabbage head. Nucleotide binding site (NBS)-encoding resistance (R) genes play a vital role in the plant immune response to various pathogens. In this study, we analyzed the expression and DNA sequence variation of 31 NBS-encoding genes in cabbage (Brassica oleracea var. capitata). These genes encoded TIR, NBS, LRR and RPW8 protein domains, all of which are known to be involved in disease resistance. RNA-seq revealed that these 31 genes were differentially expressed in leaf, root, silique, and stem tissues. Furthermore, qPCR analyses revealed that several of these genes were more highly expressed in resistant compared to susceptible cabbage lines, including Bol003711, Bol010135, Bol010559, Bol022784, Bol029866, Bol042121, Bol031422, Bol040045 and Bol042095. Further analysis of these genes promises to yield both practical benefits, such as molecular markers for marker-assisted breeding, and fundamental insights to the mechanisms of resistance to black rot in cabbage.

Incorporating pleiotropic quantitative trait loci in dissection of complex traits: seed yield in rapeseed as an example

A comprehensive linkage atlas for seed yield in rapeseed. Most agronomic traits of interest for crop improvement (including seed yield) are highly complex quantitative traits controlled by numerous genetic loci, which brings challenges for comprehensively capturing associated markers/genes. We propose that multiple trait interactions underlie complex traits such as seed yield, and that considering these component traits and their interactions can dissect individual quantitative trait loci (QTL) effects more effectively and improve yield predictions. Using a segregating rapeseed (Brassica napus) population, we analyzed a large set of trait data generated in 19 independent experiments to investigate correlations between seed yield and other complex traits, and further identified QTL in this population with a SNP-based genetic bin map. A total of 1904 consensus QTL accounting for 22 traits, including 80 QTL directly affecting seed yield, were anchored to the B. napus reference sequence. Through trait association analysis and QTL meta-analysis, we identified a total of 525 indivisible QTL that either directly or indirectly contributed to seed yield, of which 295 QTL were detected across multiple environments. A majority (81.5%) of the 525 QTL were pleiotropic. By considering associations between traits, we identified 25 yield-related QTL previously ignored due to contrasting genetic effects, as well as 31 QTL with minor complementary effects. Implementation of the 525 QTL in genomic prediction models improved seed yield prediction accuracy. Dissecting the genetic and phenotypic interrelationships underlying complex quantitative traits using this method will provide valuable insights for genomics-based crop improvement.

Fine Mapping and Transcriptome Analysis Reveal Candidate Genes Associated with Hybrid Lethality in Cabbage (Brassica Oleracea)

Hybrid lethality is a deleterious phenotype that is vital to species evolution. We previously reported hybrid lethality in cabbage (Brassica oleracea) and performed preliminary mapping of related genes. In the present study, the fine mapping of hybrid lethal genes revealed that BoHL1 was located on chromosome C1 between BoHLTO124 and BoHLTO130, with an interval of 101 kb. BoHL2 was confirmed to be between insertion-deletion (InDels) markers HL234 and HL235 on C4, with a marker interval of 70 kb. Twenty-eight and nine annotated genes were found within the two intervals of BoHL1 and BoHL2, respectively. We also applied RNA-Seq to analyze hybrid lethality in cabbage. In the region of BoHL1, seven differentially expressed genes (DEGs) and five resistance (R)-related genes (two in common, i.e., Bo1g153320 and Bo1g153380) were found, whereas in the region of BoHL2, two DEGs and four R-related genes (two in common, i.e., Bo4g173780 and Bo4g173810) were found. Along with studies in which R genes were frequently involved in hybrid lethality in other plants, these interesting R-DEGs may be good candidates associated with hybrid lethality. We also used SNP/InDel analyses and quantitative real-time PCR to confirm the results. This work provides new insight into the mechanisms of hybrid lethality in cabbage.

Genetics and fine mapping of a purple leaf gene, BoPr, in ornamental kale (Brassica oleracea L. var. acephala)

Background

Due to its variegated and colorful leaves, ornamental kale (Brassica oleracea L. var. acephala) has become a popular ornamental plant. In this study, we report the fine mapping and analysis of a candidate purple leaf gene using a backcross population and an F₂ population derived from two parental lines: W1827 (with white leaves) and P1835 (with purple leaves).

Results

Genetic analysis indicated that the purple leaf trait is controlled by a single dominant gene, which we named BoPr. Using markers developed based on the reference genome ‘02–12’, the BoPr gene was preliminarily mapped to a 280-kb interval of chromosome C09, with flanking markers M17 and BoID4714 at genetic distances of 4.3 cM and 1.5 cM, respectively. The recombination rate within this interval is almost 12 times higher than the usual level, which could be caused by assembly error for reference genome ‘02–12’ at this interval. Primers were designed based on ‘TO1000’, another B. oleracea reference genome. Among the newly designed InDel markers, BRID485 and BRID490 were found to be the closest to BoPr, flanking the gene at genetic distances of 0.1 cM and 0.2 cM, respectively; the interval between the two markers is 44.8 kb (reference genome ‘TO1000’). Seven annotated genes are located within the 44.8 kb genomic region, of which only Bo9g058630 shows high homology to AT5G42800 (dihydroflavonol reductase), which was identified as a candidate gene for BoPr. Blast analysis revealed that this 44.8 kb interval is located on an unanchored scaffold (Scaffold000035_P2) of ‘02–12’, confirming the existence of assembly error at the interval between M17 and BoID4714 for reference genome ‘02–12’.

Conclusions

This study identified a candidate gene for BoPr and lays a foundation for the cloning and functional analysis of this gene.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3613-x) contains supplementary material, which is available to authorized users.

GBS-based single dosage markers for linkage and QTL mapping allow gene mining for yield-related traits in sugarcane

BACKGROUND

Sugarcane (Saccharum spp.) is predominantly an autopolyploid plant with a variable ploidy level, frequent aneuploidy and a large genome that hampers investigation of its organization. Genetic architecture studies are important for identifying genomic regions associated with traits of interest. However, due to the genetic complexity of sugarcane, the practical applications of genomic tools have been notably delayed in this crop, in contrast to other crops that have already advanced to marker-assisted selection (MAS) and genomic selection. High-throughput next-generation sequencing (NGS) technologies have opened new opportunities for discovering molecular markers, especially single nucleotide polymorphisms (SNPs) and insertion-deletion (indels), at the genome-wide level. The objectives of this study were to (i) establish a pipeline for identifying variants from genotyping-by-sequencing (GBS) data in sugarcane, (ii) construct an integrated genetic map with GBS-based markers plus target region amplification polymorphisms and microsatellites, (iii) detect QTLs related to yield component traits, and (iv) perform annotation of the sequences that originated the associated markers with mapped QTLs to search putative candidate genes.

RESULTS

We used four pseudo-references to align the GBS reads. Depending on the reference, from 3,433 to 15,906 high-quality markers were discovered, and half of them segregated as single-dose markers (SDMs) on average. In addition to 7,049 non-redundant SDMs from GBS, 629 gel-based markers were used in a subsequent linkage analysis. Of 7,678 SDMs, 993 were mapped. These markers were distributed throughout 223 linkage groups, which were clustered in 18 homo(eo)logous groups (HGs), with a cumulative map length of 3,682.04 cM and an average marker density of 3.70 cM. We performed QTL mapping of four traits and found seven QTLs. Our results suggest the presence of a stable QTL across locations. Furthermore, QTLs to soluble solid content (BRIX) and fiber content (FIB) traits had markers linked to putative candidate genes.

CONCLUSIONS

This study is the first to report the use of GBS for large-scale variant discovery and genotyping of a mapping population in sugarcane, providing several insights regarding the use of NGS data in a polyploid, non-model species. The use of GBS generated a large number of markers and still enabled ploidy and allelic dosage estimation. Moreover, we were able to identify seven QTLs, two of which had great potential for validation and future use for molecular breeding in sugarcane.

Fine Mapping and Whole-Genome Resequencing Identify the Seed Coat Color Gene in Brassica rapa

A yellow seed coat is a desirable agronomic trait in the seeds of oilseed-type Brassica crops. In this study, we identified a candidate gene for seed coat color in Dahuang, a landrace of Brassica rapa. A previous study of Dahuang mapped the seed coat color gene Brsc1 to a 2.8-Mb interval on chromosome A9 of B. rapa. In the present study, the density of the linkage map for Brsc1 was increased by adding simple sequence repeat (SSR) markers, and the candidate region for Brsc1 was narrowed to 1.04 Mb. In addition, whole-genome resequencing with bulked segregant analysis (BSA) was conducted to identify candidate intervals for Brsc1. A genome-wide comparison of SNP profiles was performed between yellow-seeded and brown-seeded bulk samples. SNP index analyses identified a major candidate interval on chromosome A9 (A09:18,255,838-18,934,000, 678 kb) containing a long overlap with the target region recovered from the fine mapping results. According to gene annotation, Bra028067 (BrTT1) is an important candidate gene for Brsc1 in the overlapping region. Quantitative reverse transcription (qRT)-PCR revealed that BrTT1 mainly functions in the seed. Point mutations and small deletions in BrTT1 were found between yellow- and brown-seeded Dahuang plants. Collectively, the expression and sequence analysis results provide preliminary evidence that BrTT1 is a candidate gene for the seed coat color trait in Dahuang.

Transcriptomics-assisted quantitative trait locus fine mapping for the rapid identification of a nodulin 26-like intrinsic protein gene regulating boron efficiency in allotetraploid rapeseed

Allotetraploid rapeseed (Brassica napus L., An An Cn Cn , 2n = 4x = 38) is extraordinarily susceptible to boron (B) deficiency, a ubiquitous problem causing severe losses in seed yield. The breeding of B-efficient rapeseed germ plasm is a cost-effective and environmentally friendly strategy for the agricultural industry; however, genes regulating B efficiency in allotetraploid rapeseed have not yet been isolated. In this research, quantitative trait locus (QTL) fine mapping and digital gene expression (DGE) profiling were combined to identify the candidate genes underlying the major-effect QTL qBEC-A3a, which regulates B efficiency. Comparative phenotype analyses of the near-isogenic lines (NILs) indicated that qBEC-A3a plays a significant role in improving B efficiency under B deficiency. Exploiting QTL fine mapping and DGE analyses revealed a nodulin 26-like intrinsic protein (NIP) gene, which encodes a likely boric acid channel. The gene co-expression network for putative B transporters also highlighted its central role in the efficiency of B uptake. An integration of whole-genome re-sequencing (WGS) with bulked segregant analysis (BSA) authenticated the emerging availability of QTL-seq for the QTL analyses in allotetraploid rapeseed. Transcriptomics-assisted QTL mapping and comparative genomics provided novel insights into the rapid identification of quantitative trait genes (QTGs) in plant species with complex genomes.

Whole-Genome Mapping Reveals Novel QTL Clusters Associated with Main Agronomic Traits of Cabbage (Brassica oleracea var. capitata L.)

We describe a comprehensive quantitative trait locus (QTL) analysis for 24 main agronomic traits of cabbage. Field experiments were performed using a 196-line double haploid population in three seasons in 2011 and 2012 to evaluate important agronomic traits related to plant type, leaf, and head traits. In total, 144 QTLs with LOD threshold >3.0 were detected for the 24 agronomic traits: 25 for four plant-type-related traits, 64 for 10 leaf-related traits, and 55 for 10 head-related traits; each QTL explained 6.0-55.7% of phenotype variation. Of the QTLs, 95 had contribution rates higher than 10%, and 51 could be detected in more than one season. Major QTLs included Ph 3.1 (max R (2) = 55.7, max LOD = 28.2) for plant height, Ll 3.2 (max R (2) = 31.7, max LOD = 13.95) for leaf length, and Htd 3.2 (max R (2) = 28.5, max LOD = 9.49) for head transverse diameter; these could all be detected in more than one season. Twelve QTL clusters were detected on eight chromosomes, and the most significant four included Indel481-scaffold18376 (3.20 Mb), with five QTLs for five traits; Indel64-scaffold35418 (2.22 Mb), six QTLs for six traits; scaffold39782-Indel84 (1.78 Mb), 11 QTLs for 11 traits; and Indel353-Indel245 (9.89 Mb), seven QTLs for six traits. Besides, most traits clustered within the same region were significantly correlated with each other. The candidate genes at these regions were also discussed. Robust QTLs and their clusters obtained in this study should prove useful for marker-assisted selection (MAS) in cabbage breeding and in furthering our understanding of the genetic control of these traits.

Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L.)

Clubroot is a devastating disease caused by Plasmodiophora brassicae and results in severe losses of yield and quality in Brassica crops. Many clubroot resistance genes and markers are available in Brassica rapa but less is known in Brassica oleracea. Here, we applied the genotyping-by-sequencing (GBS) technique to construct a high-resolution genetic map and identify clubroot resistance (CR) genes. A total of 43,821 SNPs were identified using GBS data for two parental lines, one resistant and one susceptible lines to clubroot, and 18,187 of them showed >5× coverage in the GBS data. Among those, 4,103 were credibly genotyped for all 78 F₂ individual plants. These markers were clustered into nine linkage groups spanning 879.9 cM with an average interval of 1.15 cM. Quantitative trait loci (QTLs) survey based on three rounds of clubroot resistance tests using F_2:3 progenies revealed two and single major QTLs for Race 2 and Race 9 of P. brassicae, respectively. The QTLs show similar locations to the previously reported CR loci for Race 4 in B. oleracea but are in different positions from any of the CR loci found in B. rapa. We utilized two reference genome sequences in this study. The high-resolution genetic map developed herein allowed us to reposition 37 and 2 misanchored scaffolds in the 02–12 and TO1000DH genome sequences, respectively. Our data also support additional positioning of two unanchored 3.3 Mb scaffolds into the 02–12 genome sequence.

Next Generation Sequencing Technologies: The Doorway to the Unexplored Genomics of Non-Model Plants

Non-model plants i.e., the species which have one or all of the characters such as long life cycle, difficulty to grow in the laboratory or poor fecundity, have been schemed out of sequencing projects earlier, due to high running cost of Sanger sequencing. Consequently, the information about their genomics and key biological processes are inadequate. However, the advent of fast and cost effective next generation sequencing (NGS) platforms in the recent past has enabled the unearthing of certain characteristic gene structures unique to these species. It has also aided in gaining insight about mechanisms underlying processes of gene expression and secondary metabolism as well as facilitated development of genomic resources for diversity characterization, evolutionary analysis and marker assisted breeding even without prior availability of genomic sequence information. In this review we explore how different Next Gen Sequencing platforms, as well as recent advances in NGS based high throughput genotyping technologies are rewarding efforts on de-novo whole genome/transcriptome sequencing, development of genome wide sequence based markers resources for improvement of non-model crops that are less costly than phenotyping.