Detection and genetic evolution analysis of porcine parvovirus type 7 (PPV7) in Fujian Province.

Porcine parvovirus (PPV) is an important pathogen causing reproductive disorders in sows, with clinical symptoms including stillbirth, mummified fetuses, embryonic dysplasia and death, and sow infertility. Porcine parvovirus 7 (PPV7) is a recently discovered type of PPV and its widespread distribution and rapid evolution has caused huge economic losses in the pig industry. To investigate the molecular epidemiology of PPV7 in Fujian Province, China, we collected 491 blood samples and 72 tissue samples from diseased pigs in large-scale pig farms across selected areas of Fujian Province from 2019 to 2022. PPV7 infection was determined using real-time quantitative PCR, and positive samples underwent whole-genome amplification, sequencing, and subsequent homology, phylogenetic, and recombination analyses. The PPV7 positive detection rate was 25.73% (145/563) in Fujian Province, among which the positive rate of blood and tissue samples was 26.47% (130/491) and 20.83% (15/72), respectively. The nucleotide sequence homology among the 29 PPV7 whole-genome sequences obtained in this study was 90.0%-97.2%, whereas that with 128 reference strains from China and other countries was 88.9%-98.1%. Six strains had partial nucleotide deletions or insertions. Phylogenetic analysis based on the whole-genome sequences classified the 29 PPV7 strains and 128 reference strains into eight subtypes (PPV7a-PPV7h), and PPV7h was the predominant subtype in Fujian Province. Recombination analysis revealed evidence of inferred recombination events in the genomes of four strains. This study provides significant insights into the molecular characteristics of PPV7 in Fujian Province and serves as a crucial foundation for further advancements in PPV7 prevention and control strategies.

Variations in the NSP4 gene of the type 2 porcine reproductive and respiratory syndrome virus isolated in China from 1996 to 2021.

Porcine reproductive and respiratory syndrome virus (PRRSV) has continuously mutated since its first isolation in China in 1996, leading to difficulties in infection prevention and control. Infections caused by PRRSV-2 strains are the main epidemic strains in China, as determined by phylogenetic analysis. In this study, we focused on the prevalence and genetic variations of the non-structural protein 4 (NSP4) from PRRSV-2 over the past 20 years in China. The fundamental biological properties of the NSP4 were predicted, and an analysis and comparison of NSP4 homology at the nucleotide and amino acid levels was conducted using 123 PRRSV-2 strains. The predicted molecular weight of the NSP4 protein was determined to be 21.1 kDa, and it was predicted to be a stable hydrophobic protein that lacks a signal peptide. NSP4 from different strains exhibited a high degree of amino acid (85.8-100%) and nucleotide sequence homology (81.0-100%). Multiple amino acid substitutions were identified in NSP4 among 15 representative PRRSV-2 strains. Phylogenetic analysis showed that the lineage 8 and 1 strains, the most prevalent strains in China, were indifferent clades with a long genetic distance. This analysis will help fully elucidate the parameters of the PRRSV NSP4 epidemic in China to lay a foundation for adequate understanding of the function of NSP4. Genetic information results from the accumulation of conserved and non-conserved sequences. The high conservation of the NSP4 gene determines the most basic life traits and functions of PRRSV. Analyzing the spatial structure of NSP4 protein and studying the genetic evolution of NSP4 not only provide the theoretical basis for how NSP4 participates in the regulation of the innate response of the host but also provide a target for genetic manipulation and a reasonable target molecule and structure for new drug molecules.

Virus-specific editing identification approach reveals the landscape of A-to-I editing and its impacts on SARS-CoV-2 characteristics and evolution.

Upon SARS-CoV-2 infection, viral intermediates specifically activate the IFN response through MDA5-mediated sensing and accordingly induce ADAR1 p150 expression, which might lead to viral A-to-I RNA editing. Here, we developed an RNA virus-specific editing identification pipeline, surveyed 7622 RNA-seq data from diverse types of samples infected with SARS-CoV-2, and constructed an atlas of A-to-I RNA editing sites in SARS-CoV-2. We found that A-to-I editing was dynamically regulated, varied between tissue and cell types, and was correlated with the intensity of innate immune response. On average, 91 editing events were deposited at viral dsRNA intermediates per sample. Moreover, editing hotspots were observed, including recoding sites in the spike gene that affect viral infectivity and antigenicity. Finally, we provided evidence that RNA editing accelerated SARS-CoV-2 evolution in humans during the epidemic. Our study highlights the ability of SARS-CoV-2 to hijack components of the host antiviral machinery to edit its genome and fuel its evolution, and also provides a framework and resource for studying viral RNA editing.

