Search

Background: Ross River virus (RRV) is Australia's most common and widespread mosquito-transmitted arbovirus and is of significant public health concern. With increasing anthropogenic impacts on wildlife and mosquito populations, it is important that we understand how RRV circulates in its endemic hotspots to determine where public health efforts should be directed. Current surveillance methods are effective in locating the virus but do not provide data on the circulation of the virus and its strains within the environment. This study examined the ability to identify single nucleotide polymorphisms (SNPs) within the variable E2/E3 region by generating full-length haplotypes from a range of mosquito trap-derived samples., Methods: A novel tiled primer amplification workflow for amplifying RRV was developed with analysis using Oxford Nanopore Technology's MinION and a custom ARTIC/InterARTIC bioinformatic protocol. By creating a range of amplicons across the whole genome, fine-scale SNP analysis was enabled by specifically targeting the variable region that was amplified as a single fragment and established haplotypes that informed spatial-temporal variation of RRV in the study site in Victoria., Results: A bioinformatic and laboratory pipeline was successfully designed and implemented on mosquito whole trap homogenates. Resulting data showed that genotyping could be conducted in real time and that whole trap consensus of the viruses (with major SNPs) could be determined in a timely manner. Minor variants were successfully detected from the variable E2/E3 region of RRV, which allowed haplotype determination within complex mosquito homogenate samples., Conclusions: The novel bioinformatic and wet laboratory methods developed here will enable fast detection and characterisation of RRV isolates. The concepts presented in this body of work are transferable to other viruses that exist as quasispecies in samples. The ability to detect minor SNPs, and thus haplotype strains, is critically important for understanding the epidemiology of viruses their natural environment., (© 2023. The Author(s).)

Key Message: Breeding target traits can be broadened to include nutritive value and plant breeder's rights traits in perennial ryegrass by using in-field regression-based spectroscopy phenotyping and genomic selection. Perennial ryegrass breeding has focused on biomass yield, but expansion into a broader set of traits is needed to benefit livestock industries whilst also providing support for intellectual property protection of cultivars. Numerous breeding objectives can be targeted simultaneously with the development of sensor-based phenomics and genomic selection (GS). Of particular interest are nutritive value (NV), which has been difficult and expensive to measure using traditional phenotyping methods, resulting in limited genetic improvement to date, and traits required to obtain varietal protection, known as plant breeder's rights (PBR) traits. In order to assess phenotyping requirements for NV improvement and potential for genetic improvement, in-field reflectance-based spectroscopy was assessed and GS evaluated in a single population for three key NV traits, captured across four timepoints. Using three prediction approaches, the possibility of targeting PBR traits using GS was evaluated for five traits recorded across three years of a breeding program. Prediction accuracy was generally low to moderate for NV traits and moderate to high for PBR traits, with heritability highly correlated with GS accuracy. NV did not show significant or consistent correlation between timepoints highlighting the need to incorporate seasonal NV into selection indexes and the value of being able to regularly monitor NV across seasons. This study has demonstrated the ability to implement GS for both NV and PBR traits in perennial ryegrass, facilitating the expansion of ryegrass breeding targets to agronomically relevant traits while ensuring necessary varietal protection is achieved., (© 2023. Crown.)

Cannabis is commercially cultivated for both therapeutic and recreational purposes in a growing number of jurisdictions. The main cannabinoids of interest are cannabidiol (CBD) and delta-9 tetrahydrocannabidiol (THC), which have applications in different therapeutic treatments. The rapid, nondestructive determination of cannabinoid levels has been achieved using near-infrared (NIR) spectroscopy coupled to high-quality compound reference data provided by liquid chromatography. However, most of the literature describes prediction models for the decarboxylated cannabinoids, e.g., THC and CBD, rather than naturally occurring analogues, tetrahydrocannabidiolic acid (THCA) and cannabidiolic acid (CBDA). The accurate prediction of these acidic cannabinoids has important implications for quality control for cultivators, manufacturers and regulatory bodies. Using high-quality liquid chromatography-mass spectroscopy (LCMS) data and NIR spectra data, we developed statistical models including principal component analysis (PCA) for data quality control, partial least squares regression (PLS-R) models to predict cannabinoid concentrations for 14 different cannabinoids and partial least squares discriminant analysis (PLS-DA) models to characterise cannabis samples into high-CBDA, high-THCA and even-ratio classes. This analysis employed two spectrometers, a scientific grade benchtop instrument (Bruker MPA II-Multi-Purpose FT-NIR Analyzer) and a handheld instrument (VIAVI MicroNIR Onsite-W). While the models from the benchtop instrument were generally more robust (99.4-100% accuracy prediction), the handheld device also performed well (83.1-100% accuracy prediction) with the added benefits of portability and speed. In addition, two cannabis inflorescence preparation methods were evaluated: finely ground and coarsely ground. The models generated from coarsely ground cannabis provided comparable predictions to that of the finely ground but represent significant timesaving in terms of sample preparation. This study demonstrates that a portable NIR handheld device paired with LCMS quantitative data can provide accurate cannabinoid predictions and potentially be of use for the rapid, high-throughput, nondestructive screening of cannabis material.

Maintaining specific and reproducible cannabinoid compositions (type and quantity) is essential for the production of cannabis-based remedies that are therapeutically effective. The current study investigates factors that determine the plant's cannabinoid profile and examines interrelationships between plant features (growth rate, phenology and biomass), inflorescence morphology (size, shape and distribution) and cannabinoid content. An examination of differences in cannabinoid profile within genotypes revealed that across the cultivation facility, cannabinoids' qualitative traits (ratios between cannabinoid quantities) remain fairly stable, while quantitative traits (the absolute amount of Δ 9 -tetrahydrocannabinol (THC), cannabidiol (CBD), cannabichromene (CBC), cannabigerol (CBG), Δ 9 -tetrahydrocannabivarin (THCV) and cannabidivarin (CBDV)) can significantly vary. The calculated broad-sense heritability values imply that cannabinoid composition will have a strong response to selection in comparison to the morphological and phenological traits of the plant and its inflorescences. Moreover, it is proposed that selection in favour of a vigorous growth rate, high-stature plants and wide inflorescences is expected to increase overall cannabinoid production. Finally, a range of physiological and phenological features was utilised for generating a successful model for the prediction of cannabinoid production. The holistic approach presented in the current study provides a better understanding of the interaction between the key features of the cannabis plant and facilitates the production of advanced plant-based medicinal substances.

