Targeted Mutagenesis of the Multicopy Myrosinase Gene Family in Allotetraploid Brassica juncea Reduces Pungency in Fresh Leaves across Environments.

Recent breeding efforts in Brassica have focused on the development of new oilseed feedstock crop for biofuels (e.g., ethanol, biodiesel, bio-jet fuel), bio-industrial uses (e.g., bio-plastics, lubricants), specialty fatty acids (e.g., erucic acid), and producing low glucosinolates levels for oilseed and feed meal production for animal consumption. We identified a novel opportunity to enhance the availability of nutritious, fresh leafy greens for human consumption. Here, we demonstrated the efficacy of disarming the 'mustard bomb' reaction in reducing pungency upon the mastication of fresh tissue-a major source of unpleasant flavor and/or odor in leafy Brassica. Using gene-specific mutagenesis via CRISPR-Cas12a, we created knockouts of all functional copies of the type-I myrosinase multigene family in tetraploid Brassica juncea. Our greenhouse and field trials demonstrate, via sensory and biochemical analyses, a stable reduction in pungency in edited plants across multiple environments. Collectively, these efforts provide a compelling path toward boosting the human consumption of nutrient-dense, fresh, leafy green vegetables.

Social network analysis of the genealogy of strawberry: retracing the wild roots of heirloom and modern cultivars.

The widely recounted story of the origin of cultivated strawberry (Fragaria × ananassa) oversimplifies the complex interspecific hybrid ancestry of the highly admixed populations from which heirloom and modern cultivars have emerged. To develop deeper insights into the three-century-long domestication history of strawberry, we reconstructed the genealogy as deeply as possible-pedigree records were assembled for 8,851 individuals, including 2,656 cultivars developed since 1775. The parents of individuals with unverified or missing pedigree records were accurately identified by applying an exclusion analysis to array-genotyped single-nucleotide polymorphisms. We identified 187 wild octoploid and 1,171 F. × ananassa founders in the genealogy, from the earliest hybrids to modern cultivars. The pedigree networks for cultivated strawberry are exceedingly complex labyrinths of ancestral interconnections formed by diverse hybrid ancestry, directional selection, migration, admixture, bottlenecks, overlapping generations, and recurrent hybridization with common ancestors that have unequally contributed allelic diversity to heirloom and modern cultivars. Fifteen to 333 ancestors were predicted to have transmitted 90% of the alleles found in country-, region-, and continent-specific populations. Using parent-offspring edges in the global pedigree network, we found that selection cycle lengths over the past 200 years of breeding have been extraordinarily long (16.0-16.9 years/generation), but decreased to a present-day range of 6.0-10.0 years/generation. Our analyses uncovered conspicuous differences in the ancestry and structure of North American and European populations, and shed light on forces that have shaped phenotypic diversity in F. × ananassa.

Author Correction: Origin and evolution of the octoploid strawberry genome.

In the version of this article originally published, author Joshua R. Puzey was incorrectly listed as having affiliation 7 (School of Plant Sciences, University of Arizona, Tucson, AZ, USA); affiliation 6 (Department of Biology, College of William and Mary, Williamsburg, VA, USA) is the correct affiliation. The error has been corrected in the HTML and PDF versions of the article.

Origin and evolution of the octoploid strawberry genome.

Cultivated strawberry emerged from the hybridization of two wild octoploid species, both descendants from the merger of four diploid progenitor species into a single nucleus more than 1 million years ago. Here we report a near-complete chromosome-scale assembly for cultivated octoploid strawberry (Fragaria × ananassa) and uncovered the origin and evolutionary processes that shaped this complex allopolyploid. We identified the extant relatives of each diploid progenitor species and provide support for the North American origin of octoploid strawberry. We examined the dynamics among the four subgenomes in octoploid strawberry and uncovered the presence of a single dominant subgenome with significantly greater gene content, gene expression abundance, and biased exchanges between homoeologous chromosomes, as compared with the other subgenomes. Pathway analysis showed that certain metabolomic and disease-resistance traits are largely controlled by the dominant subgenome. These findings and the reference genome should serve as a powerful platform for future evolutionary studies and enable molecular breeding in strawberry.

Domestication of Temperate and Coastal Hybrids with Distinct Ancestral Gene Selection in Octoploid Strawberry.

