Construction of MicroRNA-Target Interaction Networks Based on MicroRNA Expression Profiles of HRV16-infected H1-HeLa Cells.

In the present study we investigated the changes in miRNA levels inhuman rhinovirus 16 (HRV16)-infected cells. A small RNA deep sequencing experiment was performed through next-generation sequencing. In total, 53 differentially expressed miRNAs were confirmed by RT-qPCR, including 37 known miRNAs and 16 novel miRNAs. Interaction networks between differentially expressed miRNAs and their targets were established by mirDIP and Navigator. The prediction results showed that QKI, NFAT5, BNC2, CELF2, LCOR, MBNL2, MTMR3, NFIB, PPARGC1A, RSBN1, TRPS1, WDR26, and ZNF148, which are associated with cellular differentiation and transcriptional regulation, were recognized by 12, 11, or 9 miRNAs. Many correlations were observed between transcriptional or post-transcriptional regulation of an miRNA and the expression levels of its target genes in HRV16-infected H1-HeLa cells.

Comparative Transcriptome and Metabolome Analysis of Resistant and Susceptible Piper Species Upon Infection by the Oomycete Phytophthora Capsici.

Phytophthora capsici is a destructive oomycete pathogen that causes devastating disease in black pepper, resulting in a significant decline in yield and economic losses. Piper nigrum (black pepper) is documented as susceptible to P. capsici, whereas its close relative Piper flaviflorum is known to be resistant. However, the molecular mechanism underlying the resistance of P. flaviflorum remains obscure. In this study, we conducted a comparative transcriptome and metabolome analysis between P. flaviflorum and P. nigrum upon P. capsici infection and found substantial differences in their gene expression profiles, with altered genes being significantly enriched in terms relating to plant-pathogen interaction, phytohormone signal transduction, and secondary metabolic pathways, including phenylpropanoid biosynthesis. Further metabolome analysis revealed the resistant P. flaviflorum to have a high background endogenous ABA reservoir and time-course-dependent accumulation of ABA and SA upon P. capsici inoculation, while the susceptible P. nigrum had a high background endogenous IAA reservoir and time-course-dependent accumulation of JA-Ile, the active form of JA. Investigation of the phenylpropanoid biosynthesis metabolome further indicated the resistant P. flaviflorum to have more accumulation of lignin precursors than the susceptible P. nigrum, resulting in a higher accumulation after inoculation. This study provides an overall characterization of biologically important pathways underlying the resistance of P. flaviflorum, which theoretically explains the advantage of using this species as rootstock for the management of oomycete pathogen in black pepper production.

Exploration of IRES Elements within the ORF of the Coxsackievirus B3 Genome.

Objective: This study aimed to identify internal ribosome entry sites (IRESs) in the open reading frame (ORF) of the Coxsackievirus B3 (CVB3) genome. Methods: The sequences of P1, P2, or P3 of the CVB3 genome or the truncated sequences from each antithymocyte globulin (ATG) to the end of the P1, P2, or P3 gene were inserted into the pEGFP-N1 vector. After transfection, possible IRES-dependent green fluorescent protein (GFP)-fused proteins were detected by anti-GFP western blotting. The sequences of possible IRESs were inserted into specific Fluc/Rluc bicistronic vectors, in which the potential IRESs were determined according to the Fluc/Rluc activity ratio. Expression of Fluc and Rluc mRNA of the bicistronic vector was detected by RT-qPCR. Results: After transfection of full length or truncated sequences of the P1, P2, or P3 plasmids, six GFP-fused protein bands in P1, six bands in P2 and nine bands in P3 were detected through western blotting. Two IRESs in VP2 (1461-1646 nt) and VP1 (2784-2983 nt) of P1; one IRES in 2C (4119-4564 nt) of P2; and two IRESs in 3C (5634-5834 nt) and 3D (6870-7087 nt) of P3 were identified according to Fluc/Rluc activity ratio. The cryptic promoter was also excluded by RT-qPCR. Conclusion: Five IRESs are present in the CVB3 coding region.

Genotype and mutation patterns of macrolide resistance genes of Mycoplasma pneumoniae from children with pneumonia in Qingdao, China, in 2019.

