Overexpression of geranyl diphosphate synthase (PmGPPS1) boosts monoterpene and diterpene production involved in the response to pine wood nematode invasion

Outbreaks of pine wood nematode (PWN; Bursaphelenchus xylophilus) represent a severe biotic epidemic for the Pinus massoniana in China. When invaded by the PWN, the resistant P. massoniana might secret abundant oleoresin terpenoid to form certain defensive fronts for survival. However, the regulatory mechanisms of this process remain unclear. Here, the geranyl diphosphate synthase (PmGPPS1) gene was identified from resistant P. massoniana. Tissue-specific expression patterns of PmGPPS1 at transcript and protein level in resistant P. massoniana were determined by quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry. Functional characteristics analysis of PmGPPS1 was performed on transgenic Nicotiana benthamiana by overexpression, as genetic transformation of P. massoniana is, so far, not possible. In summary, we identified and functionally characterized PmGPPS1 from the resistant P. massoniana following PWN inoculation. Tissue-specific expression patterns and localization of PmGPPS1 indicated that it may play a positive role involved in the metabolic and defensive processes of oleoresin terpenes production in response to PWN attack. Furthermore, overexpression of PmGPPS1 may enhance the production of monoterpene, among which limonene reduced the survival of PWN in vitro. In addition, PmGPPS1 upregulated the expression level of key genes involved in mevalonic acid (MVA) pathway, the methylerythritol phosphate (MEP) pathway and gibberellins (GAs) biosynthesis to boost the growth and development of tobacco through a feedback regulation mechanism. Our results offered new insights into the pivotal role of the PmGPPS1 involved in terpene-based defense mechanisms responding to the PWN invasion in resistant P. massoniana and provided a new metabolic engineering scenario to improve monoterpene production in tobacco.

Comparative Transcriptome Analysis of Flower Senescence of Camellia lutchuensis

Background

Flower senescence is the last stage of flower development and affects the ornamental and economic value of flower plants. There is still less known on flower senescence of the ornamental plant Camellia lutchuensis, a precious species of Camellia with significant commercial application value.

Methods

Transcriptome sequencing was used to investigate the flower senescence in five developmental stages of C. lutchuensis.

Results

By Illumina HiSeq sequencing, we generated approximately 101.16 Gb clean data and 46649 differentially expressed unigenes. Based on the different expression pattern, differentially expressed unigenes were classified into 10 Sub Class. And Sub Class 9 including 8252 unigenes, was highly expressed in the flower senescent stage, suggesting it had a potential regulatory relationship of flower senescence. First, we found that ethylene biosynthesis genes ACSs, ACOs, receptor ETR genes and signaling genes EINs, ERFs all upregulated during flower senescence, suggesting ethylene might play a key role in the flower senescence of C. lutchuensis. Furthermore, reactive oxygen species (ROS) production related genes peroxidase (POD), lipase (LIP), polyphenoloxidase (PPO), and ROS scavenging related genes glutathione S-transferase (GST), glutathione reductase (GR) and superoxide dismutase (SOD) were induced in senescent stage, suggesting ROS might be involved in the flower senescence. Besides, the expression of monoterpenoid and isoflavonoid biosynthesis genes, transcription factors (WRKY, NAC, MYB and C₂H₂), senescence-associated gene SAG20 also were increased during flower senescence.

Conclusion

In C. lutchuensis, ethylene pathway might be the key to regulate flower senescence, and ROS signal might play a role in the flower senescence.

Analysis of Neutral B-Meson Decays into Two Muons

Branching fraction and effective lifetime measurements of the rare decay B_{s}^{0}→μ^{+}μ^{-} and searches for the decays B^{0}→μ^{+}μ^{-} and B_{s}^{0}→μ^{+}μ^{-}γ are reported using proton-proton collision data collected with the LHCb detector at center-of-mass energies of 7, 8, and 13 TeV, corresponding to a luminosity of 9  fb^{-1}. The branching fraction B(B_{s}^{0}→μ^{+}μ^{-})=(3.09_{-0.43-0.11}^{+0.46+0.15})×10^{-9} and the effective lifetime τ(B_{s}^{0}→μ^{+}μ^{-})=2.07±0.29±0.03  ps are measured, where the first uncertainty is statistical and the second systematic. No significant signal for B^{0}→μ^{+}μ^{-} and B_{s}^{0}→μ^{+}μ^{-}γ decays is found and upper limits B(B^{0}→μ^{+}μ^{-})<2.6×10^{-10} and B(B_{s}^{0}→μ^{+}μ^{-}γ)<2.0×10^{-9} at the 95% C.L. are determined, where the latter is limited to the range m_{μμ}>4.9  GeV/c^{2}. The results are in agreement with the standard model expectations.