Outbreaks of avian influenza virus (AIV) from wild waterfowl into the poultry industry is of upmost significance and is an ongoing and constant threat to the industry. Accurate surveillance of AIV in wild waterfowl is critical in understanding viral diversity in the natural reservoir. Current surveillance methods for AIV involve collection of samples and transportation to a laboratory for molecular diagnostics. Processing of samples using this approach takes more than three days and may limit testing locations to those with practical access to laboratories. In potential outbreak situations, response times are critical, and delays have implications in terms of the spread of the virus that leads to increased economic cost. This study used nanopore sequencing technology for in-field sequencing and subtype characterisation of AIV strains collected from wild bird faeces and poultry. A custom in-field virus screening and sequencing protocol, including a targeted offline bioinformatic pipeline, was developed to accurately subtype AIV. Due to the lack of optimal diagnostic MinION packages for Australian AIV strains the bioinformatic pipeline was specifically targeted to confidently subtype local strains. The method presented eliminates the transportation of samples, dependence on internet access and delivers critical diagnostic information in a timely manner., (© 2022. The Author(s).)

In recent decades with the reacknowledgment of the medicinal properties of Cannabis sativa L. (cannabis) plants, there is an increased demand for high performing cultivars that can deliver quality products for various applications. However, scientific knowledge that can facilitate the generation of advanced cannabis cultivars is scarce. In order to improve cannabis breeding and optimize cultivation techniques, the current study aimed to examine the morphological attributes of cannabis inflorescences using novel image analysis practices. The investigated plant population comprises 478 plants ascribed to 119 genotypes of high-THC or blended THC-CBD ratio that was cultivated under a controlled environment facility. Following harvest, all plants were manually processed and an image of the trimmed and refined inflorescences extracted from each plant was captured. Image analysis was then performed using in-house custom-made software which extracted 8 morphological features (such as size, shape and perimeter) for each of the 127,000 extracted inflorescences. Our findings suggest that environmental factors play an important role in the determination of inflorescences' morphology. Therefore, further studies that focus on genotype X environment interactions are required in order to generate inflorescences with desired characteristics. An examination of the intra-plant inflorescences weight distribution revealed that processing 75% of the plant's largest inflorescences will gain 90% of its overall yield weight. Therefore, for the optimization of post-harvest tasks, it is suggested to evaluate if the benefits from extracting and processing the plant's smaller inflorescences outweigh its operational costs. To advance selection efficacy for breeding purposes, a prediction equation for forecasting the plant's production biomass through width measurements of specific inflorescences, formed under the current experimental methodology, was generated. Thus, it is anticipated that findings from the current study will contribute to the field of medicinal cannabis by improving targeted breeding programs, advancing crop productivity and enhancing the efficacy of post-harvest procedures., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Naim-Feil, Breen, Pembleton, Spooner, Spangenberg and Cogan.)

Heterosis is defined as increased performance of the F1 hybrid relative to its parents. In the current study, a cohort of populations and parents were created to evaluate and understand heterosis across generations (i.e., F1 to F3) in lentil, a self-pollinated annual diploid (2n = 2× = 14) crop species. Lentil plants were evaluated for heterotic traits in terms of plant height, biomass fresh weight, seed number, yield per plant and 100 grain weight. A total of 47 selected lentil genotypes were cross hybridized to generate 72 F1 hybrids. The F1 hybrids from the top five crosses exhibited between 31%-62% heterosis for seed number with reference to the better parent. The five best performing heterotic crosses were selected with a negative control for evaluation at the subsequent F2 generation and only the tails of the distribution taken forward to be assessed in the F3 generation as a sub selection. Overall, heterosis decreases across the subsequent generations for all traits studied. However, some individual genotypes were identified at the F2 and sub-selected F3 generations with higher levels of heterosis than the best F1 mean value (hybrid mimics). The phenotypic data for the selected F2 and sub selected F3 hybrids were analysed, and the study suggested that 100 grain weight was the biggest driver of yield followed by seed number. A genetic diversity analysis of all the F1 parents failed to correlate genetic distance and divergence among parents with heterotic F1's. Therefore, genetic distance was not a key factor to determine heterosis in lentil. The study highlights the challenges associated with different breeding systems for heterosis (i.e., F1 hybrid-based breeding systems and/or via hybrid mimics) but demonstrates the potential significant gains that could be achieved in lentil productivity., Competing Interests: The authors have declared that no competing interests exist.

Cannabis sativa L. produces unique phytocannabinoids, which are used for their pharmaceutical benefits. To date, there are no reports of in vivo engineering targeting the cannabinoid biosynthesis genes to greater elucidate the role each of these genes play in synthesis of these medically important compounds. Reported here is the first modulation of cannabinoid biosynthesis genes using RNAi via agroinfiltration. Vacuum infiltrated leaf segments of the Cannbio-2 C. sativa strain, transfected with different RNAi constructs corresponding to THCAS, CBDAS , and CBCAS gene sequences, showed significant downregulation of all cannabinoid biosynthesis genes using real-time quantitative PCR. Using RNAi, significant off-targeting occurs resulting in the downregulation of highly homologous transcripts. Significant ( p  < 0.05) downregulation was observed for THCAS (92%), CBDAS (97%), and CBCAS (70%) using pRNAi-GG- CBDAS-UNIVERSAL . Significant ( p  < 0.05) upregulation of CBCAS (76%) and non-significant upregulation of THCAS (13%) were observed when transfected with pRNAi-GG- CBCAS , suggesting the related gene's ability to synthesize multiple cannabinoids. Using this approach, increased understanding of the relationship between cannabinoid biosynthesis genes can be further elucidated. This RNAi approach enables functional genomics screens for further reverse genetic studies as well as the development of designer cannabis strains with over-expression and/or downregulation of targeted cannabinoid biosynthesis genes. Functional genomics screens, such as these, will further provide insights into gene regulation of cannabinoid biosynthesis in Cannabis ., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Matchett-Oates, Spangenberg and Cogan.)

