Integrative transcriptomics analysis and experimental validation reveal immunomodulatory patterns in keratoconus.

Keratoconus is a progressive disorder of the cornea and is typically considered a noninflammatory disease. However, increasing evidence indicates that immune disorders play an essential role in keratoconus progression, but the immune-related etiology remains elusive. Here, we comprehensively utilized bioinformatics approaches and experimental methods to explore the potential immunoregulatory mechanism of keratoconus progression. Transcriptomics data containing two keratoconus patient groups was derived from the public dataset GSE151631. The intersection of genes and known immunological genes was used to obtain differentially expressed immune-related genes. We utilized various protein clustering algorithms to screen out and validated the hub immune-related genes, and further explored their potential biological functions via gene annotation and pathway enrichment analyses. Moreover, the underlying immune landscape and drug targets were predicted by immune cell infiltration analysis and drug-gene interaction analysis. Furthermore, keratoconus-related immunoregulatory competitive endogenous RNA networks were constructed and experimentally validated. After filtering and experimental validation, nine keratoconus-associated immune-related genes were credible. Infiltrated monocytes might play an essential role in the progression of keratoconus. Moreover, eleven intersecting drugs targeting four genes, CCR2, CCR5, F2RL1, and ADORA1, were considered as potential druggable molecular targets for keratoconus. Furthermore, in the competitive endogenous RNA network, we identified several lncRNAs and miRNAs as critical noncoding RNAs regulating the hub genes. Overall, our data indicated that the immunomodulatory patterns had undergone changes in the pathogenesis of keratoconus, which might facilitate the understanding of keratoconus-related immune processes and provide novel insights into developing new immunotherapies for keratoconus.

Genome-Wide Identification of GmSPS Gene Family in Soybean and Expression Analysis in Response to Cold Stress.

Sucrose metabolism plays a critical role in development, stress response, and yield formation of plants. Sucrose phosphate synthase (SPS) is the key rate-limiting enzyme in the sucrose synthesis pathway. To date, genome-wide survey and comprehensive analysis of the SPS gene family in soybean (Glycine max) have yet to be performed. In this study, seven genes encoding SPS were identified in soybean genome. The structural characteristics, phylogenetics, tissue expression patterns, and cold stress response of these GmSPSs were investigated. A comparative phylogenetic analysis of SPS proteins in soybean, Medicago truncatula, Medicago sativa, Lotus japonicus, Arabidopsis, and rice revealed four families. GmSPSs were clustered into three families from A to C, and have undergone five segmental duplication events under purifying selection. All GmSPS genes had various expression patterns in different tissues, and family A members GmSPS13/17 were highly expressed in nodules. Remarkably, all GmSPS promoters contain multiple low-temperature-responsive elements such as potential binding sites of inducer of CBF expression 1 (ICE1), the central regulator in cold response. qRT-PCR proved that these GmSPS genes, especially GmSPS8/18, were induced by cold treatment in soybean leaves, and the expression pattern of GmICE1 under cold treatment was similar to that of GmSPS8/18. Further transient expression analysis in Nicotiana benthamiana and electrophoretic mobility shift assay (EMSA) indicated that GmSPS8 and GmSPS18 transcriptions were directly activated by GmICE1. Taken together, our findings may aid in future efforts to clarify the potential roles of GmSPS genes in response to cold stress in soybean.

RNA-Seq and Comparative Transcriptomic Analyses of Asian Soybean Rust Resistant and Susceptible Soybean Genotypes Provide Insights into Identifying Disease Resistance Genes.

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is one of the most destructive foliar diseases that affect soybeans. Developing resistant cultivars is the most cost-effective, environmentally friendly, and easy strategy for controlling the disease. However, the current understanding of the mechanisms underlying soybean resistance to P. pachyrhizi remains limited, which poses a significant challenge in devising effective control strategies. In this study, comparative transcriptomic profiling using one resistant genotype and one susceptible genotype was performed under infected and control conditions to understand the regulatory network operating between soybean and P. pachyrhizi. RNA-Seq analysis identified a total of 6540 differentially expressed genes (DEGs), which were shared by all four genotypes. The DEGs are involved in defense responses, stress responses, stimulus responses, flavonoid metabolism, and biosynthesis after infection with P. pachyrhizi. A total of 25,377 genes were divided into 33 modules using weighted gene co-expression network analysis (WGCNA). Two modules were significantly associated with pathogen defense. The DEGs were mainly enriched in RNA processing, plant-type hypersensitive response, negative regulation of cell growth, and a programmed cell death process. In conclusion, these results will provide an important resource for mining resistant genes to P. pachyrhizi infection and valuable resources to potentially pyramid quantitative resistance loci for improving soybean germplasm.

