Kushen root, from the woody legume Sophora flavescens, is a traditional Chinese medicine that is a key ingredient in several promising cancer treatments. This activity is attributed in part to two quinolizidine alkaloids (QAs), oxymatrine and matrine, that have a variety of therapeutic activities in vitro. Genetic selection is needed to adapt S. flavescens for cultivation and to improve productivity and product quality. Genetic diversity of S. flavescens was investigated using genotyping-by-sequencing (GBS) on 85 plants grown from seeds collected from 9 provinces of China. DArTSeq provided over 10,000 single nucleotide polymorphism (SNP) markers, 1636 of which were used in phylogenetic analysis to reveal clear regional differences for S. flavescens. One accession from each region was selected for PCR-sequencing to identify gene-specific SNPs in the first two QA pathway genes, lysine decarboxylase (LDC) and copper amine oxidase (CAO). To obtain SfCAO sequence for primer design we used a targeted transcript capture and assembly strategy using publicly available RNA sequencing data. Partial gene sequence analysis of SfCAO revealed two recently duplicated genes, SfCAO1 and SfCAO2, in contrast to the single gene found in the QA-producing legume Lupinus angustifolius. We demonstrate high efficiency converting SNPs to Kompetitive Allele Specific PCR (KASP) markers developing 27 new KASP markers, 17 from DArTSeq data, 7 for SfLDC, and 3 for SfCAO1. To complement this genetic diversity analysis a field trial site has been established in South Australia, providing access to diverse S. flavescens material for morphological, transcriptomic, and QA metabolite analysis. Analysis of dissected flower buds revealed that anthesis occurs before buds fully open suggesting a potential for S. flavescens to be an inbreeding species, however this is not supported by the relatively high level of heterozygosity observed. Two plants from the field trial site were analysed by quantitative real-time PCR and levels of matrine and oxymatrine were assessed in a variety of tissues. We are now in a strong position to select diverse plants for crosses to accelerate the process of genetic selection needed to adapt kushen to cultivation and improve productivity and product quality.

Genetic loci, particularly those with an effect in the independent panel, could be utilised to further reduce LMA expression when used with favourable combinations of genes known to affect LMA. Late maturity α-amylase (LMA) is a grain quality defect involving elevated α-amylase within the aleurone of wheat (Triticum aestivum L.) grains. The genes known to affect expression are the reduced height genes Rht-B1 (chromosome 4B) and Rht-D1 (chromosome 4D), and an ent-copalyl diphosphate synthase gene (LMA-1) on chromosome 7B. Other minor effect loci have been reported, but these are poorly characterised and further genetic understanding is needed. In this study, twelve F4-derived populations were created through single seed descent, genotyped and evaluated for LMA. LMA-1 haplotype C and the Rht-D1b allele substantially reduced LMA expression. The alternative dwarfing genes Rht13 and Rht18 had no significant effect on LMA expression. Additional quantitative trait loci (QTL) were mapped at 16 positions in the wheat genome. Effects on LMA expression were detected for four of these QTL in a large independent panel of Australian wheat lines. The QTL detected in mapping populations and confirmed in the large independent panel provide further opportunity for selection against LMA, especially if combined with Rht-D1b and/or favourable haplotypes of LMA-1.

MicroRNAs (miRNAs) are small non-coding RNAs that play critical roles in plant development and abiotic stress responses. The miRNA transcriptome (miRNAome) under water deficit stress has been investigated in many plant species, but is poorly characterised in durum wheat (Triticum turgidum L. ssp. durum). Water stress during early reproductive stages can result in significant yield loss in durum wheat and this study describes genotypic differences in the miRNAome between water deficit tolerant and sensitive durum genotypes. Small RNA libraries (96 in total) were constructed from flag leaf and developing head tissues of four durum genotypes, with or without water stress to identify differentially abundant miRNAs. Illumina sequencing detected 110 conserved miRNAs and 159 novel candidate miRNA hairpins with 66 conserved miRNAs and five novel miRNA hairpins differentially abundant under water deficit stress. Ten miRNAs (seven conserved, three novel) were validated through qPCR. Several conserved and novel miRNAs showed unambiguous inverted regulatory profiles between the durum genotypes. Several miRNAs also showed differential abundance between two tissue types regardless of treatment. Predicted mRNA targets (130) of four novel durum miRNAs were characterised using Gene Ontology (GO) which revealed functions common to stress responses and plant development. Negative correlation was observed between several target genes and the corresponding miRNA under water stress. For the first time, we present a comprehensive study of the durum miRNAome under water deficit stress. The identification of differentially abundant miRNAs provides molecular evidence that miRNAs are potential determinants of water stress tolerance in durum wheat. GO analysis of predicted targets contributes to the understanding of genotypic physiological responses leading to stress tolerance capacity. Further functional analysis of specific stress responsive miRNAs and their interaction with targets is ongoing and will assist in developing future durum wheat varieties with enhanced water deficit stress tolerance.

