Single-worm long-read sequencing reveals genome diversity in free-living nematodes.

Obtaining sufficient genetic material from a limited biological source is currently the primary operational bottleneck in studies investigating biodiversity and genome evolution. In this study, we employed multiple displacement amplification (MDA) and Smartseq2 to amplify nanograms of genomic DNA and mRNA, respectively, from individual Caenorhabditis elegans. Although reduced genome coverage was observed in repetitive regions, we produced assemblies covering 98% of the reference genome using long-read sequences generated with Oxford Nanopore Technologies (ONT). Annotation with the sequenced transcriptome coupled with the available assembly revealed that gene predictions were more accurate, complete and contained far fewer false positives than de novo transcriptome assembly approaches. We sampled and sequenced the genomes and transcriptomes of 13 nematodes from early-branching species in Chromadoria, Dorylaimia and Enoplia. The basal Chromadoria and Enoplia species had larger genome sizes, ranging from 136.6 to 738.8 Mb, compared with those in the other clades. Nine mitogenomes were fully assembled, and displayed a complete lack of synteny to other species. Phylogenomic analyses based on the new annotations revealed strong support for Enoplia as sister to the rest of Nematoda. Our result demonstrates the robustness of MDA in combination with ONT, paving the way for the study of genome diversity in the phylum Nematoda and beyond.

The Aphelenchoides genomes reveal substantial horizontal gene transfers in the last common ancestor of free-living and major plant-parasitic nematodes.

Aphelenchoides besseyi is a plant-parasitic nematode (PPN) in the family Aphelenchoididae capable of infecting more than 200 plant species. A. besseyi is also a species complex with strains exhibiting varying pathogenicity to plants. We present the genome and annotations of six Aphelenchoides species, four of which belonged to the A. besseyi species complex. Most Aphelenchoides genomes have a size of 44.7-47.4 Mb and are among the smallest in clade IV, with the exception of A. fujianensis, which has a size of 143.8 Mb and is one of the largest. Phylogenomic analysis successfully delimited the species complex into A. oryzae and A. pseudobesseyi and revealed a reduction of transposon elements in the last common ancestor of Aphelenchoides. Synteny analyses between reference genomes indicated that three chromosomes in A. besseyi were derived from fission and fusion events. A systematic identification of horizontal gene transfer (HGT) genes across 27 representative nematodes allowed us to identify two major episodes of acquisition corresponding to the last common ancestor of clade IV or major PPNs, respectively. These genes were mostly lost and differentially retained between clades or strains. Most HGT events were acquired from bacteria, followed by fungi, and also from plants; plant HGT was especially prevalent in Bursaphelenchus mucronatus. Our results comprehensively improve the understanding of HGT in nematodes.

Translocation of Fungicides and Their Efficacy in Controlling Phellinus noxius, the Cause of Brown Root Rot Disease.

Brown root rot disease (BRRD), caused by Phellinus noxius, is an important tree disease in tropical and subtropical areas. To improve chemical control of BRRD and deter emergence of fungicide resistance in P. noxius, this study investigated control efficacies and systemic activities of fungicides with different modes of action. Fourteen fungicides with 11 different modes of action were tested for inhibitory effects in vitro on 39 P. noxius isolates from Taiwan, Hong Kong, Malaysia, Australia, and Pacific Islands. Cyproconazole, epoxiconazole, and tebuconazole (Fungicide Resistance Action Committee [FRAC] 3, target-site G1) inhibited colony growth of P. noxius by 99.9 to 100% at 10 ppm and 97.7 to 99.8% at 1 ppm. The other effective fungicide was cyprodinil + fludioxonil (FRAC 9 + 12, target-site D1 + E2), which showed growth inhibition of 96.9% at 10 ppm and 88.6% at 1 ppm. Acropetal translocation of six selected fungicides was evaluated in bishop wood (Bischofia javanica) seedlings by immersion of the root tips in each fungicide at 100 ppm, followed by liquid or gas chromatography tandem mass spectrometry analyses of consecutive segments of root, stem, and leaf tissues at 7 and 21 days posttreatment. Bidirectional translocation of the fungicides was also evaluated by stem injection of fungicide stock solutions. Cyproconazole and tebuconazole were the most readily absorbed by roots and efficiently transported acropetally. Greenhouse experiments suggested that cyproconazole, tebuconazole, and epoxiconazole have a slightly higher potential for controlling BRRD than mepronil, prochloraz, and cyprodinil + fludioxonil. Because all tested fungicides lacked basipetal translocation, soil drenching should be considered instead of trunk injection for their use in BRRD control.

