Topographically Distinguished Microbiome Taxonomy and Stress-Response Genes of Royal Belum Rainforest and Raja Muda Musa Peat Swamp Revealed through Metagenomic Inquisition.

Soil ecosystems are home to a diverse range of microorganisms, but they are only partially understood because no single-cell sequencing or whole-community sequencing provides a complete picture of these complex communities. Using one of such metagenomics approaches, we succeeded in monitoring the microbial diversity and stress-response gene in the soil samples. This study aims to test whether known differences in taxonomic diversity and composition are reflected in functional gene profiles by implementing whole gene sequencing (WGS) metagenomic analysis of geographically dispersed soils from two distinct pristine forests. The study was commenced by sequencing three rainforest soil samples and three peat swamp soil samples. Soil richness effects were assessed by exploring the changes in specific functional gene abundances to elucidate physiological constraints acting on different soil systems and identify variance in functional pathways relevant to soil biogeochemical cycling. Proteobacteria shows abundances of microbial diversity for 52.15% in Royal Belum Reserved Forest and 48.28% in Raja Musa; 177 out of 1,391,841 and 449 out of 3,586,577 protein coding represent acidic stress-response genes for Royal Belum and Raja Musa, respectively. Raja Musa indicates pH 2.5, which is extremely acidic. The analysis of the taxonomic community showed that Royal Belum soils are dominated by bacteria (98% in Sungai Kooi (SK), 98% in Sungai Papan (SP), and 98% in Sungai Ruok (SR), Archaea (0.9% in SK, 0.9% in SP, and 1% in SR), and the remaining were classed under Eukaryota and viruses. Likewise, the soils of Raja Muda Musa are also dominated by bacteria (95% in Raja Musa 1 (RM1), 98% in Raja Musa 2 (RM2), and 96% in Raja Musa 3 (RM3)), followed by Archaea (4% in RM1, 1% in RM2, and 3% in RM3), and the remaining were classed under Eukaryota and viruses. This study revealed that RBFR (Royal Belum Foresr Reserve) and RMFR (Raja Musa Forest Reserve) metagenomes contained abundant stress-related genes assigned to various stress-response pathways, many of which did not show any difference among samples from both sites. Our findings indicate that the structure and functional potential of the microbial community will be altered by future environmental potential as the first glimpse of both the taxonomic and functional composition of soil microbial communities.

The molecular epidemiology of multiple zoonotic origins of SARS-CoV-2.

Understanding the circumstances that lead to pandemics is important for their prevention. We analyzed the genomic diversity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) early in the coronavirus disease 2019 (COVID-19) pandemic. We show that SARS-CoV-2 genomic diversity before February 2020 likely comprised only two distinct viral lineages, denoted "A" and "B." Phylodynamic rooting methods, coupled with epidemic simulations, reveal that these lineages were the result of at least two separate cross-species transmission events into humans. The first zoonotic transmission likely involved lineage B viruses around 18 November 2019 (23 October to 8 December), and the separate introduction of lineage A likely occurred within weeks of this event. These findings indicate that it is unlikely that SARS-CoV-2 circulated widely in humans before November 2019 and define the narrow window between when SARS-CoV-2 first jumped into humans and when the first cases of COVID-19 were reported. As with other coronaviruses, SARS-CoV-2 emergence likely resulted from multiple zoonotic events.

Proteomic analysis revealed the biofilm-degradation abilities of the bacteriophage UPMK_1 and UPMK_2 against Methicillin-resistant Staphylococcus aureus.

OBJECTIVE: The degradation activity of two bacteriophages UPMK_1 and UPMK_2 against methicillin-resistant Staphylococcus aureus phages were examined using gel zymography. METHODS: The analysis was done using BLASTP to detect peptides catalytic domains. Many peptides that are related to several phage proteins were revealed. RESULTS: UPMK_1 and UPMK_2 custom sequence database were used for peptide identification. The biofilm-degrading proteins in the bacteriophage UPMK_2 revealed the same lytic activity towards polysaccharide intercellular adhesin-dependent and independent of Methicillin-resistant Staphylococcus aureus (MRSA) biofilm producers in comparison to UPMK_1, which had lytic activity restricted solely to its host. CONCLUSION: Both bacteriophage enzymes were involved in MRSA biofilm degradation during phage infection and they have promising enzybiotics properties against MRSA biofilm formation.

Transcriptome analysis of Rafflesia cantleyi flower stages reveals insights into the regulation of senescence.