The cumulative action of combinations of alleles at several loci on the wheat genome is associated with different levels of resistance to late maturity α-amylase in bread wheat. Resistance to late maturity α-amylase (LMA) in bread wheat (Triticum aestivum L.) involves a complex interaction between the genotype and the environment. Unfortunately, the incidence and severity of LMA expression is difficult to predict and once the trait has been triggered an unacceptably low falling number, high grain α-amylase may be the inevitable consequence. Wheat varieties with different levels of resistance to LMA have been identified but whilst some genetic loci have been reported, the mechanisms involved in resistance and the interaction between resistance loci requires further research. This investigation was focused on mapping resistance loci in populations derived by inter-crossing resistant wheat varieties or crossing resistant lines with a very susceptible line and then mapping quantitative trait loci. In addition to the previously reported locus on chromosome 7B for which a candidate gene has been proposed, loci were mapped on chromosomes 1B, 2A, 2B, 3A, 3B, 4A, 6A and 7D. These loci have limited effects on their own but have a cumulative effect in combination with each other. Further research will be required to determine the nature of the causal genes at these loci, to develop diagnostic markers and determine how the genes fit into the pathway that leads to the induction of α-AMY1 transcription in the aleurone of developing wheat grains. Depending on the target environmental conditions, different combinations of alleles may be required to achieve a low risk of LMA expression.

In many plant species, self-incompatibility systems limit self-pollination and mating among relatives. This helps maintain genetic diversity in natural populations but imposes constraints in agriculture and plant breeding. In almond [Prunus dulcis (Mill.) D.A. Webb], the specificity of self-incompatibility is mainly determined by stylar ribonuclease (S-RNase) and S-haplotype-specific F-box (SFB) proteins, both encoded within a complex locus, S. Prior to this research, a nearly complete sequence was available for one S-locus haplotype. Here, we report complete sequences for four haplotypes and partial sequences for 11 haplotypes. Haplotypes vary in sequences of genes (particularly S-RNase and SFB), distances between genes and numbers and positions of long terminal repeat transposons. Haplotype variation outside of the S-RNase and SFB genes may help maintain functionally important associations between S-RNase and SFB alleles. Fluorescence-based assays were developed to distinguish among some S-RNase alleles. With three-dimensional modelling of five S-RNase proteins, conserved active sites were identified and variation was observed in electrostatic potential and in the numbers, characteristics and positions of secondary structural elements, loop anchoring points and glycosylation sites. A hypervariable region on the protein surface and differences in the number, location and types of glycosylation sites may contribute to determining S-RNase specificity.

Novel QTL for salinity tolerance traits have been detected using non-destructive and destructive phenotyping in bread wheat and were shown to be linked to improvements in yield in saline fields. Soil salinity is a major limitation to cereal production. Breeding new salt-tolerant cultivars has the potential to improve cereal crop yields. In this study, a doubled haploid bread wheat mapping population, derived from the bi-parental cross of Excalibur × Kukri, was grown in a glasshouse under control and salinity treatments and evaluated using high-throughput non-destructive imaging technology. Quantitative trait locus (QTL) analysis of this population detected multiple QTL under salt and control treatments. Of these, six QTL were detected in the salt treatment including one for maintenance of shoot growth under salinity (QG(1-5).asl-7A), one for leaf Na+ exclusion (QNa.asl-7A) and four for leaf K+ accumulation (QK.asl-2B.1, QK.asl-2B.2, QK.asl-5A and QK:Na.asl-6A). The beneficial allele for QG(1-5).asl-7A (the maintenance of shoot growth under salinity) was present in six out of 44 mainly Australian bread and durum wheat cultivars. The effect of each QTL allele on grain yield was tested in a range of salinity concentrations at three field sites across 2 years. In six out of nine field trials with different levels of salinity stress, lines with alleles for Na+ exclusion and/or K+ maintenance at three QTL (QNa.asl-7A, QK.asl-2B.2 and QK:Na.asl-6A) excluded more Na+ or accumulated more K+ compared to lines without these alleles. Importantly, the QK.asl-2B.2 allele for higher K+ accumulation was found to be associated with higher grain yield at all field sites. Several alleles at other QTL were associated with higher grain yields at selected field sites.