Comparative genomic and transcriptomic analyses of trans-kingdom pathogen Fusarium solani species complex reveal degrees of compartmentalization.

BACKGROUND: The Fusarium solani species complex (FSSC) comprises fungal pathogens responsible for mortality in a diverse range of animals and plants, but their genome diversity and transcriptome responses in animal pathogenicity remain to be elucidated. We sequenced, assembled and annotated six chromosome-level FSSC clade 3 genomes of aquatic animal and plant host origins. We established a pathosystem and investigated the expression data of F. falciforme and F. keratoplasticum in Chinese softshell turtle (Pelodiscus sinensis) host. RESULTS: Comparative analyses between the FSSC genomes revealed a spectrum of conservation patterns in chromosomes categorised into three compartments: core, fast-core (FC), and lineage-specific (LS). LS chromosomes contribute to variations in genomes size, with up to 42.2% of variations between F. vanettenii strains. Each chromosome compartment varied in structural architectures, with FC and LS chromosomes contain higher proportions of repetitive elements with genes enriched in functions related to pathogenicity and niche expansion. We identified differences in both selection in the coding sequences and DNA methylation levels between genome features and chromosome compartments which suggest a multi-speed evolution that can be traced back to the last common ancestor of Fusarium. We further demonstrated that F. falciforme and F. keratoplasticum are opportunistic pathogens by inoculating P. sinensis eggs and identified differentially expressed genes also associated with plant pathogenicity. These included the most upregulated genes encoding the CFEM (Common in Fungal Extracellular Membrane) domain. CONCLUSIONS: The high-quality genome assemblies provided new insights into the evolution of FSSC chromosomes, which also serve as a resource for studies of fungal genome evolution and pathogenesis. This study also establishes an animal model for fungal pathogens of trans-kingdom hosts.

Utilisation of Oxford Nanopore sequencing to generate six complete gastropod mitochondrial genomes as part of a biodiversity curriculum.

High-throughput sequencing has enabled genome skimming approaches to produce complete mitochondrial genomes (mitogenomes) for species identification and phylogenomics purposes. In particular, the portable sequencing device from Oxford Nanopore Technologies (ONT) has the potential to facilitate hands-on training from sampling to sequencing and interpretation of mitogenomes. In this study, we present the results from sampling and sequencing of six gastropod mitogenomes (Aplysia argus, Cellana orientalis, Cellana toreuma, Conus ebraeus, Conus miles and Tylothais aculeata) from a graduate level biodiversity course. The students were able to produce mitogenomes from sampling to annotation using existing protocols and programs. Approximately 4 Gb of sequence was produced from 16 Flongle and one MinION flow cells, averaging 235 Mb and N50 = 4.4 kb per flow cell. Five of the six 14.1-18 kb mitogenomes were circlised containing all 13 core protein coding genes. Additional Illumina sequencing revealed that the ONT assemblies spanned over highly AT rich sequences in the control region that were otherwise missing in Illumina-assembled mitogenomes, but still contained a base error of one every 70.8-346.7 bp under the fast mode basecalling with the majority occurring at homopolymer regions. Our findings suggest that the portable MinION device can be used to rapidly produce low-cost mitogenomes onsite and tailored to genomics-based training in biodiversity research.

Genome Assembly and Evolutionary Analysis of the Mandarin Duck Aix galericulata Reveal Strong Genome Conservation among Ducks.

The mandarin duck, Aix galericulata, is popular in East Asian cultures and displays exaggerated sexual dimorphism, especially in feather traits during breeding seasons. We generated and annotated the first mandarin duck de novo assembly, which was 1.08 Gb in size and encoded 16,615 proteins. Using a phylogenomic approach calibrated with fossils and molecular divergences, we inferred that the last common ancestor of ducks occurred 13.3-26.7 Ma. The majority of the mandarin duck genome repetitive sequences belonged to the chicken repeat 1 (CR1) retroposon CR1-J2_Pass, which underwent a duck lineage-specific burst. Synteny analyses among ducks revealed infrequent chromosomal rearrangements in which breaks were enriched in LINE retrotransposons and DNA transposons. The calculation of the dN/dS ratio revealed that the majority of duck genes were under strong purifying selection. The expanded gene families in the mandarin duck are primarily involved in olfactory perception as well as the development and morphogenesis of feather and branching structures. This new reference genome will improve our understanding of the morphological and physiological characteristics of ducks and provide a valuable resource for functional genomics studies to investigate the feather traits of the mandarin duck.

