Chromosome-level assembly of the Phytophthora agathidicida genome reveals adaptation in effector gene families.

Phytophthora species are notorious plant pathogens, with some causing devastating tree diseases that threaten the survival of their host species. One such example is Phytophthora agathidicida, the causal agent of kauri dieback - a root and trunk rot disease that kills the ancient, iconic and culturally significant tree species, Agathis australis (New Zealand kauri). A deeper understanding of how Phytophthora pathogens infect their hosts and cause disease is critical for the development of effective treatments. Such an understanding can be gained by interrogating pathogen genomes for effector genes, which are involved in virulence or pathogenicity. Although genome sequencing has become more affordable, the complete assembly of Phytophthora genomes has been problematic, particularly for those with a high abundance of repetitive sequences. Therefore, effector genes located in repetitive regions could be truncated or missed in a fragmented genome assembly. Using a combination of long-read PacBio sequences, chromatin conformation capture (Hi-C) and Illumina short reads, we assembled the P. agathidicida genome into ten complete chromosomes, with a genome size of 57 Mb including 34% repeats. This is the first Phytophthora genome assembled to chromosome level and it reveals a high level of syntenic conservation with the complete genome of Peronospora effusa, the only other completely assembled genome sequence of an oomycete. All P. agathidicida chromosomes have clearly defined centromeres and contain candidate effector genes such as RXLRs and CRNs, but in different proportions, reflecting the presence of gene family clusters. Candidate effector genes are predominantly found in gene-poor, repeat-rich regions of the genome, and in some cases showed a high degree of duplication. Analysis of candidate RXLR effector genes that occur in multicopy gene families indicated half of them were not expressed in planta. Candidate CRN effector gene families showed evidence of transposon-mediated recombination leading to new combinations of protein domains, both within and between chromosomes. Further analysis of this complete genome assembly will help inform new methods of disease control against P. agathidicida and other Phytophthora species, ultimately helping decipher how Phytophthora pathogens have evolved to shape their effector repertoires and how they might adapt in the future.

A single fungal strain was the unexpected cause of a mass aspergillosis outbreak in the world's largest and only flightless parrot.

Kakapo are a critically endangered species of parrots restricted to a few islands off the coast of New Zealand. Kakapo are very closely monitored, especially during nesting seasons. In 2019, during a highly successful nesting season, an outbreak of aspergillosis affected 21 individuals and led to the deaths of 9, leaving a population of only 211 kakapo. In monitoring this outbreak, cultures of aspergillus were grown, and genome sequenced. These sequences demonstrate that, very unusually for an aspergillus outbreak, a single strain of aspergillus caused the outbreak. This strain was found on two islands, but only one had an outbreak of aspergillosis; indicating that the strain was necessary, but not sufficient, to cause disease. Our analysis provides an understanding of the 2019 outbreak and provides potential ways to manage such events in the future.

Cross-kingdom transcriptomic trends in the evolution of hybrid gene expression.

Hybridization is a route to speciation that occurs widely across the eukaryote tree of life. The success of allopolyploids (hybrid species with increased ploidy) and homoploid hybrids (with unchanged ploidy) is well documented. However, their formation and establishment is not straightforward, with a suite of near-instantaneous and longer term biological repercussions faced by the new species. Central to these challenges is the rewiring of gene regulatory networks following the merger of distinct genomes inherited from both parental species. Research on the evolution of hybrid gene expression has largely involved studies on a single hybrid species or a few gene families. Here, we present the first standardized transcriptome-wide study exploring the fates of genes following hybridization across three kingdoms: animals, plants and fungi. Within each kingdom, we pair an allopolyploid system with a closely related homoploid hybrid to decouple the influence of increased ploidy from genome merger. Genome merger, not changes in ploidy, has the greatest effect on posthybridization expression patterns across all study systems. Strikingly, we find that differentially expressed genes in parent species preferentially switch to more similar expression in hybrids across all kingdoms, likely as a consequence of regulatory trans-acting cross-talk within the hybrid nucleus. We also highlight the prevalence of gene loss or silencing among extremely differentially expressed genes in hybrid species across all kingdoms. These shared patterns suggest that the evolutionary process of hybridization leads to common high-level expression outcomes, regardless of the particular species or kingdom.

Chromosome-level genomes provide insights into genome evolution, organization and size in Epichloe fungi.

