Structure and function of Mycobacterium tuberculosis EfpA as a lipid transporter and its inhibition by BRD-8000.3.

EfpA, the first major facilitator superfamily (MFS) protein identified in Mycobacterium tuberculosis (Mtb), is an essential efflux pump implicated in resistance to multiple drugs. EfpA-inhibitors have been developed to kill drug-tolerant Mtb. However, the biological function of EfpA has not yet been elucidated. Here, we present the cryo-EM structures of EfpA complexed with lipids or the inhibitor BRD-8000.3 at resolutions of 2.9 Å and 3.4 Å, respectively. Unexpectedly, EfpA forms an antiparallel dimer. Functional studies reveal that EfpA is a lipid transporter and BRD-8000.3 inhibits its lipid transport activity. Intriguingly, the mutation V319F, known to confer resistance to BRD-8000.3, alters the expression level and oligomeric state of EfpA. Based on our results and the observation of other antiparallel dimers in the MFS family, we propose an antiparallel-function model of EfpA. Collectively, our work provides structural and functional insights into EfpA's role in lipid transport and drug resistance, which would accelerate the development of antibiotics against this promising drug target.

ZmPTOX1, a plastid terminal oxidase, contributes to redox homeostasis during seed development and germination.

Maize plastid terminal oxidase1 (ZmPTOX1) plays a pivotal role in seed development by upholding redox balance within seed plastids. This study focuses on characterizing the white kernel mutant 3735 (wk3735) mutant, which yields pale-yellow seeds characterized by heightened protein but reduced carotenoid levels, along with delayed germination compared to wild-type (WT) seeds. We successfully cloned and identified the target gene ZmPTOX1, responsible for encoding maize PTOX-a versatile plastoquinol oxidase and redox sensor located in plastid membranes. While PTOX's established role involves regulating redox states and participating in carotenoid metabolism in Arabidopsis leaves and tomato fruits, our investigation marks the first exploration of its function in storage organs lacking a photosynthetic system. Through our research, we validated the existence of plastid-localized ZmPTOX1, existing as a homomultimer, and established its interaction with ferredoxin-NADP+ oxidoreductase 1 (ZmFNR1), a crucial component of the electron transport chain (ETC). This interaction contributes to the maintenance of redox equilibrium within plastids. Our findings indicate a propensity for excessive accumulation of reactive oxygen species (ROS) in wk3735 seeds. Beyond its known role in carotenoids' antioxidant properties, ZmPTOX1 also impacts ROS homeostasis owing to its oxidizing function. Altogether, our results underscore the critical involvement of ZmPTOX1 in governing seed development and germination by preserving redox balance within the seed plastids.

Population-level gene expression can repeatedly link genes to functions in maize.

Transcriptome-wide association studies (TWAS) can provide single gene resolution for candidate genes in plants, complementing genome-wide association studies (GWAS) but efforts in plants have been met with, at best, mixed success. We generated expression data from 693 maize genotypes, measured in a common field experiment, sampled over a 2-h period to minimize diurnal and environmental effects, using full-length RNA-seq to maximize the accurate estimation of transcript abundance. TWAS could identify roughly 10 times as many genes likely to play a role in flowering time regulation as GWAS conducted data from the same experiment. TWAS using mature leaf tissue identified known true-positive flowering time genes known to act in the shoot apical meristem, and trait data from a new environment enabled the identification of additional flowering time genes without the need for new expression data. eQTL analysis of TWAS-tagged genes identified at least one additional known maize flowering time gene through trans-eQTL interactions. Collectively these results suggest the gene expression resource described here can link genes to functions across different plant phenotypes expressed in a range of tissues and scored in different experiments.

Structural and inhibitor sensitivity analysis of influenza B-like viral neuraminidases derived from Asiatic toad and spiny eel.