HCV Infection Alters Salivary Gland Histology and Saliva Composition

Hepatitis C virus (HCV) infection is the most common blood-borne chronic infection in the United States. Chronic lymphocytic sialadenitis and sicca syndrome have been reported in chronic HCV infection. Up to 55% of these patients may have xerostomia; the mechanisms of the xerostomia and salivary gland (SG) hypofunction remain controversial. The objectives of this project are to establish if xerostomia associates with SG and HCV infection and to characterize the structural changes in SG and saliva composition. Eighteen HCV-infected patients with xerostomia were evaluated for SG dysfunction; 6 of these patients (patients 1–6) were further evaluated for SG histopathological changes and changes in saliva composition. The techniques used include clinical and laboratory assessment, SG ultrasonography, histological evaluation, sialochemical and proteomics analysis, and RNA in situ hybridization. All the HCV patients had low saliva flow, chronic sialadenitis, and SG fibrosis and lacked Sjögren syndrome (SS) characteristic autoantibodies. Further evaluation of a subgroup of 6 HCV patients (patients 1–6) demonstrated diffuse lymphocytic infiltrates that are predominantly CD8+ T cells with a significant increase in the number of inflammatory cells. Alcian Blue/periodic acid–Schiff staining showed significant changes in the ratio and intensity of the acinar secretory units of the HCV patients’ minor SG. The submandibular glands showed significant ultrasonographic abnormalities in the parenchyma relative to the parotid glands. Significant changes were also observed in the concentration of sodium and mucin 5b. Although no significant correlation was observed between the lymphocytic infiltrates and the years of HCV chronic infection, a positive correlation was observed between HCV RNA–positive epithelial cells and the years of HCV infection. Consistent with the low saliva flow and xerostomia, patients showed changes in several markers of SG acinar and ductal function. Changes in the composition of the saliva suggest that HCV infection can cause xerostomia by mechanisms distinct from SS.

The genome of oil-Camellia and population genomics analysis provide insights into seed oil domestication

BACKGROUND

As a perennial crop, oil-Camellia possesses a long domestication history and produces high-quality seed oil that is beneficial to human health. Camellia oleifera Abel. is a sister species to the tea plant, which is extensively cultivated for edible oil production. However, the molecular mechanism of the domestication of oil-Camellia is still limited due to the lack of sufficient genomic information.

RESULTS

To elucidate the genetic and genomic basis of evolution and domestication, here we report a chromosome-scale reference genome of wild oil-Camellia (2.95 Gb), together with transcriptome sequencing data of 221 cultivars. The oil-Camellia genome, assembled by an integrative approach of multiple sequencing technologies, consists of a large proportion of repetitive elements (76.1%) and high heterozygosity (2.52%). We construct a genetic map of high-density corrected markers by sequencing the controlled-pollination hybrids. Genome-wide association studies reveal a subset of artificially selected genes that are involved in the oil biosynthesis and phytohormone pathways. Particularly, we identify the elite alleles of genes encoding sugar-dependent triacylglycerol lipase 1, β-ketoacyl-acyl carrier protein synthase III, and stearoyl-acyl carrier protein desaturases; these alleles play important roles in enhancing the yield and quality of seed oil during oil-Camellia domestication.

CONCLUSIONS

We generate a chromosome-scale reference genome for oil-Camellia plants and demonstrate that the artificial selection of elite alleles of genes involved in oil biosynthesis contributes to oil-Camellia domestication.

The Inverted U-Shaped Association of Caffeine Intake with Serum Uric Acid in U.S. Adults

BACKGROUND AND AIMS

Caffeine is a worldwide popularly consumed constituent in foods that can exert physiological effects. However, previous researches about the relationship between caffeine intake and serum uric acid (SUA) were limited and controversial. Therefore, we sought to investigate that relationship in U.S. adults.

METHODS