A Case Study of the Morphological and Molecular Variation within a Ciliate Genus: Taxonomic Descriptions of Three Dysteria Species (Ciliophora, Cyrtophoria), with the Establishment of a New Species.

Three Dysteria species, D. crassipes Claparède & Lachmann, 1859; D. brasiliensis Faria et al., 1922; and D. paracrassipes n. sp., were collected from subtropical coastal waters of the East China Sea, near Ningbo, China. The three species were studied based on their living morphology, infraciliature, and molecular data. The new species D. paracrassipes n. sp. is very similar to D. crassipes in most morphological features except the preoral kinety, which is double-rowed in the new species (vs. single-rowed in D. crassipes). The difference in the small ribosomal subunit sequences (SSU rDNA) between these two species is 56 bases, supporting the establishment of the new species. The Ningbo population of D. crassipes is highly similar in morphology to other known populations. Nevertheless, the SSU rDNA sequences of these populations are very different, indicating high genetic diversity and potentially cryptic species. Dysteria brasiliensis is cosmopolitan with many described populations worldwide and four deposited SSU rDNA sequences. The present work supplies morphological and molecular information from five subtropical populations of D. brasiliensis that bear identical molecular sequences but show significant morphological differences. The findings of this study provide an opportunity to improve understanding of the morphological and genetic diversity of ciliates.

Identification and functional characterization of interleukin-12 receptor beta 1 and 2 in grass carp (Ctenopharyngodon idella).

Interleukin 12 (IL-12) binds its receptor complex of IL-12 receptor beta 1 (IL-12Rß1) and IL-12Rß2 to transduce cellular signaling in mammals. In teleosts, the function of Il-12 is drawing increasing attention, but molecular and functional features of Il-12 receptors remain obscure. Especially, the existence of multiple Il-12 isoforms in some fish species elicits the requirement to clarify their receptors. In this study, we isolated three cDNA sequences as Il-12 receptor candidates from grass carp, entitled as grass carp Il-12rß1 (gcIl-12rß1), gcIl-12rß2a and gcIl-12rß2b. In silico analysis showed that gcIl-12rß1 and gcIl-12rß2a shared the conserved gene locus and similar structure characteristics with their orthologues of zebrafish, frog, chicken, mouse and human, respectively. However, the Il-12rß2b of grass carp and zebrafish was similar to IL-27Ra in non-fish species. Further locally installed BLAST and gene synteny analysis uncovered three gcIl-12 receptors being single copied genes. Tissue distribution assay revealed that gcil12rß1 and gcil12rß2a transcripts were predominantly expressed in head kidney, differing from the even distribution of gcil12rß2b transcripts in all detected tissues. Subsequently, the binding ability and antagonistic effects of recombinant extracellular region of gcIl-12rß1 with recombinant grass carp Il-12 (rgcIl-12) isoforms were explored, providing functional evidence of the newly cloned gcIl-12rß1 being genuine orthologues of mammalian IL-12Rß1. Moreover, our data showed that gcIl-12rß1 and gcIl-12rß2a but not gcIl-12rß1 and gcIl-12rß2b mediated the effects of rgcIl-12 isoforms on ifn-γ promoter activity, thereby revealing Il-12 receptor signaling in fish. These results identified grass carp Il-12 receptors, thereby advancing our understanding of Il-12 isoform signaling in fish.

Long non­coding RNA CTBP1­AS2 upregulates USP22 to promote pancreatic carcinoma progression by sponging miR­141­3p.

