GRAbB: Selective Assembly of Genomic Regions, a New Niche for Genomic Research

Phylogenomics of transcriptionally active AP2/ERF and bHLH transcription factors and study of their promoter regions in Nothapodytes nimmoniana (J.Graham) Mabb

Nothapodytes nimmoniana is a medicinally important plant producing anticancer monoterpene indole alkaloid (MIA), camptothecin (CPT). The CPT is synthesised through the strictosidine intermediate following the MIA pathway; however, transcriptional regulation of CPT pathway is still elusive in N. nimmoniana. Biosynthesis of MIA is regulated by various transcription factors (TFs) belonging to AP2/ERF, bHLH, MYB, and WRKY families. The present study identified transcriptionally active full-length 105 AP2/ERF and 68 bHLH family TFs from the N. nimmoniana. AP2/ERF TFs were divided into three subfamilies along with a soloist, while bHLH TFs were divided into 10 subfamilies according to their phylogenetic similarities. Three group IXa ERFs, Nn-ERF22, Nn-ERF29, and Nn-ERF41, one subfamily IVa TF Nn-bHLH7, and three subfamilies IIIe Nn-bHLH33, Nn-bHLH51, and Nn-bHLH52 clustered with the TFs regulating alkaloid biosynthesis in Catharanthus roseus, tomato, tobacco, and Artemisia annua. Expression of these TFs in N. nimmoniana was higher in roots, which is a primary CPT accumulating tissue. Moreover, genome skimming approach was used to reconstruct the promoter regions of candidate ERF genes to identify the cis-regulatory elements. The presence of G-boxes and other jasmonic acid-responsive elements in the promoter suggests the regulation of ERFs by bHLHs. The present study effectively generated and used genomics resource for characterisation of regulatory TFs from non-model medicinal plant.

Evolutionary mechanisms underlying the diversification of nuclear factor of activated T cells across vertebrates

The mechanisms of immunity linked to biological evolution are crucial for understanding animal morphogenesis, organogenesis, and biodiversity. The nuclear factor of activated T cells (NFAT) family consists of five members (NFATc1-c4, 5) with different functions in the immune system. However, the evolutionary dynamics of NFATs in vertebrates has not been explored. Herein, we investigated the origin and mechanisms underlying the diversification of NFATs by comparing the gene, transcript and protein sequences, and chromosome information. We defined an ancestral origin of NFATs during the bilaterian development, dated approximately 650 million years ago, where NFAT5 and NFATc1-c4 were derived independently. The conserved parallel evolution of NFATs in multiple species was probably attributed to their innate nature. Conversely, frequent gene duplications and chromosomal rearrangements in the recently evolved taxa have suggested their roles in the adaptive immune evolution. A significant correlation was observed between the chromosome rearrangements with gene duplications and the structural fixation changes in vertebrate NFATs, suggesting their role in NFAT diversification. Remarkably, a conserved gene structure around NFAT genes with vertebrate evolutionary-related breaking points indicated the inheritance of NFATs with their neighboring genes as a unit. The close relationship between NFAT diversification and vertebrate immune evolution was suggested.

Discovery of Early-Branching Wolbachia Reveals Functional Enrichment on Horizontally Transferred Genes

Wolbachia is a widespread endosymbiont of insects and filarial nematodes that profoundly influences host biology. Wolbachia has also been reported in rhizosphere hosts, where its diversity and function remain poorly characterized. The discovery that plant-parasitic nematodes (PPNs) host Wolbachia strains with unknown roles is of interest evolutionarily, ecologically, and for agriculture as a potential target for developing new biological controls. The goal of this study was to screen communities for PPN endosymbionts and analyze genes and genomic patterns that might indicate their role. Genome assemblies revealed 1 out of 16 sampled sites had nematode communities hosting a Wolbachia strain, designated wTex, that has highly diverged as one of the early supergroup L strains. Genome features, gene repertoires, and absence of known genes for cytoplasmic incompatibility, riboflavin, biotin, and other biosynthetic functions placed wTex between mutualist C + D strains and reproductive parasite A + B strains. Functional terms enriched in group L included protoporphyrinogen IX, thiamine, lysine, fatty acid, and cellular amino acid biosynthesis, while dN/dS analysis suggested the strongest purifying selection on arginine and lysine metabolism, and vitamin B6, heme, and zinc ion binding, suggesting these as candidate roles in PPN Wolbachia. Higher dN/dS pathways between group L, wPni from aphids, wFol from springtails, and wCfeT from cat fleas suggested distinct functional changes characterizing these early Wolbachia host transitions. PPN Wolbachia had several putative horizontally transferred genes, including a lysine biosynthesis operon like that of the mitochondrial symbiont Midichloria, a spirochete-like thiamine synthesis operon shared only with wCfeT, an ATP/ADP carrier important in Rickettsia, and a eukaryote-like gene that may mediate plant systemic acquired resistance through the lysine-to-pipecolic acid system. The Discovery of group L-like variants from global rhizosphere databases suggests diverse PPN Wolbachia strains remain to be discovered. These findings support the hypothesis of plant-specialization as key to shaping early Wolbachia evolution and present new functional hypotheses, demonstrating promise for future genomics-based rhizosphere screens.

SLAG: A Program for Seeded Local Assembly of Genes in Complex Genomes

