Comparative Genomic Analysis of 31 Phytophthora Genomes Reveals Genome Plasticity and Horizontal Gene Transfer.

Phytophthora species are oomycete plant pathogens that cause great economic and ecological impacts. The Phytophthora genus includes over 180 known species, infecting a wide range of plant hosts, including crops, trees, and ornamentals. We sequenced the genomes of 31 individual Phytophthora species and 24 individual transcriptomes to study genetic relationships across the genus. De novo genome assemblies revealed variation in genome sizes, numbers of predicted genes, and in repetitive element content across the Phytophthora genus. A genus-wide comparison evaluated orthologous groups of genes. Predicted effector gene counts varied across Phytophthora species by effector family, genome size, and plant host range. Predicted numbers of apoplastic effectors increased as the host range of Phytophthora species increased. Predicted numbers of cytoplasmic effectors also increased with host range but leveled off or decreased in Phytophthora species that have enormous host ranges. With extensive sequencing across the Phytophthora genus, we now have the genomic resources to evaluate horizontal gene transfer events across the oomycetes. Using a machine-learning approach to identify horizontally transferred genes with bacterial or fungal origin, we identified 44 candidates over 36 Phytophthora species genomes. Phylogenetic reconstruction indicates that the transfers of most of these 44 candidates happened in parallel to major advances in the evolution of the oomycetes and Phytophthora spp. We conclude that the 31 genomes presented here are essential for investigating genus-wide genomic associations in genus Phytophthora. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Long transposon-rich centromeres in an oomycete reveal divergence of centromere features in Stramenopila-Alveolata-Rhizaria lineages.

Centromeres are chromosomal regions that serve as platforms for kinetochore assembly and spindle attachments, ensuring accurate chromosome segregation during cell division. Despite functional conservation, centromere DNA sequences are diverse and often repetitive, making them challenging to assemble and identify. Here, we describe centromeres in an oomycete Phytophthora sojae by combining long-read sequencing-based genome assembly and chromatin immunoprecipitation for the centromeric histone CENP-A followed by high-throughput sequencing (ChIP-seq). P. sojae centromeres cluster at a single focus at different life stages and during nuclear division. We report an improved genome assembly of the P. sojae reference strain, which enabled identification of 15 enriched CENP-A binding regions as putative centromeres. By focusing on a subset of these regions, we demonstrate that centromeres in P. sojae are regional, spanning 211 to 356 kb. Most of these regions are transposon-rich, poorly transcribed, and lack the histone modification H3K4me2 but are embedded within regions with the heterochromatin marks H3K9me3 and H3K27me3. Strikingly, we discovered a Copia-like transposon (CoLT) that is highly enriched in the CENP-A chromatin. Similar clustered elements are also found in oomycete relatives of P. sojae, and may be applied as a criterion for prediction of oomycete centromeres. This work reveals a divergence of centromere features in oomycetes as compared to other organisms in the Stramenopila-Alveolata-Rhizaria (SAR) supergroup including diatoms and Plasmodium falciparum that have relatively short and simple regional centromeres. Identification of P. sojae centromeres in turn also advances the genome assembly.

Pan-tissue transcriptome analysis of long noncoding RNAs in the American beaver Castor canadensis.

BACKGROUND: Long noncoding RNAs (lncRNAs) have roles in gene regulation, epigenetics, and molecular scaffolding and it is hypothesized that they underlie some mammalian evolutionary adaptations. However, for many mammalian species, the absence of a genome assembly precludes the comprehensive identification of lncRNAs. The genome of the American beaver (Castor canadensis) has recently been sequenced, setting the stage for the systematic identification of beaver lncRNAs and the characterization of their expression in various tissues. The objective of this study was to discover and profile polyadenylated lncRNAs in the beaver using high-throughput short-read sequencing of RNA from sixteen beaver tissues and to annotate the resulting lncRNAs based on their potential for orthology with known lncRNAs in other species. RESULTS: Using de novo transcriptome assembly, we found 9528 potential lncRNA contigs and 187 high-confidence lncRNA contigs. Of the high-confidence lncRNA contigs, 147 have no known orthologs (and thus are putative novel lncRNAs) and 40 have mammalian orthologs. The novel lncRNAs mapped to the Oregon State University (OSU) reference beaver genome with greater than 90% sequence identity. While the novel lncRNAs were on average shorter than their annotated counterparts, they were similar to the annotated lncRNAs in terms of the relationships between contig length and minimum free energy (MFE) and between coverage and contig length. We identified beaver orthologs of known lncRNAs such as XIST, MEG3, TINCR, and NIPBL-DT. We profiled the expression of the 187 high-confidence lncRNAs across 16 beaver tissues (whole blood, brain, lung, liver, heart, stomach, intestine, skeletal muscle, kidney, spleen, ovary, placenta, castor gland, tail, toe-webbing, and tongue) and identified both tissue-specific and ubiquitous lncRNAs. CONCLUSIONS: To our knowledge this is the first report of systematic identification of lncRNAs and their expression atlas in beaver. LncRNAs-both novel and those with known orthologs-are expressed in each of the beaver tissues that we analyzed. For some beaver lncRNAs with known orthologs, the tissue-specific expression patterns were phylogenetically conserved. The lncRNA sequence data files and raw sequence files are available via the web supplement and the NCBI Sequence Read Archive, respectively.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest.

