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

The power of model-to-crop translation illustrated by reducing seed loss from pod shatter in oilseed rape

Elucidation of key regulators in Arabidopsis fruit patterning has facilitated knowledge-translation into crop species to address yield loss caused by premature seed dispersal (pod shatter). In the 1980s, plant scientists descended on a small weed Arabidopsis thaliana (thale cress) and developed it into a powerful model system to study plant biology. The massive advances in genetics and genomics since then have allowed us to obtain incredibly detailed knowledge on specific biological processes of Arabidopsis growth and development, its genome sequence and the function of many of the individual genes. This wealth of information provides immense potential for translation into crops to improve their performance and address issues of global importance such as food security. Here, we describe how fundamental insight into the genetic mechanism by which seed dispersal occurs in members of the Brassicaceae family can be exploited to reduce seed loss in oilseed rape (Brassica napus). We demonstrate that by exploiting data on gene function in model species, it is possible to adjust the pod-opening process in oilseed rape, thereby significantly increasing yield. Specifically, we identified mutations in multiple paralogues of the INDEHISCENT and GA4 genes in B. napus and have overcome genetic redundancy by combining mutant alleles. Finally, we present novel software for the analysis of pod shatter data that is applicable to any crop for which seed dispersal is a serious problem. These findings highlight the tremendous potential of fundamental research in guiding strategies for crop improvement.

CRISPR/Cas9-Mediated Multiplex Genome Editing of JAGGED Gene in Brassica napus L

Pod shattering resistance is an essential component to achieving a high yield, which is a substantial objective in polyploid rapeseed cultivation. Previous studies have suggested that the Arabidopsis JAGGED (JAG) gene is a key factor implicated in the regulatory web of dehiscence fruit. However, its role in controlling pod shattering resistance in oilseed rape is still unknown. In this study, multiplex genome editing was carried out by the CRISPR/Cas9 system on five homoeologs (BnJAG.A02, BnJAG.C02, BnJAG.C06, BnJAG.A07, and BnJAG.A08) of the JAG gene. Knockout mutagenesis of all homoeologs drastically affected the development of the lateral organs in organizing pod shape and size. The cylindrical body of the pod comprised a number of undifferentiated cells like a callus, without distinctive valves, replum, septum, and valve margins. Pseudoseeds were produced, which were divided into two halves with an incomplete layer of cells (probably septum) that separated the undifferentiated cells. These mutants were not capable of generating any productive seeds for further generations. However, one mutant line was identified in which only a BnJAG.A08-NUB-Like paralog of the JAG gene was mutated. Knockout mutagenesis in BnJAG.A08-NUB gene caused significant changes in the pod dehiscence zone. The replum region of the mutant was increased to a great extent, resulting in enlarged cell size, bumpy fruit, and reduced length compared with the wild type. A higher replum-valve joint area may have increased the resistance to pod shattering by ~2-fold in JAG mutants compared with wild type. Our results offer a basis for understanding variations in Brassica napus fruit by mutating JAG genes and providing a way forward for other Brassicaceae species.

Identification of Markers Associated with Yield Traits and Morphological Features in Maize (Zea mays L.)

Association mapping is a powerful approach to detect associations between traits of interest and genetic markers based on linkage disequilibrium in molecular plant breeding. The aim of this study was the identification of single nucleotide polymorphisms (SNPs) and SilicoDArT markers associated with yield traits and morphological features in maize. Plant material constituted inbred lines. The field experiment with inbred lines was established on 10 m² plots in a set of complete random blocks in three replicates. We observed 22 quantitative traits. Association mapping was performed in this study using a method based on the mixed linear model with the population structure estimated by eigenanalysis (principal component analysis applied to all markers) and modeled by random effects. As a result of mapping, 969 markers (346 SNPs and 623 SilocoDArT) were selected from 49,911 identified polymorphic molecular markers, which were significantly associated with the analyzed morphological features and yield structure traits. Markers associated with five or six traits were selected during further analyses, including SilicoDArT 4591115 (anthocyanin coloration of anthers, length of main axis above the highest lateral branch, cob length, number of grains per cob, weight of fresh grains per cob and weight of fresh grains per cob at 15% moisture), SilicoDArT 7059939 (anthocyanin coloration of glumes of cob, time of anthesis—50% of flowering plants, time of silk emergence—50% of flowering plants, anthocyanin coloration of anthers and cob diameter), SilicoDArT 5587991 (anthocyanin coloration of glumes of cob, time of anthesis—50% of flowering plants, anthocyanin coloration of anthers, curvature of lateral branches and number of rows of grain). The two genetic similarity dendrograms between the inbred lines were constructed based on all significant SNPs and SilicoDArT markers. On both dendrograms lines clustered according to the kernel structure (flint, dent) and origin. The selected markers may be useful in predicting hybrid formulas in a heterosis culture. The present study demonstrated that molecular SNP and Silico DArT markers could be used in this species to group lines in terms of origin and lines with incomplete origin data. They can also be useful in maize in predicting the hybrid formula and can find applications in the selection of parental components for heterosis crossings.

Identification of SNP loci and candidate genes related to four important fatty acid composition in Brassica napus using genome wide association study

Rapeseed oil (canola, Brassica napus L.) is an important healthy vegetable oil throughout the world, the nutritional and economical value of which largely depends on its seed fatty acid composition. In this study, based on 201,187 SNP markers developed from the SLAF-seq (specific locus amplified fragment sequencing), a genome wide association study of four important fatty acid content traits (erucic acid, oleic acid, linoleic acid and linolenic acid) in a panel of 300 inbred lines of rapeseed in two environments (JXAU and JXRIS) was carried out. A total of 148 SNP loci significantly associated with these traits were detected by MLM model analysis respectively, and 30 SNP loci on A08 and C03 chromosomes were detected in three traits of erucic acid, oleic acid and linoleic acid contents simultaneously. Furthermore, 108 highly favorable alleles for increasing oleic acid and linoleic acid content, also for decreasing erucic acid content simultaneously were observed. By a basic local alignment search tool (BLAST) search with in a distance of 100 Kb around these significantly SNP-trait associations, we identified 20 orthologs of the functional candidate genes related to fatty acid biosynthesis, including the known vital fatty acid biosynthesis genes of BnaA.FAE1 and BnaC. FAE1 on the A08 and C03 chromosomes, and other potential candidate genes involving in the fatty acid biosynthesis pathway, such as the orthologs genes of FAD2, LACS09, KCS17, CER4, TT16 and ACBP5. This study lays a basis for uncovering the genetic variations and the improvement of fatty acid composition in B. napus.

Female choice for related males in wild red-backed toadlets (Pseudophryne coriacea)

Mate choice for genetic benefits is assumed to be widespread in nature, yet very few studies have comprehensively examined relationships between female mate choice and male genetic quality in wild populations. Here, we use exhaustive sampling and single nucleotide polymorphisms to provide a partial test of the “good genes as heterozygosity” hypothesis and the “genetic compatibility” hypothesis in an entire population of terrestrial breeding red-backed toadlets, Pseudophryne coriacea. We found that successful males did not display higher heterozygosity, despite a positive relationship between male heterozygosity and offspring heterozygosity. Rather, in the larger of 2 breeding events, we found that successful males were more genetically similar to their mate than expected under random mating, indicating that females can use pre- or post-copulatory mate choice mechanisms to bias paternity toward more related males. These findings provide no support for the good genes as heterozygosity hypothesis but lend support to the genetic compatibility hypothesis. A complete test of this hypothesis will now require evaluating how parental genetic similarity impacts offspring fitness. Terrestrial toadlets show a high degree of site fidelity, high levels of genetic structuring between populations, and frequently hybridize with sister species. As such, female mate choice for related males may be an adaptive strategy to reduce outbreeding depression. Our findings provide the first population-wide evidence for non-random preferential inbreeding in a wild amphibian. We argue that such reproductive patterns may be common in amphibians because extreme genetic differentiation within meta-populations creates an inherently high risk of outbreeding depression.

Stable Quantitative Resistance Loci to Blackleg Disease in Canola (Brassica napus L.) Over Continents