Although finished genomes have become more common, there is still a need for assemblies of individual genes or chromosomal regions when only unassembled reads are available. slag (Seeded Local Assembly of Genes) fulfils this need by performing iterative local assembly based on cycles of matching‐read retrieval with blast and assembly with cap3, phrap, spades, canu or unicycler. The target sequence can be nucleotide or protein. Read fragmentation allows slag to use phrap or cap3 to assemble long reads at lower coverage (e.g., 5×) than is possible with canu or unicycler. In simple, nonrepetitive genomes, a slag assembly can cover a whole chromosome, but in complex genomes the growth of target‐matching contigs is limited as additional reads are consumed by consensus contigs consisting of repetitive elements. Apart from genomic complexity, contig length and correctness depend on read length and accuracy. With pyrosequencing or Illumina reads, slag‐assembled contigs are accurate enough to allow design of PCR primers, while contigs assembled from Oxford Nanopore or pre‐HiFi Pacific Biosciences long reads are generally only accurate enough to design baiting sequences for further targeted sequencing. In an application with real reads, slag successfully extended sequences for four wheat genes, which were verified by cloning and Sanger sequencing of overlapping amplicons. slag is a robust alternative to atram2 for local assemblies, especially for read sets with less than 20× coverage. slag is freely available at https://github.com/cfcrane/SLAG.

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

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

Should we hail the Red King? Evolutionary consequences of a mutualistic lifestyle in genomes of lichenized ascomycetes

The Red Queen dynamic is often brought into play for antagonistic relationships. However, the coevolutionary effects of mutualistic interactions, which predict slower evolution for interacting organisms (Red King), have been investigated to a lesser extent. Lichens are a stable, mutualistic relationship of fungi and cyanobacteria and/or algae, which originated several times independently during the evolution of fungi. Therefore, they represent a suitable system to investigate the coevolutionary effect of mutualism on the fungal genome. We measured substitution rates and selective pressure of about 2000 protein-coding genes (plus the rDNA region) in two different classes of Ascomycota, each consisting of closely related lineages of lichenized and non-lichenized fungi. Our results show that independent lichenized clades are characterized by significantly slower rates for both synonymous and non-synonymous substitutions. We hypothesize that this evolutionary pattern is connected to the lichen life cycle (longer generation time of lichenized fungi) rather than a result of different selection strengths, which is described as the main driver for the Red Kind dynamic. This first empirical evidence of slower evolution in lichens provides an important insight on how biotic cooperative interactions are able to shape the evolution of symbiotic organisms.

Mitochondrial Genome of Nonmodel Marine Metazoans by Next-Generation Sequencing (NGS)

Mitochondrial genomes (mtgenome) represent an important source of information for addressing fundamental evolutionary, phylogeographic, systematic, and ecological questions in marine organisms. In the last two decades the advent of high-throughput next-generation sequencing (NGS) has provided an unprecedented possibility to access large amount of genomic data and, as such, there has been a rapid growth in mtgenome resources and studies. In particular, NGS strategies represent a great advantage for investigating nonmodel marine organisms for which no or limited genomic resources are available. Here, we describe a routinely used standardized protocol to obtain mtgenome of nonmodel marine organisms by NGS. The protocol is composed of five main steps, including DNA extraction, DNA fragmentation, library preparation, high-throughput sequencing, and bioinformatic analyses. Each of the first three steps is followed by size/quality and concentration validations. The advantages of the described protocol rely on the assumption that no a priori information on mtgenome of the studied organism is needed and on its versatility as researchers may choose several kits for DNA extraction and library preparation and adopt different methods for DNA fragmentation depending on their needs, experience, and suppliers.

LAMP Assay for Distinguishing Fusarium oxysporum and Fusarium commune in Lotus (Nelumbo nucifera) Rhizomes

Yields of edible rhizome from cultivation of the perennial hydrophyte lotus (Nelumbo nucifera) can be severely reduced by rhizome rot disease caused by Fusarium species. There is a lack of rapid field-applicable methods for detection of these pathogens on lotus plants displaying symptoms of rhizome rot. Fusarium commune (91%) and Fusarium oxysporum (9%) were identified at different frequencies from lotus samples showing symptoms of rhizome rot. Because these two species can cause different severity of disease and their morphology is similar, molecular diagnostic-based methods to detect these two species were developed. Based on the comparison of the mitochondrial genome of the two species, three specific DNA loci targets were found. The designed primer sets for conventional PCR, quantitative PCR, and loop-mediated isothermal amplification (LAMP) precisely distinguished the above two species when isolated from lotus and other plants. The LAMP detection limits were 10 pg/μl and 1 pg/μl of total DNA for F. commune and F. oxysporum, respectively. We also carried out field-mimicked experiments on lotus seedlings and rhizomes (including inoculated samples and field-diseased samples), and the results indicated that the LAMP primer sets and the supporting portable methods are suitable for rapid diagnosis of the lotus disease in the field. The LAMP-based detection method will aid in the rapid identification of whether F. oxysporum or F. commune is infecting lotus plants with symptoms of rhizome rot and can facilitate efficient pesticide use and prevent disease spread through vegetative propagation of Fusarium-infected lotus rhizomes.

Salicylic acid fights against Fusarium wilt by inhibiting target of rapamycin signaling pathway in Fusarium oxysporum

•

Isolating and sequencing the genome of F. oxysporum from potato tubers with dry rot symptoms.

•

SA efficiently arrests hyphal growth, sporular production and pathogenicity of F. oxysporum.

•

SA inhibits the activity of FoTORC1 via activating FoSNF1 in F. oxysporum.