Garden strawberry ( × Duchesne ex Rozier) arose from spontaneous hybridization of distinct octoploid species 300 yr ago. Since its discovery in the 1700s, migration and selection restructured the genetic diversity of early hybrids to produce elite fruit-bearing groups. Breeders' understanding of the genetic architecture of domesticated populations is incomplete. To resolve the impacts of domestication on strawberry genetic diversity, we analyzed genome-wide DNA profiles of 1300 octoploid individuals (1814-present), including wild species, historic varieties, and the University of California germplasm collection. Commercially important California genotypes, adapted to mild coastal climates and accounting for a large fraction of global production, have diverged from temperate cultivars originating in eastern North America and Europe. Whereas temperate cultivars were shown to have selected North American Miller ssp. ancestral diversity at higher frequencies, coastal breeding increased selection of (L.) Miller (beach strawberry) alleles in . × , in addition to photoperiod-insensitive flowering alleles from nonancestral (S.Watson) Staudt ssp. , underscoring the role of continued adaptive introgressions in the domestication of artificial hybrids. Selection for mass production traits in coastal climates over the last 20 to 30 yr has restructured domesticated strawberry diversity on a scale similar to the first 200 yr of breeding; coastal × has diverged further from temperate × than the latter from their wild progenitors. Selection signatures indicate that strawberry domestication targeted genes regulating hormone-mediated fruit expansion, providing a blueprint for genetic factors underlying elite phenotypes.

Opening the file drawer: Unexpected insights from a chytrid infection experiment.

Infection experiments are critical for understanding wildlife disease dynamics. Although infection experiments are typically designed to reduce complexity, disease outcomes still result from complex interactions between host, pathogen, and environmental factors. Cryptic variation across factors can lead to decreased repeatability of infection experiments within and between research groups and hinder research progress. Furthermore, studies with unexpected results are often relegated to the "file drawer" and potential insights gained from these experimental outcomes are lost. Here, we report unexpected results from an infection experiment studying the response of two differentially-susceptible but related frogs (American Bullfrog Rana catesbeiana and the Mountain yellow-legged frog Rana muscosa) to the amphibian-killing chytrid fungus (Batrachochytrium dendrobatidis, Bd). Despite well-documented differences in susceptibility between species, we found no evidence for antibody-mediated immune response and no Bd-related mortality in either species. Additionally, during the study, the sham-inoculated R. catesbeiana control group became unexpectedly Bd-positive. We used a custom genotyping assay to demonstrate that the aberrantly-infected R. catesbeiana carried a Bd genotype distinct from the inoculation genotype. Thus R. catesbeiana individuals were acquired with low-intensity infections that could not be detected with qPCR. In the Bd-inoculated R. catesbeiana treatment group, the inoculated genotype appeared to out-compete the cryptic infection. Thus, our results provide insight into Bd coinfection dynamics, a phenomenon that is increasingly relevant as different pathogen strains are moved around the globe. Our experiment highlights how unexpected experimental outcomes can serve as both cautionary tales and opportunities to explore unanswered research questions. We use our results as a case study to highlight common sources of anomalous results for infection experiments. We argue that understanding these factors will aid researchers in the design, execution, and interpretation of experiments to understand wildlife disease processes.

Genome-Wide Association Mapping Uncovers Fw1, a Dominant Gene Conferring Resistance to Fusarium Wilt in Strawberry.

Fusarium wilt, a soil-borne disease caused by the fungal pathogen Fusarium oxysporum f. sp. fragariae, threatens strawberry (Fragaria × ananassa) production worldwide. The spread of the pathogen, coupled with disruptive changes in soil fumigation practices, have greatly increased disease pressure and the importance of developing resistant cultivars. While resistant and susceptible cultivars have been reported, a limited number of germplasm accessions have been analyzed, and contradictory conclusions have been reached in earlier studies to elucidate the underlying genetic basis of resistance. Here, we report the discovery of Fw1, a dominant gene conferring resistance to Fusarium wilt in strawberry. The Fw1 locus was uncovered in a genome-wide association study of 565 historically and commercially important strawberry accessions genotyped with 14,408 SNP markers. Fourteen SNPs in linkage disequilibrium with Fw1 physically mapped to a 2.3 Mb segment on chromosome 2 in a diploid F. vesca reference genome. Fw1 and 11 tightly linked GWAS-significant SNPs mapped to linkage group 2C in octoploid segregating populations. The most significant SNP explained 85% of the phenotypic variability and predicted resistance in 97% of the accessions tested-broad-sense heritability was 0.96. Several disease resistance and defense-related gene homologs, including a small cluster of genes encoding nucleotide-binding leucine-rich-repeat proteins, were identified in the 0.7 Mb genomic segment predicted to harbor Fw1 DNA variants and candidate genes identified in the present study should facilitate the development of high-throughput genotyping assays for accurately predicting Fusarium wilt phenotypes and applying marker-assisted selection.

