Resveratrol analogues and metabolites selectively inhibit human and rat 11β-hydroxysteroid dehydrogenase 1 as the therapeutic drugs: structure-activity relationship and molecular dynamics analysis

Resveratrol is converted to various metabolites by gut microbiota. Human and rat liver 11β-hydroxysteroid dehydrogenase 1 (11β-HSD1) are critical for glucocorticoid activation, while 11β-HSD2 in the kidney does the opposite reaction. It is still uncertain whether resveratrol and its analogues selectively inhibit 11β-HSD1. In this study, the inhibitory strength, mode of action, structure-activity relationship (SAR), and docking analysis of resveratrol analogues on human, rat, and mouse 11β-HSD1 and 11β-HSD2 were performed. The inhibitory strength of these chemicals on human 11β-HSD1 was dihydropinosylvin (6.91 μM) > lunularin (45.44 μM) > pinostilbene (46.82 μM) > resveratrol (171.1 μM) > pinosylvin (193.8 μM) > others. The inhibitory strength of inhibiting rat 11β-HSD1 was pinostilbene (9.67 μM) > lunularin (17.39 μM) > dihydropinosylvin (19.83 μM) > dihydroresveratrol (23.07 μM) > dihydroxystilbene (27.84 μM) > others and dihydropinosylvin (85.09 μM) and pinostilbene (>100 μM) inhibited mouse 11β-HSD1. All chemicals did not inhibit human, rat, and mouse 11β-HSD2. It was found that dihydropinosylvin, lunularin, and pinostilbene were competitive inhibitors of human 11β-HSD1 and that pinostilbene, lunularin, dihydropinosylvin, dihydropinosylvin and dihydroxystilbene were mixed inhibitors of rat 11β-HSD1. Docking analysis showed that they bind to the steroid-binding site of human and rat 11β-HSD1. In conclusion, resveratrol and its analogues can selectively inhibit human and rat 11β-HSD1, and mouse 11β-HSD1 is insensitive to the inhibition of resveratrol analogues.

An interpretable knee replacement risk assessment system for osteoarthritis patients

Objective

Knee osteoarthritis (OA) is a complex disease with heterogeneous representations. Although it is modifiable to prevention and early treatment, there still lacks a reliable and accurate prognostic tool. Hence, we aim to develop a quantitative and self-administrable knee replacement (KR) risk stratification system for knee osteoarthritis (KOA) patients with clinical features.

Method

A total of 14 baseline features were extracted from 9592 cases in the Osteoarthritis Initiative (OAI) cohort. A survival model was constructed using the Random Survival Forests algorithm. The prediction performance was evaluated with the concordance index (C-index) and average receiver operating characteristic curve (AUC). A three-class KR risk stratification system was built to differentiate three distinct KR-free survival groups. Thereafter, Shapley Additive Explanations (SHAP) was introduced for model explanation.

Results

KR incidence was accurately predicted by the model with a C-index of 0.770 (±0.0215) and an average AUC of 0.807 (±0.0181) with 14 clinical features. Three distinct survival groups were observed from the ten-point KR risk stratification system with a four-year KR rate of 0.79%, 5.78%, and 16.2% from the low, medium, and high-risk groups respectively. KR is mainly caused by pain medication use, age, surgery history, diabetes, and a high body mass index, as revealed by SHAP.

Conclusion

A self-administrable and interpretable KR survival model was developed, underscoring a KR risk scoring system to stratify KOA patients. It will encourage regular self-assessments within the community and facilitate personalised healthcare for both primary and secondary prevention of KOA.

A long noncoding RNA functions in pumpkin fruit development through S-adenosyl-L-methionine synthetase

Long noncoding RNAs (lncRNAs) play important roles in various biological processes. However, the regulatory roles of lncRNAs underlying fruit development have not been extensively studied. The pumpkin (Cucurbita spp.) is a preferred model for understanding the molecular mechanisms regulating fruit development because of its variable shape and size and large inferior ovary. Here, we performed strand-specific transcriptome sequencing on pumpkin (Cucurbita maxima "Rimu") fruits at 6 developmental stages and identified 5,425 reliably expressed lncRNAs. Among the 332 lncRNAs that were differentially expressed during fruit development, the lncRNA MSTRG.44863.1 was identified as a negative regulator of pumpkin fruit development. MSTRG.44863.1 showed a relatively high expression level and an obvious period-specific expression pattern. Transient overexpression and silencing of MSTRG.44863.1 significantly increased and decreased the content of 1-aminocyclopropane carboxylic acid (a precursor of ethylene) and ethylene production, respectively. RNA pull-down and microscale thermophoresis assays further revealed that MSTRG.44863.1 can interact with S-adenosyl-L-methionine synthetase (SAMS), an enzyme in the ethylene synthesis pathway. Considering that ethylene negatively regulates fruit development, these results indicate that MSTRG.44863.1 plays an important role in the regulation of pumpkin fruit development, possibly through interacting with SAMS and affecting ethylene synthesis. Overall, our findings provide a rich resource for further study of fruit-related lncRNAs while offering insights into the regulation of fruit development in plants.

Mitochondria targeted biomimetic platform for chemo/photodynamic combination therapy against osteosarcoma