The hemibiotrophic fungus, Leptosphaeria maculans is the most devastating pathogen, causing blackleg disease in canola (Brassica napus L). To study the genomic regions involved in quantitative resistance (QR), 259-276 DH lines from Darmor-bzh/Yudal (DYDH) population were assessed for resistance to blackleg under shade house and field conditions across 3 years. In different experiments, the broad sense heritability varied from 43 to 95%. A total of 27 significant quantitative trait loci (QTL) for QR were detected on 12 chromosomes and explained between 2.14 and 10.13% of the genotypic variance. Of the significant QTL, at least seven were repeatedly detected across different experiments on chromosomes A02, A07, A09, A10, C01, and C09. Resistance alleles were mainly contributed by 'Darmor-bzh' but 'Yudal' also contributed few of them. Our results suggest that plant maturity and plant height may have a pleiotropic effect on QR in our conditions. We confirmed that Rlm9 which is present in 'Darmor-bzh' is not effective to confer resistance in our Australian field conditions. Comparative mapping showed that several R genes coding for nucleotide-binding leucine-rich repeat (LRR) receptors map in close proximity (within 200 Kb) of the significant trait-marker associations on the reference 'Darmor-bzh' genome assembly. More importantly, eight significant QTL regions were detected across diverse growing environments: Australia, France, and United Kingdom. These stable QTL identified herein can be utilized for enhancing QR in elite canola germplasm via marker- assisted or genomic selection strategies.

Molecular Evidence for Two Domestication Events in the Pea Crop

Pea, one of the founder crops from the Near East, has two wild species: Pisum sativum subsp. elatius, with a wide distribution centered in the Mediterranean, and P. fulvum, which is restricted to Syria, Lebanon, Israel, Palestine and Jordan. Using genome wide analysis of 11,343 polymorphic single nucleotide polymorphisms (SNPs) on a set of wild P. elatius (134) and P. fulvum (20) and 74 domesticated accessions (64 P. sativum landraces and 10 P. abyssinicum), we demonstrated that domesticated P. sativum and the Ethiopian pea (P. abyssinicum) were derived from different P. elatius genepools. Therefore, pea has at least two domestication events. The analysis does not support a hybrid origin of P. abyssinicum, which was likely introduced into Ethiopia and Yemen followed by eco-geographic adaptation. Both P. sativum and P. abyssinicum share traits that are typical of domestication, such as non-dormant seeds. Non-dormant seeds were also found in several wild P. elatius accessions which could be the result of crop to wild introgression or natural variation that may have been present during pea domestication. A sub-group of P. elatius overlaps with P. sativum landraces. This may be a consequence of bidirectional gene-flow or may suggest that this group of P. elatius is the closest extant wild relative of P. sativum.

Novel gene Sen2 conferring broad-spectrum resistance to Synchytrium endobioticum mapped to potato chromosome XI

Key message Sen2 gene for potato wart resistance, located on chromosome XI in a locus distinct from Sen1 , provides resistance against eight wart pathotypes, including the virulent ones important in Europe. Synchytrium endobioticum causes potato wart disease imposing severe losses in potato production, and as a quarantine pathogen in many countries, it results in lost trade markets and land for potato cultivation. The resistance to S. endobioticum pathotype 1(D1) is widespread in potato cultivars but new virulent pathotypes appear and the problem re-emerges. To characterize and map a new gene for resistance to potato wart, we used diploid F1 potato population from a cross of potato clone resistant to S. endobioticum pathotype 1(D1) and virulent pathotypes: 2(G1), 6(O1), 8(F1), 18(T1), 2(Ch1), 3(M1) and 39(P1) with a potato clone resistant to pathotype 1(D1) only. The 176 progeny clones were tested for resistance to eight wart pathotypes with a modified Glynne-Lemmerzahl method. Bimodal distributions and co-segregation of resistance in the population show that a single resistance gene, Sen2, underlies the resistance to eight pathotypes. Resistance to pathotype 1(D1) was additionally conferred by the locus Sen1 inherited from both parents. Sen2 was mapped to chromosome XI using DArTseq markers. The genetic and physical distances between Sen1 and Sen2 loci were indirectly estimated at 63 cM and 32 Mbp, respectively. We developed PCR markers co-segregating with the Sen2 locus that can be applied in marker-assisted selection of potatoes resistant to eight important pathotypes of S. endobioticum. Wide spectrum of the Sen2 resistance may be an indication of durability which can be enhanced by the pyramiding of the Sen2 and Sen1 loci as in 61 clones selected within this study.

The Site Frequency/Dosage Spectrum of Autopolyploid Populations

The Site Frequency Spectrum (SFS) and the heterozygosity of allelic variants are among the most important summary statistics for population genetic analysis of diploid organisms. We discuss the generalization of these statistics to populations of autopolyploid organisms in terms of the joint Site Frequency/Dosage Spectrum and its expected value for autopolyploid populations that follow the standard neutral model. Based on these results, we present estimators of nucleotide variability from High-Throughput Sequencing (HTS) data of autopolyploids and discuss potential issues related to sequencing errors and variant calling. We use these estimators to generalize Tajima's D and other SFS-based neutrality tests to HTS data from autopolyploid organisms. Finally, we discuss how these approaches fail when the number of individuals is small. In fact, in autopolyploids there are many possible deviations from the Hardy–Weinberg equilibrium, each reflected in a different shape of the individual dosage distribution. The SFS from small samples is often dominated by the shape of these deviations of the dosage distribution from its Hardy–Weinberg expectations.

Characterization of genetic diversity in Turkish common bean gene pool using phenotypic and whole-genome DArTseq-generated silicoDArT marker information

Turkey presents a great diversity of common bean landraces in farmers' fields. We collected 183 common bean accessions from 19 different Turkish geographic regions and 5 scarlet runner bean accessions to investigate their genetic diversity and population structure using phenotypic information (growth habit, and seed weight, flower color, bracteole shape and size, pod shape and leaf shape and color), geographic provenance and 12,557 silicoDArT markers. A total of 24.14% markers were found novel. For the entire population (188 accessions), the expected heterozygosity was 0.078 and overall gene diversity, Fst and Fis were 0.14, 0.55 and 1, respectively. Using marker information, model-based structure, principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA) algorithms clustered the 188 accessions into two main populations A (predominant) and B, and 5 unclassified genotypes, representing 3 meaningful heterotic groups for breeding purposes. Phenotypic information clearly distinguished these populations; population A and B, respectively, were bigger (>40g/100 seeds) and smaller (<40g/100 seeds) seed-sized. The unclassified population was pure and only contained climbing genotypes with 100 seed weight 2-3 times greater than populations A and B. Clustering was mainly based on A: seed weight, B: growth habit, C: geographical provinces and D: flower color. Mean kinship was generally low, but population B was more diverse than population A. Overall, a useful level of gene and genotypic diversity was observed in this work and can be used by the scientific community in breeding efforts to develop superior common bean strains.

Contrasting gene flow at different spatial scales revealed by genotyping-by-sequencing in Isocladus armatus, a massively colour polymorphic New Zealand marine isopod

Understanding how genetic diversity is maintained within populations is central to evolutionary biology. Research on colour polymorphism (CP), which typically has a genetic basis, can shed light on this issue. However, because gene flow can homogenise genetic variation, understanding population connectivity is critical in examining the maintenance of polymorphisms. In this study we assess the utility of genotyping-by-sequencing to resolve gene flow, and provide a preliminary investigation into the genetic basis of CP in Isocladus armatus, an endemic New Zealand marine isopod. Analysis of the genetic variation in 4,000 single nucleotide polymorphisms (SNPs) within and among populations and colour morphs revealed large differences in gene flow across two spatial scales. Marine isopods, which lack a pelagic larval phase, are typically assumed to exhibit greater population structuring than marine invertebrates possessing a biphasic life cycle. However, we found high gene flow rates and no genetic subdivision between two North Island populations situated 8 km apart. This suggests that I. armatus is capable of substantial dispersal along coastlines. In contrast, we identified a strong genetic disjunction between North and South Island populations. This result is similar to those reported in other New Zealand marine species, and is congruent with the presence of a geophysical barrier to dispersal down the east coast of New Zealand. We also found some support for a genetic basis to colouration evidenced by positive FST outlier tests, with two SNPs in particular showing strong association to the expression of a striped morph. Our study provides one of the first population genomic studies of a marine organism in New Zealand, and suggests that genotyping-by-sequencing can be a good alternative to more traditional investigations based on traditional markers such as microsatellites. Our study provides a foundation for further development of a highly tractable system for research on the evolutionary maintenance of CP.

High-Density Mapping of an Adult-Plant Stripe Rust Resistance Gene YrBai in Wheat Landrace Baidatou Using the Whole Genome DArTseq and SNP Analysis