Long non­coding RNAs (lncRNAs) feature prominently in pancreatic carcinoma progression. The present study aimed to clarify the biological functions, clinical significance and underlying mechanism of lncRNA CTBP1 antisense RNA 2 (CTBP1­AS2) in pancreatic carcinoma. Reverse transcription­quantitative PCR was performed to assess the expression levels of CTBP1­AS2, microRNA (miR)­141­3p and ubiquitin­specific protease 22 (USP22) mRNA in pancreatic carcinoma tissues and cell lines. Western blotting was used to examine USP22 protein expression in pancreatic carcinoma cell lines. Loss­of­function experiments were used to analyze the regulatory effects of CTBP1­AS2 on proliferation, apoptosis, migration and invasion of pancreatic carcinoma cells. Dual­luciferase reporter assay was used to examine the binding relationship between CTBP1­AS2 and miR­141­3p, as well as between miR­141­3p and USP22. It was demonstrated that CTBP1­AS2 expression was markedly increased in pancreatic carcinoma tissues and cell lines. High CTBP1­AS2 expression was associated with advanced clinical stage and lymph node metastasis of patients. Functional experiments confirmed that knocking down CTBP1­AS2 significantly inhibited pancreatic carcinoma cell proliferation, migration and invasion, and promoted cell apoptosis. In terms of mechanism, it was found that CTBP1­AS2 adsorbed miR­141­3p as a molecular sponge to upregulate the expression level of USP22. In conclusion, lncRNA CTBP1­AS2 may be involved in pancreatic carcinoma progression by regulating miR­141­3p and USP22 expressions; in addition, CTBP1­AS2 may be a diagnostic biomarker and treatment target for pancreatic carcinoma.

RNALocate v2.0: an updated resource for RNA subcellular localization with increased coverage and annotation.

Resolving the spatial distribution of the transcriptome at a subcellular level can increase our understanding of biology and diseases. To facilitate studies of biological functions and molecular mechanisms in the transcriptome, we updated RNALocate, a resource for RNA subcellular localization analysis that is freely accessible at http://www.rnalocate.org/ or http://www.rna-society.org/rnalocate/. Compared to RNALocate v1.0, the new features in version 2.0 include (i) expansion of the data sources and the coverage of species; (ii) incorporation and integration of RNA-seq datasets containing information about subcellular localization; (iii) addition and reorganization of RNA information (RNA subcellular localization conditions and descriptive figures for method, RNA homology information, RNA interaction and ncRNA disease information) and (iv) three additional prediction tools: DM3Loc, iLoc-lncRNA and iLoc-mRNA. Overall, RNALocate v2.0 provides a comprehensive RNA subcellular localization resource for researchers to deconvolute the highly complex architecture of the cell.

Physiological characteristics and miRNA sequencing of two root zones with contrasting ammonium assimilation patterns in Populus.

BACKGROUND: The net ammonium fluxes differ among the different root zones of Populus, but the physiological and microRNA regulatory mechanisms are unclear. OBJECTIVE: To elucidate the physiological and miRNA regulatory mechanisms, we investigated the two root zones displaying significant differences in net NH4+ effluxes of P. × canescens. METHODS: Populus plantlets were cultivated with 500 µM NH4Cl for 10 days. Six plants were randomly selected to determine the net NH4+ fluxes using a noninvasive microtest technique. High-throughput sequencing were used to determine the dynamic expression profile of miRNA among the different root zones of Populus. RESULTS: Net NH4+ efflux in zone I (from 0 to 40 mm from the root apex) was - 19.64 pmol cm-2 s-1 and in zone II (from 40 to 80 mm) it was - 43.96 pmol cm-2 s-1. The expression of eleven miRNAs was significantly upregulated, whereas fifteen miRNAs were downregulated. Moreover, eighty-eight target genes of the significantly differentially expressed miRNAs were identified in root zone II compared with zone I. Particularly, ptc-miR171a/b/e and their target, SCL6, were found to be important for the difference in net NH4+ effluxes in the two root zones. Moreover, the expression of the target of ptc-miR169d, NFYA3 was upregulated in root zone II compared with root zone I, contributing to increased NH4+ efflux and decreased NH4+ assimilation in root zone II. CONCLUSION: These results indicate that miRNAs regulate the expression levels of their target genes and thus play key roles in net NH4+ fluxes and NH4+ assimilation in different poplar root zones.

CircN4BP2L2 promotes colorectal cancer growth and metastasis through regulation of the miR-340-5p/CXCR4 axis.