Clinical treatment for osteosarcoma (OS) is still lacking effective means, and no significant progress in OS treatment have been made in recent years. Single chemotherapy has serious side effects and can produce drug resistance easily, resulting poor therapeutic effect. As a modern and non-invasive treatment form, photodynamic therapy (PDT) is widely used to treat diverse cancers. Chemotherapy in combination with PDT is a particularly efficient antitumor method that could overcome the defects of monotherapies. Since mitochondria is a key subcellular organelle involved in cell apoptosis regulation, targeting tumor cells mitochondria for drug delivery has become an important entry point for anti-tumor therapy. Herein, we rationally designed a core-shell structured biomimetic nanoplatform, i.e., D@SLNP@OSM-IR780, to achieve tumor homologous targeting and mitochondria targeted drug release for chemotherapy combined with PDT against OS. Upon 808 nm laser irradiation, D@SLNP@OSM-IR780 exhibited excellent photo-cytotoxicity in vitro. The excellent targeting effect of D@SLNP@OSM-IR780 in tumor tissues produced a tumor inhibition rate of 98.9% in vivo. We further indicated that synergistic chemo-photodynamic effect induced by D@SLNP@OSM-IR780 could activate mitochondria-mediated apoptosis pathway, along with host immune response and potential photothermal effect. On the whole, D@SLNP@OSM-IR780 is revealed to be a promising platform for OS targeted combination therapeutics.

Development and validation of a nomogram and risk stratification system to predict overall survival after surgical repair for pediatric patients with medulloblastoma based on easily accessible variables

OBJECTIVE

This study aimed to develop and validate a nomogram and risk stratification system for the overall survival of pediatric patients with medulloblastoma after surgical repair.

PATIENTS AND METHODS

In this multicenter, retrospective study, consecutive patients who underwent surgery for medulloblastoma at Shanghai Children's Medical Center and the First Affiliated Hospital of Fujian Medical University from 2010 to 2022 formed the training and external validation datasets, respectively. Univariable and multivariable Cox regression analyses were performed to identify variables associated with mortality in the training dataset. A nomogram prediction model was developed based on independent variables in the multivariable Cox regression analysis to predict the 1-, 3-, and 5-year overall survival. The area under receiver operating characteristic curve (AUC) and calibration curve were used to evaluate the discrimination and calibration of the nomogram. A risk stratification system based on the median risk score was also established to divide patients into two risk groups.

RESULTS

In the training dataset, Cox regression analyses identified tumor size, brainstem involvement and chemotherapy as independent predictors for overall survival. The AUC of the nomogram was 0.75 at 1 year, 0. 75 at 3 years, 0.77 at 5 years in the training dataset, 0.74 at 1 year, 0.70 at 3 years, and 0.70 at 5 years in the validation dataset. The calibration curve for the probability of 1-, 3-, and 5-year survival showed good agreement between the nomogram prediction and actual observation in the training and validation datasets. The risk stratification system could perfectly classify patients into two risk groups, and the overall survival in the two groups had a good division.

CONCLUSIONS

This low-cost, convenient, and noninvasive nomogram can be translated into clinical practice as a tool for risk stratification and individualized prognosis prediction for children with medulloblastoma.

Biomimetic lipid nanoparticles for homologous-targeting and enhanced photodynamic therapy against glioma

Biomimetic nano-platforms have attracted extensive attention due to their good biocompatibility, low immunogenicity, and homologous targeting to lesions. In this study, glioma cell membranes are used to encapsulate indocyanine green (ICG) loaded nanoparticles (SLNP/ICG), termed as SLNP/ICG@M for targeted photodynamic therapy (PDT) against glioma. Cell membrane modification significantly enhances cellular uptake of SLNP/ICG@M in homologous glioma cells in vitro and tumor distribution in vivo. Furthermore, SLNP/ICG@M can stimulate glioma cells to generate plentiful reactive oxygen species (ROS) under NIR irradiation, finally producing excellent photo-cytotoxicity and the optimal tumor growth inhibition with a tumor suppression rate of 93.2%. We also confirm that SLNP/ICG@M combined with NIR irradiation could activate mitochondria mediated apoptosis pathway, and the increased proliferation of CD4+ T cells and CD8+ T cells accompanied by immune activation further enhances PDT effect of SLNP/ICG@M. Herein, SLNP/ICG@M is a promising biomimetic nano drug delivery system for glioma targeted PDT therapy.

McAPRR2: The Key Regulator of Domesticated Pericarp Color in Bitter Gourd

Pericarp color is a crucial commercial trait influencing consumer preferences for bitter gourds. However, until now, the gene responsible for this trait has remained unidentified. In this study, we identified a gene (McAPRR2) controlling pericarp color via a genome-wide association study (GWAS) utilizing the resequencing data of 106 bitter gourd accessions. McAPRR2 exhibits three primary haplotypes: Hap1 is a wild type with a green pericarp, Hap2 is a SA (South Asian) and SEA (Southeast Asia) type with a green pericarp, and Hap3 is primarily a SEA type with a light green pericarp. The McAPRR2 haplotype is significantly correlated with both pericarp color and ecological type. Importantly, McAPRR2 with the light green pericarp demonstrated premature termination due to a 15 bp sequence insertion. The phylogenetic tree clustered according to pericarp color and ecological type, using SNPs located in the McAPRR2 gene and its promoter. High πwild/SEA and πSA/SEA values indicate high nucleotide diversity between wild and SEA types and between SA and SEA types in the McAPRR2 gene. The haplotypes, phylogenetic tree, and nucleotide diversity of McAPRR2 suggest that McAPRR2 has undergone domestication selection. This study identifies McAPRR2 as the key gene determining pericarp color in bitter gourds and introduces a novel insight that McAPRR2 is subject to domestication selection.