Stripe rust, caused by the biotrophic fungus Puccinia striiformis f. sp. tritici (Pst), is one of the most widespread and destructive wheat diseases worldwide. Growing resistant cultivars is an effective approach for controlling this disease. However, because host resistance genes were easily overcome by new virulent Pst races, there is a continuous demand for identifying new effective wheat stripe rust resistance genes and develop closely linked markers for marker-assisted selection (MAS). Baidatou, an old Chinese wheat landrace, has been grown for several decades in Longnan region, Gansu Province, where stripe rust epidemics are frequent and severe. In our previous study, a single dominant gene YrBai in Baidatou was identified to control the adult-plant resistance (APR) to Chinese prevalent Pst race CYR33. And the gene was located on wheat chromosome 6DS by four polymorphic simple sequence repeat (SSR) and two sequence-related amplified polymorphism (SRAP) markers, with the genetic distances of two closely linked markers 3.6 and 5.4 cM, respectively. To further confirm the APR gene in Baidatou and construct the high-density map for the resistance gene, adult plants of F1, F2, F3, and F5:6 populations derived from the cross Mingxian169/Baidatou and two parents were inoculated with CYR33 at Yangling field, Shaanxi Province during 2014-2015, 2015-2016, and 2016-2017 crop seasons, respectively. The field evaluation results indicated that a single dominant gene confers the APR to Pst race CYR33 in Baidatou. 92 F3 lines and parents were sequenced using DArTseq technology based on wheat GBS1.0 platform, and 31 genetic maps consisted of 2,131 polymorphic SilicoDArT and 952 SNP markers spanning 4,293.94 cM were constructed. Using polymorphic SilicoDArT, SNP markers and infection types (ITs) data of F3 lines, the gene YrBai was further located in 0.8 cM region on wheat chromosome 6D. These closely linked markers developed in this study should be useful for MAS for Baidatou in crop improvement and map-based clone this gene.

Genome-Wide Association Study of Cadmium Accumulation at the Seedling Stage in Rapeseed (Brassica napus L.)

Cadmium is a potentially toxic heavy metal to human health. Rapeseed (Brassica napus L.), a vegetable and oilseed crop, might also be a Cd hyperaccumulator, but there is little information on this trait in rapeseed. We evaluated Cd accumulation in different oilseed accessions and employed a genome-wide association study to identify quantitative trait loci (QTLs) related to Cd accumulation. A total of 419 B. napus accessions and inbred lines were genotyped with a 60K Illumina Infinium SNP array of Brassica. Wide genotypic variations in Cd concentration and translocation were found. Twenty-five QTLs integrated with 98 single-nucleotide polymorphisms (SNPs) located at 15 chromosomes were associated with Cd accumulation traits. These QTLs explained 3.49-7.57% of the phenotypic variation observed. Thirty-two candidate genes were identified in these genomic regions, and they were 0.33-497.97 kb away from the SNPs. We found orthologs of Arabidopsis thaliana located near the significant SNPs on the B. napus genome, including NRAMP6 (natural resistance-associated macrophage protein 6), IRT1 (iron-regulated transporter 1), CAD1 (cadmium-sensitive 1), and PCS2 (phytochelatin synthase 2). Of them, four candidate genes were verified by qRT-PCR, the expression levels of which were significantly higher after exposure to Cd than in the controls. Our results might facilitate the study of the genetic basis of Cd accumulation and the cloning of candidate Cd accumulation genes, which could be used to help reduce Cd levels in edible plant parts and/or create more efficient hyperaccumulators.

The effect of INDEHISCENT point mutations on silique shatter resistance in oilseed rape (Brassica napus)

This study elucidates the influence of indehiscent mutations on rapeseed silique shatter resistance. A phenotype with enlarged replum-valve joint area and altered cell dimensions in the dehiscence zone is described. Silique shattering is a major factor reducing the yield stability of oilseed rape (Brassica napus). Attempts to improve shatter resistance often include the use of mutations in target genes identified from Arabidopsis (Arabidopsis thaliana). A variety of phenotyping methods assessing the level of shatter resistance were previously described. However, a comparative and comprehensive evaluation of the methods has not yet been undertaken. We verified the increase of shatter resistance in indehiscent double knock-down mutants obtained by TILLING with a systematic approach comparing three independent phenotyping methods. A positive correlation of silique length and shatter resistance was observed and accounted for in the analyses. Microscopic studies ruled out the influence of different lignification patterns. Instead, we propose a model to explain increased shattering resistance of indehiscent rapeseed mutants by altered cell shapes and sizes within the contact surfaces of replum and valves.

Populations of doubled haploids for genetic mapping in hexaploid winter triticale

To create a framework for genetic dissection of hexaploid triticale, six populations of doubled haploid (DH) lines were developed from pairwise hybrids of high-yielding winter triticale cultivars. The six populations comprise between 97 and 231 genotyped DH lines each, totaling 957 DH lines. A consensus genetic map spans 4593.9 cM is composed of 1576 unique DArT markers. The maps reveal several structural rearrangements in triticale genomes. In preliminary tests of the populations and maps, markers specific to wheat segments of the engineered rye chromosome 1R (RM1B) were identified. Example QTL mapping of days to heading in cv. Krakowiak revealed loci on chromosomes 2BL and 2R responsible for extended vernalization requirement, and candidate genes were identified. The material is available to all parties interested in triticale genetics.

VARIETAL IDENTIFICATION IN HOUSEHOLD SURVEYS: RESULTS FROM THREE HOUSEHOLD-BASED METHODS AGAINST THE BENCHMARK OF DNA FINGERPRINTING IN SOUTHERN ETHIOPIA

Accurate crop varietal identification is the backbone of any high-quality assessment of outcomes and impacts. Sweetpotato (Ipomoea batatas) varieties have important nutritional differences, and there is a strong interest to identify nutritionally superior varieties for dissemination. In agricultural household surveys, such information is often collected based on the farmer's self-report. In this article, we present the results of a data capture experiment on sweet potato varietal identification in southern Ethiopia. Three household-based methods of identifying varietal adoption are tested against the benchmark of DNA fingerprinting: (A) Elicitation from farmers with basic questions for the most widely planted variety; (B) Farmer elicitation on five sweet potato phenotypic attributes by showing a visual-aid protocol; and (C) Enumerator recording observations on five sweet potato phenotypic attributes using a visual-aid protocol and visiting the field. In total, 20% of farmers identified a variety as improved when in fact it was local and 19% identified a variety as local when it was in fact improved. The variety names given by farmers delivered inconsistent and inaccurate varietal identities. Visual-aid protocols employed in methods B and C were better than those in method A, but greatly underestimated the adoption estimates given by the DNA fingerprinting method. Our results suggest that estimating the adoption of improved varieties with methods based on farmer self-reports is questionable and point towards a wider use of DNA fingerprinting in adoption and impact assessments.

Fresh is best: Accurate SNP genotyping from koala scats

Maintaining genetic diversity is a crucial component in conserving threatened species. For the iconic Australian koala, there is little genetic information on wild populations that is not either skewed by biased sampling methods (e.g., sampling effort skewed toward urban areas) or of limited usefulness due to low numbers of microsatellites used. The ability to genotype DNA extracted from koala scats using next‐generation sequencing technology will not only help resolve location sample bias but also improve the accuracy and scope of genetic analyses (e.g., neutral vs. adaptive genetic diversity, inbreeding, and effective population size). Here, we present the successful SNP genotyping (1272 SNP loci) of koala DNA extracted from scat, using a proprietary DArTseq^™ protocol. We compare genotype results from two‐day‐old scat DNA and 14‐day‐old scat DNA to a blood DNA template, to test accuracy of scat genotyping. We find that DNA from fresher scat results in fewer loci with missing information than DNA from older scat; however, 14‐day‐old scat can still provide useful genetic information, depending on the research question. We also find that a subset of 209 conserved loci can accurately identify individual koalas, even from older scat samples. In addition, we find that DNA sequences identified from scat samples through the DArTseq^™ process can provide genetic identification of koala diet species, bacterial and viral pathogens, and parasitic organisms.

Genomic prediction using DArT-Seq technology for yellowtail kingfish Seriola lalandi

Background

Genomic prediction using Diversity Arrays Technology (DArT) genotype by sequencing platform has not been reported in yellowtail kingfish (Seriola lalandi). The principal aim of this study was to address this knowledge gap and to assess predictive ability of genomic Best Linear Unbiased Prediction (gBLUP) for traits of commercial importance in a yellowtail kingfish population comprising 752 individuals that had DNA sequence and phenotypic records for growth traits (body weight, fork length and condition index). The gBLUP method was used due to its computational efficiency and it showed similar predictive performance to other approaches, especially for traits whose variation is of polygenic nature, such as body traits analysed in this study. The accuracy or predictive ability of the gBLUP model was estimated for three growth traits: body weight, folk length and condition index.