Colorectal cancer (CRC) is the third leading cause of cancer-related death worldwide. Dysregulation of circular RNAs (circRNAs) appears to be a critical factor in CRC progression. However, mechanistic studies delineating the role of circRNAs in CRC remain limited. In this study, qRT-PCR and western blot assays were used to measure the expression of genes and proteins. Migration, invasion, proliferation, and apoptosis were examined by wound-healing, transwell, CCK-8, colony formation, and flow cytometry assays, respectively. Molecular interactions were validated by a dual-luciferase report system. A xenograft animal model was established to examine in vivo tumor growth and lung metastasis. Our data indicated that circN4BP2L2 expression was increased in CRC tissues and cell lines. Notably, inhibition of circN4BP2L2 effectively inhibited proliferation, migration, and invasion of LoVo cells, and inhibited tumor growth and metastasis in vivo, whereas the forced expression of circN4BP2L2 facilitated the proliferation, migration, and invasion of HT-29 cells. Mechanistic studies revealed that circN4BP2L2 acted as a molecular sponge of miR-340-5p to competitively promote CXCR4 expression. Furthermore, inhibition of miR-340-5p reversed the anti-cancer effects of circN4BP2L2 or CXCR4 silencing. Our data indicated an oncogenic role of circN4BP2L2 in CRC via regulation of the miR-340-5p/CXCR4 axis, which may be a promising biomarker and target for CRC treatment.

Sequence complementarity between human noncoding RNAs and SARS-CoV-2 genes: What are the implications for human health?

OBJECTIVES: To investigate in silico the presence of nucleotide sequence complementarity between the RNA genome of Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2) and human non-coding (nc)RNA genes. METHODS: The FASTA sequence (NC_045512.2) of each of the 11 SARS-CoV-2 isolate Wuhan-Hu-1 genes was retrieved from NCBI.nlm.nih.gov/gene and the Ensembl.org library interrogated for any base-pair match with human ncRNA genes. SARS-CoV-2 gene-matched human ncRNAs were screened for functional activity using bioinformatic analysis. Finally, associations between identified ncRNAs and human diseases were searched in GWAS databases. RESULTS: A total of 252 matches were found between the nucleotide sequence of SARS-CoV-2 genes and human ncRNAs. With the exception of two small nuclear RNAs, all of them were long non-coding (lnc)RNAs expressed mainly in testis and central nervous system under physiological conditions. The percentage of alignment ranged from 91.30% to 100% with a mean nucleotide alignment length of 17.5 ± 2.4. Thirty-three (13.09%) of them contained predicted R-loop forming sequences, but none of these intersected the complementary sequences of SARS-CoV-2. However, in 31 cases matches fell on ncRNA regulatory sites, whose adjacent coding genes are mostly involved in cancer, immunological and neurological pathways. Similarly, several polymorphic variants of detected non-coding genes have been associated with neuropsychiatric and proliferative disorders. CONCLUSION: This pivotal in silico study shows that SARS-CoV-2 genes have Watson-Crick nucleotide complementarity to human ncRNA sequences, potentially disrupting ncRNA epigenetic control of target genes. It remains to be elucidated whether this could result in the development of human disease in the long term.

Differential Regulation of Human Surfactant Protein A Genes, SFTPA1 and SFTPA2, and Their Corresponding Variants.

The human SFTPA1 and SFTPA2 genes encode the surfactant protein A1 (SP-A1) and SP-A2, respectively, and they have been identified with significant genetic and epigenetic variability including sequence, deletion/insertions, and splice variants. The surfactant proteins, SP-A1 and SP-A2, and their corresponding variants play important roles in several processes of innate immunity as well in surfactant-related functions as reviewed elsewhere [1]. The levels of SP-A have been shown to differ among individuals both under baseline conditions and in response to various agents or disease states. Moreover, a number of agents have been shown to differentially regulate SFTPA1 and SFTPA2 transcripts. The focus in this review is on the differential regulation of SFTPA1 and SFTPA2 with primary focus on the role of 5' and 3' untranslated regions (UTRs) and flanking sequences on this differential regulation as well molecules that may mediate the differential regulation.

Genome-wide promoter analysis, homology modeling and protein interaction network of Dehydration Responsive Element Binding (DREB) gene family in Solanum tuberosum.

