Search

Permanent heterozygous loci, such as sex- or mating-compatibility regions, often display suppression of recombination and signals of genomic degeneration. In Basidiomycota, two distinct loci confer mating compatibility. These loci encode homeodomain (HD) transcription factors and pheromone receptor (Pra)-ligand allele pairs. To date, an analysis of genome level mating-type (MAT) loci is lacking for obligate biotrophic basidiomycetes in the Pucciniales, an order containing serious agricultural plant pathogens. Here, we focus on four species of Puccinia that infect oat and wheat, including P. coronata f. sp. avenae, P. graminis f. sp. tritici, P. triticina and P. striiformis f. sp. tritici. MAT loci are located on two separate chromosomes supporting previous hypotheses of a tetrapolar mating compatibility system in the Pucciniales. The HD genes are multiallelic in all four species while the PR locus appears biallelic, except for P. graminis f. sp. tritici, which potentially has multiple alleles. HD loci are largely conserved in their macrosynteny, both within and between species, without strong signals of recombination suppression. Regions proximal to the PR locus, however, displayed signs of recombination suppression and genomic degeneration in the three species with a biallelic PR locus. Our observations support a link between recombination suppression, genomic degeneration, and allele diversity of MAT loci that is consistent with recent mathematical modelling and simulations. Finally, we confirm that MAT genes are expressed during the asexual infection cycle, and we propose that this may support regulating nuclear maintenance and pairing during infection and spore formation. Our study provides insights into the evolution of MAT loci of key pathogenic Puccinia species. Understanding mating compatibility can help predict possible combinations of nuclear pairs, generated by sexual reproduction or somatic recombination, and the potential evolution of new virulent isolates of these important plant pathogens.

Abstract Background Long-read sequencing platforms offered by Oxford Nanopore Technologies (ONT) allow native DNA containing epigenetic modifications to be directly sequenced, but can be limited by lower per-base accuracies. A key step post-sequencing is basecalling, the process of converting raw electrical signals produced by the sequencing device into nucleotide sequences. This is challenging as current basecallers are primarily based on mixtures of model species for training. Here we utilise both ONT PromethION and higher accuracy PacBio Sequel II HiFi sequencing on two plants, Phebalium stellatum and Xanthorrhoea johnsonii, to train species-specific basecaller models with the aim of improving per-base accuracy. We investigate sequencing accuracies achieved by ONT basecallers and assess accuracy gains by training single-species and species-specific basecaller models. We also evaluate accuracy gains from ONT’s improved flowcells (R10.4, FLO-PRO112) and sequencing kits (SQK-LSK112). For the truth dataset for both model training and accuracy assessment, we developed highly accurate, contiguous diploid reference genomes with PacBio Sequel II HiFi reads. Results Basecalling with ONT Guppy 5 and 6 super-accurate gave almost identical results, attaining read accuracies of 91.96% and 94.15%. Guppy’s plant-specific model gave highly mixed results, attaining read accuracies of 91.47% and 96.18%. Species-specific basecalling models improved read accuracy, attaining 93.24% and 95.16% read accuracies. R10.4 sequencing kits also improve sequencing accuracy, attaining read accuracies of 95.46% (super-accurate) and 96.87% (species-specific). Conclusions The use of a single mixed-species basecaller model, such as ONT Guppy super-accurate, may be reducing the accuracy of nanopore sequencing, due to conflicting genome biology within the training dataset and study species. Training of single-species and genome-specific basecaller models improves read accuracy. Studies that aim to do large-scale long-read genotyping would primarily benefit from training their own basecalling models. Such studies could use sequencing accuracy gains and improving bioinformatics tools to improve study outcomes.

Abstract Infectious fungal diseases can have devastating effects on wildlife health and a detailed understanding of the evolution of related emerging fungal pathogen along with the ability to detect them in the wild is considered indispensable for effective management strategies. Several fungi from the genera Nannizziopsis and Paranannizziopsis are emerging pathogens of reptiles and have been observed to cause disease in a wide range of taxa. Nannizziopsis barbatae has become a particularly important pathogen of Australian reptiles with an increasing number of herpetofauna being reported with cases of infection from across the country. Here, we present the mitochondrial genome sequences and phylogenetic analysis for seven species in this group of fungi uncovering new information on the evolutionary relationship of these emerging pathogens. From this analysis, we designed a species‐specific qPCR assay for the rapid detection of N. barbatae and demonstrate its application in a wild urban population of a dragon lizard.

