CLAW: An automated Snakemake workflow for the assembly of chloroplast genomes from long-read data.

Chloroplasts are photosynthetic organelles in algal and plant cells that contain their own genome. Chloroplast genomes are commonly used in evolutionary studies and taxonomic identification and are increasingly becoming a target for crop improvement studies. As DNA sequencing becomes more affordable, researchers are collecting vast swathes of high-quality whole-genome sequence data from laboratory and field settings alike. Whole tissue read libraries sequenced with the primary goal of understanding the nuclear genome will inadvertently contain many reads derived from the chloroplast genome. These whole-genome, whole-tissue read libraries can additionally be used to assemble chloroplast genomes with little to no extra cost. While several tools exist that make use of short-read second generation and third-generation long-read sequencing data for chloroplast genome assembly, these tools may have complex installation steps, inadequate error reporting, poor expandability, and/or lack scalability. Here, we present CLAW (Chloroplast Long-read Assembly Workflow), an easy to install, customise, and use Snakemake tool to assemble chloroplast genomes from chloroplast long-reads found in whole-genome read libraries (https://github.com/aaronphillips7493/CLAW). Using 19 publicly available reference chloroplast genome assemblies and long-read libraries from algal, monocot and eudicot species, we show that CLAW can rapidly produce chloroplast genome assemblies with high similarity to the reference assemblies. CLAW was designed such that users have complete control over parameterisation, allowing individuals to optimise CLAW to their specific use cases. We expect that CLAW will provide researchers (with varying levels of bioinformatics expertise) with an additional resource useful for contributing to the growing number of publicly available chloroplast genome assemblies.

Fructilactobacillus cliffordii sp. nov., Fructilactobacillus hinvesii sp. nov., Fructilactobacillus myrtifloralis sp. nov., Fructilactobacillus carniphilus sp. nov. and Fructobacillus americanaquae sp. nov., five novel lactic acid bacteria isolated from insects or flowers of Kangaroo Island, South Australia.

Six strains, KI11_D11T, KI4_B1, KI11_C11T, KI16_H9T, KI4_A6T and KI3_B9T, were isolated from insects and flowers on Kangaroo Island, South Australia. On the basis of 16S rRNA gene phylogeny, strains KI11_D11T, KI4_B1, KI11_C11T, KI16_H9T, KI4_A6T were found to be closely related to Fructilactobacillus ixorae Ru20-1T. Due to the lack of a whole genome sequence for this species, whole genome sequencing of Fructilactobacillus ixorae Ru20-1T was undertaken. KI3_B9T was found to be closely related to Fructobacillus tropaeoli F214-1T. Utilizing core gene phylogenetics and whole genome analyses, such as determination of AAI, ANI and dDDH, we propose that these six isolates represent five novel species with the names Fructilactobacillus cliffordii (KI11_D11T= LMG 32130T = NBRC 114988T), Fructilactobacillus hinvesii (KI11_C11T = LMG 32129T = NBRC 114987T), Fructilactobacillus myrtifloralis (KI16_H9T= LMG 32131T = NBRC 114989T) Fructilactobacillus carniphilus (KI4_A6T = LMG 32127T = NBRC 114985T) and Fructobacillus americanaquae (KI3_B9T = LMG 32124T = NBRC 114983T). Chemotaxonomic analyses detected no fructophilic characters for these strains of member of the genus Fructilactobacillus. KI3_B9T was found to be obligately fructophilic, similarly to its phylogenetic neighbours in the genus Fructobacillus. This study represents the first isolation, to our knowledge, of novel species in the family Lactobacillaceae from the Australian wild.

A chromosome-level genome assembly of Plantago ovata.

Plantago ovata is cultivated for production of its seed husk (psyllium). When wet, the husk transforms into a mucilage with properties suitable for pharmaceutical industries, utilised in supplements for controlling blood cholesterol levels, and food industries for making gluten-free products. There has been limited success in improving husk quantity and quality through breeding approaches, partly due to the lack of a reference genome. Here we constructed the first chromosome-scale reference assembly of P. ovata using a combination of 5.98 million PacBio and 636.5 million Hi-C reads. We also used corrected PacBio reads to estimate genome size and transcripts to generate gene models. The final assembly covers ~ 500 Mb with 99.3% gene set completeness. A total of 97% of the sequences are anchored to four chromosomes with an N50 of ~ 128.87 Mb. The P. ovata genome contains 61.90% repeats, where 40.04% are long terminal repeats. We identified 41,820 protein-coding genes, 411 non-coding RNAs, 108 ribosomal RNAs, and 1295 transfer RNAs. This genome will provide a resource for plant breeding programs to, for example, reduce agronomic constraints such as seed shattering, increase psyllium yield and quality, and overcome crop disease susceptibility.

Apilactobacillus apisilvae sp. nov., Nicolia spurrieriana gen. nov. sp. nov., Bombilactobacillus folatiphilus sp. nov. and Bombilactobacillus thymidiniphilus sp. nov., four new lactic acid bacterial isolates from stingless bees Tetragonula carbonaria and Austroplebeia australis.

Four strains, SG5_A10T, SGEP1_A5T, SG4_D2T, and SG4_A1T, were isolated from the honey or homogenate of Australian stingless bee species Tetragonula carbonaria and Austroplebeia australis. Based on 16S rRNA gene phylogeny, core gene phylogenetics, whole genome analyses such as determination of amino acid identity (AAI), cAAI of conserved genes, average nucleotide identity (ANI), and digital DNA-DNA hybridization (dDDH), chemotaxonomic analyses, and the novel isolation sources and unique geography, we propose three new species and one genus with the names Apilactobacillus apisilvae sp. nov. (SG5_A10T = LMG 32133T = NBRC 114991T), Bombilactobacillus thymidiniphilus sp. nov. (SG4_A1T = LMG 32125T = NBRC 114984T), Bombilactobacillus folatiphilus sp. nov. (SG4_D2T = LMG 32126T = NBRC 115004T) and Nicolia spurrieriana sp. nov. (SGEP1_A5T = LMG 32134T = NBRC 114992T). Three out of the four strains were found to be fructophilic, where SG5_A10T and SGEP1_A5T belong to obligately fructophilic lactic acid bacteria, and SG4_D2T representing a new type denoted here as kinetically fructophilic. This study represents the first published lactic acid bacterial species associated with the unique niche of Australian stingless bees.

The first long-read nuclear genome assembly of Oryza australiensis, a wild rice from northern Australia.

Oryza australiensis is a wild rice native to monsoonal northern Australia. The International Oryza Map Alignment Project emphasises its significance as the sole representative of the EE genome clade. Assembly of the O. australiensis genome has previously been challenging due to its high Long Terminal Repeat (LTR) retrotransposon (RT) content. Oxford Nanopore long reads were combined with Illumina short reads to generate a high-quality ~ 858 Mbp genome assembly within 850 contigs with 46× long read coverage. Reference-guided scaffolding increased genome contiguity, placing 88.2% of contigs into 12 pseudomolecules. After alignment to the Oryza sativa cv. Nipponbare genome, we observed several structural variations. PacBio Iso-Seq data were generated for five distinct tissues to improve the functional annotation of 34,587 protein-coding genes and 42,329 transcripts. We also report SNV numbers for three additional O. australiensis genotypes based on Illumina re-sequencing. Although genetic similarity reflected geographical separation, the density of SNVs also correlated with our previous report on variations in salinity tolerance. This genome re-confirms the genetic remoteness of the O. australiensis lineage within the O. officinalis genome complex. Assembly of a high-quality genome for O. australiensis provides an important resource for the discovery of critical genes involved in development and stress tolerance.

Identifying the genetic control of salinity tolerance in the bread wheat landrace Mocho de Espiga Branca.

Salinity tolerance in bread wheat is frequently reported to be associated with low leaf sodium (Na+) concentrations. However, the Portuguese landrace, Mocho de Espiga Branca, accumulates significantly higher leaf Na+ but has comparable salinity tolerance to commercial bread wheat cultivars. To determine the genetic loci associated with the salinity tolerance of this landrace, an F2 mapping population was developed by crossing Mocho de Espiga Branca with the Australian cultivar Gladius. The population was phenotyped for 19 salinity tolerance subtraits using both non-destructive and destructive techniques. Genotyping was performed using genotyping-by-sequencing (GBS). Genomic regions associated with salinity tolerance were detected on chromosomes 1A, 1D, 4B and 5A for the subtraits of relative and absolute growth rate (RGR, AGR respectively), and on chromosome 2A, 2B, 4D and 5D for Na+, potassium (K+) and chloride (Cl-) accumulation. Candidate genes that encode proteins associated with salinity tolerance were identified within the loci including Na+/H+ antiporters, K+ channels, H+-ATPase, calcineurin B-like proteins (CBLs), CBL-interacting protein kinases (CIPKs), calcium dependent protein kinases (CDPKs) and calcium-transporting ATPase. This study provides a new insight into the genetic control of salinity tolerance in a Na+ accumulating bread wheat to assist with the future development of salt tolerant cultivars.

DArT-based evaluation of soybean germplasm from Polish Gene Bank.

OBJECTIVE: Soybean is an important plant used for food, feed and many industrial purposes. Interest in soybean breeding is growing in Central Europe, including Poland. A very large number of soybean accessions are stored in gene banks, but less than 1% of them have been used for breeding. Here, we present genotypic data as well as phenotypic data on plant and seed performance, including seed chlorophyll fluorescence traits, and on yield components within a collection of soybean accessions that are conserved in the Polish Gene Bank at the Plant Breeding and Acclimatization Institute-National Research Institute. RESULTS: The materials used consisted of sub-collections: 79 Polish genotypes, including old traditional cultivars, 24 Canadian, 21 American, 21 Swedish and 31 from Central and Eastern European Countries, 9 from France and 6 from Japan. In total, 9602 high quality SNPs were derived from DArTseq, a method utilising GBS technology. GWAS, performed with the BLINK model, revealed that a total of 41 significant SNPs were mapped for days to flowering, flower colour, plant height, days to pod formation, 100 seed weight, pod colour, seeds and hilum colour and steady-state chlorophyll fluorescence under light (Ft_Lss). This is the first report about the diversity of traditional old Polish soybean cultivars.

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.

gRNA validation for wheat genome editing with the CRISPR-Cas9 system.

BACKGROUND: The CRISPR-Cas9 system is a powerful and versatile tool for crop genome editing. However, achieving highly efficient and specific editing in polyploid species can be a challenge. The efficiency and specificity of the CRISPR-Cas9 system depends critically on the gRNA used. Here, we assessed the activities and specificities of seven gRNAs targeting 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) in hexaploid wheat protoplasts. EPSPS is the biological target of the widely used herbicide glyphosate. RESULTS: The seven gRNAs differed substantially in their on-target activities, with mean indel frequencies ranging from 0% to approximately 20%. There was no obvious correlation between experimentally determined and in silico predicted on-target gRNA activity. The presence of a single mismatch within the seed region of the guide sequence greatly reduced but did not abolish gRNA activity, whereas the presence of an additional mismatch, or the absence of a PAM, all but abolished gRNA activity. Large insertions (≥20 bp) of DNA vector-derived sequence were detected at frequencies up to 8.5% of total indels. One of the gRNAs exhibited several properties that make it potentially suitable for the development of non-transgenic glyphosate resistant wheat. CONCLUSIONS: We have established a rapid and reliable method for gRNA validation in hexaploid wheat protoplasts. The method can be used to identify gRNAs that have favourable properties. Our approach is particularly suited to polyploid species, but should be applicable to any plant species amenable to protoplast transformation.

Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases.