Single-molecule sequencing and optical mapping yields an improved genome of woodland strawberry (Fragaria vesca) with chromosome-scale contiguity.

Background: Although draft genomes are available for most agronomically important plant species, the majority are incomplete, highly fragmented, and often riddled with assembly and scaffolding errors. These assembly issues hinder advances in tool development for functional genomics and systems biology. Findings: Here we utilized a robust, cost-effective approach to produce high-quality reference genomes. We report a near-complete genome of diploid woodland strawberry (Fragaria vesca) using single-molecule real-time sequencing from Pacific Biosciences (PacBio). This assembly has a contig N50 length of ∼7.9 million base pairs (Mb), representing a ∼300-fold improvement of the previous version. The vast majority (>99.8%) of the assembly was anchored to 7 pseudomolecules using 2 sets of optical maps from Bionano Genomics. We obtained ∼24.96 Mb of sequence not present in the previous version of the F. vesca genome and produced an improved annotation that includes 1496 new genes. Comparative syntenic analyses uncovered numerous, large-scale scaffolding errors present in each chromosome in the previously published version of the F. vesca genome. Conclusions: Our results highlight the need to improve existing short-read based reference genomes. Furthermore, we demonstrate how genome quality impacts commonly used analyses for addressing both fundamental and applied biological questions.

Unlocking the story in the swab: A new genotyping assay for the amphibian chytrid fungus Batrachochytrium dendrobatidis.

One of the most devastating emerging pathogens of wildlife is the chytrid fungus, Batrachochytrium dendrobatidis (Bd), which affects hundreds of amphibian species around the world. Genomic data from pure Bd cultures have advanced our understanding of Bd phylogenetics, genomic architecture and mechanisms of virulence. However, pure cultures are laborious to obtain and whole-genome sequencing is comparatively expensive, so relatively few isolates have been genetically characterized. Thus, we still know little about the genetic diversity of Bd in natural systems. The most common noninvasive method of sampling Bd from natural populations is to swab amphibian skin. Hundreds of thousands of swabs have been collected from amphibians around the world, but Bd DNA collected via swabs is often low in quality and/or quantity. In this study, we developed a custom Bd genotyping assay using the Fluidigm Access Array platform to amplify 192 carefully selected regions of the Bd genome. We obtained robust sequence data for pure Bd cultures and field-collected skin swabs. This new assay has the power to accurately discriminate among the major Bd clades, recovering the basic tree topology previously revealed using whole-genome data. Additionally, we established a critical value for initial Bd load for swab samples (150 Bd genomic equivalents) above which our assay performs well. By leveraging advances in microfluidic multiplex PCR technology and the globally distributed resource of amphibian swab samples, noninvasive skin swabs can now be used to address critical spatial and temporal questions about Bd and its effects on declining amphibian populations.

Mountain Yellow-legged Frogs (Rana muscosa) did not Produce Detectable Antibodies in Immunization Experiments with Batrachochytrium dendrobatidis.

Chytridiomycosis is a devastating infectious disease of amphibians caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). A growing number of studies have examined the role of amphibian adaptive immunity in response to this pathogen, with varying degrees of immune activation reported. Here we present immunologic data for the mountain yellow-legged frog, Rana muscosa, and the Sierra Nevada yellow-legged frog, Rana sierrae, which are two endangered and ecologically important species experiencing Bd-inflicted declines. Previous studies on these species that examined transcriptional response during Bd infection, and the effective of immunization, provided little evidence of immune activation to Bd. However, the studies did not directly assay immune effectors in the frog hosts. We performed experiments to examine antibody production, which is a hallmark of systemic adaptive immune activation. We used controlled laboratory experiments and enzyme-linked immunosorbent assays to examine the antibody response to Bd immunization and live Bd exposure. Rana muscosa and R. sierrae individuals did not produce detectable antibodies with the capacity to bind to denatured Bd antigens under our experimental conditions. While we cannot rule out antibody response to Bd in these species, our results suggest weak, poor, or inefficient production of antibodies to denatured Bd antigens. Our findings are consistent with susceptibility to chytridiomycosis in these species and suggest additional work is needed to characterize the potential for adaptive immunity.

Genomic Correlates of Virulence Attenuation in the Deadly Amphibian Chytrid Fungus, Batrachochytrium dendrobatidis.