Rafflesia is a unique plant species existing as a single flower and produces the largest flower in the world. While Rafflesia buds take up to 21 months to develop, its flowers bloom and wither within about a week. In this study, transcriptome analysis was carried out to shed light on the molecular mechanism of senescence in Rafflesia. A total of 53.3 million high quality reads were obtained from two Rafflesia cantleyi flower developmental stages and assembled to generate 64,152 unigenes. Analysis of this dataset showed that 5,166 unigenes were differentially expressed, in which 1,073 unigenes were identified as genes involved in flower senescence. Results revealed that as the flowers progress to senescence, more genes related to flower senescence were significantly over-represented compared to those related to plant growth and development. Senescence of the R. cantleyi flower activates senescence-associated genes in the transcription activity (members of the transcription factor families MYB, bHLH, NAC, and WRKY), nutrient remobilization (autophagy-related protein and transporter genes), and redox regulation (CATALASE). Most of the senescence-related genes were found to be differentially regulated, perhaps for the fine-tuning of various responses in the senescing R. cantleyi flower. Additionally, pathway analysis showed the activation of genes such as ETHYLENE RECEPTOR, ETHYLENE-INSENSITIVE 2, ETHYLENE-INSENSITIVE 3, and ETHYLENE-RESPONSIVE TRANSCRIPTION FACTOR, indicating the possible involvement of the ethylene hormone response pathway in the regulation of R. cantleyi senescence. Our results provide a model of the molecular mechanism underlying R. cantleyi flower senescence, and contribute essential information towards further understanding the biology of the Rafflesiaceae family.

Cellular and molecular profiles of anterior nervous system regeneration in Diopatra claparedii Grube, 1878 (Annelida, Polychaeta).

The polychaete Diopatra claparedii Grube, 1878 is among those organisms successfully carrying out full body regeneration, including the whole nervous system. Thus, D. claparedii potentially can be regarded for the nervous system regeneration (NSR) study. However, data on the property of its nervous system and the NSR profile are still lacking. In this study, we investigated the morphology of D. claparedii anterior nervous system (ANS) and examined the cellular and molecular profiles on its early anterior NSR. The nervous system of D. claparedii consists of a symmetry brain with nerves branching off, circumpharyngeal connectives that connect the brain and nerve cord as well as obvious segmental ganglia. Moreover, we identified changes in the cellular condition of the ganglionic cells in the regenerating tissue, such as the accumulation of lysosomes and lipofuscins, elongated mitochondria and multiple nucleoli. Furthermore, mRNA of tissues at two regenerating stages, as well as intact tissue (non-regenerating), were sequenced with Illumina sequencer. We identified from these tissues 37,248 sequences, 18 differential expressed proteins of which upregulated were involved in NSR with noelin-like isoform X2 turned up to be the highest being expressed. Our results highlight the cellular and molecular changes during early phase of NSR, thus providing essential insights on regeneration within Annelida and understanding the neurodegenerative diseases.

Transcriptome landscape of Rafflesia cantleyi floral buds reveals insights into the roles of transcription factors and phytohormones in flower development.

Rafflesia possesses unique biological features and known primarily for producing the world's largest and existing as a single flower. However, to date, little is known about key regulators participating in Rafflesia flower development. In order to further understand the molecular mechanism that regulates Rafflesia cantleyi flower development, RNA-seq data from three developmental stages of floral bud, representing the floral organ primordia initiation, floral organ differentiation, and floral bud outgrowth, were analysed. A total of 89,890 transcripts were assembled of which up to 35% could be annotated based on homology search. Advanced transcriptome analysis using K-mean clustering on the differentially expressed genes (DEGs) was able to identify 12 expression clusters that reflect major trends and key transitional states, which correlate to specific developmental stages. Through this, comparative gene expression analysis of different floral bud stages identified various transcription factors related to flower development. The members of WRKY, NAC, bHLH, and MYB families are the most represented among the DEGs, suggesting their important function in flower development. Furthermore, pathway enrichment analysis also revealed DEGs that are involved in various phytohormone signal transduction events such as auxin and auxin transport, cytokinin and gibberellin biosynthesis. Results of this study imply that transcription factors and phytohormone signalling pathways play major role in Rafflesia floral bud development. This study provides an invaluable resource for molecular studies of the flower development process in Rafflesia and other plant species.