Comparing newly obtained and previously known nucleotide and amino-acid sequences underpins modern biological research. BLAST is a well-established tool for such comparisons but is challenging to use on new data sets. We combined a user-centric design philosophy with sustainable software development approaches to create Sequenceserver, a tool for running BLAST and visually inspecting BLAST results for biological interpretation. Sequenceserver uses simple algorithms to prevent potential analysis errors and provides flexible text-based and visual outputs to support researcher productivity. Our software can be rapidly installed for use by individuals or on shared servers.

CRISPR/Cas9-mediated knockout of Ms1 enables the rapid generation of male-sterile hexaploid wheat lines for use in hybrid seed production.

The development and adoption of hybrid seed technology have led to dramatic increases in agricultural productivity. However, it has been a challenge to develop a commercially viable platform for the production of hybrid wheat (Triticum aestivum) seed due to wheat's strong inbreeding habit. Recently, a novel platform for commercial hybrid seed production was described. This hybridization platform utilizes nuclear male sterility to force outcrossing and has been applied to maize and rice. With the recent molecular identification of the wheat male fertility gene Ms1, it is now possible to extend the use of this novel hybridization platform to wheat. In this report, we used the CRISPR/Cas9 system to generate heritable, targeted mutations in Ms1. The introduction of biallelic frameshift mutations into Ms1 resulted in complete male sterility in wheat cultivars Fielder and Gladius, and several of the selected male-sterile lines were potentially non-transgenic. Our study demonstrates the utility of the CRISPR/Cas9 system for the rapid generation of male sterility in commercial wheat cultivars. This represents an important step towards capturing heterosis to improve wheat yields, through the production and use of hybrid seed on an industrial scale.

