HvSL1 and HvMADS16 promote stamen identity to restrict multiple ovary formation in barley.

Correct floral development is the result of a sophisticated balance of molecular cues. Floral mutants provide insight into the main genetic determinants that integrate these cues, as well as providing opportunities to assess functional variation across species. In this study, we characterize the barley (Hordeum vulgare) multiovary mutants mov2.g and mov1, and propose causative gene sequences: a C2H2 zinc-finger gene HvSL1 and a B-class gene HvMADS16, respectively. In the absence of HvSL1, florets lack stamens but exhibit functional supernumerary carpels, resulting in multiple grains per floret. Deletion of HvMADS16 in mov1 causes homeotic conversion of lodicules and stamens into bract-like organs and carpels that contain non-functional ovules. Based on developmental, genetic, and molecular data, we propose a model by which stamen specification in barley is defined by HvSL1 acting upstream of HvMADS16. The present work identifies strong conservation of stamen formation pathways with other cereals, but also reveals intriguing species-specific differences. The findings lay the foundation for a better understanding of floral architecture in Triticeae, a key target for crop improvement.

Comparative transcriptome analysis reveals molecular regulation of salt tolerance in two contrasting chickpea genotypes.

Salinity is a major abiotic stress that causes substantial agricultural losses worldwide. Chickpea (Cicer arietinum L.) is an important legume crop but is salt-sensitive. Previous physiological and genetic studies revealed the contrasting response of two desi chickpea varieties, salt-sensitive Rupali and salt-tolerant Genesis836, to salt stress. To understand the complex molecular regulation of salt tolerance mechanisms in these two chickpea genotypes, we examined the leaf transcriptome repertoire of Rupali and Genesis836 in control and salt-stressed conditions. Using linear models, we identified categories of differentially expressed genes (DEGs) describing the genotypic differences: salt-responsive DEGs in Rupali (1,604) and Genesis836 (1,751) with 907 and 1,054 DEGs unique to Rupali and Genesis836, respectively, salt responsive DEGs (3,376), genotype-dependent DEGs (4,170), and genotype-dependent salt-responsive DEGs (122). Functional DEG annotation revealed that the salt treatment affected genes involved in ion transport, osmotic adjustment, photosynthesis, energy generation, stress and hormone signalling, and regulatory pathways. Our results showed that while Genesis836 and Rupali have similar primary salt response mechanisms (common salt-responsive DEGs), their contrasting salt response is attributed to the differential expression of genes primarily involved in ion transport and photosynthesis. Interestingly, variant calling between the two genotypes identified SNPs/InDels in 768 Genesis836 and 701 Rupali salt-responsive DEGs with 1,741 variants identified in Genesis836 and 1,449 variants identified in Rupali. In addition, the presence of premature stop codons was detected in 35 genes in Rupali. This study provides valuable insights into the molecular regulation underpinning the physiological basis of salt tolerance in two chickpea genotypes and offers potential candidate genes for the improvement of salt tolerance in chickpeas.

Global transcriptome and gene co-expression network analyses reveal regulatory and non-additive effects of drought and heat stress in grapevine.

Despite frequent co-occurrence of drought and heat stress, the molecular mechanisms governing plant responses to these stresses in combination have not often been studied. This is particularly evident in non-model, perennial plants. We conducted large scale physiological and transcriptome analyses to identify genes and pathways associated with grapevine response to drought and/or heat stress during stress progression and recovery. We identified gene clusters with expression correlated to leaf temperature and water stress and five hub genes for the combined stress co-expression network. Several differentially expressed genes were common to the individual and combined stresses, but the majority were unique to the individual or combined stress treatments. These included heat-shock proteins, mitogen-activated kinases, sugar metabolizing enzymes, and transcription factors, while phenylpropanoid biosynthesis and histone modifying genes were unique to the combined stress treatment. Following physiological recovery, differentially expressed genes were found only in plants under heat stress, both alone and combined with drought. Taken collectively, our results suggest that the effect of the combined stress on physiology and gene expression is more severe than that of individual stresses, but not simply additive, and that epigenetic chromatin modifications may play an important role in grapevine responses to combined drought and heat stress.

