Evolution of the Toxb Gene in Pyrenophora tritici-repentis and Related Species

Pyrenophora tritici-repentis (tan spot) is a destructive foliar pathogen of wheat with global impact. This ascomycete fungus possesses a highly plastic open pangenome shaped by the gain and loss of effector genes. This study investigated the allelic variations in the chlorosis-encoding gene ToxB across 422 isolates representing all identified pathotypes and worldwide origins. To gain better insights into ToxB evolution, we examined its presence and variability in other Pyrenophora spp. A ToxB haplotype network was constructed, revealing the evolutionary relationships of this gene (20 haplotypes) across four Pyrenophora species. Notably, toxb, the homolog of ToxB, was detected for the first time in the barley pathogen Pyrenophora teres. The ToxB/toxb genes display evidence of selection that is characterized by loss of function, duplication, and diverse mutations. Within the ToxB/toxb open reading frame, 72 mutations were identified, including 14 synonymous, 55 nonsynonymous, and 3 indel mutations. Remarkably, a, ∼5.6-kb Copia-like retrotransposon, named Copia-1_Ptr, was found inserted in the toxb gene of a race 3 isolate. This insert disrupted the ToxB gene's function, a first case of effector gene disruption by a transposable element in P. tritici-repentis. Additionally, a microsatellite with 25 nucleotide repeats (0 to 10) in the upstream region of ToxB suggested a potential mechanism influencing ToxB expression and regulation. Exploring ToxB-like protein distribution in other ascomycetes revealed the presence of ToxB-like proteins in 19 additional species, including the Leotiomycetes class for the first time. The presence/absence pattern of ToxB-like proteins defied species relatedness compared with a phylogenetic tree, suggesting a past horizontal gene transfer event during the evolution of the ToxB gene. [Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Agriculture and Agri-Food. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A Revised Nomenclature for ToxA Haplotypes Across Multiple Fungal Species

ToxA is one of the most studied proteinaceous necrotrophic effectors produced by plant pathogens. It has been identified in four pathogens (Pyrenophora tritici-repentis, Parastagonospora nodorum, Parastagonospora pseudonodorum [formerly Parastagonospora avenaria f. sp. tritici], and Bipolaris sorokiniana) causing leaf spot diseases on cereals worldwide. To date, 24 different ToxA haplotypes have been identified. Some P. tritici-repentis and related species also express ToxB, another small protein necrotrophic effector. We present here a revised and standardized nomenclature for these effectors, which could be extended to other poly-haplotypic genes found across multiple species.

A global pangenome for the wheat fungal pathogen Pyrenophora tritici-repentis and prediction of effector protein structural homology

The adaptive potential of plant fungal pathogens is largely governed by the gene content of a species, consisting of core and accessory genes across the pathogen isolate repertoire. To approximate the complete gene repertoire of a globally significant crop fungal pathogen, a pan genomic analysis was undertaken for Pyrenophora tritici-repentis (Ptr), the causal agent of tan (or yellow) spot disease in wheat. In this study, 15 new Ptr genomes were sequenced, assembled and annotated, including isolates from three races not previously sequenced. Together with 11 previously published Ptr genomes, a pangenome for 26 Ptr isolates from Australia, Europe, North Africa and America, representing nearly all known races, revealed a conserved core-gene content of 57 % and presents a new Ptr resource for searching natural homologues (orthologues not acquired by horizontal transfer from another species) using remote protein structural homology. Here, we identify for the first time a non-synonymous mutation in the Ptr necrotrophic effector gene ToxB, multiple copies of the inactive toxb within an isolate, a distant natural Pyrenophora homologue of a known Parastagonopora nodorum necrotrophic effector (SnTox3), and clear genomic break points for the ToxA effector horizontal transfer region. This comprehensive genomic analysis of Ptr races includes nine isolates sequenced via long read technologies. Accordingly, these resources provide a more complete representation of the species, and serve as a resource to monitor variations potentially involved in pathogenicity.

Immunomic Investigation of Holocyclotoxins to Produce the First Protective Anti-Venom Vaccine Against the Australian Paralysis Tick, Ixodes holocyclus