Abstract Background Most animals and plants have more than one set of chromosomes and package these haplotypes into a single nucleus within each cell. In contrast, many fungal species carry multiple haploid nuclei per cell. Rust fungi are such species with two nuclei (karyons) that contain a full set of haploid chromosomes each. The physical separation of haplotypes in dikaryons means that, unlike in diploids, Hi-C chromatin contacts between haplotypes are false-positive signals. Results We generate the first chromosome-scale, fully-phased assembly for the dikaryotic leaf rust fungus Puccinia triticina and compare Nanopore MinION and PacBio HiFi sequence-based assemblies. We show that false-positive Hi-C contacts between haplotypes are predominantly caused by phase switches rather than by collapsed regions or Hi-C read mis-mappings. We introduce a method for phasing of dikaryotic genomes into the two haplotypes using Hi-C contact graphs, including a phase switch correction step. In the HiFi assembly, relatively few phase switches occur, and these are predominantly located at haplotig boundaries and can be readily corrected. In contrast, phase switches are widespread throughout the Nanopore assembly. We show that haploid genome read coverage of 30–40 times using HiFi sequencing is required for phasing of the leaf rust genome, with 0.7% heterozygosity, and that HiFi sequencing resolves genomic regions with low heterozygosity that are otherwise collapsed in the Nanopore assembly. Conclusions This first Hi-C based phasing pipeline for dikaryons and comparison of long-read sequencing technologies will inform future genome assembly and haplotype phasing projects in other non-haploid organisms.

Massively parallel, second-generation short-read DNA sequencing has become an integral tool in biology for genomic studies. Offering highly accurate base-pair resolution at the most competitive price, the technology has become widespread. However, high-throughput generation of multiplexed DNA libraries can be costly and cumbersome. Here, we present a cost-conscious protocol for generating multiplexed short-read DNA libraries using a bead-linked transposome from Illumina. We prepare libraries in high-throughput with small reaction volumes that use 1/50th the amount of transposome compared to Illumina DNA Prep tagmentation protocols. By reducing transposome usage and optimising the protocol to circumvent magnetic bead-based clean-ups between steps, we reduce costs, labour time and DNA input requirements. Developing our own dual index primers further reduced costs and enables up to nine 96-well microplate combinations. This facilitates efficient usage of large-scale sequencing platforms, such as the Illumina NovaSeq 6000, which offers up to three terabases of sequencing per S4 flow cell. The protocol presented substantially reduces the cost per library by approximately 1/20th compared to conventional Illumina methods.

Abstract Background Silencing of transposable elements (TEs) is essential for maintaining genome stability. Plants use small RNAs (sRNAs) to direct DNA methylation to TEs (RNA-directed DNA methylation; RdDM). Similar mechanisms of epigenetic silencing in the fungal kingdom have remained elusive. Results We use sRNA sequencing and methylation data to gain insight into epigenetics in the dikaryotic fungus Puccinia graminis f. sp. tritici (Pgt), which causes the devastating stem rust disease on wheat. We use Hi-C data to define the Pgt centromeres and show that they are repeat-rich regions (~250 kb) that are highly diverse in sequence between haplotypes and, like in plants, are enriched for young TEs. DNA cytosine methylation is particularly active at centromeres but also associated with genome-wide control of young TE insertions. Strikingly, over 90% of Pgt sRNAs and several RNAi genes are differentially expressed during infection. Pgt induces waves of functionally diversified sRNAs during infection. The early wave sRNAs are predominantly 21 nts with a 5′ uracil derived from genes. In contrast, the late wave sRNAs are mainly 22-nt sRNAs with a 5′ adenine and are strongly induced from centromeric regions. TEs that overlap with late wave sRNAs are more likely to be methylated, both inside and outside the centromeres, and methylated TEs exhibit a silencing effect on nearby genes. Conclusions We conclude that rust fungi use an epigenetic silencing pathway that might have similarity with RdDM in plants. The Pgt RNAi machinery and sRNAs are under tight temporal control throughout infection and might ensure genome stability during sporulation.

Early career researchers (ECRs) are important stakeholders leading efforts to catalyze systemic change in research culture and practice. Here, we summarize the outputs from a virtual unconventional conference (unconference), which brought together 54 invited experts from 20 countries with extensive experience in ECR initiatives designed to improve the culture and practice of science. Together, we drafted 2 sets of recommendations for (1) ECRs directly involved in initiatives or activities to change research culture and practice; and (2) stakeholders who wish to support ECRs in these efforts. Importantly, these points apply to ECRs working to promote change on a systemic level, not only those improving aspects of their own work. In both sets of recommendations, we underline the importance of incentivizing and providing time and resources for systems-level science improvement activities, including ECRs in organizational decision-making processes, and working to dismantle structural barriers to participation for marginalized groups. We further highlight obstacles that ECRs face when working to promote reform, as well as proposed solutions and examples of current best practices. The abstract and recommendations for stakeholders are available in Dutch, German, Greek (abstract only), Italian, Japanese, Polish, Portuguese, Spanish, and Serbian.