Wheat landraces, wild relatives and other 'exotic' accessions are important sources of new favorable alleles. The use of those exotic alleles is facilitated by having access to information on the association of specific genomic regions with desirable traits. Here, we conducted a genome-wide association study (GWAS) using a wheat panel that includes landraces, synthetic hexaploids and other exotic wheat accessions to identify loci that contribute to increases in grain yield in southern Australia. The 568 accessions were grown in the field during the 2014 and 2015 seasons and measured for plant height, maturity, spike length, spike number, grain yield, plant biomass, HI and TGW. We used the 90K SNP array and two GWAS approaches (GAPIT and QTCAT) to identify loci associated with the different traits. We identified 17 loci with GAPIT and 25 with QTCAT. Ten of these loci were associated with known genes that are routinely employed in marker assisted selection such as Ppd-D1 for maturity and Rht-D1 for plant height and seven of those were detected with both methods. We identified one locus for yield per se in 2014 on chromosome 6B with QTCAT and three in 2015, on chromosomes 4B and 5A with GAPIT and 6B with QTCAT. The 6B loci corresponded to the same region in both years. The favorable haplotypes for yield at the 5A and 6B loci are widespread in Australian accessions with 112 out of 153 carrying the favorable haplotype at the 5A locus and 136 out of 146 carrying the favorable haplotype at the 6A locus, while the favorable haplotype at 4B is only present in 65 out of 149 Australian accessions. The low number of yield QTL in our study corroborate with other GWAS for yield in wheat, where most of the identified loci have very small effects.

The cellulose synthase-like gene HvCslF6, which is essential for (1,3;1,4)-β-glucan biosynthesis in barley, collocates with quantitative trait loci (QTL) for grain (1,3;1,4)-β-glucan concentration in several populations, including CDC Bold × TR251. Here, an alanine-to-threonine substitution (caused by the only non-synonymous difference between the CDC Bold and TR251 HvCslF6 alleles) was mapped to a position within HvCSLF6 that seems unlikely to affect enzyme stability or function. Consistent with this, transient expression of full-length HvCslF6 cDNAs from CDC Bold and TR251 in Nicotianabenthamiana led to accumulation of similar amounts of (1,3;1,4)-β-glucan accumulation. Monitoring of HvCslF6 transcripts throughout grain development revealed a significant difference late in grain development (more than 30 days after pollination), with TR251 [the parent with higher grain (1,3;1,4)-β-glucan] exhibiting higher transcript levels than CDC Bold. A similar difference was observed between Beka and Logan, the parents of another population in which a QTL had been mapped in the HvCslF6 region. Sequencing of a putative promoter region of HvCslF6 revealed numerous polymorphisms between CDC Bold and TR251, but none between Beka and Logan. While the results of this work indicate that naturally occurring quantitative differences in (1,3;1,4)-β-glucan accumulation may be due to cis-regulated differences in HvCslF6 expression, these could not be attributed to any specific DNA sequence polymorphism. Nevertheless, information on HvCslF6 sequence polymorphism was used to develop molecular markers that could be used in barley breeding to select for the desired [low or high (1,3;1,4)-β-glucan] allele of the QTL.

The cereal cyst nematode resistance locus Rha2 was mapped to a 978 kbp region on the long arm of barley chromosome 2H. Three candidate genes are discussed. The cereal cyst nematode (CCN) Heterodera avenae is a soil-borne obligate parasite that can cause severe damage to cereals. This research involved fine mapping of Rha2, a CCN resistance locus on chromosome 2H of barley. Rha2 was previously mapped relative to restriction fragment length polymorphisms (RFLPs) in two mapping populations. Anchoring of flanking RFLP clone sequences to the barley genome assembly defined an interval of 5077 kbp. Genotyping-by-sequencing of resistant and susceptible materials led to the discovery of potentially useful single nucleotide polymorphisms (SNPs). Assays were designed for these SNPs and applied to mapping populations. This narrowed the region of interest to 3966 kbp. Further fine mapping was pursued by crossing and backcrossing the resistant cultivar Sloop SA to its susceptible ancestor Sloop. Evaluation of F2 progeny confirmed that the resistance segregates as a single dominant gene. Genotyping of 9003 BC2F2 progeny identified recombinants. Evaluation of recombinant BC2F3 progeny narrowed the region of interest to 978 kbp. Two of the SNPs within this region proved to be diagnostic of CCN resistance across a wide range of barley germplasm. Fluorescence-based and gel-based assays were developed for these SNPs for use in marker-assisted selection. Within the candidate region of the reference genome, there are nine high-confidence predicted genes. Three of these, one that encodes RAR1 (a cysteine- and histidine-rich domain-containing protein), one that is predicted to encode an acetylglutamate kinase and one that is predicted to encode a tonoplast intrinsic protein, are discussed as candidate genes for CCN resistance.

