Development of genomic simple sequence repeat markers for Glycyrrhiza lepidota and cross-amplification of other Glycyrrhiza species

Background

Licorice (Glycyrrhiza spp. L.) is used as a natural sweetener and medicinal herb in European and Asian countries. Molecular studies have been conducted to find differences between wild and cultivated species because most wild species are highly resistant to abiotic and biotic stresses compared with their cultivated species. However, few molecular markers have been developed for studying the genetic diversity and population structure of licorice species and to identify differences between cultivars. Thus, the present study aimed to develop a set of genomic simple sequence repeat (SSR) markers for molecular studies of these species.

Methods

In the present study, we developed polymorphic SSR markers based on whole-genomesequence data of Glycyrrhiza lepidota. Then, based on the sequence information, the polymorphic SSR markers were developed. The SSR markers were applied to 23 Glycyrrhiza individual plants. We also evaluated the phylogenetic relationships and interspecies transferability among samples.

Results

The genetic diversity analysis using these markers identified 2–23 alleles, and the major allele frequency, observed heterozygosity, genetic diversity, and polymorphism information content were 0.11–0.91, 0–0.90, 0.17–0.94, and 0.15–0.93, respectively. Interspecies transferability values were 93.5%, 91.6%, and 91.1% for G. echinata, G. glabra, and G. uralensis, respectively. Phylogenetic analysis clustered cultivated (group 1) and wild (group 2) species into three and two subgroups, respectively. The reported markers represent a valuable resource for the genetic characteri z ation of Glycyrrhiza spp. for theanalysis of its genetic variability, and as a tool for licorice transferability. This is the first intraspecific study in a collection of Glycyrrhiza spp. germplasm using SSR markers.

Preoperative N stage evaluation in advanced gastric cancer patients using multidetector CT: can the sum of the diameters of metastatic LNs be used for N stage evaluation?

AIM

To compare the diagnostic performance of total counts of metastatic lymph nodes (LN-sum) and conventional multidetector (MD) computed tomography (CT) staging in the nodal evaluation of advanced gastric cancer (AGC) patients.

MATERIALS AND METHODS

In total, 127 consecutive patients who underwent preoperative MDCT and gastrectomy for AGC were identified. Metastatic LNs on MDCT were defined as LNs with a short axis ≥8 mm, marked or heterogeneous enhancement, and morphological features (central necrosis, round shape, clustering). The sensitivity, specificity, accuracy, positive predictive value (PPV), and negative predictive value (NPV) of the N-stage using LN-sum and conventional MDCT staging were generated and compared. In addition, metastatic LN counts between the MDCT and the histopathological examinations and correlation between LN-sum and histopathological nodal status were analysed.

RESULTS

The total counts of metastatic LNs on MDCT was significantly smaller than those detected in histopathological assessments (p<0.0001). LN-sum showed significant correlation with the pathological N stage and the number of metastatic LNs (rho=0.69, 0.73, p<0.0001). The areas under the receiver operating characteristic curve were 0.896, and 0.835, for N stage ≥N2 and N3, with cut-off values of 12.5 and 23.5 mm, respectively. LN-sum provided better diagnostic performance than conventional MDCT staging for discriminating N0-2 versus N3; sensitivity, accuracy, PPV and NPV of LN-sum were significantly higher (80.4 versus 52.2%, 81.1 versus 68.5%, 71.2 versus 57.1%, and 88 versus 74.1%).

CONCLUSION

LN-sum may be sufficiently useful in assessing the N3 stage of AGC and may help to plan appropriate therapy for AGC patients.

ABSCISIC ACID-INSENSITIVE 3 is involved in brassinosteroid-mediated regulation of flowering in plants

ABSCISIC ACID-INSENSITIVE 3 (ABI3) is one of the essential transcription factors of ABSCISIC ACID (ABA) signaling, functioning in seed germination, early seedling development, and abiotic stress tolerance. A recent study showed that epigenetic repression of ABI3 by brassinosteroid (BR)-activated BRI1 EMS SUPPRESSOR1 (BES1)-TOPLESS (TPL)HISTONE DEACETYLASE 19 (HDA19) repressor complex is a critical event for promoting seed germination and early seedling development. However, other physiological roles of the repression of ABI3 and ABA responses by BES1-mediated BR signaling pathways remain elusive. Here, we show that BES1-mediated suppression of ABI3 promotes floral transition and ABI3 acts as a negative regulator for flowering. Ectopic expression of ABI3 specifically compromised the early flowering phenotype of bes1-D and induced severe late-flowering phenotypes in wild-type Arabidopsis and Solanum lycopersicum plants. Both spatiotemporal expression patterns and global transcriptome analysis of ABI3-overexpressing plants supported the biological roles of ABI3 in the negative regulation of floral transition and reproduction. Finally, we confirmed that the loss of function of ABI3 induced early-flowering phenotypes in both long- and short-day conditions. In conclusion, our data suggest that BES1-mediated regulation of ABI3 is important in the reproductive phase transition of plants.