ABSTRACT The kingdom Fungi is highly diverse in morphology and ecosystem function. Yet fungi are challenging to characterize as they can be difficult to culture and morphologically indistinct. Overall, their description and analysis lag far behind other microbes such as bacteria. Classification of species via high-throughput sequencing is increasingly becoming the norm for pathogen detection, microbiome studies, and environmental monitoring. With the rapid development of sequencing technologies, however, standardized procedures for taxonomic assignment of long sequence reads have not yet been well established. Focusing on nanopore sequencing technology, we compared classification and community composition analysis pipelines using shotgun and amplicon sequencing data generated from mock communities comprising 43 fungal species. We show that regardless of the sequencing methodology used, the highest accuracy of species identification was achieved by sequence alignment against a fungal-specific database. During the assessment of classification algorithms, we found that applying cutoffs to the query coverage of each read or contig significantly improved the classification accuracy and community composition analysis without major data loss. We also generated draft genome assemblies for three fungal species from nanopore data which were absent from genome databases. Our study improves sequence-based classification and estimation of relative sequence abundance using real fungal community data and provides a practical guide for the design of metagenomics analyses focusing on fungi. IMPORTANCE Our study is unique in that it provides an in-depth comparative study of a real-life complex fungal community analyzed with multiple long- and short-read sequencing approaches. These technologies and their application are currently of great interest to diverse biologists as they seek to characterize the community compositions of microbiomes. Although great progress has been made on bacterial community compositions, microbial eukaryotes such as fungi clearly lag behind. Our study provides a detailed breakdown of strategies to improve species identification with immediate relevance to real-world studies. We find that real-life data sets do not always behave as expected, distinct from reports based on simulated data sets.

Identification of the causative infectious agent is essential in the management of infectious diseases, with the ideal diagnostic method being rapid, accurate, and informative, while remaining cost-effective. Traditional diagnostic techniques rely on culturing and cell propagation to isolate and identify the causative pathogen. These techniques are limited by the ability and the time required to grow or propagate an agent in vitro and the facts that identification based on morphological traits are non-specific, insensitive, and reliant on technical expertise. The evolution of next-generation sequencing has revolutionized genomic studies to generate more data at a cheaper cost. These are divided into short- and long-read sequencing technologies, depending on the length of reads generated during sequencing runs. Long-read sequencing also called third-generation sequencing emerged commercially through the instruments released by Pacific Biosciences and Oxford Nanopore Technologies, although relying on different sequencing chemistries, with the first one being more accurate both platforms can generate ultra-long sequence reads. Long-read sequencing is capable of entirely spanning previously established genomic identification regions or potentially small whole genomes, drastically improving the accuracy of the identification of pathogens directly from clinical samples. Long-read sequencing may also provide additional important clinical information, such as antimicrobial resistance profiles and epidemiological data from a single sequencing run. While initial applications of long-read sequencing in clinical diagnosis showed that it could be a promising diagnostic technique, it also has highlighted the need for further optimization. In this review, we show the potential long-read sequencing has in clinical diagnosis of fungal infections and discuss the pros and cons of its implementation.

Summary: Background: A cornerstone of Australia's ability to control COVID-19 has been effective border control with an extensive supervised quarantine programme. However, a rapid recrudescence of COVID-19 was observed in the state of Victoria in June, 2020. We aim to describe the genomic findings that located the source of this second wave and show the role of genomic epidemiology in the successful elimination of COVID-19 for a second time in Australia. Methods: In this observational, genomic epidemiological study, we did genomic sequencing of all laboratory-confirmed cases of COVID-19 diagnosed in Victoria, Australia between Jan 25, 2020, and Jan 31, 2021. We did phylogenetic analyses, genomic cluster discovery, and integrated results with epidemiological data (detailed information on demographics, risk factors, and exposure) collected via interview by the Victorian Government Department of Health. Genomic transmission networks were used to group multiple genomic clusters when epidemiological and genomic data suggested they arose from a single importation event and diversified within Victoria. To identify transmission of emergent lineages between Victoria and other states or territories in Australia, all publicly available SARS-CoV-2 sequences uploaded before Feb 11, 2021, were obtained from the national sequence sharing programme AusTrakka, and epidemiological data were obtained from the submitting laboratories. We did phylodynamic analyses to estimate the growth rate, doubling time, and number of days from the first local infection to the collection of the first sequenced genome for the dominant local cluster, and compared our growth estimates to previously published estimates from a similar growth phase of lineage B.1.1.7 (also known as the Alpha variant) in the UK. Findings: Between Jan 25, 2020, and Jan 31, 2021, there were 20 451 laboratory-confirmed cases of COVID-19 in Victoria, Australia, of which 15 431 were submitted for sequencing, and 11 711 met all quality control metrics and were included in our analysis. We identified 595 genomic clusters, with a median of five cases per cluster (IQR 2–11). Overall, samples from 11 503 (98·2%) of 11 711 cases clustered with another sample in Victoria, either within a genomic cluster or transmission network. Genomic analysis revealed that 10 426 cases, including 10 416 (98·4%) of 10 584 locally acquired cases, diagnosed during the second wave (between June and October, 2020) were derived from a single incursion from hotel quarantine, with the outbreak lineage (transmission network G, lineage D.2) rapidly detected in other Australian states and territories. Phylodynamic analyses indicated that the epidemic growth rate of the outbreak lineage in Victoria during the initial growth phase (samples collected between June 4 and July 9, 2020; 47·4 putative transmission events, per branch, per year [1/years; 95% credible interval 26·0–85·0]), was similar to that of other reported variants, such as B.1.1.7 in the UK (mean approximately 71·5 1/years). Strict interventions were implemented, and the outbreak lineage has not been detected in Australia since Oct 29, 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread. Interpretation: Our study highlights how rapid escalation of clonal outbreaks can occur from a single incursion. However, strict quarantine measures and decisive public health responses to emergent cases are effective, even with high epidemic growth rates. Real-time genomic surveillance can alter the way in which public health agencies view and respond to COVID-19 outbreaks. Funding: The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