Extensive sampling of Saccharomyces cerevisiae in Taiwan reveals ecology and evolution of predomesticated lineages.

The ecology and genetic diversity of the model yeast Saccharomyces cerevisiae before human domestication remain poorly understood. Taiwan is regarded as part of this yeast's geographic birthplace, where the most divergent natural lineage was discovered. Here, we extensively sampled the broadleaf forests across this continental island to probe the ancestral species' diversity. We found that S. cerevisiae is distributed ubiquitously at low abundance in the forests. Whole-genome sequencing of 121 isolates revealed nine distinct lineages that diverged from Asian lineages during the Pleistocene, when a transient continental shelf land bridge connected Taiwan to other major landmasses. Three lineages are endemic to Taiwan and six are widespread in Asia, making this region a focal biodiversity hotspot. Both ancient and recent admixture events were detected between the natural lineages, and a genetic ancestry component associated with isolates from fruits was detected in most admixed isolates. Collectively, Taiwanese isolates harbor genetic diversity comparable to that of the whole Asia continent, and different lineages have coexisted at a fine spatial scale even on the same tree. Patterns of variations within each lineage revealed that S. cerevisiae is highly clonal and predominantly reproduces asexually in nature. We identified different selection patterns shaping the coding sequences of natural lineages and found fewer gene family expansion and contractions that contrast with domesticated lineages. This study establishes that S. cerevisiae has rich natural diversity sheltered from human influences, making it a powerful model system in microbial ecology.