Influenza virus neuraminidase (NA) is an important target for antiviral development because it plays a crucial role in releasing newly assembled viruses. Two unique influenza-like virus genomes were recently reported in the Wuhan Asiatic toad and Wuhan spiny eel. Their NA genes appear to be highly divergent from all known influenza NAs, raising key questions as to whether the Asiatic toad influenza-like virus NA (tNA) and spiny eel NA (eNA) have canonical NA activities and structures and whether they show sensitivity to NA inhibitors (NAIs). Here, we found that both tNA and eNA have neuraminidase activities. A detailed structural analysis revealed that tNA and eNA present similar overall structures to currently known NAs, with a conserved calcium binding site. Inhibition assays indicated that tNA is resistant to NAIs, while eNA is still sensitive to NAIs. E119 is conserved in canonical NAs. The P119E substitution in tNA can restore sensitivity to NAIs, and, in contrast, the E119P substitution in eNA decreased its sensitivity to NAIs. The structures of NA-inhibitor complexes further provide a detailed insight into NA-inhibitor interactions at the atomic level. Moreover, tNA and eNA have unique N-glycosylation sites compared with canonical NAs. Collectively, the structural features, NA activities, and sensitivities to NAIs suggest that fish- and amphibian-derived influenza-like viruses may circulate in these vertebrates. More attention should be paid to these influenza-like viruses because their NA molecules may play roles in the emergence of NAI resistance.

Flowering time regulator qFT13-3 involved in soybean adaptation to high latitudes.

Soybean is a short-day plant that typically flowers earlier when exposed to short-day conditions. However, the identification of genes associated with earlier flowering time but without a yield penalty is rare. In this study, we conducted genome-wide association studies (GWAS) using two re-sequencing datasets that included 113 wild soybeans (G. soja) and 1192 cultivated soybeans (G. max), respectively, and simultaneously identified a candidate flowering gene, qFT13-3, which encodes a protein homologous to the pseudo-response regulator (PRR) transcription factor. We identified four major haplotypes of qFT13-3 in the natural population, with haplotype H4 (qFT13-3H4) being lost during domestication, while qFT13-3H1 underwent natural and artificial selection, increasing in proportion from 4.5% in G. soja to 43.8% in landrace and to 81.9% in improve cultivars. Notably, most cultivars harbouring qFT13-3H1 were located in high-latitude regions. Knockout of qFT13-3 accelerated flowering and maturity time under long-day conditions, indicating that qFT13-3 functions as a flowering inhibitor. Our results also showed that qFT13-3 directly downregulates the expression of GmELF3b-2 which is a component of the circadian clock evening complex. Field trials revealed that the qft13-3 mutants shorten the maturity period by 11 days without a concomitant penalty on yield. Collectively, qFT13-3 can be utilized for the breeding of high-yield cultivars with a short maturity time suitable for high latitudes.

Polarity-Evolution Control and Luminescence Regulation in Multiple-Site Hydrogen-Bonded Organic Frameworks.

Hydrogen-bonded organic frameworks (HOFs) have shown great potential in separation, sensing and host-guest chemistry, however, the pre-design of HOFs remains challenging due to the uncertainty of solvents' participation in framework formation. Herein, the polarity-evolution-controlled framework/luminescence regulation is demonstrated based on multiple-site hydrogen-bonded organic frameworks. Several distinct HOFs were prepared by changing bonding modes of building units via the evolution of electrostatic forces induced by various solvent polarities. High-polar solvents with strong electrostatic attraction to surrounding units showed the tendency to form cage structures, while low-polar solvents with weak electrostatic attraction only occupy hydrogen-bond sites, conducive to the channel formation. Furthermore, the conformation of optical building unit can be adjusted by affecting the solvent polarity, generating different luminescence outputs. These results pave the way for the rational design of ideal HOFs with on-demand framework regulation and luminescence properties.

Deciphering the genetic basis of resistance to soybean cyst nematode combining IBD and association mapping.