OBJECTIVES: This study assessed the incidence and resistance of Mycoplasma pneumoniae (MP) in children in Qingdao, China, in 2019. METHODS: We detected MP infection in 78 pharyngeal swabs from children with pneumonia by qPCR. The RepMP4 element in the P1 adhesin gene, domain V of the 23S rRNA gene, and the L4/L22 ribosomal proteins were amplified by nested PCR. Evolutionary analysis was conducted based on the P1 gene sequence. Resistance mutations in domain V of the 23S rRNA gene and L4/L22 ribosomal proteins were analysed. RESULTS: The incidence of MP infection in children with pneumonia was 59.0% (46/78). The mean duration of MP infection was longer than that of non-MP infection. According to P1 gene sequencing of 21 samples, 12 (57.1%) were type 1 and 9 (42.9%) were type 2. Drug resistance mutations A2063G in domain V of 23S rRNA gene and T508C in L22 were identified from all sequenced MP. However, mutations at positions 2064 and 2617 were not found in this study. C162A mutation appeared in most type 2 samples. A430G mutation appeared in one type 1 sample and in several type 2 samples. T279C mutation in L22 was mostly found in type 2 samples. CONCLUSION: The incidence of MP infection was 59.0% in children with pneumonia in Qingdao in 2019. Type 1 MP infection was slightly more common than type 2, indicating that the genotype of MP is gradually shifting from type 1 to type 2. Macrolide resistance mutation A2063G could be detected in all sequenced MP.

Differential expression of the Hsf family in cassava under biotic and abiotic stresses.

Heat shock transcription factors (Hsfs) are important regulators of biotic and abiotic stress responses in plants. Currently, the Hsf gene family is not well understood in cassava, an important tropical crop. In the present study, 32 MeHsf genes were identified from the cassava genome database, which were divided into three types based on functional domain and motif distribution analyses. Analysis of the differential expression of the genes belonging to the Hsf family in cassava was carried out based on published cassava transcriptome data from tissues/organs (leaf blade, leaf midvein, lateral buds, organized embryogenic structures, friable embryogenic callus, fibrous roots, storage roots, stem, petiole, shoot apical meristem, and root apical meristem) under abiotic stress (cold, drought) or biotic stress (mealybugs. cassava brown streak disease, cassava bacterial blight). The results show the expression diversity of cassava Hsfs genes in various tissues/organs. The transcript levels of MeHsfB3a, MeHsfA6a, MeHsfA2a, and MeHsfA9b were upregulated by abiotic and biotic stresses, such as cold, drought, cassava bacterial blight, cassava brown streak disease, and mealybugs, indicating their potential roles in mediating the response of cassava plants to environment stresses. Further interaction network and co-expression analyses suggests that Hsf genes may interact with Hsp70 family members to resist environmental stresses in cassava. These results provide valuable information for future studies of the functional characterization of the MeHsf gene family.

Construction of a high-density linkage map and QTL mapping for important agronomic traits in Stylosanthes guianensis (Aubl.) Sw.

Stylosanthes guianensis (Aubl.) Sw. is an economically important pasture and forage legume in tropical regions of the world. Genetic improvement of the crop can be enhanced through marker-assisted breeding. However, neither single nucleotide polymorphism (SNP) markers nor SNP-based genetic linkage map has been previously reported. In this study, a high-quality genetic linkage map of 2572 SNP markers for S. guianensis is generated using amplified-fragment single nucleotide polymorphism and methylation (AFSM) approach. The genetic map has 10 linkage groups (LGs), which spanned 2226.6 cM, with an average genetic distance of 0.87 cM between adjacent markers. Genetic mapping of quantitative trait loci (QTLs) for important agronomic traits such as yield-related and nutritional or quality-related traits was performed using F2 progeny of a cross between a male-sterile female parent TPRC1979 and male parent TPRCR273 with contrasting phenotypes for morphological and physiological traits. A total of 30 QTLs for 8 yield-related traits and 18 QTLs for 4 nutritional or quality-related traits are mapped on the linkage map. Both the high-quality genetic linkage map and the QTL mapping for important agronomic traits described here will provide valuable genetic resources for marker-assisted selection for S. guianensis.

Modeling the Transmission of Middle East Respirator Syndrome Corona Virus in the Republic of Korea.

The 2015 epidemic of Middle East respiratory syndrome (MERS) in the Republic of Korea has been the largest outbreak outside Middle East. This epidemic had caused 185 laboratory-confirmed cases and 36 deaths in the Republic of Korea until September 2, 2015, which attracted public's attention. Based on the detailed data of patients released by World Health Organization (WHO) and actual propagation of the epidemic, we construct two dynamical models to simulate the propagation processes from May 20 to June 8 and from June 9 to July 10, 2015, respectively and find that the basic reproduction number R0 reaches up to 4.422. The numerical analysis shows that the reasons of the outbreak spread quickly are lack of self-protection sense and targeted control measures. Through partial correction analysis, the parameters ß1 and γ have strong correlations with R0, i.e., the infectivity and proportion of the asymptomatic infected cases have much influence on the spread of disease. By sensitivity analysis, strengthening self-protection ability of susceptible and quickly isolating or monitoring close contacts are effective measures to control the disease.