BACKGROUND: Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies. RESULTS: Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications. CONCLUSIONS: Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Transcriptome analysis of responses in Brachypodium distachyon overexpressing the BdbZIP26 transcription factor.

BACKGROUND: Biotic and abiotic stresses are the major cause of reduced growth, persistence, and yield in agriculture. Over the past decade, RNA-Sequencing and the use of transgenics with altered expression of stress related genes have been utilized to gain a better understanding of the molecular mechanisms leading to salt tolerance in a variety of species. Identification of transcription factors that, when overexpressed in plants, improve multiple stress tolerance may be valuable for crop improvement, but sometimes overexpression leads to deleterious effects during normal plant growth. RESULTS: Brachypodium constitutively expressing the BdbZIP26:GFP gene showed reduced stature compared to wild type plants (WT). RNA-Seq analysis comparing WT and bZIP26 transgenic plants revealed 7772 differentially expressed genes (DEGs). Of these DEGs, 987 of the DEGs were differentially expressed in all three transgenic lines. Many of these DEGs are similar to those often observed in response to abiotic and biotic stress, including signaling proteins such as kinases/phosphatases, calcium/calmodulin related proteins, oxidases/reductases, hormone production and signaling, transcription factors, as well as disease responsive proteins. Interestingly, there were many DEGs associated with protein turnover including ubiquitin-related proteins, F-Box and U-box related proteins, membrane proteins, and ribosomal synthesis proteins. Transgenic and control plants were exposed to salinity stress. Many of the DEGs between the WT and transgenic lines under control conditions were also found to be differentially expressed in WT in response to salinity stress. This suggests that the over-expression of the transcription factor is placing the plant in a state of stress, which may contribute to the plants diminished stature. CONCLUSION: The constitutive expression of BdbZIP26:GFP had an overall negative effect on plant growth and resulted in stunted plants compared to WT plants under control conditions, and a similar response to WT plants under salt stress conditions. The results of gene expression analysis suggest that the transgenic plants are in a constant state of stress, and that they are trying to allocate resources to survive.

Transcriptome analysis of the model grass Lolium temulentum exposed to green leaf volatiles.

BACKGROUND: Forage and turf grasses are routinely cut and grazed upon throughout their lifecycle. When grasses are cut or damaged, they rapidly release a volatile chemical cocktail called green leaf volatiles (GLV). Previously we have shown that mechanical wounding or exposure to GLV released from cut grass, activated a Lt 46 kDa mitogen-activated protein kinase (MAPK) within 3 min and a 44 kDa MAPK within 15-20 min in the model grass species Lolium temulentum (Lt). Currently very little is known concerning the perception, signaling or molecular responses associated with wound stress in grasses. Since GLV are released during wounding, we wanted to investigate what genes and signaling pathways would be induced in undamaged plants exposed to GLV. RESULTS: RNA-Seq generated transcriptome of Lolium plants exposed to GLV identified 4308 up- and 2794 down-regulated distinct differentially-expressed sequences (DES). Gene Ontology analysis revealed a strong emphasis on signaling, response to stimulus and stress related categories. Transcription factors and kinases comprise over 13% of the total DES found in the up-regulated dataset. The analysis showed a strong initial burst within the first hour of GLV exposure with over 60% of the up-regulated DES being induced. Specifically sequences annotated for enzymes involved in the biosynthesis of jasmonic acid and other plant hormones, mitogen-activated protein kinases and WRKY transcription factors were identified. Interestingly, eleven DES for ferric reductase oxidase, an enzyme involved in iron uptake and transport, were exclusively found in the down-regulated dataset. Twelve DES of interest were selected for qRT-PCR analysis; all displayed a rapid induction one hour after GLV exposure and were also strongly induced by mechanical wounding. CONCLUSION: The information gained from the analysis of this transcriptome and previous studies suggests that GLV released from cut grasses transiently primes an undamaged plant's wound stress pathways for potential oncoming damage, and may have a dual role for inter- as well as intra-plant signaling.