KEY MESSAGE: IBD analysis clarified the dynamics of chromosomal recombination during the ZP pedigree breeding process and identified ten genomic regions resistant to SCN race3 combining association mapping. Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most devastating pathogens for soybean production worldwide. The cultivar Zhongpin03-5373 (ZP), derived from SCN-resistant progenitor parents, Peking, PI 437654 and Huipizhi Heidou, is an elite line with high resistance to SCN race3. In the current study, a pedigree variation map was generated for ZP and its ten progenitors using 3,025,264 high-quality SNPs identified from an average of 16.2 × re-sequencing for each genome. Through identity by decent (IBD) tracking, we showed the dynamic change of genome and detected important IBD fragments, which revealed the comprehensively artificial selection of important traits during ZP breeding process. A total of 2,353 IBD fragments related to SCN resistance including SCN-resistant genes rhg1, rhg4 and NSFRAN07 were identified based on the resistant-related genetic paths. Moreover, 23 genomic regions underlying resistance to SCN race3 were identified by genome-wide association study (GWAS) in 481 re-sequenced cultivated soybeans. Ten common loci were found by both IBD tracking and GWAS analysis. Haplotype analysis of 16 potential candidate genes suggested a causative SNP (C/T, - 1065) located in the promoter of Glyma.08G096500 and encoding a predicted TIFY5b-related protein on chr8 was highly correlated with SCN race3 resistance. Our results more thoroughly elucidated the dynamics of genomic fragments during ZP pedigree breeding and the genetic basis of SCN resistance, which will provide useful information for gene cloning and the development of resistant soybean cultivars using a marker-assisted selection approach.

The elite variations in germplasms for soybean breeding.

The genetic base of soybean cultivars (Glycine max (L.) Merr.) has been narrowed through selective domestication and specific breeding improvement, similar to other crops. This presents challenges in breeding new cultivars with improved yield and quality, reduced adaptability to climate change, and increased susceptibility to diseases. On the other hand, the vast collection of soybean germplasms offers a potential source of genetic variations to address those challenges, but it has yet to be fully leveraged. In recent decades, rapidly improved high-throughput genotyping technologies have accelerated the harness of elite variations in soybean germplasm and provided the important information for solving the problem of a narrowed genetic base in breeding. In this review, we will overview the situation of maintenance and utilization of soybean germplasms, various solutions provided for different needs in terms of the number of molecular markers, and the omics-based high-throughput strategies that have been used or can be used to identify elite alleles. We will also provide an overall genetic information generated from soybean germplasms in yield, quality traits, and pest resistance for molecular breeding.

Time series canopy phenotyping enables the identification of genetic variants controlling dynamic phenotypes in soybean.

Advances in plant phenotyping technologies are dramatically reducing the marginal costs of collecting multiple phenotypic measurements across several time points. Yet, most current approaches and best statistical practices implemented to link genetic and phenotypic variation in plants have been developed in an era of single-time-point data. Here, we used time-series phenotypic data collected with an unmanned aircraft system for a large panel of soybean (Glycine max (L.) Merr.) varieties to identify previously uncharacterized loci. Specifically, we focused on the dissection of canopy coverage (CC) variation from this rich data set. We also inferred the speed of canopy closure, an additional dimension of CC, from the time-series data, as it may represent an important trait for weed control. Genome-wide association studies (GWASs) identified 35 loci exhibiting dynamic associations with CC across developmental stages. The time-series data enabled the identification of 10 known flowering time and plant height quantitative trait loci (QTLs) detected in previous studies of adult plants and the identification of novel QTLs influencing CC. These novel QTLs were disproportionately likely to act earlier in development, which may explain why they were missed in previous single-time-point studies. Moreover, this time-series data set contributed to the high accuracy of the GWASs, which we evaluated by permutation tests, as evidenced by the repeated identification of loci across multiple time points. Two novel loci showed evidence of adaptive selection during domestication, with different genotypes/haplotypes favored in different geographic regions. In summary, the time-series data, with soybean CC as an example, improved the accuracy and statistical power to dissect the genetic basis of traits and offered a promising opportunity for crop breeding with quantitative growth curves.

TWAS results are complementary to and less affected by linkage disequilibrium than GWAS.

A genome-wide association study (GWAS) is used to identify genetic markers associated with phenotypic variation. In contrast, a transcriptome-wide association study (TWAS) detects associations between gene expression levels and phenotypic variation. It has previously been shown that in the cross-pollinated species, maize (Zea mays), GWAS, and TWAS identify complementary sets of trait-associated genes, many of which exhibit characteristics of true positives. Here, we extend this conclusion to the self-pollinated species, Arabidopsis thaliana and soybean (Glycine max). Linkage disequilibrium (LD) can result in the identification, via GWAS, of false-positive associations. In all three analyzed plant species, most trait-associated genes identified via TWAS are well separated physically from other candidate genes. Hence, TWAS is less affected by LD than is GWAS, demonstrating that TWAS is particularly well suited for association studies in genomes with slow rates of LD decay, such as soybean. TWAS is reasonably robust to the plant organs/tissues used to determine expression levels. In summary, this study confirms that TWAS is a promising approach for accurate gene-level association mapping in plants that is complementary to GWAS, and established that TWAS can exhibit substantial advantages relative to GWAS in species with slow rates of LD decay.