Modeling the transmission dynamics of Ebola virus disease in Liberia.

Ebola virus disease (EVD) has erupted many times in some zones since it was first found in 1976. The 2014 EVD outbreak in West Africa is the largest ever, which has caused a large number of deaths and the most serious country is Liberia during the outbreak period. Based on the data released by World Health Organization and the actual transmission situations, we investigate the impact of different transmission routes on the EVD outbreak in Liberia and estimate the basic reproduction number R0 = 2.012 in the absence of effective control measures. Through sensitivity and uncertainty analysis, we reveal that the transmission coefficients of suspected and probable cases have stronger correlations on the basic reproduction number. Furthermore, we study the influence of control measures (isolation and safe burial measures) on EVD outbreak. It is found that if combined control measures are taken, the basic reproduction number will be less than one and thus EVD in Liberia may be well contained. The obtained results may provide new guidance to prevent and control the spread of disease.

Transcriptome sequencing and analysis of rubber tree (Hevea brasiliensis Muell.) to discover putative genes associated with tapping panel dryness (TPD).

BACKGROUND: Tapping panel dryness (TPD) involves in the partial or complete cessation of latex flow thus seriously affect latex production in rubber tree (Hevea brasiliensis). Numerous studies have been conducted to define its origin and nature, but the molecular nature and mechanism of TPD occurrence remains unknown. This study is committed to de novo sequencing and comparative analysis of the transcriptomes of healthy (H) and TPD-affected (T) rubber trees to identify the genes and pathways related to the TPD. RESULTS: Total raw reads of 34,632,012 and 35,913,020 bp were obtained from H and T library, respectively using Illumina Hiseq 2000 sequencing technology. De novo assemblies yielded 141,456 and 169,285 contigs, and 96,070 and 112,243 unigenes from H and T library, respectively. Among 73597 genes, 22577 genes were identified as differential expressed genes between H and T library via comparative transcript profiling. A majority of genes involved in natural rubber biosynthesis and jasmonate synthesis with most potential relevance in TPD occurrence were found to be differentially expressed. CONCLUSIONS: In TPD-affected trees, the expression of most genes related to the latex biosynthesis and jasmonate synthesis was severely inhibited and is probably the direct cause of the TPD. These new de novo transcriptome data sets provide a significant resource for the discovery of genes related to TPD and improve our understanding of the occurrence and maintainace of TPD.

Pinus taeda phenylpropenal double-bond reductase: purification, cDNA cloning, heterologous expression in Escherichia coli, and subcellular localization in P. taeda.

A phenylpropenal double-bond reductase (PPDBR) was obtained from cell suspension cultures of loblolly pine (Pinus taeda L.). Following trypsin digestion and amino acid sequencing, the cDNA encoding this protein was subsequently cloned, with the functional recombinant protein expressed in Escherichia coli and characterized. PPDBR readily converted both dehydrodiconiferyl and coniferyl aldehydes into dihydrodehydrodiconiferyl and dihydroconiferyl aldehydes, when NADPH was added as cofactor. However, it was unable to reduce directly either the double bond of dehydrodiconiferyl or coniferyl alcohols in the presence of NADPH. During this reductive step, the corresponding 4-proR hydrogen was abstracted from [4R-3H]-NADPH during hydride transfer. This is thus the first report of a double-bond reductase involved in phenylpropanoid metabolism, and which is presumed to be involved in plant defense. In situ mRNA hybridization indicated that the PPDBR transcripts in P. taeda stem sections were localized to the vascular cambium, as well as to radial and axial parenchyma cell types. Additionally, using P. taeda cell suspension culture crude protein extracts, dehydrodiconiferyl and coniferyl alcohols could be dehydrogenated to afford dehydrodiconiferyl and coniferyl aldehydes. Furthermore, these same extracts were able to convert dihydrodehydrodiconiferyl and dihydroconiferyl aldehydes into the corresponding alcohols. Taken together, these results indicate that in the crude extracts dehydrodiconiferyl and coniferyl alcohols can be converted to dihydrodehydrodiconiferyl and dihydroconiferyl alcohols through a three-step process, i.e. by initial phenylpropenol oxidation, then sequential PPDBR and phenylpropanal reductions, respectively.