In this cross-sectional study, the total sample of 7888 selected participants (3838 males and 4050 females) were identified from the National Health and Nutritional Examination Surveys (NHANES) 2015-2018. All subjects were tested for serum uric acid levels (μmmol/L), and their daily caffeine intakes (mg/d) were obtained by an average of two 24-hour dietary recalls. Multivariate linear regression models were used to evaluate the association between two variables in total subjects and subgroup analyses. Generalized additive models with smooth curve fittings were also performed.

RESULTS

Multivariate regression analyses showed caffeine intake was negatively correlated with SUA after adjustment of other confounders. The subgroup analyses stratified by gender showed the negative correlation of caffeine intake with SUA was statistically significant in males but not in females. Furthermore, we observed a nonlinear inverse association of caffeine intake with SUA (P nonlinear <0.001) in the generalized additive model, followed by an inverted U-shaped curve (inflection point: 60.5mg/d) for all participants. This inverted U-shaped relationship between them could also be found in both genders, individuals aged below 60 years old, overweight (BMI of 25 to 30), and Non-Hispanic White individuals.

CONCLUSIONS

This study indicated that caffeine intake exhibited an inverse correlation with SUA, especially in males. In addition, this inverse relationship was nonlinear, which followed an inverted U-shaped curve.

The complete chloroplast genome of Camellia fluviatilis (Theaceae), a wild oil-Camellia species

Camellia fluviatilis is an important shrub producing edible seed oil, which is widely cultivated in South China. In this study, the complete chloroplast genome was sequenced and analyzed based on the Illumina HiSeq platform. The results showed that the complete chloroplast genome is 157,041 bp with 37.29% GC content, including a large single copy (LSC) region of 86,718 bp, a small single copy (SSC) region of 18,293 bp, and a pair of inverted repeats (IRs) regions of 26,015 bp. There are 128 genes in the chloroplast genome of C. fluviatilis, including 83 protein-coding genes, 8 ribosomal RNAs, and 37 transfer RNAs. Phylogenetic analysis revealed that C. fluviatilis is closely related to C. lanceoleosa, indicating that both belong to the Sect. Paracamellia Sealy.

Global Burden of Disease from Major Air Pollution Sources (GBD MAPS): A Global Approach

Ambient fine particulate matter (particles <2.5 μm in aerodynamic diameter [PM2.5]) is the world's leading environmental health risk factor. Reducing the PM2.5 disease burden requires specific strategies that target dominant sources across multiple spatial scales. The Global Burden of Disease from Major Air Pollution Sources (GBD MAPS) project provides a contemporary and comprehensive evaluation of contributions to the ambient PM2.5 disease burden from source sectors and fuels across 21 regions, 204 countries, and 200 subnational areas. We first derived quantitative contributions from 24 emission sensitivity simulations using an updated global atmospheric chemistry-transport model, input with a newly developed detailed anthropogenic emissions dataset that includes emissions specific to source sector and fuels. These simulation results were integrated with newly available high-resolution satellite-derived PM2.5 exposure estimates and disease-specific concentration-response relationships consistent with the GBD project to quantify contributions of specific source sector and fuel to the ambient PM2.5 disease burden across all regions, countries, and subnational areas. To improve the transparency and reproducibility of this and future work, we publicly provided the global atmospheric chemistry-transport model source code, emissions dataset and emissions model source code, analysis scripts, and source sensitivity results, and further described the emissions dataset and source contribution results in two publications.

We found that nearly 1.05 million (95% uncertainty interval [UI]: 0.74-1.36 million) deaths worldwide (27.3% of the total mortality attributable to PM2.5) would be avoidable by eliminating fossil fuel combustion, with coal contributing over half of that burden. Residential (19.2%; 736,000 deaths [95% UI: 521,000-955,000]), industrial (11.7%; 448,000 deaths [95% UI: 318,000-582,000]), and energy (10.2%; 391,000 deaths [95% UI: 277,000-507,000]) sector emissions are among the dominant global sources Uncertainty in these estimates reflects those of the input datasets. Regions with the largest anthropogenic contributions generally have the highest numbers of attributable deaths, which clearly demonstrates the importance of reducing these emissions to realize reductions in global air pollution and its disease burden.

The Flavonoid Biosynthesis Network in Plants