Brassinosteroids facilitate xylem differentiation and wood formation in tomato

BR signaling pathways facilitate xylem differentiation and wood formation by fine tuning SlBZR1/SlBZR2-mediated gene expression networks involved in plant secondary growth. Brassinosteroid (BR) signaling and BR crosstalk with diverse signaling cues are involved in the pleiotropic regulation of plant growth and development. Recent studies reported the critical roles of BR biosynthesis and signaling in vascular bundle development and plant secondary growth; however, the molecular bases of these roles are unclear. Here, we performed comparative physiological and anatomical analyses of shoot morphological growth in a cultivated wild-type tomato (Solanum lycopersicum cv. BGA) and a BR biosynthetic mutant [Micro Tom (MT)]. We observed that the canonical BR signaling pathway was essential for xylem differentiation and sequential wood formation by facilitating plant secondary growth. The gradual retardation of xylem development phenotypes during shoot vegetative growth in the BR-deficient MT tomato mutant recovered completely in response to exogenous BR treatment or genetic complementation of the BR biosynthetic DWARF (D) gene. By contrast, overexpression of the tomato Glycogen synthase kinase 3 (SlGSK3) or CRISPR-Cas9 (CR)-mediated knockout of the tomato Brassinosteroid-insensitive 1 (SlBRI1) impaired BR signaling and resulted in severely defective xylem differentiation and secondary growth. Genetic modulation of the transcriptional activity of the tomato Brassinazole-resistant 1/2 (SlBZR1/SlBZR2) confirmed the positive roles of BR signaling pathways for xylem differentiation and secondary growth. Our data indicate that BR signaling pathways directly promote xylem differentiation and wood formation by canonical BR-activated SlBZR1/SlBZR2.

Aptamer-Based Pathogen Monitoring for Salmonella enterica ser. Typhimurium

Salmonella enterica ser. Typhimurium is a foodborne pathogen that causes salmonellosis. Symptoms of salmonellosis include fever, diarrhea, and gastroenteritis. Conventional culture methods for detecting foodborne bacterial pathogens require long incubation time, expensive immunoassay methods, and sample enrichment steps. The objective of this study was to develop an aptamer-based sandwich assay to detect S. enterica ser. Typhimurium. First, S. enterica ser. Typhimurium specific binding aptamers (S11 and S24) were obtained by whole-cell SELEX for high sensitivity and specificity detection using live S. enterica ser. Typhimurium. S11 and S24 aptamers were able to capture S. enterica ser. Typhimurium selectively and distinguish it from other species of Salmonella (S. enterica ser. Typhimurium, S. enterica ser. Choleraesuis, S. enterica ser. Dublin, and S. enterica ser. Enteritidis) and food-borne bacterial pathogens (Escherichia coli K12, Listeria monocytogenes, Shigella sonnei and Staphylococcus aureus) with KD values of 4.41×10-12 M and 3.75×10-11 M, respectively. S. enterica ser. Typhimurium. aptamer-based sandwich assay exhibited a linear response for sensing S. enterica ser. Typhimurium. cells at concentration ranging from 2×101 to 2×105 CFU/mL. This aptamer-based sandwich assay can be used for on-site detection of S. enterica ser. Typhimurium.

Comparative transcriptome analysis identified candidate genes involved in mycelium browning in Lentinula edodes

Background

Lentinula edodes is one of the most popular edible mushroom species in the world and contains useful medicinal components, such as lentinan. The light-induced formation of brown film on the vegetative mycelial tissues of L. edodes is an important process for ensuring the quantity and quality of this edible mushroom. To understand the molecular mechanisms underlying this critical developmental process in L. edodes, we characterized the morphological phenotypic changes in a strain, Chamaram, associated with abnormal brown film formation and compared its genome-wide transcriptional features.

Results