Stress-Induced Volatile Emissions and Signalling in Inter-Plant Communication.

The sessile plant has developed mechanisms to survive the "rough and tumble" of its natural surroundings, aided by its evolved innate immune system. Precise perception and rapid response to stress stimuli confer a fitness edge to the plant against its competitors, guaranteeing greater chances of survival and productivity. Plants can "eavesdrop" on volatile chemical cues from their stressed neighbours and have adapted to use these airborne signals to prepare for impending danger without having to experience the actual stress themselves. The role of volatile organic compounds (VOCs) in plant-plant communication has gained significant attention over the past decade, particularly with regard to the potential of VOCs to prime non-stressed plants for more robust defence responses to future stress challenges. The ecological relevance of such interactions under various environmental stresses has been much debated, and there is a nascent understanding of the mechanisms involved. This review discusses the significance of VOC-mediated inter-plant interactions under both biotic and abiotic stresses and highlights the potential to manipulate outcomes in agricultural systems for sustainable crop protection via enhanced defence. The need to integrate physiological, biochemical, and molecular approaches in understanding the underlying mechanisms and signalling pathways involved in volatile signalling is emphasised.

Nested association mapping-based GWAS for grain yield and related traits in wheat grown under diverse Australian environments.

KEY MESSAGE: Utilising a nested association mapping (NAM) population-based GWAS, 98 stable marker-trait associations with 127 alleles unique to the exotic parents were detected for grain yield and related traits in wheat. Grain yield, thousand-grain weight, screenings and hectolitre weight are important wheat yield traits. An understanding of their genetic basis is crucial for improving grain yield in breeding programmes. Nested association mapping (NAM) populations are useful resources for the dissection of the genetic basis of complex traits such as grain yield and related traits in wheat. Coupled with phenotypic data collected from multiple environments, NAM populations have the power to detect quantitative trait loci and their multiple alleles, providing germplasm that can be incorporated into breeding programmes. In this study, we evaluated a large-scale wheat NAM population with two recurrent parents in unbalanced trials in nine diverse Australian field environments over three years. By applying a single-stage factor analytical linear mixed model (FALMM) to the NAM multi-environment trials (MET) data and conducting a genome-wide association study (GWAS), we detected 98 stable marker-trait associations (MTAs) with their multiple alleles. 74 MTAs had 127 alleles that were derived from the exotic parents and were absent in either of the two recurrent parents. The exotic alleles had favourable effects on 46 MTAs of the 74 MTAs, for grain yield, thousand-grain weight, screenings and hectolitre weight. Two NAM RILs with consistently high yield in multiple environments were also identified, highlighting the potential of the NAM population in supporting plant breeding through provision of germplasm that can be readily incorporated into breeding programmes. The identified beneficial exotic alleles introgressed into the NAM population provide potential target alleles for the genetic improvement of wheat and further studies aimed at pinpointing the underlying genes.

The complete genome sequence of elite bread wheat cultivar, "Sonmez".

High-yielding crop varieties will become critical in meeting the future food demand in the face of worsening weather extremes and threatening biotic stressors. The bread wheat cultivar Sonmez-2001 is a registered variety that is notable for its performance under low-irrigation conditions, which further improves upon irrigation. Additionally, Sonmez-2001 is resilient against certain biotic stressors, particularly soil-borne pathogens. Here, we provide a reference-guided whole genome sequence of Sonmez-2001, assembled into 21 chromosomes of the A, B and D genomes and totaling 13.3 gigabase-pairs in length. Additionally, a de novo assembly of an additional 1.05 gigabase-pairs was generated that represents either Sonmez-specific sequences or sequences that considerably diverged between Sonmez and Chinese Spring. Within this de novo assembly, we identified 35 gene models, of which 11 were high-confidence, that may contribute to the favorable traits of this high-performing variety. We identified up to 24 million sequence variants, of which up to 2.4% reside in coding sequences, that can be used to develop molecular markers that should be of immediate use to the cereal community.

