The causal mutation leading to sweetness in modern white lupin cultivars

Lupins are high-protein crops that are rapidly gaining interest as hardy alternatives to soybean; however, they accumulate antinutritional alkaloids of the quinolizidine type (QAs). Lupin domestication was enabled by the discovery of genetic loci conferring low QA levels (sweetness), but the precise identity of the underlying genes remains uncertain. We show that pauper, the most common sweet locus in white lupin, encodes an acetyltransferase (AT) unexpectedly involved in the early QA pathway. In pauper plants, a single-nucleotide polymorphism (SNP) strongly impairs AT activity, causing pathway blockage. We corroborate our hypothesis by replicating the pauper chemotype in narrow-leafed lupin via mutagenesis. Our work adds a new dimension to QA biosynthesis and establishes the identity of a lupin sweet gene for the first time, thus facilitating lupin breeding and enabling domestication of other QA-containing legumes.

FLOWERING LOCUS T indel variants confer vernalization-independent and photoperiod-insensitive flowering of yellow lupin (Lupinus luteus L.)

Ongoing climate change has considerably reduced the seasonal window for crop vernalization, concurrently expanding cultivation area into northern latitudes with long-day photoperiod. To address these changes, cool season legume breeders need to understand molecular control of vernalization and photoperiod. A key floral transition gene integrating signals from these pathways is the Flowering locus T (FT). Here, a recently domesticated grain legume, yellow lupin (Lupinus luteus L.), was explored for potential involvement of FT homologues in abolition of vernalization and photoperiod requirements. Two FTa (LlutFTa1a and LlutFTa1b) and FTc (LlutFTc1 and LlutFTc2) homologues were identified and sequenced for two contrasting parents of a reference recombinant inbred line (RIL) population, an early-flowering cultivar Wodjil and a late-flowering wild-type P28213. Large deletions were detected in the 5' promoter regions of three FT homologues. Quantitative trait loci were identified for flowering time and vernalization response in the RIL population and in a diverse panel of wild and domesticated accessions. A 2227 bp deletion found in the LlutFTc1 promoter was linked with early phenology and vernalization independence, whereas LlutFTa1a and LlutFTc2 indels with photoperiod responsiveness. Comparative mapping highlighted convergence of FTc1 indel evolution in two Old World lupin species, addressing both artificial selection during domestication and natural adaptation to short season environmental conditions. We concluded that rapid flowering in yellow lupin is associated with the de-repression of the LlutFTc1 homologue from the juvenile phase, putatively due to the elimination of all binding sites in the promoter region for the AGAMOUS-like 15 transcription factor.

A Trimethylguanosine Synthase1-like (TGS1) homologue is implicated in vernalisation and flowering time control

KEY MESSAGE

A plant-specific Trimethylguanosine Synthase1-like homologue was identified as a candidate gene for the efl mutation in narrow-leafed lupin, which alters phenology by reducing vernalisation requirement. The vernalisation pathway is a key component of flowering time control in plants from temperate regions but is not well understood in the legume family. Here we examined vernalisation control in the temperate grain legume species, narrow-leafed lupin (Lupinus angustifolius L.), and discovered a candidate gene for an ethylene imine mutation (efl). The efl mutation changes phenology from late to mid-season flowering and additionally causes transformation from obligate to facultative vernalisation requirement. The efl locus was mapped to pseudochromosome NLL-10 in a recombinant inbred line (RIL) mapping population developed by accelerated single seed descent. Candidate genes were identified in the reference genome, and a diverse panel of narrow-leafed lupins was screened to validate mutations specific to accessions with efl. A non-synonymous SNP mutation within an S-adenosyl-L-methionine-dependent methyltransferase protein domain of a Trimethylguanosine Synthase1-like (TGS1) orthologue was identified as the candidate mutation giving rise to efl. This mutation caused substitution of an amino acid within an established motif at a position that is otherwise highly conserved in several plant families and was perfectly correlated with the efl phenotype in F2 and F6 genetic population and a panel of diverse accessions, including the original efl mutant. Expression of the TGS1 homologue did not differ between wild-type and efl genotypes, supporting altered functional activity of the gene product. This is the first time a TGS1 orthologue has been associated with vernalisation response and flowering time control in any plant species.

Phenotypic characterisation and linkage mapping of domestication syndrome traits in yellow lupin (Lupinus luteus L.)

The transformation of wild plants into domesticated crops usually modifies a common set of characters referred to as 'domestication syndrome' traits such as the loss of pod shattering/seed dehiscence, loss of seed dormancy, reduced anti-nutritional compounds and changes in growth habit, phenology, flower and seed colour. Understanding the genetic control of domestication syndrome traits facilitates the efficient transfer of useful traits from wild progenitors into crops through crossing and selection. Domesticated forms of yellow lupin (Lupinus luteus L.) possess many domestication syndrome traits, while their genetic control remains a mystery. This study aimed to reveal the genetic control of yellow lupin domestication traits. This involved phenotypic characterisation of those traits, defining the genomic regions controlling domestication traits on a linkage map and performing a comparative genomic analysis of yellow lupin with its better-understood relatives, narrow-leafed lupin (L. angustifolius L.) and white lupin (L. albus L.). We phenotyped an F₉ recombinant inbred line (RIL) population of a wide cross between Wodjil (domesticated) × P28213 (wild). Vernalisation responsiveness, alkaloid content, flower and seed colour in yellow lupin were each found to be controlled by single loci on linkage groups YL-21, YL-06, YL-03 and YL-38, respectively. Aligning the genomes of yellow with narrow-leafed lupin and white lupin revealed well-conserved synteny between these sister species (76% and 71%, respectively). This genomic comparison revealed that one of the key domestication traits, vernalisation-responsive flowering, mapped to a region of conserved synteny with the vernalisation-responsive flowering time Ku locus of narrow-leafed lupin, which has previously been shown to be controlled by an FT homologue. In contrast, the loci controlling alkaloid content were each found at non-syntenic regions among the three species. This provides a first glimpse into the molecular control of flowering time in yellow lupin and demonstrates both the power and the limitation of synteny as a tool for gene discovery in lupins.

The First Genetic Map for a Psoraleoid Legume (Bituminaria bituminosa) Reveals Highly Conserved Synteny with Phaseoloid Legumes

We present the first genetic map of tedera (Bituminaria bituminosa (L.) C.H. Stirton), a drought-tolerant forage legume from the Canary Islands with useful pharmaceutical properties. It is also the first genetic map for any species in the tribe Psoraleeae (Fabaceae). The map comprises 2042 genotyping-by-sequencing (GBS) markers distributed across 10 linkage groups, consistent with the haploid chromosome count for this species (n = 10). Sequence tags from the markers were used to find homologous matches in the genome sequences of the closely related species in the Phaseoleae tribe: soybean, common bean, and cowpea. No tedera linkage groups align in their entirety to chromosomes in any of these phaseoloid species, but there are long stretches of collinearity that could be used in tedera research for gene discovery purposes using the better-resourced phaseoloid species. Using Ks analysis of a tedera transcriptome against five legume genomes provides an estimated divergence time of 17.4 million years between tedera and soybean. Genomic information and resources developed here will be invaluable for breeding tedera varieties for forage and pharmaceutical purposes.

The first genetic map for yellow lupin enables genetic dissection of adaptation traits in an orphan grain legume crop

BACKGROUND

Yellow lupin (Lupinus luteus L.) is a promising grain legume for productive and sustainable crop rotations. It has the advantages of high tolerance to soil acidity and excellent seed quality, but its current yield potential is poor, especially in low rainfall environments. Key adaptation traits such as phenology and enhanced stress tolerance are often complex and controlled by several genes. Genomic-enabled technologies may help to improve our basic understanding of these traits and to provide selective markers in breeding. However, in yellow lupin there are very limited genomic resources to support research and no published information is available on the genetic control of adaptation traits.