Epichloe fungi are endophytes of cool season grasses, both wild species and commercial cultivars, where they may exhibit mutualistic or pathogenic lifestyles. The Epichloe-grass symbiosis is of great interest to agricultural research for the fungal bioprotective properties conferred to host grasses but also serves as an ideal system to study the evolution of fungal plant-pathogens in natural environments. Here, we assembled and annotated gapless chromosome-level genomes of two pathogenic Epichloe sibling species. Both genomes have a bipartite genome organization, with blocks of highly syntenic gene-rich regions separated by blocks of AT-rich DNA. The AT-rich regions show an extensive signature of RIP (repeat-induced point mutation) and the expansion of this compartment accounts for the large difference in genome size between the two species. This study reveals how the rapid evolution of repeat structure can drive divergence between closely related taxa and highlights the evolutionary role of dynamic compartments in fungal genomes.

Keratin intermediate filament chains in the European common wall lizard (Podarcis muralis) and a potential keratin filament crosslinker.

On the basis of sequence homology with mammalian α-keratins, and on the criteria that the coiled-coil segments and central linker in the rod domain of these molecules must have conserved lengths if they are to assemble into viable intermediate filaments, a total of 28 Type I and Type II keratin intermediate filament chains (KIF) have been identified from the genome of the European common wall lizard (Podarcis muralis). Using the same criteria this number may be compared to 33 found here in the green anole lizard (Anole carolinensis) and 25 in the tuatara (Sphenodon punctatus). The Type I and Type II KIF genes in the wall lizard fall in clusters on chromosomes 13 and 2 respectively. Although some differences occur in the terminal domains in the KIF chains of the two lizards and tuatara, the similarities between key indicator residues - cysteine, glycine and proline - are significant. The terminal domains of the KIF chains in the wall lizard also contain sequence repeats commonly based on glycine and large apolar residues and would permit the fine tuning of physical properties when incorporated within the intermediate filaments. The H1 domain in the Type II chain is conserved across the lizards, tuatara and mammals, and has been related to its role in assembly at the 2-4 molecule level. A KIF-like chain (K80) with an extensive tail domain comprised of multiple tandem repeats has been identified as having a potential filament-crosslinking role.

Keratin intermediate filament chains in tuatara (Sphenodon punctatus): A comparison of tuatara and human sequences.

Determination of the sequences of the keratin intermediate filament chains in tuatara has shown that these are closely akin to the α-keratins in human and other vertebrates, especially in the central, coiled-coil rod region. The domain lengths within the rod are preserved exactly, both Type I and Type II chains have been recognised, and sequence identity/homology exists between their respective chains. Nonetheless, there are characteristic differences in amino acid composition and sequence between their respective head (N-terminal) domains and their tail (C-terminal) domains, though some similarities are retained. Further, there is evidence of tandem repeats of a variety of lengths in the tuatara heads and tails indicative of sequence duplication events. These are not evident in human α-keratins and would therefore have implications for the physical attributes of the tissues in the two species. Multiple families of keratin-associated proteins that are ultra-high cysteine-rich or glycine + tyrosine-rich in human and other species do not have direct equivalents in the tuatara. Although high-sulphur proteins are present in tuatara the cysteine residue contents are significantly lower than in human. Further, no sequence homologies between the HS proteins in the two species have been found, thereby casting considerable doubt as to whether any matrix-forming high-sulphur proteins exist in tuatara. These observations may be correlated with the numerous cysteine-rich ß-keratins (corneous ß-proteins) that are present in tuatara, but which are conspicuously absent in mammals.

Regulation of host-infection ability in the grass-symbiotic fungus Epichloë festucae by histone H3K9 and H3K36 methyltransferases.

Recent studies have identified key genes that control the symbiotic interaction between Epichloë festucae and Lolium perenne. Here we report on the identification of specific E. festucae genes that control host infection. Deletion of setB, which encodes a homologue of the H3K36 histone methyltransferase Set2/KMT3, reduced histone H3K36 trimethylation and led to severe defects in colony growth and hyphal development. The E. festucae ΔclrD mutant, which lacks the gene encoding the homologue of the H3K9 methyltransferase KMT1, displays similar developmental defects. Both mutants are completely defective in their ability to infect L. perenne. Alleles that complement the culture and plant phenotypes of both mutants also complement the histone methylation defects. Co-inoculation of either ΔsetB or ΔclrD with the wild-type strain enables these mutants to colonize the host. However, successful colonization by the mutants resulted in death or stunting of the host plant. Transcriptome analysis at the early infection stage identified four fungal candidate genes, three of which encode small-secreted proteins, that are differentially regulated in these mutants compared to wild type. Deletion of crbA, which encodes a putative carbohydrate binding protein, resulted in significantly reduced host infection rates by E. festucae.