Comparative analysis of nucleus-encoded plastid-targeting proteins in Rafflesia cantleyi against photosynthetic and non-photosynthetic representatives reveals orthologous systems with potentially divergent functions.

Parasitic plants are known to discard photosynthesis thus leading to the deletion or loss of the plastid genes. Despite plastid genome reduction in non-photosynthetic plants, some nucleus-encoded proteins are transported back to the plastid to carry out specific functions. In this work, we study such proteins in Rafflesia cantleyi, a member of the holoparasitic genus well-known for producing the largest single flower in the world. Our analyses of three transcriptome datasets, two holoparasites (R. cantleyi and Phelipanche aegyptiaca) and one photosynthetic plant (Arabidopsis thaliana), suggest that holoparasites, such as R. cantleyi, retain some common plastid associated processes such as biosynthesis of amino acids and lipids, but are missing photosynthesis components that can be extensions of these pathways. The reconstruction of two selected biosynthetic pathways involving plastids correlates the trend of plastid retention to pathway complexity - transcriptome evidence for R. cantleyi suggests alternate mechanisms in regulating the plastidial heme and terpenoid backbone biosynthesis pathways. The evolution to holoparasitism from autotrophy trends towards devolving the plastid genes to the nuclear genome despite the functional sites remaining in the plastid, or maintaining non-photosynthetic processes in the plastid, before the eventual loss of the plastid and any site dependent functions.

Unravelling the adaptation strategies employed by Glaciozyma antarctica PI12 on Antarctic sea ice.

Glaciozyma antarctica PI12, is a psychrophilic yeast isolated from Antarctic sea. In this work, Expressed Sequence Tags (EST) from cells exposed to three different temperatures; 15 °C, 0 °C and -12 °C were generated to identify genes associated with cold adaptation. A total of 5376 clones from each library were randomly picked and sequenced. Comparative analyses from the resulting ESTs in each condition identified several groups of genes required for cold adaptation. Additionally, 319 unique transcripts that encoded uncharacterised functions were identified in the -12 °C library and are currently unique to G. antarctica. Gene expression analysis using RT-qPCR revealed two of the unknown genes to be up-regulated at -12 °C compared to 0 °C and 15 °C. These findings further contribute to the collective knowledge into G. antarctica cold adaptation and as a resource for understanding the ecological and physiological tolerance of psychrophilic microbes in general.

Genome Sequence of Serratia marcescens subsp. sakuensis Strain K27, a Marine Bacterium Isolated from Sponge (Haliclona amboinensis).

Serratia marcescens subsp. sakuensis strain K27 was isolated from sponge (Haliclona amboinensis). The genome of this strain consists of 5,325,727 bp, with 5,140 open reading frames (ORFs), 3 rRNAs, and 67 tRNAs. It contains genes for the production of amylases, lipases, and proteases. Gene clusters for the biosynthesis of nonribosomal peptides and thiopeptide were also identified.

The Glaciozyma antarctica genome reveals an array of systems that provide sustained responses towards temperature variations in a persistently cold habitat.

Extremely low temperatures present various challenges to life that include ice formation and effects on metabolic capacity. Psyhcrophilic microorganisms typically have an array of mechanisms to enable survival in cold temperatures. In this study, we sequenced and analysed the genome of a psychrophilic yeast isolated in the Antarctic region, Glaciozyma antarctica. The genome annotation identified 7857 protein coding sequences. From the genome sequence analysis we were able to identify genes that encoded for proteins known to be associated with cold survival, in addition to annotating genes that are unique to G. antarctica. For genes that are known to be involved in cold adaptation such as anti-freeze proteins (AFPs), our gene expression analysis revealed that they were differentially transcribed over time and in response to different temperatures. This indicated the presence of an array of adaptation systems that can respond to a changing but persistent cold environment. We were also able to validate the activity of all the AFPs annotated where the recombinant AFPs demonstrated anti-freeze capacity. This work is an important foundation for further collective exploration into psychrophilic microbiology where among other potential, the genes unique to this species may represent a pool of novel mechanisms for cold survival.

RNA-seq data of Oryza sativa cultivar Kuku Belang under PEG treatment.