First Report of Broad Bean Wilt Virus 2 Infecting Balsam (Impatiens balsamina) in China

Balsam (Impatiens balsamina L.) is an ornamental plant cultivated extensively in China and elsewhere, but it has also been used as a medicinal plant for thousands of years (Qian et al., 2023). In 2022, an examination of 10 garden-grown I. balsamina plants in Chaoyang, Beijing, China revealed eight plants with blotches, mosaic symptoms, and deformed leaves (Fig. S1A). Total RNA was extracted from the symptomatic leaf tissue of these eight plants using the TRIzol reagent (Invitrogen, USA). Four RNA preparations (high quality and quantity) were combined for the small RNA sequencing analysis (TIANGEN Biotech Co., China). A total of 16,509,586 clean reads (18-30 nt) were obtained and assembled into larger contigs using Velvet 1.0.5. A search of the National Center for Biotechnology Information non-redundant database using BLASTX indicated 72, 24, and 19 contigs were homologous to broad bean wilt virus 2 (BBWV2), cucumber mosaic virus (CMV), and impatiens cryptic virus 1 (ICV1) sequences (Zheng et al., 2022), respectively. To verify the next-generation sequencing data, the following three sets of primer pairs were designed according to the contig sequences of these three viruses: CMV-F:5'-ATGGACAAATCTGAATCAACCAGTGC-3'/CMV-R: 5'-CCGTAAGCTGGATGGACAACC-3'; BBWV2-F:5'-CAATTTGGACAACTACAATTTGCC-3'/ BBWV2-R: 5'-GCTGAGTCTAAATCCCATCTATC-3'; and ICV1-F: 5'-CGCACAACT CTACAAT GACATGGTC-3'/ICV1-R: 5'-AGTTCCATCGTCCAGTAGGCG-3'. The primers were used to amplify CMV, BBWV2, and ICV1 sequences by reverse transcription-polymerase chain reaction (RT-PCR), with individual RNA preparations serving as the template. The CMV, BBWV2, and ICV1 target sequences were amplified from eight, four, and four samples, respectively (Fig. S1B). To evaluate virus infectivity, Nicotiana benthamiana seedlings were inoculated using a leaf tissue extract prepared from an infected I. balsamina plant. At 7 days post-inoculation, disease symptoms were detected on N. benthamiana systemic leaves (e.g., deformation and apical necrosis) (Fig. S1C). Confirmation tests involving RT-PCR indicated the N. benthamiana plants were infected with BBWV2 and CMV, but not with ICV1 (Fig. S1D). To obtain the complete BBWV2 genome sequence (RNA1 and RNA2), virus-specific PCR primers (Table S1) were designed to produce the terminal sequences via 5' and 3' rapid amplification of cDNA ends (RACE), which was completed using the SMARTer RACE 5'/3' Kit (Clontech, China). The RNA1 and RNA2 sequences comprised 5,957 nt (GenBank: OQ857921) and 3,614 nt (GenBank: OQ857922), respectively. The BLAST analyses revealed RNA1 and RNA2 were similar to sequences in other BBWV2 isolates (sequence identities of 78.88% to 95.15% and 80.83% to 91.51%, respectively). Using the neighbor-joining method and MEGA 7.0, the phylogenetic relationships between the BBWV2 isolated in this study and other isolates were determined on the basis of the full-length RNA1 and RNA2 sequences (Kumar et al., 2016). According to the RNA1 and RNA2 sequences, the BBWV2 isolated in this study was most closely related to the BBWV2 isolate from Gynura procumbens (GenBank: KX686589) and the BBWV2 isolate from Nicotiana tabacum (GenBank: KX650868), respectively (Fig. S1E). To the best of our knowledge, this is the first report of I. balsamina naturally infected with BBWV2 in China. The study findings may be useful for detecting BBWV2 in I. balsamina and for diagnosing and managing the associated disease. The authors declare no conflict of interest. Yanhong Qiu and Haijun Zhang contributed equally to this paper. Funding: This research was supported by the Beijing Academy of Agriculture and Forestry Foundation, China (KYCX202305, QNJJ202131, and KJCX20230214). References: Qian H.Q., et al. 2023. J Ethnopharmacol. 303. Zheng Y., et al. 2022. Arch Virol. 167: 2099-2102. Kumar et al. 2016. Mol Biol Evol. 33: 1870-1874.

Whole-Genome Sequencing of Fusarium oxysporum f. sp. cucumerinum Strain Race-4 Infecting Cucumber in China

Fusarium oxysporum f. sp. cucumerinum, which causes root and vascular wilting, is one of the most devastating diseases infecting cucumber. Here, we report the first genome resource with high-quality assembly for F. oxysporum f. sp. cucumerinum strain Race-4, which is primarily endemic to China. The genome was 59.11 Mb in size and consisted of 48 scaffolds with an N50 of 3.87 Mb using PacBio long reads (301.77×) sequencing, and encodes 14,898 proteins from analyzing RNA-seq data. Gene annotations identified pathogen-host interaction genes, fungal virulence factors, secreted proteins, transcription factors, and secondary metabolite biosynthesis gene. Moreover, functional genes reported in previous studies were also identified in the genome of Race-4. These genes and genome resource may play important roles in understanding F. oxysporum f. sp. cucumerinum-cucumber interactions and will be useful for further research.