•

Transgenic potato plants with interference of FoTOR1 and FoSAH1 genes prevent the occurrence of Fusarium wilt.

•

Providing insights SA into controlling various fungal diseases by targeting the SNF1-TORC1 pathway of pathogens.

Introduction

Objectives

Methods

Results

Conclusion

A bioinformatic platform to integrate target capture and whole genome sequences of various read depths for phylogenomics

The increasing availability of short‐read whole genome sequencing (WGS) provides unprecedented opportunities to study ecological and evolutionary processes. Although loci of interest can be extracted from WGS data and combined with target sequence data, this requires suitable bioinformatic workflows. Here, we test different assembly and locus extraction strategies and implement them into secapr, a pipeline that processes short‐read data into multilocus alignments for phylogenetics and molecular ecology analyses. We integrate the processing of data from low‐coverage WGS (<30×) and target sequence capture into a flexible framework, while optimizing de novo contig assembly and loci extraction. Specifically, we test different assembly strategies by contrasting their ability to recover loci from targeted butterfly protein‐coding genes, using four data sets: a WGS data set across different average coverages (10×, 5× and 2×) and a data set for which these loci were enriched prior to sequencing via target sequence capture. Using the resulting de novo contigs, we account for potential errors within contigs and infer phylogenetic trees to evaluate the ability of each assembly strategy to recover species relationships. We demonstrate that choosing multiple sizes of kmer simultaneously for assembly results in the highest yield of extracted loci from de novo assembled contigs, while data sets derived from sequencing read depths as low as 5× recovers the expected species relationships in phylogenetic trees. By making the tested assembly approaches available in the secapr pipeline, we hope to inspire future studies to incorporate complementary data and make an informed choice on the optimal assembly strategy.

Insights into the genome and secretome of Fusarium metavorans DSM105788 by cultivation on agro-residual biomass and synthetic nutrient sources

BACKGROUND

The transition to a biobased economy involving the depolymerization and fermentation of renewable agro-industrial sources is a challenge that can only be met by achieving the efficient hydrolysis of biomass to monosaccharides. In nature, lignocellulosic biomass is mainly decomposed by fungi. We recently identified six efficient cellulose degraders by screening fungi from Vietnam.

RESULTS

We characterized a high-performance cellulase-producing strain, with an activity of 0.06 U/mg, which was identified as a member of the Fusarium solani species complex linkage 6 (Fusarium metavorans), isolated from mangrove wood (FW16.1, deposited as DSM105788). The genome, representing nine potential chromosomes, was sequenced using PacBio and Illumina technology. In-depth secretome analysis using six different synthetic and artificial cellulose substrates and two agro-industrial waste products identified 500 proteins, including 135 enzymes assigned to five different carbohydrate-active enzyme (CAZyme) classes. The F. metavorans enzyme cocktail was tested for saccharification activity on pre-treated sugarcane bagasse, as well as untreated sugarcane bagasse and maize leaves, where it was complemented with the commercial enzyme mixture Accellerase 1500. In the untreated sugarcane bagasse and maize leaves, initial cell wall degradation was observed in the presence of at least 196 µg/mL of the in-house cocktail. Increasing the dose to 336 µg/mL facilitated the saccharification of untreated sugarcane biomass, but had no further effect on the pre-treated biomass.

CONCLUSION

Our results show that F. metavorans DSM105788 is a promising alternative pre-treatment for the degradation of agro-industrial lignocellulosic materials. The enzyme cocktail promotes the debranching of biopolymers surrounding the cellulose fibers and releases reduced sugars without process disadvantages or loss of carbohydrates.

The complete molecular sequence of chloroplast genome of Lablab purpureus (L.) Sweet

The complete molecular sequence of chloroplast genome of Lablab purpureus (L.) Sweet was firstly assembled and characterized using Illumina sequencing technology. It is 151916 bp in length, with a GC content of 35.4%, and has a typical quadrant structure, including a large single-copy region (LSC), a pair of inverted repeat regions (IRs) and a small single-copy region (SSC), the sequence length is 81132, 53244, 17540 bp, respectively. There are 131 genes in the L. purpureus chloroplast genome, including 84 encoding protein genes, 8 rRNA genes, and 38 tRNA genes. Phylogenetic analysis showed that L. purpureus clustered into a large evolutionary clade with three Vigna species.

Impact of parasitic lifestyle and different types of centromere organization on chromosome and genome evolution in the plant genus Cuscuta

The parasitic genus Cuscuta (Convolvulaceae) is exceptional among plants with respect to centromere organization, including both monocentric and holocentric chromosomes, and substantial variation in genome size and chromosome number. We investigated 12 species representing the diversity of the genus in a phylogenetic context to reveal the molecular and evolutionary processes leading to diversification of their genomes. We measured genome sizes and investigated karyotypes and centromere organization using molecular cytogenetic techniques. We also performed low-pass whole genome sequencing and comparative analysis of repetitive DNA composition. A remarkable 102-fold variation in genome sizes (342-34 734 Mbp/1C) was detected for monocentric Cuscuta species, while genomes of holocentric species were of moderate sizes (533-1545 Mbp/1C). The genome size variation was primarily driven by the differential accumulation of LTR-retrotransposons and satellite DNA. The transition to holocentric chromosomes in the subgenus Cuscuta was associated with loss of histone H2A phosphorylation and elimination of centromeric retrotransposons. In addition, basic chromosome number of holocentric species (x = 7) was smaller than in monocentrics (x = 15 or 16). We demonstrated that the transition to holocentricity in Cuscuta was accompanied by significant changes in epigenetic marks, chromosome number and the repetitive DNA sequence composition.