The Antibiotic Resistant Target Seeker (ARTS), an exploration engine for antibiotic cluster prioritization and novel drug target discovery.

With the rise of multi-drug resistant pathogens and the decline in number of potential new antibiotics in development there is a fervent need to reinvigorate the natural products discovery pipeline. Most antibiotics are derived from secondary metabolites produced by microorganisms and plants. To avoid suicide, an antibiotic producer harbors resistance genes often found within the same biosynthetic gene cluster (BGC) responsible for manufacturing the antibiotic. Existing mining tools are excellent at detecting BGCs or resistant genes in general, but provide little help in prioritizing and identifying gene clusters for compounds active against specific and novel targets. Here we introduce the 'Antibiotic Resistant Target Seeker' (ARTS) available at https://arts.ziemertlab.com. ARTS allows for specific and efficient genome mining for antibiotics with interesting and novel targets. The aim of this web server is to automate the screening of large amounts of sequence data and to focus on the most promising strains that produce antibiotics with new modes of action. ARTS integrates target directed genome mining methods, antibiotic gene cluster predictions and 'essential gene screening' to provide an interactive page for rapid identification of known and putative targets in BGCs.

Whole Genome Sequences of the Raspberry and Strawberry Pathogens Phytophthora rubi and P. fragariae.

Phytophthora rubi and P. fragariae are two closely related oomycete plant pathogens that exhibit strong morphological and physiological similarities but are specialized to infect different hosts of economic importance, namely, raspberry and strawberry. Here, we report the draft genome sequences of these two Phytophthora species as a first step toward understanding the genomic processes underlying plant host adaptation in these pathogens.

Essential role for the response regulator PmrA in Coxiella burnetii type 4B secretion and colonization of mammalian host cells.

Successful host cell colonization by the Q fever pathogen, Coxiella burnetii, requires translocation of effector proteins into the host cytosol by a Dot/Icm type 4B secretion system (T4BSS). In Legionella pneumophila, the two-component system (TCS) PmrAB regulates the Dot/Icm T4BSS and several additional physiological processes associated with pathogenesis. Because PmrA consensus regulatory elements are associated with some dot/icm and substrate genes, a similar role for PmrA in regulation of the C. burnetii T4BSS has been proposed. Here, we constructed a C. burnetii pmrA deletion mutant to directly probe PmrA-mediated gene regulation. Compared to wild-type bacteria, C. burnetii ΔpmrA exhibited severe intracellular growth defects that coincided with failed secretion of effector proteins. Luciferase gene reporter assays demonstrated PmrA-dependent expression of 5 of 7 dot/icm operons and 9 of 11 effector-encoding genes with a predicted upstream PmrA regulatory element. Mutational analysis verified consensus sequence nucleotides required for PmrA-directed transcription. RNA sequencing and whole bacterial cell mass spectrometry of wild-type C. burnetii and the ΔpmrA mutant uncovered new components of the PmrA regulon, including several genes lacking PmrA motifs that encoded Dot/Icm substrates. Collectively, our results indicate that the PmrAB TCS is a critical virulence factor that regulates C. burnetii Dot/Icm secretion. The presence of PmrA-responsive genes lacking PmrA regulatory elements also suggests that the PmrAB TCS controls expression of regulatory systems associated with the production of additional C. burnetii proteins involved in host cell parasitism.

RNA sequencing least shrew (Cryptotis parva) brainstem and gut transcripts following administration of a selective substance P neurokinin NK1 receptor agonist and antagonist expands genomics resources for emesis research.