RESULTS

We aimed to address these deficiencies by developing the first linkage map for yellow lupin and conducting quantitative trait locus (QTL) analysis of yield under well-watered (WW) and water-deficit (WT) conditions. Two next-generation sequencing marker approaches - genotyping-by-sequencing (GBS) and Diversity Array Technology (DArT) sequencing - were employed to genotype a recombinant inbred line (RIL) population developed from a bi-parental cross between wild and domesticated parents. A total of 2,458 filtered single nucleotide polymorphism (SNP) and presence / absence variation (PAV) markers were used to develop a genetic map comprising 40 linkage groups, the first reported for this species. A number of significant QTLs controlling total biomass and 100-seed weight under two water (WW and WD) regimes were found on linkage groups YL-03, YL-09 and YL-26 that together explained 9 and 28% of total phenotypic variability. QTLs associated with length of the reproductive phase and time to flower were found on YL-01, YL-21, YL-35 and YL-40 that together explained a total of 12 and 44% of total phenotypic variation.

CONCLUSION

These genomic resources and the QTL information offer significant potential for use in marker-assisted selection in yellow lupin.

Prediction of pathogenicity genes involved in adaptation to a lupin host in the fungal pathogens Botrytis cinerea and Sclerotinia sclerotiorum via comparative genomics

Background

Narrow-leafed lupin is an emerging crop of significance in agriculture, livestock feed and human health food. However, its susceptibility to various diseases is a major obstacle towards increased adoption. Sclerotinia sclerotiorum and Botrytis cinerea – both necrotrophs with broad host-ranges - are reported among the top 10 lupin pathogens. Whole-genome sequencing and comparative genomics are useful tools to discover genes responsible for interactions between pathogens and their hosts.

Results

Genomes were assembled for one isolate of B. cinerea and two isolates of S. sclerotiorum, which were isolated from either narrow-leafed or pearl lupin species. Comparative genomics analysis between lupin-derived isolates and others isolated from alternate hosts was used to predict between 94 to 98 effector gene candidates from among their respective non-conserved gene contents.

Conclusions

Detection of minor differences between relatively recently-diverged isolates, originating from distinct regions and with hosts, may highlight novel or recent gene mutations and losses resulting from host adaptation in broad host-range fungal pathogens.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5774-2) contains supplementary material, which is available to authorized users.

INDEL variation in the regulatory region of the major flowering time gene LanFTc1 is associated with vernalization response and flowering time in narrow-leafed lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) cultivation was transformed by 2 dominant vernalization-insensitive, early flowering time loci known as Ku and Julius (Jul), which allowed expansion into shorter season environments. However, reliance on these loci has limited genetic and phenotypic diversity for environmental adaptation in cultivated lupin. We recently predicted that a 1,423-bp deletion in the cis-regulatory region of LanFTc1, a FLOWERING LOCUS T (FT) homologue, derepressed expression of LanFTc1 and was the underlying cause of the Ku phenotype. Here, we surveyed diverse germplasm for LanFTc1 cis-regulatory variation and identified 2 further deletions of 1,208 and 5,162 bp in the 5' regulatory region, which overlap the 1,423-bp deletion. Additionally, we confirmed that no other polymorphisms were perfectly associated with vernalization responsiveness. Phenotyping and gene expression analyses revealed that Jul accessions possessed the 5,162-bp deletion and that the Jul and Ku deletions were equally capable of removing vernalization requirement and up-regulating gene expression. The 1,208-bp deletion was associated with intermediate phenology, vernalization responsiveness, and gene expression and therefore may be useful for expanding agronomic adaptation of lupin. This insertion/deletion series may also help resolve how the vernalization response is mediated at the molecular level in legumes.

Nondestructive Phenomic Tools for the Prediction of Heat and Drought Tolerance at Anthesis in Brassica Species

Oilseed Brassica species are vulnerable to heat and drought stress, especially in the early reproductive stage. We evaluated plant imaging of whole plant and flower tissue, leaf stomatal conductance, leaf and bud temperature, photochemical reflectance index, quantum yield of photosynthesis, and leaf gas exchange for their suitability to detect tolerance to heat (H) and/or drought (D) stress treatments in 12 Brassica genotypes (G). A replicated factorial experiment was set up with 7 d of stress treatment from the beginning of anthesis with various levels of three factors H, D, and G. Most phenomics tools detected plant stress as indicated by significant main effects of H, D, and H×D. Whole plant volume was highly correlated with fresh weight changes, suggesting that whole plant imaging may be a useful surrogate for fresh weight in future studies. Vcmax, the maximum carboxylation rate of photosynthesis, increased rapidly on day 1 in H and H+D treatments, and there were significant interactions of G×H and G×D. Vcmax of genotypes on day 1 in H and H+D treatments was positively correlated with their harvested seed yield. Vcmax on day 1 and day 3 were clustered with seed yield in H and H+D treatments as shown in the heatmaps of genotypic correlations. TPU, the rate of triose phosphate use, also showed significant positive genotypic correlations with seed yield in H+D treatments. Flower volume showed significant interactions of G×H and G×D on day 7, and flower volume of genotypes on day 7 in H was positively correlated with their harvested seed yield. There were few interactions of G×H or G×D for leaf stomatal conductance, leaf and bud temperature, photochemical reflectance index, and quantum yield of photosynthesis. Vcmax, TPU, and volume of flowers are potential nondestructive phenomic traits for heat or combined heat and drought stress tolerance screening in Brassica germplasm.

Adapting legume crops to climate change using genomic approaches

Our agricultural system and hence food security is threatened by combination of events, such as increasing population, the impacts of climate change, and the need to a more sustainable development. Evolutionary adaptation may help some species to overcome environmental changes through new selection pressures driven by climate change. However, success of evolutionary adaptation is dependent on various factors, one of which is the extent of genetic variation available within species. Genomic approaches provide an exceptional opportunity to identify genetic variation that can be employed in crop improvement programs. In this review, we illustrate some of the routinely used genomics-based methods as well as recent breakthroughs, which facilitate assessment of genetic variation and discovery of adaptive genes in legumes. Although additional information is needed, the current utility of selection tools indicate a robust ability to utilize existing variation among legumes to address the challenges of climate uncertainty.

The western Mediterranean region provided the founder population of domesticated narrow-leafed lupin

This study revealed that the western Mediterranean provided the founder population for domesticated narrow-leafed lupin and that genetic diversity decreased significantly during narrow-leafed lupin domestication. The evolutionary history of plants during domestication profoundly shaped the genome structure and genetic diversity of today's crops. Advances in next-generation sequencing technologies allow unprecedented opportunities to understand genome evolution in minor crops, which constitute the majority of plant domestications. A diverse set of 231 wild and domesticated narrow-leafed lupin (Lupinus angustifolius L.) accessions were subjected to genotyping-by-sequencing using diversity arrays technology. Phylogenetic, genome-wide divergence and linkage disequilibrium analyses were applied to identify the founder population of domesticated narrow-leafed lupin and the genome-wide effect of domestication on its genome. We found wild western Mediterranean population as the founder of domesticated narrow-leafed lupin. Domestication was associated with an almost threefold reduction in genome diversity in domesticated accessions compared to their wild relatives. Selective sweep analysis identified no significant footprints of selection around domestication loci. A genome-wide association study identified single nucleotide polymorphism markers associated with pod dehiscence. This new understanding of the genomic consequences of narrow-leafed lupin domestication along with molecular marker tools developed here will assist plant breeders more effectively access wild genetic diversity for crop improvement.

Exploring the genetic and adaptive diversity of a pan-Mediterranean crop wild relative: narrow-leafed lupin

KEY MESSAGE