Extraordinary Sequence Diversity and Promiscuity of Centromeric Satellites in the Legume Tribe Fabeae

Satellite repeats are major sequence constituents of centromeres in many plant and animal species. Within a species, a single family of satellite sequences typically occupies centromeres of all chromosomes and is absent from other parts of the genome. Due to their common origin, sequence similarities exist among the centromere-specific satellites in related species. Here, we report a remarkably different pattern of centromere evolution in the plant tribe Fabeae, which includes genera Pisum, Lathyrus, Vicia, and Lens. By immunoprecipitation of centromeric chromatin with CENH3 antibodies, we identified and characterized a large and diverse set of 64 families of centromeric satellites in 14 species. These families differed in their nucleotide sequence, monomer length (33-2,979 bp), and abundance in individual species. Most families were species-specific, and most species possessed multiple (2-12) satellites in their centromeres. Some of the repeats that were shared by several species exhibited promiscuous patterns of centromere association, being located within CENH3 chromatin in some species, but apart from the centromeres in others. Moreover, FISH experiments revealed that the same family could assume centromeric and noncentromeric positions even within a single species. Taken together, these findings suggest that Fabeae centromeres are not shaped by the coevolution of a single centromeric satellite with its interacting CENH3 proteins, as proposed by the centromere drive model. This conclusion is also supported by the absence of pervasive adaptive evolution of CENH3 sequences retrieved from Fabeae species.

Detecting Introgression Between Members of the Fusarium fujikuroi and F. oxysporum Species Complexes by Comparative Mitogenomics

The Fusarium fujikuroi species complex (FFSC) and F. oxysporum species complex (FOSC) are two related groups of plant pathogens causing a wide diversity of diseases in agricultural crops world wide. The aims of this study are (1) to clarify the phylogeny of the FFSC, (2) to identify potential deviation from tree-like evolution, (3) to explore the value of using mitogenomes for these kinds of analyses, and (4) to better understand mitogenome evolution. In total, we have sequenced 24 species from the FFSC and a representative set of recently analyzed FOSC strains was chosen, while F. redolens was used as outgroup for the two species complexes. A species tree was constructed based on the concatenated alignment of seven nuclear genes and the mitogenome, which was contrasted to individual gene trees to identify potential conflicts. These comparisons indicated conflicts especially within the previously described African clade of the FFSC. Furthermore, the analysis of the mitogenomes revealed the presence of a variant of the large variable (LV) region in FFSC which was previously only reported for FOSC. The distribution of this variant and the results of sequence comparisons indicate horizontal genetic transfer between members of the two species complexes, most probably through introgression. In addition, a duplication of atp9 was found inside an intron of cob, which suggests that even highly conserved mitochondrial genes can have paralogs. Paralogization in turn may lead to inaccurate single gene phylogenies. In conclusion, mitochondrial genomes provide a robust basis for phylogeny. Comparative phylogenetic analysis indicated that gene flow among and between members of FFSC and FOSC has played an important role in the evolutionary history of these two groups. Since mitogenomes show greater levels of conservation and synteny than nuclear regions, they are more likely to be compatible for recombination than nuclear regions. Therefore, mitogenomes can be used as indicators to detect interspecies gene flow.

Population Genomic Analysis Reveals a Highly Conserved Mitochondrial Genome in Fusarium asiaticum

Fusarium asiaticum is one of the pivotal members of the Fusarium graminearum species complex (FGSC) causing Fusarium head blight (FHB) on wheat, barley and rice in large parts of Asia. Besides resulting in yield losses, FHB also causes the accumulation of mycotoxins such as nivalenol (NIV) and deoxynivalenol (DON). The aim of this study was to conduct population studies on F. asiaticum from Southern China through mitochondrial genome analyses. All strains were isolated from wheat or rice from several geographic areas in seven provinces in Southern China. Based on geographic location and host, 210 isolates were selected for next generation sequencing, and their mitogenomes were assembled by GRAbB and annotated to explore the mitochondrial genome variability of F. asiaticum. The F. asiaticum mitogenome proves extremely conserved and variation is mainly caused by absence/presence of introns harboring homing endonuclease genes. These variations could be utilized to develop molecular markers for track and trace of migrations within and between populations. This study illustrates how mitochondrial introns can be used as markers for population genetic analysis. SNP analysis demonstrate the occurrence of mitochondrial recombination in F. asiaticum as was previously found for F. oxysporum and implied for F. graminearum. Furthermore, varying degrees of genetic diversity and recombination showed a high association with different geographic regions as well as with cropping systems. The mitogenome of F. graminearum showed a much higher SNP diversity while the interspecies intron variation showed no evidence of gene flow between the two closely related and sexual compatible species.

Direct visualization of degradation microcompartments at the ER membrane

Endoplasmic reticulum-associated degradation (ERAD) is an essential process that removes misfolded proteins from the ER, preventing cellular dysfunction and disease. While most of the key components of ERAD are known, their specific localization remains a mystery. This study uses in situ cryo-electron tomography to directly visualize the ERAD machinery within the native cellular environment. Proteasomes and Cdc48, the complexes that extract and degrade ER proteins, cluster together in non–membrane-bound cytosolic microcompartments that contact ribosome-free patches on the ER membrane. This discrete molecular organization may facilitate efficient ERAD. Structural analysis reveals that proteasomes directly engage ER-localized substrates, providing evidence for a noncanonical “direct ERAD” pathway. In addition, live-cell fluorescence microscopy suggests that these ER-associated proteasome clusters form by liquid–liquid phase separation.