[Analysis of genotypes on 850 newborns with SLC26A4 single-allele mutation and the phenotypes of those with second variant]

Objective: To clarify the phenotypes of the newborns with SLC26A4 single-allele mutation in deafness genetic screening and second variant; to analyze the SLC26A4 genotype and hearing phenotype. Methods: 850 newborns born in Beijing from April 2015 to December 2019 were included and there were 468 males and 382 females. They received genetic deafness screening for 9 or 15 variants, with the result of SLC26A4 single-allele mutation. Firstly, three step deafness gene sequencing was adopted in this work, i.e., the first step was "SLC26A4 gene whole exons and splice sites" sequencing; the second step was "SLC26A4 gene promoter, FOXI1 gene and KCNJ10 gene whole exons" sequencing; and the third step was detection for "SLC26A4 gene copy number variation". Secondly, we collected the results of newborn hearing screening for all patients with the second mutation found in the three step test, and conducted audiological examinations, such as acoustic immittance, auditory brainstem response and auditory steady state response. Thirdly, for novel/VUS mutations, we searched the international deafness gene database or software, such as DVD, ClinVar and Mutation Taster, to predict the pathogenicity of mutations according to the ACMG guideline. Lastly, we analyzed the relationship between genotype and phenotype of newborns with SLC26A4 single allele mutation. Results: Among 850 cases, the median age of diagnosis was 4 months. In the first step, 850 cases were sequenced. A total of 32 cases (3.76%, 32/850) of a second variants were detected, including 18 cases (2.12%, 18/850) with identified pathogenic variants; 832 cases were sequenced and 8 cases of KCNJ10 gene missense variants were detected among the second step. No missense mutations in the FOXI1 gene and abnormal SLC26A4 gene promoter were detected; the third step sequencing results were all negative. Genotypes and hearing phenotypes included 18 cases combined with the second clear pathogenic variant, 16 cases (16/18) referred newborn hearing screening and 2 cases (2/18) passed in both ears; degree of hearing loss consisted of 18 profound ears (18/36), 13 severe ears (13/36) and 5 moderate ears (5/36); audiogram patterns comprised 17 high frequency drop ears (17/36), 14 flat ears (14/36), 3 undistinguished ears (3/36), and 2 U shaped ears (2/36); 11 cases underwent imaging examination, all of which were bilateral enlarged vestibular aqueduct. As for 22 cases of other genotypes, all passed neonatal hearing screening and the hearing diagnosis was normal, including 9 cases with VUS or possibly novel benign variants, 8 cases with KCNJ10 double gene heterozygous variants, and 5 cases with double heterozygous variants. Conclusions: The probability of individuals with SLC26A4 single-allele variant who merge with a second pathogenic variant is 2.12%, all of which are SNV, which can provide scientific basis for the genetic diagnosis and genetic counseling of SLC26A4 variants. Those who have merged with second pathogenic variant are all diagnosed with sensorineural hearing loss. Patients with KCNJ10 gene mutations do not manifest hearing loss during the infancy, suggesting the need for further follow-up.

Genetic diversity analysis and variety identification using SSR and SNP markers in melon

Melon is an important horticultural crop with a pleasant aromatic flavor and abundance of health-promoting substances. Numerous melon varieties have been cultivated worldwide in recent years, but the high number of varieties and the high similarity between them poses a major challenge for variety evaluation, discrimination, as well as innovation in breeding. Recently, simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs), two robust molecular markers, have been utilized as a rapid and reliable method for variety identification. To elucidate the genetic structure and diversity of melon varieties, we screened out 136 perfect SSRs and 164 perfect SNPs from the resequencing data of 149 accessions, including the most representative lines worldwide. This study established the DNA fingerprint of 259 widely-cultivated melon varieties in China using Target-seq technology. All melon varieties were classified into five subgruops, including ssp. agrestis, ssp. melo, muskmelon and two subgroups of foreign individuals. Compared with ssp. melo, the ssp. agrestis varieties might be exposed to a high risk of genetic erosion due to their extremely narrow genetic background. Increasing the gene exchange between ssp. melo and ssp. agrestis is therefore necessary in the breeding procedure. In addition, analysis of the DNA fingerprints of the 259 melon varieties showed a good linear correlation (R² = 0.9722) between the SSR genotyping and SNP genotyping methods in variety identification. The pedigree analysis based on the DNA fingerprint of 'Jingyu' and 'Jingmi' series melon varieties was consistent with their breeding history. Based on the SNP index analysis, ssp. agrestis had low gene exchange with ssp. melo in chromosome 4, 7, 10, 11and 12, two specific SNP loci were verified to distinguish ssp. agrestis and ssp. melon varieties. Finally, 23 SSRs and 40 SNPs were selected as the core sets of markers for application in variety identification, which could be efficiently applied to variety authentication, variety monitoring, as well as the protection of intellectual property rights in melon.

Telomere-to-telomere genome assembly of bitter melon (Momordica charantia L. var. abbreviata Ser.) reveals fruit development, composition and ripening genetic characteristics