Emerging infectious diseasespose a significant threat to global health, but predicting disease outcomes for particular species can be complicated when pathogen virulence varies across space, time, or hosts. The pathogenic chytrid fungus Batrachochytrium dendrobatidis (Bd) has caused worldwide declines in frog populations. Not only do Bd isolates from wild populations vary in virulence, but virulence shifts can occur over short timescales when Bd is maintained in the laboratory. We leveraged changes in Bd virulence over multiple generations of passage to better understand mechanisms of pathogen virulence. We conducted whole-genome resequencing of two samples of the same Bd isolate, differing only in passage history, to identify genomic processes associated with virulence attenuation. The isolate with shorter passage history (and greater virulence) had greater chromosome copy numbers than the isolate maintained in culture for longer, suggesting that virulence attenuation may be associated with loss of chromosome copies. Our results suggest that genomic processes proposed as mechanisms for rapid evolution in Bd are correlated with virulence attenuation in laboratory culture within a single lineage of Bd. Moreover, these genomic processes can occur over extremely short timescales. On a practical level, our results underscore the importance of immediately cryo-archiving new Bd isolates and using fresh isolates, rather than samples cultured in the laboratory for long periods, for laboratory infection experiments. Finally, when attempting to predict disease outcomes for this ecologically important pathogen, it is critical to consider existing variation in virulence among isolates and the potential for shifts in virulence over short timescales.

Complex history of the amphibian-killing chytrid fungus revealed with genome resequencing data.

Understanding the evolutionary history of microbial pathogens is critical for mitigating the impacts of emerging infectious diseases on economically and ecologically important host species. We used a genome resequencing approach to resolve the evolutionary history of an important microbial pathogen, the chytrid Batrachochytrium dendrobatidis (Bd), which has been implicated in amphibian declines worldwide. We sequenced the genomes of 29 isolates of Bd from around the world, with an emphasis on North, Central, and South America because of the devastating effect that Bd has had on amphibian populations in the New World. We found a substantial amount of evolutionary complexity in Bd with deep phylogenetic diversity that predates observed global amphibian declines. By investigating the entire genome, we found that even the most recently evolved Bd clade (termed the global panzootic lineage) contained more genetic variation than previously reported. We also found dramatic differences among isolates and among genomic regions in chromosomal copy number and patterns of heterozygosity, suggesting complex and heterogeneous genome dynamics. Finally, we report evidence for selection acting on the Bd genome, supporting the hypothesis that protease genes are important in evolutionary transitions in this group. Bd is considered an emerging pathogen because of its recent effects on amphibians, but our data indicate that it has a complex evolutionary history that predates recent disease outbreaks. Therefore, it is important to consider the contemporary effects of Bd in a broader evolutionary context and identify specific mechanisms that may have led to shifts in virulence in this system.

Substrate-specific gene expression in Batrachochytrium dendrobatidis, the chytrid pathogen of amphibians.

Determining the mechanisms of host-pathogen interaction is critical for understanding and mitigating infectious disease. Mechanisms of fungal pathogenicity are of particular interest given the recent outbreaks of fungal diseases in wildlife populations. Our study focuses on Batrachochytrium dendrobatidis (Bd), the chytrid pathogen responsible for amphibian declines around the world. Previous studies have hypothesized a role for several specific families of secreted proteases as pathogenicity factors in Bd, but the expression of these genes has only been evaluated in laboratory growth conditions. Here we conduct a genome-wide study of Bd gene expression under two different nutrient conditions. We compare Bd gene expression profiles in standard laboratory growth media and in pulverized host tissue (i.e., frog skin). A large proportion of genes in the Bd genome show increased expression when grown in host tissue, indicating the importance of studying pathogens on host substrate. A number of gene classes show particularly high levels of expression in host tissue, including three families of secreted proteases (metallo-, serine- and aspartyl-proteases), adhesion genes, lipase-3 encoding genes, and a group of phylogenetically unusual crinkler-like effectors. We discuss the roles of these different genes as putative pathogenicity factors and discuss what they can teach us about Bd's metabolic targets, host invasion, and pathogenesis.

Only skin deep: shared genetic response to the deadly chytrid fungus in susceptible frog species.

Amphibian populations around the world are threatened by an emerging infectious pathogen, the chytrid fungus Batrachochytrium dendrobatidis (Bd). How can a fungal skin infection kill such a broad range of amphibian hosts? And do different host species have a similar response to Bd infection? Here, we use a genomics approach to understand the genetic response of multiple susceptible frog species to Bd infection. We characterize the transcriptomes of two closely related endangered frog species (Rana muscosa and Rana sierrae) and analyse whole genome expression profiles from frogs in controlled Bd infection experiments. We integrate the Rana results with a comparable data set from a more distantly related susceptible species (Silurana tropicalis). We demonstrate that Bd-infected frogs show massive disruption of skin function and show no evidence of a robust immune response. The genetic response to infection is shared across the focal susceptible species, suggesting a common effect of Bd on susceptible frogs.