Cyst nematodes induce host-plant root cells to form syncytia from which the nematodes feed. Comprehensive histological investigation of these feeding sites is complicated by their variable shape and their positions deep within root tissue. Using tissue clearing and confocal microscopy, we examined thick (up to 150 μm) sections of wheat roots infected by cereal cyst nematodes (Heterodera avenae). This approach provided clear views of feeding sites and surrounding tissues, with resolution sufficient to reveal spatial relationships among nematodes, syncytia and host vascular tissues at the cellular level. Regions of metaxylem vessels near syncytia were found to have deviated from classical developmental patterns. Xylem vessel elements in these regions had failed to elongate but had undergone radial expansion, becoming short and plump rather than long and cylindrical. Further investigation revealed that vessel elements cease to elongate shortly after infection and that they later experience delays in secondary thickening (lignification) of their outer cell walls. Some of these elements were eventually incorporated into syncytial feeding sites. By interfering with a developmental program that normally leads to programmed cell death, H. avenae may permit xylem vessel elements to remain alive for later exploitation by the parasite.

Molecular markers and knowledge of traits associated with heat tolerance are likely to provide breeders with a more efficient means of selecting wheat varieties able to maintain grain size after heat waves during early grain filling.

In crop plant genetics, linkage maps provide the basis for the mapping of loci that affect important traits and for the selection of markers to be applied in crop improvement. In outcrossing species such as almond (Prunus dulcis Mill. D. A. Webb), application of a double pseudotestcross mapping approach to the F1 progeny of a biparental cross leads to the construction of a linkage map for each parent. Here, we report on the application of genotyping by sequencing to discover and map single nucleotide polymorphisms in the almond cultivars "Nonpareil" and "Lauranne." Allele-specific marker assays were developed for 309 tag pairs. Application of these assays to 231 Nonpareil × Lauranne F1 progeny provided robust linkage maps for each parent. Analysis of phenotypic data for shell hardness demonstrated the utility of these maps for quantitative trait locus mapping. Comparison of these maps to the peach genome assembly confirmed high synteny and collinearity between the peach and almond genomes. The marker assays were applied to progeny from several other Nonpareil crosses, providing the basis for a composite linkage map of Nonpareil. Applications of the assays to a panel of almond clones and a panel of rootstocks used for almond production demonstrated the broad applicability of the markers and provide subsets of markers that could be used to discriminate among accessions. The sequence-based linkage maps and single nucleotide polymorphism assays presented here could be useful resources for the genetic analysis and genetic improvement of almond.

QTL for tan spot resistance were mapped on wheat chromosomes 1A and 2A. Lines were developed with resistance alleles at these loci and at the tsn1 locus on chromosome 5B. These lines expressed significantly higher resistance than the parent with tsn1 only. Tan spot (syn. yellow spot and yellow leaf spot) caused by Pyrenophora tritici-repentis is an important foliar disease of wheat in Australia. Few resistance genes have been mapped in Australian germplasm and only one, known as tsn1 located on chromosome 5B, is known in Australian breeding programs. This gene confers insensitivity to the fungal effector ToxA. The main aim of this study was to map novel resistance loci in two populations: Calingiri/Wyalkatchem, which is fixed for the ToxA-insensitivity allele tsn1, and IGW2574/Annuello, which is fixed for the ToxA-sensitivity allele Tsn1. A second aim was to combine new loci with tsn1 to develop lines with improved resistance. Tan spot severity was evaluated at various growth stages and in multiple environments. Symptom severity traits exhibited quantitative variation. The most significant quantitative trait loci (QTL) were detected on chromosomes 2A and 1A. The QTL on 2A explained up to 29.2% of the genotypic variation in the Calingiri/Wyalkatchem population with the resistance allele contributed by Wyalkatchem. The QTL on 1A explained up to 28.1% of the genotypic variation in the IGW2574/Annuello population with the resistance allele contributed by Annuello. The resistance alleles at both QTL were successfully combined with tsn1 to develop lines that express significantly better resistance at both seedling and adult plant stages than Calingiri which has tsn1 only.

