Complete Genome Sequence Resource of Pectobacterium colocasium Strain F1-1 that Causes Soft Rot Disease of Taro.

Pectobacterium colocasium is a recently named, narrow-host-range phytopathogenic bacterium causing soft rot of taro (Colocasium esculenta). It is found on the Chinese mainland and the island of Taiwan. Taro is a domesticated crop with a long history of cultivation in Taiwan and the Pacific islands. However, not much was known about Pectobacterium spp. from taro, especially from the islands in the Pacific. Herein, we report a high-quality, completely annotated genome sequence of P. colosacium strain F1-1. The 4,816,345 bp genome, which was assembled with Illumina and Nanopore reads with 217× and 311× coverage, respectively, consists of one chromosome and no plasmid. This completely circularized genome will aid future studies in comparative genomics, evolution, and pathogenicity of P. colocasium. This genome resource will also be helpful for developing strategies to control P. colocasium in taro.[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.

Whole genome resequencing and complementation tests reveal candidate loci contributing to bacterial wilt (Ralstonia sp.) resistance in tomato.

Tomato (Solanum lycopersicum) is one of the most economically important vegetable crops worldwide. Bacterial wilt (BW), caused by the Ralstonia solanacearum species complex, has been reported as the second most important plant pathogenic bacteria worldwide, and likely the most destructive. Extensive research has identified two major loci, Bwr-6 and Bwr-12, that contribute to resistance to BW in tomato; however, these loci do not completely explain resistance. Segregation of resistance in two populations that were homozygous dominant or heterozygous for all Bwr-6 and Bwr-12 associated molecular markers suggested the action of one or two resistance loci in addition to these two major QTLs. We utilized whole genome sequence data analysis and pairwise comparison of six BW resistant and nine BW susceptible tomato lines to identify candidate genes that, in addition to Bwr-6 and Bwr-12, contributed to resistance. Through this approach we found 27,046 SNPs and 5975 indels specific to the six resistant lines, affecting 385 genes. One sequence variant on chromosome 3 captured by marker Bwr3.2dCAPS located in the Asc (Solyc03g114600.4.1) gene had significant association with resistance, but it did not completely explain the resistance phenotype. The SNP associated with Bwr3.2dCAPS was located within the resistance gene Asc which was inside the previously identified Bwr-3 locus. This study provides a foundation for further investigations into new loci distributed throughout the tomato genome that could contribute to BW resistance and into the role of resistance genes that may act against multiple pathogens.

DNA-free CRISPR-Cas9 gene editing of wild tetraploid tomato Solanum peruvianum using protoplast regeneration.

Wild tomatoes (Solanum peruvianum) are important genomic resources for tomato research and breeding. Development of a foreign DNA-free clustered regularly interspaced short palindromic repeat (CRISPR)-Cas delivery system has potential to mitigate public concern about genetically modified organisms. Here, we established a DNA-free CRISPR-Cas9 genome editing system based on an optimized protoplast regeneration protocol of S. peruvianum, an important resource for tomato introgression breeding. We generated mutants for genes involved in small interfering RNAs biogenesis, RNA-DEPENDENT RNA POLYMERASE 6 (SpRDR6), and SUPPRESSOR OF GENE SILENCING 3 (SpSGS3); pathogen-related peptide precursors, PATHOGENESIS-RELATED PROTEIN-1 (SpPR-1) and PROSYSTEMIN (SpProSys); and fungal resistance (MILDEW RESISTANT LOCUS O, SpMlo1) using diploid or tetraploid protoplasts derived from in vitro-grown shoots. The ploidy level of these regenerants was not affected by PEG-Ca2+-mediated transfection, CRISPR reagents, or the target genes. By karyotyping and whole genome sequencing analysis, we confirmed that CRISPR-Cas9 editing did not introduce chromosomal changes or unintended genome editing sites. All mutated genes in both diploid and tetraploid regenerants were heritable in the next generation. spsgs3 null T0 regenerants and sprdr6 null T1 progeny had wiry, sterile phenotypes in both diploid and tetraploid lines. The sterility of the spsgs3 null mutant was partially rescued, and fruits were obtained by grafting to wild-type (WT) stock and pollination with WT pollen. The resulting seeds contained the mutated alleles. Tomato yellow leaf curl virus proliferated at higher levels in spsgs3 and sprdr6 mutants than in the WT. Therefore, this protoplast regeneration technique should greatly facilitate tomato polyploidization and enable the use of CRISPR-Cas for S. peruvianum domestication and tomato breeding.

Efficient and Economical Targeted Insertion in Plant Genomes via Protoplast Regeneration.

Versatile genome editing can be facilitated by the insertion of DNA sequences into specific locations. Current protocols involving CRISPR and Cas proteins rely on low efficiency homology-directed repair or non-homologous end joining with modified double-stranded DNA oligonucleotides as donors. Our simple protocol eliminates the need for expensive equipment, chemical and enzymatic donor DNA modification, or plasmid construction by using polyethylene glycol-calcium to deliver non-modified single-stranded DNA oligonucleotides and CRISPR-Cas9 ribonucleoprotein into protoplasts. Plants regenerated via edited protoplasts achieved targeted insertion frequencies of up to 50% in Nicotiana benthamiana and 13.6% in rapid cycling Brassica oleracea without antibiotic selection. Using a 60 nt donor containing 27 nt in each homologous arm, 6/22 regenerated N. benthamiana plants showed targeted insertions, and one contained a precise insertion of a 6 bp HindIII site. The inserted sequences were transmitted to the next generation and invite the possibility of future exploration of versatile genome editing by targeted DNA insertion in plants.

Comparative Genomic Analysis Uncovered Evolution of Pathogenicity Factors, Horizontal Gene Transfer Events, and Heavy Metal Resistance Traits in Citrus Canker Bacterium Xanthomonas citri subsp. citri.

Background: Worldwide citrus production is severely threatened by Asiatic citrus canker which is caused by the proteobacterium Xanthomonas citri subsp. citri. Foliar sprays of copper-based bactericides are frequently used to control plant bacterial diseases. Despite the sequencing of many X. citri strains, the genome diversity and distribution of genes responsible for metal resistance in X. citri subsp. citri strains from orchards with different management practices in Taiwan are not well understood. Results: The genomes of three X. citri subsp. citri strains including one copper-resistant strain collected from farms with different management regimes in Taiwan were sequenced by Illumina and Nanopore sequencing and assembled into complete circular chromosomes and plasmids. CRISPR spoligotyping and phylogenomic analysis indicated that the three strains were located in the same phylogenetic lineages and shared ∼3,000 core-genes with published X. citri subsp. citri strains. These strains differed mainly in the CRISPR repeats and pathogenicity-related plasmid-borne transcription activator-like effector (TALE)-encoding pthA genes. The copper-resistant strain has a unique, large copper resistance plasmid due to an unusual ∼40 kbp inverted repeat. Each repeat contains a complete set of the gene cluster responsible for copper and heavy metal resistance. Conversely, the copper sensitive strains carry no metal resistance genes in the plasmid. Through comparative analysis, the origin and evolution of the metal resistance clusters was resolved. Conclusion: Chromosomes remained constant among three strains collected in Taiwan, but plasmids likely played an important role in maintaining pathogenicity and developing bacterial fitness in the field. The evolution of pathogenicity factors and horizontal gene transfer events were observed in the three strains. These data suggest that agricultural management practices could be a potential trigger for the evolution of citrus canker pathogens. The decrease in the number of CRISPR repeats and pthA genes might be the result of adaptation to a less stressful environment. The metal resistance genes in the copper resistant X. citri strain likely originated from the Mauritian strain not the local copper-resistant X. euvesicatoria strain. This study highlights the importance of plasmids as 'vehicles' for exchanging genetic elements between plant pathogenic bacteria and contributing to bacterial adaptation to the environment.

Author Correction: The Apostasia genome and the evolution of orchids.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Small RNA profiling in Pinus pinaster reveals the transcriptome of developing seeds and highlights differences between zygotic and somatic embryos.

Regulation of seed development by small non-coding RNAs (sRNAs) is an important mechanism controlling a crucial phase of the life cycle of seed plants. In this work, sRNAs from seed tissues (zygotic embryos and megagametophytes) and from somatic embryos of Pinus pinaster were analysed to identify putative regulators of seed/embryo development in conifers. In total, sixteen sRNA libraries covering several developmental stages were sequenced. We show that embryos and megagametophytes express a large population of 21-nt sRNAs and that substantial amounts of 24-nt sRNAs were also detected, especially in somatic embryos. A total of 215 conserved miRNAs, one third of which are conifer-specific, and 212 high-confidence novel miRNAs were annotated. MIR159, MIR171 and MIR394 families were found in embryos, but were greatly reduced in megagametophytes. Other families, like MIR397 and MIR408, predominated in somatic embryos and megagametophytes, suggesting their expression in somatic embryos is associated with in vitro conditions. Analysis of the predicted miRNA targets suggests that miRNA functions are relevant in several processes including transporter activity at the cotyledon-forming stage, and sulfur metabolism across several developmental stages. An important resource for studying conifer embryogenesis is made available here, which may also provide insightful clues for improving clonal propagation via somatic embryogenesis.

Reference genome and comparative genome analysis for the WHO reference strain for Mycobacterium bovis BCG Danish, the present tuberculosis vaccine.

BACKGROUND: Mycobacterium bovis bacillus Calmette-Guérin (M. bovis BCG) is the only vaccine available against tuberculosis (TB). In an effort to standardize the vaccine production, three substrains, i.e. BCG Danish 1331, Tokyo 172-1 and Russia BCG-1 were established as the WHO reference strains. Both for BCG Tokyo 172-1 as Russia BCG-1, reference genomes exist, not for BCG Danish. In this study, we set out to determine the completely assembled genome sequence for BCG Danish and to establish a workflow for genome characterization of engineering-derived vaccine candidate strains. RESULTS: By combining second (Illumina) and third (PacBio) generation sequencing in an integrated genome analysis workflow for BCG, we could construct the completely assembled genome sequence of BCG Danish 1331 (07/270) (and an engineered derivative that is studied as an improved vaccine candidate, a SapM KO), including the resolution of the analytically challenging long duplication regions. We report the presence of a DU1-like duplication in BCG Danish 1331, while this tandem duplication was previously thought to be exclusively restricted to BCG Pasteur. Furthermore, comparative genome analyses of publicly available data for BCG substrains showed the absence of a DU1 in certain BCG Pasteur substrains and the presence of a DU1-like duplication in some BCG China substrains. By integrating publicly available data, we provide an update to the genome features of the commonly used BCG strains. CONCLUSIONS: We demonstrate how this analysis workflow enables the resolution of genome duplications and of the genome of engineered derivatives of the BCG Danish vaccine strain. The BCG Danish WHO reference genome will serve as a reference for future engineered strains and the established workflow can be used to enhance BCG vaccine standardization.

Draft Genome Assembly of the False Spider Mite Brevipalpus yothersi.

The false spider mite Brevipalpus yothersi infests a broad host plant range and has become one of the most economically important species within the genus Brevipalpus. This phytophagous mite inflicts damage by both feeding on plants and transmitting plant viruses. Here, we report the first draft genome sequence of the false spider mite, which is also the first plant virus mite vector to be sequenced. The ∼72 Mb genome (sequenced at 42× coverage) encodes ∼16,000 predicted protein-coding genes.

Functional and evolutionary genomic inferences in Populus through genome and population sequencing of American and European aspen.

The Populus genus is one of the major plant model systems, but genomic resources have thus far primarily been available for poplar species, and primarily Populus trichocarpa (Torr. & Gray), which was the first tree with a whole-genome assembly. To further advance evolutionary and functional genomic analyses in Populus, we produced genome assemblies and population genetics resources of two aspen species, Populus tremula L. and Populus tremuloides Michx. The two aspen species have distributions spanning the Northern Hemisphere, where they are keystone species supporting a wide variety of dependent communities and produce a diverse array of secondary metabolites. Our analyses show that the two aspens share a similar genome structure and a highly conserved gene content with P. trichocarpa but display substantially higher levels of heterozygosity. Based on population resequencing data, we observed widespread positive and negative selection acting on both coding and noncoding regions. Furthermore, patterns of genetic diversity and molecular evolution in aspen are influenced by a number of features, such as expression level, coexpression network connectivity, and regulatory variation. To maximize the community utility of these resources, we have integrated all presented data within the PopGenIE web resource (PopGenIE.org).

IMA Genome-F 8: Draft genome of Cercospora zeina, Fusarium pininemorale, Hawksworthiomyces lignivorus, Huntiella decipiens and Ophiostoma ips.

The genomes of Cercospora zeina, Fusarium pininemorale, Hawksworthiomyces lignivorus, Huntiella decipiens, and Ophiostoma ips are presented in this genome announcement. Three of these genomes are from plant pathogens and otherwise economically important fungal species. Fusarium pininemorale and H. decipiens are not known to cause significant disease but are closely related to species of economic importance. The genome sizes range from 25.99 Mb in the case of O. ips to 4.82 Mb for H. lignivorus. These genomes include the first reports of a genome from the genus Hawksworthiomyces. The availability of these genome data will allow the resolution of longstanding questions regarding the taxonomy of these species. In addition these genome sequences through comparative studies with closely related organisms will increase our understanding of how these species or close relatives cause disease.

Complete Genome Sequence of Bacillus cereus C1L, a Plant Growth-Promoting Rhizobacterium from the Rhizosphere of Formosa Lily in Taiwan.

Bacillus cereus C1L, a plant growth-promoting rhizobacterium, provides protection against fungal pathogens in monocot plants. To gain new insights into the biocontrol mechanisms used by this rhizobacterium, we determined the complete genome sequence of B. cereus C1L. One chromosome and three plasmids were identified with a total size of ~6.0 Mb.

The Apostasia genome and the evolution of orchids.

Constituting approximately 10% of flowering plant species, orchids (Orchidaceae) display unique flower morphologies, possess an extraordinary diversity in lifestyle, and have successfully colonized almost every habitat on Earth. Here we report the draft genome sequence of Apostasia shenzhenica, a representative of one of two genera that form a sister lineage to the rest of the Orchidaceae, providing a reference for inferring the genome content and structure of the most recent common ancestor of all extant orchids and improving our understanding of their origins and evolution. In addition, we present transcriptome data for representatives of Vanilloideae, Cypripedioideae and Orchidoideae, and novel third-generation genome data for two species of Epidendroideae, covering all five orchid subfamilies. A. shenzhenica shows clear evidence of a whole-genome duplication, which is shared by all orchids and occurred shortly before their divergence. Comparisons between A. shenzhenica and other orchids and angiosperms also permitted the reconstruction of an ancestral orchid gene toolkit. We identify new gene families, gene family expansions and contractions, and changes within MADS-box gene classes, which control a diverse suite of developmental processes, during orchid evolution. This study sheds new light on the genetic mechanisms underpinning key orchid innovations, including the development of the labellum and gynostemium, pollinia, and seeds without endosperm, as well as the evolution of epiphytism; reveals relationships between the Orchidaceae subfamilies; and helps clarify the evolutionary history of orchids within the angiosperms.

Systems genetics reveals a transcriptional network associated with susceptibility in the maize-grey leaf spot pathosystem.

We used a systems genetics approach to elucidate the molecular mechanisms of the responses of maize to grey leaf spot (GLS) disease caused by Cercospora zeina, a threat to maize production globally. Expression analysis of earleaf samples in a subtropical maize recombinant inbred line population (CML444 × SC Malawi) subjected in the field to C. zeina infection allowed detection of 20 206 expression quantitative trait loci (eQTLs). Four trans-eQTL hotspots coincided with GLS disease QTLs mapped in the same field experiment. Co-expression network analysis identified three expression modules correlated with GLS disease scores. The module (GY-s) most highly correlated with susceptibility (r = 0.71; 179 genes) was enriched for the glyoxylate pathway, lipid metabolism, diterpenoid biosynthesis and responses to pathogen molecules such as chitin. The GY-s module was enriched for genes with trans-eQTLs in hotspots on chromosomes 9 and 10, which also coincided with phenotypic QTLs for susceptibility to GLS. This transcriptional network has significant overlap with the GLS susceptibility response of maize line B73, and may reflect pathogen manipulation for nutrient acquisition and/or unsuccessful defence responses, such as kauralexin production by the diterpenoid biosynthesis pathway. The co-expression module that correlated best with resistance (TQ-r; 1498 genes) was enriched for genes with trans-eQTLs in hotspots coinciding with GLS resistance QTLs on chromosome 9. Jasmonate responses were implicated in resistance to GLS through co-expression of COI1 and enrichment of genes with the Gene Ontology term 'cullin-RING ubiquitin ligase complex' in the TQ-r module. Consistent with this, JAZ repressor expression was highly correlated with the severity of GLS disease in the GY-s susceptibility network.

The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea.

Seagrasses colonized the sea on at least three independent occasions to form the basis of one of the most productive and widespread coastal ecosystems on the planet. Here we report the genome of Zostera marina (L.), the first, to our knowledge, marine angiosperm to be fully sequenced. This reveals unique insights into the genomic losses and gains involved in achieving the structural and physiological adaptations required for its marine lifestyle, arguably the most severe habitat shift ever accomplished by flowering plants. Key angiosperm innovations that were lost include the entire repertoire of stomatal genes, genes involved in the synthesis of terpenoids and ethylene signalling, and genes for ultraviolet protection and phytochromes for far-red sensing. Seagrasses have also regained functions enabling them to adjust to full salinity. Their cell walls contain all of the polysaccharides typical of land plants, but also contain polyanionic, low-methylated pectins and sulfated galactans, a feature shared with the cell walls of all macroalgae and that is important for ion homoeostasis, nutrient uptake and O2/CO2 exchange through leaf epidermal cells. The Z. marina genome resource will markedly advance a wide range of functional ecological studies from adaptation of marine ecosystems under climate warming, to unravelling the mechanisms of osmoregulation under high salinities that may further inform our understanding of the evolution of salt tolerance in crop plants.

Comparative in silico analysis of SSRs in coding regions of high confidence predicted genes in Norway spruce (Picea abies) and Loblolly pine (Pinus taeda).

BACKGROUND: Microsatellites or simple sequence repeats (SSRs) are DNA sequences consisting of 1-6 bp tandem repeat motifs present in the genome. SSRs are considered to be one of the most powerful tools in genetic studies. We carried out a comparative study of perfect SSR loci belonging to class I (≥20) and class II (≥12 and <20 bp) types located in coding regions of high confidence genes in Picea abies and Pinus taeda. SSRLocator was used to retrieve SSRs from the full length CDS of predicted genes in both species. RESULTS: Trimers were the most abundant motifs in class I followed by hexamers in Picea abies, while trimers and hexamers were equally abundant in Pinus taeda class I SSRs. Hexamers were most frequent within class II SSRs followed by trimers, in both species. Although the frequency of genes containing SSRs was slightly higher in Pinus taeda, SSR counts per Mbp for class I was similar in both species (P-value = 0.22); while for class II SSRs, it was significantly higher in Picea abies (P-value = 0.00009). AT-rich motifs were higher in abundance than the GC-rich motifs, within class II SSRs in both the species (P-values = 10(-9) and 0). With reference to class I SSRs, AT-rich and GC-rich motifs were detected with equal frequency in Pinus taeda (P-value = 0.24); while in Picea abies, GC-rich motifs were detected with higher frequency than the AT-rich motifs (P-value = 0.0005). CONCLUSIONS: Our study gives a comparative overview of the genome SSRs composition based on high confidence genes in the two recently sequenced and economically important conifers and, also provides information on functional molecular markers that can be applied in genetic studies in Pinus and Picea species.

Serendipitous Meta-Transcriptomics: The Fungal Community of Norway Spruce (Picea abies).

After performing de novo transcript assembly of >1 billion RNA-Sequencing reads obtained from 22 samples of different Norway spruce (Picea abies) tissues that were not surface sterilized, we found that assembled sequences captured a mix of plant, lichen, and fungal transcripts. The latter were likely expressed by endophytic and epiphytic symbionts, indicating that these organisms were present, alive, and metabolically active. Here, we show that these serendipitously sequenced transcripts need not be considered merely as contamination, as is common, but that they provide insight into the plant's phyllosphere. Notably, we could classify these transcripts as originating predominantly from Dothideomycetes and Leotiomycetes species, with functional annotation of gene families indicating active growth and metabolism, with particular regards to glucose intake and processing, as well as gene regulation.

The Plant Genome Integrative Explorer Resource: PlantGenIE.org.

Accessing and exploring large-scale genomics data sets remains a significant challenge to researchers without specialist bioinformatics training. We present the integrated PlantGenIE.org platform for exploration of Populus, conifer and Arabidopsis genomics data, which includes expression networks and associated visualization tools. Standard features of a model organism database are provided, including genome browsers, gene list annotation, Blast homology searches and gene information pages. Community annotation updating is supported via integration of WebApollo. We have produced an RNA-sequencing (RNA-Seq) expression atlas for Populus tremula and have integrated these data within the expression tools. An updated version of the ComPlEx resource for performing comparative plant expression analyses of gene coexpression network conservation between species has also been integrated. The PlantGenIE.org platform provides intuitive access to large-scale and genome-wide genomics data from model forest tree species, facilitating both community contributions to annotation improvement and tools supporting use of the included data resources to inform biological insight.