Genetic control of compatibility in crosses between wheat and its wild or cultivated relatives.

In the recent years, the agricultural world has been progressing towards integrated crop protection, in the context of sustainable and reasoned agriculture to improve food security and quality, and to preserve the environment through reduced uses of water, pesticides, fungicides or fertilisers. For this purpose, one possible issue is to cross-elite varieties widely used in fields for crop productions with exotic or wild genetic resources in order to introduce new diversity for genes or alleles of agronomical interest to accelerate the development of new improved cultivars. However, crossing ability (or crossability) often depends on genetic background of the recipient varieties or of the donor, which hampers a larger use of wild resources in breeding programmes of many crops. In this review, we tried to provide a comprehensive summary of genetic factors controlling crossing ability between Triticeae species with a special focus on the crossability between wheat (Triticum aestivum L.) and rye (Secale cereale), which lead to the creation of Triticale (x Triticosecale Wittm.). We also discussed potential applications of newly identified genes or markers associated with crossability for accelerating wheat and Triticale improvement by application of modern genomics technologies in breeding programmes.

A targeted bioinformatics approach identifies highly variable cell surface proteins that are unique to Glomeromycotina.

Diversity in arbuscular mycorrhizal fungi (AMF) contributes to biodiversity and resilience in natural environments and healthy agricultural systems. Functional complementarity exists among species of AMF in symbiosis with their plant hosts, but the molecular basis of this is not known. We hypothesise this is in part due to the difficulties that current sequence assembly methodologies have assembling sequences for intrinsically disordered proteins (IDPs) due to their low sequence complexity. IDPs are potential candidates for functional complementarity because they often exist as extended (non-globular) proteins providing additional amino acids for molecular interactions. Rhizophagus irregularis arabinogalactan-protein-like proteins (AGLs) are small secreted IDPs with no known orthologues in AMF or other fungi. We developed a targeted bioinformatics approach to identify highly variable AGLs/IDPs in RNA-sequence datasets. The approach includes a modified multiple k-mer assembly approach (Oases) to identify candidate sequences, followed by targeted sequence capture and assembly (mirabait-mira). All AMF species analysed, including the ancestral family Paraglomeraceae, have small families of proteins rich in disorder promoting amino acids such as proline and glycine, or glycine and asparagine. Glycine- and asparagine-rich proteins also were found in Geosiphon pyriformis (an obligate symbiont of a cyanobacterium), from the same subphylum (Glomeromycotina) as AMF. The sequence diversity of AGLs likely translates to functional diversity, based on predicted physical properties of tandem repeats (elastic, amyloid, or interchangeable) and their broad pI ranges. We envisage that AGLs/IDPs could contribute to functional complementarity in AMF through processes such as self-recognition, retention of nutrients, soil stability, and water movement.

The genetics of vigour-related traits in chickpea (Cicer arietinum L.): insights from genomic data.