Tsp-1 is involved in DNA stability through Tgf-ß1 activation domain in cone photoreceptor 661 W cells.

Thrombospondin-1 (Tsp-1), a matricellular protein, could protect retinal neurons from endogenous or exogenous insults; however, its underlying mechanism remains unclear. Thus, this study aimed to investigate Tsp-1-mediated neuron-protection effect in retinal cells. Our data showed that Tsp-1 downregulation would aggravate UV irradiation-induced DNA damage in 661 W cells and cone photoreceptor cells. The increasing levels of poly (ADP ribose) polymer (PAR) and γ-H2AX in Tsp-1-silenced 661 W cells indicate severe DNA single-strand breaks (SSBs) and double-strand breaks (DSBs). By utilizing an error-prone substrate, Tsp-1 silencing significantly increased deleted DNA end joining in 661 W cells with spontaneous DNA damage (SDD). Moreover, Tsp-1 is indirectly involved in DNA stability in 661 W cells as UV treatment caused a significant Tsp-1 decreasing in cytoplasm, but no obvious Tsp-1 alteration in cell nuclear of 661 W cells. Furthermore, our data indicate that Tgf-ß1 activation domain in Tsp-1 plays a critical role in DNA stability in 661 W cells through expressing mutated exogenous Tsp-1 and Tgf-ß inhibitor, LSKL. Therefore, this study provides new insights into the mechanism of the neuroprotective action positively mediated by Tsp-1, which might be a therapeutic target for the treatment of retinal pathology.

Kruppel-like factor 2 acts as a tumor suppressor in human retinoblastoma.

Krüppel-like factor 2 (KLF2) belongs to the KLF family of zinc-finger transcription factors and mediates the occurrence and progression of various cancers. However, little is known about its expression pattern and biological role in retinoblastoma (RB). In the present study, we showed that KLF2 was markedly downregulated in human RB tissue compared with retina. KLF2 overexpression significantly inhibited RB cell proliferation and decreased proliferating cell nuclear antigen (PCNA) expression. Subsequently, we confirmed that KLF2 arrested cells at the G1-S phase transition, accompanied by the upregulation of p21 and downregulation of CyclinD1, as well as the activation of mitochondria-mediated apoptosis in RB cells. In addition, KLF2 overexpression contributed to suppressing RB cell migration and invasion by downregulating matrix metallopeptidase 9 (MMP9). On the contrary, KLF2 downregulation promoted RB cells proliferation, migration and invasion. Notably, the KLF2 expression pattern was opposite to that of C-X-C chemokine receptor 4 (CXCR4) in the two RB cell lines, KLF2 overexpression significantly decreased CXCR4 expression, silencing KLF2 had the opposite effect. Furthermore, dual-luciferase reporter and chromatin immunoprecipitation (ChIP) assays confirmed that KLF2 directly bound to the CXCR4 promoter and negatively regulated its expression in RB cells. Collectively, our results suggested that KLF2 function as a tumor suppressor in RB and may represent a potential therapeutic target for RB.

Comparison of DNA stability and its related genes of neurons derived from induced pluripotent stem cells and primary retinal neurons.