This first pan-Mediterranean analysis of genetic diversity in wild narrow-leafed lupin revealed strong East-West genetic differentiation of populations, an historic eastward migration, and signatures of genetic adaptation to climatic variables. Most grain crops suffer from a narrow genetic base, which limits their potential for adapting to new challenges such as increased stresses associated with climate change. Plant breeders are returning to the wild ancestors of crops and their close relatives to broaden the genetic base of their crops. Understanding the genetic adaptation of these wild relatives will help plant breeders most effectively use available wild diversity. Here, we took narrow-leafed lupin (Lupinus angustifolius L.) as a model to understand adaptation in a wild crop ancestor. A set of 142 wild accessions of narrow-leafed lupin from across the Mediterranean basin were subjected to genotyping-by-sequencing using Diversity Arrays Technology. Phylogenetic, linkage disequilibrium and demographic analyses were employed to explore the history of narrow-leafed lupin within the Mediterranean region. We found strong genetic differentiation between accessions from the western and eastern Mediterranean, evidence of an historic West to East migration, and that eastern Mediterranean narrow-leafed lupin experienced a severe and recent genetic bottleneck. We showed that these two populations differ for flowering time as a result of local adaptation, with the West flowering late while the East flowers early. A genome-wide association study identified single nucleotide polymorphism markers associated with climatic adaptation. Resolving the origin of wild narrow-leafed lupin and how its migration has induced adaptation to specific regions of the Mediterranean serves as a useful resource not only for developing narrow-leafed lupin cultivars with greater resilience to a changing climate, but also as a model which can be applied to other legumes.

A transcriptome-SNP-derived linkage map of Apios americana (potato bean) provides insights about genome re-organization and synteny conservation in the phaseoloid legumes

KEY MESSAGE

We report a linkage map for Apios americana and describe synteny with selected warm-season legumes. A translocation event in common bean and soybean is confirmed against Apios and Vigna species. Apios (Apios americana; "apios"), a tuberous perennial legume in the Phaseoleae tribe, was widely used as a food by Native Americans. Work in the last 40 years has led to several improved breeding lines. Aspects of the pollination biology (complex floral structure and tripping mechanism) have made controlled crosses difficult, and the previous reports indicated that the plant is likely primarily an outcrosser. We used a pseudo-testcross strategy to construct a genetic map specific to the maternal parent. The map was built using single-nucleotide polymorphism markers identified by comparing the expressed sequences of individuals in the mapping population against a de novo maternal reference transcriptome assembly. The apios map consists of 11 linkage groups and 1121 recombinationally distinct loci, covering ~ 938.6 cM. By sequencing the transcriptomes of all potential pollen parents, we were able to identify the probable pollen donors and to discover new aspects of the pollination biology in apios. No selfing was observed, but multiple pollen parents were seen within individual pods. Comparisons with genome sequences in other species in the Phaseoleae showed extended synteny for most apios linkage groups. This synteny supports the robustness of the map, and also sheds light on the history of the Phaseoleae, as apios is relatively early diverging in this tribe. We detected a translocation event that separates apios and two Vigna species from Phaseolus vulgaris and Glycine max. This apios mapping work provides a general protocol for sequencing-based construction of high-density linkage maps in outcrossing species with heterogeneous pollen parents.

Drought-Tolerant Brassica rapa Shows Rapid Expression of Gene Networks for General Stress Responses and Programmed Cell Death Under Simulated Drought Stress

Production of oilseed rape/canola (Brassica napus) is increasingly threatened by dry conditions while the demand for vegetable oil is increasing. Brassica rapa is a genetically diverse ancestor of B. napus, and is readily crossed with B. napus. Recently, we reported promising levels of drought tolerance in a wild type of B. rapa which could be a source of drought tolerance for B. napus. We analysed global gene expression by messenger RNA sequencing in seedlings of the drought-tolerant and a drought-sensitive genotype of B. rapa under simulated drought stress and control conditions. A subset of stress-response genes were validated by reverse transcription quantitative PCR. Gene ontology enrichment analysis and pathway enrichment analysis revealed major differences between the two genotypes in the mode and onset of stress responses in the first 12 h of treatment. Drought-tolerant plants reacted uniquely and rapidly by upregulating genes associated with jasmonic acid and salicylic acid metabolism, as well as genes known to cause endoplasmic reticulum stress and induction of programmed cell death. Conversely, active responses in drought-sensitive plants were delayed until 8 or 12 h after stress application. The results may help to identify biomarkers for selection of breeding materials with potentially improved drought tolerance.

A comprehensive draft genome sequence for lupin (Lupinus angustifolius), an emerging health food: insights into plant-microbe interactions and legume evolution

Lupins are important grain legume crops that form a critical part of sustainable farming systems, reducing fertilizer use and providing disease breaks. It has a basal phylogenetic position relative to other crop and model legumes and a high speciation rate. Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is gaining popularity as a health food, which is high in protein and dietary fibre but low in starch and gluten-free. We report the draft genome assembly (609 Mb) of NLL cultivar Tanjil, which has captured >98% of the gene content, sequences of additional lines and a dense genetic map. Lupins are unique among legumes and differ from most other land plants in that they do not form mycorrhizal associations. Remarkably, we find that NLL has lost all mycorrhiza-specific genes, but has retained genes commonly required for mycorrhization and nodulation. In addition, the genome also provided candidate genes for key disease resistance and domestication traits. We also find evidence of a whole-genome triplication at around 25 million years ago in the genistoid lineage leading to Lupinus. Our results will support detailed studies of legume evolution and accelerate lupin breeding programmes.

The loss of vernalization requirement in narrow-leafed lupin is associated with a deletion in the promoter and de-repressed expression of a Flowering Locus T (FT) homologue

Adaptation of Lupinus angustifolius (narrow-leafed lupin) to cropping in southern Australian and northern Europe was transformed by a dominant mutation (Ku) that removed vernalization requirement for flowering. The Ku mutation is now widely used in lupin breeding to confer early flowering and maturity. We report here the identity of the Ku mutation. We used a range of genetic, genomic and gene expression approaches to determine whether Flowering Locus T (FT) homologues are associated with the Ku locus. One of four FT homologues present in the narrow-leafed lupin genome, LanFTc1, perfectly co-segregated with the Ku locus in a reference mapping population. Expression of LanFTc1 in the ku (late-flowering) parent was strongly induced by vernalization, in contrast to the Ku (early-flowering) parent, which showed constitutively high LanFTc1 expression. Co-segregation of this expression phenotype with the LanFTc1 genotype indicated that the Ku mutation impairs cis-regulation of LanFTc1. Sequencing of LanFTc1 revealed a 1.4-kb deletion in the promoter region, which was perfectly predictive of vernalization response in 216 wild and domesticated accessions. Linkage disequilibrium rapidly decayed around LanFTc1, suggesting that this deletion caused the loss of vernalization response. This is the first time a legume FTc subclade gene has been implicated in the vernalization response.

Expansion of the phosphatidylethanolamine binding protein family in legumes: a case study of Lupinus angustifolius L. FLOWERING LOCUS T homologs, LanFTc1 and LanFTc2

BACKGROUND

The Arabidopsis FLOWERING LOCUS T (FT) gene, a member of the phosphatidylethanolamine binding protein (PEBP) family, is a major controller of flowering in response to photoperiod, vernalization and light quality. In legumes, FT evolved into three, functionally diversified clades, FTa, FTb and FTc. A milestone achievement in narrow-leafed lupin (Lupinus angustifolius L.) domestication was the loss of vernalization responsiveness at the Ku locus. Recently, one of two existing L. angustifolius homologs of FTc, LanFTc1, was revealed to be the gene underlying Ku. It is the first recorded involvement of an FTc homologue in vernalization. The evolutionary basis of this phenomenon in lupin has not yet been deciphered.