Drought stress is the main abiotic factor affecting rice production. Rain-fed upland rice which is grown on unbounded fields and totally dependent on rainfall for moisture is more prone to drought stress compared to rice from other ecosystems. However, upland rice has adapted to this limited water condition, thus are more drought tolerant than rice from other ecosystems. We performed the first transcriptome sequencing of drought tolerant indica upland rice cultivar Kuku Belang to identify differentially expressed genes related to drought tolerance mechanism. Raw reads for non-treated and PEG-treated Oryza sativa subspecies indica cv. Kuku Belang were deposited in the NCBI SRA database with accession number SRP074520 (https://www.ncbi.nlm.nih.gov/sra?term=SRP074520).

RNA-seq data from different developmental stages of Rafflesia cantleyi floral buds.

Rafflesia cantleyi, known as one of the world's largest flowers, is a specialised holoparasite due to dramatic morphological modifications. It possesses highly reduced vegetative structure and only appears as a flower for sexual reproduction. Moreover, it has an unusual life cycle in that its floral bud development takes up to nine months. In order to fully understand the highly modified floral organ structure and long life cycle of R. cantleyi, we used Illumina sequencing technology (HiSeq) for sequence generation followed by de novo assembly of sequence reads. We obtained the RNA-seq data from three different stages of floral bud, representing the early, mid and advanced developmental stages. These data are available via BioProject accession number PRJNA378435. More than 10.3 Gb raw sequence data were generated, corresponding to 102,203,042 raw reads. Following removal of low-quality reads and trimming of adapter sequences, a total of 91,638,836 reads were obtained. De novo assembly of these sequences using Trinity resulted in 89,690 unique transcripts with an N50 of 1653 bp. The obtained transcriptomic data will be useful for further study to understand the molecular interactions that result in R. cantleyi floral development.

Expression and characterization of thermostable glycogen branching enzyme from Geobacillus mahadia Geo-05.

The glycogen branching enzyme (EC 2.4.1.18), which catalyses the formation of α-1,6-glycosidic branch points in glycogen structure, is often used to enhance the nutritional value and quality of food and beverages. In order to be applicable in industries, enzymes that are stable and active at high temperature are much desired. Using genome mining, the nucleotide sequence of the branching enzyme gene (glgB) was extracted from the Geobacillus mahadia Geo-05 genome sequence provided by the Malaysia Genome Institute. The size of the gene is 2013 bp, and the theoretical molecular weight of the protein is 78.43 kDa. The gene sequence was then used to predict the thermostability, function and the three dimensional structure of the enzyme. The gene was cloned and overexpressed in E. coli to verify the predicted result experimentally. The purified enzyme was used to study the effect of temperature and pH on enzyme activity and stability, and the inhibitory effect by metal ion on enzyme activity. This thermostable glycogen branching enzyme was found to be most active at 55 °C, and the half-life at 60 °C and 70 °C was 24 h and 5 h, respectively. From this research, a thermostable glycogen branching enzyme was successfully isolated from Geobacillus mahadia Geo-05 by genome mining together with molecular biology technique.

Perigone Lobe Transcriptome Analysis Provides Insights into Rafflesia cantleyi Flower Development.

Rafflesia is a biologically enigmatic species that is very rare in occurrence and possesses an extraordinary morphology. This parasitic plant produces a gigantic flower up to one metre in diameter with no leaves, stem or roots. However, little is known about the floral biology of this species especially at the molecular level. In an effort to address this issue, we have generated and characterised the transcriptome of the Rafflesia cantleyi flower, and performed a comparison with the transcriptome of its floral bud to predict genes that are expressed and regulated during flower development. Approximately 40 million sequencing reads were generated and assembled de novo into 18,053 transcripts with an average length of 641 bp. Of these, more than 79% of the transcripts had significant matches to annotated sequences in the public protein database. A total of 11,756 and 7,891 transcripts were assigned to Gene Ontology categories and clusters of orthologous groups respectively. In addition, 6,019 transcripts could be mapped to 129 pathways in Kyoto Encyclopaedia of Genes and Genomes Pathway database. Digital abundance analysis identified 52 transcripts with very high expression in the flower transcriptome of R. cantleyi. Subsequently, analysis of differential expression between developing flower and the floral bud revealed a set of 105 transcripts with potential role in flower development. Our work presents a deep transcriptome resource analysis for the developing flower of R. cantleyi. Genes potentially involved in the growth and development of the R. cantleyi flower were identified and provide insights into biological processes that occur during flower development.

Full genome SNP-based phylogenetic analysis reveals the origin and global spread of Brucella melitensis.

BACKGROUND: Brucellosis is an important zoonotic disease that affects both humans and animals. We sequenced the full genome and characterised the genetic diversity of two Brucella melitensis isolates from Malaysia and the Philippines. In addition, we performed a comparative whole-genome single nucleotide polymorphism (SNP) analysis of B. melitensis strains collected from around the world, to investigate the potential origin and the history of the global spread of B. melitensis. RESULTS: Single sequencing runs of each genome resulted in draft genome sequences of MY1483/09 and Phil1136/12, which covered 99.85% and 99.92% of the complete genome sequences, respectively. The B. melitensis genome sequences, and two B. abortus strains used as the outgroup strains, yielded a total of 13,728 SNP sites. Phylogenetic analysis using whole-genome SNPs and geographical distribution of the isolates revealed spatial clustering of the B. melitensis isolates into five genotypes, I, II, III, IV and V. The Mediterranean strains, identified as genotype I, occupied the basal node of the phylogenetic tree, suggesting that B. melitensis may have originated from the Mediterranean regions. All of the Asian B. melitensis strains clustered into genotype II with the SEA strains, including the two isolates sequenced in this study, forming a distinct clade denoted here as genotype IId. Genotypes III, IV and V of B. melitensis demonstrated a restricted geographical distribution, with genotype III representing the African lineage, genotype IV representing the European lineage and genotype V representing the American lineage. CONCLUSION: We showed that SNPs retrieved from the B. melitensis draft full genomes were sufficient to resolve the interspecies relationships between B. melitensis strains and to discriminate between the vaccine and endemic strains. Phylogeographic reconstruction of the history of B. melitensis global spread at a finer scale by using whole-genome SNP analyses supported the origin of all B. melitensis strains from the Mediterranean region. The possible global distribution of B. melitensis following the ancient trade routes was also consistent with whole-genome SNP phylogeny. The whole genome SNP phylogenetics analysis, hence is a powerful tool for intraspecies discrimination of closely related species.

Draft Genome Sequence of Pasteurella multocida subsp. multocida Strain PMTB, Isolated from a Buffalo.

Pasteurella multocida serotypes B:2 and E:2 are the main causative agents of ruminant hemorrhagic septicemia in Asia and Africa, respectively. Pasteurella multocida strain PMTB was isolated from a buffalo with hemorrhagic septicemia and has been determined to be serotype B:2. Here we report the draft genome sequence of strain PMTB.

EmaxDB: Availability of a first draft genome sequence for the apicomplexan Eimeria maxima.

Apicomplexan parasites are serious pathogens of animals and man that cause diseases including coccidiosis, malaria and toxoplasmosis. The importance of these parasites has prompted the establishment of genomic resources in support of developing effective control strategies. For the Eimeria species resources have developed most rapidly for the reference Eimeria tenella Houghton strain (http://www.genedb.org/Homepage/Etenella). The value of these resources can be enhanced by comparison with related parasites. The well characterised immunogenicity and genetic diversity associated with Eimeria maxima promote its use in genetics-led studies on coccidiosis and recommended its selection for sequencing. Using a combination of sequencing technologies a first draft assembly and annotation has been produced for an E. maxima Houghton strain-derived clone (EmaxDB; http://www.genomemalaysia.gov.my/emaxdb/). The assembly of a draft genome sequence for E. maxima provides a resource for comparative studies with Eimeria and related parasites as demonstrated here through the identification of genes predicted to encode microneme proteins in E. maxima.

Metabolic routes affecting rubber biosynthesis in Hevea brasiliensis latex.

The cytosolic mevalonate (MVA) pathway in Hevea brasiliensis latex is the conventionally accepted pathway which provides isopentenyl diphosphate (IPP) for cis-polyisoprene (rubber) biosynthesis. However, the plastidic 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway may be an alternative source of IPP since its more recent discovery in plants. Quantitative RT-PCR (qRT-PCR) expression profiles of genes from both pathways in latex showed that subcellular compartmentalization of IPP for cis-polyisoprene synthesis is related to the degree of plastidic carotenoid synthesis. From this, the occurrence of two schemes of IPP partitioning and utilization within one species is proposed whereby the supply of IPP for cis-polyisoprene from the MEP pathway is related to carotenoid production in latex. Subsequently, a set of latex unique gene transcripts was sequenced and assembled and they were then mapped to IPP-requiring pathways. Up to eight such pathways, including cis-polyisoprene biosynthesis, were identified. Our findings on pre- and post-IPP metabolic routes form an important aspect of a pathway knowledge-driven approach to enhancing cis-polyisoprene biosynthesis in transgenic rubber trees.