KEY MESSAGE: QTL controlling vigour and related traits were identified in a chickpea RIL population and validated in diverse sets of germplasm. Robust KASP markers were developed for marker-assisted selection. To understand the genetic constitution of vigour in chickpea (Cicer arietinum L.), genomic data from a bi-parental population and multiple diversity panels were used to identify QTL, sequence-level haplotypes and genetic markers associated with vigour-related traits in Australian environments. Using 182 Recombinant Inbred Lines (RILs) derived from a cross between two desi varieties, Rupali and Genesis836, vigour QTL independent of flowering time were identified on chromosomes (Ca) 1, 3 and 4 with genotypic variance explained (GVE) ranging from 7.1 to 28.8%. Haplotype analysis, association analysis and graphical genotyping of whole-genome re-sequencing data of two diversity panels consisting of Australian and Indian genotypes and an ICRISAT Chickpea Reference Set revealed a deletion in the FTa1-FTa2-FTc gene cluster of Ca3 significantly associated with vigour and flowering time. Across the RIL population and diversity panels, the impact of the deletion was consistent for vigour but not flowering time. Vigour-related QTL on Ca4 co-located with a QTL for seed size in Rupali/Genesis836 (GVE = 61.3%). Using SNPs from this region, we developed and validated gene-based KASP markers across different panels. Two markers were developed for a gene on Ca1, myo -inositol monophosphatase (CaIMP), previously proposed to control seed size, seed germination and seedling growth in chickpea. While associated with vigour in the diversity panels, neither the markers nor broader haplotype linked to CaIMP was polymorphic in Rupali/Genesis836. Importantly, vigour appears to be controlled by different sets of QTL across time and with components which are independent from phenology.

Proton-pumping pyrophosphatase homeolog expression is a dynamic trait in bread wheat (Triticum aestivum).

Proton-pumping pyrophosphatases (H+-PPases) have been shown to enhance biomass and yield. However, to date, there has been little work towards identify genes encoding H+-PPases in bread wheat (Triticum aestivum) (TaVPs) and limited knowledge on how the expression of these genes varies across different growth stages and tissue types. In this study, the IWGSC database was used to identify two novel TaVP genes, TaVP4 and TaVP5, and elucidate the complete homeolog sequences of the three known TaVP genes, bringing the total number of bread wheat TaVPs from 9 to 15. Gene expression levels of each TaVP homeolog were assessed using quantitative real-time PCR (qRT-PCR) in four diverse wheat varieties in terms of phenotypic traits related to high vacuolar pyrophosphatase expression. Homeolog expression was analyzed across multiple tissue types and developmental stages. Expression levels of the TaVP homeologs were found to vary significantly between varieties, tissues and plant developmental stages. During early development (Z10 and Z13), expressions of TaVP1 and TaVP2 homeologs were higher in shoot tissue than root tissue, with both shoot and root expression increasing in later developmental stages (Z22). TaVP2-D was expressed in all varieties and tissue types and was the most highly expressed homeolog at all developmental stages. Expression of the TaVP3 homeologs was restricted to developing grain (Z75), while TaVP4 homeolog expression was higher at Z22 than earlier developmental stages. Variation in TaVP4B was detected among varieties at Z22 and Z75, with Buck Atlantico (high biomass) and Scout (elite Australian cultivar) having the highest levels of expression. These findings offer a comprehensive overview of the bread wheat H+-PPase family and identify variation in TaVP homeolog expression that will be of use to improve the growth, yield, and abiotic stress tolerance of bread wheat.

HvLEAFY controls the early stages of floral organ specification and inhibits the formation of multiple ovaries in barley.

Transition to the reproductive phase, inflorescence formation and flower development are crucial elements that ensure maximum reproductive success in a plant's life cycle. To understand the regulatory mechanisms underlying correct flower development in barley (Hordeum vulgare), we characterized the multiovary 5 (mov5.o) mutant. This mutant develops abnormal flowers that exhibit mosaic floral organs typified by multiple carpels at the total or partial expense of stamens. Genetic mapping positioned mov5 on the long arm of chromosome 2H, incorporating a region that encodes HvLFY, the barley orthologue of LEAFY from Arabidopsis. Sequencing revealed that, in mov5.o plants, HvLFY contains a single amino acid substitution in a highly conserved proline residue. CRISPR-mediated knockout of HvLFY replicated the mov5.o phenotype, suggesting that HvLFYmov5 represents a loss of function allele. In heterologous assays, the HvLFYmov5 polymorphism influenced protein-protein interactions and affinity for a putative binding site in the promoter of HvMADS58, a C-class MADS-box gene. Moreover, molecular analysis indicated that HvLFY interacts with HvUFO and regulates the expression of floral homeotic genes including HvMADS2, HvMADS4 and HvMADS16. Other distinct changes in expression differ from those reported in the rice LFY mutants apo2/rfl, suggesting that LFY function in the grasses is modulated in a species-specific manner. This pathway provides a key entry point for the study of LFY function and multiple ovary formation in barley, as well as cereal species in general.