In the present study, we performed genome-wide transcriptome analyses of different vegetative mycelium growth phenotypes, namely, early white, normal brown, and defective dark yellow partial brown films phenotypes which were exposed to different light conditions. The analysis revealed the identification of clusters of genes specific to the light-induced brown film phenotypes. These genes were significantly associated with light sensing via photoreceptors such as FMN- and FAD-bindings, signal transduction by kinases and GPCRs, melanogenesis via activation of tyrosinases, and cell wall degradation by glucanases, chitinases, and laccases, which suggests these processes are involved in the formation of mycelial browning in L. edodes. Interestingly, hydrophobin genes such as SC1 and SC3 exhibited divergent expression levels in the normal and abnormal brown mycelial films, indicating the ability of these genes to act in fruiting body initiation and formation of dikaryotic mycelia. Furthermore, we identified the up-regulation of glycoside hydrolase domain-containing genes in the normal brown film but not in the abnormal film phenotype, suggesting that cell wall degradation in the normal brown film phenotype is crucial in the developmental processes related to the initiation and formation of fruiting bodies.

Conclusions

This study systematically analysed the expression patterns of light-induced browning-related genes in L. edodes. Our findings provide information for further investigations of browning formation mechanisms in L. edodes and a foundation for future L. edodes breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5509-4) contains supplementary material, which is available to authorized users.

Development of Polymorphic Simple Sequence Repeat Markers using High-Throughput Sequencing in Button Mushroom (Agaricus bisporus)

The white button mushroom (Agaricus bisporus) is one of the most widely cultivated species of edible mushroom. Despite its economic importance, relatively little is known about the genetic diversity of this species. Illumina paired-end sequencing produced 43,871,558 clean reads and 69,174 contigs were generated from five offspring. These contigs were subsequently assembled into 57,594 unigenes. The unigenes were annotated with reference genome in which 6,559 unigenes were associated with clusters, indicating orthologous genes. Gene ontology classification assigned many unigenes. Based on genome data of the five offspring, 44 polymorphic simple sequence repeat (SSR) markers were developed. The major allele frequency ranged from 0.42 to 0.92. The number of genotypes and the number of alleles ranged from 1 to 4, and from 2 to 4, respectively. The observed heterozygosity and the expected heterozygosity ranged from 0.00 to 1.00, and from 0.15 to 0.64, respectively. The polymorphic information content value ranged from 0.14 to 0.57. The genetic distances and UPGMA clustering discriminated offspring strains. The SSR markers developed in this study can be applied in polymorphism analyses of button mushroom and for cultivar discrimination.

The complete mitochondrial genome of the edible and phytopathogenic fungus Desarmillaria tabescens

Desarmillaria tabescens is one of the most important edible, medicinal, and phytopathogenic basidiomycetes. The complete mitochondrial genome of this species was determined using next-generation sequencing technology. This mitogenome is a circular molecule of 93,439 bp with a GC content of 29.28% and contains 15 protein-coding, two rRNA (rnl and rns), and 24 tRNA genes. Phylogenetic analysis revealed that D. tabescens is genetically closest to Agrocybe aegerita. Desarmillaria tabescens mitogenome can contribute to our understanding of the phylogeny and evolution of this species.

TaF: a web platform for taxonomic profile-based fungal gene prediction

INTRODUCTION

The accurate prediction and annotation of gene structures from the genome sequence of an organism enable genome-wide functional analyses to obtain insight into the biological properties of an organism.

OBJECTIVES

We recently developed a highly accurate filamentous fungal gene prediction pipeline and web platform called TaF. TaF is a homology-based gene predictor employing large-scale taxonomic profiling to search for close relatives in genome queries.

METHODS

TaF pipeline consists of four processing steps; (1) taxonomic profiling to search for close relatives to query, (2) generation of hints for determining exon-intron boundaries from orthologous protein sequence data of the profiled species, (3) gene prediction by combination of ab inito and evidence-based prediction methods, and (4) homology search for gene models.

RESULTS

TaF generates extrinsic evidence that suggests possible exon-intron boundaries based on orthologous protein sequence data, thus reducing false-positive predictions of gene structure based on distantly related orthologs data. In particular, the gene prediction method using taxonomic profiling shows very high accuracy, including high sensitivity and specificity for gene models, suggesting a new approach for homology-based gene prediction from newly sequenced or uncharacterized fungal genomes, with the potential to improve the quality of gene prediction.

CONCLUSION

TaF will be a useful tool for fungal genome-wide analyses, including the identification of targeted genes associated with a trait, transcriptome profiling, comparative genomics, and evolutionary analysis.

Translational control of phloem development by RNA G-quadruplex-JULGI determines plant sink strength