Flavonoids are an important class of secondary metabolites widely found in plants, contributing to plant growth and development and having prominent applications in food and medicine. The biosynthesis of flavonoids has long been the focus of intense research in plant biology. Flavonoids are derived from the phenylpropanoid metabolic pathway, and have a basic structure that comprises a C15 benzene ring structure of C6-C3-C6. Over recent decades, a considerable number of studies have been directed at elucidating the mechanisms involved in flavonoid biosynthesis in plants. In this review, we systematically summarize the flavonoid biosynthetic pathway. We further assemble an exhaustive map of flavonoid biosynthesis in plants comprising eight branches (stilbene, aurone, flavone, isoflavone, flavonol, phlobaphene, proanthocyanidin, and anthocyanin biosynthesis) and four important intermediate metabolites (chalcone, flavanone, dihydroflavonol, and leucoanthocyanidin). This review affords a comprehensive overview of the current knowledge regarding flavonoid biosynthesis, and provides the theoretical basis for further elucidating the pathways involved in the biosynthesis of flavonoids, which will aid in better understanding their functions and potential uses.

Transcriptomic Analysis Reveals Key Candidate Genes Related to Seed Abortion in Chinese Jujube (Ziziphus jujuba Mill.)

Background

Seed abortion is a common phenomenon in Chinese jujube that seriously hinders the process of cross-breeding. However, the molecular mechanisms of seed abortion remain unclear in jujube.

Methods

Here, we performed transcriptome sequencing using eight flower and fruit tissues at different developmental stages in Ziziphus jujuba Mill. ‘Zhongqiusucui’ to identify key genes related to seed abortion. Histological analysis revealed a critical developmental process of embryo abortion after fertilization.

Results

Comparisons of gene expression revealed a total of 14,012 differentially expressed genes. Functional enrichment analyses of differentially expressed genes between various sample types uncovered several important biological processes, such as embryo development, cellular metabolism, and stress response, that were potentially involved in the regulation of seed abortion. Furthermore, gene co-expression network analysis revealed a suite of potential key genes related to ovule and seed development. We focused on three types of candidate genes, agamous subfamily genes, plant ATP-binding cassette subfamily G transporters, and metacaspase enzymes, and showed that the expression profiles of some members were associated with embryo abortion.

Conclusion

This work generates a comprehensive gene expression data source for unraveling the molecular mechanisms of seed abortion and aids future cross-breeding efforts in jujube.

D6 protein kinase in root xylem benefiting resistance to Fusarium reveals infection and defense mechanisms in tung trees

Fusarium oxysporum, a global soil-borne pathogen, causes severe disease in various cultivated plants. The mechanism underlying infection and resistance remains largely elusive. Vernicia fordii, known as the tung tree, suffers from disease caused by F. oxysporum f. sp. fordiis (Fof-1), while its sister species V. montana displays high resistance to Fof-1. To investigate the process of infection and resistance ability, we demonstrated that Fof-1 can penetrate the epidermis of root hairs and then centripetally invade the cortex and phloem in both species. Furthermore, Fof-1 spread upwards through the root xylem in susceptible V. fordii trees, whereas it failed to infect the root xylem in resistant V. montana trees. We found that D6 PROTEIN KINASE LIKE 2 (VmD6PKL2) was specifically expressed in the lateral root xylem and was induced after Fof-1 infection in resistant trees. Transgenic analysis in Arabidopsis and tomato revealed that VmD6PKL2 significantly enhanced resistance in both species, whereas the d6pkl2 mutant displayed reduced resistance against Fof-1. Additionally, VmD6PKL2 was identified to interact directly with synaptotagmin (VmSYT3), which is specifically expressed in the root xylem and mediates the negative regulation responding to Fof-1. Our data suggested that VmD6PKL2 could act as a resistance gene against Fof-1 through suppression of VmSYT3-mediated negative regulation in the lateral root xylem of the resistant species. These findings provide novel insight into Fusarium wilt resistance in plants.

Characterization and phylogenetic analysis of the complete chloroplast genome of Camellia chrysanthoides (Theaceae)