Development of an Australian Bread Wheat Nested Association Mapping Population, a New Genetic Diversity Resource for Breeding under Dry and Hot Climates.

Genetic diversity, knowledge of the genetic architecture of the traits of interest and efficient means of transferring the desired genetic diversity into the relevant genetic background are prerequisites for plant breeding. Exotic germplasm is a rich source of genetic diversity; however, they harbor undesirable traits that limit their suitability for modern agriculture. Nested association mapping (NAM) populations are valuable genetic resources that enable incorporation of genetic diversity, dissection of complex traits and providing germplasm to breeding programs. We developed the OzNAM by crossing and backcrossing 73 diverse exotic parents to two Australian elite varieties Gladius and Scout. The NAM parents were genotyped using the iSelect wheat 90K Infinium SNP array, and the progeny were genotyped using a custom targeted genotyping-by-sequencing assay based on molecular inversion probes designed to target 12,179 SNPs chosen from the iSelect wheat 90K Infinium SNP array of the parents. In total, 3535 BC1F4:6 RILs from 125 families with 21 to 76 lines per family were genotyped and we found 4964 polymorphic and multi-allelic haplotype markers that spanned the whole genome. A subset of 530 lines from 28 families were evaluated in multi-environment trials over three years. To demonstrate the utility of the population in QTL mapping, we chose to map QTL for maturity and plant height using the RTM-GWAS approach and we identified novel and known QTL for maturity and plant height.

A Major Locus on Wheat Chromosome 7B Associated With Late-Maturity α-Amylase Encodes a Putative ent-Copalyl Diphosphate Synthase.

Many wheat varieties have the potential to develop unacceptably high levels of α-amylase in the grains if exposed to a cool temperature shock or simply cool temperature during the early to middle stages of grain filling. This phenomenon is referred to as late maturity α-amylase (LMA). The enzyme persists in the grain until harvest and may result in wheat with a low Falling Number that does not meet receival and export specifications. Resistance to LMA is therefore a valuable target for wheat breeders and wheat industries in general. Genetic evidence implicating a locus on the long arm of chromosome 7B in variation in LMA phenotype was confirmed in this investigation. Through intensive fine-mapping an ent-copalyl diphosphate synthase (CPS), hitherto named LMA-1, was identified as the likely candidate gene associated with variation in LMA phenotype. Single Nucleotide Polymorphisms (SNPs) within the LMA-1 coding sequence of Chinese Spring, Maringa and Halberd result in either prematurely terminated or functionally altered proteins that are associated with useful levels of resistance to LMA. LMA-1 transcripts detected in de-embryonated grain tissue from around 15 days after anthesis, several days before the synthesis of α-amylase, were low in the resistant varieties Chinese Spring and Maringa compared with LMA susceptible genotype Spica. This was associated with a dramatic reduction in the concentrations of intermediates in the gibberellin biosynthesis pathway such as GA19, evidence that LMA-1 was functioning as CPS in the gibberellin biosynthesis pathway. A survey of a large collection of Australian and international wheat varieties distinguished 9 major haplotypes at the LMA-1 locus. Generally, within classes, there was notable variation for LMA phenotype and evidence for genotypes whose resistance is presumed to be due to genetic loci located elsewhere on the wheat genome. Further investigation is required to characterize the sequence of steps between LMA-1 and α-amylase synthesis as well as to gain a better understanding of the role and potential impact of other genetic loci. Diagnostic markers for sources of resistance and SNP variation reported in this study should assist breeders to deploy resistance associated with LMA-1 variants in breeding programs.