Effects of Rht-B1 and Ppd-D1 loci on pollinator traits in wheat.

KEY MESSAGE: Elite wheat pollinators are critical for successful hybrid breeding. We identified Rht-B1 and Ppd-D1 loci affecting multiple pollinator traits and therefore represent major targets for improving hybrid seed production. Hybrid breeding has a great potential to significantly boost wheat yields. Ideal male pollinators would be taller in stature, contain many spikelets well-spaced along the spike and exhibit high extrusion of large anthers. Most importantly, flowering time would match with that of the female parent. Available genetic resources for developing an elite wheat pollinator are limited, and the genetic basis for many of these traits is largely unknown. Here, we report on the genetic analysis of pollinator traits using biparental mapping populations. We identified two anther extrusion QTLs of medium effect, one on chromosome 1BL and the other on 4BS coinciding with the semi-dwarfing Rht-B1 locus. The effect of Rht-B1 alleles on anther extrusion is genotype dependent, while tall plant Rht-B1a allele is consistently associated with large anthers. Multiple QTLs were identified at the Ppd-D1 locus for anther length, spikelet number and spike length, with the photoperiod-sensitive Ppd-D1b allele associated with favourable pollinator traits in the populations studied. We also demonstrated that homeoloci, Rht-D1 and Ppd-B1, influence anther length among other traits. These results suggest that combinations of Rht-B1 and Ppd-D1 alleles control multiple pollinator traits and should be major targets of hybrid wheat breeding programs.

