In spite of the growth of molecular ecology, systematics and next-generation sequencing, the discovery and analysis of diversity is not currently integrated with building the tree-of-life. Tropical arthropod ecologists are well placed to accelerate this process if all specimens obtained through mass-trapping, many of which will be new species, could be incorporated routinely into phylogeny reconstruction. Here we test a shotgun sequencing approach, whereby mitochondrial genomes are assembled from complex ecological mixtures through mitochondrial metagenomics, and demonstrate how the approach overcomes many of the taxonomic impediments to the study of biodiversity. DNA from approximately 500 beetle specimens, originating from a single rainforest canopy fogging sample from Borneo, was pooled and shotgun sequenced, followed by de novo assembly of complete and partial mitogenomes for 175 species. The phylogenetic tree obtained from this local sample was highly similar to that from existing mitogenomes selected for global coverage of major lineages of Coleoptera. When all sequences were combined only minor topological changes were induced against this reference set, indicating an increasingly stable estimate of coleopteran phylogeny, while the ecological sample expanded the tip-level representation of several lineages. Robust trees generated from ecological samples now enable an evolutionary framework for ecology. Meanwhile, the inclusion of uncharacterized samples in the tree-of-life rapidly expands taxon and biogeographic representation of lineages without morphological identification. Mitogenomes from shotgun sequencing of unsorted environmental samples and their associated metadata, placed robustly into the phylogenetic tree, constitute novel DNA "superbarcodes" for testing hypotheses regarding global patterns of diversity.
Backcross breeding is the most commonly used method for incorporating a blast resistance gene into a rice cultivar. Linkage between the resistance gene and undesirable units can persist for many generations of backcrossing. Marker-assisted backcrossing (MABC) along with marker-assisted selection (MAS) contributes immensely to overcome the main limitation of the conventional breeding and accelerates recurrent parent genome (RPG) recovery. The MABC approach was employed to incorporate (a) blast resistance gene(s) from the donor parent Pongsu Seribu 1, the blast-resistant local variety in Malaysia, into the genetic background of MR219, a popular high-yielding rice variety that is blast susceptible, to develop a blast-resistant MR219 improved variety. In this perspective, the recurrent parent genome recovery was analyzed in early generations of backcrossing using simple sequence repeat (SSR) markers. Out of 375 SSR markers, 70 markers were found polymorphic between the parents, and these markers were used to evaluate the plants in subsequent generations. Background analysis revealed that the extent of RPG recovery ranged from 75.40% to 91.3% and from 80.40% to 96.70% in BC1F1 and BC2F1 generations, respectively. In this study, the recurrent parent genome content in the selected BC2F2 lines ranged from 92.7% to 97.7%. The average proportion of the recurrent parent in the selected improved line was 95.98%. MAS allowed identification of the plants that are more similar to the recurrent parent for the loci evaluated in backcross generations. The application of MAS with the MABC breeding program accelerated the recovery of the RP genome, reducing the number of generations and the time for incorporating resistance against rice blast.
Enterovirus 71 (EV 71) is a causative agent of mild Hand Foot and Mouth Disease but is capable of causing severe complications in the CNS in young children. Reverse genetics technology is currently widely used to study the pathogenesis of the virus. The aim of this work was to determine and evaluate the factors which can contribute to infectivity of EV 71 RNA transcripts in vitro. Two strategies, overlapping RT-PCR and long distance RT-PCR, were employed to obtain the full-length genome cDNA clones of the virus. The length of the poly(A) tail and the presence of non-viral 3'-terminal sequences were studied in regard to their effects on infectivity of the in vitro RNA transcripts of EV 71 in cell culture. The data revealed that only cDNA clones obtained after long distance RT-PCR were infectious. No differences were observed in virus titres after transfection with in vitro RNA harbouring a poly(A) tail of 18 or 30 adenines in length, irrespective of the non-viral sequences at the 3'-terminus.
Southeast Asia houses various culturally and linguistically diverse ethnic groups. In Malaysia, where the Malay, Chinese, and Indian ethnic groups form the majority, there exist minority groups such as the "negritos" who are believed to be descendants of the earliest settlers of Southeast Asia. Here we report patterns of genetic substructure and admixture in two Malaysian negrito populations (Jehai and Kensiu), using ~50,000 genome-wide single-nucleotide polymorphism (SNP) data. We found traces of recent admixture in both the negrito populations, particularly in the Jehai, with the Malay through principal component analysis and STRUCTURE analysis software, which suggested that the admixture was as recent as one generation ago. We also identified significantly differentiated nonsynonymous SNPs and haplotype blocks related to intracellular transport, metabolic processes, and detection of stimulus. These results highlight the different levels of admixture experienced by the two Malaysian negritos. Delineating admixture and differentiated genomic regions should be of importance in designing and interpretation of molecular anthropology and disease association studies.
Bacterial heart rot of pineapple reported in Hawaii in 2003 and reoccurring in 2006 was caused by an undetermined species of Dickeya. Classification of the bacterial strains isolated from infected pineapple to one of the recognized Dickeya species and their phylogenetic relationships with Dickeya were determined by a multilocus sequence analysis (MLSA), based on the partial gene sequences of dnaA, dnaJ, dnaX, gyrB and recN. Individual and concatenated gene phylogenies revealed that the strains form a clade with reference Dickeya sp. isolated from pineapple in Malaysia and are closely related to D. zeae; however, previous DNA-DNA reassociation values suggest that these strains do not meet the genomic threshold for consideration in D. zeae, and require further taxonomic analysis. An analysis of the markers used in this MLSA determined that recN was the best overall marker for resolution of species within Dickeya. Differential intraspecies resolution was observed with the other markers, suggesting that marker selection is important for defining relationships within a clade. Phylogenies produced with gene sequences from the sequenced genomes of strains D. dadantii Ech586, D. dadantii Ech703 and D. zeae Ech1591 did not place the sequenced strains with members of other well-characterized members of their respective species. The average nucleotide identity (ANI) and tetranucleotide frequencies determined for the sequenced strains corroborated the results of the MLSA that D. dadantii Ech586 and D. dadantii Ech703 should be reclassified as Dickeya zeae Ech586 and Dickeya paradisiaca Ech703, respectively, whereas D. zeae Ech1591 should be reclassified as Dickeya chrysanthemi Ech1591.
The ratio of genetic diversity on X chromosomes relative to autosomes in organisms with XX/XY sex chromosomes could provide fundamental insight into the process of genome evolution. Here we report this ratio for 24 cynomolgus monkeys (Macaca fascicularis) originating in Indonesia, Malaysia, and the Philippines. The average X/A diversity ratios in these samples was 0.34 and 0.20 in the Indonesian-Malaysian and Philippine populations, respectively, considerably lower than the null expectation of 0.75. A Philippine population supposed to derive from an ancestral population by founding events showed a significantly lower ratio than the parental population, suggesting a demographic effect for the reduction. Taking sex-specific mutation rate bias and demographic effect into account, expected X/A diversity ratios generated by computer simulations roughly agreed with the observed data in the intergenic regions. In contrast, silent sites in genic regions on X chromosomes showed strong reduction in genetic diversity and the observed X/A diversity ratio in the genic regions cannot be explained by mutation rate bias and demography, indicating that natural selection also reduces the level of polymorphism near genes. Whole-genome analysis of a female cynomolgus monkey also supported the notion of stronger reduction of genetic diversity near genes on the X chromosome.
Machine learning techniques are becoming useful as an alternative approach to conventional medical diagnosis or prognosis as they are good for handling noisy and incomplete data, and significant results can be attained despite a small sample size. Traditionally, clinicians make prognostic decisions based on clinicopathologic markers. However, it is not easy for the most skilful clinician to come out with an accurate prognosis by using these markers alone. Thus, there is a need to use genomic markers to improve the accuracy of prognosis. The main aim of this research is to apply a hybrid of feature selection and machine learning methods in oral cancer prognosis based on the parameters of the correlation of clinicopathologic and genomic markers.
Background. The gram-negative Neisseria is associated with two of the most potent human epidemic diseases: meningococcal meningitis and gonorrhoea. In both cases, disease is caused by bacteria colonizing human mucosal membrane surfaces. Overall, the genus shows great diversity and genetic variation mainly due to its ability to acquire and incorporate genetic material from a diverse range of sources through horizontal gene transfer. Although a number of databases exist for the Neisseria genomes, they are mostly focused on the pathogenic species. In this present study we present the freely available NeisseriaBase, a database dedicated to the genus Neisseria encompassing the complete and draft genomes of 15 pathogenic and commensal Neisseria species. Methods. The genomic data were retrieved from National Center for Biotechnology Information (NCBI) and annotated using the RAST server which were then stored into the MySQL database. The protein-coding genes were further analyzed to obtain information such as calculation of GC content (%), predicted hydrophobicity and molecular weight (Da) using in-house Perl scripts. The web application was developed following the secure four-tier web application architecture: (1) client workstation, (2) web server, (3) application server, and (4) database server. The web interface was constructed using PHP, JavaScript, jQuery, AJAX and CSS, utilizing the model-view-controller (MVC) framework. The in-house developed bioinformatics tools implemented in NeisseraBase were developed using Python, Perl, BioPerl and R languages. Results. Currently, NeisseriaBase houses 603,500 Coding Sequences (CDSs), 16,071 RNAs and 13,119 tRNA genes from 227 Neisseria genomes. The database is equipped with interactive web interfaces. Incorporation of the JBrowse genome browser in the database enables fast and smooth browsing of Neisseria genomes. NeisseriaBase includes the standard BLAST program to facilitate homology searching, and for Virulence Factor Database (VFDB) specific homology searches, the VFDB BLAST is also incorporated into the database. In addition, NeisseriaBase is equipped with in-house designed tools such as the Pairwise Genome Comparison tool (PGC) for comparative genomic analysis and the Pathogenomics Profiling Tool (PathoProT) for the comparative pathogenomics analysis of Neisseria strains. Discussion. This user-friendly database not only provides access to a host of genomic resources on Neisseria but also enables high-quality comparative genome analysis, which is crucial for the expanding scientific community interested in Neisseria research. This database is freely available at http://neisseria.um.edu.my.
A unique feature of influenza A virus (IAV) life cycle is replication of the viral genome in the host cell nucleus. The nuclear import of IAV genome is an indispensable step in establishing virus infection. IAV nucleoprotein (NP) is known to mediate the nuclear import of viral genome via its nuclear localization signals. Here, we demonstrate that cellular heat shock protein 40 (Hsp40/DnaJB1) facilitates the nuclear import of incoming IAV viral ribonucleoproteins (vRNPs) and is important for efficient IAV replication. Hsp40 was found to interact with NP component of IAV RNPs during early stages of infection. This interaction is mediated by the J domain of Hsp40 and N-terminal region of NP. Drug or RNAi mediated inhibition of Hsp40 resulted in reduced nuclear import of IAV RNPs, diminished viral polymerase function and attenuates overall viral replication. Hsp40 was also found to be required for efficient association between NP and importin alpha, which is crucial for IAV RNP nuclear translocation. These studies demonstrate an important role for cellular chaperone Hsp40/DnaJB1 in influenza A virus life cycle by assisting nuclear trafficking of viral ribonucleoproteins.
Briskly evolving phytopathogens are dire threats to our food supplies and threaten global food security. From the recent advances made toward high-throughput sequencing technologies, understanding of pathogenesis and effector biology, and plant innate immunity, translation of these means into new control tools is being introduced to develop durable disease resistance. Effectoromics as a powerful genetic tool for uncovering effector-target genes, both susceptibility genes and executor resistance genes in effector-assisted breeding, open up new avenues to improve resistance. TALENs (Transcription Activator-Like Effector Nucleases), engineered nucleases and CRISPR (Clustered Regulatory Interspaced Short Palindromic Repeats)/Cas9 systems are breakthrough and powerful techniques for genome editing, providing efficient mechanisms for targeted crop protection strategies in disease resistance programs. In this review, major advances in plant disease management to confer durable disease resistance and novel strategies for boosting plant innate immunity are highlighted.
Extraintestinal pathogenic Escherichia coli (ExPEC) that causes extraintestinal infections often harbor plasmids encoding fitness traits such as resistance and virulence determinants that are of clinical importance. We determined the complete nucleotide sequence of plasmid pEC302/04 from a multidrug-resistant E. coli EC302/04 which was isolated from the tracheal aspirate of a patient in Malaysia. In addition, we also performed comparative sequence analyses of 18 related IncFIIA plasmids to determine the phylogenetic relationship and diversity of these plasmids. The 140,232 bp pEC302/04 is a multireplicon plasmid that bears three replication systems (FII, FIA, and FIB) with subtype of F2:A1:B1. The plasmid is self-transmissible with a complete transfer region. pEC302/04 also carries antibiotic resistance genes such as bla TEM-1 and a class I integron containing sul1, cml and aadA resistance genes, conferring multidrug resistance (MDR) to its host, E. coli EC302/04. Besides, two iron acquisition systems (SitABCD and IutA-IucABCD) which are the conserved virulence determinants of ExPEC-colicin V or B and M (ColV/ColBM)-producing plasmids were identified in pEC302/04. Multiple toxin-antitoxin (TA)-based addiction systems (i.e., PemI/PemK, VagC/VagD, CcdA/CcdB, and Hok/Sok) and a plasmid partitioning system, ParAB, and PsiAB, which are important for plasmid maintenance were also found. Comparative plasmid analysis revealed only one conserved gene, the repA1 as the core genome, showing that there is an extensive diversity among the IncFIIA plasmids. The phylogenetic relationship of 18 IncF plasmids based on the core regions revealed that ColV/ColBM-plasmids and non-ColV/ColBM plasmids were separated into two distinct groups. These plasmids, which carry highly diverse genetic contents, are also mosaic in nature. The atypical combination of genetic materials, i.e., the MDR- and ColV/ColBM-plasmid-virulence encoding regions in a single ExPEC plasmid is rare but of clinical importance. Such phenomenon is bothersome when the plasmids are transmissible, facilitating the spread of virulence and resistance plasmids among pathogenic bacteria. Notably, certain TA systems are more commonly found in particular ExPEC plasmid types, indicating the possible relationships between certain TA systems and ExPEC pathogenesis.
Background. Two non-tuberculous mycobacterial strains, UM_3 and UM_11, were isolated from the trunk wash of captive elephants in Malaysia. As they appeared to be identical phenotypes, they were investigated further by conventional and whole genome sequence-based methods of strain differentiation. Methods. Multiphasic investigations on the isolates included species identification with hsp65 PCR-sequencing, conventional biochemical tests, rapid biochemical profiling using API strips and the Biolog Phenotype Microarray analysis, protein profiling with liquid chromatography-mass spectrometry, repetitive sequence-based PCR typing and whole genome sequencing followed by phylogenomic analyses. Results. The isolates were shown to be possibly novel slow-growing schotochromogens with highly similar biological and genotypic characteristics. Both strains have a genome size of 5.2 Mbp, G+C content of 68.8%, one rRNA operon and 52 tRNAs each. They qualified for classification into the same species with their average nucleotide identity of 99.98% and tetranucleotide correlation coefficient of 0.99999. At the subspecies level, both strains showed 98.8% band similarity in the Diversilab automated repetitive sequence-based PCR typing system, 96.2% similarity in protein profiles obtained by liquid chromatography mass spectrometry, and a genomic distance that is close to zero in the phylogenomic tree constructed with conserved orthologs. Detailed epidemiological tracking revealed that the elephants shared a common habitat eight years apart, thus, strengthening the possibility of a clonal relationship between the two strains.
Pandoraea thiooxydans DSM 25325(T) is a thiosulfate-oxidizing bacterium isolated from rhizosphere soils of a sesame plant. Here, we present the first complete genome of P. thiooxydans DSM 25325(T). Several genes involved in thiosulfate oxidation and biodegradation of aromatic compounds were identified.
Dickeya solani is an emerging pathogen that causes soft rot and blackleg diseases in several crops including Solanum tuberosum, but little is known about its genomic diversity and evolution.
The whole mitochondrial genome of the pest fruit fly Bactrocera arecae was obtained from next-generation sequencing of genomic DNA. It had a total length of 15,900 bp, consisting of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and a non-coding region (A + T-rich control region). The control region (952 bp) was flanked by rrnS and trnI genes. The start codons included 6 ATG, 3 ATT and 1 each of ATA, ATC, GTG and TCG. Eight TAA, two TAG, one incomplete TA and two incomplete T stop codons were represented in the protein-coding genes. The cloverleaf structure for trnS1 lacked the D-loop, and that of trnN and trnF lacked the TΨC-loop. Molecular phylogeny based on 13 protein-coding genes was concordant with 37 mitochondrial genes, with B. arecae having closest genetic affinity to B. tryoni. The subgenus Bactrocera of Dacini tribe and the Dacinae subfamily (Dacini and Ceratitidini tribes) were monophyletic. The whole mitogenome of B. arecae will serve as a useful dataset for studying the genetics, systematics and phylogenetic relationships of the many species of Bactrocera genus in particular, and tephritid fruit flies in general.
Genome-wide association studies have identified 20 genomic regions associated with risk of epithelial ovarian cancer (EOC), but many additional risk variants may exist. Here, we evaluated associations between common genetic variants [single nucleotide polymorphisms (SNPs) and indels] in DNA repair genes and EOC risk. We genotyped 2896 common variants at 143 gene loci in DNA samples from 15 397 patients with invasive EOC and controls. We found evidence of associations with EOC risk for variants at FANCA, EXO1, E2F4, E2F2, CREB5 and CHEK2 genes (P ≤ 0.001). The strongest risk association was for CHEK2 SNP rs17507066 with serous EOC (P = 4.74 x 10(-7)). Additional genotyping and imputation of genotypes from the 1000 genomes project identified a slightly more significant association for CHEK2 SNP rs6005807 (r (2) with rs17507066 = 0.84, odds ratio (OR) 1.17, 95% CI 1.11-1.24, P = 1.1×10(-7)). We identified 293 variants in the region with likelihood ratios of less than 1:100 for representing the causal variant. Functional annotation identified 25 candidate SNPs that alter transcription factor binding sites within regulatory elements active in EOC precursor tissues. In The Cancer Genome Atlas dataset, CHEK2 gene expression was significantly higher in primary EOCs compared to normal fallopian tube tissues (P = 3.72×10(-8)). We also identified an association between genotypes of the candidate causal SNP rs12166475 (r (2) = 0.99 with rs6005807) and CHEK2 expression (P = 2.70×10(-8)). These data suggest that common variants at 22q12.1 are associated with risk of serous EOC and CHEK2 as a plausible target susceptibility gene.
A vaccine against human enterovirus 71 (EV-A71) is urgently needed to combat outbreaks of EV-A71 and in particular, the serious neurological complications that manifest during these outbreaks. In this study, an EV-A71 virus-like-particle (VLP) based on a B5 subgenogroup (EV-A71-B5 VLP) was generated using an insect cell/baculovirus platform. Biochemical analysis demonstrated that the purified VLP had a highly native procapsid structure and initial studies in vivo demonstrated that the VLPs were immunogenic in mice. The impact of VLP immunization on infection was examined in non-human primates using a VLP prime-boost strategy prior to EV-A71 challenge. Rhesus macaques were immunized on day 0 and day 21 with VLPs (100 μg/dose) containing adjuvant or with adjuvant alone (controls), and were challenged with EV-A71 on day 42. Complete blood counts, serum chemistry, magnetic resonance imaging (MRI) scans, and histopathology results were mostly normal in vaccinated and control animals after virus challenge demonstrating that the fatal EV-A71-B3 clinical isolate used in this study was not highly virulent in rhesus macaques. Viral genome and/or infectious virus were detected in blood, spleen or brain of two of three control animals, but not in any specimens from the vaccinated animals, indicating that VLP immunization prevented systemic spread of EV-A71 in rhesus macaques. High levels of IgM and IgG were detected in VLP-vaccinated animals and these responses were highly specific for EV-A71 particles and capsid proteins. Serum from vaccinated animals also exhibited similar neutralizing activity against different subgenogroups of EV-A71 demonstrating that the VLPs induced cross-neutralizing antibodies. In conclusion, our EV-A71-B5 VLP is safe, highly immunogenic, and prevents systemic EV-A71-B3 infection in nonhuman primates making it a viable attractive vaccine candidate for EV-A71.
Viral infections are a common cause of morbidity worldwide. The emergence of Coronavirus Disease 2019 (COVID-19) has led to more attention to viral infections and finding novel therapeutics. The CRISPR-Cas9 system has been recently proposed as a potential therapeutic tool for the treatment of viral diseases. Here, we review the research progress in the use of CRISPR-Cas technology for treating viral infections, as well as the strategies for improving the delivery of this gene-editing tool in vivo. Key challenges that hinder the widespread clinical application of CRISPR-Cas9 technology are also discussed, and several possible directions for future research are proposed.
The bone morphogenic protein (BMP) family is a member of the TGF-beta superfamily and plays a crucial role during the onset of gut inflammation and arthritis diseases. Recent studies have reported a connection with the gut-joint axis; however, the genetic players are still less explored. Meanwhile, BDMC33 is a newly synthesized anti-inflammatory drug candidate. Therefore, in our present study, we analysed the genome-wide features of the BMP family as well as the role of BMP members in gut-associated arthritis in an inflammatory state and the ability of BDMC33 to attenuate this inflammation. Firstly, genome-wide analyses were performed on the BMP family in the zebrafish genome, employing several in silico techniques. Afterwards, the effects of curcumin analogues on BMP gene expression in zebrafish larvae induced with TNBS (0.78 mg/mL) were determined using real time-qPCR. A total of 38 identified BMP proteins were revealed to be clustered in five major clades and contain TGF beta and TGF beta pro peptide domains. Furthermore, BDMC33 suppressed the expression of four selected BMP genes in the TNBS-induced larvae, where the highest gene suppression was in the BMP2a gene (an eight-fold decrement), followed by BMP7b (four-fold decrement), BMP4 (four-fold decrement), and BMP6 (three-fold decrement). Therefore, this study reveals the role of BMPs in gut-associated arthritis and proves the ability of BDMC33 to act as a potential anti-inflammatory drug for suppressing TNBS-induced BMP genes in zebrafish larvae.
Nosocomial outbreaks caused by carbapenem-resistant Acinetobacter baumannii (CRAB) strains are rapidly emerging worldwide and are cause for concern. Herein, we aimed to describe the genomic characteristics of CRAB strains isolated from two hospitals in China in 2023. The A. baumannii isolates were mainly collected from the ICU and isolated from the sputum (71.43%, 15/21), followed by urine (14.29%, 3/21). Twenty-one A. baumannii strains possessed a multidrug-resistant (MDR) profile, and whole-genome sequencing showed that they all carried blaOXA-23. Based on the Pasteur multilocus sequence typing (MLST) scheme, all strains were typed into a sequence type 2 (ST2). Based on the Oxford MLST scheme, six strains belonged to ST540, three of which were ST208, and four strains were assigned to ST784. Kaptive showed most of the strains (38.10%, 8/21) contained KL93. As for the lipoolygosaccharide (OC locus) type, OCL1c and OCL1d were identified, accounting for 33.33% (7/21) and 66.67% (14/21), respectively. Based on the BacWGSTdb server, we found that the strains belonging to ST540 and ST784 were all collected from China. However, the ST938 strains were isolated from Malaysia and Thailand. Comparative genomics analysis showed that the AB10 strain had a closed relationship with SXAB10-SXAB13 strains, suggesting the transmission happened in these two hospitals and other hospital in China. In addition, the 4300STDY7045869 strain, which was collected from Thailand, possessed near genetic relationship with our isolates in this study, suggesting the possible spread among various countries. Additionally, 3-237 single nucleotide polymorphisms were observed among these strains. In conclusion, this study conducted a genome-based study for A. baumannii strains collected from two hospitals in China and revealed their epidemiological and molecular features. Clone spreading occurred in these two hospitals. Hence, there is an urgent need for increased surveillance in hospitals and other clinical settings to prevent and control CRAB spreading.