Improving Nitrogen Use Efficiency Through Overexpression of Alanine Aminotransferase in Rice, Wheat, and Barley.

Nitrogen is an essential nutrient for plants, but crop plants are inefficient in the acquisition and utilization of applied nitrogen. This often results in producers over applying nitrogen fertilizers, which can negatively impact the environment. The development of crop plants with more efficient nitrogen usage is, therefore, an important research goal in achieving greater agricultural sustainability. We utilized genetically modified rice lines over-expressing a barley alanine aminotransferase (HvAlaAT) to help characterize pathways which lead to more efficient use of nitrogen. Under the control of a stress-inducible promoter OsAnt1, OsAnt1:HvAlaAT lines have increased above-ground biomass with little change to both nitrate and ammonium uptake rates. Based on metabolic profiles, carbon metabolites, particularly those involved in glycolysis and the tricarboxylic acid (TCA) cycle, were significantly altered in roots of OsAnt1:HvAlaAT lines, suggesting higher metabolic turnover. Moreover, transcriptomic data revealed that genes involved in glycolysis and TCA cycle were upregulated. These observations suggest that higher activity of these two processes could result in higher energy production, driving higher nitrogen assimilation, consequently increasing biomass production. Other potential mechanisms contributing to a nitrogen-use efficient phenotype include involvements of phytohormonal responses and an alteration in secondary metabolism. We also conducted basic growth studies to evaluate the effect of the OsAnt1:HvAlaAT transgene in barley and wheat, which the transgenic crop plants increased seed production under controlled environmental conditions. This study provides comprehensive profiling of genetic and metabolic responses to the over-expression of AlaAT and unravels several components and pathways which contribute to its nitrogen-use efficient phenotype.

The wheat Seven in absentia gene is associated with increases in biomass and yield in hot climates.

Wheat (Triticum aestivum L.) productivity is severely reduced by high temperatures. Breeding of heat-tolerant cultivars can be achieved by identifying genes controlling physiological and agronomical traits when high temperatures occur and using these to select superior genotypes, but no gene underlying genetic variation for heat tolerance has previously been described. We advanced the positional cloning of qYDH.3BL, a quantitative trait locus (QTL) on bread wheat chromosome 3B associated with increased yield in hot and dry climates. The delimited genomic region contained 12 putative genes and a sequence variant in the promoter region of one gene, Seven in absentia, TaSINA. This was associated with the QTL's effects on early vigour, root growth, plant biomass, and yield components in two distinct wheat populations grown under various growth conditions. Near isogenic lines carrying the positive allele at qYDH.3BL underexpressed TaSINA and had increased vigour and water use efficiency early in development, as well as increased biomass, grain number, and grain weight following heat stress. A survey of worldwide distribution indicated that the positive allele became widespread from the 1950s through the CIMMYT wheat breeding programme but, to date, has been selected only in breeding programmes in Mexico and Australia.

Ph2 encodes the mismatch repair protein MSH7-3D that inhibits wheat homoeologous recombination.

Meiotic recombination is a critical process for plant breeding, as it creates novel allele combinations that can be exploited for crop improvement. In wheat, a complex allohexaploid that has a diploid-like behaviour, meiotic recombination between homoeologous or alien chromosomes is suppressed through the action of several loci. Here, we report positional cloning of Pairing homoeologous 2 (Ph2) and functional validation of the wheat DNA mismatch repair protein MSH7-3D as a key inhibitor of homoeologous recombination, thus solving a half-century-old question. Similar to ph2 mutant phenotype, we show that mutating MSH7-3D induces a substantial increase in homoeologous recombination (up to 5.5 fold) in wheat-wild relative hybrids, which is also associated with a reduction in homologous recombination. These data reveal a role for MSH7-3D in meiotic stabilisation of allopolyploidy and provides an opportunity to improve wheat's genetic diversity through alien gene introgression, a major bottleneck facing crop improvement.