The Rlnn1 locus, which resides on chromosome 7A of bread wheat (Triticum aestivum L.) confers moderate resistance against the root lesion nematode Pratylenchus neglectus. Prior to this research, the exact linkage relationships of Rlnn1 with other loci on chromosome 7A were not clear and there were no simple codominant markers available for selection of Rlnn1 in wheat breeding. The objectives of the research reported here were to (1) develop an improved genetic map of the Rlnn1 region of chromosome 7A and (2) develop molecular markers that could be used in marker-assisted selection to improve resistance of wheat against P. neglectus.

Para rubber trees (Hevea brasiliensis) are the largest major source of natural rubber in the world. Its major pathogens are Phytophthora spp., Corynespora cassiicola, and Colletotrichum spp. A rubber diversity panel of 116 clones using over 12,000 single nucleotide polymorphisms (SNPs) from DArTSeq genotyping revealed clear phylogenetic differences in clones that originated from different geographical regions of the world. An integrated linkage map constructed with an F1 progeny of 86 from an interspecific cross between H. brasiliensis and H. benthamiana using 23,978 markers [10,323 SNPs and 13,655 SilicoDArTs] spanned 3947.83 cM with 0.83 cM average marker-interval. The genome scaffolds that were anchored to the linkage map, covering 1.44 Gb of H. brasiliensis reference genome, revealed a high level of collinearity between the genetic map and reference genome. Association analysis identified 12 SNPs significantly associated with the resistance against Phytophthora, Corynespora, and Colletotrichum in six linkage groups: 2, 6, 12, 14, 17, and 18. Kompetitive Allele-Specific PCR marker assays were developed for those 12 SNPs, screened with 178 individuals, and detected clear separation between two genotypes. Within the proximity to those SNPs, 41 potentially key genes that have previously been reported to associate with plant disease resistance were predicted with high confidence.

Many research and diagnostic applications rely upon the assay of individual single nucleotide polymorphisms (SNPs). Thus, methods to improve the speed and efficiency for single-marker SNP genotyping are highly desirable. Here, we describe the method of temperature-switch PCR (TSP), a biphasic four-primer PCR system with a universal primer design that permits amplification of the target locus in the first phase of thermal cycling before switching to the detection of the alleles. TSP can simplify assay design for a range of commonly used single-marker SNP genotyping methods, and reduce the requirement for individual assay optimization and operator expertise in the deployment of SNP assays.

Cellulose is a fundamentally important component of cell walls of higher plants. It provides a scaffold that allows the development and growth of the plant to occur in an ordered fashion. Cellulose also provides mechanical strength, which is crucial for both normal development and to enable the plant to withstand both abiotic and biotic stresses. We quantified the cellulose concentration in the culm of 288 two--rowed and 288 six--rowed spring type barley accessions that were part of the USDA funded barley Coordinated Agricultural Project (CAP) program in the USA. When the population structure of these accessions was analysed we identified six distinct populations, four of which we considered to be comprised of a sufficient number of accessions to be suitable for genome-wide association studies (GWAS). These lines had been genotyped with 3072 SNPs so we combined the trait and genetic data to carry out GWAS. The analysis allowed us to identify regions of the genome containing significant associations between molecular markers and cellulose concentration data, including one region cross-validated in multiple populations. To identify candidate genes we assembled the gene content of these regions and used these to query a comprehensive RNA-seq based gene expression atlas. This provided us with gene annotations and associated expression data across multiple tissues, which allowed us to formulate a supported list of candidate genes that regulate cellulose biosynthesis. Several regions identified by our analysis contain genes that are co-expressed with cellulose synthase A (HvCesA) across a range of tissues and developmental stages. These genes are involved in both primary and secondary cell wall development. In addition, genes that have been previously linked with cellulose synthesis by biochemical methods, such as HvCOBRA, a gene of unknown function, were also associated with cellulose levels in the association panel. Our analyses provide new insights into the genes that contribute to cellulose content in cereal culms and to a greater understanding of the interactions between them.