Maintaining DNA stability in induced pluripotent stem cells (iPSCs) and iPSCs-derived neurons is a challenge in their clinical application. In the present study, we compared DNA stability between primary retinal neurons and differentiated neurons. We found that the basal level of γ-H2AX phosphorylation, a specific marker of DNA breaks, was notably higher (~26-folds) in human iPSCs compared to iPSCs-derived neurons. However, iPSCs-derived neurons are more sensitive to UV treatment compared to primary rat retinal neurons (postnatal Day 1). UV treatment induced a significantly decreasing in the cell viability of iPSCs-derived neurons by ~76.1%, whereas ~20.8% in primary retinal neurons. After analyzing the expression levels of genes involved in DNA stability, such as Brca1, Ligase IV, Ku80, and Mre11, we found that Ku80 and its heterodimeric partner, Ku70 were positive in iPSCs-derived neurons. However, both Ku80 and Ku70 are not expressed in primary retinal neurons and cerebellar neurons. Similarly, both Ku80 and Ku70 are also expressed in 3D retinal organoids from human embryonic stem cells (ESCs), except for a few Map2-negative cells and the hyaloid vessels of mice E12.5 retinas. Hence, Ku80, and Ku70 are specifically expressed in stem cell-derived neurons. Moreover, using the Ku80 inhibitor Compound L, our data showed that Ku80 promotes the DNA stability and cell viability of iPSCs-derived neurons. Thus, our results demonstrated that iPSCs-, ESCs-derived neurons have specific characteristics of DNA stability. This study provides new insights into the neural differentiation of stem cells but might also warrant the future clinical application of stem cells in neurodegenerative diseases.

Comparison of the Response to the CXCR4 Antagonist AMD3100 during the Development of Retinal Organoids Derived from ES Cells and Zebrafish Retina.

Retinal organoids generated from human embryonic stem cells or iPSCs recreate the key structural and functional features of mammalian retinal tissue in vitro. However, the differences in the development of retinal organoids and normal retina in vivo are not well defined. Thus, in the present study, we analyzed the development of retinal organoids and zebrafish retina after inhibition of CXCR4, a key role in neurogenesis and optic nerve development, with the antagonist AMD3100. Our data indicated that CXCR4 was mainly expressed in ganglion cells in retinal organoids and was rarely expressed in amacrine or photoreceptor cells. AMD3100 treatment reduced the retinal organoid generation ratio, impaired differentiation, and induced morphological changes. Ganglion cells, amacrine cells, and photoreceptors were decreased and abnormal locations were observed in organoids treated with AMD3100. Neuronal axon outgrowth was also damaged in retinal organoids. Similarly, a decrease of ganglion cells, amacrine cells, and photoreceptors and the distribution of neural outgrowth was induced by AMD3100 treatment in zebrafish retina. However, abnormal photoreceptor ensembles induced by AMD3100 treatment in the organoids were not detected in zebrafish retina. Therefore, our study suggests that although retinal organoids might provide a reliable model for reproducing a retinal developmental model, there is a difference between the organoids and the retina in vivo.

Gata3 Silencing Is Involved in Neuronal Differentiation and Its Abnormal Expression Impedes Neural Activity in Adult Retinal Neurocytes.

GATA binding protein 3 (Gata3), a zinc-finger transcription factor, plays an important role in neural development. However, its expression and bioactivity in the retina remain unclear. In the present study, our data indicated that Gata3 maintains the precursor state of 661W cells, and Gata3 silencing induces cell differentiation. The expression of Nestin, a marker of precursor cells, was significantly decreased in parallel, whereas the expression of Map2, a marker of differentiated neurons, was significantly increased following the decrease in Gata3. Neurite outgrowth was increased by 2.78-fold in Gata3-silenced cells. Moreover, Gata3 expression generally paralleled that of Nestin in developing mouse retinas. Both Gata3 and Nestin were expressed in the retina at postnatal day 1 and silenced in the adult mouse retina. Exogenous Gata3 significantly inhibited the neural activity of primary retinal neurocytes (postnatal day 1) by decreasing synaptophysin levels, neurite outgrowth, and cell viability. Furthermore, in vivo, exogenous Gata3 significantly induced apoptosis and the contraction of retinal outlay filaments and decreased the a- and b-waves in adult mouse intravitreal injected with AAV-Re-Gata3-T2A-GFP. Thus, Gata3 silencing promotes neuronal differentiation and neurite outgrowth. Its abnormal expression impedes neural activity in adult retinal neurocytes. This study provides new insights into Gata3 bioactivity in retinal neurocytes.

Brca1 Is Regulated by the Transcription Factor Gata3, and Its Silencing Promotes Neural Differentiation in Retinal Neurons.