To promote the biochemical reactions of life, cells can compartmentalize molecular interaction partners together within separated non–membrane-bound regions. It is unknown whether this strategy is used to facilitate protein degradation at specific locations within the cell. Leveraging in situ cryo-electron tomography to image the native molecular landscape of the unicellular alga Chlamydomonas reinhardtii, we discovered that the cytosolic protein degradation machinery is concentrated within ∼200-nm foci that contact specialized patches of endoplasmic reticulum (ER) membrane away from the ER–Golgi interface. These non–membrane-bound microcompartments exclude ribosomes and consist of a core of densely clustered 26S proteasomes surrounded by a loose cloud of Cdc48. Active proteasomes in the microcompartments directly engage with putative substrate at the ER membrane, a function canonically assigned to Cdc48. Live-cell fluorescence microscopy revealed that the proteasome clusters are dynamic, with frequent assembly and fusion events. We propose that the microcompartments perform ER-associated degradation, colocalizing the degradation machinery at specific ER hot spots to enable efficient protein quality control.

STRAIN: an R package for multi-locus sequence typing from whole genome sequencing data

Background

Multi-locus sequence typing (MLST) is a standard typing technique used to associate a sequence type (ST) to a bacterial isolate. When the output of whole genome sequencing (WGS) of a sample is available the ST can be assigned directly processing the read-set.

Current approaches employ reads mapping (SRST2) against the MLST loci, k-mer distribution (stringMLST), selective assembly (GRAbB) or whole genome assembly (BIGSdb) followed by BLASTn sequence query.

Here we present STRAIN (ST Reduced Assembly IdentificatioN), an R package that implements a hybrid strategy between assembly and mapping of the reads to assign the ST to an isolate starting from its read-sets.

Results

Analysis of 540 publicly accessible Illumina read sets showed STRAIN to be more accurate at correct allele assignment and new alleles identification compared to SRTS2, stringMLST and GRAbB. STRAIN assigned correctly 3666 out of 3780 alleles (capability to identify correct alleles 97%) and, when presented with samples containing new alleles, identified them in 3730 out of 3780 STs (capability to identify new alleles 98.7%) of the cases. On the same dataset the other tested tools achieved lower capability to identify correct alleles (from 28.5 to 96.9%) and lower capability to identify new alleles (from 1.1 to 97.1%).

Conclusions

STRAIN is a new accurate method to assign the alleles and ST to an isolate by processing the raw reads output of WGS. STRAIN is also able to retrieve new allele sequences if present. Capability to identify correct and new STs/alleles, evaluated on a benchmark dataset, are higher than other existing methods. STRAIN is designed for single allele typing as well as MLST. Its implementation in R makes allele and ST assignment simple, direct and prompt to be integrated in wider pipeline of downstream bioinformatics analyses.

Tracking outbreak populations of the pepper weevil Anthonomus eugenii (Coleoptera; Curculionidae) using complete mitochondrial genomes

The pepper weevil, Anthonomus eugenii, is a major pest on Capsicum species. Apart from natural spread, there is a risk of spread via international pepper trade. In the Netherlands, a pepper weevil outbreak occurred in 2012 and affected six greenhouses producing different sweet pepper varieties. The following year, a pepper weevil outbreak occurred in Italy. To trace the origin of the Dutch outbreak and to establish if the Dutch and Italian outbreaks were linked, we determined the mitogenomes of A. eugenii specimens collected at outbreak locations, and compared these with specimens from the native area, and other areas where the pest was introduced either by natural dispersal or via trade. The circular 17,257 bp A. eugenii mitogenome comprises thirteen mitochondrial genes typically found in insect species. Intra-species variation of these mitochondrial genes revealed four main mitochondrial lineages encompassing 41 haplotypes. The highest diversity was observed for specimens from its presumed native area (i.e. Mexico). The Dutch outbreak specimens represented three highly similar haplotypes, suggesting a single introduction of the pest. The major Dutch haplotype was also found in two specimens from New Jersey. As the Netherlands does not have pepper trade with New Jersey, it is likely that the specimens sampled in New Jersey and those sampled in the Netherlands originate from a shared source that was not included in this study. In addition, our analysis shows that the Italian and Dutch outbreaks were not linked. The mitochondrial genome is a useful tool to trace outbreak populations and the methodology presented in this paper could prove valuable for other invasive pest species, such as the African fruit moth Thaumatotibia leucotreta and emerald ash borer Agrilus planipennis.

riboSeed: leveraging prokaryotic genomic architecture to assemble across ribosomal regions

The vast majority of bacterial genome sequencing has been performed using Illumina short reads. Because of the inherent difficulty of resolving repeated regions with short reads alone, only ∼10% of sequencing projects have resulted in a closed genome. The most common repeated regions are those coding for ribosomal operons (rDNAs), which occur in a bacterial genome between 1 and 15 times, and are typically used as sequence markers to classify and identify bacteria. Here, we exploit the genomic context in which rDNAs occur across taxa to improve assembly of these regions relative to de novo sequencing by using the conserved nature of rDNAs across taxa and the uniqueness of their flanking regions within a genome. We describe a method to construct targeted pseudocontigs generated by iteratively assembling reads that map to a reference genome’s rDNAs. These pseudocontigs are then used to more accurately assemble the newly sequenced chromosome. We show that this method, implemented as riboSeed, correctly bridges across adjacent contigs in bacterial genome assembly and, when used in conjunction with other genome polishing tools, can assist in closure of a genome.