DAWN: a resource for yielding insights into the diversity among wheat genomes.

BACKGROUND: Democratising the growing body of whole genome sequencing data available for Triticum aestivum (bread wheat) has been impeded by the lack of a genome reference and the large computational requirements for analysing these data sets. RESULTS: DAWN (Diversity Among Wheat geNomes) integrates data from the T. aestivum Chinese Spring (CS) IWGSC RefSeq v1.0 genome with public WGS and exome data from 17 and 62 accessions respectively, enabling researchers and breeders alike to investigate genotypic differences between wheat accessions at the level of whole chromosomes down to individual genes. CONCLUSIONS: Using DAWN we show that it is possible to visualise small and large chromosomal deletions, identify haplotypes at a glance and spot the consequences of selective breeding. DAWN allows us to detect the break points of alien introgression segments brought into an accession when transferring desired genes. Furthermore, we can find possible explanations for reduced recombination in parts of a chromosome, we can predict regions with linkage drag, and also look at diversity in centromeric regions.

POTAGE: A Visualisation Tool for Speeding up Gene Discovery in Wheat.

POPSEQ Ordered Triticum aestivum Gene Expression (POTAGE) is a web application which accelerates the process of identifying candidate genes for quantitative trait loci (QTL) in hexaploid wheat. This is achieved by leveraging several of the most commonly used data sets in wheat research. These include the Chromosome Survey Sequences, their order along the chromosomes determined by the population sequencing (POPSEQ) approach, the gene predictions and RNA-Seq expression data. POTAGE aggregates those data sets and provides an intuitive interface for biologists to explore the expression of the predicted genes and their functional annotation in a chromosomal context. The interface accelerates some of the laborious and repetitive tasks commonly undertaken in the process of identifying and prioritising genes which may underlie QTL. We illustrate the utility of POTAGE by showing how a short-list of candidate genes can quickly be identified for a QTL linked to pre-harvest sprouting - a major cause of quality and yield loss in wheat production. The candidate genes identified using POTAGE included TaMKK3, which was recently reported as a causal gene for seed dormancy in wheat, and a mutation in its barley ortholog has been shown to reduce pre-harvest sprouting. POTAGE is available at http://crobiad.agwine.adelaide.edu.au/potage .