Previous studies have indicated that Brca1 (Breast cancer suppressor gene 1) plays an important role in neural development and degenerative diseases. However, the bioactivity and regulatory mechanism of Brca1 expression in retinal neurocytes remain unclear. In the present study, our data indicated that Brca1 maintains the state of neuronal precursor cells. Brca1 silencing induces differentiation in 661W cells. Nestin, a marker of precursor cells, was significantly decreased in parallel with Brca1 silencing in 661W cells, whereas Map2 (Microtubule associated protein 2), a marker of differentiated neurons, was significantly increased. Neurite outgrowth was increased by ~4.0-fold in Brca1-silenced cells. Moreover, DNA affinity purification assays and ChIP assays demonstrated that Gata3 (GATA binding protein 3) regulates Brca1 transcription in 661W cells. Silencing or overexpressing Gata3 could significantly regulate the expression of Brca1 and affect its promoter inducibility. Furthermore, the expression of Gata3 generally occurred in parallel with that of Brca1 in developing mouse retinas. Both Gata3 and Brca1 are expressed in the neonatal mouse retina but are developmentally silenced with age. Exogenous Gata3 significantly inhibited neural activity by decreasing synaptophysin and neurite outgrowth. Thus, this study demonstrated that Brca1 is transcriptionally regulated by Gata3. Brca1/Gata3 silencing is involved in neuronal differentiation and maturation.

Identification of Quantitative Trait Locus and Candidate Genes for Drought Tolerance in a Soybean Recombinant Inbred Line Population.

With global warming and regional decreases in precipitation, drought has become a problem worldwide. As the number of arid regions in the world is increasing, drought has become a major factor leading to significant crop yield reductions and food crises. Soybean is a crop that is relatively sensitive to drought. It is also a crop that requires more water during growth and development. The aim of this study was to identify the quantitative trait locus (QTL) that affects drought tolerance in soybean by using a recombinant inbred line (RIL) population from a cross between the drought-tolerant cultivar 'Jindou21' and the drought-sensitive cultivar 'Zhongdou33'. Nine agronomic and physiological traits were identified under drought and well-watered conditions. Genetic maps were constructed with 923,420 polymorphic single nucleotide polymorphism (SNP) markers distributed on 20 chromosomes at an average genetic distance of 0.57 centimorgan (cM) between markers. A total of five QTLs with a logarithm of odds (LOD) value of 4.035-8.681 were identified on five chromosomes. Under well-watered conditions and drought-stress conditions, one QTL related to the main stem node number was located on chromosome 16, accounting for 17.177% of the phenotypic variation. Nine candidate genes for drought resistance were screened from this QTL, namely Glyma.16G036700, Glyma.16G036400, Glyma.16G036600, Glyma.16G036800, Glyma.13G312700, Glyma.13G312800, Glyma.16G042900, Glyma.16G043200, and Glyma.15G100700. These genes were annotated as NAC transport factor, GATA transport factor, and BTB/POZ-MATH proteins. This result can be used for molecular marker-assisted selection and provide a reference for breeding for drought tolerance in soybean.

Overexpression of GmMYB14 improves high-density yield and drought tolerance of soybean through regulating plant architecture mediated by the brassinosteroid pathway.

MYB transcription factors (TFs) have been reported to regulate the biosynthesis of secondary metabolites, as well as to mediate plant adaption to abiotic stresses, including drought. However, the roles of MYB TFs in regulating plant architecture and yield potential remain poorly understood. Here, we studied the roles of the dehydration-inducible GmMYB14 gene in regulating plant architecture, high-density yield and drought tolerance through the brassinosteroid (BR) pathway in soybean. GmMYB14 was shown to localize to nucleus and has a transactivation activity. Stable GmMYB14-overexpressing (GmMYB14-OX) transgenic soybean plants displayed a semi-dwarfism and compact plant architecture associated with decreased cell size, resulting in a decrease in plant height, internode length, leaf area, leaf petiole length and leaf petiole angle, and improved yield in high density under field conditions. Results of the transcriptome sequencing suggested the involvement of BRs in regulating GmMYB14-OX plant architecture. Indeed, GmMYB14-OX plants showed reduced endogenous BR contents, while exogenous application of brassinolide could partly rescue the phenotype of GmMYB14-OX plants. Furthermore, GmMYB14 was shown to directly bind to the promoter of GmBEN1 and up-regulate its expression, leading to reduced BR content in GmMYB14-OX plants. GmMYB14-OX plants also displayed improved drought tolerance under field conditions. GmBEN1 expression was also up-regulated in the leaves of GmMYB14-OX plants under polyethylene glycol treatment, indicating that the GmBEN1-mediated reduction in BR level under stress also contributed to drought/osmotic stress tolerance of the transgenic plants. Our findings provided a strategy for stably increasing high-density yield and drought tolerance in soybean using a single TF-encoding gene.