Dehydration Responsive Element Binding (DREB) regulates the expression of numerous stress-responsive genes, and hence plays a pivotal role in abiotic stress responses and tolerance in plants. The study aimed to develop a complete overview of the cis-acting regulatory elements (CAREs) present in S. tuberosum DREB gene promoters. A total of one hundred and four (104) cis-regulatory elements (CREs) were identified from 2.5kbp upstream of the start codon (ATG). The in-silico promoter analysis revealed variable sets of cis-elements and functional diversity with the predominance of light-responsive (30%), development-related (20%), abiotic stress-responsive (14%), and hormone-responsive (12%) elements in StDREBs. Among them, two light-responsive elements (Box-4 and G-box) were predicted in 64 and 61 StDREB genes, respectively. Two development-related motifs (AAGAA-motif and as-1) were abundant in StDREB gene promoters. Most of the DREB genes contained one or more Myeloblastosis (MYB) and Myelocytometosis (MYC) elements associated with abiotic stress responses. Hormone-responsive element i.e. ABRE was found in 59 out of 66 StDREB genes, which implied their role in dehydration and salinity stress. Moreover, six proteins were chosen corresponding to A1-A6 StDREB subgroups for secondary structure analysis and three-dimensional protein modeling followed by model validation through PROCHECK server by Ramachandran Plot. The predicted models demonstrated >90% of the residues in the favorable region, which further ensured their reliability. The present study also anticipated pocket binding sites and disordered regions (DRs) to gain insights into the structural flexibility and functional annotation of StDREB proteins. The protein association network determined the interaction of six selected StDREB proteins with potato proteins encoded by other gene families such as MYB and NAC, suggesting their similar functional roles in biological and molecular pathways. Overall, our results provide fundamental information for future functional analysis to understand the precise molecular mechanisms of the DREB gene family in S. tuberosum.

CircRNA75 and CircRNA72 Function as the Sponge of MicroRNA-200 to Suppress Coelomocyte Apoptosis Via Targeting Tollip in Apostichopus japonicus.

Circular RNAs (circRNAs) act as essential regulators in many biological processes, especially in mammalian immune response. Nonetheless, the functions and mechanisms of circRNAs in the invertebrate immune system are largely unclarified. In our previous work, 261 differentially expressed circRNAs potentially related to the development of Apostichopus japonicus skin ulceration syndrome (SUS), which is a major problem restricting the sea cucumber breeding industry, were identified by genome-wide screening. In this study, via miRanda analysis, both circRNA75 and circrRNA72 were shown to share the miR-200 binding site, a key microRNA in the SUS. The two circRNAs were verified to be increased significantly in LPS-exposed primary coelomocytes, similar to the results of circRNA-seq in sea cucumber under Vibrio splendidus-challenged conditions. A dual-luciferase assay indicated that both circRNA75 and circRNA72 could bind miR-200 in vivo, in which circRNA75 had four binding sites of miR-200 and only one for circRNA72. Furthermore, we found that miR-200 could bind the 3'-UTR of Toll interacting protein (Tollip) to negatively mediate the expression of Tollip. Silencing Tollip increased primary coelomocyte apoptosis. Consistently, inference of circRNA75 and circRNA72 could also downregulate Tollip expression, thereby increasing the apoptosis of primary coelomocytes, which could be blocked by miR-200 inhibitor treatment. Moreover, the rate of si-circRNA75-downregulated Tollip expression was higher than that of si-circRNA72 under an equivalent amount. CircRNA75 and circRNA72 suppressed coelomocyte apoptosis by sponging miR-200 to promote Tollip expression. The ability of circRNA to adsorb miRNA might be positively related to the number of binding sites for miRNA.

Negative feedback for DARS2-Fis complex by ATP-DnaA supports the cell cycle-coordinated regulation for chromosome replication.

In Escherichia coli, the replication initiator DnaA oscillates between an ATP- and an ADP-bound state in a cell cycle-dependent manner, supporting regulation for chromosome replication. ATP-DnaA cooperatively assembles on the replication origin using clusters of low-affinity DnaA-binding sites. After initiation, DnaA-bound ATP is hydrolyzed, producing initiation-inactive ADP-DnaA. For the next round of initiation, ADP-DnaA binds to the chromosomal locus DARS2, which promotes the release of ADP, yielding the apo-DnaA to regain the initiation activity through ATP binding. This DnaA reactivation by DARS2 depends on site-specific binding of IHF (integration host factor) and Fis proteins and IHF binding to DARS2 occurs specifically during pre-initiation. Here, we reveal that Fis binds to an essential region in DARS2 specifically during pre-initiation. Further analyses demonstrate that ATP-DnaA, but not ADP-DnaA, oligomerizes on a cluster of low-affinity DnaA-binding sites overlapping the Fis-binding region, which competitively inhibits Fis binding and hence the DARS2 activity. DiaA (DnaA initiator-associating protein) stimulating ATP-DnaA assembly enhances the dissociation of Fis. These observations lead to a negative feedback model where the activity of DARS2 is repressed around the time of initiation by the elevated ATP-DnaA level and is stimulated following initiation when the ATP-DnaA level is reduced.