Chromosome-Level Reference Genome of Venturia effusa, Causative Agent of Pecan Scab.

Pecan scab, caused by Venturia effusa, is a devastating disease of pecan (Carya illinoinensis), which results in economic losses on susceptible cultivars throughout the southeastern United States. To enhance our understanding of pathogenicity in V. effusa, we have generated a complete telomere-to-telomere reference genome of V. effusa isolate FRT5LL7-Albino. By combining Illumina MiSeq and Oxford Nanopore MinION data, we assembled a 45.2-Mb genome represented by 20 chromosomes and containing 10,820 putative genes, of which 7,619 have at least one functional annotation. The likely causative mutation of the albino phenotype was identified as a single base insertion and a resulting frameshift in the gene encoding the polyketide synthase ALM1. This genome represents the first full chromosome-level assembly of any Venturia sp.

Repeat elements organise 3D genome structure and mediate transcription in the filamentous fungus Epichloë festucae.

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism.

accuMUlate: a mutation caller designed for mutation accumulation experiments.

Summary: Mutation accumulation (MA) is the most widely used method for directly studying the effects of mutation. By sequencing whole genomes from MA lines, researchers can directly study the rate and molecular spectra of spontaneous mutations and use these results to understand how mutation contributes to biological processes. At present there is no software designed specifically for identifying mutations from MA lines. Here we describe accuMUlate, a probabilistic mutation caller that reflects the design of a typical MA experiment while being flexible enough to accommodate properties unique to any particular experiment. Availability and implementation accuMUlate is available from https://github.com/dwinter/accuMUlate. Supplementary information: Supplementary data are available at Bioinformatics online.

Interspecific hybridization causes long-term phylogenetic discordance between nuclear and mitochondrial genomes in freshwater fishes.

Classification, phylogeography and the testing of evolutionary hypotheses rely on correct estimation of species phylogeny. Early molecular phylogenies often relied on mtDNA alone, which acts as a single linkage group with one history. Over the last decade, the use of multiple nuclear sequences has often revealed conflict among gene trees. This observation can be attributed to hybridization, lineage sorting, paralogy or selection. Here, we use 54 groups of fishes from 48 studies to estimate the degree of concordance between mitochondrial and nuclear gene trees in two ecological grades of fishes: marine and freshwater. We test the hypothesis that freshwater fish phylogenies should, on average, show more discordance because of their higher propensity for hybridization in the past. In keeping with this idea, concordance between mitochondrial and nuclear gene trees (as measured by proportion of components shared) is on average 50% higher in marine fishes. We discuss why this difference almost certainly results from introgression caused by greater historical hybridization among lineages in freshwater groups, and further emphasize the need to use multiple nuclear genes, and identify conflict among them, in estimation of species phylogeny.

Estimating error models for whole genome sequencing using mixtures of Dirichlet-multinomial distributions.

MOTIVATION: Accurate identification of genotypes is an essential part of the analysis of genomic data, including in identification of sequence polymorphisms, linking mutations with disease and determining mutation rates. Biological and technical processes that adversely affect genotyping include copy-number-variation, paralogous sequences, library preparation, sequencing error and reference-mapping biases, among others. RESULTS: We modeled the read depth for all data as a mixture of Dirichlet-multinomial distributions, resulting in significant improvements over previously used models. In most cases the best model was comprised of two distributions. The major-component distribution is similar to a binomial distribution with low error and low reference bias. The minor-component distribution is overdispersed with higher error and reference bias. We also found that sites fitting the minor component are enriched for copy number variants and low complexity regions, which can produce erroneous genotype calls. By removing sites that do not fit the major component, we can improve the accuracy of genotype calls. AVAILABILITY AND IMPLEMENTATION: Methods and data files are available at https://github.com/CartwrightLab/WuEtAl2017/ (doi:10.5281/zenodo.256858). CONTACT: cartwright@asu.edu. SUPPLEMENTARY INFORMATION: Supplementary data is available at Bioinformatics online.

Low Base-Substitution Mutation Rate in the Germline Genome of the Ciliate Tetrahymena thermophil.