Genome-Wide Investigation and Expression Profiling Under Abiotic Stresses of a Soybean Unknown Function (DUF21) and Cystathionine-ß-Synthase (CBS) Domain-Containing Protein Family.

Cystathionine-ß-synthase (CBS) domain-containing proteins (CDCPs) constitute a large family in plants, and members of this family have been implicated in a variety of biological processes. However, the precise functions and the underlying mechanisms of most members of this family in plants remain to be elucidated. CBSDUF proteins belong to the CDCP superfamily, which contains one domain of unknown function (DUF21) and an N terminus that is adjacent to two intracellular CBS domains. In this study, a comprehensive genome database analysis of soybean was performed to investigate the role(s) of these CBSDUFs and to explore their nomenclature, classification, chromosomal distribution, exon-intron organization, protein structure, and phylogenetic relationships; the analysis identified a total of 18 putative CBSDUF genes. Using specific protein domains and phylogenetic analysis, the CBSDUF gene family was subdivided into eight groups. The soybean CBSDUF genes showed an uneven distribution on 12 chromosomes of Glycine max. RNA-seq transcriptome data from different tissues in public databases revealed tissue-specific and differential expression profiles of the GmCBSDUFs, and qPCR analysis revealed that certain groups of soybean CBSDUFs are likely involved in specific stress responses. In addition, GmCBSDUF3 transgenic Arabidopsis was subjected to phenotypic analysis under NaCl, PEG, and ABA stress treatments. The overexpression of GmCBSDUF3 could enhance tolerance to drought and salt stress in Arabidopsis. This study presents a first comprehensive look at soybean CBSDUF proteins and provides valuable resources for functionally elucidating this protein subgroup within the CBS domain-containing protein family.

Genome-Wide Identification and Characterization of Soybean GmLOR Gene Family and Expression Analysis in Response to Abiotic Stresses.

The LOR (LURP-one related) family genes encode proteins containing a conserved LOR domain. Several members of the LOR family genes are required for defense against Hyaloperonospora parasitica (Hpa) in Arabidopsis. However, there are few reports of LOR genes in response to abiotic stresses in plants. In this study, a genome-wide survey and expression levels in response to abiotic stresses of 36 LOR genes from Glycine max were conducted. The results indicated that the GmLOR gene family was divided into eight subgroups, distributed on 14 chromosomes. A majority of members contained three extremely conservative motifs. There were four pairs of tandem duplicated GmLORs and nineteen pairs of segmental duplicated genes identified, which led to the expansion of the number of GmLOR genes. The expansion patterns of the GmLOR family were mainly segmental duplication. A heatmap of soybean LOR family genes showed that 36 GmLOR genes exhibited various expression patterns in different tissues. The cis-acting elements in promoter regions of GmLORs include abiotic stress-responsive elements, such as dehydration-responsive elements and drought-inducible elements. Real-time quantitative PCR was used to detect the expression level of GmLOR genes, and most of them were expressed in the leaf or root except that GmLOR6 was induced by osmotic and salt stresses. Moreover, GmLOR4/10/14/19 were significantly upregulated after PEG and salt treatments, indicating important roles in the improvement of plant tolerance to abiotic stress. Overall, our study provides a foundation for future investigations of GmLOR gene functions in soybean.

Identification of soybean drought-tolerant genotypes and loci correlated with agronomic traits contributes new candidate genes for breeding.