Background: Cannabis has been used worldwide for centuries for industrial, recreational and medicinal use, however, to date no successful attempts at editing genes involved in cannabinoid biosynthesis have been reported. This study proposes and develops an in silico best practices approach for the design and implementation of genome editing technologies in cannabis to target all genes involved in cannabinoid biosynthesis., Results: A large dataset of reference genomes was accessed and mined to determine copy number variation and associated SNP variants for optimum target edit sites for genotype independent editing. Copy number variance and highly polymorphic gene sequences exist in the genome making genome editing using CRISPR, Zinc Fingers and TALENs technically difficult. Evaluation of allele or additional gene copies was determined through nucleotide and amino acid alignments with comparative sequence analysis performed. From determined gene copy number and presence of SNPs, multiple online CRISPR design tools were used to design sgRNA targeting every gene, accompanying allele and homologs throughout all involved pathways to create knockouts for further investigation. Universal sgRNA were designed for highly homologous sequences using MultiTargeter and visualised using Sequencher, creating unique sgRNA avoiding SNP and shared nucleotide locations targeting optimal edit sites., Conclusions: Using this framework, the approach has wider applications to all plant species regardless of ploidy number or highly homologous gene sequences., Significance Statement: Using this framework, a best-practice approach to genome editing is possible in all plant species, including cannabis, delivering a comprehensive in silico evaluation of the cannabinoid pathway diversity from a large set of whole genome sequences. Identification of SNP variants across all genes could improve genome editing potentially leading to novel applications across multiple disciplines, including agriculture and medicine., Competing Interests: The authors have declared that no competing interests exist.

Background: For millennia, drug-type cannabis strains were extensively used for various medicinal, ritual, and inebriant applications. However, cannabis prohibition during the last century led to cultivation and breeding activities being conducted under clandestine conditions, while scientific development of the crop ceased. Recently, the potential of medicinal cannabis has been reacknowledged and the now expanding industry requires optimal and scientifically characterized varieties. However, scientific knowledge that can propel this advancement is sorely lacking. To address this issue, the current study aims to provide a better understanding of key physiological and phenological traits that can facilitate the breeding of advanced cultivars., Results: A diverse population of 121 genotypes of high-THC or balanced THC-CBD ratio was cultivated under a controlled environment facility and 13 plant parameters were measured. No physiological association across genotypes attributed to the same vernacular classification was observed. Floral bud dry weight was found to be positively associated with plant height and stem diameter but not with days to maturation. Furthermore, the heritability of both plant height and days to maturation was relatively high, but for plant height it decreased during the vegetative growth phase. To advance breeding efficacy, a prediction equation for forecasting floral bud dry weight was generated, driven by parameters that can be detected during the vegetative growth phase solely., Conclusions: Our findings suggest that selection for taller and fast-growing genotypes is likely to lead to an increase in floral bud productivity. It was also found that the final plant height and stem diameter are determined by 5 independent factors that can be used to maximize productivity through cultivation adjustments. The proposed prediction equation can facilitate the selection of prolific genotypes without the completion of a full cultivation cycle. Future studies that will associate genome-wide variation with plants morphological traits and cannabinoid profile will enable precise and accelerated breeding through genomic selection approaches.

Soil salinity is a major abiotic stress, limiting lentil productivity worldwide. Understanding the genetic basis of salt tolerance is vital to develop tolerant varieties. A diversity panel consisting of 276 lentil accessions was screened in a previous study through traditional and image-based approaches to quantify growth under salt stress. Genotyping was performed using two contrasting methods, targeted (tGBS) and transcriptome (GBS-t) genotyping-by-sequencing, to evaluate the most appropriate methodology. tGBS revealed the highest number of single-base variants (SNPs) (c. 56,349), and markers were more evenly distributed across the genome compared to GBS-t. A genome-wide association study (GWAS) was conducted using a mixed linear model. Significant marker-trait associations were observed on Chromosome 2 as well as Chromosome 4, and a range of candidate genes was identified from the reference genome, the most plausible being potassium transporters, which are known to be involved in salt tolerance in related species. Detailed mineral composition performed on salt-treated and control plant tissues revealed the salt tolerance mechanism in lentil, in which tolerant accessions do not transport Na + ions around the plant instead localize within the root tissues. The pedigree analysis identified two parental accessions that could have been the key sources of tolerance in this dataset.

Molecular characterization of genetically modified plants can provide crucial information for the development of detection and identification methods, to comply with traceability, and labeling requirements prior to commercialization. Detailed description of the genetic modification was previously a challenging step in the safety assessment, since it required the use of laborious and time-consuming techniques. In this study an accurate, simple, and fast method was developed for molecular characterization of genetically modified (GM) plants, following a user-friendly workflow for researchers with limited bioinformatic capabilities. Three GM events from a diverse array of crop species-perennial ryegrass, white clover, and canola-were used to test the approach that exploits long-read sequencing by the MinION device, from Oxford Nanopore Technologies. The method delivered a higher degree of resolution of the transgenic events within the host genome than has previously been possible with the standard Illumina short-range sequencing strategies. The flanking sequences, copy number, and presence of backbone sequences, and overall transgene insertion structure were determined for each of the plant genomes, with the additional identification of moderate-sized secondary insertions that would have previously been missed. The proposed workflow takes only about 1 week from DNA extraction to analyzed result, and the method will complement the existing approaches for molecular characterization of GM plants, since it makes the process faster, simpler, and more cost-effective., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Giraldo, Shinozuka, Spangenberg, Smith and Cogan.)

Cannabis is a diploid species (2 n = 20), the estimated haploid genome sizes of the female and male plants using flow cytometry are 818 and 843 Mb respectively. Although the genome of Cannabis has been sequenced (from hemp, wild and high-THC strains), all assemblies have significant gaps. In addition, there are inconsistencies in the chromosome numbering which limits their use. A new comprehensive draft genome sequence assembly (∼900 Mb) has been generated from the medicinal cannabis strain Cannbio-2, that produces a balanced ratio of cannabidiol and delta-9-tetrahydrocannabinol using long-read sequencing. The assembly was subsequently analysed for completeness by ordering the contigs into chromosome-scale pseudomolecules using a reference genome assembly approach, annotated and compared to other existing reference genome assemblies. The Cannbio-2 genome sequence assembly was found to be the most complete genome sequence available based on nucleotides assembled and BUSCO evaluation in Cannabis sativa with a comprehensive genome annotation. The new draft genome sequence is an advancement in Cannabis genomics permitting pan-genome analysis, genomic selection as well as genome editing., Competing Interests: The authors declare that they have no competing interests., (© The Author(s) 2020.)