Momordica charantia L. var. abbreviata Ser. (Mca), known as bitter gourd or bitter melon, is a Momordica variety with medicinal value and belongs to the Cucurbitaceae family. In view of the lack of genomic information on bitter gourd and other Momordica species and to promote Mca genomic research, we assembled a 295.6-Mb telomere-to-telomere (T2T) high-quality Mca genome with six gap-free chromosomes after Hi-C correction. This genome is anchored to 11 chromosomes, which is consistent with the karyotype information, and comprises 98 contigs (N50 of 25.4 Mb) and 95 scaffolds (N50 of 25.4 Mb). The Mca genome harbors 19 895 protein-coding genes, of which 45.59% constitute predicted repeat sequences. Synteny analysis revealed variations involved in fruit quality during the divergence of bitter gourd. In addition, assay for transposase-accessible chromatin by high-throughput sequencing and metabolic analysis showed that momordicosides and other substances are characteristic of Mca fruit pulp. A combined transcriptomic and metabolomic analysis revealed the mechanisms of pigment accumulation and cucurbitacin biosynthesis in Mca fruit peels, providing fundamental molecular information for further research on Mca fruit ripening. This report provides a new genetic resource for Momordica genomic studies and contributes additional insights into Cucurbitaceae phylogeny.

Genome-Wide Identification of Powdery Mildew Responsive Long Non-Coding RNAs in Cucurbita pepo

Cucurbita pepo L. is an essential economic vegetable crop worldwide, and its production is severely affected by powdery mildew (PM). However, our understanding of the molecular mechanism of PM resistance in C. pepo is very limited. Long non-coding RNAs (lncRNAs) play an important role in regulating plant responses to biotic stress. Here, we systematically identified 2,363 reliably expressed lncRNAs from the leaves of PM-susceptible (PS) and PM-resistant (PR) C. pepo. The C. pepo lncRNAs are shorter in length and expressed at a lower level than the protein-coding transcripts. Among the 2,363 lncRNAs, a total of 113 and 146 PM-responsive lncRNAs were identified in PS and PR, respectively. Six PM-responsive lncRNAs were predicted as potential precursors of microRNAs (miRNAs). In addition, 58 PM-responsive lncRNAs were predicted as targets of miRNAs and one PM-responsive lncRNA was predicted as an endogenous target mimic (eTM). Furthermore, a total of 5,200 potential cis target genes and 5,625 potential trans target genes were predicted for PM-responsive lncRNAs. Functional enrichment analysis showed that these potential target genes are involved in different biological processes, such as the plant-pathogen interaction pathway, MAPK signaling pathway, and plant hormone signal transduction pathway. Taken together, this study provides a comprehensive view of C. pepo lncRNAs and explores the putative functions of PM-responsive lncRNAs, thus laying the foundation for further study of the regulatory mechanisms of lncRNAs responding to PM.

First report of tomato spotted wilt virus infecting water dropwort in China

Water dropwort (Oenanthe javanica) is an aquatic perennial plant that has been cultivated in many regions in Asia for thousands of years. In China, it is an economically important vegetable that has been consumed as food, while also being used as a folk remedy to alleviate diseases (Liu et al., 2021). In 2021, during a disease survey of a greenhouse in Beijing, China, chlorotic spots were detected on many water dropwort plants (Fig. S1A). Twenty-seven water dropwort samples were collected for the extraction of total RNA using the TRIzol reagent (Invitrogen, USA). High-quality RNA samples from three water dropwort plants were combined and used as the template for constructing a single small RNA library (BGI-Shenzhen Company, China). The Velvet 1.0.5 software was used to assemble the clean reads (18 to 28 nt) into larger contigs, which were then compared with the nucleotide sequences in the National Center database using the BLASTn algorithm. Thirty-eight contigs matched sequences in the tomato spotted wilt virus (TSWV) genome. No other viruses were detected. Twenty-seven leaf samples were analyzed in an enzyme-linked immunosorbent assay (ELISA) with anti-TSWV antibody (Agdia, USA), which revealed 17 positive reaction. Two sets of primer pairs targeting different parts of the S RNA (Table S1) was used to verify the TSWV infection on water dropwort by reverse transcription (RT)-PCR followed by Sanger sequencing (BGI-Shenzhen, China). The TSWV target sequences were amplified from 17 samples, which was consistent with the ELISA results. The sequenced 861-bp PCR product shared 99.8% nucleotide sequence identity with TSWV isolate MR-01 (MG593199), while the 441-bp amplicon shared a 99.2% nucleotide sequence identity with MR-01 (MG593199). To obtain the whole genome sequence of TSWV (S, M, and L RNA sequences), specific RT-PCR primers were designed (Table S1) and used to amplify their respective fragments from one representative sample (TSWV-water dropwort). The amplified products were inserted into PCE2TA/Blunt-Zero vector (Vazyme Biotech Co., Ltd, China) and then sequenced (BGI-Shenzhen, China). The S, M, and L RNA sequences were determined to be 2,952 nt (accession no. OM154969), 4,776 nt (accession no. OM154970), and 8,914 nt (accession no. OM154971), respectively. BLASTn analysis demonstrated that the whole genome sequence was highly conserved. The nucleotide identities between this isolate and other TSWV isolates ranged from 98.6% to 99.6% (S RNA), 98.9% to 99.2% (M RNA), and 97.3% to 98.7% (L RNA). Using MEGA 7.0, the phylogenetic relationships of TSWV were determined on the basis of the S, M, and L RNA full-length sequences (Kumar et al., 2016). In the S RNA derived phylogenetic tree, the water dropwort isolate was closely related to the MR-01 isolate from the USA (MG593199). In the M RNA and L RNA derived phylogenetic trees, the water dropwort isolate formed a branch with only a TSWV isolate from eggplant. Additionally, the M and L RNA sequences were most similar to sequences in TSWV isolates from China and Korea, respectively (Fig. S1B). To the best of our knowledge, this is the first report of water dropwort as a natural host for TSWV in China and the second report worldwide since the first finding in the Korea (Kil et al. 2020). TSWV has caused serious problems on many crops in the world, and the infection of TSWV on water dropwort in a greenhouse should not be looked lightly. Firstly, the virus can be passed on from generation to generation in infected water dropwort due to the vegetative propagation mode of the plant in production, thus threaten the production of this vegetable crop. In addition, infected water dropwort may serve as a reservoir for the virus, thus potentially posing a threat for causing TSWV spread in the affected greenhouses. The author(s) declare no conflict of interest. Funding: This research was supported by the Beijing Academy of Agriculture and Forestry Foundation, China (QNJJ202131, KJCX20200212, and KJCX20200113). References: Kil et al. 2020. Plant Pathol. J. 36: 67-75 Kumar et al. 2016. Mol Biol Evol, 33: 1870-1874 Liu et al. 2021. Horticulture Research. 8:1-17.