KEY MESSAGE: Drought tolerance level of 136 soybean genotypes, the correlations among traits were evaluated, and several important drought-tolerant genotypes, traits, SNPs and genes were possibly useful for soybean genetic breeding. Drought is an adverse environmental factor affecting crops growth, development, and yield. Promising genotypes and genes with improved tolerance to drought are probably effective ways to alleviate the situation. In this study, our main task was to determine drought tolerance level of 136 soybean genotypes, the correlations among physiological and agronomic traits under drought, and drought-tolerant single nucleotide polymorphism (SNPs) and genes. In this study, twenty-six varieties were identified as excellent tolerant genotypes to stress among which S14, S93 and S135 with high drought-tolerant index (DTI > 1.3) and yield (Y > 300 kg). Fourteen varieties were identified as drought-sensitive genotypes, such as S25, S45 and S58, with low drought-tolerant index (DTI < 0.5). 422 SNPs and 302 genes correlated with seed number per plant (SNPP), maturity (M), number of seeds per pod (NSPP), node number of main stem (NNMS), Stem diameter (SD) and pull stem (PS) were detected under well-watered and drought conditions by genome-wide association study (GWAS). Among them, we found SNPs (Chr 3:1758920-1958934) between drought-tolerant and sensitive genotypes. 13 genes (Glyma.03G017800, Glyma.03G018000, Glyma.03G018200, Glyma.03G018400, Glyma.03G018500, Glyma.03G018600, Glyma.03G018700, Glyma.03G018800, Glyma.03G018900, Glyma.03G019000, Glyma.03G019100, Glyma.03G019200, Glyma.03G019300) correlated with NNMS were detected. By qRT-PCR, the expression level of Glyma.03G018000 and Glyma.03G018900 in drought-tolerant varieties was significantly increased, but low or no expression in sensitive varieties under drought stress. This study provides important drought-tolerant genotypes, traits, SNPs and potential genes, possibly useful for soybean genetic breeding.

Genome-wide survey of soybean papain-like cysteine proteases and their expression analysis in root nodule symbiosis.

BACKGROUND: Plant papain-like cysteine proteases (PLCPs) are a large class of proteolytic enzymes and play important roles in root nodule symbiosis (RNS), while the whole-genome studies of PLCP family genes in legume are quite limited, and the roles of Glycine max PLCPs (GmPLCPs) in nodulation, nodule development and senescence are not fully understood. RESULTS: In the present study, we identified 97 GmPLCPs and performed a genome-wide survey to explore the expansion of soybean PLCP family genes and their relationships to RNS. Nineteen paralogous pairs of genomic segments, consisting of 77 GmPLCPs, formed by whole-genome duplication (WGD) events were identified, showing a high degree of complexity in duplication. Phylogenetic analysis among different species showed that the lineage differentiation of GmPLCPs occurred after family expansion, and large tandem repeat segment were specifically in soybean. The expression patterns of GmPLCPs in symbiosis-related tissues and nodules identified RNS-related GmPLCPs and provided insights into their putative symbiotic functions in soybean. The symbiotic function analyses showed that a RNS-related GmPLCP gene (Glyma.04G190700) really participate in nodulation and nodule development. CONCLUSIONS: Our findings improved our understanding of the functional diversity of legume PLCP family genes, and provided insights into the putative roles of the legume PLCPs in nodulation, nodule development and senescence.

CRISPR/Cas9-mediated targeted mutagenesis of GmSPL9 genes alters plant architecture in soybean.

BACKGROUND: The plant architecture has significant effects on grain yield of various crops, including soybean (Glycine max), but the knowledge on optimization of plant architecture in order to increase yield potential is still limited. Recently, CRISPR/Cas9 system has revolutionized genome editing, and has been widely utilized to edit the genomes of a diverse range of crop plants. RESULTS: In the present study, we employed the CRISPR/Cas9 system to mutate four genes encoding SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors of the SPL9 family in soybean. These four GmSPL9 genes are negatively regulated by GmmiR156b, a target for the improvement of soybean plant architecture and yields. The soybean Williams 82 was transformed with the binary CRISPR/Cas9 plasmid, assembled with four sgRNA expression cassettes driven by the Arabidopsis thaliana U3 or U6 promoter, targeting different sites of these four SPL9 genes via Agrobacterium tumefaciens-mediated transformation. A 1-bp deletion was detected in one target site of the GmSPL9a and one target site of the GmSPL9b, respectively, by DNA sequencing analysis of two T0-generation plants. T2-generation spl9a and spl9b homozygous single mutants exhibited no obvious phenotype changes; but the T2 double homozygous mutant spl9a/spl9b possessed shorter plastochron length. In T4 generation, higher-order mutant plants carrying various combinations of mutations showed increased node number on the main stem and branch number, consequently increased total node number per plants at different levels. In addition, the expression levels of the examined GmSPL9 genes were higher in the spl9b-1 single mutant than wild-type plants, which might suggest a feedback regulation on the expression of the investigated GmSPL9 genes in soybean. CONCLUSIONS: Our results showed that CRISPR/Cas9-mediated targeted mutagenesis of four GmSPL9 genes in different combinations altered plant architecture in soybean. The findings demonstrated that GmSPL9a, GmSPL9b, GmSPL9c and GmSPL9 function as redundant transcription factors in regulating plant architecture in soybean.