RESULTS

Bacterial artificial chromosome (BAC) clones carrying LanFTc1 and LanFTc2 genes were localized in different mitotic chromosomes and constituted sequence-specific landmarks for linkage groups NLL-10 and NLL-17. BAC-derived superscaffolds containing LanFTc genes revealed clear microsyntenic patterns to genome sequences of nine legume species. Superscaffold-1 carrying LanFTc1 aligned to regions encoding one or more FT-like genes whereas superscaffold-2 mapped to a region lacking such a homolog. Comparative mapping of the L. angustifolius genome assembly anchored to linkage map localized superscaffold-1 in the middle of a 15 cM conserved, collinear region. In contrast, superscaffold-2 was found at the edge of a 20 cM syntenic block containing highly disrupted collinearity at the LanFTc2 locus. 118 PEBP-family full-length homologs were identified in 10 legume genomes. Bayesian phylogenetic inference provided novel evidence supporting the hypothesis that whole-genome and tandem duplications contributed to expansion of PEBP-family genes in legumes. Duplicated genes were subjected to strong purifying selection. Promoter analysis of FT genes revealed no statistically significant sequence similarity between duplicated copies; only RE-alpha and CCAAT-box motifs were found at conserved positions and orientations.

CONCLUSIONS

Numerous lineage-specific duplications occurred during the evolution of legume PEBP-family genes. Whole-genome duplications resulted in the origin of subclades FTa, FTb and FTc and in the multiplication of FTa and FTb copy number. LanFTc1 is located in the region conserved among all main lineages of Papilionoideae. LanFTc1 is a direct descendant of ancestral FTc, whereas LanFTc2 appeared by subsequent duplication.

Root trait diversity, molecular marker diversity, and trait-marker associations in a core collection of Lupinus angustifolius

Narrow-leafed lupin (Lupinus angustifolius L.) is the predominant grain legume crop in southern Australia, contributing half of the total grain legume production of Australia. Its yield in Australia is hampered by a range of subsoil constraints. The adaptation of lupin genotypes to subsoil constraints may be improved by selecting for optimal root traits from new and exotic germplasm sources. We assessed root trait diversity and genetic diversity of a core collection of narrow-leafed lupin (111 accessions) using 191 Diversity Arrays Technology (DArT) markers. The genetic relationship among accessions was determined using the admixture model in STRUCTURE. Thirty-eight root-associated traits were characterized, with 21 having coefficient of variation values >0.5. Principal coordinate analysis and cluster analysis of the DArT markers revealed broad diversity among the accessions. An ad hoc statistics calculation resulted in 10 distinct populations with significant differences among and within them (P < 0.001). The mixed linear model test in TASSEL showed a significant association between all root traits and some DArT markers, with the numbers of markers associated with an individual trait ranging from 2 to 13. The percentage of phenotypic variation explained by any one marker ranged from 6.4 to 21.8%, with 15 associations explaining >10% of phenotypic variation. The genetic variation values ranged from 0 to 7994, with 23 associations having values >240. Root traits such as deeper roots and lateral root proliferation at depth would be useful for this species for improved adaptation to drier soil conditions. This study offers opportunities for discovering useful root traits that can be used to increase the yield of Australian cultivars across variable environmental conditions.

Centromere Locations in Brassica A and C Genomes Revealed Through Half-Tetrad Analysis

Locating centromeres on genome sequences can be challenging. The high density of repetitive elements in these regions makes sequence assembly problematic, especially when using short-read sequencing technologies. It can also be difficult to distinguish between active and recently extinct centromeres through sequence analysis. An effective solution is to identify genetically active centromeres (functional in meiosis) by half-tetrad analysis. This genetic approach involves detecting heterozygosity along chromosomes in segregating populations derived from gametes (half-tetrads). Unreduced gametes produced by first division restitution mechanisms comprise complete sets of nonsister chromatids. Along these chromatids, heterozygosity is maximal at the centromeres, and homologous recombination events result in homozygosity toward the telomeres. We genotyped populations of half-tetrad-derived individuals (from Brassica interspecific hybrids) using a high-density array of physically anchored SNP markers (Illumina Brassica 60K Infinium array). Mapping the distribution of heterozygosity in these half-tetrad individuals allowed the genetic mapping of all 19 centromeres of the Brassica A and C genomes to the reference Brassica napus genome. Gene and transposable element density across the B. napus genome were also assessed and corresponded well to previously reported genetic map positions. Known centromere-specific sequences were located in the reference genome, but mostly matched unanchored sequences, suggesting that the core centromeric regions may not yet be assembled into the pseudochromosomes of the reference genome. The increasing availability of genetic markers physically anchored to reference genomes greatly simplifies the genetic and physical mapping of centromeres using half-tetrad analysis. We discuss possible applications of this approach, including in species where half-tetrads are currently difficult to isolate.

Characterization and mapping of LanrBo: a locus conferring anthracnose resistance in narrow-leafed lupin (Lupinus angustifolius L.)

A novel and highly effective source of anthracnose resistance in narrow-leafed lupin was identified. Resistance was shown to be governed by a single dominant locus. Molecular markers have been developed, which can be used for selecting resistant genotypes in lupin breeding. A screening for anthracnose resistance of a set of plant genetic resources of narrow-leafed lupin (Lupinus angustifolius L.) identified the breeding line Bo7212 as being highly resistant to anthracnose (Colletotrichum lupini). Segregation analysis indicated that the resistance of Bo7212 is inherited by a single dominant locus. The corresponding resistance gene was given the designation LanrBo. Previously published molecular anchor markers allowed us to locate LanrBo on linkage group NLL-11 of narrow-leafed lupin. Using information from RNAseq data obtained with inoculated resistant vs. susceptible lupin entries as well as EST-sequence information from the model genome Lotus japonicus, additional SNP and EST markers linked to LanrBo were derived. A bracket of two LanrBo-flanking markers allows for precise marker-assisted selection of the novel resistance gene in narrow-leafed lupin breeding programs.

Microspore culture reveals complex meiotic behaviour in a trigenomic Brassica hybrid

BACKGROUND

Development of synthetic allohexaploid Brassica (2n = AABBCC) would be beneficial for agriculture, as allelic contributions from three genomes could increase hybrid vigour and broaden adaptation. Microspore culture of a near-allohexaploid hybrid derived from the cross (B. napus × B. carinata) × B. juncea was undertaken in order to assess the frequency and distribution of homologous and homoeologous crossovers in this trigenomic hybrid. SNP and SSR molecular markers were used to detect inheritance of A, B and C genome alleles in microspore-derived (MD) progeny. SNP allele copy number was also assessed. The MD progeny were also compared to progeny derived by self-pollination and open-pollination for fertility (estimated by self-pollinated seed set and pollen viability) and DNA ploidy (measured by flow cytometry).

RESULTS

In the trigenomic hybrid, homologous chromosome pairs A(j)-A(n), B(j)-B(c) and C(n)-C(c) had similar meiotic crossover frequencies and segregation to that previously observed in established Brassica species, as demonstrated by marker haplotype analysis of the MD population. Homoeologous pairing between chromosomes A1-C1, A2-C2 and A7-C6 was detected at frequencies of 12-18 %, with other homoeologous chromosome regions associating from 8 % (A3-C3) to 0-1 % (A8-C8, A8-C9) of the time. Copy number analysis revealed eight instances of additional chromosomes and 20 instances of chromosomes present in one copy in somatically doubled MD progeny. Presence of chromosome A6 was positively correlated with self-pollinated seed set and pollen viability in the MD population. Many MD progeny were unable to produce self-pollinated seed (76 %) or viable pollen (53 %), although one MD plant produced 198 self-pollinated seeds. Average fertility was significantly lower in progeny obtained by microspore culture than progeny obtained by self-pollination or open-pollination, after excluding MD progeny which had not undergone chromosome doubling.