Mycena genomes resolve the evolution of fungal bioluminescence.

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species (Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ∼150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

Evidence of Extensive Intraspecific Noncoding Reshuffling in a 169-kb Mitochondrial Genome of a Basidiomycetous Fungus.

Comparative genomics of fungal mitochondrial genomes (mitogenomes) have revealed a remarkable pattern of rearrangement between and within major phyla owing to horizontal gene transfer and recombination. The role of recombination was exemplified at a finer evolutionary time scale in basidiomycetes group of fungi as they display a diversity of mitochondrial DNA inheritance patterns. Here, we assembled mitogenomes of six species from the Hymenochaetales order of basidiomycetes and examined 59 mitogenomes from 2 genetic lineages of Phellinus noxius. Gene order is largely collinear, while intergene regions are major determinants of mitogenome size variation. Substantial sequence divergence was found in shared introns consistent with high horizontal gene transfer frequency observed in yeasts, but we also identified a rare case where an intron was retained in five species since speciation. In contrast to the hyperdiversity observed in nuclear genomes of Phellinus noxius, mitogenomes' intraspecific polymorphisms at protein-coding sequences are extremely low. Phylogeny network based on introns revealed turnover as well as exchange of introns between two lineages. Strikingly, some strains harbor a mosaic origin of introns from both lineages. Analysis of intergenic sequence indicated substantial differences between and within lineages, and an expansion may be ongoing as a result of exchange between distal intergenes. These findings suggest that the evolution in mitochondrial DNAs is usually lineage specific but chimeric mitotypes are frequently observed, thus capturing the possible evolutionary processes shaping mitogenomes in a basidiomycete. The large mitogenome sizes reported in various basidiomycetes appear to be a result of interspecific reshuffling of intergenes.

Comparative and population genomic landscape of Phellinus noxius: A hypervariable fungus causing root rot in trees.

The order Hymenochaetales of white rot fungi contain some of the most aggressive wood decayers causing tree deaths around the world. Despite their ecological importance and the impact of diseases they cause, little is known about the evolution and transmission patterns of these pathogens. Here, we sequenced and undertook comparative genomic analyses of Hymenochaetales genomes using brown root rot fungus Phellinus noxius, wood-decomposing fungus Phellinus lamaensis, laminated root rot fungus Phellinus sulphurascens and trunk pathogen Porodaedalea pini. Many gene families of lignin-degrading enzymes were identified from these fungi, reflecting their ability as white rot fungi. Comparing against distant fungi highlighted the expansion of 1,3-beta-glucan synthases in P. noxius, which may account for its fast-growing attribute. We identified 13 linkage groups conserved within Agaricomycetes, suggesting the evolution of stable karyotypes. We determined that P. noxius has a bipolar heterothallic mating system, with unusual highly expanded ~60 kb A locus as a result of accumulating gene transposition. We investigated the population genomics of 60 P. noxius isolates across multiple islands of the Asia Pacific region. Whole-genome sequencing showed this multinucleate species contains abundant poly-allelic single nucleotide polymorphisms with atypical allele frequencies. Different patterns of intra-isolate polymorphism reflect mono-/heterokaryotic states which are both prevalent in nature. We have shown two genetically separated lineages with one spanning across many islands despite the geographical barriers. Both populations possess extraordinary genetic diversity and show contrasting evolutionary scenarios. These results provide a framework to further investigate the genetic basis underlying the fitness and virulence of white rot fungi.

The Genetic Structure of Phellinus noxius and Dissemination Pattern of Brown Root Rot Disease in Taiwan.

Since the 1990s, brown root rot caused by Phellinus noxius (Corner) Cunningham has become a major tree disease in Taiwan. This fungal pathogen can infect more than 200 hardwood and softwood tree species, causing gradual to fast decline of the trees. For effective control, we must determine how the pathogen is disseminated and how the new infection center of brown root rot is established. We performed Illumina sequencing and de novo assembly of a single basidiospore isolate Daxi42 and obtained a draft genome of ~40 Mb. By comparing the 12,217 simple sequence repeat (SSR) regions in Daxi42 with the low-coverage Illumina sequencing data for four additional P. noxius isolates, we identified 154 SSR regions with potential polymorphisms. A set of 13 polymorphic SSR markers were then developed and used to analyze 329 P. noxius isolates collected from 73 tree species from urban/agricultural areas in 14 cities/counties all around Taiwan from 1989 to 2012. The results revealed a high proportion (~98%) of distinct multilocus genotypes (MLGs) and that none of the 329 isolates were genome-wide homozygous, which supports a possible predominant outcrossing reproductive mode in P. noxius. The diverse MLGs exist as discrete patches, so brown root rot was most likely caused by multiple clones rather than a single predominant strain. The isolates collected from diseased trees near each other tend to have similar genotype(s), which indicates that P. noxius may spread to adjacent trees via root-to-root contact. Analyses based on Bayesian clustering, FST statistics, analysis of molecular variance, and isolation by distance all suggest a low degree of population differentiation and little to no barrier to gene flow throughout the P. noxius population in Taiwan. We discuss the involvement of basidiospore dispersal in disease dissemination.

Bond cutting in K-doped tris(8-hydroxyquinoline) aluminium.

A series of Al 2p, K 2p, O 1s and N 1s core-level spectra have been used to characterize the interaction between potassium (K) and tris(8-hydroxyquinoline) aluminium (Alq(3)) molecules in the K-doped Alq(3) layer. All core-level spectra were tuned to be very surface sensitive in selecting various photon energies provided by the wide-range beamline at the National Synchrotron Radiation Research Center, Taiwan. A critical K concentration (x = 2.4) exists in the K-doped Alq(3) layer, below which the K-doped atoms generate a strained environment near the O and N atoms within 8-quinolinoline ligands. This creates new O 1s and N 1s components on the lower binding-energy side. Above the critical K coverage, the K-doped atoms attach the O atoms in the Al-O-C bonds next to the phenoxide ring and replace Al-O-C bonds by forming K-O-C bonds. An Alq(3) molecule is disassembled into Alq(2) and Kq by bond cutting and bond formation. The Alq(2) molecule can be further dissociated into Alq, or even Al, through subsequent formations of Kq.

Bond cutting in Cs-doped tris(8-hydroxyquinoline) aluminium.

A series of Cs 4d and Al 2p spectra associated with valence-band and cut-off spectra have been used to characterize the interaction between caesium and tris(8-hydroxyquinoline) aluminium (Alq(3)) molecules in a Cs-doped Alq(3) layer. The Cs 4d and Al 2p spectra were tuned to be very surface sensitive by selecting a photon energy of 120 eV at the National Synchrotron Radiation Research Center, Taiwan. A critical Cs concentration exists, above which a new Al 2p signal appears next to the Al 2p peak of Alq(3) in the lower binding-energy side. The Al 2p signal was analyzed and assigned as being contributed from a mixture of Alq(2), Alq and Al. Experimental data supported the observation that bond cutting of Alq(3) by the doped Cs atoms occurred at high Cs doping concentration.