Genetic factors associated with favourable pollinator traits in the wheat cultivar Piko.

Hybrid breeding in wheat has the potential to boost yields. An efficient hybrid seed production system requires elite pollinators; however, such germplasm is limited among modern cultivars. Piko, a winter wheat (Triticum aestivum L.) cultivar, has been identified as a superior pollinator and has been used in Europe. Piko has favourable pollinator traits for anther extrusion, anther length, pollen mass and hybrid seed set. However, the genetic factors responsible for Piko's favourable traits are largely unknown. Here, we report on the genetic analysis of a Piko-derived F2 mapping population. We confirmed that Piko's Rht-D1a allele for tall stature is associated with large anthers and high anther extrusion. However, Rht-D1 was not found to be associated with anther filament length, confirmed by near isogenic lines. Piko's photoperiod sensitive Ppd-B1b allele shows an association with increased spike length, more spikelets and spike architecture traits, while the insensitive Ppd-B1a allele is linked with high anther extrusion and larger anthers. We identified an anther extrusion quantitative trait locus (QTL) on chromosome 6A that showed significantly biased transmission of the favourable Piko allele amongst F2 progenies. The Piko allele is completely absent in the distal 6AS region and the central 6A region revealed a significantly lower ratio (<8%) of F2 with homozygous Piko alleles. Our study provided further evidence for the effects of Rht-D1 and Ppd-B1 loci on multiple pollinator traits and a novel anther extrusion QTL that exhibits segregation distortion.

Hybrid breeding in wheat: how shaping floral biology can offer new perspectives.

Hybrid breeding in wheat (Triticum aestivum L.) has the potential to deliver major yield increases. This is a requisite to guarantee food security for increasing population demands and to counterbalance the effects of extreme environmental conditions. Successful hybrid breeding in wheat relies on forced outcrossing while preventing self-pollination. To achieve this, research has been directed towards identifying and improving fertility control systems. To maximise cross-pollination and seed set, however, fertility control systems need to be complemented by breeding phenotypically distinct male and female lines. This review summarises existing and novel male sterility systems for wheat hybridisation. We also consider the genetic resources that can be used to alter wheat's floral development and spike morphology, with a focus on the genetic variation already available. Exploiting these resources can lead to enhanced outcrossing, a key requirement in the progress towards hybrid wheat breeding.

Dual-locus DNA metabarcoding reveals southern hairy-nosed wombats (Lasiorhinus latifrons Owen) have a summer diet dominated by toxic invasive plants.

Habitat degradation and summer droughts severely restrict feeding options for the endangered southern hairy-nosed wombat (SHNW; Lasiorhinus latifrons). We reconstructed SHNW summer diets by DNA metabarcoding from feces. We initially validated rbcL and ndhJ diet reconstructions using autopsied and captive animals. Subsequent diet reconstructions of wild wombats broadly reflected vegetative ground cover, implying local rather than long-range foraging. Diets were all dominated by alien invasives. Chemical analysis of alien food revealed Carrichtera annua contains high levels of glucosinolates. Clinical examination (7 animals) and autopsy (12 animals) revealed that the most degraded site also contained most individuals showing signs of glucosinolate poisoning. We infer that dietary poisoning through the ingestion of alien invasives may have contributed to the recent population crashes in the region. In floristically diverse sites, individuals appear to be able to manage glucosinolate intake by avoidance or episodic feeding but this strategy is less tractable in the most degraded sites. We conclude that recovery of the most affected populations may require effective Carrichtera management and interim supplementary feeding. More generally, we argue that protection against population decline by poisoning in territorial herbivores requires knowledge of their diet and of those food plants containing toxic principles.