The least shrew is among the subset of animals that are capable of vomiting and therefore serves as a valuable research model for investigating the biochemistry, molecular biology, pharmacology, and genomics of emesis. Both nausea and vomiting are associated with a variety of illnesses (bacterial/viral infections, bulimia, exposure to toxins, gall bladder disease), conditions (pregnancy, motion sickness, emotional stress, overeating) and reactions to drugs (chemotherapeutics, opiates). The severe discomfort and intense fear associated with the stressful symptoms of nausea and emesis are the major reason for patient non-compliance when being treated with cancer chemotherapeutics. Increased understanding of the physiology, pharmacology and pathophysiology underlying vomiting and nausea can accelerate progress for developing new antiemetics. As a major animal model for emesis, expanding genomic knowledge associated with emesis in the least shrew will further enhance the laboratory utility of this model. A key question is which genes mediate emesis, and are they expressed in response to emetics/antiemetics. To elucidate the mediators of emesis, in particular emetic receptors, their downstream signaling pathways, as well as the shared emetic signals, we carried out an RNA sequencing study focused on the central and peripheral emetic loci, the brainstem and gut. Thus, we sequenced RNA extracted from brainstem and gut tissues from different groups of least shrews treated with either a neurokinin NK1 receptor selective emetic agonist, GR73632 (5 mg/kg, i.p.), its corresponding selective antagonist netupitant (5 mg/kg, i.p.), a combination of these two agents, versus their corresponding vehicle-pretreated controls and drug naïve animals. The resulting sequences were processed using a de novo transcriptome assembly and used it to identify orthologs within human, dog, mouse, and ferret gene sets. We compared the least shrew to human and a veterinary species (dog) that may be treated with vomit-inducing chemotherapeutics, and the ferret, another well-established model organism for emesis research. The mouse was included because it does not vomit. In total, we identified a final set of 16,720 least shrew orthologs. We employed comparative genomics analyses as well as gene ontology enrichment, KEGG pathway enrichment and phenotype enrichment to better understand the molecular biology of genes implicated in vomiting.

SSR-based genetic maps of Miscanthus sinensis and M. sacchariflorus, and their comparison to sorghum.

We present SSR-based genetic maps from a cross between Miscanthus sacchariflorus Robustus and M. sinensis, the progenitors of the promising cellulosic biofuel feedstock Miscanthus × giganteus. cDNA-derived SSR markers were mapped by the two-way pseudo-testcross model due to the high heterozygosity of each parental species. A total of 261 loci were mapped in M. sacchariflorus, spanning 40 linkage groups and 1,998.8 cM, covering an estimated 72.7% of the genome. For M. sinensis, a total of 303 loci were mapped, forming 23 linkage groups and 2,238.3 cM, covering 84.9% of the genome. The use of cDNA-derived SSR loci permitted alignment of the Miscanthus linkage groups to the sorghum chromosomes, revealing a whole genome duplication affecting the Miscanthus lineage after the divergence of subtribes Sorghinae and Saccharinae, as well as traces of the pan-cereal whole genome duplication. While the present maps provide for many early research needs in this emerging crop, additional markers are also needed to improve map density and to further characterize the structural changes of the Miscanthus genome since its divergence from sorghum and Saccharum.

Chromosome-level genome assembly of Mentha longifolia L. reveals gene organization underlying disease resistance and essential oil traits.

Mentha longifolia (L.) Huds., a wild, diploid mint species, has been developed as a model for mint genetic and genomic research to aid breeding efforts that target Verticillium wilt disease resistance and essential oil monoterpene composition. Here, we present a near-complete, chromosome-scale mint genome assembly for M. longifolia USDA accession CMEN 585. This new assembly is an update of a previously published genome draft, with dramatic improvements. A total of 42,107 protein-coding genes were annotated and placed on 12 chromosomal scaffolds. One hundred fifty-three genes contained conserved sequence domains consistent with nucleotide binding site-leucine-rich-repeat plant disease resistance genes. Homologs of genes implicated in Verticillium wilt resistance in other plant species were also identified. Multiple paralogs of genes putatively involved in p-menthane monoterpenoid biosynthesis were identified and several cases of gene clustering documented. Heterologous expression of candidate genes, purification of recombinant target proteins, and subsequent enzyme assays allowed us to identify the genes underlying the pathway that leads to the most abundant monoterpenoid volatiles. The bioinformatic and functional analyses presented here are laying the groundwork for using marker-assisted selection in improving disease resistance and essential oil traits in mints.

Transgenic Soybeans Expressing Phosphatidylinositol-3-Phosphate-Binding Proteins Show Enhanced Resistance Against the Oomycete Pathogen Phytophthora sojae.