Results

The prediction accuracy was moderate to high (0.44 to 0.69) for growth-related traits. The predictive ability for body weight increased by 17.0% (from 0.69 to 0.83) when missing genotype was imputed. Within population prediction using five-fold across validation approach showed that the gBLUP model performed well for growth traits (weight, length and condition factor), with the coefficient of determination (R²) from linear regression analysis ranging from 0.49 to 0.71.

Conclusions

Collectively our results demonstrated, for the first time in yellowtail kingfish, the potential application of genomic selection for growth-related traits in the future breeding program for this species, S. lalandi.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4493-4) contains supplementary material, which is available to authorized users.

Speed breeding is a powerful tool to accelerate crop research and breeding

The growing human population and a changing environment have raised significant concern for global food security, with the current improvement rate of several important crops inadequate to meet future demand 1 . This slow improvement rate is attributed partly to the long generation times of crop plants. Here, we present a method called 'speed breeding', which greatly shortens generation time and accelerates breeding and research programmes. Speed breeding can be used to achieve up to 6 generations per year for spring wheat (Triticum aestivum), durum wheat (T. durum), barley (Hordeum vulgare), chickpea (Cicer arietinum) and pea (Pisum sativum), and 4 generations for canola (Brassica napus), instead of 2-3 under normal glasshouse conditions. We demonstrate that speed breeding in fully enclosed, controlled-environment growth chambers can accelerate plant development for research purposes, including phenotyping of adult plant traits, mutant studies and transformation. The use of supplemental lighting in a glasshouse environment allows rapid generation cycling through single seed descent (SSD) and potential for adaptation to larger-scale crop improvement programs. Cost saving through light-emitting diode (LED) supplemental lighting is also outlined. We envisage great potential for integrating speed breeding with other modern crop breeding technologies, including high-throughput genotyping, genome editing and genomic selection, accelerating the rate of crop improvement.

A High-Density Integrated DArTseq SNP-Based Genetic Map of Pisum fulvum and Identification of QTLs Controlling Rust Resistance

Pisum fulvum, a wild relative of pea is an important source of allelic diversity to improve the genetic resistance of cultivated species against fungal diseases of economic importance like the pea rust caused by Uromyces pisi. To unravel the genetic control underlying resistance to this fungal disease, a recombinant inbred line (RIL) population was generated from a cross between two P. fulvum accessions, IFPI3260 and IFPI3251, and genotyped using Diversity Arrays Technology. A total of 9,569 high-quality DArT-Seq and 8,514 SNPs markers were generated. Finally, a total of 12,058 markers were assembled into seven linkage groups, equivalent to the number of haploid chromosomes of P. fulvum and P. sativum. The newly constructed integrated genetic linkage map of P. fulvum covered an accumulated distance of 1,877.45 cM, an average density of 1.19 markers cM^-1 and an average distance between adjacent markers of 1.85 cM. The composite interval mapping revealed three QTLs distributed over two linkage groups that were associated with the percentage of rust disease severity (DS%). QTLs UpDSII and UpDSIV were located in the LGs II and IV respectively and were consistently identified both in adult plants over 3 years at the field (Córdoba, Spain) and in seedling plants under controlled conditions. Whenever they were detected, their contribution to the total phenotypic variance varied between 19.8 and 29.2. A third QTL (UpDSIV.2) was also located in the LGIVand was environmentally specific as was only detected for DS % in seedlings under controlled conditions. It accounted more than 14% of the phenotypic variation studied. Taking together the data obtained in the study, it could be concluded that the expression of resistance to fungal diseases in P. fulvum originates from the resistant parent IFPI3260.

Genome of wild olive and the evolution of oil biosynthesis

Here we present the genome sequence and annotation of the wild olive tree (Olea europaea var. sylvestris), called oleaster, which is considered an ancestor of cultivated olive trees. More than 50,000 protein-coding genes were predicted, a majority of which could be anchored to 23 pseudochromosomes obtained through a newly constructed genetic map. The oleaster genome contains signatures of two Oleaceae lineage-specific paleopolyploidy events, dated at ∼28 and ∼59 Mya. These events contributed to the expansion and neofunctionalization of genes and gene families that play important roles in oil biosynthesis. The functional divergence of oil biosynthesis pathway genes, such as FAD2, SACPD, EAR, and ACPTE, following duplication, has been responsible for the differential accumulation of oleic and linoleic acids produced in olive compared with sesame, a closely related oil crop. Duplicated oleaster FAD2 genes are regulated by an siRNA derived from a transposable element-rich region, leading to suppressed levels of FAD2 gene expression. Additionally, neofunctionalization of members of the SACPD gene family has led to increased expression of SACPD2, 3, 5, and 7, consequently resulting in an increased desaturation of steric acid. Taken together, decreased FAD2 expression and increased SACPD expression likely explain the accumulation of exceptionally high levels of oleic acid in olive. The oleaster genome thus provides important insights into the evolution of oil biosynthesis and will be a valuable resource for oil crop genomics.

Why does inbreeding reduce male paternity? Effects on sexually selected traits

Mating with relatives has often been shown to negatively affect offspring fitness (inbreeding depression). There is considerable evidence for inbreeding depression due to effects on naturally selected traits, particularly those expressed early in life, but there is less evidence of it for sexually selected traits. This is surprising because sexually selected traits are expected to exhibit strong inbreeding depression. Here, we experimentally created inbred and outbred male mosquitofish (Gambusia holbrooki). Inbred males were the offspring of matings between full siblings. We then investigated how inbreeding influenced a number of sexually selected male traits, specifically: attractiveness, sperm number and velocity, as well as sperm competitiveness based on a male's share of paternity. We found no inbreeding depression for male attractiveness or sperm traits. There was, however, evidence that lower heterozygosity decreased paternity due to reduced sperm competitiveness. Our results add to the growing evidence that competitive interactions exacerbate the negative effects of the increased homozygosity that arises when there is inbreeding.

A High Density Genetic Map Derived from RAD Sequencing and Its Application in QTL Analysis of Yield-Related Traits in Vigna unguiculata

Cowpea [Vigna unguiculata (L.) Walp.] is an annual legume of economic importance and widely grown in the semi-arid tropics. However, high-density genetic maps of cowpea are still lacking. Here, we identified 34,868 SNPs (single nucleotide polymorphisms) that were distributed in the cowpea genome based on the RAD sequencing (restriction-site associated DNA sequencing) technique using a population of 170 individuals (two cowpea parents and 168 F_2:3 progenies). Of these, 17,996 reliable SNPs were allotted to 11 consensus linkage groups (LGs). The length of the genetic map was 1,194.25 cM in total with a mean distance of 0.066 cM/SNP marker locus. Using this map and the F_2:3 population, combined with the CIM (composite interval mapping) method, eleven quantitative trait loci (QTL) of yield-related trait were detected on seven LGs (LG4, 5, 6, 7, 9, 10, and 11) in cowpea. These QTL explained 0.05-17.32% of the total phenotypic variation. Among these, four QTL were for pod length, four QTL for thousand-grain weight (TGW), two QTL for grain number per pod, and one QTL for carpopodium length. Our results will provide a foundation for understanding genes related to grain yield in the cowpea and genus Vigna.

Genetics of dispersal

Dispersal is a process of central importance for the ecological and evolutionary dynamics of populations and communities, because of its diverse consequences for gene flow and demography. It is subject to evolutionary change, which begs the question, what is the genetic basis of this potentially complex trait? To address this question, we (i) review the empirical literature on the genetic basis of dispersal, (ii) explore how theoretical investigations of the evolution of dispersal have represented the genetics of dispersal, and (iii) discuss how the genetic basis of dispersal influences theoretical predictions of the evolution of dispersal and potential consequences.

Dispersal has a detectable genetic basis in many organisms, from bacteria to plants and animals. Generally, there is evidence for significant genetic variation for dispersal or dispersal‐related phenotypes or evidence for the micro‐evolution of dispersal in natural populations. Dispersal is typically the outcome of several interacting traits, and this complexity is reflected in its genetic architecture: while some genes of moderate to large effect can influence certain aspects of dispersal, dispersal traits are typically polygenic. Correlations among dispersal traits as well as between dispersal traits and other traits under selection are common, and the genetic basis of dispersal can be highly environment‐dependent.