Venom producing animals are ubiquitously disseminated among vertebrates and invertebrates such as fish, snakes, scorpions, spiders, and ticks. Of the ~890 tick species worldwide, 27 have been confirmed to cause paralysis in mammalian hosts. The Australian paralysis tick (Ixodes holocyclus) is the most potent paralyzing tick species known. It is an indigenous three host tick species that secretes potent neurotoxins known as holocyclotoxins (HTs). Holocyclotoxins cause a severe and harmful toxicosis leading to a rapid flaccid paralysis which can result in death of susceptible hosts such as dogs. Antivenins are generally polyclonal antibody treatments developed in sheep, horses or camels to administer following bites from venomous creatures. Currently, the methods to prevent or treat tick paralysis relies upon chemical acaricide preventative treatments or prompt removal of all ticks attached to the host followed by the administration of a commercial tick-antiserum (TAS) respectively. However, these methods have several drawbacks such as poor efficacies, non-standardized dosages, adverse effects and are expensive to administer. Recently the I. holocyclus tick transcriptome from salivary glands and viscera reported a large family of 19 holocyclotoxins at 38-99% peptide sequence identities. A pilot trial demonstrated that correct folding of holocyclotoxins is needed to induce protection from paralysis. The immunogenicity of the holocyclotoxins were measured using commercial tick antiserum selecting HT2, HT4, HT8 and HT11 for inclusion into the novel cocktail vaccine. A further 4 HTs (HT1, HT12, HT14 and HT17) were added to the cocktail vaccine to ensure that the sequence variation among the HT protein family was encompassed in the formulation. A second trial comparing the cocktail of 8 HTs to a placebo group demonstrated complete protection from tick challenge. Here we report the first successful anti-venom vaccine protecting dogs from tick paralysis.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

Expansion and Conservation of Biosynthetic Gene Clusters in Pathogenic Pyrenophora spp

Pyrenophora is a fungal genus responsible for a number of major cereal diseases. Although fungi produce many specialised or secondary metabolites for defence and interacting with the surrounding environment, the repertoire of specialised metabolites (SM) within Pyrenophora pathogenic species remains mostly uncharted. In this study, an in-depth comparative analysis of the P. teres f. teres, P teres f. maculata and P. tritici-repentis potential to produce SMs, based on in silico predicted biosynthetic gene clusters (BGCs), was conducted using genome assemblies from PacBio DNA reads. Conservation of BGCs between the Pyrenophora species included type I polyketide synthases, terpene synthases and the first reporting of a type III polyketide synthase in P teres f. maculata. P. teres isolates exhibited substantial expansion of non-ribosomal peptide synthases relative to P. tritici-repentis, hallmarked by the presence of tailoring cis-acting nitrogen methyltransferase domains. P. teres isolates also possessed unique non-ribosomal peptide synthase (NRPS)-indole and indole BGCs, while a P. tritici-repentis phytotoxin BGC for triticone production was absent in P. teres. These differences highlight diversification between the pathogens that reflects their different evolutionary histories, host adaption and lifestyles.

Genomic distribution of a novel Pyrenophora tritici-repentis ToxA insertion element

The ToxA effector is a major virulence gene of Pyrenophora tritici-repentis (Ptr), a necrotrophic fungus and the causal agent of tan spot disease of wheat. ToxA and co-located genes are believed to be the result of a recent horizontally transferred highly conserved 14kb region a major pathogenic event for Ptr. Since this event, monitoring isolates for pathogenic changes has become important to help understand the underlying mechanisms in play. Here we examined ToxA in 100 Ptr isolates from Australia, Europe, North and South America and the Middle East, and uncovered in isolates from Denmark, Germany and New Zealand a new variation, a novel 166 bp insertion element (PtrHp1) which can form a perfectly matched 59 bp inverted repeat hairpin structure located downstream of the ToxA coding sequence in the 3’ UTR exon. A wider examination revealed PtrHp1 elements to be distributed throughout the genome. Analysis of genomes from Australia and North America had 50–112 perfect copies that often overlap other genes. The hairpin element appears to be unique to Ptr and the lack of ancient origins in other species suggests that PtrHp1 emerged after Ptr speciation. Furthermore, the ToxA UTR insertion site is identical for different isolates, which suggests a single insertion event occurred after the ToxA horizontal transfer. In vitro and in planta-detached leaf assays found that the PtrHp1 element insertion had no effect on ToxA expression. However, variation in the expression of ToxA was detected between the Ptr isolates from different demographic locations, which appears to be unrelated to the presence of the element. We envision that this discovery may contribute towards future understanding of the possible role of hairpin elements in Ptr.

Gene-enriched draft genome of the cattle tick Rhipicephalus microplus: assembly by the hybrid Pacific Biosciences/Illumina approach enabled analysis of the highly repetitive genome

The genome of the cattle tick Rhipicephalus microplus, an ectoparasite with global distribution, is estimated to be 7.1Gbp in length and consists of approximately 70% repetitive DNA. We report the draft assembly of a tick genome that utilized a hybrid sequencing and assembly approach to capture the repetitive fractions of the genome. Our hybrid approach produced an assembly consisting of 2.0Gbp represented in 195,170 scaffolds with a N50 of 60,284bp. The Rmi v2.0 assembly is 51.46% repetitive with a large fraction of unclassified repeats, short interspersed elements, long interspersed elements and long terminal repeats. We identified 38,827 putative R. microplus gene loci, of which 24,758 were protein coding genes (≥100 amino acids). OrthoMCL comparative analysis against 11 selected species including insects and vertebrates identified 10,835 and 3,423 protein coding gene loci that are unique to R. microplus or common to both R. microplus and Ixodes scapularis ticks, respectively. We identified 191 microRNA loci, of which 168 have similarity to known miRNAs and 23 represent novel miRNA families. We identified the genomic loci of several highly divergent R. microplus esterases with sequence similarity to acetylcholinesterase. Additionally we report the finding of a novel cytochrome P450 CYP41 homolog that shows similar protein folding structures to known CYP41 proteins known to be involved in acaricide resistance.