Oomycete and fungal pathogens cause billions of dollars of damage to crops worldwide annually. Therefore, there remains a need for broad-spectrum resistance genes, especially ones that target pathogens but do not interfere with colonization by beneficial microbes. Motivated by evidence suggesting that phosphatidylinositol-3-phosphate (PI3P) may be involved in the delivery of some oomycete and fungal virulence effector proteins, we created stable transgenic soybean plants that express and secrete two different PI3P-binding proteins, GmPH1 and VAM7, in an effort to interfere with effector delivery and confer resistance. Soybean plants expressing the two PI3P-binding proteins exhibited reduced infection by the oomycete pathogen Phytophthora sojae compared to control lines. Measurements of nodulation by nitrogen-fixing mutualistic bacterium Bradyrhizobium japonicum, which does not produce PI3P, revealed that the two lines with the highest levels of GmPH1 transcripts exhibited reductions in nodulation and in benefits from nodulation. Transcriptome and plant hormone measurements were made of soybean lines with the highest transcript levels of GmPH1 and VAM7, as well as controls, following P. sojae- or mock-inoculation. The results revealed increased levels of infection-associated transcripts in the transgenic lines, compared to controls, even prior to P. sojae infection, suggesting that the plants were primed for increased defense. The lines with reduced nodulation exhibited elevated levels of jasmonate-isoleucine and of transcripts of a JAR1 ortholog encoding jasmonate-isoleucine synthetase. However, lines expressing VAM7 transgenes exhibited normal nodulation and no increases in jasmonate-isoleucine. Overall, together with previously published data from cacao and from P. sojae transformants, the data suggest that secretion of PI3P-binding proteins may confer disease resistance through a variety of mechanisms.

Transcriptome Analysis of Lolium temulentum Exposed to a Combination of Drought and Heat Stress.

Drought and heat are two major stresses predicted to increase in the future due to climate change. Plants exposed to multiple stressors elicit unique responses from those observed under individual stresses. A comparative transcriptome analysis of Lolium temulentum exposed to drought plus heat and non-stressed control plants revealed 20,221 unique up-regulated and 17,034 unique down-regulated differentially regulated transcripts. Gene ontology analysis revealed a strong emphasis on transcriptional regulation, protein folding, cell cycle/parts, organelles, binding, transport, signaling, oxidoreductase, and antioxidant activity. Differentially expressed genes (DEGs) encoding for transcriptional control proteins such as basic leucine zipper, APETALA2/Ethylene Responsive Factor, NAC, and WRKY transcription factors, and Zinc Finger (CCCH type and others) proteins were more often up-regulated, while DEGs encoding Basic Helix-Loop-Helix, MYB and GATA transcription factors, and C2H2 type Zinc Finger proteins were more often down-regulated. The DEGs encoding heat shock transcription factors were only up-regulated. Of the hormones, auxin-related DEGs were the most prevalent, encoding for auxin response factors, binding proteins, and efflux/influx carriers. Gibberellin-, cytokinin- and ABA-related DEGs were also prevalent, with fewer DEGs related to jasmonates and brassinosteroids. Knowledge of genes/pathways that grasses use to respond to the combination of heat/drought will be useful in developing multi-stress resistant grasses.

A draft phased assembly of the diploid Cascade hop (Humulus lupulus) genome.

Hop (Humulus lupulus L. var Lupulus) is a diploid, dioecious plant with a history of cultivation spanning more than one thousand years. Hop cones are valued for their use in brewing and contain compounds of therapeutic interest including xanthohumol. Efforts to determine how biochemical pathways responsible for desirable traits are regulated have been challenged by the large (2.8 Gb), repetitive, and heterozygous genome of hop. We present a draft haplotype-phased assembly of the Cascade cultivar genome. Our draft assembly and annotation of the Cascade genome is the most extensive representation of the hop genome to date. PacBio long-read sequences from hop were assembled with FALCON and partially phased with FALCON-Unzip. Comparative analysis of haplotype sequences provides insight into selective pressures that have driven evolution in hop. We discovered genes with greater sequence divergence enriched for stress-response, growth, and flowering functions in the draft phased assembly. With improved resolution of long terminal retrotransposons (LTRs) due to long-read sequencing, we found that hop is over 70% repetitive. We identified a homolog of cannabidiolic acid synthase (CBDAS) that is expressed in multiple tissues. The approaches we developed to analyze the draft phased assembly serve to deepen our understanding of the genomic landscape of hop and may have broader applicability to the study of other large, complex genomes.