By contrast, models have historically considered a highly simplified genetic architecture of dispersal. It is only recently that models have started to consider multiple loci influencing dispersal, as well as non‐additive effects such as dominance and epistasis, showing that the genetic basis of dispersal can influence evolutionary rates and outcomes, especially under non‐equilibrium conditions. For example, the number of loci controlling dispersal can influence projected rates of dispersal evolution during range shifts and corresponding demographic impacts. Incorporating more realism in the genetic architecture of dispersal is thus necessary to enable models to move beyond the purely theoretical towards making more useful predictions of evolutionary and ecological dynamics under current and future environmental conditions. To inform these advances, empirical studies need to answer outstanding questions concerning whether specific genes underlie dispersal variation, the genetic architecture of context‐dependent dispersal phenotypes and behaviours, and correlations among dispersal and other traits.

Assessment of Genetic Diversity and Structure of Large Garlic (Allium sativum) Germplasm Bank, by Diversity Arrays Technology "Genotyping-by-Sequencing" Platform (DArTseq)

Garlic (Allium sativum) is used worldwide in cooking and industry, including pharmacology/medicine and cosmetics, for its interesting properties. Identifying redundancies in germplasm blanks to generate core collections is a major concern, mostly in large stocks, in order to reduce space and maintenance costs. Yet, similar appearance and phenotypic plasticity of garlic varieties hinder their morphological classification. Molecular studies are challenging, due to the large and expected complex genome of this species, with asexual reproduction. Classical molecular markers, like isozymes, RAPD, SSR, or AFLP, are not convenient to generate germplasm core-collections for this species. The recent emergence of high-throughput genotyping-by-sequencing (GBS) approaches, like DArTseq, allow to overcome such limitations to characterize and protect genetic diversity. Therefore, such technology was used in this work to: (i) assess genetic diversity and structure of a large garlic-germplasm bank (417 accessions); (ii) create a core collection; (iii) relate genotype to agronomical features; and (iv) describe a cost-effective method to manage genetic diversity in garlic-germplasm banks. Hierarchical-cluster analysis, principal-coordinates analysis and STRUCTURE showed general consistency, generating three main garlic-groups, mostly determined by variety and geographical origin. In addition, high-resolution genotyping identified 286 unique and 131 redundant accessions, used to select a reduced size germplasm-bank core collection. This demonstrates that DArTseq is a cost-effective method to analyze species with large and expected complex genomes, like garlic. To the best of our knowledge, this is the first report of high-throughput genotyping of a large garlic germplasm. This is particularly interesting for garlic adaptation and improvement, to fight biotic and abiotic stresses, in the current context of climate change and global warming.

Physical Map of the Short Arm of Bread Wheat Chromosome 3D