Blackleg disease causes yield losses in canola (Brassica napus L.). To identify resistance genes and genomic regions, genome-wide association studies (GWAS) of 585 diverse winter and spring canola accessions were performed using imputed whole-genome sequence (WGS) and transcriptome genotype-by-sequencing (GBSt). Blackleg disease phenotypes were collected across three years in six trials. GWAS were performed in several ways and their respective power was judged by the number of significant single nucleotide polymorphisms (SNP), the false discovery rate (FDR), and the percentage of SNP that validated in additional field trials in two subsequent years. WGS GWAS with 1,234,708 million SNP detected a larger number of significant SNP, achieved a lower FDR and a higher validation rate than GBSt with 64,072 SNP. A meta-analysis combining survival and average internal infection resulted in lower FDR but also lower validation rates. The meta-analysis GWAS identified 79 genomic regions (674 SNP) conferring potential resistance to L. maculans. While several GWAS signals localised in regions of known Rlm genes, fifty-three new potential resistance regions were detected. Seventeen regions had underlying genes with putative functions related to disease defence or stress response in Arabidopsis thaliana. This study provides insight into the genetic architecture and potential molecular mechanisms underlying canola L. maculans resistance.

Plant seeds have long been promoted as a production platform for novel fatty acids such as the ω3 long-chain (≥ C 20 ) polyunsaturated fatty acids eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) commonly found in fish oil. In this article we describe the creation of a canola ( Brassica napus ) variety producing fish oil-like levels of DHA in the seed. This was achieved by the introduction of a microalgal/yeast transgenic pathway of seven consecutive enzymatic steps which converted the native substrate oleic acid to α-linolenic acid and, subsequently, to EPA, docosapentaenoic acid (DPA) and DHA. This paper describes construct design and evaluation, plant transformation, event selection, field testing in a wide range of environments, and oil profile stability of the transgenic seed. The stable, high-performing event NS-B50027-4 produced fish oil-like levels of DHA (9-11%) in open field trials of T 3 to T 7 generation plants in several locations in Australia and Canada. This study also describes the highest seed DHA levels reported thus far and is one of the first examples of a deregulated genetically modified crop with clear health benefits to the consumer., (Copyright © 2020 Petrie, Zhou, Leonforte, McAllister, Shrestha, Kennedy, Belide, Buzza, Gororo, Gao, Lester, Mansour, Mulder, Liu, Tian, Silva, Cogan, Nichols, Green, de Feyter, Devine and Singh.)