Transcriptome Analysis of Wounding in the Model Grass Lolium temulentum.

For forage and turf grasses, wounding is a predominant stress that often results in extensive loss of vegetative tissues followed by rapid regrowth. Currently, little is known concerning the perception, signaling, or molecular responses associated with wound stress in forage- and turf-related grasses. A transcriptome analysis of Lolium temulentum plants subjected to severe wounding revealed 9413 upregulated and 7704 downregulated, distinct, differentially expressed genes (DEGs). Categories related to signaling, transcription, and response to stimuli were enriched in the upregulated DEGs. Specifically, sequences annotated as enzymes involved in hormone biosynthesis/action and cell wall modifications, mitogen-activated protein kinases, WRKY transcription factors, proteinase inhibitors, and pathogen defense-related DEGs were identified. Surprisingly, DEGs related to heat shock and chaperones were more prevalent in the downregulated DEGs when compared with the upregulated DEGs. This wound transcriptome analysis is the first step in identifying the molecular components and pathways used by grasses in response to wounding. The information gained from the analysis will provide a valuable molecular resource that will be used to develop approaches that can improve the recovery, regrowth, and long-term fitness of forage and turf grasses before/after cutting or grazing.

Clusters of polymorphic transmembrane genes control resistance to schistosomes in snail vectors.

Schistosomiasis is a debilitating parasitic disease infecting hundreds of millions of people. Schistosomes use aquatic snails as intermediate hosts. A promising avenue for disease control involves leveraging innate host mechanisms to reduce snail vectorial capacity. In a genome-wide association study of Biomphalaria glabrata snails, we identify genomic region PTC2 which exhibits the largest known correlation with susceptibility to parasite infection (>15 fold effect). Using new genome assemblies with substantially higher contiguity than the Biomphalaria reference genome, we show that PTC2 haplotypes are exceptionally divergent in structure and sequence. This variation includes multi-kilobase indels containing entire genes, and orthologs for which most amino acid residues are polymorphic. RNA-Seq annotation reveals that most of these genes encode single-pass transmembrane proteins, as seen in another resistance region in the same species. Such groups of hyperdiverse snail proteins may mediate host-parasite interaction at the cell surface, offering promising targets for blocking the transmission of schistosomiasis.

Schistosomiasis is a widespread parasitic disease, affecting over 200 million people in tropical countries. It is caused by schistosome worms, which are carried by freshwater snails. These snails release worm larvae into the water, where they can infect humans ­ for example, after bathing or swimming. Treatment options for schistosomiasis are limited. Eliminating the freshwater snails is one way to control the disease, but this is not always effective in the long term and the chemicals used can also harm other animals in the water. Another way to manage schistosomiasis could be to stop the worms from infecting their snail host by breaking the parasites' life cycle without killing the snails. It is already known that some snails are naturally resistant to infection by some strains of schistosomes. Since this immunity is also inherited by the offspring of resistant snails, there is likely a genetic mechanism behind it. However, very little else is known about any genes that might be involved. Tennessen et al. therefore set out to identify what genes were responsible for schistosome resistance and how they worked. The experiments used a large laboratory colony of snails, whose susceptibility to schistosome infection varied among individual animals. To determine the genes behind this variation, Tennessen et al. first searched for areas of DNA that also differed between the immune and infected snails. Comparing genetic sequences across over 1,000 snails revealed a distinct region of DNA that had a large effect on how likely they were to be infected. This section of DNA turned out to be highly diverse, with different snails carrying varying numbers and different forms of the genes within this region. Many of these genes appear to encode proteins found on the surface of snail cells, which could affect whether snails and worms can recognize each other when they come into contact. This in turn could determine whether or not the worms can infect their hosts. These results shed new light on how the snails that carry schistosomes may be able to resist infections. In the future, this knowledge could be key to controlling schistosomiasis, either by releasing genetically engineered, immune snails into the wild (thus making it harder for the parasites to reproduce) or by using the snails' mechanism of resistance to design better drug therapies.

Genomic Investigation of the Strawberry Pathogen Phytophthora fragariae Indicates Pathogenicity Is Associated With Transcriptional Variation in Three Key Races.