Mutation is the ultimate source of all genetic variation and is, therefore, central to evolutionary change. Previous work on Paramecium tetraurelia found an unusually low germline base-substitution mutation rate in this ciliate. Here, we tested the generality of this result among ciliates using Tetrahymena thermophila. We sequenced the genomes of 10 lines of T. thermophila that had each undergone approximately 1,000 generations of mutation accumulation (MA). We applied an existing mutation-calling pipeline and developed a new probabilistic mutation detection approach that directly models the design of an MA experiment and accommodates the noise introduced by mismapped reads. Our probabilistic mutation-calling method provides a straightforward way of estimating the number of sites at which a mutation could have been called if one was present, providing the denominator for our mutation rate calculations. From these methods, we find that T. thermophila has a germline base-substitution mutation rate of 7.61 × 10 - 12 per-site, per cell division, which is consistent with the low base-substitution mutation rate in P. tetraurelia. Over the course of the evolution experiment, genomic exclusion lines derived from the MA lines experienced a fitness decline that cannot be accounted for by germline base-substitution mutations alone, suggesting that other genetic or epigenetic factors must be involved. Because selection can only operate to reduce mutation rates based upon the "visible" mutational load, asexual reproduction with a transcriptionally silent germline may allow ciliates to evolve extremely low germline mutation rates.

Whole Genome Sequencing of Field Isolates Reveals Extensive Genetic Diversity in Plasmodium vivax from Colombia.

Plasmodium vivax is the most prevalent malarial species in South America and exerts a substantial burden on the populations it affects. The control and eventual elimination of P. vivax are global health priorities. Genomic research contributes to this objective by improving our understanding of the biology of P. vivax and through the development of new genetic markers that can be used to monitor efforts to reduce malaria transmission. Here we analyze whole-genome data from eight field samples from a region in Cordóba, Colombia where malaria is endemic. We find considerable genetic diversity within this population, a result that contrasts with earlier studies suggesting that P. vivax had limited diversity in the Americas. We also identify a selective sweep around a substitution known to confer resistance to sulphadoxine-pyrimethamine (SP). This is the first observation of a selective sweep for SP resistance in this species. These results indicate that P. vivax has been exposed to SP pressure even when the drug is not in use as a first line treatment for patients afflicted by this parasite. We identify multiple non-synonymous substitutions in three other genes known to be involved with drug resistance in Plasmodium species. Finally, we found extensive microsatellite polymorphisms. Using this information we developed 18 polymorphic and easy to score microsatellite loci that can be used in epidemiological investigations in South America.

MMOD: an R library for the calculation of population differentiation statistics.

MMOD is a library for the R programming language that allows the calculation of the population differentiation measures D(est), G″(ST) and φ'(ST). R provides a powerful environment in which to conduct and record population genetic analyses but, at present, no R libraries provide functions for the calculation of these statistics from standard population genetic files. In addition to the calculation of differentiation measures, mmod can produce parametric bootstrap and jackknife samples of data sets for further analysis. By integrating with and complimenting the existing libraries adegenet and pegas, mmod extends the power of R as a population genetic platform.

Multigene phylogeny of the southern bull-kelp genus Durvillaea (Phaeophyceae: Fucales).

Durvillaea (southern bull-kelp) is an economically and ecologically important brown algal genus that dominates many exposed, rocky coasts in the cold-temperate Southern Hemisphere. Of its five currently-recognized species, four are non-buoyant and restricted to the south-western Pacific, whereas one is both buoyant and widely distributed. Durvillaea has had an unsettled taxonomic history. Although its position within the brown algae (Phaeophyceae) has now been largely resolved through the use of molecular techniques, the taxonomic status of several Durvillaea species/morphotypes remains unresolved. Previous molecular phylogenetic studies of phaeophycean taxa have included few Durvillaea samples, and have consequently paid little or no attention to variation within this genus. The current study presents phylogenetic analyses of four genetic markers (mitchondrial: COI; chloroplast: rbcL; and nuclear: 18S and 28S) to resolve phylogenetic relationships within Durvillaea. Results support the monophyly of solid-bladed taxa D. willana, D. potatorum, and D. sp. A (an undescribed species from the Antipodes Islands), whereas the widespread, buoyant D. antarctica is paraphyletic, with solid-bladed D.chathamensis placed sister to a D. antarctica clade from northern NZ but within D. antarctica sensu lato. The phylogenetic and ecological diversity detected within D. antarctica indicate that it is a species complex of five deeply divergent clades. Under a phylogenetic species concept, Durvillaea can be interpreted as a complex of nine distinct evolutionary lineages, only one of which has an intercontinental distribution ('subantarctic'D. antarctica).