Bread wheat ( L.) is one of the most important crops worldwide. Although a reference genome sequence would represent a valuable resource for wheat improvement through genomics-assisted breeding and gene cloning, its generation has long been hampered by its allohexaploidy, high repeat content, and large size. As a part of a project coordinated by the International Wheat Genome Sequencing Consortium (IWGSC), a physical map of the short arm of wheat chromosome 3D (3DS) was prepared to facilitate reference genome assembly and positional gene cloning. It comprises 869 contigs with a cumulative length of 274.5 Mbp and represents 85.5% of the estimated chromosome arm size. Eighty-six Mbp of survey sequences from chromosome arm 3DS were assigned in silico to physical map contigs via next-generation sequencing of bacterial artificial chromosome pools, thus providing a high-density framework for physical map ordering along the chromosome arm. About 60% of the physical map was anchored in this single experiment. Finally, 1393 high-confidence genes were anchored to the physical map. Comparisons of gene space of the chromosome arm 3DS with genomes of closely related species [ (L.) P.Beauv., rice ( L.), and sorghum [ (L.) Moench] and homeologous wheat chromosomes provided information about gene movement on the chromosome arm.

In-depth genome characterization of a Brazilian common bean core collection using DArTseq high-density SNP genotyping

BACKGROUND

Common bean is a legume of social and nutritional importance as a food crop, cultivated worldwide especially in developing countries, accounting for an important source of income for small farmers. The availability of the complete sequences of the two common bean genomes has dramatically accelerated and has enabled new experimental strategies to be applied for genetic research. DArTseq has been widely used as a method of SNP genotyping allowing comprehensive genome coverage with genetic applications in common bean breeding programs.

RESULTS

Using this technology, 6286 SNPs (1 SNP/86.5 Kbp) were genotyped in genic (43.3%) and non-genic regions (56.7%). Genetic subdivision associated to the common bean gene pools (K = 2) and related to grain types (K = 3 and K = 5) were reported. A total of 83% and 91% of all SNPs were polymorphic within the Andean and Mesoamerican gene pools, respectively, and 26% were able to differentiate the gene pools. Genetic diversity analysis revealed an average H E of 0.442 for the whole collection, 0.102 for Andean and 0.168 for Mesoamerican gene pools (F ST  = 0.747 between gene pools), 0.440 for the group of cultivars and lines, and 0.448 for the group of landrace accessions (F ST  = 0.002 between cultivar/line and landrace groups). The SNP effects were predicted with predominance of impact on non-coding regions (77.8%). SNPs under selection were identified within gene pools comparing landrace and cultivar/line germplasm groups (Andean: 18; Mesoamerican: 69) and between the gene pools (59 SNPs), predominantly on chromosomes 1 and 9. The LD extension estimate corrected for population structure and relatedness (r2SV) was ~ 88 kbp, while for the Andean gene pool was ~ 395 kbp, and for the Mesoamerican was ~ 130 kbp.

CONCLUSIONS

For common bean, DArTseq provides an efficient and cost-effective strategy of generating SNPs for large-scale genome-wide studies. The DArTseq resulted in an operational panel of 560 polymorphic SNPs in linkage equilibrium, providing high genome coverage. This SNP set could be used in genotyping platforms with many applications, such as population genetics, phylogeny relation between common bean varieties and support to molecular breeding approaches.

Genome-Wide SNP Markers Based on SLAF-Seq Uncover Breeding Traces in Rapeseed (Brassica napus L.)

Single Nucleotide Polymorphisms (SNPs) are the most abundant and richest form of genomic polymorphism, and hence make highly favorable markers for genetic map construction and genome-wide association studies. In this study, a total of 300 rapeseed accessions (278 representative of Chinese germplasm, plus 22 outgroup accessions of different origins and ecotypes) were collected and sequenced using Specific-Locus Amplified Fragment Sequencing (SLAF-seq) technology, obtaining 660.25M reads with an average sequencing depth of 6.27 × and a mean Q30 of 85.96%. Based on the 238,711 polymorphic SLAF tags a total of 1,197,282 SNPs were discovered, and a subset of 201,817 SNPs with minor allele frequency >0.05 and integrity >0.8 were selected. Of these, 30,877 were designated SNP "hotspots," and 41 SNP-rich genomic regions could be delineated, with 100 genes associated with plant resistance, vernalization response, and signal transduction detected in these regions. Subsequent analysis of genetic diversity, linkage disequilibrium (LD), and population structure in the 300 accessions was carried out based on the 201,817 SNPs. Nine subpopulations were observed based on the population structure analysis. Hierarchical clustering and principal component analysis divided the 300 varieties roughly in accordance with their ecotype origins. However, spring-type varieties were intermingled with semi-winter type varieties, indicating frequent hybridization between spring and semi-winter ecotypes in China. In addition, LD decay across the whole genome averaged 299 kb when r² = 0.1, but the LD decay in the A genome (43 kb) was much shorter than in the C genome (1,455 kb), supporting the targeted introgression of the A genome from progenitor species B. rapa into Chinese rapeseed. This study also lays the foundation for genetic analysis of important agronomic traits using this rapeseed population.

Investigation of the Genetic Diversity and Quantitative Trait Loci Accounting for Important Agronomic and Seed Quality Traits in Brassica carinata

Brassica carinata (BBCC) is an allotetraploid in Brassicas with unique alleles for agronomic traits and has huge potential as source for biodiesel production. To investigate the genome-wide molecular diversity, population structure and linkage disequilibrium (LD) pattern in this species, we genotyped a panel of 81 accessions of B. carinata with genotyping by sequencing approach DArTseq, generating a total of 54,510 polymorphic markers. Two subpopulations were exhibited in the B. carinata accessions. The average distance of LD decay (r2 = 0.1) in B subgenome (0.25 Mb) was shorter than that of C subgenome (0.40 Mb). Genome-wide association analysis (GWAS) identified a total of seven markers significantly associated with five seed quality traits in two experiments. To further identify the quantitative trait loci (QTL) for important agronomic and seed quality traits, we phenotyped a doubled haploid (DH) mapping population derived from the "YW" cross between two parents (Y-BcDH64 and W-BcDH76) representing from the two subpopulations. The YW DH population and its parents were grown in three contrasting environments; spring (Hezheng and Xining, China), semi-winter (Wuhan, China), and spring (Wagga Wagga, Australia) across 5 years for QTL mapping. Genetic bases of phenotypic variation in seed yield and its seven related traits, and six seed quality traits were determined. A total of 282 consensus QTL accounting for these traits were identified including nine major QTL for flowering time, oleic acid, linolenic acid, pod number of main inflorescence, and seed weight. Of these, 109 and 134 QTL were specific to spring and semi-winter environment, respectively, while 39 consensus QTL were identified in both contrasting environments. Two QTL identified for linolenic acid (B3) and erucic acid (C7) were validated in the diverse lines used for GWAS. A total of 25 QTL accounting for flowering time, erucic acid, and oleic acid were aligned to the homologous QTL or candidate gene regions in the C genome of B. napus. These results would not only provide insights for genetic improvement of this species, but will also identify useful genetic variation hidden in the Cc subgenome of B. carinata to improve canola cultivars.

Identification and Mapping of Adult Plant Resistance Loci to Leaf Rust and Stripe Rust in Common Wheat Cultivar Kundan

Leaf rust (LR) and stripe rust (YR) are important diseases of wheat worldwide. We used 148 recombinant inbred lines (RIL) from the cross of Avocet × Kundan for determining and mapping the genetic basis of adult plant resistance (APR) loci. The population was phenotyped LR and YR for three seasons in field trials conducted in Mexico and genotyped with the diversity arrays technology sequencing (DArT-Seq) and simple sequence repeat markers. The final genetic map was constructed using 2,937 polymorphic markers with an average distance of 1.29 centimorgans between markers. Inclusive composite interval mapping identified two co-located APR quantitative trait loci (QTL) for LR and YR, two LR QTL, and three YR QTL. The co-located resistance QTL on chromosome 1BL corresponded to the pleiotropic APR gene Lr46/Yr29. QLr.cim-2BL, QYr.cim-2AL, and QYr.cim-5AS could be identified as new resistance loci in this population. Lr46/Yr29 contributed 49.5 to 65.1 and 49.2 to 66.1% of LR and YR variations, respectively. The additive interaction between detected QTL showed that LR severities for RIL combining four QTL ranged between 5.3 and 25.8%, whereas the lowest YR severities were for RIL carrying QTL on chromosomes 1BL + 2AL + 6AL. The high-density DArT-Seq markers across chromosomes can be used in fine mapping of the targeted loci and development SNP markers.

The in silico identification and characterization of a bread wheat/Triticum militinae introgression line

The capacity of the bread wheat (Triticum aestivum) genome to tolerate introgression from related genomes can be exploited for wheat improvement. A resistance to powdery mildew expressed by a derivative of the cross-bread wheat cv. Tähti × T. militinae (Tm) is known to be due to the incorporation of a Tm segment into the long arm of chromosome 4A. Here, a newly developed in silico method termed rearrangement identification and characterization (RICh) has been applied to characterize the introgression. A virtual gene order, assembled using the GenomeZipper approach, was obtained for the native copy of chromosome 4A; it incorporated 570 4A DArTseq markers to produce a zipper comprising 2132 loci. A comparison between the native and introgressed forms of the 4AL chromosome arm showed that the introgressed region is located at the distal part of the arm. The Tm segment, derived from chromosome 7G, harbours 131 homoeologs of the 357 genes present on the corresponding region of Chinese Spring 4AL. The estimated number of Tm genes transferred along with the disease resistance gene was 169. Characterizing the introgression's position, gene content and internal gene order should not only facilitate gene isolation, but may also be informative with respect to chromatin structure and behaviour studies.

Genotyping-by-Sequencing and Its Application to Oat Genomic Research

Genotyping-by-sequencing (GBS) has emerged as a useful genomic approach for sampling genome-wide genetic variation, performing genome-wide association mapping, and conducting genomic selection. It is a combined one-step process of SNP marker discovery and genotyping through genome reduction with restriction enzymes and SNP calling with or without a sequenced genome. This approach has the advantage of being rapid, high throughput, cost effective, and applicable to organisms without sequenced genomes. It has been increasingly applied to generate SNP genotype data for plant genetic and genomic studies. To facilitate a wider GBS application, particularly in oat genetic and genomic research, we describe the GBS approach, review the current applications of GBS in plant species, and highlight some applications of GBS to oat research. We also discuss issues in various applications of GBS and provide some perspectives in GBS research. Recent developments of bioinformatics pipelines in high-quality SNP discovery for polyploid crops will enhance the application of GBS to oat genetic and genomic research.

Genome-Wide SNP Markers Based on SLAF-Seq Uncover Breeding Traces in Rapeseed (Brassica napus L.)

Single Nucleotide Polymorphisms (SNPs) are the most abundant and richest form of genomic polymorphism, and hence make highly favorable markers for genetic map construction and genome-wide association studies. In this study, a total of 300 rapeseed accessions (278 representative of Chinese germplasm, plus 22 outgroup accessions of different origins and ecotypes) were collected and sequenced using Specific-Locus Amplified Fragment Sequencing (SLAF-seq) technology, obtaining 660.25M reads with an average sequencing depth of 6.27 × and a mean Q30 of 85.96%. Based on the 238,711 polymorphic SLAF tags a total of 1,197,282 SNPs were discovered, and a subset of 201,817 SNPs with minor allele frequency >0.05 and integrity >0.8 were selected. Of these, 30,877 were designated SNP "hotspots," and 41 SNP-rich genomic regions could be delineated, with 100 genes associated with plant resistance, vernalization response, and signal transduction detected in these regions. Subsequent analysis of genetic diversity, linkage disequilibrium (LD), and population structure in the 300 accessions was carried out based on the 201,817 SNPs. Nine subpopulations were observed based on the population structure analysis. Hierarchical clustering and principal component analysis divided the 300 varieties roughly in accordance with their ecotype origins. However, spring-type varieties were intermingled with semi-winter type varieties, indicating frequent hybridization between spring and semi-winter ecotypes in China. In addition, LD decay across the whole genome averaged 299 kb when r² = 0.1, but the LD decay in the A genome (43 kb) was much shorter than in the C genome (1,455 kb), supporting the targeted introgression of the A genome from progenitor species B. rapa into Chinese rapeseed. This study also lays the foundation for genetic analysis of important agronomic traits using this rapeseed population.

A Major Locus for Manganese Tolerance Maps on Chromosome A09 in a Doubled Haploid Population of Brassica napus L

Soil acidity poses a major threat to productivity of several crops; mainly due to the prevalence of toxic levels of Al³⁺ and Mn²⁺. Crop productivity could be harnessed on acid soils via the development of plant varieties tolerant to phytotoxic levels of these cations. In this study, we investigated the extent of natural variation for Mn²⁺ tolerance among ten parental lines of the Australian and International canola mapping populations. Response to Mn²⁺ toxicity was measured on the bases of cotyledon chlorosis, shoot biomass, and leaf area in nutrient solution under control (9 μM of MnCl₂⋅4H₂O) and Mn treatment (125 μM of MnCl₂⋅4H₂O). Among parental lines, we selected Darmor-bzh and Yudal that showed significant and contrasting variation in Mn²⁺ tolerance to understand genetic control and identify the quantitative trait loci (QTL) underlying Mn²⁺ tolerance. We evaluated parental lines and their doubled haploid (DH) progenies (196 lines) derived from an F₁ cross, Darmor-bzh/Yudal for Mn²⁺ tolerance. Mn²⁺-tolerant genotypes had significantly higher shoot biomass and leaf area compared to Mn²⁺-sensitive genotypes. A genetic linkage map based on 7,805 DArTseq markers corresponding to 2,094 unique loci was constructed and further utilized for QTL identification. A major locus, BnMn²⁺.A09 was further mapped with a SNP marker, Bn-A09-p29012402 (LOD score of 34.6) accounting for most of the variation in Mn²⁺ tolerance on chromosome A09. This is the first report on the genomic localization of a Mn²⁺ tolerance locus in B. napus. Additionally, an ortholog of A. thaliana encoding for cation efflux facilitator transporter was located within 3,991 bp from significant SNP marker associated with BnMn²⁺.A09. A suite of genome sequence based markers (DArTseq and Illumina Infinium SNPs) flanking the BnMn²⁺.A09 locus would provide an invaluable tool for various molecular breeding applications to improve canola production and profitability on Mn²⁺ toxic soils.

Transcriptome Analyses Reveal Candidate Pod Shattering-Associated Genes Involved in the Pod Ventral Sutures of Common Vetch (Vicia sativa L.)

The seed dispersion caused by pod shattering is a form of propagation used by many wild species. Loss of seeds from pod shattering is frequent in the common vetch (Vicia sativa L.), an important self-pollinating annual forage legume. However, pod shattering is one of the most important defects that limits the reproduction of the vetch in the field and the usage as a leguminous forage crop. To better understand the vetch pod shattering mechanism, we used high-throughput RNA sequencing to assess the global changes in the transcriptomes of the pod ventral sutures of shattering-susceptible and shattering-resistant vetch accessions screened from 541 vetch germplasms. A total of 1,285 significantly differentially expressed unigenes (DEGs) were detected, including 575 up-regulated unigenes and 710 down-regulated unigenes. Analyses of Gene Ontology and KEGG metabolic enrichment pathways of 1,285 DEGs indicated that 22 DEGs encoding cell wall modifications and hydrolases associated with pod shattering were highly expressed in shattering-susceptible accessions. These genes were mainly enriched in "hydrolase activity," "cytoplasm," and "carbohydrate metabolic process" systems. These cell wall modifications and hydrolases genes included β-glucosidase and endo-polygalacturonase, which work together to break down the glycosidic bonds of pectin and cellulose, and to promote the dissolution and disappearance of the cell wall in the ventral suture of the pod and make the pod more susceptible to shattering. We demonstrated the differences in gene transcription levels between the shattering-susceptible and shattering-resistant vetch accessions for the first time and our results provided valuable information for the identifying and characterizing of pod shattering regulation networks in vetch. This information may facilitate the future identification of pod shattering-related genes and their underlying molecular mechanisms in the common vetch.

Molecular Diversity Analysis and Genetic Mapping of Pod Shatter Resistance Loci in Brassica carinata L

Seed lost due to easy pod dehiscence at maturity (pod shatter) is a major problem in several members of Brassicaceae family. We investigated the level of pod shatter resistance in Ethiopian mustard (Brassica carinata) and identified quantitative trait loci (QTL) for targeted introgression of this trait in Ethiopian mustard and its close relatives of the genus Brassica. A set of 83 accessions of B. carinata, collected from the Australian Grains Genebank, was evaluated for pod shatter resistance based on pod rupture energy (RE). In comparison to B. napus (RE = 2.16 mJ), B. carinata accessions had higher RE values (2.53 to 20.82 mJ). A genetic linkage map of an F₂ population from two contrasting B. carinata selections, BC73526 (shatter resistant with high RE) and BC73524 (shatter prone with low RE) comprising 300 individuals, was constructed using a set of 6,464 high quality DArTseq markers and subsequently used for QTL analysis. Genetic analysis of the F₂ and F_2:3 derived lines revealed five statistically significant QTL (LOD ≥ 3) that are linked with pod shatter resistance on chromosomes B1, B3, B8, and C5. Herein, we report for the first time, identification of genetic loci associated with pod shatter resistance in B. carinata. These characterized accessions would be useful in Brassica breeding programs for introgression of pod shatter resistance alleles in to elite breeding lines. Molecular markers would assist marker-assisted selection for tracing the introgression of resistant alleles. Our results suggest that the value of the germplasm collections can be harnessed through genetic and genomics tools.

Identification of BnaYUCCA6 as a candidate gene for branch angle in Brassica napus by QTL-seq

Oilseed rape (Brassica napus L.) is one of the most important oil crops in China as well as worldwide. Branch angle as a plant architecture component trait plays an important role for high density planting and yield performance. In this study, bulked segregant analysis (BSA) combined with next generation sequencing technology was used to fine map QTL for branch angle. A major QTL, designated as branch angle 1 (ba1) was identified on A06 and further validated by Indel marker-based classical QTL mapping in an F₂ population. Eighty-two genes were identified in the ba1 region. Among these genes, BnaA0639380D is a homolog of AtYUCCA6. Sequence comparison of BnaA0639380D from small- and big-branch angle oilseed rape lines identified six SNPs and four amino acid variation in the promoter and coding region, respectively. The expression level of BnaA0639380D is significantly higher in the small branch angle line Purler than in the big branch angle line Huyou19, suggesting that the genomic mutations may result in reduced activity of BnaA0639380D in Huyou19. Phytohormone determination showed that the IAA content in Purler was also obviously increased. Taken together, our results suggested BnaA0639380D is a possible candidate gene for branch angle in oilseed rape.

Sex-linked markers in the North American green frog (Rana clamitans) developed using DArTseq provide early insight into sex chromosome evolution

BACKGROUND

The extent to which sex reversal is associated with transitions in sex determining systems (XX-XY, ZZ-ZW, etc.) or abnormal sexual differentiation is predominantly unexplored in amphibians. This is in large part because most amphibian taxa have homomorphic sex chromosomes, which has traditionally made it challenging to identify discordance between phenotypic and genetic sex in amphibians, despite all amphibians having a genetic component to sex determination. Recent advances in molecular techniques such as genome complexity reduction and high throughput sequencing present a valuable avenue for furthering our understanding of sex determination in amphibians and other taxa with homomorphic sex chromosomes like many fish and reptiles.

RESULTS

We use DArTseq as a novel approach to identify sex-linked markers in the North American green frog (Rana clamitans melanota) using lab-reared tadpoles as well as wild-caught adults from seven ponds either in undeveloped, forested habitats or suburban ponds known to be subject to contamination by anthropogenic chemicals. The DArTseq methodology identified 13 sex-linked SNP loci and eight presence-absence loci associated with males, indicating an XX-XY system. Both alleles from a single locus show partial high sequence homology to Dmrt1, a gene linked to sex determination and differentiation throughout Metazoa. Two other loci have sequence similarities to regions of the chimpanzee and human X-chromosome as well as the chicken Z-chromosome. Several loci also show geographic variation in sex-linkage, possibly indicating sex chromosome recombination. While all loci are statistically sex-linked, they show varying degrees of female heterozygosity and male homozygosity, providing further evidence that some markers are on regions of the sex chromosomes undergoing higher rates of recombination and therefore further apart from the putative sex determining locus.

CONCLUSION

The ease of the DArTseq platform provides a useful avenue for future research on sex reversal and sex chromosome evolution in vertebrates, particularly for non-model species with homomorphic or cryptic or nascent sex chromosomes.

Who's your mama? Riverine hybridisation of threatened freshwater Trout Cod and Murray Cod

Rates of hybridization and introgression are increasing dramatically worldwide because of translocations, restocking of organisms and habitat modifications; thus, determining whether hybridization is occuring after reintroducing extirpated congeneric species is commensurately important for conservation. Restocking programs are sometimes criticized because of the genetic consequences of hatchery-bred fish breeding with wild populations. These concerns are important to conservation restocking programs, including those from the Australian freshwater fish family, Percichthyidae. Two of the better known Australian Percichthyidae are the Murray Cod, Maccullochella peelii and Trout Cod, Maccullochella macquariensis which were formerly widespread over the Murray Darling Basin. In much of the Murrumbidgee River, Trout Cod and Murray Cod were sympatric until the late 1970s when Trout Cod were extirpated. Here we use genetic single nucleotide polymorphism (SNP) data together with mitochondrial sequences to examine hybridization and introgression between Murray Cod and Trout Cod in the upper Murrumbidgee River and consider implications for restocking programs. We have confirmed restocked riverine Trout Cod reproducing, but only as inter-specific matings, in the wild. We detected hybrid Trout Cod-Murray Cod in the Upper Murrumbidgee, recording the first hybrid larvae in the wild. Although hybrid larvae, juveniles and adults have been recorded in hatcheries and impoundments, and hybrid adults have been recorded in rivers previously, this is the first time fertile F1 have been recorded in a wild riverine population. The F1 backcrosses with Murray cod have also been found to be fertile. All backcrosses noted were with pure Murray Cod. Such introgression has not been recorded previously in these two species, and the imbalance in hybridization direction may have important implications for restocking programs.

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.

A Genome-Wide Association Study Reveals New Loci for Resistance to Clubroot Disease in Brassica napus

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. However, the yield and quality of rapeseed were largely decreased by clubroot (Plasmodiophora brassicae Woronin). Therefore, it is of great importance for screening more resistant germplasms or genes and improving the resistance to P. brassicae in rapeseed breeding. In this study, a massive resistant identification for a natural global population was conducted in two environments with race/pathotype 4 of P. brassicae which was the most predominant in China, and a wide range of phenotypic variation was found in the population. In addition, a genome-wide association study of 472 accessions for clubroot resistance (CR) was performed with 60K Brassica Infinium SNP arrays for the first time. In total, nine QTLs were detected, seven of which were novel through integrative analysis. Furthermore, additive effects in genetic control of CR in rapeseed among the above loci were found. By bioinformatic analyses, the candidate genes of these loci were predicted, which indicated that TIR-NBS gene family might play an important role in CR. It is believable that the results presented in our study could provide valuable information for understanding the genetic mechanism and molecular regulation of CR.

A Genome-Wide Association Study Reveals New Loci for Resistance to Clubroot Disease in Brassica napus

Rapeseed (Brassica napus L.) is one of the most important oil crops in the world. However, the yield and quality of rapeseed were largely decreased by clubroot (Plasmodiophora brassicae Woronin). Therefore, it is of great importance for screening more resistant germplasms or genes and improving the resistance to P. brassicae in rapeseed breeding. In this study, a massive resistant identification for a natural global population was conducted in two environments with race/pathotype 4 of P. brassicae which was the most predominant in China, and a wide range of phenotypic variation was found in the population. In addition, a genome-wide association study of 472 accessions for clubroot resistance (CR) was performed with 60K Brassica Infinium SNP arrays for the first time. In total, nine QTLs were detected, seven of which were novel through integrative analysis. Furthermore, additive effects in genetic control of CR in rapeseed among the above loci were found. By bioinformatic analyses, the candidate genes of these loci were predicted, which indicated that TIR-NBS gene family might play an important role in CR. It is believable that the results presented in our study could provide valuable information for understanding the genetic mechanism and molecular regulation of CR.

Genome-Wide Association Study Dissecting the Genetic Architecture Underlying the Branch Angle Trait in Rapeseed (Brassica napus L.)

The rapeseed branch angle is an important morphological trait because an adequate branch angle enables more efficient light capture under high planting densities. Here, we report that the average angle of the five top branches provides a reliable representation of the average angle of all branches. Statistical analyses revealed a significantly positive correlation between the branch angle and multiple plant-type and yield-related traits. The 60 K Brassica Infinium^® single nucleotide polymorphism (SNP) array was utilized to genotype an association panel with 520 diverse accessions. A genome-wide association study was performed to determine the genetic architecture of branch angle, and 56 loci were identified as being significantly associated with the branch angle trait via three models, including a robust, novel, nonparametric Anderson-Darling (A-D) test. Moreover, these loci explained 51.1% of the phenotypic variation when a simple additive model was applied. Within the linkage disequilibrium (LD) decay ranges of 53 loci, we observed plausible candidates orthologous to documented Arabidopsis genes, such as LAZY1, SGR2, SGR4, SGR8, SGR9, PIN3, PIN7, CRK5, TIR1, and APD7. These results provide insight into the genetic basis of the branch angle trait in rapeseed and might facilitate marker-based breeding for improvements in plant architecture.

Genomic selection to resistance to Stenocarpella maydis in maize lines using DArTseq markers

BACKGROUND

The identification of lines resistant to ear diseases is of great importance in maize breeding because such diseases directly interfere with kernel quality and yield. Among these diseases, ear rot disease is widely relevant due to significant decrease in grain yield. Ear rot may be caused by the fungus Stenocarpella maydi; however, little information about genetic resistance to this pathogen is available in maize, mainly related to candidate genes in genome. In order to exploit this genome information we used 23.154 Dart-seq markers in 238 lines and apply genome-wide selection to select resistance genotypes. We divide the lines into clusters to identify groups related to resistance to Stenocarpella maydi and use Bayesian stochastic search variable approach and rr-BLUP methods to comparate their selection results.

RESULTS

Through a principal component analysis (PCA) and hierarchical clustering, it was observed that the three main genetic groups (Stiff Stalk Synthetic, Non-Stiff Stalk Synthetic and Tropical) were clustered in a consistent manner, and information on the resistance sources could be obtained according to the line of origin where populations derived from genetic subgroup Suwan presenting higher levels of resistance. The ridge regression best linear unbiased prediction (rr-BLUP) and Bayesian stochastic search variable (BSSV) models presented equivalent abilities regarding predictive processes.

CONCLUSION

Our work showed that is possible to select maize lines presenting a high resistance to Stenocarpella maydis. This claim is based on the acceptable level of predictive accuracy obtained by Genome-wide Selection (GWS) using different models. Furthermore, the lines related to background Suwan present a higher level of resistance than lines related to other groups.

Molecular-cytogenetic characterization of C-genome chromosome substitution lines in Brassica juncea (L.) Czern and Coss

C genome chromosome substitution lines of B. juncea constitute a key genetic resource for increased genetic diversity and hybrid performance. C genome chromosome substitution lines were found in the progenies of derived B. juncea (2n = 36; AABB), synthesized through hybridization between B. napus and B. carinata. These were originally recognized based on the morphology and genomic in situ hybridization. Genotyping using the Brassica Illumina 60K Infinium SNP array confirmed the presence of C genome chromosomes in a large number of derived B. juncea genotypes. Three whole chromosome substitutions and 13 major C genome fragment substitutions were identified. Fragment substitutions were primarily terminal, but intercalary substitution(s) involving chromosome C1 and C2 were identified in three genotypes. The size of substituted C genome fragments varied from 0.04 Mbp (C1) to 64.85 Mbp (C3). In terms of proportions, these ranged from 0.10 % (C1) to 100 % (C1, C3 and C7) of the substituted chromosome. SSR genotyping with B genome specific primers suggested that substituting C genome chromosome(s) are likely to have replaced B genome chromosome(s). C1 was the most common substituting chromosome while substituted B chromosome seemed random. Study of the phenotypic traits underlined the importance of the substitution lines (especially of chromosome C1) for conferring superior trait performance (main shoot length and pods on the main shoot). High heterosis was also indicated in hybrid combinations of substitution lines with natural B. juncea. These substitution genotypes constituted a valuable resource for targeted gene transfer and QTL identification.

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.

To Have and to Hold: Selection for Seed and Fruit Retention During Crop Domestication

Crop domestication provides a useful model system to characterize the molecular and developmental bases of morphological variation in plants. Among the most universal changes resulting from selection during crop domestication is the loss of seed and fruit dispersal mechanisms, which greatly facilitates harvesting efficiency. In this review, we consider the molecular genetic and developmental bases of the loss of seed shattering and fruit dispersal in six major crop plant families, three of which are primarily associated with seed crops (Poaceae, Brassicaceae, Fabaceae) and three of which are associated with fleshy-fruited crops (Solanaceae, Rosaceae, Rutaceae). We find that the developmental basis of the loss of seed/fruit dispersal is conserved in a number of independently domesticated crops, indicating the widespread occurrence of developmentally convergent evolution in response to human selection. With regard to the molecular genetic approaches used to characterize the basis of this trait, traditional biparental quantitative trait loci mapping remains the most commonly used strategy; however, recent advances in next-generation sequencing technologies are now providing new avenues to map and characterize loss of shattering/dispersal alleles. We anticipate that continued application of these approaches, together with candidate gene analyses informed by known shattering candidate genes from other crops, will lead to a rapid expansion of our understanding of this critical domestication trait.

Identification QTLs Controlling Genes for Se Uptake in Lentil Seeds

Lentil (Lens culinaris Medik.) is an excellent source of protein and carbohydrates and is also rich in essential trace elements for the human diet. Selenium (Se) is an essential micronutrient for human health and nutrition, providing protection against several diseases and regulating important biological systems. Dietary intake of 55 μg of Se per day is recommended for adults, with inadequate Se intake causing significant health problems. The objective of this study was to identify and map quantitative trait loci (QTL) of genes controlling Se accumulation in lentil seeds using a population of 96 recombinant inbred lines (RILs) developed from the cross “PI 320937” × “Eston” grown in three different environments for two years (2012 and 2013). Se concentration in seed varied between 119 and 883 μg/kg. A linkage map consisting of 1,784 markers (4 SSRs, and 1,780 SNPs) was developed. The map spanned a total length of 4,060.6 cM, consisting of 7 linkage groups (LGs) with an average distance of 2.3 cM between adjacent markers. Four QTL regions and 36 putative QTL markers, with LOD scores ranging from 3.00 to 4.97, distributed across two linkage groups (LG2 and LG5) were associated with seed Se concentration, explaining 6.3–16.9% of the phenotypic variation.