The oomycete Phytophthora fragariae is a highly destructive pathogen of cultivated strawberry (Fragaria × ananassa), causing the root rotting disease, "red core". The host-pathogen interaction has a well described gene-for-gene resistance relationship, but to date neither candidate avirulence nor resistance genes have been identified. We sequenced a set of American, Canadian, and United Kingdom isolates of known race type, along with three representatives of the closely related pathogen of the raspberry (Rubus idaeus), P. rubi, and found a clear population structure, with a high degree of nucleotide divergence seen between some race types and abundant private variation associated with race types 4 and 5. In contrast, between isolates defined as United Kingdom races 1, 2, and 3 (UK1-2-3) there was no evidence of gene loss or gain; or the presence of insertions/deletions (INDELs) or Single Nucleotide Polymorphisms (SNPs) within or in proximity to putative pathogenicity genes could be found associated with race variation. Transcriptomic analysis of representative UK1-2-3 isolates revealed abundant expression variation in key effector family genes associated with pathogen race; however, further long read sequencing did not reveal any long range polymorphisms to be associated with avirulence to race UK2 or UK3 resistance, suggesting either control in trans or other stable forms of epigenetic modification modulating gene expression. This work reveals the combined power of population resequencing to uncover race structure in pathosystems and in planta transcriptomic analysis to identify candidate avirulence genes. This work has implications for the identification of putative avirulence genes in the absence of associated expression data and points toward the need for detailed molecular characterisation of mechanisms of effector regulation and silencing in oomycete plant pathogens.

Not so Dangerous After All? Venom Composition and Potency of the Pholcid (Daddy Long-Leg) Spider Physocyclus mexicanus.

Pholcid spiders (Araneae: Pholcidae), officially "cellar spiders" but popularly known as "daddy long-legs," are renown for the potential of deadly toxic venom, even though venom composition and potency has never formally been studied. Here we detail the venom composition of male Physocyclus mexicanus using proteomic analyses and venom-gland transcriptomes ("venomics"). We also analyze the venom's potency on insects, and assemble available evidence regarding mammalian toxicity. The majority of the venom (51% of tryptic polypeptides and 62% of unique tryptic peptides) consists of proteins homologous to known venom toxins including enzymes (astacin metalloproteases, serine proteases and metalloendopeptidases, particularly neprilysins) and venom peptide neurotoxins. We identify 17 new groups of peptides (U1-17-PHTX) most of which are homologs of known venom peptides and are predicted to have an inhibitor cysteine knot fold; of these, 13 are confirmed in the proteome. Neprilysins (M13 peptidases), and astacins (M12 peptidases) are the most abundant venom proteins, respectively representing 15 and 11% of the individual proteins and 32 and 20% of the tryptic peptides detected in crude venom. Comparative evidence suggests that the neprilysin gene family is expressed in venoms across a range of spider taxa, but has undergone an expansion in the venoms of pholcids and may play a central functional role in these spiders. Bioassays of crude venoms on crickets resulted in an effective paralytic dose of 3.9 µg/g, which is comparable to that of crude venoms of Plectreurys tristis and other Synspermiata taxa. However, crickets exhibit flaccid paralysis and regions of darkening that are not observed after P. tristis envenomation. Documented bites on humans make clear that while these spiders can bite, the typical result is a mild sting with no long-lasting effects. Together, the evidence we present indicates pholcid venoms are a source of interesting new peptides and proteins, and effects of bites on humans and other mammals are inconsequential.

A taxonomic monograph of Ipomoea integrated across phylogenetic scales.

Taxonomic monographs have the potential to make a unique contribution to the understanding of global biodiversity. However, such studies, now rare, are often considered too daunting to undertake within a realistic time frame, especially as the world's collections have doubled in size in recent times. Here, we report a global-scale monographic study of morning glories (Ipomoea) that integrated DNA barcodes and high-throughput sequencing with the morphological study of herbarium specimens. Our approach overhauled the taxonomy of this megadiverse group, described 63 new species and uncovered significant increases in net diversification rates comparable to the most iconic evolutionary radiations in the plant kingdom. Finally, we show that more than 60 species of Ipomoea, including sweet potato, independently evolved storage roots in pre-human times, indicating that the storage root is not solely a product of human domestication but a trait that predisposed the species for cultivation. This study demonstrates how the world's natural history collections can contribute to global challenges in the Anthropocene.