Open and reproducible research practices increase the reusability and impact of scientific research. The reproducibility of research results is influenced by many factors, most of which can be addressed by improved education and training. Here we describe how workshops developed by the Reproducibility for Everyone (R4E) initiative can be customized to provide researchers at all career stages and across most disciplines with education and training in reproducible research practices. The R4E initiative, which is led by volunteers, has reached more than 3000 researchers worldwide to date, and all workshop materials, including accompanying resources, are available under a CC-BY 4.0 license at https://www.repro4everyone.org/.

Macrophomina phaseolina is a soil-borne phytopathogenic fungus that causes charcoal rot in several plant species, including sorghum. We constructed a draft genome of M. phaseolina isolate BRIP 70780a from sorghum, using long-read native DNA from MinION sequencing, which was error-corrected using short-read Illumina MiSeq reads. The draft genome, consisting of 22 contigs with an N50 of 4,257,441 bp, 99.3% complete benchmarking universal single-copy orthologs, and 14,471 genes, is a valuable resource to aid future studies in population genomics and molecular diagnostic marker development for rapid detection of the pathogen.

Strain Ug99 of the wheat stem rust pathogen, Puccinia graminis f. sp. tritici, is a major threat to global food security. Here, the authors present genomic analyses supporting that Ug99 arose as a result of non-sexual genetic exchange between dikaryotic ancestors.

Fungal diseases of plants are responsible for major losses in agriculture, highlighting the need for rapid and accurate identification of plant pathogens. Disease outcomes are often defined not only by the main pathogen but are influenced by diverse microbial communities known as the microbiome at sites of infection. Here we present the first use of whole genome shot-gun sequencing with a portable DNA sequencing device as a method for the detection of fungal pathogens from wheat (Triticum aestivum) in a standard molecular biology laboratory. The data revealed that our method is robust and applicable to the diagnosis of fungal diseases including wheat stripe rust (caused by Puccinia striiformis f. sp. tritici), Septoria tritici blotch (caused by Zymoseptoria tritici), and yellow leaf spot (caused by Pyrenophora tritici repentis). We also identified the bacterial genus Pseudomonas co-present with Puccinia and Zymoseptoria but not Pyrenophora infections. One limitation of the method is the over-representation of redundant wheat genome sequences from samples. This could be addressed by long-range amplicon-based sequencing approaches in future studies, which specifically target nonhost organisms. Our work outlines a new approach for detection of a broad range of plant pathogens and associated microbes using a portable sequencer in a standard laboratory, providing the basis for future development of an on-site disease monitoring system.

Abstract Background Chloroplasts are organelles that conduct photosynthesis in plant and algal cells. The information chloroplast genome contained is widely used in agriculture and studies of evolution and ecology. Correctly assembling chloroplast genomes can be challenging because the chloroplast genome contains a pair of long inverted repeats (10–30 kb). Typically, it is simply assumed that the gross structure of the chloroplast genome matches the most commonly observed structure of two single-copy regions separated by a pair of inverted repeats. The advent of long-read sequencing technologies should remove the need to make this assumption by providing sufficient information to completely span the inverted repeat regions. Yet, long-reads tend to have higher error rates than short-reads, and relatively little is known about the best way to combine long- and short-reads to obtain the most accurate chloroplast genome assemblies. Using Eucalyptus pauciflora, the snow gum, as a test case, we evaluated the effect of multiple parameters, such as different coverage of long-(Oxford nanopore) and short-(Illumina) reads, different long-read lengths, different assembly pipelines, with a view to determining the most accurate and efficient approach to chloroplast genome assembly. Results Hybrid assemblies combining at least 20x coverage of both long-reads and short-reads generated a single contig spanning the entire chloroplast genome with few or no detectable errors. Short-read-only assemblies generated three contigs (the long single copy, short single copy and inverted repeat regions) of the chloroplast genome. These contigs contained few single-base errors but tended to exclude several bases at the beginning or end of each contig. Long-read-only assemblies tended to create multiple contigs with a much higher single-base error rate. The chloroplast genome of Eucalyptus pauciflora is 159,942 bp, contains 131 genes of known function. Conclusions Our results suggest that very accurate assemblies of chloroplast genomes can be achieved using a combination of at least 20x coverage of long- and short-reads respectively, provided that the long-reads contain at least ~5x coverage of reads longer than the inverted repeat region. We show that further increases in coverage give little or no improvement in accuracy, and that hybrid assemblies are more accurate than long-read-only or short-read-only assemblies.