Assembly of Ancient Mitochondrial Genomes Without a Closely Related Reference Sequence

Recent methodological advances have transformed the field of ancient DNA (aDNA). Basic bioinformatics skills are becoming essential requirements to process and analyze the sheer amounts of data generated by current aDNA studies and in biomedical research in general. This chapter is intended as a practical guide to the assembly of ancient mitochondrial genomes, directly from genomic DNA-derived next-generation sequencing (NGS) data, specifically in the absence of closely related reference genomes. In a hands-on tutorial suitable for readers with little to no prior bioinformatics experience, we reconstruct the mitochondrial genome of a woolly mammoth deposited ~45,000 years ago. We introduce key software tools and outline general strategies for mitogenome assembly, including the critical quality assessment of assembly results without a reference genome.

First steps towards mitochondrial pan-genomics: detailed analysis of Fusarium graminearum mitogenomes

There is a gradual shift from representing a species' genome by a single reference genome sequence to a pan-genome representation. Pan-genomes are the abstract representations of the genomes of all the strains that are present in the population or species. In this study, we employed a pan-genomic approach to analyze the intraspecific mitochondrial genome diversity of Fusarium graminearum. We present an improved reference mitochondrial genome for F. graminearum with an intron-exon annotation that was verified using RNA-seq data. Each of the 24 studied isolates had a distinct mitochondrial sequence. Length variation in the F. graminearum mitogenome was found to be largely due to variation of intron regions (99.98%). The "intronless" mitogenome length was found to be quite stable and could be informative when comparing species. The coding regions showed high conservation, while the variability of intergenic regions was highest. However, the most important variable parts are the intron regions, because they contain approximately half of the variable sites, make up more than half of the mitogenome, and show presence/absence variation. Furthermore, our analyses show that the mitogenome of F. graminearum is recombining, as was previously shown in F. oxysporum, indicating that mitogenome recombination is a common phenomenon in Fusarium. The majority of mitochondrial introns in F. graminearum belongs to group I introns, which are associated with homing endonuclease genes (HEGs). Mitochondrial introns containing HE genes may spread within populations through homing, where the endonuclease recognizes and cleaves the recognition site in the target gene. After cleavage of the "host" gene, it is replaced by the gene copy containing the intron with HEG. We propose to use introns unique to a population for tracking the spread of the given population, because introns can spread through vertical inheritance, recombination as well as via horizontal transfer. We demonstrate how pooled sequencing of strains can be used for mining mitogenome data. The usage of pooled sequencing offers a scalable solution for population analysis and for species level comparisons studies. This study may serve as a basis for future mitochondrial genome variability studies and representations.

The linear mitochondrial genome of the quarantine chytrid Synchytrium endobioticum; insights into the evolution and recent history of an obligate biotrophic plant pathogen

BACKGROUND

Chytridiomycota species (chytrids) belong to a basal lineage in the fungal kingdom. Inhabiting terrestrial and aquatic environments, most are free-living saprophytes but several species cause important diseases: e.g. Batrachochytrium dendrobatidis, responsible for worldwide amphibian decline; and Synchytrium endobioticum, causing potato wart disease. S. endobioticum has an obligate biotrophic lifestyle and isolates can be further characterized as pathotypes based on their virulence on a differential set of potato cultivars. Quarantine measures have been implemented globally to control the disease and prevent its spread. We used a comparative approach using chytrid mitogenomes to determine taxonomical relationships and to gain insights into the evolution and recent history of introductions of this plant pathogen.

RESULTS

We assembled and annotated the complete mitochondrial genome of 30 S. endobioticum isolates and generated mitochondrial genomes for five additional chytrid species. The mitochondrial genome of S. endobioticum is linear with terminal inverted repeats which was validated by tailing and PCR amplifying the telomeric ends. Surprisingly, no conservation in organisation and orientation of mitochondrial genes was observed among the Chytridiomycota except for S. endobioticum and its sister species Synchytrium microbalum. However, the mitochondrial genome of S. microbalum is circular and comprises only a third of the 72.9 Kbp found for S. endobioticum suggesting recent linearization and expansion. Four mitochondrial lineages were identified in the S. endobioticum mitochondrial genomes. Several pathotypes occur in different lineages, suggesting that these have emerged independently. In addition, variations for polymorphic sites in the mitochondrial genome of individual isolates were observed demonstrating that S. endobioticum isolates represent a community of different genotypes. Such communities were shown to be complex and stable over time, but we also demonstrate that the use of semi-resistant potato cultivars triggers a rapid shift in the mitochondrial haplotype associated with increased virulence.

CONCLUSIONS

Mitochondrial genomic variation shows that S. endobioticum has been introduced into Europe multiple times, that several pathotypes emerged multiple times, and that isolates represent communities of different genotypes. Our study represents the most comprehensive dataset of chytrid mitogenomes, which provides new insights into the extraordinary dynamics and evolution of mitochondrial genomes involving linearization, expansion and reshuffling.

Kollector: transcript-informed, targeted de novo assembly of gene loci

Motivation