Hydration-Facilitated Coordination Tuning of Metal-Organic Frameworks toward Water-Responsive Fluorescence and Proton Conduction.

Developing smart stimuli-responsive metal-organic frameworks (MOFs) with diversified induced readable signals is highly desirable; however, reported multimode responsive MOFs are always achieved under strong environmental stimulations, making it difficult to keep MOF structures stable for practical applications. Herein, we reported a hydration-facilitated coordination tuning strategy to achieve the dual-mode water response in fluorescence and proton conduction from a single MOF. The designed MOF permitted reversible single-crystal transformation via the controllable hydration effect on metal nodes. The change in coordination modes leads to the regulation on conformations of optical ligands, contributing to the switch of fluorescence emissions. Moreover, the hydration effect adds additional hydrogen-bond sites in channels and optimizes hydrogen-bond networks, abruptly enhancing the proton conductivity by ∼20 times. These results pave new avenues for the exploitation of smart MOFs with multimode responsive behavior for on-demand sensing/detection applications.

Keystone microbes affect the evolution and ecological coexistence of the community via species/strain specificity.

AIM: Microbial communities exhibit different diversity and fluctuations in the ecological functions due to time and environmental migration. Despite a long history of research and a plethora of data, the factors determining the biodiversity and stability of ecosystems is still elusive. METHODS AND RESULTS: Here, the Chinese Xiaoqu fermentation system was used as a template to explore the mechanism in which the species specificity and strain in the initial phase affect the community structure and metabolites in the subsequent micro-ecosystem. The micro-ecosystem has been applied for hundreds of years, and the main metabolic function can be reproduced and traced. CONCLUSIONS: The result proved that Rhizopus spp. is a keystone microbe with a species/strain specificity affecting the trophic interaction niche and function of modules in the complex community through glucose. The fungal community was demonstrated to have a high sealing and stability, while the bacterial community was generally found to change the community structure, physiological function, and interaction relationship, producing strains with connector functions to adapt to fluctuations. SIGNIFICANCE AND IMPACT OF THE STUDY: This study shows that the taxonomic level of key microbial strains can be changed to affect the evolution of coexistence and functional realisation of the community.

Rivaroxaban versus Vitamin K Antagonists (warfarin) based on the triple therapy for left ventricular thrombus after ST-Elevation myocardial infarction.

BACKGROUND: Left ventricular thrombus (LVT) can complicate ST-Elevation myocardial infarction (STEMI) and is associated with poor outcomes. Conventional triple anticoagulation [Vitamin K Antagonists (VKA) plus dual-antiplatelet therapy (DAPT)] is the first-line therapy for LVT after STEMI. In patients with LVT following STEMI, contemporary data of triple therapy with rivaroxaban are lacking. METHODS: We conducted a retrospective cohort study involving 1335 STEMI patients who underwent primary percutaneous coronary intervention (PCI). Among patients who developed LVT after STEMI, we observed differences in efficacy between rivaroxaban plus DAPT therapy and VKA plus DAPT. The time of LVT resolution was also evaluated, as well as net clinical adverse events, and rates of bleeding events. RESULTS: In 1335 patients with STEMI, a total of 77 (5.7%) developed LVT over the follow-up period (median 25.0 months). Of the patients diagnosed with LVT, 31 patients were started on triple therapy with VKA, 33 patients on triple therapy with rivaroxaban. There was a consistent similarity in LVT resolution with rivaroxaban application compared to VKA application during the follow-up period [HR (log-rank test) 1.57(95% CI 0.89-2.77), p = 0.096; Adjusted HR 1.70(95% CI 0.90-3.22), p = 0.104]. Triple therapy with rivaroxaban showed quicker resolution than with VKA (6 months: p = 0.049; 12 months: p = 0.044; 18 months: p = 0.045). Similar risks of ISTH bleeding were not significantly different between the 2 groups [VKA 9.7% vs Rivaroxaban 6.1%, Adjusted HR 0.48 (95% CI 0.73-3.20); p = 0.444)]. Fewer net adverse clinical events (NACE) were observed in the rivaroxaban group [VKA 58.1% vs Rivaroxaban 24.2%; HR (log-rank test) 0.31(95% CI 0.14-0.68), p = 0.003; Adjusted HR 0.23(95% CI 0.09-0.57), p = 0.001]. CONCLUSION: In the observational study, triple therapy with rivaroxaban has similar and quicker LVT resolution in patients with LVT after STEMI, compared with triple therapy with VKA, and perhaps was associated with a better clinical benefit. Larger sample sizes and randomized controlled trials are needed to confirm this observation.