[Clinical features of children with Cunninghamella spp. infection: a case report and literature review]

We report a case of mucormycosis induced by Cunninghamella spp. infection in a ten-year-old girl with acute lymphoblastic leukemia, who developed fever and respiratory symptoms after chemotherapy and was diagnosed with invasive fungal disease. Peripheral blood DNA sequences were analyzed using metagenomic next-generation sequencing (mNGS), and by comparison with the Pathogens Metagenomics Database (PMDB), we identified Cunninghamella spp. with sequence number 514 as the pathogen. The patient was treated with amphotericin B combined with posaconazole and showed a favorable response. We searched Pubmed, Embase, CNKI, and Wanfang database for reports of cases of Cunninghamella spp. infection in children and retrieved 22 reported cases (including 12 males) with a median age of 13.5 (3-18) years. In these 22 cases, hematological malignancy was the most common underlying condition (19/22), and most of patients experienced an acute onset and rapid progression with respiratory symptoms (14/20) and fever (16/20) as the most common symptoms. CT imaging often showed unilateral lesions with varying imaging findings, including pulmonary nodules or masses, infiltrative changes, and pleural effusion. Definite diagnoses were established in 18 of the cases, and 4 had probable diagnoses; the lungs and skin were the most frequent organs compromised by the infection. A definite diagnosis of Cunninghamella spp. infection still relied on histopathological examination and fungal culture, but the molecular techniques including PCR and mNGS had shown potentials in the diagnosis. Almost all the cases received antifungal treatment after diagnosis (21/22), and 13 patients also underwent surgeries. Death occurred in 9 (42%) of the cases at a median of 19 (4-54) days after onset of the signs or symptoms. The patients receiving antifungal therapy combined with surgery had a high survival rate (9/13, 69%) than those with antifungal therapy alone (3/8, 37%). Invasive fungal disease is a common complication in immunoco-mpromised patients, but Cunninghamella spp. infection is rare and has a high mortality rate. In cases highly suspected of this disease, active diagnosis and early treatment are critical to improve the survival outcomes of the patients.

The vegetable SNP database: An integrated resource for plant breeders and scientists