Acetylcholinesterase 1 in populations of organophosphate-resistant North American strains of the cattle tick, Rhipicephalus microplus (Acari: Ixodidae)

Rhipicephalus microplus, the cattle fever tick, is a global economic problem to the cattle industry due to direct infestation of cattle and pathogens transmitted during feeding. Cattle fever tick outbreaks continue to occur along the Mexico-US border even though the tick has been eradicated from the USA. The organophosphate (OP) coumaphos targets acetylcholinesterase (AChE) and is the approved acaricide for eradicating cattle fever tick outbreaks. There is evidence for coumaphos resistance developing in cattle ticks in Mexico, and OP-resistant R. microplus ticks were discovered in outbreak populations of Texas in 2005. The molecular basis of coumaphos resistance is not known, and our study was established to gather further information on whether AChE1 is involved in the resistance mechanism. We also sought information on allele diversity in tick populations with different levels of coumaphos resistance. The overarching project goal was to define OP resistance-associated gene mutations such that a DNA-based diagnostic assay could be developed to assist the management of resistance. Three different AChE transcripts have been reported in R. microplus, and supporting genomic and transcriptomic data are available at CattleTickBase. Here, we report the complete R. microplus AChE1 gene ascertained by sequencing a bacterial artificial chromosome clone containing the entire coding region and the flanking 5' and 3' regions. We also report AChE1 sequences of larval ticks from R. microplus strains having different sensitivities to OP. To accomplish this, we sequenced a 669-bp region of the AChE1 gene corresponding to a 223 amino acid region of exon 2 to assess alleles in seven strains of R. microplus with varying OP resistance phenotypes. We identified 72 AChE1 sequence variants, 2 of which are strongly associated with OP-resistant phenotypes. Esterase-like sequences from the R. microplus transcriptome RmiTr Version 1.0 were compared to the available sequence databases to identify other transcripts with similarity to AChE1.

Genome sequences of six wheat-infecting fusarium species isolates

Fusarium pathogens represent a major constraint to wheat and barley production worldwide. To facilitate future comparative studies of Fusarium species that are pathogenic to wheat, the genome sequences of four Fusarium pseudograminearum isolates, a single Fusarium acuminatum isolate, and an organism from the Fusarium incarnatum-F. equiseti species complex are reported.

The complexity of Rhipicephalus (Boophilus) microplus genome characterised through detailed analysis of two BAC clones

BACKGROUND

Rhipicephalus (Boophilus) microplus (Rmi) a major cattle ectoparasite and tick borne disease vector, impacts on animal welfare and industry productivity. In arthropod research there is an absence of a complete Chelicerate genome, which includes ticks, mites, spiders, scorpions and crustaceans. Model arthropod genomes such as Drosophila and Anopheles are too taxonomically distant for a reference in tick genomic sequence analysis. This study focuses on the de-novo assembly of two R. microplus BAC sequences from the understudied R microplus genome. Based on available R. microplus sequenced resources and comparative analysis, tick genomic structure and functional predictions identify complex gene structures and genomic targets expressed during tick-cattle interaction.

RESULTS

In our BAC analyses we have assembled, using the correct positioning of BAC end sequences and transcript sequences, two challenging genomic regions. Cot DNA fractions compared to the BAC sequences confirmed a highly repetitive BAC sequence BM-012-E08 and a low repetitive BAC sequence BM-005-G14 which was gene rich and contained short interspersed elements (SINEs). Based directly on the BAC and Cot data comparisons, the genome wide frequency of the SINE Ruka element was estimated. Using a conservative approach to the assembly of the highly repetitive BM-012-E08, the sequence was de-convoluted into three repeat units, each unit containing an 18S, 5.8S and 28S ribosomal RNA (rRNA) encoding gene sequence (rDNA), related internal transcribed spacer and complex intergenic region.In the low repetitive BM-005-G14, a novel gene complex was found between to 2 genes on the same strand. Nested in the second intron of a large 9 Kb papilin gene was a helicase gene. This helicase overlapped in two exonic regions with the papilin. Both these genes were shown expressed in different tick life stage important in ectoparasite interaction with the host. Tick specific sequence differences were also determined for the papilin gene and the protein binding sites of the 18S subunit in a comparison to Bos taurus.