CONCLUSIONS

Based on SNP data analysis of the microspore-derived progeny, crossover frequency per chromosome in the allohexaploid hybrid was found to be similar to that in established Brassica species, suggesting that the higher chromosome number did not significantly disrupt cellular regulation of meiosis. SNP allele copy number analysis revealed the occurrence not only of homoeologous duplication/deletion events but also other cryptic duplications and deletions that may have been the result of mitotic instability. Microspore culture simplified the assessment of chromosome behaviour in the allohexaploid hybrid but yielded progeny with lower fertility and a greater range of ploidy levels compared to progeny obtained by self- or open-pollination.

Using the Animal Model to Accelerate Response to Selection in a Self-Pollinating Crop

We used the animal model in S0 (F1) recurrent selection in a self-pollinating crop including, for the first time, phenotypic and relationship records from self progeny, in addition to cross progeny, in the pedigree. We tested the model in Pisum sativum, the autogamous annual species used by Mendel to demonstrate the particulate nature of inheritance. Resistance to ascochyta blight (Didymella pinodes complex) in segregating S0 cross progeny was assessed by best linear unbiased prediction over two cycles of selection. Genotypic concurrence across cycles was provided by pure-line ancestors. From cycle 1, 102/959 S0 plants were selected, and their S1 self progeny were intercrossed and selfed to produce 430 S0 and 575 S2 individuals that were evaluated in cycle 2. The analysis was improved by including all genetic relationships (with crossing and selfing in the pedigree), additive and nonadditive genetic covariances between cycles, fixed effects (cycles and spatial linear trends), and other random effects. Narrow-sense heritability for ascochyta blight resistance was 0.305 and 0.352 in cycles 1 and 2, respectively, calculated from variance components in the full model. The fitted correlation of predicted breeding values across cycles was 0.82. Average accuracy of predicted breeding values was 0.851 for S2 progeny of S1 parent plants and 0.805 for S0 progeny tested in cycle 2, and 0.878 for S1 parent plants for which no records were available. The forecasted response to selection was 11.2% in the next cycle with 20% S0 selection proportion. This is the first application of the animal model to cyclic selection in heterozygous populations of selfing plants. The method can be used in genomic selection, and for traits measured on S0-derived bulks such as grain yield.

Dissecting Molecular Evolution in the Highly Diverse Plant Clade Caryophyllales Using Transcriptome Sequencing

Many phylogenomic studies based on transcriptomes have been limited to “single-copy” genes due to methodological challenges in homology and orthology inferences. Only a relatively small number of studies have explored analyses beyond reconstructing species relationships. We sampled 69 transcriptomes in the hyperdiverse plant clade Caryophyllales and 27 outgroups from annotated genomes across eudicots. Using a combined similarity- and phylogenetic tree-based approach, we recovered 10,960 homolog groups, where each was represented by at least eight ingroup taxa. By decomposing these homolog trees, and taking gene duplications into account, we obtained 17,273 ortholog groups, where each was represented by at least ten ingroup taxa. We reconstructed the species phylogeny using a 1,122-gene data set with a gene occupancy of 92.1%. From the homolog trees, we found that both synonymous and nonsynonymous substitution rates in herbaceous lineages are up to three times as fast as in their woody relatives. This is the first time such a pattern has been shown across thousands of nuclear genes with dense taxon sampling. We also pinpointed regions of the Caryophyllales tree that were characterized by relatively high frequencies of gene duplication, including three previously unrecognized whole-genome duplications. By further combining information from homolog tree topology and synonymous distance between paralog pairs, phylogenetic locations for 13 putative genome duplication events were identified. Genes that experienced the greatest gene family expansion were concentrated among those involved in signal transduction and oxidoreduction, including a cytochrome P450 gene that encodes a key enzyme in the betalain synthesis pathway. Our approach demonstrates a new approach for functional phylogenomic analysis in nonmodel species that is based on homolog groups in addition to inferred ortholog groups.

Transcriptome sequencing of different narrow-leafed lupin tissue types provides a comprehensive uni-gene assembly and extensive gene-based molecular markers

Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is an important grain legume crop that is valuable for sustainable farming and is becoming recognized as a human health food. NLL breeding is directed at improving grain production, disease resistance, drought tolerance and health benefits. However, genetic and genomic studies have been hindered by a lack of extensive genomic resources for the species. Here, the generation, de novo assembly and annotation of transcriptome datasets derived from five different NLL tissue types of the reference accession cv. Tanjil are described. The Tanjil transcriptome was compared to transcriptomes of an early domesticated cv. Unicrop, a wild accession P27255, as well as accession 83A:476, together being the founding parents of two recombinant inbred line (RIL) populations. In silico predictions for transcriptome-derived gene-based length and SNP polymorphic markers were conducted and corroborated using a survey assembly sequence for NLL cv. Tanjil. This yielded extensive indel and SNP polymorphic markers for the two RIL populations. A total of 335 transcriptome-derived markers and 66 BAC-end sequence-derived markers were evaluated, and 275 polymorphic markers were selected to genotype the reference NLL 83A:476 × P27255 RIL population. This significantly improved the completeness, marker density and quality of the reference NLL genetic map.

Multiple polyploidy events in the early radiation of nodulating and nonnodulating legumes

Unresolved questions about evolution of the large and diverse legume family include the timing of polyploidy (whole-genome duplication; WGDs) relative to the origin of the major lineages within the Fabaceae and to the origin of symbiotic nitrogen fixation. Previous work has established that a WGD affects most lineages in the Papilionoideae and occurred sometime after the divergence of the papilionoid and mimosoid clades, but the exact timing has been unknown. The history of WGD has also not been established for legume lineages outside the Papilionoideae. We investigated the presence and timing of WGDs in the legumes by querying thousands of phylogenetic trees constructed from transcriptome and genome data from 20 diverse legumes and 17 outgroup species. The timing of duplications in the gene trees indicates that the papilionoid WGD occurred in the common ancestor of all papilionoids. The earliest diverging lineages of the Papilionoideae include both nodulating taxa, such as the genistoids (e.g., lupin), dalbergioids (e.g., peanut), phaseoloids (e.g., beans), and galegoids (=Hologalegina, e.g., clovers), and clades with nonnodulating taxa including Xanthocercis and Cladrastis (evaluated in this study). We also found evidence for several independent WGDs near the base of other major legume lineages, including the Mimosoideae-Cassiinae-Caesalpinieae (MCC), Detarieae, and Cercideae clades. Nodulation is found in the MCC and papilionoid clades, both of which experienced ancestral WGDs. However, there are numerous nonnodulating lineages in both clades, making it unclear whether the phylogenetic distribution of nodulation is due to independent gains or a single origin followed by multiple losses.

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

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

Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea

Background

Brassica oleracea is a valuable vegetable species that has contributed to human health and nutrition for hundreds of years and comprises multiple distinct cultivar groups with diverse morphological and phytochemical attributes. In addition to this phenotypic wealth, B. oleracea offers unique insights into polyploid evolution, as it results from multiple ancestral polyploidy events and a final Brassiceae-specific triplication event. Further, B. oleracea represents one of the diploid genomes that formed the economically important allopolyploid oilseed, Brassica napus. A deeper understanding of B. oleracea genome architecture provides a foundation for crop improvement strategies throughout the Brassica genus.

Results

We generate an assembly representing 75% of the predicted B. oleracea genome using a hybrid Illumina/Roche 454 approach. Two dense genetic maps are generated to anchor almost 92% of the assembled scaffolds to nine pseudo-chromosomes. Over 50,000 genes are annotated and 40% of the genome predicted to be repetitive, thus contributing to the increased genome size of B. oleracea compared to its close relative B. rapa. A snapshot of both the leaf transcriptome and methylome allows comparisons to be made across the triplicated sub-genomes, which resulted from the most recent Brassiceae-specific polyploidy event.