Rapid advancements in long-read sequencing technologies have transformed read lengths from bps to Mbps, which has enabled chromosome-scale genome assemblies. However, read lengths are now becoming limited by the extraction of pure high-molecular weight DNA suitable for long-read sequencing, which is particularly challenging in plants and fungi. To overcome this, we present a protocol collection; high-molecular weight DNA extraction, clean-up and size selection for long-read sequencing. We optimised a gentle magnetic bead based high-molecular weight DNA extraction, which is presented here in detail. The protocol circumvents spin columns and high-centrifugation, to limit DNA fragmentation. The protocol is scalable based on tissue input, which can be used on many species of plants, fungi, reptiles and bacteria. It is also cost effective compared to kit-based protocols and hence applicable at scale in low resource settings. An optional sorbitol wash is listed and is highly recommended for plant and fungal tissues. To further remove any remaining contaminants such as phenols and polysaccharides, optional DNA clean-up and size selection strategies are given. This protocol collection is suitable for all common long-read sequencing platforms, such as technologies offered by PacBio and Oxford Nanopore. Using these protocols, sequencing on the Oxford Nanopore MinION can achieve read length N50 values of 30-50 kb, with reads exceeding 200 kb and outputs ranging from 15-30 Gbp. This has been routinely achieved with various plant, fungi, animal and bacteria samples.

Abstract Austropuccinia psidiiAustropuccinia psidiiA. psidiiA. psidii

Background European brown hares (Lepus europaeus) and European rabbits (Oryctolagus cuniculus) are invasive pest species in Australia, with rabbits having a substantially larger environmental impact than hares. As their spatial distribution in Australia partially overlaps, we conducted a comparative microbiome study to determine how the composition of gastrointestinal microbiota varies between these species, since this may indicate species differences in diet, physiology, and other internal and external factors. Methods We analysed the faecal microbiome of nine wild hares and twelve wild rabbits from a sympatric periurban reserve in Canberra, Australia, using a 16S rRNA amplicon-based sequencing approach. Additionally, we compared the concordance between results from Illumina and Nanopore sequencing platforms. Results We identified significantly more variation in faecal microbiome composition between individual rabbits compared to hares, despite both species occupying a similar habitat. The faecal microbiome in both species was dominated by the phyla Firmicutes and Bacteroidetes, typical of many vertebrates. Many phyla, including Actinobacteria, Proteobacteria and Patescibacteria, were shared between rabbits and hares. In contrast, bacteria from phylum Verrucomicrobia were present only in rabbits, while phyla Lentisphaerae and Synergistetes were represented only in hares. We did not identify phylum Spirochaetes in Australian hares; this phylum was previously shown to be present at high relative abundance in European hare faecal samples. These differences in the composition of faecal microbiota may be indicative of less discriminate foraging behaviour in rabbits, which in turn may enable them to adapt quicker to new environments, and may reflect the severe environmental impacts that this species has in Australia.

The rice XA21 and XA3 pattern receptor kinases, derived from Oryza longistaminata and an Oryza. sativa japonica cultivar Wase Aikoku 3, respectively, confer resistance to strains of the Gram-negative bacterium Xanthomonas oryzae pv. oryzae (Xoo), the causal agent of rice bacterial blight disease. Previously, we showed that transfer of Xa21 to the model rice cultivar Kitaake enhances resistance to Xoo. In this manuscript we demonstrate that Kitaake expressing Xa3 confers resistance to Xoo strain PXO79 and that the stress-related marker genes PR10b and KO5 are upregulated in Xoo-infected Xa3 rice leaves. We also show that rice somatic embryogenesis receptor kinase 2 (OsSERK2) positively regulates XA3-mediated immunity in Kitaake. We found that overexpression of XA21 binding protein 15 (XB15) and XB24, two negative regulators of XA21-mediated immunity, do not affect XA3-mediated immunity in the Kitaake genetic background. Our results indicate that the rice immune receptors XA21 and XA3 employ both shared and distinct signaling components in their response to Xoo. The results are important to further understand pathogen-associated molecular pattern (PAMP)-triggered immunity in rice. Furthermore, the presence of Kitaake rice carrying Xa3 will facilitate genetic research to study the XA3-mediated immunity.

ABSTRACT Most known examples of horizontal gene transfer (HGT) between eukaryotes are ancient. These events are identified primarily using phylogenetic methods on coding regions alone. Only rarely are there examples of HGT where noncoding DNA is also reported. The gene encoding the wheat virulence protein ToxA and the surrounding 14 kb is one of these rare examples. ToxA has been horizontally transferred between three fungal wheat pathogens (Parastagonospora nodorum, Pyrenophora tritici-repentis, and Bipolaris sorokiniana) as part of a conserved ∼14 kb element which contains coding and noncoding regions. Here we used long-read sequencing to define the extent of HGT between these three fungal species. Construction of near-chromosomal-level assemblies enabled identification of terminal inverted repeats on either end of the 14 kb region, typical of a type II DNA transposon. This is the first description of ToxA with complete transposon features, which we call ToxhAT. In all three species, ToxhAT resides in a large (140-to-250 kb) transposon-rich genomic island which is absent in isolates that do not carry the gene (annotated here as toxa−). We demonstrate that the horizontal transfer of ToxhAT between P. tritici-repentis and P. nodorum occurred as part of a large (∼80 kb) HGT which is now undergoing extensive decay. In B. sorokiniana, in contrast, ToxhAT and its resident genomic island are mobile within the genome. Together, these data provide insight into the noncoding regions that facilitate HGT between eukaryotes and into the genomic processes which mask the extent of HGT between these species. IMPORTANCE This work dissects the tripartite horizontal transfer of ToxA, a gene that has a direct negative impact on global wheat yields. Defining the extent of horizontally transferred DNA is important because it can provide clues to the mechanisms that facilitate HGT. Our analysis of ToxA and its surrounding 14 kb suggests that this gene was horizontally transferred in two independent events, with one event likely facilitated by a type II DNA transposon. These horizontal transfer events are now in various processes of decay in each species due to the repeated insertion of new transposons and subsequent rounds of targeted mutation by a fungal genome defense mechanism known as repeat induced point mutation. This work highlights the role that HGT plays in the evolution of host adaptation in eukaryotic pathogens. It also increases the growing body of evidence indicating that transposons facilitate adaptive HGT events between fungi present in similar environments and hosts.

