Circular RNAs (circRNAs) are a class of non-coding RNAs that are formed via back-splicing. These circRNAs are predominantly studied for their roles as regulators of various biological processes. Notably, emerging evidence demonstrates that host circRNAs can be differentially expressed (DE) upon infection with pathogens (e.g., influenza and coronaviruses), suggesting a role for circRNAs in regulating host innate immune responses. However, investigations on the role of circRNAs during pathogenic infections are limited by the knowledge and skills required to carry out the necessary bioinformatic analysis to identify DE circRNAs from RNA sequencing (RNA-seq) data. Bioinformatics prediction and identification of circRNAs is crucial before any verification, and functional studies using costly and time-consuming wet-lab techniques. To solve this issue, a step-by-step protocol of in silico prediction and characterization of circRNAs using RNA-seq data is provided in this manuscript. The protocol can be divided into four steps: 1) Prediction and quantification of DE circRNAs via the CIRIquant pipeline; 2) Annotation via circBase and characterization of DE circRNAs; 3) CircRNA-miRNA interaction prediction through Circr pipeline; 4) functional enrichment analysis of circRNA parental genes using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). This pipeline will be useful in driving future in vitro and in vivo research to further unravel the role of circRNAs in host-pathogen interactions.
To examine the impact of regional and seasonal variations on the incidence and severity of sheath rot, a major seed-borne disease of rice caused by Sarocladium oryzae, data on incidence and severity were collected from 27 selected fields in the Gazipur, Rangpur, Bogra, Chittagong, Comilla, Gopalgonj, Jessore, Manikgonj, and Bhola districts of Bangladesh in rain-fed and irrigated conditions. Cultural variability of 29 pathogen isolates obtained from 8 different locations was studied on potato dextrose agar (PDA) and genetic variability was determined by DNA fingerprinting using variable number tandem repeat-polymerase chain reaction markers. Overall, disease incidence and severity were higher in irrigated rice. Disease incidence and severity were highest in the Bhola district in rain-fed rice and lowest in irrigated rice. Mycelial growth of 29 representative isolates was found to vary on PDA and the isolates were divided into 6 groups. The range of the overall size of conidia of the selected isolates was 2.40-7.20 x 1.20-2.40 μm. Analysis of the DNA fingerprint types of the 29 isolates of S. oryzae, obtained from the amplification reactions, revealed 10 fingerprinting types (FPTs) that were 80% similar. FPT-1 was the largest group and included 13 isolates (44.8%), while FPT-2 was the third largest group and included 3 isolates. Each of FPT-3, 4, 5, and 6 included only 1 isolate. We observed no relationship between cultural and genetic groupings.
The rice blast fungus Magnaporthe oryzae is a serious pathogen that jeopardises the world's most important food-security crop. Ten common Malaysian rice varieties were examined for their morphological, physiological and genomic responses to this rice blast pathogen. qPCR quantification was used to assess the growth of the pathogen population in resistant and susceptible rice varieties. The chlorophyll content and photosynthesis were also measured to further understand the disruptive effects that M. oryzae has on infected plants of these varieties. Real-time PCR was used to explore the differential expression of eight blast resistance genes among the ten local varieties. Blast disease has destructive effects on the growth of rice, and the findings of our study provide evidence that the Pikh, Pi9, Pi21, and Osw45 genes are involved in defence responses in the leaves of Malaysian rice at 31 h after inoculation with M. oryzae pathotype P7.2. Both the chlorophyll content and photosynthesis were reduced, but the levels of Pikh gene expression remained constant in susceptible varieties, with a developed pathogen population and mild or severe symptoms. The Pi9, Pi21, and Osw45 genes, however, were simultaneously upregulated in infected rice plants. Therefore, the presence of the Pikh, Pi9, Pi21, and Osw45 genes in the germplasm is useful for improving the resistance of rice varieties.
Transposable elements (TEs) are agents of genetic variability in phytopathogens as they are a source of adaptive evolution through genome diversification. Although many studies have uncovered information on TEs, the exact mechanism behind TE-induced changes within the genome remains poorly understood. Furthermore, convergent trends towards bigger genomes, emergence of novel genes and gain or loss of genes implicate a TE-regulated genome plasticity of fungal phytopathogens. TEs are able to alter gene expression by revamping the cis-regulatory elements or recruiting epigenetic control. Recent findings show that TEs recruit epigenetic control on the expression of effector genes as part of the coordinated infection strategy. In addition to genome plasticity and diversity, fungal pathogenicity is an area of economic concern. A survey of TE distribution suggests that their proximity to pathogenicity genes TEs may act as sites for emergence of novel pathogenicity factors via nucleotide changes and expansion or reduction of the gene family. Through a systematic survey of literature, we were able to conclude that the role of TEs in fungi is wide: ranging from genome plasticity, pathogenicity to adaptive behavior in evolution. This review also identifies the gaps in knowledge that requires further elucidation for a better understanding of TEs' contribution to genome architecture and versatility.
The detection and identification of two endocytobiotic bacterial strains, one affiliated to the "Candidatus Caedibacter acanthamoebae"/"Ca. Paracaedimonas acanthamoeba", and another to the endosymbiont of Acanthamoeba UWC8 and "Ca. Jidaibacter acanthamoeba" are described. For endocytobiont screening, we developed a PCR method with a set of broad-range bacterial 16S rRNA primers to substitute the commonly used but technically demanding fluorescent in situ hybridization technique. Our PCR test alone without sequencing failed to discriminate the endocytobiont-containing and endocytobiont-free Acanthamoeba sp. due to the presence of mismatched primers to host mitochondrial DNA. We highlighted the need to perform bacterial primer checking against the Acanthamoeba genome to avoid false positive detection in PCR. Although the genetic aspect of "Ca. Caedibacter acanthamoebae"/"Ca. Paracaedimonas acanthamoeba" and the endosymbiont of Acanthamoeba UWC8/"Ca. Jidaibacter acanthamoeba" are well studied, knowledge pertaining to their morphologies are quite vague. Hence, we used transmission electron microscopy to examine our endocytobionts which are affiliated to previously described intracellular bacteria of Acanthamoeba sp. We used good-quality TEM images for the localization and the fate of the current endocytobionts inside different life stages of the hosts. Furthermore, to the best of our knowledge, our TEM findings are the first to provide morphological evidence for the clearance of defective Acanthamoeba endocytobionts via an autophagic-like process.
Black pepper is an important commodity crop in Malaysia that generates millions of annual revenue for the country. However, black pepper yield is affected by slow decline disease caused by a soil-borne fungus Fusarium solani. RNA sequencing transcriptomics approach has been employed in this study to explore the differential gene expression in susceptible Piper nigrum L. and resistant Piper colubrinum Link. Gene expression comparative analysis of the two pepper species has yielded 2,361 differentially expressed genes (DEGs). Among them, higher expression of 1,426 DEGs was detected in resistant plant. These DEGs practically demonstrated the major branches of plant-pathogen interaction pathway (Path: ko04626). We selected five groups of defence-related DEGs for downstream qRT-PCR analysis. Cf-9, the gene responsible for recognizing fungal avirulence protein activity was found inexpressible in susceptible plant. However, this gene exhibited promising expression in resistant plant. Inactivation of Cf-9 could be the factor that causes susceptible plant fail in recognition of F. solani and subsequently delay activation of adaptive response to fungal invasion. This vital study advance the understanding of pepper plant defence in response to F. solani and aid in identifying potential solution to manage slow decline disease in black pepper cultivation.
Zika virus (ZIKV) causes no-to-mild symptoms or severe neurological disorders. To investigate the importance of viral and host genetic variations in determining ZIKV infection outcomes, we created three full-length infectious cDNA clones as bacterial artificial chromosomes for each of three spatiotemporally distinct and genetically divergent ZIKVs: MR-766 (Uganda, 1947), P6-740 (Malaysia, 1966), and PRVABC-59 (Puerto Rico, 2015). Using the three molecularly cloned ZIKVs, together with 13 ZIKV region-specific polyclonal antibodies covering nearly the entire viral protein-coding region, we made three conceptual advances: (i) We created a comprehensive genome-wide portrait of ZIKV gene products and their related species, with several previously undescribed gene products identified in the case of all three molecularly cloned ZIKVs. (ii) We found that ZIKV has a broad cell tropism in vitro, being capable of establishing productive infection in 16 of 17 animal cell lines from 12 different species, although its growth kinetics varied depending on both the specific virus strain and host cell line. More importantly, we identified one ZIKV-non-susceptible bovine cell line that has a block in viral entry but fully supports the subsequent post-entry steps. (iii) We showed that in mice, the three molecularly cloned ZIKVs differ in their neuropathogenicity, depending on the particular combination of viral and host genetic backgrounds, as well as in the presence or absence of type I/II interferon signaling. Overall, our findings demonstrate the impact of viral and host genetic variations on the replication kinetics and neuropathogenicity of ZIKV and provide multiple avenues for developing and testing medical countermeasures against ZIKV.
Two dominant species of wild small rodents trapped in Novosibirsk region, South-Western Siberia, Russia differed in their susceptibility to the tick-borne encephalitis virus (TBEV) infection. TBEV RNA average detection rate for Northern red-backed vole Myodes rutilus (Pallas, 1779) (82.2 ± 5.8% blood samples and 63.1 ± 2.7% organ samples) significantly exceeded the corresponding values for the striped field mouse Apodemus agrarius (Pallas, 1771) (47.0 ± 8.7% blood and 24.5 ± 2.8% organ samples) (p <0.001). Innate immunity may be one of possible reasons of the differences. Th1 cytokine gene expression distinguished between M. rutilus (12.5 ± 8.5%) and A. agrarius (66.6 ± 11.4%), whereas Th2 cytokine frequencies were statistically similar (81.8 ± 12.2% and 100.0%, respectively). Polarization indexes (PI) of the innate immunity calculated as ratio of Th2 to Th1 cytokine RNA detection rates for both M. rutilus (6.5) and A. agrarius (1.5) suggested Th2 mainly humoral immune response against persistent TBEV in natural mammalian hosts. Therefore, the TBEV-induced antibodies were analyzed by ELISA and hemagglutination inhibition (HI) tests. The TBEV-specific antibodies were detected in 74.8 ± 4.3% sera of M. rutilus and 67.3 ± 6.8% of A. agrarius. Among them HI antibodies were found in 4.8 ± 2.1% of the same analyzed sera of M. rutilus and in 6.0 ± 3.4% blood samples of A. agrarius only. To model the TBEV persistence both M. rutilus and A. agrarius were infected with the suspensions of the TBEV-infected ticks with further observations during 4 subsequent months. Detection rate of the TBEV RNA and antigen E remained high during the whole period, however, pathogenic for laboratory suckling mice virus was isolated up to 8 days postinfection. At late stages of the persistent infection (1-4 months) the TBEV RNA detection rate in northern red-backed voles remained high 70.6 ± 7.9% whereas in striped field mice significantly declined to 26.7 ± 9.2% (p .05) but Th1 cytokine mRNA detection rates were different (44.4 ± 12.5% and 85.7 ± 9.7%, respectively) (p
The macaque parasite Plasmodium knowlesi is a significant concern in Malaysia where cases of human infection are increasing. Parasites infecting humans originate from genetically distinct subpopulations associated with the long-tailed (Macaca fascicularis (Mf)) or pig-tailed macaques (Macaca nemestrina (Mn)). We used a new high-quality reference genome to re-evaluate previously described subpopulations among human and macaque isolates from Malaysian-Borneo and Peninsular-Malaysia. Nuclear genomes were dimorphic, as expected, but new evidence of chromosomal-segment exchanges between subpopulations was found. A large segment on chromosome 8 originating from the Mn subpopulation and containing genes encoding proteins expressed in mosquito-borne parasite stages, was found in Mf genotypes. By contrast, non-recombining organelle genomes partitioned into 3 deeply branched lineages, unlinked with nuclear genomic dimorphism. Subpopulations which diverged in isolation have re-connected, possibly due to deforestation and disruption of wild macaque habitats. The resulting genomic mosaics reveal traits selected by host-vector-parasite interactions in a setting of ecological transition.
Influenza A viruses (IAVs) have a remarkable tropism in their ability to circulate in both mammalian and avian species. The IAV NS1 protein is a multifunctional virulence factor that inhibits the type I interferon host response through a myriad of mechanisms. How NS1 has evolved to enable this remarkable property across species and its specific impact in the overall replication, pathogenicity, and host preference remain unknown. Here we analyze the NS1 evolutionary landscape and host tropism using a barcoded library of recombinant IAVs. Results show a surprisingly great variety of NS1 phenotypes according to their ability to replicate in different hosts. The IAV NS1 genes appear to have taken diverse and random evolutionary pathways within their multiple phylogenetic lineages. In summary, the high evolutionary plasticity of this viral protein underscores the ability of IAVs to adapt to multiple hosts and aids in our understanding of its global prevalence.
Lipid biosynthesis produces glycerol, which is important in fueling turgor pressure necessary for germination and penetration of plant host by fungi. As the relationship between pathogenicity and the lipid biosynthetic pathway is not fully understood, we have elucidated the role of the fatty acid synthase beta subunit dehydratase (FAS1) gene in lipid biosynthesis. The FAS1 gene was silenced through homologous double crossover in Magnaporthe oryzae strain S6 to study the effect on lipid biosynthesis. The vegetative growth of Δfas1 mutants show the highest drop on oleic acid (between 10 and 50%), while the mycelial dry weight of mutants dropped significantly on all media. Conidiation of FAS1 mutants show a ~10- and ~5-fold reduction on oatmeal and Potato Dextrose Agar (PDA), respectively. Mutants formed mycelium that were mildly pigmented, indicating that the deletion of FAS1 may have affected melanin biosynthesis. Biochemical and gene expression studies concluded that the fatty acid degradation pathway might have been interrupted by FAS1 deletion. FAS1 mutants showed no enzyme activity on glucose or olive oil, suggesting that the mutants may lack functional peroxisomes and be defective in β-oxidation of fatty acids, hence explaining the reduced lipid deposits in the spores.
Chitinases play a vital role during the pathogenic invasion and immunosuppression in various organisms including invertebrates and vertebrates. In this study, we have investigated the participation of MrChit-3 (Macrobrachium rosenbergii Chitinase-3) during host-pathogenic interaction in freshwater prawn, M. rosenbergii. Quantitative real-time PCR analysis showed that the expression of MrChit-3 was up-regulated during bacterial, viral and laminarin challenge. Moreover, to understand the antimicrobial role of the GH18 domain, a putative membrane-targeting antimicrobial peptide (MrVG) was identified from the GH18 domain region of the protein and it was chemically synthesized. Physico-chemical features of the GH18 derived antimicrobial peptide (AMP) was assessed by various in silico tools and the antimicrobial property of the peptide was confirmed from in vitro studies. The membrane targeting mechanism of the peptide was determined by flow cytometry (FACS) and scanning electron microscope (SEM) analysis. Interestingly, the peptide was able to inhibit the growth of a chitinolytic fungal pathogen, Aspergillus niger, which was isolated from the shells of M. rosenbergii. The toxicity studies such as hemolysis activity on human blood erythrocytes and cell viability assay with primary kidney cells, HEK293 of MrVG revealed that the peptide was not involved in inducing any toxicity.
At present, very little is known about how Burkholderia pseudomallei (B. pseudomallei) interacts with its host to elicit melioidosis symptoms. We established a murine acute-phase melioidosis model and used DNA microarray technology to investigate the global host/pathogen interaction. We compared the transcriptome of infected liver and spleen with uninfected tissues over an infection period of 42 hr to identify genes whose expression is altered in response to an acute infection.
BACKGROUND: Human metapneumovirus (HMPV) is established as one of the causative agents of respiratory tract infections. To date, there are limited reports that describe the effect of HMPV genotypes and/or viral load on disease pathogenesis in adults. This study aims to determine the role of HMPV genetic diversity and nasopharyngeal viral load on symptom severity in outpatient adults with acute respiratory tract infections.
METHODS: Severity of common cold symptoms of patients from a teaching hospital was assessed by a four-category scale and summed to obtain the total symptom severity score (TSSS). Association between the fusion and glycoprotein genes diversity, viral load (quantified using an improved RT-qPCR assay), and symptom severity were analyzed using bivariate and linear regression analyses.
RESULTS: Among 81/3706 HMPV-positive patients, there were no significant differences in terms of demographics, number of days elapsed between symptom onset and clinic visit, respiratory symptoms manifestation and severity between different HMPV genotypes/sub-lineages. Surprisingly, elderly patients (≥65 years old) had lower severity of symptoms (indicated by TSSS) than young and middle age adults (p = 0.008). Nasopharyngeal viral load did not correlate with nor predict symptom severity of HMPV infection. Interestingly, at 3-5 days after symptom onset, genotype A-infected patients had higher viral load compared to genotype B (4.4 vs. 3.3 log10 RNA copies/μl) (p = 0.003).
CONCLUSIONS: Overall, HMPV genetic diversity and viral load did not impact symptom severity in adults with acute respiratory tract infections. Differences in viral load dynamics over time between genotypes may have important implications on viral transmission.
Study site: Primary Care Clinic, University of Malaya Medical Centre (UMMC), Kuala Lumpur, Malaysia
Interplay between EBV infection and acquired genetic alterations during nasopharyngeal carcinoma (NPC) development remains vague. Here we report a comprehensive genomic analysis of 70 NPCs, combining whole-genome sequencing (WGS) of microdissected tumor cells with EBV oncogene expression to reveal multiple aspects of cellular-viral co-operation in tumorigenesis. Genomic aberrations along with EBV-encoded LMP1 expression underpin constitutive NF-κB activation in 90% of NPCs. A similar spectrum of somatic aberrations and viral gene expression undermine innate immunity in 79% of cases and adaptive immunity in 47% of cases; mechanisms by which NPC may evade immune surveillance despite its pro-inflammatory phenotype. Additionally, genomic changes impairing TGFBR2 promote oncogenesis and stabilize EBV infection in tumor cells. Fine-mapping of CDKN2A/CDKN2B deletion breakpoints reveals homozygous MTAP deletions in 32-34% of NPCs that confer marked sensitivity to MAT2A inhibition. Our work concludes that NPC is a homogeneously NF-κB-driven and immune-protected, yet potentially druggable, cancer.