Conclusions

Differential expression of the triplicated syntelogs and cytosine methylation levels across the sub-genomes suggest residual marks of the genome dominance that led to the current genome architecture. Although cytosine methylation does not correlate with individual gene dominance, the independent methylation patterns of triplicated copies suggest epigenetic mechanisms play a role in the functional diversification of duplicate genes.

High-resolution molecular karyotyping uncovers pairing between ancestrally related Brassica chromosomes

How do chromosomal regions with differing degrees of homology and homeology interact at meiosis? We provide a novel analytical method based on simple genetics principles which can help to answer this important question. This method interrogates high-throughput molecular marker data in order to infer chromosome behavior at meiosis in interspecific hybrids. We validated this method using high-resolution molecular marker karyotyping in two experimental Brassica populations derived from interspecific crosses among B. juncea, B. napus and B. carinata, using a single nucleotide polymorphism chip. This method of analysis successfully identified meiotic interactions between chromosomes sharing different degrees of similarity: full-length homologs; full-length homeologs; large sections of primary homeologs; and small sections of secondary homeologs. This analytical method can be applied to any allopolyploid species or fertile interspecific hybrid in order to detect meiotic associations. This genetic information can then be used to identify which genomic regions share functional homeology (i.e., retain enough similarity to allow pairing and segregation at meiosis). When applied to interspecific hybrids for which reference genome sequences are available, the question of how differing degrees of homology and homeology affect meiotic interactions may finally be resolved.

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

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

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

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

Delayed water loss and temperature rise in floral buds compared with leaves of Brassica rapa subjected to a transient water stress during reproductive development

Leaf canopy temperature has been proposed as a rapid selection tool for drought tolerance among crop genotypes. However, floral bud temperature may be a better indicator of drought tolerance than leaf temperature in grain crops. In this study, we examined whether the floral bud and leaves of Brassica rapa L. had similar stomatal characteristics and showed similar water loss during a drying cycle. We also compared the leaf and bud temperatures when the plants were exposed to a 10-day transient water stress during reproductive development that affected flower development, increased flower abortion, increased pod abortion and reduced yield by an average of 85%. The water loss of detached leaves and floral buds showed that the stomata on the leaves closed before those of the floral buds as the leaf water potential decreased. Consistent with the water loss studies, the temperature of the intact bud showed a delayed increase during the drying process compared with the leaves. This suggested that floral bud temperature could be a useful indicator of the water status of the reproductive organs of B. rapa.

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

Background

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

Results

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

Conclusions

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

Evidence from genome-wide simple sequence repeat markers for a polyphyletic origin and secondary centers of genetic diversity of Brassica juncea in China and India

The oilseed Brassica juncea is an important crop with a long history of cultivation in India and China. Previous studies have suggested a polyphyletic origin of B. juncea and more than one migration from the primary to secondary centers of diversity. We investigated molecular genetic diversity based on 99 simple sequence repeat markers in 119 oilseed B. juncea varieties from China, India, Europe, and Australia to test whether molecular differentiation follows Vavilov's proposal of secondary centers of diversity in India and China. Two distinct groups were identified by markers in the A genome, and the same two groups were confirmed by markers in the B genome. Group 1 included accessions from central and western India, in addition to those from eastern China. Group 2 included accessions from central and western China, as well as those from northern and eastern India. European and Australian accessions were found only in Group 2. Chinese accessions had higher allelic diversity per accession (Group 1) and more private alleles per accession (Groups 1 and 2) than those from India. The marker data and geographic distribution of Groups 1 and 2 were consistent with two independent migrations of B. juncea from its center of origin in the Middle East and neighboring regions along trade routes to western China and northern India, followed by regional adaptation. Group 1 migrated further south and west in India, and further east in China, than Group 2. Group 2 showed diverse agroecological adaptation, with yellow-seeded spring-sown types in central and western China and brown-seeded autumn-sown types in India.

How a phosphorus-acquisition strategy based on carboxylate exudation powers the success and agronomic potential of lupines (Lupinus, Fabaceae)

Lupines (Lupinus species; Fabaceae) are an ancient crop with great potential to be developed further for high-protein feed and food, cover crops, and phytoremediation. Being legumes, they are capable of symbiotically fixing atmospheric nitrogen. However, Lupinus species appear to be nonmycorrhizal or weakly mycorrhizal at most; instead some produce cluster roots, which release vast amounts of phosphate-mobilizing carboxylates (inorganic anions). Other lupines produce cluster-like roots, which function in a similar manner, and some release large amounts of carboxylates without specialized roots. These traits associated with nutrient acquisition make lupines ideally suited for either impoverished soils or soils with large amounts of phosphorus that is poorly available for most plants, e.g., acidic or alkaline soils. Here we explore how common the nonmycorrhizal phosphorus-acquisition strategy based on exudation of carboxylates is in the genus Lupinus, concluding it is very likely more widespread than generally acknowledged. This trait may partly account for the role of lupines as pioneers or invasive species, but also makes them suitable crop plants while we reach "peak phosphorus".

The Brassica napus blackleg resistance gene LepR3 encodes a receptor-like protein triggered by the Leptosphaeria maculans effector AVRLM1

LepR3, found in the Brassica napus cv 'Surpass 400', provides race-specific resistance to the fungal pathogen Leptosphaeria maculans, which was overcome after great devastation in Australia in 2004. We investigated the LepR3 locus to identify the genetic basis of this resistance interaction. We employed a map-based cloning strategy, exploiting collinearity with the Arabidopsis thaliana and Brassica rapa genomes to enrich the map and locate a candidate gene. We also investigated the interaction of LepR3 with the L. maculans avirulence gene AvrLm1 using transgenics. LepR3 was found to encode a receptor-like protein (RLP). We also demonstrated that avirulence towards LepR3 is conferred by AvrLm1, which is responsible for both the Rlm1 and LepR3-dependent resistance responses in B. napus. LepR3 is the first functional B. napus disease resistance gene to be cloned. AvrLm1's interaction with two independent resistance loci, Rlm1 and LepR3, highlights the need to consider redundant phenotypes in 'gene-for-gene' interactions and offers an explanation as to why LepR3 was overcome so rapidly in parts of Australia.

Domestication bottlenecks limit genetic diversity and constrain adaptation in narrow-leafed lupin (Lupinus angustifolius L.)

In contrast to most widespread broad-acre crops, the narrow-leafed lupin (Lupinus angustifolius L.) was domesticated very recently, in breeding programmes isolated in both space and time. Whereas domestication was initiated in Central Europe in the early twentieth century, the crop was subsequently industrialized in Australia, which now dominates world production. To investigate the ramifications of these bottlenecks, the genetic diversity of wild (n = 1,248) and domesticated populations (n = 95) was characterized using diversity arrays technology, and adaptation studied using G × E trials (n = 31) comprising all Australian cultivars released from 1967 to 2004 (n = 23). Principal coordinates analysis demonstrates extremely limited genetic diversity in European and Australian breeding material compared to wild stocks. AMMI analysis indicates that G × E interaction is a minor, albeit significant effect, dominated by strong responses to local, Western Australian (WA) optima. Over time Australian cultivars have become increasingly responsive to warm, intermediate rainfall environments in the northern WA grainbelt, but much less so to cool vegetative phase eastern environments, which have considerably more yield potential. G × E interaction is well explained by phenology, and its interaction with seasonal climate, as a result of varying vernalization responses. Yield differences are minimized when vegetative phase temperatures fully satisfy the vernalization requirement (typical of eastern Australia), and maximized when they do not (typical of WA). In breeding for WA optima, the vernalization response has been eliminated and there has been strong selection for terminal drought avoidance through early phenology, which limits yield potential in longer season eastern environments. Conversely, vernalization-responsive cultivars are more yield-responsive in the east, where low temperatures moderately extend the vegetative phase. The confounding of phenology and vernalization response limits adaptation in narrow-leafed lupin, isolates breeding programmes, and should be eliminated by widening the flowering time range in a vernalization-unresponsive background. Concomitantly, breeding strategies that will widen the genetic base of the breeding pool in an ongoing manner should be initiated.

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

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines ‘Lynx-037DH’ and ‘Monty-028DH’. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.