Peer-reviewed journal publication is the main means for academic researchers in the life sciences to create a permanent public record of their work. These publications are also the de facto currency for career progress, with a strong link between journal brand recognition and perceived value. The current peer-review process can lead to long delays between submission and publication, with cycles of rejection, revision, and resubmission causing redundant peer review. This situation creates unique challenges for early career researchers (ECRs), who rely heavily on timely publication of their work to gain recognition for their efforts. Today, ECRs face a changing academic landscape, including the increased interdisciplinarity of life sciences research, expansion of the researcher population, and consequent shifts in employer and funding demands. The publication of preprints, publicly available scientific manuscripts posted on dedicated preprint servers prior to journal-managed peer review, can play a key role in addressing these ECR challenges. Preprinting benefits include rapid dissemination of academic work, open access, establishing priority or concurrence, receiving feedback, and facilitating collaborations. Although there is a growing appreciation for and adoption of preprints, a minority of all articles in life sciences and medicine are preprinted. The current low rate of preprint submissions in life sciences and ECR concerns regarding preprinting need to be addressed. We provide a perspective from an interdisciplinary group of ECRs on the value of preprints and advocate their wide adoption to advance knowledge and facilitate career development.

Gram-negative bacteria continuously pinch off portions of their outer membrane, releasing membrane vesicles. These outer membrane vesicles (OMVs) are involved in multiple processes including cell-to-cell communication, biofilm formation, stress tolerance, horizontal gene transfer, and virulence. OMVs are also known modulators of the mammalian immune response. Despite the well-documented role of OMVs in mammalian-bacterial communication, their interaction with plants is not well studied. To examine whether OMVs of plant pathogens modulate the plant immune response, we purified OMVs from four different plant pathogens and used them to treat Arabidopsis thaliana. OMVs rapidly induced a reactive oxygen species burst, medium alkalinization, and defense gene expression in A. thaliana leaf discs, cell cultures, and seedlings, respectively. Western blot analysis revealed that EF-Tu is present in OMVs and that it serves as an elicitor of the plant immune response in this form. Our results further show that the immune coreceptors BAK1 and SOBIR1 mediate OMV perception and response. Taken together, our results demonstrate that plants can detect and respond to OMV-associated molecules by activation of their immune system, revealing a new facet of plant-bacterial interactions.

Rice (Oryza sativa) plants expressing the XA21 cell-surface receptor kinase are resistant to Xanthomonas oryzae pv. oryzae (Xoo) infection. We previously demonstrated that expressing a chimeric protein containing the ELONGATION FACTOR Tu RECEPTOR (EFR) ectodomain and the XA21 endodomain (EFR:XA21) in rice does not confer robust resistance to Xoo. To test if the XA21 ectodomain is required for Xoo resistance, we produced transgenic rice lines expressing a chimeric protein consisting of the XA21 ectodomain and EFR endodomain (XA21:EFR) and inoculated these lines with Xoo. We also tested if the XA21:EFR rice plants respond to a synthetic sulfated 21 amino acid derivative (RaxX21-sY) of the activator of XA21-mediated immunity, RaxX. We found that five independently transformed XA21:EFR rice lines displayed resistance to Xoo as measured by lesion length analysis, and showed that five lines share characteristic markers of the XA21 defense response (generation of reactive oxygen species and defense response gene expression) after treatment with RaxX21-sY. Our results indicate that expression of the XA21:EFR chimeric receptor in rice confers resistance to Xoo. These results suggest that the endodomain of the EFR and XA21 immune receptors are interchangeable and the XA21 ectodomain is the key determinant conferring robust resistance to Xoo.

ABSTRACT Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae, is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenae. IMPORTANCE Disease management strategies for oat crown rust are challenged by the rapid evolution of Puccinia coronata f. sp. avenae, which renders resistance genes in oat varieties ineffective. Despite the economic importance of understanding P. coronata f. sp. avenae, resources to study the molecular mechanisms underpinning pathogenicity and the emergence of new virulence traits are lacking. Such limitations are partly due to the obligate biotrophic lifestyle of P. coronata f. sp. avenae as well as the dikaryotic nature of the genome, features that are also shared with other important rust pathogens. This study reports the first release of a haplotype-phased genome assembly for a dikaryotic fungal species and demonstrates the amenability of using emerging technologies to investigate genetic diversity in populations of P. coronata f. sp. avenae.