CONCLUSION

In the absence of a sequenced reference genome we have assembled two complex BAC sequences, characterised novel gene structure that was confirmed by gene expression and sequencing analyses. This is the first report to provide evidence for 2 eukaryotic genes with exon regions that overlap on the same strand, the first to describe Rhipicephalinae papilin, and the first to report the complete ribosomal DNA repeated unit sequence structure for ticks. The Cot data estimation of genome wide sequence frequency means this research will underpin future efforts for genome sequencing and assembly of the R. microplus genome.

Aligning a new reference genetic map of Lupinus angustifolius with the genome sequence of the model legume, Lotus japonicus

We have developed a dense reference genetic map of Lupinus angustifolius (2n = 40) based on a set of 106 publicly available recombinant inbred lines derived from a cross between domesticated and wild parental lines. The map comprised 1090 loci in 20 linkage groups and three small clusters, drawing together data from several previous mapping publications plus almost 200 new markers, of which 63 were gene-based markers. A total of 171 mainly gene-based, sequence-tagged site loci served as bridging points for comparing the Lu. angustifolius genome with the genome sequence of the model legume, Lotus japonicus via BLASTn homology searching. Comparative analysis indicated that the genomes of Lu. angustifolius and Lo. japonicus are highly diverged structurally but with significant regions of conserved synteny including the region of the Lu. angustifolius genome containing the pod-shatter resistance gene, lentus. We discuss the potential of synteny analysis for identifying candidate genes for domestication traits in Lu. angustifolius and in improving our understanding of Fabaceae genome evolution.

Genomic analysis of Campylobacter fetus subspecies: identification of candidate virulence determinants and diagnostic assay targets

Background

Campylobacter fetus subspecies venerealis is the causative agent of bovine genital campylobacteriosis, asymptomatic in bulls the disease is spread to female cattle causing extensive reproductive loss. The microbiological and molecular differentiation of C. fetus subsp. venerealis from C. fetus subsp. fetus is extremely difficult. This study describes the analysis of the available C. fetus subsp. venerealis AZUL-94 strain genome (~75–80%) to identify elements exclusively found in C. fetus subsp. venerealis strains as potential diagnostic targets and the characterisation of subspecies virulence genes.

Results

Eighty Kb of genomic sequence (22 contigs) was identified as unique to C. fetus subsp. venerealis AZUL-94 and consisted of type IV secretory pathway components, putative plasmid genes and hypothetical proteins. Of the 9 PCR assays developed to target C. fetus subsp. venerealis type IV secretion system genes, 4 of these were specific for C. fetus subsp. venerealis biovar venerealis and did not detect C. fetus subsp. venerealis biovar intermedius. Two assays were specific for C. fetus subsp. venerealis AZUL-94 strain, with a further single assay specific for the AZUL-94 strain and C. fetus subsp. venerealis biovar intermedius (and not the remaining C. fetus subsp. venerealis biovar venerealis strains tested). C. fetus subsp. fetus and C. fetus subsp. venerealis were found to share most common Campylobacter virulence factors such as SAP, chemotaxis, flagellar biosynthesis, 2-component systems and cytolethal distending toxin subunits (A, B, C). We did not however, identify in C. fetus the full complement of bacterial adherence candidates commonly found in other Campylobacter spp.

Conclusion

The comparison of the available C. fetus subsp. venerealis genome sequence with the C. fetus subsp. fetus genome identified 80 kb of unique C. fetus subsp. venerealis AZUL94 sequence, with subsequent PCR confirmation demonstrating inconsistent amplification of these targets in all other C. fetus subsp. venerealis strains and biovars tested. The assays developed here highlight the complexity of targeting strain specific virulence genes for field studies for the molecular identification and epidemiology of C. fetus.

The first gene-based map of Lupinus angustifolius L.-location of domestication genes and conserved synteny with Medicago truncatula

We report the first gene-based linkage map of Lupinus angustifolius (narrow-leafed lupin) and its comparison to the partially sequenced genome of Medicago truncatula. The map comprises 382 loci in 20 major linkage groups, two triplets, three pairs and 11 unlinked loci and is 1,846 cM in length. The map was generated from the segregation of 163 RFLP markers, 135 gene-based PCR markers, 75 AFLP and 4 AFLP-derived SCAR markers in a mapping population of 93 recombinant inbred lines, derived from a cross between domesticated and wild-type parents. This enabled the mapping of five major genes controlling key domestication traits in L. angustifolius. Using marker sequence data, the L. angustifolius genetic map was compared to the partially completed M. truncatula genome sequence. We found evidence of conserved synteny in some regions of the genome despite the wide evolutionary distance between these legume species. We also found new evidence of widespread duplication within the L. angustifolius genome.