Single nucleotide polymorphisms (SNPs) are widely used in genetic research and molecular breeding. To date, the genomes of many vegetable crops have been assembled, and hundreds of core germplasms for each vegetable have been sequenced. However, these data are not currently easily accessible because they are stored on different public databases. Therefore, a vegetable crop SNP database should be developed that hosts SNPs demonstrated to have a high success rate in genotyping for genetic research (herein, "alpha SNPs"). We constructed a database (VegSNPDB, http://www.vegsnpdb.cn/) containing the sequence data of 2032 germplasms from 16 vegetable crop species. VegSNPDB hosts 118,725,944 SNPs of which 4,877,305 were alpha SNPs. SNPs can be searched by chromosome number, position, SNP type, genetic population, or specific individuals, as well as the values of MAF, PIC, and heterozygosity. We hope that VegSNPDB will become an important SNP database for the vegetable research community.

Dual RNA-seq Reveals the Global Transcriptome Dynamics of Ralstonia solanacearum and Pepper (Capsicum annuum) Hypocotyls During Bacterial Wilt Pathogenesis

Bacterial wilt, caused by Ralstonia solanacearum, is a serious disease in pepper. However, the interaction between the pathogen and pepper remains largely unknown. This study aimed to gain insights into determinants of pepper susceptibility and R. solanacearum pathogenesis. We assembled the complete genome of R. solanacearum strain Rs-SY1 and identified 5,106 predicted genes, including 84 type III effectors (T3E). RNA-seq was used to identify differentially expressed genes (DEGs) in susceptible pepper CM334 at 1 and 5 days postinoculation (dpi) with R. solanacearum. Dual RNA-seq was used to simultaneously capture transcriptome changes in the host and pathogen at 3 and 7 dpi. A total of 1,400, 3,335, 2,878, and 4,484 DEGs of pepper (PDEGs) were identified in the CM334 hypocotyls at 1, 3, 5, and 7 dpi, respectively. Functional enrichment of the PDEGs suggests that inducing ethylene production, suppression of photosynthesis, downregulation of polysaccharide metabolism, and weakening of cell wall defenses may contribute to successful infection by R. solanacearum. When comparing in planta and nutrient agar growth of the R. solanacearum, 218 and 1,042 DEGs of R. solanacearum (RDEGs) were detected at 3 and 7 dpi, respectively. Additional analysis of the RDEGs suggested that enhanced starch and sucrose metabolism, and upregulation of virulence factors may promote R. solanacearum colonization. Strikingly, 26 R. solanacearum genes were found to have similar DEG patterns during a variety of host-R. solanacearum interactions. This study provides a foundation for a better understanding of the transcriptional changes during pepper-R. solanacearum interactions and will aid in the discovery of potential susceptibility and virulence factors.

First report of squash leaf curl China virus infecting tomato in China

Squash leaf curl China virus (SLCCNV) is a species in the genus Begomovirus that possess a bipartite genome. It is transmitted by the whitefly species Bemisia tabaci and infects cucurbit crops in various parts of the Old World (Wu et al., 2020). In 2020, tomato plants with curled, distorted and yellow leaves were found in a greenhouse in Shouguang, Shandong Province, China (Fig. S1). Leaves with these symptoms were collected from 11 plants and the total RNA was extracted with TRIzol reagent (Invitrogen, USA). Five RNA extracts of the highest quality were combined and a small RNA library was generated by the company (BGI-Shenzhen, China). About 22,338,920 clean reads (18-28nt) were acquired and assembled into larger contigs with the software Velvet 1.0.5. These were further compared against nucleotide sequences in the National Center for Biotechnology Information (NCBI) databases with BLASTn searches. Not unexpectedly, there were many assembled contigs that had high identities (90%-100% identities) with known tomato-infecting viruses, including 241 contigs matching tomato chlorosis virus, 26 contigs matching southern tomato virus, and 4 contigs matching tomato yellow leaf curl virus. However, 12 contigs had high identities (90%-100%) with the genomic DNA-A of SLCCNV, while 9 other contigs had high identities (90%-100%) with the genomic DNA-B of SLCCNV. To verify the presence of SLCCNV in tomato plants, two sets of primer pairs were designed according to the specific contigs assembled from derived small interfering RNAs (vsiRNAs). The primer pairs A742-F/A742-R (5'-GTAATACGAGCATCCGCACGGTAG-3'/5'-CGTGGAGGGCGAC AAACAGCTAACG-3') and B539-F/B539-R (5'-GCTACTTTCAAGGACGAAGAAGAGG-3'/5'-CG ACATAGATTTCTGGTCGGTGGGC-3') directed the amplification of 742 bp and 539 bp for DNA-A and DNA-B fragments, respectively, from the total genomic DNA of the 11 tomato samples. The DNA-A and DNA-B of SLCCNV were both detected from all of the tomato samples. After sequencing, the 742 bp PCR products shared 100% nucleotide sequence identity with the DNA-A of SLCCNV isolate GDXW (MW389919), whereas the PCR-amplified 539 bp fragments shared 100% nucleotide sequence identity with the DNA-B of SLCCNV isolate GDXW (MW389920). The full-length of DNA-A and DNA-B components were amplified with back-to-back primers A-F/A-R (Wu et al., 2020) and B-F/B-R (5'-GATAAACACGTCTCATTGCACCGC-3'/5'-GAGACGTGTTTATCAATATGGA CG-3'), respectively. The amplified fragments were further cloned into the PCE2TA/Blunt-Zero vector (Vazyme Biotech Co., China). After sequencing, the complete sequence of DNA-A was 2736 nt in length (MZ682117), while the DNA-B was 2718 nt in length (OK236348). The phylogenetic relationships of the DNA-A and DNA-B components were determined using MEGA 7 based on the full-length sequences of DNA-A and DNA-B, respectively (Kumar et al., 2016). Results showed that the DNA-A formed an independent cluster and was mostly related to the GDHY (MW389917) in the phylogenetic tree constructed using the neighbor-joining (NJ) method, while the DNA-B formed an independent cluster and was mostly related to the SLCCNV isolate BLDG (MW389928) and isolate GDBL (MW389922) (Fig. S2). The nt identities of DNA-A were also calculated with SDT v1.2 by comparison with other begomovirus sequences from the initial BLASTn analysis (Muhire et al., 2014), showing that the virus shared 99.4% sequence identity with SLCCNV isolate GDHY (MW389917). According to the current demarcation threshold for begomoviruses, recommended by the International Committee on Taxonomy of Viruses (ICTV) (91% nt identity) (Brown et al., 2015), this virus identified from tomato is a distinct strain of SLCCNV, designated SLCCNV-SDSG. To the best of our knowledge, this is the first report of a natural infection of SLCCNV on tomato in China. SLCCNV has caused serious problems in cucurbit production in some areas, so it will be important to investigate if tomato plays a role in the disease biology by serving as a reservoir host. The author(s) declare no conflict of interest. Funding: The funding for this research was supported by the Beijing Academy of Agriculture and Forestry Foundation, China (QNJJ202131, QNJJ201915, KJCX20200113). References: Brown et al. 2015. Arch Virol 160: 1593-1619 Kumar et al. 2016. Mol Biol Evol, 33: 1870-1874 Muhire et al. 2014. Plos One, 9 Wu et al. 2020. J Integr Agr, 19: 570-577.

Cucurbitaceae genome evolution, gene function and molecular breeding

The Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationship, and functional genes associated with important agronomic traits (e.g., fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

Dietary chitooligosaccharide supplementation alleviates intestinal barrier damage, and oxidative and immunological stress in lipopolysaccharide-challenged laying hens

This study aimed to investigate the effects of chitooligosaccharide (COS) on intestinal barrier, antioxidant capacity, and immunity of lipopolysaccharide (LPS)-challenged laying hens. A total of 360 Hy-line Brown laying hens (80-wk-old) were randomly divided into 5 groups with 6 replicates of 12 birds. Hens were fed a corn-soybean meal basal diet supplemented with different COS levels (0; 5; 10; 15; 20 mg/kg) for 8 wk. The results showed that 15 mg/kg COS administration elevated albumen height and Haugh unit (P < 0.05), and numerically optimized productive performance (P > 0.05), therefore, the dosage of 15 mg/kg was chosen for the subsequent experiment. Thereafter, 12 birds from non-supplemented group were randomly selected and assigned into 2 groups, and birds in each group were administered (1.5 mg/kg BW, i.p.) with saline (control group) or LPS (challenge group). Another 6 hens from 15 mg/kg COS-supplemented group were selected and injected with LPS in the same way. Compared with the control group, LPS-challenged birds exhibited elevated circulating diamine oxidase activity, and reduced jejunal villus height and ratio of villus height to crypt depth, and these indices were reversed to control levels by COS (P < 0.05). Also, LPS increased malondialdehyde accumulation and reduced several antioxidant enzyme activities in the intestinal mucosa (P < 0.05). Additionally, LPS increased jejunal secretory IgA and interferon-γ (IFN-γ), and ileal secretory IgA, IgM, and interleukin-1β (IL-1β) concentrations, whereas COS reduced jejunal IFN-γ and IL-1β, and ileal IgM levels (P < 0.05). Moreover, LPS down-regulated mRNA abundance of jejunal occludin and claudin 2, and upregulated expression of jejunal nuclear factor erythroid-2 related factor 2, superoxide dismutase 1, and IFN-γ as well as ileal IL-1β (P < 0.05). Besides, COS increased jejunal occludin and ileal claudin 2, nuclear factor erythroid-2 related factor 2, and heme oxygenase-1 expression, and decreased jejunal IFN-γ and IL-1β abundance (P < 0.05). These results suggested that COS could alleviate LPS-induced intestinal barrier impairment, and oxidative and immunological stress in laying hens.

[Analysis of epidemiological characteristics of human coronavirus in hospitalized children with respiratory tract infection in Hebei region]

Objective: To investigate the epidemiological characteristics of human coronavirus (HCoV) in hospitalized children with respiratory tract infection in Hebei region, providing evidence for the diagnosis and prevention of children with respiratory tract infection. Methods: A retrospective study was conducted on 1 062 HCoV positive children hospitalized for respiratory tract infection in Children's Hospital of Hebei Province from January 2015 to December 2020, aged from 33 days to 14 years, with a median age of 2 years. 27 932 (60.9%) were males and 17 944(39.1%) were females. And the gender, ages, seasonal distribution, HCoV-positive rates, co-detection distribution and clinical diagnosis of HCoV positive cases were analyzed by SPSS 25.0. Enumeration data were expressed by frequency and percentage; categorical variable were compared by the Pearson χ²test. Results: The overall HCoV-positive rate was 2.31% (1 062/45 876), which was 2.37% (662/27 932) in male children and 2.23% (400/17 944) in female children. There was no statistically significant difference between genders (χ²=0.916, P=0.339). Children at age groups<1 years (2.44%) and 1-<3 years (2.63%) had higher HCoV-positive rates than those at age groups 3-<5 years (1.97%) and ≥5 years (1.38%) (χ²=27.332,P<0.01). The HCoV-positive rates from 2015 to 2018 were 2.13%, 2.45%, 2.28% and 2.23%. The HCoV-positive rate of 2019 (1.71%) was significantly lower than in 2016 (χ²=12.05, P<0.01), 2017 (χ²=7.34, P=0.01) and 2018 (χ²=6.78, P=0.01), but there was no significant difference compared with 2015 (χ²=2.84, P=0.09). The HCoV-positive rate of 2020 (3.37%) was significantly higher than in 2015 (χ²=13.636, P<0.01), 2016 (χ²=11.099, P<0.01), 2017 (χ²=15.482, P<0.01), 2018(χ²=18.601, P<0.01) and 2019(χ²=45.580, P<0.01). The positive rate was highest in spring (March to May) in 2015 and 2017 to 2018. February to April and July to September of 2016 were the peak periods of positive detection. No obvious seasonal change was observed in 2019 and the HCoV-positive rate of 2020 was extremely low from January to July, following significantly increased from August to December. 26.37% (280/1 062) of HCoV were co-detected with other respiratory pathogens and the most frequently identified mixed detection was RSV. Three or more pathogens were detected in 7.34% (78/1 062) of the HCoV-positive samples. Bronchopneumonia and bronchiolitis were more frequently observed in the single HCoV positive (61.89% and 16.75%) children compared to co-detected children(34.29% and 9.64%)(χ²=63.394 and 8.228, P<0.01). However, compared to those with HCoV mono-detection, co-detected children were more likely to have severe pneumonia (4.6% and 47.14%) (χ²=280.171, P<0.01). Conclusions: HCoV is one of the respiratory pathogens in children in Hebei region and more prevalent in spring. The susceptible population of HCoV is mainly children under the age of 3 years old. HCoV often co-detects with other respiratory pathogens, and the co-infection is one of the risk factors of severe pneumonia in children with respiratory infection.