ABSTRACT A long-standing biological question is how evolution has shaped the genomic architecture of dikaryotic fungi. To answer this, high-quality genomic resources that enable haplotype comparisons are essential. Short-read genome assemblies for dikaryotic fungi are highly fragmented and lack haplotype-specific information due to the high heterozygosity and repeat content of these genomes. Here, we present a diploid-aware assembly of the wheat stripe rust fungus Puccinia striiformis f. sp. tritici based on long reads using the FALCON-Unzip assembler. Transcriptome sequencing data sets were used to infer high-quality gene models and identify virulence genes involved in plant infection referred to as effectors. This represents the most complete Puccinia striiformis f. sp. tritici genome assembly to date (83 Mb, 156 contigs, N50 of 1.5 Mb) and provides phased haplotype information for over 92% of the genome. Comparisons of the phase blocks revealed high interhaplotype diversity of over 6%. More than 25% of all genes lack a clear allelic counterpart. When we investigated genome features that potentially promote the rapid evolution of virulence, we found that candidate effector genes are spatially associated with conserved genes commonly found in basidiomycetes. Yet, candidate effectors that lack an allelic counterpart are more distant from conserved genes than allelic candidate effectors and are less likely to be evolutionarily conserved within the P. striiformis species complex and Pucciniales. In summary, this haplotype-phased assembly enabled us to discover novel genome features of a dikaryotic plant-pathogenic fungus previously hidden in collapsed and fragmented genome assemblies. IMPORTANCE Current representations of eukaryotic microbial genomes are haploid, hiding the genomic diversity intrinsic to diploid and polyploid life forms. This hidden diversity contributes to the organism’s evolutionary potential and ability to adapt to stress conditions. Yet, it is challenging to provide haplotype-specific information at a whole-genome level. Here, we take advantage of long-read DNA sequencing technology and a tailored-assembly algorithm to disentangle the two haploid genomes of a dikaryotic pathogenic wheat rust fungus. The two genomes display high levels of nucleotide and structural variations, which lead to allelic variation and the presence of genes lacking allelic counterparts. Nonallelic candidate effector genes, which likely encode important pathogenicity factors, display distinct genome localization patterns and are less likely to be evolutionary conserved than those which are present as allelic pairs. This genomic diversity may promote rapid host adaptation and/or be related to the age of the sequenced isolate since last meiosis.

The plasma membrane-localized BRI1-ASSOCIATED KINASE1 (BAK1) functions as a co-receptor with several receptor kinases including the brassinosteroid (BR) receptor BRASSINOSTEROID-INSENSITIVE 1 (BRI1), which is involved in growth, and the receptors for bacterial flagellin and EF-Tu, FLAGELLIN-SENSING 2 (FLS2) and EF-TU RECEPTOR (EFR), respectively, which are involved in immunity. BAK1 is a dual specificity protein kinase that can autophosphorylate on serine, threonine and tyrosine residues. It was previously reported that phosphorylation of Tyr-610 in the carboxy-terminal domain of BAK1 is required for its function in BR signaling and immunity. However, the functional role of Tyr-610 in vivo has recently come under scrutiny. Therefore, we have generated new BAK1 (Y610F) transgenic plants for functional studies. We first produced transgenic Arabidopsis lines expressing BAK1 (Y610F)-Flag in the homozygous bak1-4 bkk1-1 double null background. In a complementary approach, we expressed untagged BAK1 and BAK1 (Y610F) in the bak1-4 null mutant. Neither BAK1 (Y610F) transgenic line had any obvious growth phenotype when compared to wild-type BAK1 expressed in the same background. In addition, the BAK1 (Y610F)-Flag plants responded similarly to plants expressing BAK1-Flag in terms of brassinolide (BL) inhibition of root elongation, and there were only minor changes in gene expression between the two transgenic lines as monitored by microarray analysis and quantitative real-time PCR. In terms of plant immunity, there were no significant differences between plants expressing BAK1 (Y610F)-Flag and BAK1-Flag in the growth of the non-pathogenic hrpA- mutant of Pseudomonas syringae pv. tomato DC3000. Furthermore, untagged BAK1 (Y610F) transgenic plants were as responsive as plants expressing BAK1 (in the bak1-4 background) and wild-type Col-0 plants toward treatment with the EF-Tu- and flagellin-derived peptide epitopes elf18- and flg22, respectively, as measured by reactive oxygen species production, mitogen-activated protein kinase activation, and seedling growth inhibition. These new results do not support any involvement of Tyr-610 phosphorylation in either BR or immune signaling.