Camellia chrysanthoides H.T.Chang 1979 is an economic species for its high ornamental and medicinal values. In the present work, the complete chloroplast (cp) genome sequence of C. chrysanthoides was assembled and characterized using Illumina paired-end sequencing data. The total cp genome of C. chrysanthoides is 157,439 bp in size, consisting of a small single copy (SSC) and a large single copy (LSC) separated by a pair of inverted repeats (IRs) with 18,265 bp, 88,162 bp, and 25,506 bp, respectively. Further annotation revealed the cp genome encoded 124 genes, including 82 protein-coding genes, 34 tRNA genes, and eight rRNA genes, and the overall GC content of the cp genome is 37.31%. Phylogenetic analysis based on 80 protein-coding genes shows that C. chrysanthoidesis closely related to C. azalea in the genus Camellia.

Branching Fraction Measurements of the Rare B_{s}^{0}→ϕμ^{+}μ^{-} and B_{s}^{0}→f_{2}^{'}(1525)μ^{+}μ^{-} Decays

The branching fraction of the rare B_{s}^{0}→ϕμ^{+}μ^{-} decay is measured using data collected by the LHCb experiment at center-of-mass energies of 7, 8, and 13 TeV, corresponding to integrated luminosities of 1, 2, and 6  fb^{-1}, respectively. The branching fraction is reported in intervals of q^{2}, the square of the dimuon invariant mass. In the q^{2} region between 1.1 and 6.0  GeV^{2}/c^{4}, the measurement is found to lie 3.6 standard deviations below a standard model prediction based on a combination of light cone sum rule and lattice QCD calculations. In addition, the first observation of the rare B_{s}^{0}→f_{2}^{'}(1525)μ^{+}μ^{-} decay is reported with a statistical significance of 9 standard deviations and its branching fraction is determined.

Interactions between poloxamer, PEOx-PPOy-PEOx, and non-ionic surfactant, sucrose monolaurate: A study on potential allergenic effect using model phospholipid membrane

Sugar-based surfactants are involved in skin related allergy cases in the past decade. Skin irritation starts with the interaction of the surfactant with the skin lipids leading to lipid emulsification and eventual barrier damage. Polymers or co-surfactants can be used to mitigate the allergenic effect but the mechanism of formulation mildness on skin remains unclear. We have used the quartz crystal microbalance (QCM) together with dissipative particle dynamics (DPD) simulation, small angle x-ray scattering (SAXS) as well as cell viability tests to decipher the interactions between poloxamers and sucrose monolaurate (SML), and how these interactions could prevent the disruption of a model supported phospholipid bilayer (SLB). Poloxamer addition to the SML solution can delay or totally prevent the disruption of the SLB depending on poloxamer type and concentration. Poloxamer P407 (Pluronic® F127) delays the onset of disruption while poloxamer P188 (Pluronic® F68) does not preserve the bilayer integrity even at high concentration of up to 15% w/w. Preservation of the SLB is likely due to the differences in the aggregates formation between SML-F127 and SML-F68 mixtures with corresponding retarded motion of SML micelles through the SML-F127 polymer matrix that improved cell viability.

Upregulation of FoxO6 in nucleus pulposus cells promotes DNA damage repair via activation of RAD51

OBJECTIVE

DNA damage is an essential risk for intervertebral disc degeneration (IDD). Here, we attempted to uncover the effect of FoxO6 and RAD51 on the DNA damage repair of nucleus pulposus (NP) cells in IDD.

PATIENTS AND METHODS

We collected the human NP tissues of different degeneration degrees and tested the collagen II, FoxO6, and RAD51 expression. Besides, the IL-1β induced NP cell model was also used to elucidate the degenerative progress in vitro. We used Chromatin immunoprecipitation (ChIP) and luciferase reporter assay to confirm whether the FoxO6 protein could enhance the RAD51 expression by binding to its promoter. The FoxO6 gene was upregulated in NP cells by vectors transfection. Immunofluorescence staining was used to measure the RAD51 and γH2AX foci formation. Besides, the typical NP cell gene expression was analyzed by RT-PCR. Cell proliferation was determined by CCK-8, and the cell cycle distribution was determined by flow cytometry.

RESULTS

Like collagen II, FoxO6 and RAD51 expression were all decreased both in the severe degenerated NP tissue and in the IL-1β treated NP cells. Upregulation of FoxO6 gene in NP cells enhanced the RAD51 expression via activating the promoter region and inhibited the DNA damage marker γH2AX formation. FoxO6 upregulation alleviated the loss of collagen II, aggrecan, SOD1, and CAT, and suppressed the increase of collagen I/X, TNF-α, and IL-1β expression, which was affected by IL-1β. Besides, FoxO6 also helped the proliferation and cell cycle of NP cells with the activation of RAD51.

CONCLUSIONS

Upregulation of FoxO6 promotes the DNA repair and maintains the typical phenotype of NP cells, via somehow the mediation of RAD51.

Observation of the Mass Difference between Neutral Charm-Meson Eigenstates

A measurement of mixing and CP violation in neutral charm mesons is performed using data reconstructed in proton-proton collisions collected by the LHCb experiment from 2016 to 2018, corresponding to an integrated luminosity of 5.4  fb^{-1}. A total of 30.6 million D^{0}→K_{S}^{0}π^{+}π^{-} decays are analyzed using a method optimized for the measurement of the mass difference between neutral charm-meson eigenstates. Allowing for CP violation in mixing and in the interference between mixing and decay, the mass and decay-width differences are measured to be x_{CP}=[3.97±0.46(stat)±0.29(syst)]×10^{-3} and y_{CP}=[4.59±1.20(stat)±0.85(syst)]×10^{-3}, respectively. The CP-violating parameters are measured as Δx=[-0.27±0.18(stat)±0.01(syst)]×10^{-3} and Δy=[0.20±0.36(stat)±0.13(syst)]×10^{-3}. This is the first observation of a nonzero mass difference in the D^{0} meson system, with a significance exceeding seven standard deviations. The data are consistent with CP symmetry and improve existing constraints on the associated parameters.

Potential biomarkers and lncRNA-mRNA regulatory networks in invasive growth hormone-secreting pituitary adenomas

PURPOSE

Growth hormone-secreting pituitary adenomas (GH-PAs) are common subtypes of functional PAs. Invasive GH-PAs play a key role in restricting poor outcomes. The transcriptional changes in GH-PAs were evaluated.

METHODS

In this study, the transcriptome analysis of six different GH-PA samples was performed. The functional roles, co-regulatory network, and chromosome location of differentially expressed (DE) genes in invasive GH-PAs were explored.

RESULTS

Bioinformatic analysis revealed 101 DE mRNAs and 70 DE long non-coding RNAs (lncRNAs) between invasive and non-invasive GH-PAs. Functional enrichment analysis showed that epithelial cell differentiation and development pathways were suppressed in invasive GH-PAs, whereas the pathways of olfactory transduction, retinol metabolism, drug metabolism-cytochrome P450, and metabolism of xenobiotics by cytochrome P450 had an active trend. In the protein-protein interaction network, 11 main communities were characterized by cell- adhesion, -motility, and -cycle; transport process; phosphorus and hormone metabolic processes. The SGK1 gene was suggested to play a role in the invasiveness of GH-PAs. Furthermore, the up-regulated genes OR51B6, OR52E4, OR52E8, OR52E6, OR52N2, MAGEA6, MAGEC1, ST8SIA6-AS1, and the down-regulated genes GAD1-AS1 and SPINT1-AS1 were identified in the competing endogenous RNA network. The RT-qPCR results further supported the aberrant expression of those genes. Finally, the enrichment of DE genes in chromosome 11p15 and 12p13 regions were detected.

CONCLUSION

Our findings provide a new perspective for studies evaluating the underlying mechanism of invasive GH-PAs.

Observation of New Resonances Decaying to J/ψK^{+} and J/ψϕ

The first observation of exotic states with a new quark content cc[over ¯]us[over ¯] decaying to the J/ψK^{+} final state is reported with high significance from an amplitude analysis of the B^{+}→J/ψϕK^{+} decay. The analysis is carried out using proton-proton collision data corresponding to a total integrated luminosity of 9  fb^{-1} collected by the LHCb experiment at center-of-mass energies of 7, 8, and 13 TeV. The most significant state, Z_{cs}(4000)^{+}, has a mass of 4003±6_{-14}^{+4}  MeV, a width of 131±15±26  MeV, and spin parity J^{P}=1^{+}, where the quoted uncertainties are statistical and systematic, respectively. A new 1^{+}  X(4685) state decaying to the J/ψϕ final state is also observed with high significance. In addition, the four previously reported J/ψϕ states are confirmed and two more exotic states, Z_{cs}(4220)^{+} and X(4630), are observed with significance exceeding 5 standard deviations.