Estimating the transmissibility of hand, foot, and mouth disease by a dynamic model.

OBJECTIVE: Hand, foot, and mouth disease (HFMD) is a serious public health issue in many countries; however, its transmissibility in county-level outbreaks remains unclear. The aim of this study is to estimate the transmissibility of HFMD epidemics on both city level and county level, for a better understanding of the transmission dynamics of HFMD epidemics. STUDY DESIGN: Simulation based on data obtained from the China Information System for Disease Control and Prevention. METHODS: The weekly number of reported HFMD cases from April 2009 to December 2017 in nine regions of Changsha City was collected. A susceptible-infectious-recovered model was used to estimate the transmissibility of HFMD. The reproduction number of reported cases during the ascending (denoted as Rasc) and descending (denoted as Rdes) period was used to describe the transmissibility of HFMD. RESULTS: The Rasc and Rdes for HFMD in Changsha was 1.44 (95% confidence interval [CI]: 1.41-1.48) and 0.71 (95% CI: 0.69-0.73), respectively. There was no statistical significance of Rasc values among nine regions (F = 1.056, P = 0.396), nor of Rdes values among nine regions (F = 1.676, P = 0.106). The average Rasc (1.53, 95% CI: 1.46-1.61) from 2009 to 2012 was higher than the one (1.37, 95% CI: 1.34-1.40) from 2013 to 2017 (t = 3.974, P < 0.001), but the average Rdes (0.67, 95% CI: 0.63-0.70) from 2009 to 2012 was lower than the one (0.74, 95% CI: 0.73-0.76) from 2013 to 2017 (t = -3.751, P < 0.001). CONCLUSIONS: The epidemic of HFMD in Changsha City is still grim, and integrated strategies should be taken for controlling and preventing HFMD.

Rabies Virus Transmission in Solid Organ Transplantation, China, 2015-2016.

We report rabies virus transmission among solid organ transplantation recipients in Changsha, China, in 2016. Two recipients were confirmed to have rabies and died. Our findings suggest that more attention should be paid to the possibility of rabies virus transmission through organ transplantation for clinical and public health reasons.

Identification and Comparative Analysis of CBS Domain-Containing Proteins in Soybean (Glycine max) and the Primary Function of GmCBS21 in Enhanced Tolerance to Low Nitrogen Stress.

Nitrogen is an important macronutrient required for plant growth, and is a limiting factor for crop productivity. Improving the nitrogen use efficiency (NUE) is therefore crucial. At present, the NUE mechanism is unclear and information on the genes associated with NUE in soybeans is lacking. cystathionine beta synthase (CBS) domain-containing proteins (CDCPs) may be implicated in abiotic stress tolerance in plants. We identified and classified a CBS domain-containing protein superfamily in soybean. A candidate gene for NUE, GmCBS21, was identified. GmCBS21 gene characteristics, the temporal expression pattern of the GmCBS21 gene, and the phenotype of GmCBS21 overexpression in transgenic Arabidopsis thaliana under low nitrogen stress were analyzed. The phenotypes suggested that the transgenic Arabidopsis thaliana seedlings performed better under the nitrogen-deficient condition. GmCBS21-overexpressing transgenic plants exhibit higher low nitrogen stress tolerance than WT plants, and this suggests its role in low nitrogen stress tolerance in plants. We conclude that GmCBS21 may serve as an excellent candidate for breeding crops with enhanced NUE and better yield.