Structural basis for activity of TRIC counter-ion channels in calcium release.

Trimeric intracellular cation (TRIC) channels are thought to provide counter-ion currents that facilitate the active release of Ca2+ from intracellular stores. TRIC activity is controlled by voltage and Ca2+ modulation, but underlying mechanisms have remained unknown. Here we describe high-resolution crystal structures of vertebrate TRIC-A and TRIC-B channels, both in Ca2+-bound and Ca2+-free states, and we analyze conductance properties in structure-inspired mutagenesis experiments. The TRIC channels are symmetric trimers, wherein we find a pore in each protomer that is gated by a highly conserved lysine residue. In the resting state, Ca2+ binding at the luminal surface of TRIC-A, on its threefold axis, stabilizes lysine blockage of the pores. During active Ca2+ release, luminal Ca2+ depletion removes inhibition to permit the lysine-bearing and voltage-sensing helix to move in response to consequent membrane hyperpolarization. Diacylglycerol is found at interprotomer interfaces, suggesting a role in metabolic control.

The Anthocyanin Accumulation Related ZmBZ1, Facilitates Seedling Salinity Stress Tolerance via ROS Scavenging.

Anthocyanins are a class of antioxidants that scavenge free radicals in cells and play an important role in promoting human health and preventing many diseases. Here, we characterized a maize Bronze gene (BZ1) from the purple colored W22 introgression line, which encodes an anthocyanin 3-O-glucosyltransferase, a key enzyme in the anthocyanin synthesis pathway. Mutation of ZmBZ1 showed bronze-colored seeds and reduced anthocyanins in seeds aleurone layer, seedlings coleoptile, and stem of mature plants by comparison with purple colored W22 (WT). Furthermore, we proved that maize BZ1 is an aleurone layer-specific expressed protein and sub-located in cell nucleus. Real-time tracing of the anthocyanins in developing seeds demonstrated that the pigment was visible from 16 DAP (day after pollination) in field condition, and first deposited in the crown part then spread all over the seed. Additionally, it was transferred along with the embryo cell activity during seed germination, from aleurone layer to cotyledon and coleoptile, as confirmed by microscopy and real-time qRT-PCR. Finally, we demonstrated that the ZmBZ1 contributes to stress tolerance, especially salinity. Further study proved that ZmBZ1 participates in reactive oxygen scavenging (ROS) by accumulating anthocyanins, thereby enhancing the tolerance to abiotic stress.

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding.

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

Natural variations in the P-type ATPase heavy metal transporter gene ZmHMA3 control cadmium accumulation in maize grains.

Cadmium (Cd) accumulation in maize grains is detrimental to human health. Developing maize varieties with low Cd content is important for safe consumption of maize grains. However, the key genes controlling maize grain Cd accumulation have not been cloned. Here, we identified one major locus for maize grain Cd accumulation (qCd1) using a genome-wide association study (GWAS) and bulked segregant RNA-seq analysis with a biparental segregating population of Jing724 (low-Cd line) and Mo17 (high-Cd line). The candidate gene ZmHMA3 was identified by fine mapping and encodes a tonoplast-localized heavy metal P-type ATPase transporter. An ethyl methane sulfonate mutant analysis and an allelism test confirmed that ZmHMA3 influences maize grain Cd accumulation. A transposon in intron 1 of ZmHMA3 is responsible for the abnormal amino acid sequence in Mo17. Based on the natural sequence variations in the ZmHMA3 gene of diverse maize lines, four PCR-based molecular markers were developed, and these were successfully used to distinguish five haplotypes with different grain Cd contents in the GWAS panel and to predict grain Cd contents of widely used maize inbred lines and hybrids. These molecular markers can be used to breed elite maize varieties with low grain Cd contents.