Next generation DNA sequencing technology delivers valuable genetic markers for the genomic orphan legume species, Bituminaria bituminosa

BACKGROUND

Bituminaria bituminosa is a perennial legume species from the Canary Islands and Mediterranean region that has potential as a drought-tolerant pasture species and as a source of pharmaceutical compounds. Three botanical varieties have previously been identified in this species: albomarginata, bituminosa and crassiuscula. B. bituminosa can be considered a genomic 'orphan' species with very few genomic resources available. New DNA sequencing technologies provide an opportunity to develop high quality molecular markers for such orphan species.

RESULTS

432,306 mRNA molecules were sampled from a leaf transcriptome of a single B. bituminosa plant using Roche 454 pyrosequencing, resulting in an average read length of 345 bp (149.1 Mbp in total). Sequences were assembled into 3,838 isotigs/contigs representing putatively unique gene transcripts. Gene ontology descriptors were identified for 3,419 sequences. Raw sequence reads containing simple sequence repeat (SSR) motifs were identified, and 240 primer pairs flanking these motifs were designed. Of 87 primer pairs developed this way, 75 (86.2%) successfully amplified primarily single fragments by PCR. Fragment analysis using 20 primer pairs in 79 accessions of B. bituminosa detected 130 alleles at 21 SSR loci. Genetic diversity analyses confirmed that variation at these SSR loci accurately reflected known taxonomic relationships in original collections of B. bituminosa and provided additional evidence that a division of the botanical variety bituminosa into two according to geographical origin (Mediterranean region and Canary Islands) may be appropriate. Evidence of cross-pollination was also found between botanical varieties within a B. bituminosa breeding programme.

CONCLUSIONS

B. bituminosa can no longer be considered a genomic orphan species, having now a large (albeit incomplete) repertoire of expressed gene sequences that can serve as a resource for future genetic studies. This experimental approach was effective in developing codominant and polymorphic SSR markers for application in diverse genetic studies. These markers have already given new insight into genetic variation in B. bituminosa, providing evidence that a division of the botanical variety bituminosa may be appropriate. This approach is commended to those seeking to develop useful markers for genomic orphan species.

Production of viable male unreduced gametes in Brassica interspecific hybrids is genotype specific and stimulated by cold temperatures

BACKGROUND

Unreduced gametes (gametes with the somatic chromosome number) may provide a pathway for evolutionary speciation via allopolyploid formation. We evaluated the effect of genotype and temperature on male unreduced gamete formation in Brassica allotetraploids and their interspecific hybrids. The frequency of unreduced gametes post-meiosis was estimated in sporads from the frequency of dyads or giant tetrads, and in pollen from the frequency of viable giant pollen compared with viable normal pollen. Giant tetrads were twice the volume of normal tetrads, and presumably resulted from pre-meiotic doubling of chromosome number. Giant pollen was defined as pollen with more than 1.5 × normal diameter, under the assumption that the doubling of DNA content in unreduced gametes would approximately double the pollen cell volume. The effect of genotype was assessed in five B. napus, two B. carinata and one B. juncea parents and in 13 interspecific hybrid combinations. The effect of temperature was assessed in a subset of genotypes in hot (day/night 30°C/20°C), warm (25°C/15°C), cool (18°C/13°C) and cold (10°C/5°C) treatments.

RESULTS

Based on estimates at the sporad stage, some interspecific hybrid genotypes produced unreduced gametes (range 0.06 to 3.29%) at more than an order of magnitude higher frequency than in the parents (range 0.00% to 0.11%). In nine hybrids that produced viable mature pollen, the frequency of viable giant pollen (range 0.2% to 33.5%) was much greater than in the parents (range 0.0% to 0.4%). Giant pollen, most likely formed from unreduced gametes, was more viable than normal pollen in hybrids. Two B. napus × B. carinata hybrids produced 9% and 23% unreduced gametes based on post-meiotic sporad observations in the cold temperature treatment, which was more than two orders of magnitude higher than in the parents.

CONCLUSIONS

These results demonstrate that sources of unreduced gametes, required for the triploid bridge hypothesis of allopolyploid evolution, are readily available in some Brassica interspecific hybrid genotypes, especially at cold temperatures.

Characterization of Brassica nigra collections using simple sequence repeat markers reveals distinct groups associated with geographical location, and frequent mislabelling of species identity

Genetic diversity of 180 Brassica nigra (L.) Kochgenotypes from 60 different accessions was evaluated using 15 simple sequence repeat markers with known locations on the Brassica A, B, and C genomes. Two lines each from Brassica juncea (L.) Czern and Brassica carinata Braunwere also included as comparator species. A total of 218 high quality alleles were used to generate a genetic distance matrix, and clustering and multidimensional scaling analyses were used to investigate genetic relationships among the accessions. Accessions from the same country of origin tended to cluster together. Surprisingly, 13 accessions declared to be B. nigra had A- and B-genome alleles and morphology consistent with them being B. juncea, which was supported by their positioning near B. juncea in the cluster analysis. Two B. nigra accessions possessed alleles associated more closely with the A genome than the B genome, and these may be Brassica rapa L. accessions. One B. nigra accession had B- and C-genome alleles and morphology consistent with it being B. carinata. The remaining 44 accessions (73%) appeared to be truly B. nigra and formed morphologically and genetically distinct groups associated with country or region of origin, notably Ethiopia, Israel, India, and Europe. Most B. nigra accessions were highly heterozygous, consistent with their obligate outcrossing habit. This study demonstrated the value of using molecular markers with known genome locations (in this case, in the Brassica A, B, and C genomes) to confirm species identity in families such as Brassicaceae where species identification based solely on morphological characters is difficult.

Genotypic effects on the frequency of homoeologous and homologous recombination in Brassica napus × B. carinata hybrids

We investigated the influence of genotype on homoeologous and homologous recombination frequency in eight different Brassica napus (AAC(n)C(n)) × B. carinata (BBC(c)C(c)) interspecific hybrids (genome composition C(n)C(c)AB). Meiotic recombination events were assessed through microsatellite marker analysis of 67 unreduced microspore-derived progeny. Thirty-four microsatellite markers amplified 83 A-, B-, C(n)- and C(c)-genome alleles at 64 loci, of which a subset of seven markers amplifying 26 alleles could be used to determine allele copy number. Hybrid genotypes varied significantly in loss of A- and B-genome alleles (P < 0.0001), which ranged from 6 to 22% between hybrid progeny sets. Allele copy number analysis revealed 19 A-C, 3 A-B and 10 B-C duplication/deletion events attributed to homoeologous recombination. Additionally, 55 deletions and 19 duplications without an accompanying dosage change in homoeologous alleles were detected. Hybrid progeny sets varied in observed frequencies of loss, gain and exchange of alleles across the A and B genomes as well as in the diploid C genome. Self-fertility in hybrid progeny decreased as the loss of B-genome loci (but not A-genome loci) increased. Hybrid genotypes with high levels of homologous and homoeologous exchange may be exploited for genetic introgressions between B. carinata and B. napus (canola), and those with low levels may be used to develop stable synthetic Brassica allopolyploids.

The mosaic of ancestral karyotype blocks in the Sinapis alba L. genome