Intensive breeding of cultivated lentil has resulted in a relatively narrow genetic base, which limits the options to increase crop productivity through selection. Assessment of genetic diversity in the wild gene pool of lentil, as well as characterization of useful and novel alleles/genes that can be introgressed into elite germplasm, presents new opportunities and pathways for germplasm enhancement, followed by successful crop improvement. In the current study, a lentil collection consisting of 467 wild and cultivated accessions that originated from 10 diverse geographical regions was assessed, to understand genetic relationships among different lentil species/subspecies. A total of 422,101 high-confidence SNP markers were identified against the reference lentil genome (cv. CDC Redberry). Phylogenetic analysis clustered the germplasm collection into four groups, namely, Lens culinaris / Lens orientalis , Lens lamottei / Lens odemensis , Lens ervoides , and Lens nigricans . A weak correlation was observed between geographical origin and genetic relationship, except for some accessions of L . culinaris and L . ervoides . Genetic distance matrices revealed a comparable level of variation within the gene pools of L . culinaris (Nei's coefficient 0.01468-0.71163), L . ervoides (Nei's coefficient 0.01807-0.71877), and L . nigricans (Nei's coefficient 0.02188-1.2219). In order to understand any genic differences at species/subspecies level, allele frequencies were calculated from a subset of 263 lentil accessions. Among all cultivated and wild lentil species, L . nigricans exhibited the greatest allelic differentiation across the genome compared to all other species/subspecies. Major differences were observed on six genomic regions with the largest being on Chromosome 1 (c. 1 Mbp). These results indicate that L. nigricans is the most distantly related to L. culinaris and additional structural variations are likely to be identified from genome sequencing studies. This would provide further insights into evolutionary relationships between cultivated and wild lentil germplasm, for germplasm improvement and introgression., (Copyright © 2020 Dissanayake, Braich, Cogan, Smith and Kaur.)

Genomic selection accelerates genetic progress in crop breeding through the prediction of future phenotypes of selection candidates based on only their genomic information. Here we report genetic correlations and genomic prediction accuracies in 22 agronomic, disease, and seed quality traits measured across multiple years (2015-2017) in replicated trials under rain-fed and irrigated conditions in Victoria, Australia. Two hundred and two spring canola lines were genotyped for 62,082 Single Nucleotide Polymorphisms (SNPs) using transcriptomic genotype-by-sequencing (GBSt). Traits were evaluated in single trait and bivariate genomic best linear unbiased prediction (GBLUP) models and cross-validation. GBLUP were also expanded to include genotype-by-environment G × E interactions. Genomic heritability varied from 0.31to 0.66. Genetic correlations were highly positive within traits across locations and years. Oil content was positively correlated with most agronomic traits. Strong, not previously documented, negative correlations were observed between average internal infection (a measure of blackleg disease) and arachidic and stearic acids. The genetic correlations between fatty acid traits followed the expected patterns based on oil biosynthesis pathways. Genomic prediction accuracy ranged from 0.29 for emergence count to 0.69 for seed yield. The incorporation of G × E translates into improved prediction accuracy by up to 6%. The genomic prediction accuracies achieved indicate that genomic selection is ready for application in canola breeding.

Potatoes are an important human food crop, but have a number of yield limiting factors, including disease susceptibility. Potato virus Y (PVY) is found worldwide, and is one of the main virus problems for potato growers. PVY is transmitted by aphids and mechanically by machinery, tools and people, and symptoms are variable across cultivars and strains, including being symptomless in some cultivars. Therefore, breeding resistant cultivars is the best way to control this virus. This study phenotypically screened 74 of the main commercial cultivars and a few other select cultivars grown in Australia, in order to identify sources of resistance to PVY. The cultivars were screened against PVY O and PVY NTN , with 23 out of 71 resistant to PVY O and 13 out of 74 resistant to PVY NTN , and all these 13 were resistant to both strains. When the phenotypic screening was compared to the results listed on the European Cultivated Potato Database, the majority of results were found to be consistent. We then evaluated three molecular markers RYSC3, M45, and STM0003 for the extreme resistance genes Ry adg and Ry sto , to validate the usefulness of the markers for marker-assisted selection (MAS) on Australian germplasm. The degree of correlation between the resistance phenotypes and the RYSC3, M45, and STM0003 markers for Ry adg and Ry sto conferred PVY resistance was determined. Three cultivars amplified the RYSC3 marker, while the M45 marker amplified the same 3 and an additional 9. Of the 12 cultivars, 11 phenotyped as resistant, but 1 was susceptible. The STM0003 marker was amplified from only 2 cultivars that both had resistant phenotypes. The RYSC3, M45, and STM0003 markers were therefore able to identify all the 13 cultivars that were resistant to both strains of PVY. Therefore, these markers will enable the identification of genotypes with resistance to PVY, and enable PVY resistant parents to be used for the development of superior progeny; these genetic markers can be used for MAS in the Australian potato breeding program.

Food security is one of major concerns for the growing global population. Modern agricultural biotechnologies, such as genetic modification, are a possible solution through enabling an increase of production, more efficient use of natural resources, and reduced environmental impacts. However, new crop varieties with altered genetic materials may be subjected to safety assessments to fulfil the regulatory requirements, prior to marketing. The aim of the assessment is to evaluate the impact of products from the new crop variety on human, animal, and the environmental health. Although, many studies on the risk assessment of genetically modified (GM) food have been published, little consideration to GM feedstuff has been given, despite that between 70 to 90% of all GM crops and their biomass are used as animal feed. In addition, in some GM plants such as forages that are only used for animal feeds, the assessment of the genetic modification may be of relevance only to livestock feeding. In this article, the regulatory framework of GM crops intended for animal feed is reviewed using the available information on GM food as the baseline. Although, the majority of techniques used for the safety assessment of GM food can be used in GM feed, many plant parts used for livestock feeding are inedible to humans. Therefore, the concentration of novel proteins in different plant tissues and level of exposure to GM feedstuff in the diet of target animals should be considered. A further development of specific methodologies for the assessment of GM crops intended for animal consumption is required, in order to provide a more accurate and standardized assessment to the GM feed safety., (Copyright © 2019 Giraldo, Shinozuka, Spangenberg, Cogan and Smith.)

Breeding schemes that utilize modern breeding methods like genomic selection (GS) and speed breeding (SB) have the potential to accelerate genetic gain for different crops. We investigated through stochastic computer simulation the advantages and disadvantages of adopting both GS and SB (SpeedGS) into commercial breeding programs for allogamous crops. In addition, we studied the effect of omitting one or two selection stages from the conventional phenotypic scheme on GS accuracy, genetic gain, and inbreeding. As an example, we simulated GS and SB for five traits (heading date, forage yield, seed yield, persistency, and quality) with different genetic architectures and heritabilities (0.7, 0.3, 0.4, 0.1, and 0.3; respectively) for a tall fescue breeding program. We developed a new method to simulate correlated traits with complex architectures of which effects can be sampled from multiple distributions, e.g. to simulate the presence of both minor and major genes. The phenotypic selection scheme required 11 years, while the proposed SpeedGS schemes required four to nine years per cycle. Generally, SpeedGS schemes resulted in higher genetic gain per year for all traits especially for traits with low heritability such as persistency. Our results showed that running more SB rounds resulted in higher genetic gain per cycle when compared to phenotypic or GS only schemes and this increase was more pronounced per year when cycle time was shortened by omitting cycle stages. While GS accuracy declined with additional SB rounds, the decline was less in round three than in round two, and it stabilized after the fourth SB round. However, more SB rounds resulted in higher inbreeding rate, which could limit long-term genetic gain. The inbreeding rate was reduced by approximately 30% when generating the initial population for each cycle through random crosses instead of generating half-sib families. Our study demonstrated a large potential for additional genetic gain from combining GS and SB. Nevertheless, methods to mitigate inbreeding should be considered for optimal utilization of these highly accelerated breeding programs., (Copyright © 2019 Jighly, Lin, Pembleton, Cogan, Spangenberg, Hayes and Daetwyler.)

Cannabinoids are the main medicinal compounds of interest in the plant Cannabis sativa, that are primarily synthesised in the glandular trichomes; found on female floral buds. The content, composition and yield of secondary metabolites (cannabinoids and terpenoids) is influenced by the plant's genetics and environment. Some initial gene expression experiments have been performed from strains of this plant species that contrasted in cannabinoid production, however the present knowledge about detailed trichome transcriptomics in this species is limited. An extensive transcriptome atlas was generated by RNA sequencing using root, shoot, flower and trichome tissues from a female plant strain (Cannbio-2) and was enhanced with the addition of vegetative and reproductive tissues from a male cannabis plant. Differential gene expression analysis identified genes preferentially expressed in different tissues. Detailed trichomics was performed from extractions specifically from glandular trichomes as well as female floral tissues at varying developmental stages, to identify stage-specific differentially expressed genes. Candidate genes involved in terpene and cannabinoid synthesis were identified and the majority were found to have an abundant expression in trichomes. The comprehensive transcriptome is a significant resource in cannabis for further research of functional genomics to improve the yield of specialised metabolites with high pharmacological value.

Trap-based surveillance strategies are widely used for monitoring of invasive insect species, aiming to detect newly arrived exotic taxa as well as track the population levels of established or endemic pests. Where these surveillance traps have low specificity and capture non-target endemic species in excess of the target pests, the need for extensive specimen sorting and identification creates a major diagnostic bottleneck. While the recent development of standardized molecular diagnostics has partly alleviated this requirement, the single specimen per reaction nature of these methods does not readily scale to the sheer number of insects trapped in surveillance programmes. Consequently, target lists are often restricted to a few high-priority pests, allowing unanticipated species to avoid detection and potentially establish populations. DNA metabarcoding has recently emerged as a method for conducting simultaneous, multi-species identification of complex mixed communities and may lend itself ideally to rapid diagnostics of bulk insect trap samples. Moreover, the high-throughput nature of recent sequencing platforms could enable the multiplexing of hundreds of diverse trap samples on a single flow cell, thereby providing the means to dramatically scale up insect surveillance in terms of both the quantity of traps that can be processed concurrently and number of pest species that can be targeted. In this review of the metabarcoding literature, we explore how DNA metabarcoding could be tailored to the detection of invasive insects in a surveillance context and highlight the unique technical and regulatory challenges that must be considered when implementing high-throughput sequencing technologies into sensitive diagnostic applications., (© The Author(s) 2019. Published by Oxford University Press.)

Despite the high accuracy of short read sequencing (SRS), there are still issues with attaining accurate single nucleotide polymorphism (SNP) genotypes at low sequencing coverage and in highly duplicated genomes due to misalignment. Long read sequencing (LRS) systems, including the Oxford Nanopore Technologies (ONT) minION, have become popular options for de novo genome assembly and structural variant characterisation. The current high error rate often requires substantial post-sequencing correction and would appear to prevent the adoption of this system for SNP genotyping, but nanopore sequencing errors are largely random. Using low coverage ONT minION sequencing for genotyping of pre-validated SNP loci was examined in 9 canola doubled haploids. The minION genotypes were compared to the Illumina sequences to determine the extent and nature of genotype discrepancies between the two systems. The significant increase in read length improved alignment to the genome and the absence of classical SRS biases results in a more even representation of the genome. Sequencing errors are present, primarily in the form of heterozygous genotypes, which can be removed in completely homozygous backgrounds but requires more advanced bioinformatics in heterozygous genomes. Developments in this technology are promising for routine genotyping in the future.

Potato is an important food crop due to its increasing consumption, and as a result, there is demand for varieties with improved production. However, the current status of breeding for improved varieties is a long process which relies heavily on phenotypic evaluation and dated molecular techniques and has little emphasis on modern genotyping approaches. Evaluation and selection before a cultivar is commercialized typically takes 10-15 years. Molecular markers have been developed for disease and pest resistance, resulting in initial marker-assisted selection in breeding. This study has evaluated and implemented a high-throughput transcriptome sequencing method for dense marker discovery in potato for the application of genomic selection. An Australian relevant collection of commercial cultivars was selected, and identification and distribution of high quality SNPs were examined using standard bioinformatic pipelines and a custom approach for the prediction of allelic dosage. As a result, a large number of SNP markers were identified and filtered to generate a high-quality subset that was then combined with historic phenotypic data to assess the approach for genomic selection. Genomic selection potential was predicted for highly heritable traits and the approach demonstrated advantages over the previously used technologies in terms of markers identified as well as costs incurred. The high-quality SNP list also provided acceptable genome coverage which demonstrates its applicability for much larger future studies. This SNP list was also annotated to provide an indication of function and will serve as a resource for the community in future studies. Genome wide marker tools will provide significant benefits for potato breeding efforts and the application of genomic selection will greatly enhance genetic progress.

Implementation of molecular biotechnology, such as transgenic technologies, in forage species can improve agricultural profitability through achievement of higher productivity, better use of resources such as soil nutrients, water, or light, and reduced environmental impact. Development of detection and quantification techniques for genetically modified plants are necessary to comply with traceability and labeling requirements prior to regulatory approval for release. Real-time PCR has been the standard method used for detection and quantification of genetically modified events, and droplet digital PCR is a recent alternative technology that offers a higher accuracy. Evaluation of both technologies was performed using a transgenic high-energy forage grass as a case study. Two methods for detection and quantification of the transgenic cassette, containing modified fructan biosynthesis genes, and a selectable marker gene, hygromycin B phosphotransferase used for transformation, were developed. Real-time PCR was assessed using two detection techniques, SYBR Green I and fluorescent probe-based methods. A range of different agricultural commodities were tested including fresh leaves, tillers, seeds, pollen, silage and hay, simulating a broad range of processed agricultural commodities that are relevant in the commercial use of genetically modified pastures. The real-time and droplet digital PCR methods were able to detect both exogenous constructs in all agricultural products. However, a higher sensitivity and repeatability in transgene detection was observed with the droplet digital PCR technology. Taking these results more broadly, it can be concluded that the droplet digital PCR technology provides the necessary resolution for quantitative analysis and detection, allowing absolute quantification of the target sequence at the required limits of detection across all jurisdictions globally. The information presented here provides guidance and resources for pasture-based biotechnology applications that are required to comply with traceability requirements.

Whole genome sequencing offers genome wide, unbiased markers, and inexpensive library preparation. With the cost of sequencing decreasing rapidly, many plant genomes of modest size are amenable to skim whole genome resequencing (skim WGR). The use of skim WGR in diverse sample sets without the use of imputation was evaluated in silico in 149 canola samples representative of global diversity. Fastq files with an average of 10x coverage of the reference genome were used to generate skim samples representing 0.25x, 0.5x, 1x, 2x, 3x, 4x, and 5x sequencing coverage. Applying a pre-defined list of SNPs versus de novo SNP discovery was evaluated. As skim WGR is expected to result in some degree of insufficient allele sampling, all skim coverage levels were filtered at a range of minimum read depths from a relaxed minimum read depth of 2 to a stringent read depth of 5, resulting in 28 list-based SNP sets. As a broad recommendation, genotyping pre-defined SNPs between 1x and 2x coverage with relatively stringent depth filtering is appropriate for a diverse sample set of canola due to a balance between marker number, sufficient accuracy, and sequencing cost, but depends on the intended application. This was experimentally examined in two sample sets with different genetic backgrounds: 1x coverage of 1,590 individuals from 84 Australian spring type four-parent crosses aimed at maximizing diversity as well as one commercial F 1 hybrid, and 2x coverage of 379 doubled haploids (DHs) derived from a subset of the four-parent crosses. To determine optimal coverage in a simpler genetic background, the DH sample sequence coverage was further down sampled in silico . The flexible and cost-effective nature of the protocol makes it highly applicable across a range of species and purposes.

Alkaloid concentration of perennial ryegrass herbage is affected by endophyte strain and host plant genotype. However, previous studies suggest that associations between host and endophyte also depends on environmental conditions, especially those affecting nutrient reserves and that water-soluble carbohydrate (WSC) concentration of perennial ryegrass plants may influence grass-endophyte associations. In this study a single transgenic event, with altered expression of fructosyltransferase genes to produce high WSC and biomass, has been crossed into a range of cultivar backgrounds with varying Epichloë endophyte strains. The effect of the association between the transgenic trait and alkaloid production was assessed and compared with transgene free control populations. In the vast-majority of comparisons there was no significant difference between alkaloid concentrations of transgenic and non-transgenic plants within the same cultivar and endophyte backgrounds. There was no significant difference between GOI+ (gene of interest positive) and GOI- (gene of interest negative) populations in Janthritrem response. Peramine concentration was not different between GOI+ and GOI- for 10 of the 12 endophytes-cultivar combinations. Cultivar Trojan infected with NEA6 and Alto with SE (standard endophyte) exhibited higher peramine and lolitrem B (only for Alto SE) concentration, in the control GOI- compared with GOI+. Similarly, cultivar Trojan infected with NEA6 and Alto with NEA3 presented higher ergovaline concentration in GOI-. Differences in alkaloid concentration may be attributable to an indirect effect in the modulation of fungal biomass. These results conclude that the presence of this transgenic insertion, does not alter the risk (toxicity) of the endophyte-grass associations. Endophyte-host interactions are complex and further research into associations with high WSC plant should be performed in a case by case basis.

Key Message: Exploitation of data from a ryegrass breeding program has enabled rapid development and implementation of genomic selection for sward-based biomass yield with a twofold-to-threefold increase in genetic gain. Genomic selection, which uses genome-wide sequence polymorphism data and quantitative genetics techniques to predict plant performance, has large potential for the improvement in pasture plants. Major factors influencing the accuracy of genomic selection include the size of reference populations, trait heritability values and the genetic diversity of breeding populations. Global diversity of the important forage species perennial ryegrass is high and so would require a large reference population in order to achieve moderate accuracies of genomic selection. However, diversity of germplasm within a breeding program is likely to be lower. In addition, de novo construction and characterisation of reference populations are a logistically complex process. Consequently, historical phenotypic records for seasonal biomass yield and heading date over a 18-year period within a commercial perennial ryegrass breeding program have been accessed, and target populations have been characterised with a high-density transcriptome-based genotyping-by-sequencing assay. Ability to predict observed phenotypic performance in each successive year was assessed by using all synthetic populations from previous years as a reference population. Moderate and high accuracies were achieved for the two traits, respectively, consistent with broad-sense heritability values. The present study represents the first demonstration and validation of genomic selection for seasonal biomass yield within a diverse commercial breeding program across multiple years. These results, supported by previous simulation studies, demonstrate the ability to predict sward-based phenotypic performance early in the process of individual plant selection, so shortening the breeding cycle, increasing the rate of genetic gain and allowing rapid adoption in ryegrass improvement programs.

Genomic prediction is becoming a popular plant breeding method to predict the genetic merit of lines. While some genomic prediction results have been reported in canola, none have been evaluated for blackleg disease. Here, we report genomic prediction for seedling emergence, survival rate, and internal infection), using 532 Spring and Winter canola lines. These lines were phenotyped in two replicated blackleg disease nurseries grown at Wickliffe and Green Lake, Victoria, Australia. A transcriptome genotyping-by-sequencing approach revealed 98,054 single nucleotide polymorphisms (SNPs) after quality control. We assessed various genomic prediction scenarios based on Genomic Best Linear Unbiased Prediction (GBLUP), BayesR and BayesRC, which can make use of prior quantitative trait loci information, via cross-validation. Clustering based on genomic relationships showed that Winter and Spring lines were genetically distinct, indicating limited gene flow between sets. Genetic correlations within traits between Spring and Winter lines ranged from 0.68 and 0.90 (mean = 0.76). Based on GBLUP in the whole population, moderate to high genomic prediction accuracies were achieved within environments (0.35-0.74) and were reduced across environments (0.28-0.58). Prediction accuracy within the Spring set ranged from 0.30-0.69, and from 0.19-0.71 within the Winter set. The BayesR model resulted in slightly lower accuracy to GBLUP. The proportion of genetic variance explained by known blackleg quantitative trait loci (QTL) was < 30%, indicating that there is a large reservoir of genetic variation in blackleg traits that remains to be discovered, but can be captured with genomic prediction. However, providing prior information of known QTL in the BayesRC method resulted in an increased prediction accuracy for survival and internal infection, particularly with Spring lines. Overall, these promising results indicate that genomic prediction will be a valuable tool to make use of all genetic variation to improve blackleg resistance in canola., (Copyright © 2018 Crop Science Society of America.)