Despite considerable advancements in sequencing and computing technologies, de novo assembly of whole eukaryotic genomes is still a time-consuming task that requires a significant amount of computational resources and expertise. A targeted assembly approach to perform local assembly of sequences of interest remains a valuable option for some applications. This is especially true for gene-centric assemblies, whose resulting sequence can be readily utilized for more focused biological research. Here we describe Kollector, an alignment-free targeted assembly pipeline that uses thousands of transcript sequences concurrently to inform the localized assembly of corresponding gene loci. Kollector robustly reconstructs introns and novel sequences within these loci, and scales well to large genomes—properties that makes it especially useful for researchers working on non-model eukaryotic organisms.

Results

We demonstrate the performance of Kollector for assembling complete or near-complete Caenorhabditis elegans and Homo sapiens gene loci from their respective, input transcripts. In a time- and memory-efficient manner, the Kollector pipeline successfully reconstructs respectively 99% and 80% (compared to 86% and 73% with standard de novo assembly techniques) of C.elegans and H.sapiens transcript targets in their corresponding genomic space using whole genome shotgun sequencing reads. We also show that Kollector outperforms both established and recently released targeted assembly tools. Finally, we demonstrate three use cases for Kollector, including comparative and cancer genomics applications.

Availability and Implementation

Kollector is implemented as a bash script, and is available at https://github.com/bcgsc/kollector

Supplementary information

Supplementary data are available at Bioinformatics online.

Genomic Understanding of an Infectious Brain Disease from the Desert

Rhinocladiella mackenziei accounts for the majority of fungal brain infections in the Middle East, and is restricted to the arid climate zone between Saudi Arabia and Pakistan. Neurotropic dissemination caused by this fungus has been reported in immunocompromised, but also immunocompetent individuals. If untreated, the infection is fatal. Outside of humans, the environmental niche of R. mackenziei is unknown, and the fungus has been only cultured from brain biopsies. In this paper, we describe the whole-genome resequencing of two R. mackenziei strains from patients in Saudi Arabia and Qatar. We assessed intraspecies variation and genetic signatures to uncover the genomic basis of the pathogenesis, and potential niche adaptations. We found that the duplicated genes (paralogs) are more susceptible to accumulating significant mutations. Comparative genomics with other filamentous ascomycetes revealed a diverse arsenal of genes likely engaged in pathogenicity, such as the degradation of aromatic compounds and iron acquisition. In addition, intracellular accumulation of trehalose and choline suggests possible adaptations to the conditions of an arid climate region. Specifically, protein family contractions were found, including short-chain dehydrogenase/reductase SDR, the cytochrome P450 (CYP) (E-class), and the G-protein β WD-40 repeat. Gene composition and metabolic potential indicate extremotolerance and hydrocarbon assimilation, suggesting a possible environmental habitat of oil-polluted desert soil.

Phylogeny of dermatophytes with genomic character evaluation of clinically distinct Trichophyton rubrum and T. violaceum

Trichophyton rubrum and T. violaceum are prevalent agents of human dermatophyte infections, the former being found on glabrous skin and nail, while the latter is confined to the scalp. The two species are phenotypically different but are highly similar phylogenetically. The taxonomy of dermatophytes is currently being reconsidered on the basis of molecular phylogeny. Molecular species definitions do not always coincide with existing concepts which are guided by ecological and clinical principles. In this article, we aim to bring phylogenetic and ecological data together in an attempt to develop new species concepts for anthropophilic dermatophytes. Focus is on the T. rubrum complex with analysis of rDNA ITS supplemented with LSU, TUB2, TEF3 and ribosomal protein L10 gene sequences. In order to explore genomic differences between T. rubrum and T. violaceum, one representative for both species was whole genome sequenced. Draft sequences were compared with currently available dermatophyte genomes. Potential virulence factors of adhesins and secreted proteases were predicted and compared phylogenetically. General phylogeny showed clear gaps between geophilic species of Arthroderma, but multilocus distances between species were often very small in the derived anthropophilic and zoophilic genus Trichophyton. Significant genome conservation between T. rubrum and T. violaceum was observed, with a high similarity at the nucleic acid level of 99.38 % identity. Trichophyton violaceum contains more paralogs than T. rubrum. About 30 adhesion genes were predicted among dermatophytes. Seventeen adhesins were common between T. rubrum and T. violaceum, while four were specific for the former and eight for the latter. Phylogenetic analysis of secreted proteases reveals considerable expansion and conservation among the analyzed species. Multilocus phylogeny and genome comparison of T. rubrum and T. violaceum underlined their close affinity. The possibility that they represent a single species exhibiting different phenotypes due to different localizations on the human body is discussed.

Epichloë hybrida, sp. nov., an emerging model system for investigating fungal allopolyploidy

Endophytes of the genus Epichloë (Clavicipitaceae, Ascomycota) frequently occur within cool-season grasses and form interactions with their hosts that range from mutualistic to antagonistic. Many Epichloë species have arisen via interspecific hybridization, resulting in species with two or three subgenomes that retain all or nearly all of their original parental genomes, a process termed allopolyploidization. Here, we characterize Epichloë hybrida, sp. nov., a mutualistic species that has increasingly become a model system for investigating allopolyploidy in fungi. The Epichloë species so far identified as the closest known relatives of the two progenitors of E. hybrida are E. festucae var. lolii and E. typhina. We confirm that the nuclear genome of E. hybrida contains two homeologs of most protein-coding genes from E. festucae and E. typhina, with genome-wide gene expression analysis indicating a slight bias in overall gene expression from the E. typhina subgenome. Mitochondrial DNA is detectable only from E. festucae, whereas ribosomal DNA is detectable only from E. typhina. Inheriting ribosomal DNA from just one parent might be expected to preferentially favor interactions with ribosomal proteins from the same parent, but we find that ribosomal protein genes from both parental subgenomes are nearly all expressed equally in E. hybrida. Finally, we provide a comprehensive set of resources for this model system that are intended to facilitate further study of fungal hybridization by other researchers.