Four simple recommendations to encourage best practices in research software.

Scientific research relies on computer software, yet software is not always developed following practices that ensure its quality and sustainability. This manuscript does not aim to propose new software development best practices, but rather to provide simple recommendations that encourage the adoption of existing best practices. Software development best practices promote better quality software, and better quality software improves the reproducibility and reusability of research. These recommendations are designed around Open Source values, and provide practical suggestions that contribute to making research software and its source code more discoverable, reusable and transparent. This manuscript is aimed at developers, but also at organisations, projects, journals and funders that can increase the quality and sustainability of research software by encouraging the adoption of these recommendations.

Best practice data life cycle approaches for the life sciences.

Throughout history, the life sciences have been revolutionised by technological advances; in our era this is manifested by advances in instrumentation for data generation, and consequently researchers now routinely handle large amounts of heterogeneous data in digital formats. The simultaneous transitions towards biology as a data science and towards a 'life cycle' view of research data pose new challenges. Researchers face a bewildering landscape of data management requirements, recommendations and regulations, without necessarily being able to access data management training or possessing a clear understanding of practical approaches that can assist in data management in their particular research domain. Here we provide an overview of best practice data life cycle approaches for researchers in the life sciences/bioinformatics space with a particular focus on 'omics' datasets and computer-based data processing and analysis. We discuss the different stages of the data life cycle and provide practical suggestions for useful tools and resources to improve data management practices.

Towards an open, collaborative, reusable framework for sharing hands-on bioinformatics training workshops.

There is a clear demand for hands-on bioinformatics training. The development of bioinformatics workshop content is both time-consuming and expensive. Therefore, enabling trainers to develop bioinformatics workshops in a way that facilitates reuse is becoming increasingly important. The most widespread practice for sharing workshop content is through making PDF, PowerPoint and Word documents available online. While this effort is to be commended, such content is usually not so easy to reuse or repurpose and does not capture all the information required for a third party to rerun a workshop. We present an open, collaborative framework for developing and maintaining, reusable and shareable hands-on training workshop content.

Development of a cloud-based Bioinformatics Training Platform.

The Bioinformatics Training Platform (BTP) has been developed to provide access to the computational infrastructure required to deliver sophisticated hands-on bioinformatics training courses. The BTP is a cloud-based solution that is in active use for delivering next-generation sequencing training to Australian researchers at geographically dispersed locations. The BTP was built to provide an easy, accessible, consistent and cost-effective approach to delivering workshops at host universities and organizations with a high demand for bioinformatics training but lacking the dedicated bioinformatics training suites required. To support broad uptake of the BTP, the platform has been made compatible with multiple cloud infrastructures. The BTP is an open-source and open-access resource. To date, 20 training workshops have been delivered to over 700 trainees at over 10 venues across Australia using the BTP.

The evolutionary origin of CIPK16: A gene involved in enhanced salt tolerance.

Calcineurin B-like protein interacting protein kinases (CIPKs) are key regulators of pre-transcriptional and post-translational responses to abiotic stress. Arabidopsis thaliana CIPK16 (AtCIPK16) was identified from a forward genetic screen as a gene that mediates lower shoot salt accumulation and improved salinity tolerance in Arabidopsis and transgenic barley. Here, we aimed to gain an understanding of the evolution of AtCIPK16, and orthologues of CIPK16 in other plant species including barley, by conducting a phylogenetic analysis of terrestrial plant species. The resulting protein sequence based phylogenetic trees revealed a single clade that included AtCIPK16 along with two segmentally duplicated CIPKs, AtCIPK5 and AtCIPK25. No monocots had proteins that fell into this clade; instead the most closely related monocot proteins formed a group basal to the entire CIPK16, 5 and 25 clade. We also found that AtCIPK16 contains a core Brassicales specific indel and a putative nuclear localisation signal, which are synapomorphic characters of CIPK16 genes. In addition, we present a model that proposes the evolution of CIPK16, 5 and 25 clade.