Intensive breeding of Brassica napus has resulted in relatively low diversity, such that B. napus would benefit from germplasm improvement schemes that sustain diversity. As such, samples representative of global germplasm pools need to be assessed for existing population structure, diversity and linkage disequilibrium (LD). Complexity reduction genotyping-by-sequencing (GBS) methods, including GBS-transcriptomics (GBS-t), enable cost-effective screening of a large number of samples, while whole genome re-sequencing (WGR) delivers the ability to generate large numbers of unbiased genomic single nucleotide polymorphisms (SNPs), and identify structural variants (SVs). Furthermore, the development of genomic tools based on whole genomes representative of global oilseed diversity and orientated by the reference genome has substantial industry relevance and will be highly beneficial for canola breeding. As recent studies have focused on European and Chinese varieties, a global diversity panel as well as a substantial number of Australian spring types were included in this study. Focusing on industry relevance, 633 varieties were initially genotyped using GBS-t to examine population structure using 61,037 SNPs. Subsequently, 149 samples representative of global diversity were selected for WGR and both data sets used for a side-by-side evaluation of diversity and LD. The WGR data was further used to develop genomic resources consisting of a list of 4,029,750 high-confidence SNPs annotated using SnpEff, and SVs in the form of 10,976 deletions and 2,556 insertions. These resources form the basis of a reliable and repeatable system allowing greater integration between canola genomics studies, with a strong focus on breeding germplasm and industry applicability.