The emergence of a plant vascular system was a prerequisite for the colonization of land; however, it is unclear how the photosynthate transporting system was established during plant evolution. Here, we identify a novel translational regulatory module for phloem development involving the zinc-finger protein JULGI (JUL) and its targets, the 5' untranslated regions (UTRs) of the SUPPRESSOR OF MAX2 1-LIKE4/5 (SMXL4/5) mRNAs, which is exclusively conserved in vascular plants. JUL directly binds and induces an RNA G-quadruplex in the 5' UTR of SMXL4/5, which are key promoters of phloem differentiation. We show that RNA G-quadruplex formation suppresses SMXL4/5 translation and restricts phloem differentiation. In turn, JUL deficiency promotes phloem formation and strikingly increases sink strength per seed. We propose that the translational regulation by the JUL/5' UTR G-quadruplex module is a major determinant of phloem establishment, thereby determining carbon allocation to sink tissues, and that this mechanism was a key invention during the emergence of vascular plants.

Genome-wide transcriptomic analysis of BR-deficient Micro-Tom reveals correlations between drought stress tolerance and brassinosteroid signaling in tomato

Brassinosteroids (BRs) are plant steroid hormones that play crucial roles in a range of growth and developmental processes. Although BR signal transduction and biosynthetic pathways have been well characterized in model plants, their biological roles in an important crop, tomato (Solanum lycopersicum), remain unknown. Here, cultivated tomato (WT) and a BR synthesis mutant, Micro-Tom (MT), were compared using physiological and transcriptomic approaches. The cultivated tomato showed higher tolerance to drought and osmotic stresses than the MT tomato. However, BR-defective phenotypes of MT, including plant growth and stomatal closure defects, were completely recovered by application of exogenous BR or complementation with a SlDWARF gene. Using genome-wide transcriptome analysis, 619 significantly differentially expressed genes (DEGs) were identified between WT and MT plants. Several DEGs were linked to known signaling networks, including those related to biotic/abiotic stress responses, lignification, cell wall development, and hormone responses. Consistent with the higher susceptibility of MT to drought stress, several gene sets involved in responses to drought and osmotic stress were differentially regulated between the WT and MT tomato plants. Our data suggest that BR signaling pathways are involved in mediating the response to abiotic stress via fine-tuning of abiotic stress-related gene networks in tomato plants.

Diversity of A mating type in Lentinula edodes and mating type preference in the cultivated strains

Diversity of A mating type in Lentinula edodes has been assessed by analysis of A mating loci in 127 strains collected from East Asia. It was discovered that hypervariable sequence region with an approximate length of 1 kb in the A mating locus, spanning 5' region of HD2-intergenic region-5' region of HD1, could represent individual A mating type as evidenced by comprehensive mating analysis. The sequence analysis revealed 27 A mating type alleles from 96 cultivated strains and 48 alleles from 31 wild strains. Twelve of them commonly appeared, leaving 63 unique A mating type alleles. It was also revealed that only A few A mating type alleles such as A1, A4, A5, and A7 were prevalent in the cultivated strains, accounting for 62.5% of all A mating types. This implies preferred selection of certain A mating types in the process of strain development and suggests potential role of A mating genes in the expression of genes governing mushroom quality. Dominant expression of an A mating gene HD1 was observed from A1 mating locus, the most prevalent A allele, in A1-containing dikaryons. However, connections between HD1 expression and A1 preference in the cultivated strains remain to be verified. The A mating type was highly diverse in the wild strains. Thirty-six unique A alleles were discovered from relatively small and confined area of mountainous region in Korean peninsula. The number will further increase because no A allele has been recurrently observed in the wild strains and thus newly discovered strain will have good chances to contain new A allele. The high diversity in small area also suggests that the A mating locus has evolved rapidly and thus its diversity will further increase.

Long-read transcriptome data for improved gene prediction in Lentinula edodes

Lentinula edodes is one of the most popular edible mushrooms in the world and contains useful medicinal components such as lentinan. The whole-genome sequence of L. edodes has been determined with the objective of discovering candidate genes associated with agronomic traits, but experimental verification of gene models with correction of gene prediction errors is lacking. To improve the accuracy of gene prediction, we produced 12.6 Gb of long-read transcriptome data of variable lengths using PacBio single-molecule real-time (SMRT) sequencing and generated 36,946 transcript clusters with an average length of 2.2 kb. Evidence-driven gene prediction on the basis of long- and short-read RNA sequencing data was performed; a total of 16,610 protein-coding genes were predicted with error correction. Of the predicted genes, 42.2% were verified to be covered by full-length transcript clusters. The raw reads have been deposited in the NCBI SRA database under accession number PRJNA396788.

Isoform Sequencing Provides a More Comprehensive View of the Panax ginseng Transcriptome