Genome-wide analyses reveal footprints of divergent selection and popping-related traits in CIMMYT's maize inbred lines.

Popcorn (Zea mays L. var. Everta) is the most ancient type of cultivated maize. However, there is little known about the genetics of popping-related traits based on genotyping-by-sequencing (GBS) technology. Here, we characterized the phenotypic variation for seven popping-related traits in maize kernels among 526 CIMMYT inbred lines (CMLs). In total, 155 083 high-quality single nucleotide polymorphism (SNP) markers were identified by a GBS approach. Several trait-associated loci were detected by genome-wide association study for color, popping expansion volume, shape, pericarp, flotation index, floury/vitreous, and protein content, explaining a majority of the observed phenotypic variance, and these were validated by a diverse panel comprising 764 tropical landrace accessions. Sixty two of the identified loci were recognized to have undergone selection. On average, there was a 55.27% frequency for alleles that promote popping in CMLs. Our work not only pinpoints previously unknown loci for popping-related traits, but also reveals that many of these loci have undergone selection. Beyond establishing a new benchmark for the genetics of popcorn, our study provides a foundation for gene discovery and breeding. It also presents evidence to investigate the role of a gradual loss of popping ability as a by-product of diversification of culinary uses throughout the evolution of teosinte-to-modern maize.

LRSCnet: Local Reference Semantic Code learning for breast tumor classification in ultrasound images.

PURPOSE: This study proposed a novel Local Reference Semantic Code (LRSC) network for automatic breast ultrasound image classification with few labeled data. METHODS: In the proposed network, the local structure extractor is firstly developed to learn the local reference which describes common local characteristics of tumors. After that, a two-stage hierarchical encoder is developed to encode the local structures of lesion into the high-level semantic code. Based on the learned semantic code, the self-matching layer is proposed for the final classification. RESULTS: In the experiment, the proposed method outperformed traditional classification methods and AUC (Area Under Curve), ACC (Accuracy), Sen (Sensitivity), Spec (Specificity), PPV (Positive Predictive Values), and NPV(Negative Predictive Values) are 0.9540, 0.9776, 0.9629, 0.93, 0.9774 and 0.9090, respectively. In addition, the proposed method also improved matching speed. CONCLUSIONS: LRSC-network is proposed for breast ultrasound images classification with few labeled data. In the proposed network, a two-stage hierarchical encoder is introduced to learn high-level semantic code. The learned code contains more effective high-level classification information and is simpler, leading to better generalization ability.

Plasma IP-10 and MCP-3 levels are highly associated with disease severity and predict the progression of COVID-19.

BACKGROUND: The outbreak of coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 was first reported in Wuhan, December 2019, and continuously poses a serious threat to public health, highlighting the urgent need of identifying biomarkers for disease severity and progression. OBJECTIVE: We sought to identify biomarkers for disease severity and progression of COVID-19. METHODS: Forty-eight cytokines in the plasma samples from 50 COVID-19 cases including 11 critically ill, 25 severe, and 14 moderate patients were measured and analyzed in combination with clinical data. RESULTS: Levels of 14 cytokines were found to be significantly elevated in COVID-19 cases and showed different expression profiles in patients with different disease severity. Moreover, expression levels of IFN-γ-induced protein 10, monocyte chemotactic protein-3, hepatocyte growth factor, monokine-induced gamma IFN, and macrophage inflammatory protein 1 alpha, which were shown to be highly associated with disease severity during disease progression, were remarkably higher in critically ill patients, followed by severe and then the moderate patients. Serial detection of the 5 cytokines in 16 cases showed that continuously high levels were associated with deteriorated progression of disease and fatal outcome. Furthermore, IFN-γ-induced protein 10 and monocyte chemotactic protein-3 were excellent predictors for the progression of COVID-19, and the combination of the 2 cytokines showed the biggest area under the curve of the receiver-operating characteristics calculations with a value of 0.99. CONCLUSIONS: In this study, we report biomarkers that are highly associated with disease severity and progression of COVID-19. These findings add to our understanding of the immunopathologic mechanisms of severe acute respiratory syndrome coronavirus 2 infection, and provide potential therapeutic targets and strategies.