Organelle-nuclear retrograde signaling regulates gene expression, but its roles in specialized cells and integration with hormonal signaling remain enigmatic. Here we show that the SAL1-PAP (3′-phosphoadenosine 5′- phosphate) retrograde pathway interacts with abscisic acid (ABA) signaling to regulate stomatal closure and seed germination in Arabidopsis. Genetically or exogenously manipulating PAP bypasses the canonical signaling components ABA Insensitive 1 (ABI1) and Open Stomata 1 (OST1); priming an alternative pathway that restores ABA-responsive gene expression, ROS bursts, ion channel function, stomatal closure and drought tolerance in ost1-2. PAP also inhibits wild type and abi1-1 seed germination by enhancing ABA sensitivity. PAP-XRN signaling interacts with ABA, ROS and Ca2+; up-regulating multiple ABA signaling components, including lowly-expressed Calcium Dependent Protein Kinases (CDPKs) capable of activating the anion channel SLAC1. Thus, PAP exhibits many secondary messenger attributes and exemplifies how retrograde signals can have broader roles in hormone signaling, allowing chloroplasts to fine-tune physiological responses.

The rice XA21 receptor kinase confers robust resistance to the bacterial pathogen Xanthomonas oryzaepv. oryzae (Xoo). We developed a detached leaf infection assay to quickly and reliably measure activation of the XA21-mediated immune response using genetic markers. We used RNA sequencing of elf18 treated EFR:XA21:GFP plants to identify candidate genes that could serve as markers for XA21 activation. From this analysis, we identified eight genes that are up-regulated in both in elf18 treated EFR:XA21:GFP rice leaves and Xoo infected XA21 rice leaves. These results provide a rapid and reliable method to assess bacterial-rice interactions.

[This corrects the article DOI: 10.1371/journal.ppat.1004809.].

Plant plasma membrane localized pattern recognition receptors (PRRs) detect extracellular pathogen-associated molecules. PRRs such as Arabidopsis EFR and rice XA21 are taxonomically restricted and are absent from most plant genomes. Here we show that rice plants expressing EFR or the chimeric receptor EFR::XA21, containing the EFR ectodomain and the XA21 intracellular domain, sense both Escherichia coli- and Xanthomonas oryzae pv. oryzae (Xoo)-derived elf18 peptides at sub-nanomolar concentrations. Treatment of EFR and EFR::XA21 rice leaf tissue with elf18 leads to MAP kinase activation, reactive oxygen production and defense gene expression. Although expression of EFR does not lead to robust enhanced resistance to fully virulent Xoo isolates, it does lead to quantitatively enhanced resistance to weakly virulent Xoo isolates. EFR interacts with OsSERK2 and the XA21 binding protein 24 (XB24), two key components of the rice XA21-mediated immune response. Rice-EFR plants silenced for OsSERK2, or overexpressing rice XB24 are compromised in elf18-induced reactive oxygen production and defense gene expression indicating that these proteins are also important for EFR-mediated signaling in transgenic rice. Taken together, our results demonstrate the potential feasibility of enhancing disease resistance in rice and possibly other monocotyledonous crop species by expression of dicotyledonous PRRs. Our results also suggest that Arabidopsis EFR utilizes at least a subset of the known endogenous rice XA21 signaling components.

Members of UAW 5810—the union for postdoctoral researchers at the University of California—describe how their union has led to improved terms and conditions for postdocs.

Pattern recognition receptors (PRRs) play an important role in detecting invading pathogens and mounting a robust defense response to restrict infection. In rice, one of the best characterized PRRs is XA21, a leucine rich repeat receptor-like kinase that confers broad-spectrum resistance to multiple strains of the bacterial pathogen Xanthomonas oryzae pv. oryzae (Xoo). In 2009 we reported that an Xoo protein, called Ax21, is secreted by a type I-secretion system and that it serves to activate XA21-mediated immunity. This report has recently been retracted. Here we present data that corrects our previous model. We first show that Ax21 secretion does not depend on the predicted type I secretion system and that it is processed by the general secretion (Sec) system. We further show that Ax21 is an outer membrane protein, secreted in association with outer membrane vesicles. Finally, we provide data showing that ax21 knockout strains do not overcome XA21-mediated immunity.

Plants rely heavily on receptor-like kinases (RLKs) for perception and integration of external and internal stimuli. The Arabidopsis regulatory leucine-rich repeat RLK (LRR-RLK) BAK1 is involved in steroid hormone responses, innate immunity, and cell death control. Here, we describe the differential regulation of three different BAK1-dependent signaling pathways by a novel allele of BAK1, bak1-5. Innate immune signaling mediated by the BAK1-dependent RKs FLS2 and EFR is severely compromised in bak1-5 mutant plants. However, bak1-5 mutants are not impaired in BR signaling or cell death control. We also show that, in contrast to the RD kinase BRI1, the non-RD kinases FLS2 and EFR have very low kinase activity, and we show that neither was able to trans-phosphorylate BAK1 in vitro. Furthermore, kinase activity for all partners is completely dispensable for the ligand-induced heteromerization of FLS2 or EFR with BAK1 in planta, revealing another pathway specific mechanistic difference. The specific suppression of FLS2- and EFR-dependent signaling in bak1-5 is not due to a differential interaction of BAK1-5 with the respective ligand-binding RK but requires BAK1-5 kinase activity. Overall our results demonstrate a phosphorylation-dependent differential control of plant growth, innate immunity, and cell death by the regulatory RLK BAK1, which may reveal key differences in the molecular mechanisms underlying the regulation of ligand-binding RD and non-RD RKs.