The organisation of the Sinapis alba genome, comprising 12 linkage groups (n = 12), was compared with the Brassicaceae ancestral karyotype (AK) genomic blocks previously described in other crucifer species. Most of the S. alba genome falls into conserved triplicated genomic blocks that closely match the AK-defined genomic blocks found in other crucifer species including the A, B, and C genomes of closely related Brassica species. In one instance, an S. alba linkage group (S05) was completely collinear with one AK chromosome (AK1), the first time this has been observed in a member of the Brassiceae tribe. However, as observed for other members of the Brassiceae tribe, ancestral genomic blocks were fragmented in the S. alba genome, supporting previously reported comparative chromosome painting describing rearrangements of the AK karyotype prior to the divergence of the Brassiceae from other crucifers. The presented data also refute previous phylogenetic reports that suggest S. alba was more closely related to Brassica nigra (B genome) than to B. rapa (A genome) and B. oleracea (C genome). A comparison of the S. alba and Arabidopsis thaliana genomes revealed many regions of conserved gene order, which will facilitate access to the rich genomic resources available in the model species A. thaliana for genetic research in the less well-resourced crop species S. alba.

An efficient high-throughput flow cytometric method for estimating DNA ploidy level in plants

We present an efficient high-throughput flow cytometric method that builds on previously published methods and permits rapid ploidy discrimination in plants. By using Brassica napus L. microspore-derived plants as an example, we describe how 192 leaf tissue samples may be processed and analyzed comfortably by one operator in 6 h from tissue sampling to ploidy determination. The technique involves placing young leaf samples in two 96-well racks, using a bead-beating procedure to release nuclei into a lysis solution, filtering the samples on 96-well filter plates, staining with propidium iodide, and then rapidly estimating DNA ploidy using a plate loader on a BD FACS-Canto II flow cytometer. Throughout the sample preparation process, multichannel pipetting allows faster and less error-prone sample handling. In two 96-well plates of samples, the histogram peaks of DNA content from flow cytometry were wellresolved in 189 of 192 samples tested (98.4%), with CV values ranging from 2.98% to 6.20% with an average CV of 4.35% (SD = 0.68%). This new method is useful in doubled haploid plant breeding programs where early discrimination of haploid and doubled haploid (i.e., diploid) plantlets can confer significantly improved operational efficiencies. We discuss how this method could be further refined including adapting the method to robotic sample processing.

Aligning a new reference genetic map of Lupinus angustifolius with the genome sequence of the model legume, Lotus japonicus

We have developed a dense reference genetic map of Lupinus angustifolius (2n = 40) based on a set of 106 publicly available recombinant inbred lines derived from a cross between domesticated and wild parental lines. The map comprised 1090 loci in 20 linkage groups and three small clusters, drawing together data from several previous mapping publications plus almost 200 new markers, of which 63 were gene-based markers. A total of 171 mainly gene-based, sequence-tagged site loci served as bridging points for comparing the Lu. angustifolius genome with the genome sequence of the model legume, Lotus japonicus via BLASTn homology searching. Comparative analysis indicated that the genomes of Lu. angustifolius and Lo. japonicus are highly diverged structurally but with significant regions of conserved synteny including the region of the Lu. angustifolius genome containing the pod-shatter resistance gene, lentus. We discuss the potential of synteny analysis for identifying candidate genes for domestication traits in Lu. angustifolius and in improving our understanding of Fabaceae genome evolution.

Microspore culture preferentially selects unreduced (2n) gametes from an interspecific hybrid of Brassica napus L. x Brassica carinata Braun

We analysed the products of male meiosis in microspore-derived progeny from a Brassica napus (AAC(n)C(n)) x Brassica carinata (BBC(c)C(c)) interspecific hybrid (ABC(n)C(c)). Genotyping at 102 microsatellite marker loci and nuclear DNA contents provided strong evidence that 26 of the 28 progeny (93%) were derived from unreduced (2n) gametes. The high level of C(n)C(c) marker heterozygosity, and parallel spindles at Anaphase II in the ABC(n)C(c) hybrid, indicated that unreduced gametes were formed by first division restitution. The frequency of dyads at the tetrad stage of pollen development (2.6%) suggested that unreduced gametes were preferentially selected in microspore culture. Segregation of marker alleles in the microspore-derived progeny was consistent with homologous recombination between C(n) and C(c) chromosomes and homoeologous recombination involving A-, B- and C-genome chromosomes during meiosis in the ABC(n)C(c) hybrid. We discuss the potential for using microspore culture of unreduced gametes in interspecific hybrids to map Brassica centromeres through half-tetrad analysis.

Twinned microspore-derived embryos of canola (Brassica napus L.) are genetically identical

Microspore culture is used extensively in several crop species to generate diverse populations of homozygous, doubled haploid lines for breeding and genetic analyses. In our canola (Brassica napus L.) doubled haploid breeding programme we regularly observe conjoined microspore-derived embryos, most commonly twins, joined either at the base of the hypocotyl or along the length of the hypocotyl axis. The aim of this study was to determine if twinned embryos were genetically identical or non-identical in order to gauge their value for breeding and linkage analysis. Microsatellite marker fingerprinting of 12 pairs of twinned embryos produced by microspore culture of heterozygous F(1) lines revealed that pairs of twins were genetically identical. Based on this finding, we recommend breeders and geneticists using microspore culture technology to retain only one embryo from each pair of twinned embryos.

Divergent patterns of allelic diversity from similar origins: the case of oilseed rape (Brassica napus L.) in China and Australia

Oilseed rape (Brassica napus) in Australia and China have similar origins, with introductions from Europe, Canada, and Japan in the mid 20th century, and there has been some interchange of germplasm between China and Australia since that time. Allelic diversity of 72 B. napus genotypes representing contemporary germplasm in Australia and China, including samples from India, Europe, and Canada, was characterized by 55 polymorphic simple sequence repeat (SSR) markers spanning the entire B. napus genome. Hierarchical clustering and two-dimensional multidimensional scaling identified a Chinese group (China-1) that was separated from "mixed group" of Australian, Chinese (China-2), European, and Canadian lines. A small group from India was distinctly separated from all other B. napus genotypes. Chinese genotypes, especially in the China-1 group, have inherited unique alleles from interspecific crossing, primarily with B. rapa, and the China-2 group has many alleles in common with Australian genotypes. The concept of "private alleles" is introduced to describe both the greater genetic diversity and the genetic distinctiveness of Chinese germplasm, compared with Australian germplasm, after 50 years of breeding from similar origins.

The first genetic and comparative map of white lupin (Lupinus albus L.): identification of QTLs for anthracnose resistance and flowering time, and a locus for alkaloid content

We report the first genetic linkage map of white lupin (Lupinus albus L.). An F8 recombinant inbred line population developed from Kiev mutant x P27174 was mapped with 220 amplified fragment length polymorphism and 105 gene-based markers. The genetic map consists of 28 main linkage groups (LGs) that varied in length from 22.7 cM to 246.5 cM and spanned a total length of 2951 cM. There were seven additional pairs and 15 unlinked markers, and 12.8% of markers showed segregation distortion at P < 0.05. Syntenic relationships between Medicago truncatula and L. albus were complex. Forty-five orthologous markers that mapped between M. truncatula and L. albus identified 17 small syntenic blocks, and each M. truncatula chromosome aligned to between one and six syntenic blocks in L. albus. Genetic mapping of three important traits: anthracnose resistance, flowering time, and alkaloid content allowed loci governing these traits to be defined. Two quantitative trait loci (QTLs) with significant effects were identified for anthracnose resistance on LG4 and LG17, and two QTLs were detected for flowering time on the top of LG1 and LG3. Alkaloid content was mapped as a Mendelian trait to LG11.