Controlled RISC loading efficiency of miR168 defined by miRNA duplex structure adjusts ARGONAUTE1 homeostasis.

Micro RNAs (miRNAs) are processed from precursor RNA molecules with precisely defined secondary stem-loop structures. ARGONAUTE1 (AGO1) is the main executor component of miRNA pathway and its expression is controlled via the auto-regulatory feedback loop activity of miR168 in plants. Previously we have shown that AGO1 loading of miR168 is strongly restricted leading to abundant cytoplasmic accumulation of AGO-unbound miR168. Here, we report, that intrinsic RNA secondary structure of MIR168a precursor not only defines the processing of miR168, but also precisely adjusts AGO1 loading efficiency determining the biologically active subset of miR168 pool. Our results show, that modification of miRNA duplex structure of MIR168a precursor fragment or expression from artificial precursors can alter the finely adjusted loading efficiency of miR168. In dcl1-9 mutant where, except for miR168, production of most miRNAs is severely reduced this mechanism ensures the elimination of unloaded AGO1 proteins via enhanced AGO1 loading of miR168. Based on this data, we propose a new competitive loading mechanism model for miR168 action: the miR168 surplus functions as a molecular buffer for controlled AGO1 loading continuously adjusting the amount of AGO1 protein in accordance with the changing size of the cellular miRNA pool.

The hyperthermophilic archaeon Thermococcus kodakarensis is resistant to pervasive negative supercoiling activity of DNA gyrase.

In all cells, DNA topoisomerases dynamically regulate DNA supercoiling allowing essential DNA processes such as transcription and replication to occur. How this complex system emerged in the course of evolution is poorly understood. Intriguingly, a single horizontal gene transfer event led to the successful establishment of bacterial gyrase in Archaea, but its emergent function remains a mystery. To better understand the challenges associated with the establishment of pervasive negative supercoiling activity, we expressed the gyrase of the bacterium Thermotoga maritima in a naïve archaeon Thermococcus kodakarensis which naturally has positively supercoiled DNA. We found that the gyrase was catalytically active in T. kodakarensis leading to strong negative supercoiling of plasmid DNA which was stably maintained over at least eighty generations. An increased sensitivity of gyrase-expressing T. kodakarensis to ciprofloxacin suggested that gyrase also modulated chromosomal topology. Accordingly, global transcriptome analyses revealed large scale gene expression deregulation and identified a subset of genes responding to the negative supercoiling activity of gyrase. Surprisingly, the artificially introduced dominant negative supercoiling activity did not have a measurable effect on T. kodakarensis growth rate. Our data suggest that gyrase can become established in Thermococcales archaea without critically interfering with DNA transaction processes.

Transcriptome analysis reveals the molecular mechanisms of response to an emergent yellow-flower disease in green Chinese prickly ash (Zanthoxylum schinifolium).

Chinese prickly ash (Zanthoxylum) is extensively used as spice and traditional medicine in eastern Asian countries. Recently, an emergent yellow-flower disease (YFD) break out in green Chinese prickly ash (Zanthoxylum schinifolium, Qinghuajiao in Chinese) at Chongqing municipality, and then leads to a sharp reduction in the yield of Qinghuajiao, and thus results in great economic losses for farmers. To address the molecular response for the emergent YFD of Qinghuajiao, we analyzed the transcriptome of 12 samples including the leaves and inflorescences of asymptomatic and symptomatic plants from three different towns at Chongqing by high-throughput RNA-Seq technique. A total of 126,550 genes and 229,643 transcripts were obtained, and 21,054 unigenes were expressed in all 12 samples. There were 56 and 164 different expressed genes (DEGs) for the AL_vs_SL (asymptomatic leaf vs symptomatic leaf) and AF_vs_SF (asymptomatic flower vs symptomatic flower) groups, respectively. The results of KEGG analysis showed that the "phenylpropanoid biosynthesis" pathway that related to plant-pathogen interaction were found in AL_vs_SL and AF_vs_SF groups, and the "Plant-pathogen interaction" found in AF_vs_SF group, implying that this Qinghuajiao YFD might cause by plant pathogen. Interestingly, we detected 33 common unigenes for the 2 groups, and almost these unigenes were up-regulated in the symptomatic plants. Moreover, most of which were homologs to virus RNA, the components of viruses, implying that this YFD was related to virus. Our results provided a primary molecular basis for the prevention and treatment of YFD of Qinghuajiao trees.