Korean ginseng (Panax ginseng C.A. Meyer) has been widely used for medicinal purposes and contains potent plant secondary metabolites, including ginsenosides. To obtain transcriptomic data that offers a more comprehensive view of functional genomics in P. ginseng, we generated genome-wide transcriptome data from four different P. ginseng tissues using PacBio isoform sequencing (Iso-Seq) technology. A total of 135,317 assembled transcripts were generated with an average length of 3.2 kb and high assembly completeness. Of those unigenes, 67.5% were predicted to be complete full-length (FL) open reading frames (ORFs) and exhibited a high gene annotation rate. Furthermore, we successfully identified unique full-length genes involved in triterpenoid saponin synthesis and plant hormonal signaling pathways, including auxin and cytokinin. Studies on the functional genomics of P. ginseng seedlings have confirmed the rapid upregulation of negative feed-back loops by auxin and cytokinin signaling cues. The conserved evolutionary mechanisms in the auxin and cytokinin canonical signaling pathways of P. ginseng are more complex than those in Arabidopsis thaliana. Our analysis also revealed a more detailed view of transcriptome-wide alternative isoforms for 88 genes. Finally, transposable elements (TEs) were also identified, suggesting transcriptional activity of TEs in P. ginseng. In conclusion, our results suggest that long-read, full-length or partial-unigene data with high-quality assemblies are invaluable resources as transcriptomic references in P. ginseng and can be used for comparative analyses in closely related medicinal plants.

Characterization of the Rosellinia necatrix Transcriptome and Genes Related to Pathogenesis by Single-Molecule mRNA Sequencing

White root rot disease, caused by the pathogen Rosellinia necatrix, is one of the world's most devastating plant fungal diseases and affects several commercially important species of fruit trees and crops. Recent global outbreaks of R. necatrix and advances in molecular techniques have both increased interest in this pathogen. However, the lack of information regarding the genomic structure and transcriptome of R. necatrix has been a barrier to the progress of functional genomic research and the control of this harmful pathogen. Here, we identified 10,616 novel full-length transcripts from the filamentous hyphal tissue of R. necatrix (KACC 40445 strain) using PacBio single-molecule sequencing technology. After annotation of the unigene sets, we selected 14 cell cycle-related genes, which are likely either positively or negatively involved in hyphal growth by cell cycle control. The expression of the selected genes was further compared between two strains that displayed different growth rates on nutritional media. Furthermore, we predicted pathogen-related effector genes and cell wall-degrading enzymes from the annotated gene sets. These results provide the most comprehensive transcriptomal resources for R. necatrix, and could facilitate functional genomics and further analyses of this important phytopathogen.

Development and Molecular Characterization of Novel Polymorphic Genomic DNA SSR Markers in Lentinula edodes

Sixteen genomic DNA simple sequence repeat (SSR) markers of Lentinula edodes were developed from 205 SSR motifs present in 46.1-Mb long L. edodes genome sequences. The number of alleles ranged from 3-14 and the major allele frequency was distributed from 0.17-0.96. The values of observed and expected heterozygosity ranged from 0.00-0.76 and 0.07-0.90, respectively. The polymorphic information content value ranged from 0.07-0.89. A dendrogram, based on 16 SSR markers clustered by the paired hierarchical clustering' method, showed that 33 shiitake cultivars could be divided into three major groups and successfully identified. These SSR markers will contribute to the efficient breeding of this species by providing diversity in shiitake varieties. Furthermore, the genomic information covered by the markers can provide a valuable resource for genetic linkage map construction, molecular mapping, and marker-assisted selection in the shiitake mushroom.

Development of 44 Novel Polymorphic SSR Markers for Determination of Shiitake Mushroom (Lentinula edodes) Cultivars

The shiitake mushroom (Lentinulaedodes) is one of the most popular edible mushrooms in the world and has attracted attention for its value in medicinal and pharmacological uses. With recent advanced research and techniques, the agricultural cultivation of the shiitake mushroom has been greatly increased, especially in East Asia. Additionally, demand for the development of new cultivars with good agricultural traits has been greatly enhanced, but the development processes are complicated and more challenging than for other edible mushrooms. In this study, we developed 44 novel polymorphic simple sequence repeat (SSR) markers for the determination of shiitake mushroom cultivars based on a whole genome sequencing database of L. edodes. These markers were found to be polymorphic and reliable when screened in 23 shiitake mushroom cultivars. For the 44 SSR markers developed in this study, the major allele frequency ranged from 0.13 to 0.94; the number of genotypes and number of alleles were each 2-11; the observed and expected heterozygosity were 0.00-1.00 and 0.10-0.90, respectively; and the polymorphic information content value ranged from 0.10 to 0.89. These new markers can be used for molecular breeding, the determination of cultivars, and other applications.