The application of genomics in crops has the ability to significantly improve genetic gain for agriculture. Many marker-dense tools have been developed, but few have seen broad adoption in plant genomics due to issues of significant variations of genome size, levels of ploidy, single nucleotide polymorphism (SNP) frequency and reproductive habit. When combined with limited breeding activities, small research communities and scant sequence resources, the suitability of popular systems is often suboptimal and routinely fails to effectively balance cost-effectiveness and sample throughput. Genotyping-by-sequencing (GBS) encompasses a range of protocols including resequencing of the transcriptome. This study describes a skim GBS-transcriptomics (GBS-t) approach developed to be broadly applicable, cost-effective and high-throughput while still assaying a significant number of SNP loci. A range of crop species with differing levels of ploidy and degree of inbreeding/outbreeding were chosen, including perennial ryegrass, a diploid outbreeding forage grass; phalaris, a putative segmental allotetraploid outbreeding forage grass; lentil, a diploid inbreeding grain legume; and canola, an allotetraploid partially outbreeding oilseed. GBS-t was validated as a simple and largely automated, cost-effective method which generates sufficient SNPs (from 89 738 to 231 977) with acceptable levels of missing data and even genome coverage from c. 3 million sequence reads per sample. GBS-t is therefore a broadly applicable system suitable for many crops, offering advantages over other systems. The correct choice of subsequent sequence analysis software is important, and the bioinformatics process should be iterative and tailored to the specific challenges posed by ploidy variation and extent of heterozygosity., (© 2017 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

The current Illumina HiSeq and MiSeq platforms can generate paired-end reads of up to 2 x 250 bp and 2 x 300 bp in length, respectively. These read lengths may be substantially longer than genomic regions of interest when a DNA sequencing library is prepared through a target enrichment-based approach. A sequencing library preparation method has been developed based on the homology-based enzymatic DNA fragment assembly scheme to allow processing of multiple PCR products within a single read. Target sequences were amplified using locus-specific PCR primers with 8 bp tags, and using the tags, homology-based enzymatic DNA assembly was performed with DNA polymerase, T7 exonuclease and T4 DNA ligase. Short PCR amplicons can hence be assembled into a single molecule, along with sequencing adapters specific to the Illumina platforms. As a proof-of-concept experiment, short PCR amplicons (57-66 bp in length) derived from genomic DNA templates of field pea and containing variable nucleotide locations were assembled and sequenced on the MiSeq platform. The results were validated with other genotyping methods. When 5 PCR amplicons were assembled, 4.3 targeted sequences (single-nucleotide polymorphisms) on average were successfully identified within each read. The utility of this for sequencing of short fragments has consequently been demonstrated., (© The Author(s) 2018. Published by Oxford University Press.)

Effective vector and arbovirus surveillance requires timely and accurate screening techniques that can be easily upscaled. Next-generation sequencing (NGS) is a high-throughput technology that has the potential to modernize vector surveillance. When combined with DNA barcoding, it is termed 'metabarcoding.' The aim of our study was to establish a metabarcoding protocol to characterize pools of mosquitoes and screen them for virus. Pools contained 100 morphologically identified individuals, including one Ross River virus (RRV) infected mosquito, with three species present at different proportions: 1, 5, 94%. Nucleic acid extracted from both crude homogenate and supernatant was used to amplify a 269-bp section of the mitochondrial cytochrome c oxidase subunit I (COI) locus. Additionally, a 67-bp region of the RRV E2 gene was amplified from synthesized cDNA to screen for RRV. Amplicon sequencing was performed using an Illumina MiSeq, and bioinformatic analysis was performed using a DNA barcode database of Victorian mosquitoes. Metabarcoding successfully detected all mosquito species and RRV in every positive sample tested. The limits of species detection were also examined by screening a pool of 1000 individuals, successfully identifying the species and RRV from a single mosquito. The primers used for amplification, number of PCR cycles and total number of individuals present all have effects on the quantification of species in mixed bulk samples. Based on the results, a number of recommendations for future metabarcoding studies are presented. Overall, metabarcoding shows great promise for providing a new alternative approach to screening large insect surveillance trap catches., (© 2017 The Authors. Molecular Ecology Resources Published by John Wiley & Sons Ltd.)

RNA-Seq methodology has been used to generate a comprehensive transcriptome sequence resource for perennial ryegrass, an important temperate pasture grass species. A total of 931 547 255 reads were obtained from libraries corresponding to 19 distinct tissue samples, including both vegetative and reproductive stages of development. Assembly of data generated a final filtered reference set of 48 713 contigs and scaffolds. The transcriptome resource will support whole genome sequence assembly, comparative genomics, implementation of genotyping-by-sequencing (GBS) methods based on transcript sampling, and identification of candidate genes for multiple biological functions.

Molecular characterisation has convincingly demonstrated some types of horizontal gene transfer in eukaryotes, but nuclear gene transfer between distantly related eukaryotic groups appears to have been rare. For angiosperms (flowering plants), nuclear gene transfer events identified to date have been confined to genes originating from prokaryotes or other plant species. In this report, evidence for ancient horizontal transfer of a fungal nuclear gene, encoding a ß-1,6-glucanase enzyme for fungal cell wall degradation, into an angiosperm lineage is presented for the first time. The gene was identified from de novo sequencing and assembly of the genome and transcriptome of perennial ryegrass, a cool-season grass species. Molecular analysis confirmed the presence of the complete gene in the genome of perennial ryegrass. No corresponding sequence was found in other plant species, apart from members of the Poeae sub-tribes Loliinae and Dactylidinae. Evidence suggests that a common ancestor of the two sub-tribes acquired the gene from a species ancestral to contemporary grass-associated fungal endophytes around 9-13 million years ago. This first report of horizontal transfer of a nuclear gene from a taxonomically distant eukaryote to modern flowering plants provides evidence for a novel adaptation mechanism in angiosperms.