A molecular perspective: biology of the emerging pathogen Batrachochytrium dendrobatidis.

Ten years after the first discovery of the chytrid pathogen Batrachochytrium dendrobatidis (Bd), the catastrophic effect of Bd on wild amphibian populations is indisputable. However, a number of persistent questions remain about Bd's origin and mechanisms of pathogenicity. Here we discuss the promise of genetic and genomic tools for answering these previously intractable questions about the biology and evolutionary history of Bd. Full genomes of 2 Bd strains have recently been sequenced, and Bd research on this species using population genetics, phylogenetics, proteomics, comparative genomics and functional genomics is already underway. We review some of the insights gleaned from the first studies using these genome-scale approaches focusing particularly on Bd's genomic architecture, patterns of global genetic variation, virulence factors and genetic interactions with hosts. Avenues of future research promise to be particularly fruitful and highlight the need for integrative studies that unite genetic, ecological and spatial data in both Bd and its amphibian hosts.

Genome-wide transcriptional response of Silurana (Xenopus) tropicalis to infection with the deadly chytrid fungus.

Emerging infectious diseases are of great concern for both wildlife and humans. Several highly virulent fungal pathogens have recently been discovered in natural populations, highlighting the need for a better understanding of fungal-vertebrate host-pathogen interactions. Because most fungal pathogens are not fatal in the absence of other predisposing conditions, host-pathogen dynamics for deadly fungal pathogens are of particular interest. The chytrid fungus Batrachochytrium dendrobatidis (hereafter Bd) infects hundreds of species of frogs in the wild. It is found worldwide and is a significant contributor to the current global amphibian decline. However, the mechanism by which Bd causes death in amphibians, and the response of the host to Bd infection, remain largely unknown. Here we use whole-genome microarrays to monitor the transcriptional responses to Bd infection in the model frog species, Silurana (Xenopus) tropicalis, which is susceptible to chytridiomycosis. To elucidate the immune response to Bd and evaluate the physiological effects of chytridiomycosis, we measured gene expression changes in several tissues (liver, skin, spleen) following exposure to Bd. We detected a strong transcriptional response for genes involved in physiological processes that can help explain some clinical symptoms of chytridiomycosis at the organismal level. However, we detected surprisingly little evidence of an immune response to Bd exposure, suggesting that this susceptible species may not be mounting efficient innate and adaptive immune responses against Bd. The weak immune response may be partially explained by the thermal conditions of the experiment, which were optimal for Bd growth. However, many immune genes exhibited decreased expression in Bd-exposed frogs compared to control frogs, suggesting a more complex effect of Bd on the immune system than simple temperature-mediated immune suppression. This study generates important baseline data for ongoing efforts to understand differences in response to Bd between susceptible and resistant frog species and the effects of chytridiomycosis in natural populations.

Maternal transfer of antibodies to eggs in Xenopus laevis.

The immune system of the African clawed frog, Xenopus laevis, includes nearly the full repertoire of lymphoid organs and immune cell types found in mammals. In contrast to the mammalian immune system, the development of the amphibian immune system occurs in the open environment. Oviparity necessitates a rapid ontogeny of the immune system. X. laevis larvae become immunocompetent about 2 weeks after fertilization of the egg. During this 2-week window, larvae cannot mount an adaptive immune response to potential pathogens and presumably must depend on innate responses. In the present study, the possibility of maternal transfer of antibodies to eggs was examined. Adult female X. laevis were injected three times at weekly intervals with the hapten-carrier complex, trinitrophenylated bovine serum albumin (TNP-BSA). The sera of immunized frogs demonstrated antibody activity to BSA, TNP-BSA, and, importantly, trinitrophenylated ovalbumin (TNP-OVA) when examined by enzyme-linked immunosorbent assay (ELISA). Reactivity to TNP-OVA confirmed that antibodies were produced against TNP. The adult female frogs were induced to lay eggs by injection of human chorionic gonadotropin. Next, membrane-free extracts of the eggs were treated with protease inhibitors in order to prevent proteolysis of proteins found in the eggs. On analysis by ELISA, it was found that TNP-specific antibodies were present in the egg extracts. This demonstrated the transfer of antigen-specific antibodies from adult females to eggs in X. laevis.