Spatial regulation by multiple Gremlin1 enhancers provides digit development with cis-regulatory robustness and evolutionary plasticity.

Precise cis-regulatory control of gene expression is essential for normal embryogenesis and tissue development. The BMP antagonist Gremlin1 (Grem1) is a key node in the signalling system that coordinately controls limb bud development. Here, we use mouse reverse genetics to identify the enhancers in the Grem1 genomic landscape and the underlying cis-regulatory logics that orchestrate the spatio-temporal Grem1 expression dynamics during limb bud development. We establish that transcript levels are controlled in an additive manner while spatial regulation requires synergistic interactions among multiple enhancers. Disrupting these interactions shows that altered spatial regulation rather than reduced Grem1 transcript levels prefigures digit fusions and loss. Two of the enhancers are evolutionary ancient and highly conserved from basal fishes to mammals. Analysing these enhancers from different species reveal the substantial spatial plasticity in Grem1 regulation in tetrapods and basal fishes, which provides insights into the fin-to-limb transition and evolutionary diversification of pentadactyl limbs.

Calculating Orthologous Protein-Coding Sequence Set Probability Using the Poisson Process.

We extend the popular Jukes-Cantor evolution model and calculate the probability of an orthologous nucleotide sequence set [a reference sequence (B1) stays with the other sequences (B-1)], where the sequence evolution [from a last common ancestral sequence (ɑ)] follows the (prospective) Poisson process with the overall event rate λ prorated among mutation types (nucleotide/codon substitution, insertion, and deletion) and sites along each sequence. The corresponding retrospective process (reversing the prospective process) facilitates developing algorithms to calculate the marginal probability [Pr(B1)] (Monte Carlo integration) and sample ɑ (given B1). We calculate probability Pr(B-1|ɑ) based on the identified events (during "ɑ→B-1") from pairwise sequence alignment to implement Pr(B-1|B1) calculation (Monte Carlo integration). Event queue sampling and probability magnifiers are used to improve the computational efficiency when the number of events is large. We finally test our procedure on both simulated and recently studied hexapod transcriptome data (Brandt et al.), where each asexual lineage pairs with its closest related sexual lineage. Rate estimates (for Phasmatodea and Zygentoma) and model comparison indicate that the asexual lineages likely mutate several times faster than their sexual relatives.

AGO CLIP-based imputation of potent siRNA sequences targeting SARS-CoV-2 with antifibrotic miRNA-like activity.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is associated with fatal pulmonary fibrosis. Small interfering RNAs (siRNAs) can be developed to induce RNA interference against SARS-CoV-2, and their susceptible target sites can be inferred by Argonaute crosslinking immunoprecipitation sequencing (AGO CLIP). Here, by reanalysing AGO CLIP data in RNA viruses, we delineated putative AGO binding in the conserved non-structural protein 12 (nsp12) region encoding RNA-dependent RNA polymerase (RdRP) in SARS-CoV-2. We utilised the inferred AGO binding to optimise the local RNA folding parameter to calculate target accessibility and predict all potent siRNA target sites in the SARS-CoV-2 genome, avoiding sequence variants. siRNAs loaded onto AGO also repressed seed (positions 2-8)-matched transcripts by acting as microRNAs (miRNAs). To utilise this, we further screened 13 potential siRNAs whose seed sequences were matched to known antifibrotic miRNAs and confirmed their miRNA-like activity. A miR-27-mimicking siRNA designed to target the nsp12 region (27/RdRP) was validated to silence a synthesised nsp12 RNA mimic in lung cell lines and function as an antifibrotic miR-27 in regulating target transcriptomes related to TGF-ß signalling. siRNA sequences with an antifibrotic miRNA-like activity that could synergistically treat COVID-19 are available online ( http://clip.korea.ac.kr/covid19 ).