Evolution within the fungal genus Verticillium is characterized by chromosomal rearrangement and gene loss

The fungal genus Verticillium contains ten species, some of which are notorious plant pathogens causing vascular wilt diseases in host plants, while others are known as saprophytes and opportunistic plant pathogens. Whereas the genome of V. dahliae, the most notorious plant pathogen of the genus, has been well characterized, evolution and speciation of other members of the genus received little attention thus far. Here, we sequenced the genomes of the nine haploid Verticillium spp. to study evolutionary trajectories of their divergence from a last common ancestor. Frequent occurrence of chromosomal rearrangement and gene family loss was identified. In addition to ∼11 000 genes that are shared at least between two species, only 200-600 species-specific genes occur. Intriguingly, these species-specific genes show different features than the shared genes.

Complete mitogenome of the high ethanol production fungus Fusarium oxysporum Mh2-2

Fusarium spp. are significantly important plant pathogens, and some of them are ethanol-producing strains. During infection and/or ethanol production, Fusarium requires a plenty of energy that is mainly provided by mitochondria. Here we report the first mitogenome from a selected Fusarium oxysporum strain mh2-2 that produces ethanol from glucose and xylose. The size of this mitogenome, 46 kb, is different from the size of any reported Fusarium mitogenome. Our results provide insight into the functions and evolution of mitochondrial genes and genomes.

van der Lee TA, Waalwijk C, van Diepeningen AD. Mitochondrial genomes reveal recombination in the presumed asexual Fusarium oxysporum species complex

BACKGROUND

The Fusarium oxysporum species complex (FOSC) contains several phylogenetic lineages. Phylogenetic studies identified two to three major clades within the FOSC. The mitochondrial sequences are highly informative phylogenetic markers, but have been mostly neglected due to technical difficulties.

RESULTS

A total of 61 complete mitogenomes of FOSC strains were de novo assembled and annotated. Length variations and intron patterns support the separation of three phylogenetic species. The variable region of the mitogenome that is typical for the genus Fusarium shows two new variants in the FOSC. The variant typical for Fusarium is found in members of all three clades, while variant 2 is found in clades 2 and 3 and variant 3 only in clade 2. The extended set of loci analyzed using a new implementation of the genealogical concordance species recognition method support the identification of three phylogenetic species within the FOSC. Comparative analysis of the mitogenomes in the FOSC revealed ongoing mitochondrial recombination within, but not between phylogenetic species.

CONCLUSIONS

The recombination indicates the presence of a parasexual cycle in F. oxysporum. The obstacles hindering the usage of the mitogenomes are resolved by using next generation sequencing and selective genome assemblers, such as GRAbB. Complete mitogenome sequences offer a stable basis and reference point for phylogenetic and population genetic studies.

A mobile pathogenicity chromosome in Fusarium oxysporum for infection of multiple cucurbit species

The genome of Fusarium oxysporum (Fo) consists of a set of eleven 'core' chromosomes, shared by most strains and responsible for housekeeping, and one or several accessory chromosomes. We sequenced a strain of Fo f.sp. radicis-cucumerinum (Forc) using PacBio SMRT sequencing. All but one of the core chromosomes were assembled into single contigs, and a chromosome that shows all the hallmarks of a pathogenicity chromosome comprised two contigs. A central part of this chromosome contains all identified candidate effector genes, including homologs of SIX6, SIX9, SIX11 and SIX 13. We show that SIX6 contributes to virulence of Forc. Through horizontal chromosome transfer (HCT) to a non-pathogenic strain, we also show that the accessory chromosome containing the SIX gene homologs is indeed a pathogenicity chromosome for cucurbit infection. Conversely, complete loss of virulence was observed in Forc016 strains that lost this chromosome. We conclude that also a non-wilt-inducing Fo pathogen relies on effector proteins for successful infection and that the Forc pathogenicity chromosome contains all the information necessary for causing root rot of cucurbits. Three out of nine HCT strains investigated have undergone large-scale chromosome alterations, reflecting the remarkable plasticity of Fo genomes.

Whole-Genome Sequencing and In Silico Analysis of Two Strains of Sporothrix globosa

Sporothrix globosa is a thermo-dimorphic fungus belonging to a pathogenic clade that also includes Sporothrix schenckii, which causes human and animal sporotrichosis. Here, we present the first genome assemblies of two S. globosa strains providing data for future comparative genomic studies in pathogenic Sporothrix species.

Draft Genome Sequence of the Dimorphic Fungus Sporothrix pallida, a Nonpathogenic Species Belonging to Sporothrix, a Genus Containing Agents of Human and Feline Sporotrichosis

Sporothrix pallida is considered to be a mostly avirulent environmental fungus, phylogenetically closely related to the well-known pathogen Sporothrix schenckii. Here, we present the first assembly of its genome, which provides a valuable resource for future comparative genomic studies between nonpathogenic and pathogenic Sporothrix spp.